Encontrar la longitud y latitud geocéntrica de un asteroide

Empecé a usar el software que viene con el libro Algoritmos astronómicos de Jean Meeus. Ahora estoy usando el código actualizado de naughter(dot)com/aa.html.

Puedo calcular la longitud heliocéntrica y la latitud de un asteroide usando los elementos orbitales. A continuación se muestra el asteroide Juno. No puedo encontrar ningún código fuente o instrucciones para convertir las coordenadas de asteroides que calculé a continuación a longitud y latitud geocéntrica. ¿Alguien puede recomendar un método para encontrar la longitud y latitud geocéntrica de un asteroide? Gracias.

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
//-------------------------------------------------------------------------

Respuestas (1)

Esto se puede lograr a través de los siguientes pasos:

  1. calcular la posición de la Tierra en el mismo marco heliocéntrico,
  2. expresar ambas posiciones (Tierra y asteroide) en coordenadas rectangulares heliocéntricas,
  3. cambiar el origen del marco para obtener la ubicación del asteroide en coordenadas rectangulares geocéntricas,
  4. convertir la posición del asteroide en latitud y longitud geocéntrica.
Gracias. Cuando tenga esa configuración en unos días, publicaré el código fuente. Espero que funcione.