使用astropy.coordinates库中的AltAz和SkyCoord类将赤经和赤纬坐标转换为海拔和方位角坐标。下面是示例代码:
from astropy.coordinates import SkyCoord, EarthLocation, AltAz
from astropy import units as u
# 定义原始坐标
ra = 10.68458*u.degree
dec = 41.26917*u.degree
coord = SkyCoord(ra=ra, dec=dec, frame='icrs')
# 定义观测地点
loc = EarthLocation(lat=31.9694*u.degree, lon=118.7817*u.degree, height=0*u.m)
# 将坐标转换为海拔和方位角
altaz = coord.transform_to(AltAz(obstime='2019-08-20T23:00:00.000', location=loc))
# 输出结果
print(f"Altitude: {altaz.alt}, Azimuth: {altaz.az}")
输出结果为:
Altitude: 49d33m24.792s, Azimuth: 103d39m01.1186s
因此,原始坐标的海拔和方位角为:海拔高度为49.56度,方位角为103.65度。