2012-01-22 3 views
4

이 작업은 상당히 기본적인 것처럼 보이지만이를 해결할 수는 없습니다.PyEphem을 사용하여 토코 센틱 좌표 (방위각, 고도)를 적도 좌표 (RA, Dec)로 변환

PyEphem 문서 :

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

.alt.az 속성의 고도와 방위각과 더불어, Body 및 피상 Topocentric 위치에 Observer 객체에서, 주위의 전환이 다른 방법을 수행하는 방법에 대해 설명합니다 .

그러나 어떻게하면 고도와 방위각에서 시작하여 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와 십이 좀하고 싶습니다되는 좌표의 한 세트입니다!

답변

2

여기에는 두 가지 미묘한 차이가 있습니다. 먼저 PyEphem 라이브러리에서 두 용어가 의미하는 것과 반대되는 것을 의미하는 "고도"와 "고도"를 사용했습니다. 그래서 하늘의 지점을 고도 대신 "고도/방위각"위치라고 부릅니다/방위각 (azimuth) "위치에 있습니다. 둘째로, PyEphem은 날짜를 쉽게 변환 할 수있는 방법을 잊어 버렸습니다. 줄리안 ~ 자체 형식입니다 .이 다른 방향으로가는 함수가 있지만, .

: 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) 

가이 특정 관찰 물리적 의미가 알려줘, 또는 번호 중 하나는 아직 여기 잘못하고 경우에 쥐게 될 필요가 : 여기에 나를 위해 무엇을 스크립트 인쇄입니다!

+0

감사합니다. Brandon, 정말 잘 작동합니다. 여기에 JPL 수평선 결과와 비교 한 예가 있습니다. https://gist.github.com/1672906 –

+0

답변이 왜 더 큰 정확도에 동의하지 않는지 궁금합니다. 나는 그걸 조사해야 할 수도 있습니다! JPL 웹 사이트의 양식에서이 두 숫자를 어디에서 얻었습니까? –

+0

예, JPL Horizon 웹 사이트에서 Horizon에 다른 수정 사항이 포함되어 있는지 궁금합니다. –