I have state vectors describing the satellite’s orbit under the influence of moon, sun, radiation, and the spherical Earth with zonal and tesseral terms. I want to generate a time series of mean orbital elements from these state vectors so that they resemble TLE data in some way. The purpose is to generate data for model training.
I’ve read many similar answers on this forum (this this, this), but they mostly deal with converting from TLE to osculating elements.
For converting from osculating (state vectors) to TLE, the general consensus seems to be using SGP4 to iteratively guess TLE elements that return the closest state vector. However, I believe this solution is not computationally feasible for my application.
I’ve also looked at the Brouver and Kozai methods described in Vallado and this article that discusses osculating to mean conversion, but the theory behind it seems to be quite complex and I’m still quite confused about these methods and not sure if I can even use them for my application.
My understanding is that TLE are essentially the average of observed state vectors. So, when given TLE is fed to SGP4, the predicted trajectory fits the observations as best as possible. However, based on this and the fact that I have essentially very closely spaced observations, why isn't it possible to numerically integrate (average) the elements over several periods to produce mean elements? I think this averaging should remove the periodic oscillations from perturbations, but I haven’t found any resources about it. Also unsure about the time window for the averaging(one period, days, etc.)
The question is, what’s the best way to go from state vectors to mean elements in my case?