衛星で観測された物体が月/太陽によって隠されるとき

Question

私は衛星の経度、緯度と高度を持っています。それから私はRa & Declと遠くのオブジェクトを持っています。私は与えられた時間にオブジェクトと月/太陽の間の角距離を知る必要があります。

import ephem 
sun = ephem.Sun() 
scraft = ephem.Observer() 
scraft.lon = lon 
scraft.lat = lat 
scraft.elevation = altit 
scraft.date = time-15019.5 # my time is in MJD 
sun.compute(scraft) 
print float(sun.dec), float(sun.ra)   

これは正しい方法ですか？ 結果を確認する方法はありますか？ （ウェブインタフェースを備えた別のツールで、ほとんどの場合数字を入力して違いがあるかどうかを確認できます）

Answer 1

遠くにあるオブジェクトの位置を表す固定オブジェクトを作成し、衛星から見たときのその位置と太陽の位置との間の分離のために：

f = ephem.FixedBody() 
f._ra = '1:23:45.0' 
f._dec = '6:78:90.0' 
f.compute(scraft.date) 
print ephem.separation(sun, f), 'degrees' 
print float(ephem.separation(sun, f)), 'radians'