지구 궤도 위치 및 회전 계산

Question

주어진 날짜와 시간에 대한 지구의 위치 (태양에 상대적인)와 축의 회전을 계산하고 싶습니다. 태양이 0,0,0 좌표에 고정되어 있다고 가정하는 것은 괜찮습니다. 예를 들어 문스톤의 중력에 의한 매우 작은 편향도 무시할 수 있습니다. 어느 정도 정도 정확하면 충분합니다.

이 작업을 수행하는 데 도움이되는 라이브러리/소스/데이터가 있습니까?

Answer 1

aa-56 코드 인 which can be downloaded from here에는 태양의 천문력이 포함되어있어 사용자의 필요를 충족시킬 수 있습니다. 고정밀 작업의 경우 JPL의 DE421과 같은보다 정확한 것을 원하지만 불리하게 커다란 계수 테이블이 있습니다. 정확도가 1도에 만족하면 극단적 인 과잉입니다.

주어진 시간에 지구의 자전은 그리니치 항성시에 의해 주어진다.

T = (JD - 2451545.0)/36525 

theta0 = 280.46061837 + 360.98564736629*(JD-2451545.0) + 
      0.000387933*T*T - T*T*T/38710000.0 

:
장 메에 우스 ''천문 알고리즘 "(! 이러한 종류의 계산 위해이 좋은 참조) 율리우스 날짜 JD 환산 theta0위한 화학식 (각도 누적 회전 각도) 을 준다

theta0 = 0 degrees mod 360은 그리니치 자오선이 천체 좌표에서 오른쪽 자리 0시에 정렬되는 순간을 나타냅니다.

Answer 2

예. 필요한 것은 Ephemeris 패키지입니다.

JPL has an online Ephemeris service 당신을 위해 계산을 할 것이므로, 당신이 웹을 사용할 수 있다면 당신은 그것을 공격 할 수 있습니다.

나는 free Ephermeris package도 찾았지만 JPL 사이트를 방문한 것으로 보입니다. 그러나 JPL의 데이터베이스를 다운로드하고 거기에서 작동하도록하는 링크가 있습니다. 따라서 오프라인으로 작업하고 JPL에서 수동으로 데이터베이스를 업데이트하는 것을 신경 쓰지 않으려면 옵션 일 수 있습니다.

Answer 3

ephem 라이브러리를 사용하여 파이썬에서

:

>>> import ephem 
>>> sun = ephem.Sun() 
>>> sun.compute(ephem.now()) 
>>> sun.hlong, sun.hlat, sun.earth_distance 
(69:41:32.6, 0:00:00.2, 0.98602390289306641) 

ephem

가 몸으로 지구의 편리한 표현을 제공하지 않지만, sun.hlong 및 sun.hlat은 지구의 태양 중심의 위도와 경도를 제공합니다. 이것은 더 잘 문서화 될 수 있습니다.

지구 회전의 경우 여기에서 원하는 값을 말할 수 있습니다. (일반적으로 우리는 하루 중 시간을 사용하지만, 나는 그것을 잡아내는 것이 일반적으로 분명하다고 생각합니다.)