我一直试图弄清楚如何将一组AltAz坐标转换为赤道坐标,到目前为止,我所能找到的是如何使用以下方法将赤道转换为AltAz(但不是相反):
c = SkyCoord('22h50m0.19315s', '+24d36m05.6984s', frame='icrs')
loc = EarthLocation(lat = 31.7581*u.deg, lon = -95.6386*u.deg, height = 147*u.m)
time = Time('1991-06-06 12:00:00')
cAltAz = c.transform_to(AltAz(obstime = time, location = loc))
但是现在我想将mpAltAz的方位角旋转一些增量,并找出新点的相应赤道坐标。
即。我想要这样的东西:
newAltAzcoordiantes = SkyCoord(alt = cAltAz.alt.deg, az = cAltAz.az.deg + x*u.deg, obstime = time, frame = 'altaz')
newAltAzcoordiantes.transform_to(ICRS)
我遇到的问题是,我似乎无法在AltAz坐标系中构建SkyCoord对象
我希望这很清楚,这是我第一次在stackoverflow上发帖。
答案 0 :(得分:0)
我对天文学知之甚少,但似乎有很多文档:
对我而言,cAltAz.icrs
有效并返回原始c
。我需要调整一堆内容才能使其在newAltAzcoordiantes
上运行(需要定义x
,停止在已有单位的属性上调用deg
,添加location
):
>>> x = 5 # why not?
>>> newAltAzcoordiantes = SkyCoord(alt = cAltAz.alt, az = cAltAz.az + x*u.deg, obstime = time, frame = 'altaz', location = loc)
>>> newAltAzcoordiantes.icrs
<SkyCoord (ICRS): (ra, dec) in deg
(341.79674062, 24.35770826)>