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.AngularRaDec;
23  import org.orekit.estimation.measurements.EstimationModifier;
24  import org.orekit.estimation.measurements.GroundStation;
25  import org.orekit.estimation.measurements.ObservableSatellite;
26  import org.orekit.frames.Frame;
27  import org.orekit.propagation.SpacecraftState;
28  import org.orekit.propagation.sampling.OrekitStepInterpolator;
29  import org.orekit.time.AbsoluteDate;
30  import org.orekit.utils.ParameterDriver;
31  
32  
33  /** Builder for {@link AngularRaDec} measurements.
34   * @author Luc Maisonobe
35   * @since 9.3
36   */
37  public class AngularRaDecBuilder extends AbstractMeasurementBuilder<AngularRaDec> {
38  
39      /** Ground station from which measurement is performed. */
40      private final GroundStation station;
41  
42      /** Reference frame in which the right ascension - declination angles are given. */
43      private final Frame referenceFrame;
44  
45      /** Satellite related to this builder.
46       * @since 12.0
47       */
48      private final ObservableSatellite satellite;
49  
50      /** Simple constructor.
51       * @param noiseSource noise source, may be null for generating perfect measurements
52       * @param station ground station from which measurement is performed
53       * @param referenceFrame Reference frame in which the right ascension - declination angles are given
54       * @param sigma theoretical standard deviation
55       * @param baseWeight base weight
56       * @param satellite satellite related to this builder
57       */
58      public AngularRaDecBuilder(final CorrelatedRandomVectorGenerator noiseSource,
59                                 final GroundStation station, final Frame referenceFrame,
60                                 final double[] sigma, final double[] baseWeight,
61                                 final ObservableSatellite satellite) {
62          super(noiseSource, sigma, baseWeight, satellite);
63          this.station        = station;
64          this.referenceFrame = referenceFrame;
65          this.satellite      = satellite;
66      }
67  
68      /** {@inheritDoc} */
69      @Override
70      public AngularRaDec build(final AbsoluteDate date, final Map<ObservableSatellite, OrekitStepInterpolator> interpolators) {
71  
72          final double[] sigma                = getTheoreticalStandardDeviation();
73          final double[] baseWeight           = getBaseWeight();
74          final SpacecraftState[] relevant    = new SpacecraftState[] { interpolators.get(satellite).getInterpolatedState(date) };
75  
76          // create a dummy measurement
77          final AngularRaDec dummy = new AngularRaDec(station, referenceFrame, relevant[0].getDate(),
78                                                      new double[] {
79                                                          0.0, 0.0
80                                                      }, sigma, baseWeight, satellite);
81          for (final EstimationModifier<AngularRaDec> modifier : getModifiers()) {
82              dummy.addModifier(modifier);
83          }
84  
85          // set a reference date for parameters missing one
86          for (final ParameterDriver driver : dummy.getParametersDrivers()) {
87              if (driver.getReferenceDate() == null) {
88                  final AbsoluteDate start = getStart();
89                  final AbsoluteDate end   = getEnd();
90                  driver.setReferenceDate(start.durationFrom(end) <= 0 ? start : end);
91              }
92          }
93  
94          // estimate the perfect value of the measurement
95          final double[] angular = dummy.estimateWithoutDerivatives(relevant).getEstimatedValue();
96  
97          // add the noise
98          final double[] noise = getNoise();
99          if (noise != null) {
100             angular[0] += noise[0];
101             angular[1] += noise[1];
102         }
103 
104         // generate measurement
105         final AngularRaDec measurement = new AngularRaDec(station, referenceFrame, relevant[0].getDate(),
106                                                           angular, sigma, baseWeight, satellite);
107         for (final EstimationModifier<AngularRaDec> modifier : getModifiers()) {
108             measurement.addModifier(modifier);
109         }
110         return measurement;
111 
112     }
113 
114 }