PyEphemで東方座標（方位角、高度）を赤道座標（RA、Dec）に変換します

Question

これはかなり基本的な作業ですが、私はそれを並べ替えることはできません。

PyEphemドキュメント：

http://rhodesmill.org/pyephem/radec.html

は.altと.az属性における仰角と方位角と、BodyとObserverオブジェクトから明らかTopocentric位置に、他の方法の周りの変換を実行する方法について説明します。

しかし、私は代わりに、高度と方位角から始め、RAと12月を取得する必要がありますか？

az = 3.30084818 #rad 
el = 0.94610742 #rad 
lat = 34.64 #deg 
lon = -103.7 #deg 
alt = 35800.26 #m 
ut = 2455822.20000367 #julian date 

ありがとう：たとえば

は、ここで私は赤道基準フレームでRAと12月を取得したいのですがそのための座標の1セットです！

Answer 1

ここには2つの微妙な点があります。まず、「高度」と「高度」を使用して、PyEphemライブラリでの2つの用語の意味を逆にすることになりました。空の場所を「標高/方位」の代わりに「標高/ /方位角 "の位置にあることを示しています.2番目に、PyEphemは日付を簡単に変換する方法を忘れているようです。 Julian 〜独自の形式です。 。

： ephemの名前が何のためにあるのかを考え出すことにより、他の方向に行くために自分自身を仕事のビット

を念頭に置いたもの規定に、私はこのスクリプトは、あなたの質問に答えるかもしれないと思います

import ephem 

az = 3.30084818 #rad 
el = 0.94610742 #rad 
lat = 34.64 #deg 
lon = -103.7 #deg 
alt = 35800.26 #m 
ut = 2455822.20000367 #julian date 

# Which Julian Date does Ephem start its own count at? 
J0 = ephem.julian_date(0) 

observer = ephem.Observer() 
observer.lon = str(lon) # str() forces deg -> rad conversion 
observer.lat = str(lat) # deg -> rad 
observer.elevation = alt 
observer.date = ut - J0 

print observer.date 
print observer.radec_of(az, el) 

この回答は、この特定の観測に対して正しいかどうかを確認しますか？

2011/9/17 16:48:00 
(9:16:24.95, -0:45:56.8) 

は、それは、この特定の観察のための物理的な意味があるなら、私に教えてください、または数字の1がまだここに間違っているとあれば微調整する必要がありますここでは私のためにどのようなスクリプトプリントです！