Can I retrieve longitude of Lunar Reconnaissance Orbiter around the Moon by POST method (URL)?

Question

I see that I can use an URL like this to retrieve position data for LRO around the Moon:

https://ssd.jpl.nasa.gov/horizons_batch.cgi?batch=1&MAKE_EPHEM=%27YES%27&TABLE_TYPE=%27VECTORS%27&OUT_UNITS=%27KM-S%27&REF_SYSTEM=%27J2000%27&VEC_LABELS=%27YES%27&CSV_FORMAT=%27YES%27&OBJ_DATA=%27YES%27&VEC_TABLE=%273%27&STEP_SIZE=%2710m%27&START_TIME=%272019-10-9%2000:00%27&STOP_TIME=%272019-10-10%2001:15%27&COMMAND=%27-85%27&CENTER=%27@301%27

But what I need is LRO longitude, to calculate when it will pass over a specific position, so I tried using Horizon to retrieve also position of a static point on Moon surface to calculate angular distance, but it dos not work, if it's static it has no ephemeris.

Maybe it could be possible using this page, but I was not able to figure out if an how:

https://wgc.jpl.nasa.gov:8443/webgeocalc/

For example this?

https://wgc.jpl.nasa.gov:8443/webgeocalc/#SurfaceInterceptFinder

But in case it works, which would be the corresponding query URL?

Answer 1

I had to set up a whole page to experiment with Horizons, and after weeks of testing I ended with this url, which you can see here exploded for better readability:

https://ssd.jpl.nasa.gov/horizons_batch.cgi?batch=1
COMMAND='-85'            <---- Target: LRO orbiter
CENTER='coord@301'       <---- Observer: Specify coordinates of point on surface with SITE_COORD
SITE_COORD='10,20,1737'  <---- Site location: Longitude, Latitude, Distance 
                               from center of target body.
TABLE_TYPE='VECTORS'     <---- Ask vectors table containig quantity "RG", which is 
                               target/observer RanGe (distance)
VECT_TABLE='6'           <---- Type of vectors table
START_TIME='2019/11/28 17:16'
STOP_TIME='2019/11/29 17:16'
STEP_SIZE='1h'    
                 ---- Other parameters (useless?) -----
OBJ_DATA='yes'
MAKE_EPHEM='yes'
REF_PLANE='ECLIPTIC'
COORD_TYPE='GEODETIC'
REF_SYSTEM='J2000'
OUT_UNITS='KM-S'
VECT_CORR='NONE'
CAL_FORMAT='CAL'
ANG_FORMAT='HMS'
APPARENT='AIRLESS'
TIME_TYPE='UTC'
TIME_DIGITS='MINUTES'
RANGE_UNITS='AU'
SUPPRESS_RANGE_RATE='no'
SKIP_DAYLT='no'
EXTRA_PREC='yes'
CSV_FORMAT='yes'
VEC_LABELS='yes'
ELM_LABELS='yes'
TP_TYPE='ABSOLUTE'
R_T_S_ONLY='NO'
CA_TABLE_TYPE='STANDARD'

Result:

*******************************************************************************
            JDTDB,            Calendar Date (TDB),                     LT,                     RG,                     RR,
**************************************************************************************************************************
$$SOE
2458816.219444444, A.D. 2019-Nov-28 17:16:00.0000,  7.629721177258104E-03,  2.287332865584860E+03,  8.814594260485629E-01,

Use this command to get a list of all supported bodies and spacecrafts (COMMAND = '*')

Use this command to get a list of all supported surface sites for Moon (CENTER='*@301')

Answer 2

The steps are summarised as follows

curl $url 2&>1 | awk -F ',' '{if(NF==12) print \$1, \$3, \$4, \$5}' | tail -n +2

1. Redirect it to a file. This will get you the Time, Px_J2000, Py_J2000, Pz_J2000.

2. Now, you need to get the J2000_To_IAU Matrix somehow and then get Time, Px_IAU, Py_IAU, Pz_IAU.

3. Get the sidereal angle at some epoch. Convert it to MCMF (Moon Centred Moon Fixed).

4. Now, you most likely need the Latitude and Longitude in ME frame. So you need PA_TO_ME constant matrix and get Px_ME, Py_ME, Pz_ME. Now getting longitude is as simple as $atan2(y,x)$

For reference frame definition: Standardized Lunar Coordinate Systems