pythonでフィッティングするファイル（Astropy）

Question

次元（180,360）のデータを持つフィットファイルがあります。これは180行で、各行に360要素あります。私がそれをプロットするとき、x軸の範囲は0から360まで、y軸の範囲は0から180までです。x軸とy軸が実際に緯度であり、経度。私は天体を使用しています。

Answer 1

ファイルと呼ばれるFITS持っていると仮定すると：'filename.fits'を：

from astropy.io import fits 
from astropy.wcs import WCS 

with fits.open('filename.fits') as hdus: 
    data = hdus[0].data # sometimes data is also saved in the extension 1 or "sci" 
    header = hdus[0].header 
    wcs = WCS(hdus[0].header) 

pixel = np.array([[10, 10], [20, 20]]) # pixel 10/10 and 20/20 
world = w.wcs_pix2world(pixcrd, 1) # convert the list of pixel to coordinates 

今、あなたはあなたがdata、headerとwcsで、あなたはその後、新しいファイルを書き込むことができますたければ好きなものをも行うことができます。

hdulist = fits.HDUList([fits.PrimaryHDU(data, header=header)]) 
hdulist.writeto('new_name.fits') 
# hdulist.writeto('filename.fits', clobber=True) # or replace the original 

本質的に同じで詳細な情報が存在するastropy-documentationをご覧ください。

• Loading WCS from a FITS file
• fits.open()
fits.HDUList.writeto()