算不同历元的坐标转换（从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