GPS- (lon, lat, timestamp) 및 (방위각, 천정)에서 하늘의 점에 대한 은하계 좌표?

Question

주어진 정보 : 옵저버는 GPS로 주어진 위도 및 위도에 있습니다. 하늘에있는 물체를 방위각 및 천정시간에서 GPS로 표시합니다.

어떻게 얻을 수 있습니까? 은하계 좌표 하늘에 포인트가 있습니까? PyEphem으로이 작업을 수행 할 수 있습니까?

나는이 거의 다 저를 얻을 생각 :

from datetime import datetime 
import ephem 
observer = ephem.Observer() 
observer.pressure = 0 
observer.lon = str(4.95569830927) 
observer.lat = str(52.3545603701) 
observer.elevation = 56.8426188165 
observer.date = ephem.Date(datetime.utcfromtimestamp(1409097608)) 
ra, dec = observer.radec_of(str(azimuth), str(zenith)) 

하지만 어떻게 얻을 수있는 은하이에서 조정합니다.

Answer 1

가

생성 할 수있는 은하 이용한 좌표 빠른 참조로 설명한 사물 변환 좌표

가

http://rhodesmill.org/pyephem/quick.html#coordinate-conversion

다음 스크립트가 임의의 방위 및 천정 각도를 취하고 은하 좌표로 변환하여, 응답이 동의를

: http://lambda.gsfc.nasa.gov/toolbox/tb_coordconv.cfm

from datetime import datetime 
import ephem 

azimuth = '45:00' 
zenith = '84:00' 
ninety = ephem.degrees('90:00') 

observer = ephem.Observer() 
observer.pressure = 0 
observer.lon = str(4.95569830927) 
observer.lat = str(52.3545603701) 
observer.elevation = 56.8426188165 
observer.date = ephem.Date(datetime.utcfromtimestamp(1409097608)) 
ra, dec = observer.radec_of(str(azimuth), ninety - ephem.degrees(str(zenith))) 
# print ra 
# print dec 

e = ephem.Equatorial(ra, dec, epoch='2000') 
g = ephem.Galactic(e) 
print g.lon 
print g.lat 

세 빠르게 메모에서 NASA의 계산기

1. PyEphem은 항상 "고도와 방위각"이 아니라 "천정 거리와 방위각"에 대해 말합니다. 따라서 스크립트는 고도를 생성하기 위해 천정 거리를 90 °에서 빼야합니다.

2. 이러한 결과와 NASA의 작은 오차는 PyEphem의 기반이되는 libastro 라이브러리에 내재되어 있습니다.

3. NASA 계산기를 직접 사용하려면 각도 (!)로 표현 된 RA가 필요합니다. 이는 PyEhem에 문구 print ephem.degrees(ra)을 표시하도록 할 수 있습니다.