1   /* Copyright 2002-2024 CS GROUP
2    * Licensed to CS GROUP (CS) under one or more
3    * contributor license agreements.  See the NOTICE file distributed with
4    * this work for additional information regarding copyright ownership.
5    * CS licenses this file to You under the Apache License, Version 2.0
6    * (the "License"); you may not use this file except in compliance with
7    * the License.  You may obtain a copy of the License at
8    *
9    *   http://www.apache.org/licenses/LICENSE-2.0
10   *
11   * Unless required by applicable law or agreed to in writing, software
12   * distributed under the License is distributed on an "AS IS" BASIS,
13   * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14   * See the License for the specific language governing permissions and
15   * limitations under the License.
16   */
17  package org.orekit.estimation.measurements.generation;
18  
19  import java.util.Map;
20  
21  import org.hipparchus.random.CorrelatedRandomVectorGenerator;
22  import org.orekit.estimation.measurements.EstimationModifier;
23  import org.orekit.estimation.measurements.GroundStation;
24  import org.orekit.estimation.measurements.ObservableSatellite;
25  import org.orekit.estimation.measurements.TurnAroundRange;
26  import org.orekit.propagation.SpacecraftState;
27  import org.orekit.propagation.sampling.OrekitStepInterpolator;
28  import org.orekit.time.AbsoluteDate;
29  import org.orekit.utils.ParameterDriver;
30  
31  
32  /** Builder for {@link TurnAroundRange} measurements.
33   * @author Luc Maisonobe
34   * @since 9.3
35   */
36  public class TurnAroundRangeBuilder extends AbstractMeasurementBuilder<TurnAroundRange> {
37  
38      /** Primary ground station from which measurement is performed. */
39      private final GroundStation primaryStation;
40  
41      /** Secondary ground station reflecting the signal. */
42      private final GroundStation secondaryStation;
43  
44      /** Satellite related to this builder.
45       * @since 12.0
46       */
47      private final ObservableSatellite satellite;
48  
49      /** Simple constructor.
50       * @param noiseSource noise source, may be null for generating perfect measurements
51       * @param primaryStation ground station from which measurement is performed
52       * @param secondaryStation ground station reflecting the signal
53       * @param sigma theoretical standard deviation
54       * @param baseWeight base weight
55       * @param satellite satellite related to this builder
56       */
57      public TurnAroundRangeBuilder(final CorrelatedRandomVectorGenerator noiseSource,
58                                    final GroundStation primaryStation, final GroundStation secondaryStation,
59                                    final double sigma, final double baseWeight,
60                                    final ObservableSatellite satellite) {
61          super(noiseSource, sigma, baseWeight, satellite);
62          this.primaryStation   = primaryStation;
63          this.secondaryStation = secondaryStation;
64          this.satellite        = satellite;
65      }
66  
67      /** {@inheritDoc} */
68      @Override
69      public TurnAroundRange build(final AbsoluteDate date, final Map<ObservableSatellite, OrekitStepInterpolator> interpolators) {
70  
71          final double sigma                  = getTheoreticalStandardDeviation()[0];
72          final double baseWeight             = getBaseWeight()[0];
73          final SpacecraftState[] relevant    = new SpacecraftState[] { interpolators.get(satellite).getInterpolatedState(date) };
74  
75          // create a dummy measurement
76          final TurnAroundRange dummy = new TurnAroundRange(primaryStation, secondaryStation, relevant[0].getDate(),
77                                                            Double.NaN, sigma, baseWeight, satellite);
78          for (final EstimationModifier<TurnAroundRange> modifier : getModifiers()) {
79              dummy.addModifier(modifier);
80          }
81  
82          // set a reference date for parameters missing one
83          for (final ParameterDriver driver : dummy.getParametersDrivers()) {
84              if (driver.getReferenceDate() == null) {
85                  final AbsoluteDate start = getStart();
86                  final AbsoluteDate end   = getEnd();
87                  driver.setReferenceDate(start.durationFrom(end) <= 0 ? start : end);
88              }
89          }
90  
91          // estimate the perfect value of the measurement
92          double range = dummy.estimateWithoutDerivatives(relevant).getEstimatedValue()[0];
93  
94          // add the noise
95          final double[] noise = getNoise();
96          if (noise != null) {
97              range += noise[0];
98          }
99  
100         // generate measurement
101         final TurnAroundRange measurement = new TurnAroundRange(primaryStation, secondaryStation, relevant[0].getDate(),
102                                                                 range, sigma, baseWeight, satellite);
103         for (final EstimationModifier<TurnAroundRange> modifier : getModifiers()) {
104             measurement.addModifier(modifier);
105         }
106         return measurement;
107 
108     }
109 
110 }