astrocut:切割fitsfile 

from astrocut import fits_cut
from astropy.io import fits
from astropy.coordinates import SkyCoord
input_files=['13co.fits']
center_coord = SkyCoord("285.9131678 5.1650797",unit='deg')
cutout_size=[120,120]
cutout_file = fits_cut(input_files,center_coord,cutout_size,extension='all',cutout_prefix="cutout_p2381x05y09", verbose=False)
plt.imshow(fits.getdata(cutout_file,1))

 

posted @ 2022-01-11 21:02  cql_astro  阅读(49)  评论(0)  编辑  收藏  举报