Evaluates TLE mean orbital elements at epoch
[mean_oe, cjd] = CL_tle_evalAvgElem(tle, type_oe [, frame, ut1_tref, tt_tref])
Evaluates the mean orbital elements at epoch numerically. The processing consists in removing short period terms from the osculating elements (in the specified frame) over a short period of time (a few orbits) The mean elements obtained should be comparable with the mean elements defined by the "eckhech" or "lydlp" models.
The orbital elements can be returned in various forms: Keplerian elements, position and velocity, etc...
In a standard situation, the orbital elements are converted to the specified reference frame. The epoch (internally in UTC) is converted to the TREF time scale. Note that if the frame is not a true equator frame, the transformation may not be as accurate as should be as the frame transformation arguments that can be provided are limited to ut1_tref and tt_tref.
It is also possible to give frame a special value: "native". It means that the output reference frame AND time scale are those natively used for TLEs: TEME and UTC respectively. No conversion is performed.
TLE structure (size N)
(string) Type of output orbital elements: "kep", "cir", "cireq", "equin" or "pv" (1x1)
(string, optional) Output frame, see above for details. Default is "ECI".
(optional) UT1-TREF [seconds]. Default is %CL_UT1_TREF (1xN or 1x1)
(optional) TT-TREF [seconds]. Default is %CL_TT_TREF (1xN or 1x1)
Mean orbital elements (of type given by "type_oe") relative to the specified frame [m, m/s, rad] (6xN)
TLE epoch (modified julian date from 1950.0, TREF time scale except if frame == "native") (1xN)
CNES - DCT/SB/MS
fnames = CL_path("tle_examples.txt", CL_home(), "all"); tle = CL_tle_load(fnames); // all "SPOT" satellites I = grep(tle.desc, "SPOT"); tle.desc(I) // cjd in TREF, mean elements relative to ECI [mean_cir, cjd] = CL_tle_evalAvgElem(tle(I), "cir") |