使用PyEphem將座標轉換爲赤道座標（方位角，高程）至赤道座標（RA，Dec）

Question

這看起來相當基本，但我無法整理出來。

的PyEphem文檔：

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

介紹瞭如何圍繞執行轉換的另一種方式，從Body和Observer對象視站心位置，與在.alt和.az屬性仰角和方位角。

但是，我應該如何從高程和方位角開始並獲得RA和Dec？

例如，下面是一組座標對此我想獲得RA和DEC赤道參照系：

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 

謝謝！

Answer 1

這裏有兩個微妙之處。首先，您碰巧使用「高程」和「高度」來表示PyEphem圖書館中兩個術語的含義相反 - 因此您要在天空中將「高程/方位角」位置而不是「高度」/azimuth「的位置;第二，似乎PyEphem已經忘記提供一種簡單的方法來將從 Julian 轉換爲自己的格式，雖然有一個函數julian_date()會走向另一個方向，但我們必須做。有些工作自己去搞清楚什麼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) 

讓我知道，如果讓身體感覺這個特殊的觀察，或者如果其中一個號碼是錯在這裏，仍然需要進行調整！