객체의 J2000 RA/DEC 좌표를 "관측 위치"로 변환하려고합니다. 즉, 굴절 효과를 포함하는 토 코어 중심의 RA/DEC 좌표로 변환하려고합니다.pyephem으로 굴절 보정 된 RA/DEC 좌표를 계산하는 방법은 무엇입니까?
from math import pi
import ephem
from datetime import datetime
ra = 20.370473492/12. * pi
dec = 40.256674958/180. * pi
tt = datetime(2016, 07, 27, 23, 30, 0)
lowell = ephem.Observer()
lowell.lon = '-111:32.1'
lowell.lat = '35:05.8'
lowell.elevation = 2198
lowell.date = tt
lowell.pressure = 1000
bd = ephem.FixedBody()
bd._ra = ra
bd._dec = dec
bd.compute(lowell)
print "Pressure: ", lowell.pressure
print "Input: ", bd._ra, bd._dec
print "Astrometric Geocentric Position: ", bd.a_ra, bd.a_dec
print "Apparent Geocentric Position ", bd.g_ra, bd.g_dec
print "Apparent Topocentric Position: ", bd.ra, bd.dec
print "Horizontal Position: ", bd.alt, bd.az
lowell.pressure = 0
bd.compute(lowell)
print "Pressure: ", lowell.pressure
print "Input: ", bd._ra, bd._dec
print "Astrometric Geocentric Position: ", bd.a_ra, bd.a_dec
print "Apparent Geocentric Position ", bd.g_ra, bd.g_dec
print "Apparent Topocentric Position: ", bd.ra, bd.dec
print "Horizontal Position: ", bd.alt, bd.az
먼저 주목할만한 점은 다음과 같습니다 : 명백한 지구를 중심으로 겉보기 topocentric 좌표가 다르지 않다 워드 프로세서 (http://rhodesmill.org/pyephem/radec.html#how-the-three-positions-differ) 다음 내가 이런 짓을. 시차와 굴절에 대한 교정이 문서에 명시된대로 이루어지면 차이가 있어야합니다.
두 번째 사항 : 압력을 0으로 설정하면 굴절이 사라지고 무엇인가가 변경되어야합니다. 그러나, 수평 좌표 만 변경됩니다 (약 11 분). 그러나 모든 ra/dec 좌표는 동일하게 유지됩니다.
무엇이 여기에 있습니까?
또는 다른 말로하면 어떻게 든 파이 펫에서 굴절 보정 RA/DEC 좌표를 얻을 수 있습니까?
가 (BTW, 나는 최신 pyepehm 버전을 사용, 3.7.6.0)
Brandon 여러분,이 명확한 설명을 해주셔서 대단히 감사합니다. 나는 아직도 문서를 약간 혼란 스럽지만, 당신의 요지는 완벽하게 유효하다. 내 특별한 경우에는 해결 방법 사용을 의미합니다. 망원경 제어 소프트웨어를 프로그래밍 중이며 RA 12 굴절 굴절에서 축 위치를 설정해야합니다. 그러나 그것은 당신이 제안한대로 작동해야합니다. –
Btw, pyephem, skyfields 등에 대한 작업에 큰 감사를드립니다. 그들은 내 작업에 매우 중요합니다. –