2013-04-16 3 views
1

저는 최근에 pyephem을 다운로드했습니다. 나는 많은 사용자 인터페이스를 좋아한다.Pyephem : 세차 운동 계산이 정확하게 세세하게 적용되지 않았습니까?

그러나 나는 과거의 어떤 시점에서 별의 적위를 추정하려고 할 때 특히 그다지 색다른 결과를 얻고 있습니다. 나는 정확히 http://rhodesmill.org/brandon/projects/pyephem-manual.html#FixedObjects,Precession,andEpochs에 설명 된 Thuban/Polaris 예제와 관련된 코드를 사용했으며 완전히 다른 결과를 얻었습니다. 나는 일을 만들기 위해 스크립트에 추가해야 할 다른 것이 있는지, pyephem에 새로운 버그가 있는지 궁금하다.

내 스크립트는 다음과 같습니다

import ephem 
import ephem.stars 

polaris = ephem.readdb("Polaris,f|M|F7,2:31:48.704,89:15:50.72,2.02,2000") 
polaris.compute()  

thuban = ephem.readdb("Thuban,f|V|A0,14:4:23.3,64:22:33,3.65,2000") 
thuban.compute() 

print 'my result for current polaris dec',polaris.dec 
print 'webpage result     89:15:50.73\n' 

print 'my result for current thuban dec',thuban.dec 
print 'webpage result     64:22:32.99\n' 

polaris.compute(epoch='-2800') 
print 'my result for polaris dec in -2800',polaris.dec 
print 'webpage result      63:33:17.63\n' 

thuban.compute(epoch='-2800') 
print 'my result for thuban dec in -2800',thuban.dec 
print 'webpage result     89:54:34.97\n' 

는 다음과 같은 출력을 생성 실행합니다. 처음 두 비교가 약간 다른 이유를 이해할 수 있습니다. (웹 페이지의 예제는 최근 과거에 생성되었습니다.) 그러나 여러분은 어느 쪽의 별에 대해서도 현재 날짜와 에포크 -2800 사이에서 편차가 변경되지 않았 음을 알 수 있습니다.

my result for current polaris dec 89:19:19.5 
webpage result     89:15:50.73 

my result for current thuban dec 64:18:39.5 
webpage result     64:22:32.99 

my result for polaris dec in -2800 89:19:19.5 
webpage result      63:33:17.63 

my result for thuban dec in -2800 64:18:39.5 
webpage result     89:54:34.97 

답변

0

.dec 속성이 자동으로 조정하는 방법 세차 작품을 보여주는 튜토리얼의 시도를 방해 한 것 몇 가지 방법입니다 : 웹 페이지의 예는 편차는 두 시대 사이의 두 별을 큰 거래를 변경 따라서 튜토리얼의 코드를주의 깊게 살펴보면 Polar와 Thuban의 "astrometric declination"에 대해 항상 묻습니다 (.a_dec). 속성 이름 앞에 a_을 추가하여 위의 코드에 약간의 철자 수정을하면 인쇄 된 답변이 자습서에 인쇄 된 답변과 일치한다는 것을 알게 될 것입니다!