Skip to content
Snippets Groups Projects
Normal view Permalink
StationData.java 8.51 KiB
Newer Older
noeljanes's avatar
First set of files
noeljanes committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 
package fr.cs.examples.estimation;

import java.io.BufferedReader;
import java.io.File;
import java.io.FileInputStream;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.PrintStream;
import java.io.UnsupportedEncodingException;
import java.net.URISyntaxException;
import java.text.ParseException;
import java.util.ArrayList;
import java.util.Collections;
import java.util.Comparator;
import java.util.HashMap;
import java.util.List;
import java.util.Locale;
import java.util.Map;
import java.util.NoSuchElementException;
import java.util.SortedSet;
import java.util.TreeSet;
import java.util.regex.Pattern;

import org.hipparchus.exception.LocalizedCoreFormats;
import org.hipparchus.geometry.euclidean.threed.Vector3D;
import org.hipparchus.linear.QRDecomposer;
import org.hipparchus.optim.nonlinear.vector.leastsquares.GaussNewtonOptimizer;
import org.hipparchus.optim.nonlinear.vector.leastsquares.LeastSquaresOptimizer;
import org.hipparchus.optim.nonlinear.vector.leastsquares.LeastSquaresProblem;
import org.hipparchus.optim.nonlinear.vector.leastsquares.LevenbergMarquardtOptimizer;
import org.hipparchus.stat.descriptive.StreamingStatistics;
import org.hipparchus.util.FastMath;
import org.hipparchus.util.Precision;
import org.orekit.attitudes.AttitudeProvider;
import org.orekit.attitudes.BodyCenterPointing;
import org.orekit.attitudes.LofOffset;
import org.orekit.attitudes.NadirPointing;
import org.orekit.attitudes.YawCompensation;
import org.orekit.attitudes.YawSteering;
import org.orekit.bodies.CelestialBody;
import org.orekit.bodies.CelestialBodyFactory;
import org.orekit.bodies.GeodeticPoint;
import org.orekit.bodies.OneAxisEllipsoid;
import org.orekit.data.DataProvidersManager;
import org.orekit.data.DirectoryCrawler;
import org.orekit.errors.OrekitException;
import org.orekit.errors.OrekitMessages;
import org.orekit.estimation.leastsquares.BatchLSEstimator;
import org.orekit.estimation.leastsquares.BatchLSObserver;
import org.orekit.estimation.measurements.AngularAzEl;
import org.orekit.estimation.measurements.EstimatedMeasurement;
import org.orekit.estimation.measurements.EstimationsProvider;
import org.orekit.estimation.measurements.GroundStation;
import org.orekit.estimation.measurements.ObservableSatellite;
import org.orekit.estimation.measurements.ObservedMeasurement;
import org.orekit.estimation.measurements.PV;
import org.orekit.estimation.measurements.Range;
import org.orekit.estimation.measurements.RangeRate;
noeljanes's avatar
Added troposphere corrections to range-rate  
noeljanes committed
59 
import org.orekit.estimation.measurements.modifiers.*;
noeljanes's avatar
First set of files
noeljanes committed
60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 
import org.orekit.forces.PolynomialParametricAcceleration;
import org.orekit.forces.drag.DragForce;
import org.orekit.forces.drag.DragSensitive;
import org.orekit.forces.drag.IsotropicDrag;
import org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel;
import org.orekit.forces.gravity.OceanTides;
import org.orekit.forces.gravity.Relativity;
import org.orekit.forces.gravity.SolidTides;
import org.orekit.forces.gravity.ThirdBodyAttraction;
import org.orekit.forces.gravity.potential.GravityFieldFactory;
import org.orekit.forces.gravity.potential.NormalizedSphericalHarmonicsProvider;
import org.orekit.forces.radiation.IsotropicRadiationSingleCoefficient;
import org.orekit.forces.radiation.RadiationSensitive;
import org.orekit.forces.radiation.SolarRadiationPressure;
import org.orekit.frames.EOPHistory;
import org.orekit.frames.Frame;
import org.orekit.frames.FramesFactory;
import org.orekit.frames.LOFType;
import org.orekit.frames.TopocentricFrame;
import org.orekit.gnss.Frequency;
import org.orekit.gnss.MeasurementType;
import org.orekit.gnss.ObservationData;
import org.orekit.gnss.ObservationDataSet;
import org.orekit.gnss.RinexLoader;
import org.orekit.gnss.SatelliteSystem;
import org.orekit.models.AtmosphericRefractionModel;
import org.orekit.models.earth.EarthITU453AtmosphereRefraction;
import org.orekit.models.earth.atmosphere.Atmosphere;
import org.orekit.models.earth.atmosphere.DTM2000;
import org.orekit.models.earth.atmosphere.data.MarshallSolarActivityFutureEstimation;
import org.orekit.models.earth.displacement.OceanLoading;
import org.orekit.models.earth.displacement.OceanLoadingCoefficientsBLQFactory;
import org.orekit.models.earth.displacement.StationDisplacement;
import org.orekit.models.earth.displacement.TidalDisplacement;
import org.orekit.models.earth.ionosphere.IonosphericModel;
import org.orekit.models.earth.ionosphere.KlobucharIonoCoefficientsLoader;
import org.orekit.models.earth.ionosphere.KlobucharIonoModel;
import org.orekit.models.earth.troposphere.DiscreteTroposphericModel;
import org.orekit.models.earth.troposphere.EstimatedTroposphericModel;
import org.orekit.models.earth.troposphere.GlobalMappingFunctionModel;
import org.orekit.models.earth.troposphere.MappingFunction;
import org.orekit.models.earth.troposphere.NiellMappingFunctionModel;
import org.orekit.models.earth.troposphere.SaastamoinenModel;
import org.orekit.orbits.CartesianOrbit;
import org.orekit.orbits.CircularOrbit;
import org.orekit.orbits.EquinoctialOrbit;
import org.orekit.orbits.KeplerianOrbit;
import org.orekit.orbits.Orbit;
import org.orekit.orbits.PositionAngle;
import org.orekit.propagation.SpacecraftState;
import org.orekit.propagation.analytical.tle.TLE;
import org.orekit.propagation.analytical.tle.TLEPropagator;
import org.orekit.propagation.conversion.DormandPrince853IntegratorBuilder;
import org.orekit.propagation.conversion.NumericalPropagatorBuilder;
import org.orekit.time.AbsoluteDate;
import org.orekit.time.DateComponents;
import org.orekit.time.TimeScalesFactory;
import org.orekit.utils.Constants;
import org.orekit.utils.IERSConventions;
import org.orekit.utils.PVCoordinates;
import org.orekit.utils.ParameterDriver;
import org.orekit.utils.ParameterDriversList;

import fr.cs.examples.KeyValueFileParser;

public class StationData {

    /** Container for stations-related data. */

    /** Ground station. */
    public final  GroundStation station;

    /** Range sigma. */
    public final double rangeSigma;

    /** Range bias (may be null if bias is fixed to zero). */
    public final Bias<Range> rangeBias;

    /** Range rate sigma. */
    public final  double rangeRateSigma;

    /** Range rate bias (may be null if bias is fixed to zero). */
    public final  Bias<RangeRate> rangeRateBias;

    /** Azimuth-elevation sigma. */
    public final double[] azElSigma;

    /** Azimuth-elevation bias (may be null if bias is fixed to zero). */
    public final Bias<AngularAzEl> azELBias;

    /** Elevation refraction correction (may be null). */
    public final AngularRadioRefractionModifier refractionCorrection;

    /** Tropospheric correction (may be null). */
    public final RangeTroposphericDelayModifier rangeTroposphericCorrection;
noeljanes's avatar
Added troposphere corrections to range-rate  
noeljanes committed
156 157 
    public final RangeRateTroposphericDelayModifier rangeRateTroposphericCorrection;
noeljanes's avatar
First set of files
noeljanes committed
158 159 160 161 162 163 164 165 166 167 
    /** Simple constructor.
     * @param station ground station
     * @param rangeSigma range sigma
     * @param rangeBias range bias (may be null if bias is fixed to zero)
     * @param rangeRateSigma range rate sigma
     * @param rangeRateBias range rate bias (may be null if bias is fixed to zero)
     * @param azElSigma azimuth-elevation sigma
     * @param azELBias azimuth-elevation bias (may be null if bias is fixed to zero)
     * @param refractionCorrection refraction correction for elevation (may be null)
     * @param rangeTroposphericCorrection tropospheric correction  for the range (may be null)
noeljanes's avatar
Added troposphere corrections to range-rate  
noeljanes committed
168 
     * @param rangeRateTroposphericCorrection
noeljanes's avatar
First set of files
noeljanes committed
169 170 171 172 173 174 
     */
    public StationData(final GroundStation station,
                       final double rangeSigma, final Bias<Range> rangeBias,
                       final double rangeRateSigma, final Bias<RangeRate> rangeRateBias,
                       final double[] azElSigma, final Bias<AngularAzEl> azELBias,
                       final AngularRadioRefractionModifier refractionCorrection,
noeljanes's avatar
Added troposphere corrections to range-rate  
noeljanes committed
175 176 177 178 179 180 181 182 183 184 185 
                       final RangeTroposphericDelayModifier rangeTroposphericCorrection, RangeRateTroposphericDelayModifier rangeRateTroposphericCorrection) {
        this.station                         = station;
        this.rangeSigma                      = rangeSigma;
        this.rangeBias                       = rangeBias;
        this.rangeRateSigma                  = rangeRateSigma;
        this.rangeRateBias                   = rangeRateBias;
        this.azElSigma                       = azElSigma.clone();
        this.azELBias                        = azELBias;
        this.refractionCorrection            = refractionCorrection;
        this.rangeTroposphericCorrection     = rangeTroposphericCorrection;
        this.rangeRateTroposphericCorrection = rangeRateTroposphericCorrection;
noeljanes's avatar
First set of files
noeljanes committed
186 187 188 
    }

}