azalt_transformto_radec

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))

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值