|||
算不同历元的坐标转换(从J2000到当前历元)我以前用以下代码。
star = ephem.FixedBody()
star._ra = ra0
star._dec = dec0
star._epoch = '2000'
star.compute(epoch=ephem.now())
print ephem.degrees(star.ra)/convert,ephem.degrees(star.dec)/convert
一直没发现有问题。直到今天发现J2000转B1950会出问题。
用NED的坐标转换测试了一下,下面的代码似乎是对的。
new = ephem.Equatorial(ra0, dec0, epoch=ephem.J2000)
old = ephem.Equatorial(new, epoch=ephem.now())
print old.ra/convert,old.dec/convert
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 09:19
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社