Поиск геоцентрической долготы и широты астероида

Я начал использовать программное обеспечение, которое идет в комплекте с книгой Джин Миус «Астрономические алгоритмы». Сейчас я использую обновленный код из naughter(dot)com/aa.html.

Я могу рассчитать гелиоцентрическую долготу и широту астероида, используя элементы орбиты. Ниже показан астероид Юнона. Я не могу найти исходный код или инструкции для преобразования координат астероида, которые я вычислил ниже, в геоцентрическую долготу и широту. Может ли кто-нибудь порекомендовать метод определения геоцентрической долготы и широты астероида? Спасибо.

double Longitude = 0.0; //West Greenwich  
double Latitude  = 51.4826; //North Greenwich  
double TZOffSet = 0.0;


//Asteroid Juno Elements;   Using epoch J2000
CAAEllipticalObjectElements elements;  

//Calculation Julian Date  2017, 9, 25 
elements.JDEquinox = JD;   

elements.T = (JD - 2451545.0) / 36525.0;

//semi-major axis of the orbit
elements.a = 2.6700912;

//eccentricity of the orbit
elements.e = 0.25498122;

//inclination of the plane of the ecliptic
elements.i = 12.98211;

//Longitude of ascending node
elements.omega = 169.91138;

//longitude of perihelion
elements.w = 248.10807; 

//-------------------------------------------------------------------------
CAAEllipticalObjectDetails details = CAAElliptical::Calculate(JD, elements, true);

1) Heliocentric Ecliptic Longitude  = 215.6670202
2) Heliocentric Ecliptic Latitude = 9.377841131
3) True Geocentric RA = 13.90415220
4) True Geocentric Declination = -3.802673727
5) True Geocentric Distance = 4.136531368
6) True Geocentric Light Time = 0.0238906126
7) Astrometric Geocentric RA = 13.90397025
8) Astrometric Geocentric Declination = -3.80217202
9) Astrometric Geocentric Distance = 4.13653445
10) Astrometric Geocentric Light Time = 0.023890630
11) Elongation = 26.74464206
12) Phase Angle  = 7.928218201
//-------------------------------------------------------------------------

Ответы (1)

Это может быть достигнуто с помощью следующих шагов:

  1. вычислить положение Земли в той же гелиоцентрической системе отсчета,
  2. выразить оба положения (Земли и астероида) в гелиоцентрических прямоугольных координатах,
  3. изменить начало кадра, чтобы получить местоположение астероида в геоцентрических прямоугольных координатах,
  4. преобразовать положение астероида в геоцентрическую широту и долготу.
Спасибо. Когда я получу эту настройку через несколько дней, я опубликую исходный код. Я надеюсь, что это работает.