from astropy.io import fits
from astropy.wcs import WCS
import numpy as np
from astropy import units as u
from astropy.coordinates import SkyCoord
from astropy.utils.data import get_pkg_data_filename
fn = get_pkg_data_filename('4682_1_1_0046965_001_0000500_0_3_20210131210820_00060_000002.fits')
print(fn)
f = fits.open(fn)
w = WCS(f[0].header)
c1 = w.pixel_to_world(1, 2048)
c2 = w.pixel_to_world(4096, 2048)
print(c1)
c1 = SkyCoord(c1, frame='fk5')
c2 = SkyCoord(c2, frame='fk5')
sep = c1.separation(c2)
s = sep.degree
print(s)
计算探测视场大小
最新推荐文章于 2024-02-23 10:37:37 发布