Dcb.java
/* Copyright 2002-2023 CS GROUP
* Licensed to CS GROUP (CS) under one or more
* contributor license agreements. See the NOTICE file distributed with
* this work for additional information regarding copyright ownership.
* CS licenses this file to You under the Apache License, Version 2.0
* (the "License"); you may not use this file except in compliance with
* the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
package org.orekit.files.sinex;
import java.util.HashMap;
import java.util.HashSet;
import org.hipparchus.util.Pair;
import org.orekit.gnss.ObservationType;
import org.orekit.time.AbsoluteDate;
import org.orekit.utils.TimeSpanMap;
/**
* Class to store DCB Solution data parsed in the SinexLoader.
* <p>
* This class is made to handle both station and satellite DCB data.
* Bias values are stored in TimeSpanMaps associated with a given pair
* of observation codes. Those TimeSpanMaps are stored in a Map, which
* associate a pair of observation code (as a HashSet of ObservationType)
* to a TimeSpanMap, encapsulated in a DCBCode object.
* </p>
* @author Louis Aucouturier
* @since 12.0
*/
public class Dcb {
/** Ensemble of observation code pairs available for the satellite. */
private HashSet<Pair<ObservationType, ObservationType>> observationSets;
/**
* Ensemble of DCBCode object, identifiable by observation code pairs,
* each containing the corresponding TimeSpanMap of biases (DCB).
*/
private HashMap<Pair<ObservationType, ObservationType>, DcbCode> dcbCodeMap;
/** Simple constructor. */
public Dcb() {
this.observationSets = new HashSet<>();
this.dcbCodeMap = new HashMap<>();
}
// Class to store the TimeSpanMap per DCB Observation Code set
private static class DcbCode {
/** TimeSpanMap containing the DCB bias values. */
private TimeSpanMap<Double> dcbMap;
/**
* Simple constructor.
*/
DcbCode() {
this.dcbMap = new TimeSpanMap<>(null);
}
/**
* Getter for the TimeSpanMap of DCB values.
*
* @return a time span map containing DCB values, for the observation code pair
* corresponding to this DCBCode object
*/
public TimeSpanMap<Double> getDcbTimeMap() {
return dcbMap;
}
}
/**
* Add the content of a DCB line to the DCBSatellite object.
* <p>
* The method check the presence of a Code pair in a map, and add
* values to the corresponding TimeSpanMap.
* </p>
* @param obs1 String corresponding to the first code used for the DCB computation
* @param obs2 String corresponding to the second code used for the DCB computation
* @param spanBegin Absolute Date corresponding to the beginning of the validity span for this bias value
* @param spanEnd Absolute Date corresponding to the end of the validity span for this bias value
* @param biasValue DCB bias value expressed in S.I. units
*/
public void addDcbLine(final String obs1, final String obs2,
final AbsoluteDate spanBegin, final AbsoluteDate spanEnd,
final double biasValue) {
// Setting a pair of observation type.
final Pair<ObservationType, ObservationType> observationPair = new Pair<>(ObservationType.valueOf(obs1), ObservationType.valueOf(obs2));
// If not present add a new DCBCode to the map, identified by the Observation Pair.
// Then add the bias value and validity period.
if (observationSets.add(observationPair)) {
dcbCodeMap.put(observationPair, new DcbCode());
}
dcbCodeMap.get(observationPair).getDcbTimeMap().addValidBetween(biasValue, spanBegin, spanEnd);
}
/**
* Get the value of the Differential Code Bias for a given observation pair
* and a at a given date.
*
* @param obs1 string corresponding to the first code used for the DCB computation
* @param obs2 string corresponding to the second code used for the DCB computation
* @param date date at which to obtain the DCB
* @return the value of the DCB in S.I. units
*/
public double getDcb(final String obs1, final String obs2, final AbsoluteDate date) {
return getDcb(ObservationType.valueOf(obs1), ObservationType.valueOf(obs2), date);
}
/**
* Get the value of the Differential Code Bias for a given observation pair
* and a at a given date.
*
* @param obs1 first observation type
* @param obs2 second observation type
* @param date date at which to obtain the DCB
* @return the value of the DCB in S.I. units
*/
public double getDcb(final ObservationType obs1, final ObservationType obs2, final AbsoluteDate date) {
return getTimeSpanMap(obs1, obs2).get(date);
}
/**
* Get all available observation code pairs for the satellite.
*
* @return HashSet(HashSet(ObservationType)) Observation code pairs obtained.
*/
public HashSet<Pair<ObservationType, ObservationType>> getAvailableObservationPairs() {
return observationSets;
}
/**
* Get the minimum valid date for a given observation pair.
*
* @param obs1 sString corresponding to the first code used for the DCB computation
* @param obs2 string corresponding to the second code used for the DCB computation
* @return minimum valid date for the observation pair
*/
public AbsoluteDate getMinimumValidDateForObservationPair(final String obs1, final String obs2) {
return getMinimumValidDateForObservationPair(ObservationType.valueOf(obs1), ObservationType.valueOf(obs2));
}
/**
* Get the minimum valid date for a given observation pair.
*
* @param obs1 first observation type
* @param obs2 second observation type
* @return minimum valid date for the observation pair
*/
public AbsoluteDate getMinimumValidDateForObservationPair(final ObservationType obs1, final ObservationType obs2) {
return getTimeSpanMap(obs1, obs2).getFirstTransition().getDate();
}
/**
* Get the maximum valid date for a given observation pair.
*
* @param obs1 string corresponding to the first code used for the DCB computation
* @param obs2 string corresponding to the second code used for the DCB computation
* @return maximum valid date for the observation pair
*/
public AbsoluteDate getMaximumValidDateForObservationPair(final String obs1, final String obs2) {
return getMaximumValidDateForObservationPair(ObservationType.valueOf(obs1), ObservationType.valueOf(obs2));
}
/**
* Get the maximum valid date for a given observation pair.
*
* @param obs1 first observation type
* @param obs2 second observation type
* @return maximum valid date for the observation pair
*/
public AbsoluteDate getMaximumValidDateForObservationPair(final ObservationType obs1, final ObservationType obs2) {
return getTimeSpanMap(obs1, obs2).getLastTransition().getDate();
}
/**
* Return the TimeSpanMap object for a given observation code pair,
* for further operation on the object directly.
*
* @param obs1 first observation type
* @param obs2 second observation type
* @return the time span map for a given observation code pair
*/
private TimeSpanMap<Double> getTimeSpanMap(final ObservationType obs1, final ObservationType obs2) {
return dcbCodeMap.get(new Pair<>(obs1, obs2)).getDcbTimeMap();
}
}