OrbMod (Orbit Modeling ) - tools for numerical modeling of orbits. Orbit propagation and fitting to observations. Orbit modeling is implemented within the perturbed two-body problem, in Cartesian, Kustaanheimo-Stiefel (KS) and Sperling-Burdet (SB) variables. The flexible set of perturbations can be used:
- central body gravity as point mass
- central body gravity in spherical harmonics expansion
- third bodies gravity as point mass
- spherical harmonics expansion of one of third bodies gravity field
- General relativity effects
- solar radiation pressure whithin 'cannon ball' model
- athmospheric drag whithin 'cannon ball' model
In this software, SPICE routines and kernels are used for coorinate frames and time scales transformation, ephemeris computation.
There is a GUI (C#) for this software.
Pre-build binaries for OrbMod as well as for UI, auxiliary libraris, all required SPICE kernels and input files are available here: https://github.com/vdmitriev90/OrbMod_bin
Here you able to specify:
- central body
- processing mode (orbit propagation or fitting orbit to observations)
- Variables
- intial conditions, its type, coordinates and time frame
- integraion setting (step, order, tolerance and so on)
Here you able to specify the set of perturbation forces:
- Non-central part of central body's gravitation (via spherical harmonics approach). NB: tidal perturbations currently is disable .
- big planets + Pluto as point masses
- 5 additional bodies, whith IDs according to SPICE nomenclature (example: 3==Earth-Moon baricenter, 301== Moon, 499 == Mars, etc.)
- Non-central part of additional body #5 attraction (via spherical harmonics approach)
- General relativity effects
- Solar radiation pressure whithin 'cannon ball' model
- athmospheric drag whithin 'cannon ball' model
Here you can manage the parameters of fitting orbit to observations:
-
Orbit fiting mode
-
(LS) - Ordinary batch least squares adjustment. It means, that full set of observations will process in each iteration
-
(batch_LS_filer) - adjustment orbit to observations using sequential least squares flter. Please, see here for detaials.
-
(pure_Kalman) - Extended Kalman filter(EKF) for detains see (O. Montenbruckand, E. Gil - 2000)
-
(hybrid_Kalman) - Software will try to estimate state vector and covariance using few first epochs via batch least squares, after that standart EKF algorithm will be used.
-
Outliers rejection - rejection epoch with bad observations
-
Settings of sequential least squares flter
-
Settings of EKF
Here you able to specify the data whichkind of data wil be store in files during processing, and rate of output streams
Here you able to:
Load astrometric observation from text file (in MPC format) to intrnal storage (item "Load observations")
View and manage amount of observations to be proceesed. You can change time frames and include/exclude all observation from specific observatory. On pictuer below are shown example for obsevations of Apophis asteroid.
Visualize result of orbit propagation or fiiting (Plot graphs) There are two examples of such graphs.
Graph of different source of accelerations for asteroid Apophis from 2005 to 2015
Graph of post fit residuals (right ascension) for observations of Apophis asteroid from 2005 to 2015, obtained in this software