azalt_transformto_radec

懒得起名字22

已于 2022-05-14 18:02:35 修改

阅读量268

点赞数

CC 4.0 BY-SA版权

文章标签： python

于 2022-05-14 17:34:27 首次发布

本文链接：https://blog.csdn.net/qq_40692494/article/details/124771341

这篇博客介绍了如何使用astropy库将天文观测的方位角(AZ)和仰角(EL)从地平坐标系转换为国际天球参考系统(ICRS)。通过读取FITS文件头信息获取观测日期、时间和地理位置，然后利用SkyCoord和EarthLocation类进行坐标转换。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

from astropy.io import fits
import astropy.units as u
from astropy.coordinates import SkyCoord,EarthLocation
from astropy.time import Time

hdulistex = fits.open('/media/zyy/Disk_D/望远镜检测/ali20220506/站点-20220505-694低轨-600ms-230407-可见光/站点-20220505-694低轨-600ms-230407-可见光_FITS/142945.fits/FITS_00000.fits')
hdr = hdulistex[0].header
az = hdr['AZ']
ele = hdr['EL']
date= hdr['DATE-OBS']
tm = hdr['TIME-OBS']
t = date+' '+tm
location = EarthLocation(lat=32.3052787780762 * u.deg, lon=80.0458297729492 * u.deg, height=5375.99992752075 * u.m)
crab = SkyCoord(obstime= Time(t), location=location,az = az*u.deg, alt=ele*u.deg,frame = 'altaz')
coord = crab.transform_to('icrs')
ra = str((coord.ra*u.deg).value)
dec = (coord.dec*u.deg).value
print(ra,dec)
print(type(ra))