In [1]:
%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
import pymap3d
from demtrack import load_gmat_track, load_csv_track, track_analysis, track_comparison, reference_lla, moon_ellipsoid

The tracking files used below are generated with GMAT and crash.script.


In [2]:
# Luna centred ICRF ephemeris from 2019-07-25
track_analysis(load_gmat_track('dslwp_luna_pos_25jul.txt'))


IMPACT POINT: Time = 2019-07-31T14:20:03.310932000, Latitude = 16.6858, Longitude = 159.5218, Altitude = -1561.12
With respect to reference: East = -0.4 m, North = 0.5 m , Up = 0.0 m

In [3]:
# Luna centred ICRF ephemeris from 2019-07-18
track_analysis(load_gmat_track('dslwp_luna_pos_18jul.txt'))


IMPACT POINT: Time = 2019-07-31T14:19:54.004249000, Latitude = 16.6785, Longitude = 159.5196, Altitude = -1654.19
With respect to reference: East = -64.6 m, North = -221.6 m , Up = -93.1 m

In [4]:
# Luna centred ICRF ephemeris from 2019-06-28
track_analysis(load_gmat_track('dslwp_luna_pos_28jun.txt'))


IMPACT POINT: Time = 2019-07-31T14:18:43.744813000, Latitude = 16.6834, Longitude = 159.5387, Altitude = -1647.23
With respect to reference: East = 490.5 m, North = -72.6 m , Up = -86.2 m

Impact location obtained by Cees Bassa with CSV file here


In [5]:
pymap3d.geodetic2enu(16.675, 159.617, 0, reference_lla[0], reference_lla[1], reference_lla[2], ell = moon_ellipsoid)


Out[5]:
(2766.3410883251972, -326.94541006886107, 1558.8876808685357)

In [6]:
# CSV file computed by Wei Mingchuan BG2BHC
track_analysis(load_csv_track('DSLWP-B1-MOON LLA Position.csv'))


IMPACT POINT: Time = 2019-07-31T14:20:15.400000000, Latitude = 16.6755, Longitude = 159.6166, Altitude = -1874.87
With respect to reference: East = 2752.9 m, North = -311.6 m , Up = -316.0 m

Comparison between track obtained in GMAT and Wei's track. For some reason it seems that Wei's track is displaced to the southeast.


In [7]:
track_comparison(load_gmat_track('dslwp_luna_pos_25jul.txt'), load_csv_track('DSLWP-B1-MOON LLA Position.csv'))