GNSSDate.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.time;

import java.io.Serializable;
import java.util.HashMap;
import java.util.List;
import java.util.Map;
import java.util.concurrent.atomic.AtomicReference;

import org.hipparchus.util.FastMath;
import org.orekit.annotation.DefaultDataContext;
import org.orekit.data.DataContext;
import org.orekit.errors.OrekitException;
import org.orekit.errors.OrekitMessages;
import org.orekit.frames.EOPEntry;
import org.orekit.gnss.SatelliteSystem;
import org.orekit.utils.Constants;
import org.orekit.utils.IERSConventions;

/** Container for date in GNSS form.
 * <p> This class can be used to handle {@link SatelliteSystem#GPS GPS},
 * {@link SatelliteSystem#GALILEO Galileo}, {@link SatelliteSystem#BEIDOU BeiDou}
 * and {@link SatelliteSystem#QZSS QZSS} dates. </p>
 * @author Luc Maisonobe (original code)
 * @author Bryan Cazabonne (generalization to all GNSS constellations)
 * @see AbsoluteDate
 */
public class GNSSDate implements Serializable, TimeStamped {

    /** Serializable UID. */
    private static final long serialVersionUID = 20221228L;

    /** Duration of a week in days. */
    private static final int WEEK_D = 7;

    /** Duration of a week in seconds. */
    private static final double WEEK_S = WEEK_D * Constants.JULIAN_DAY;

    /** Reference date for ensuring continuity across GNSS week rollover.
     * @since 9.3.1
     */
    private static AtomicReference<DateComponents> rolloverReference = new AtomicReference<DateComponents>(null);

    /** Week number since the GNSS reference epoch. */
    private final int weekNumber;

    /** Number of seconds since week start. */
    private final double secondsInWeek;

    /** Satellite system to consider. */
    private final SatelliteSystem system;

    /** Corresponding date. */
    private final transient AbsoluteDate date;

    /** Build an instance corresponding to a GNSS date.
     * <p>
     * GNSS dates are provided as a week number starting at
     * the GNSS reference epoch and as a number of seconds
     * since week start.
     * </p>
     * <p>
     * Many interfaces provide week number modulo the constellation week cycle. In order to cope with
     * this, when the week number is smaller than the week cycle, this constructor assumes a modulo operation
     * has been performed and it will fix the week number according to the reference date set up for
     * handling rollover (see {@link #setRolloverReference(DateComponents) setRolloverReference(reference)}).
     * If the week number is equal to the week cycle or larger, it will be used without any correction.
     * </p>
     *
     * <p>This method uses the {@link DataContext#getDefault() default data context}.
     *
     * @param weekNumber week number
     * @param secondsInWeek number of seconds since week start
     * @param system satellite system to consider
     * @see #GNSSDate(int, double, SatelliteSystem, TimeScales)
     * @since 12.0
     */
    @DefaultDataContext
    public GNSSDate(final int weekNumber, final double secondsInWeek, final SatelliteSystem system) {
        this(weekNumber, secondsInWeek, system, DataContext.getDefault().getTimeScales());
    }

    /**
     * Build an instance corresponding to a GNSS date.
     * <p>
     * GNSS dates are provided as a week number starting at the GNSS reference epoch and
     * as a number of seconds since week start.
     * </p>
     * <p>
     * Many interfaces provide week number modulo the constellation week cycle. In order
     * to cope with this, when the week number is smaller than the week cycle, this
     * constructor assumes a modulo operation has been performed and it will fix the week
     * number according to the reference date set up for handling rollover (see {@link
     * #setRolloverReference(DateComponents) setRolloverReference(reference)}). If the
     * week number is equal to the week cycle or larger, it will be used without any
     * correction.
     * </p>
     *
     * @param weekNumber    week number
     * @param secondsInWeek number of seconds since week start
     * @param system        satellite system to consider
     * @param timeScales    the set of time scales. Used to retrieve the appropriate time
     *                      scale for the given {@code system}.
     * @since 12.0
     */
    public GNSSDate(final int weekNumber, final double secondsInWeek,
                    final SatelliteSystem system, final TimeScales timeScales) {

        final int day = (int) FastMath.floor(secondsInWeek / Constants.JULIAN_DAY);
        final double secondsInDay = secondsInWeek - day * Constants.JULIAN_DAY;

        int w = weekNumber;
        DateComponents dc = new DateComponents(getWeekReferenceDateComponents(system), weekNumber * 7 + day);
        final int cycleW = GNSSDateType.getRollOverWeek(system);
        if (weekNumber < cycleW) {

            DateComponents reference = rolloverReference.get();
            if (reference == null) {
                // lazy setting of a default reference, using end of EOP entries
                final UT1Scale       ut1       = timeScales.getUT1(IERSConventions.IERS_2010, true);
                final List<EOPEntry> eop       = ut1.getEOPHistory().getEntries();
                final int            lastMJD   = eop.get(eop.size() - 1).getMjd();
                reference = new DateComponents(DateComponents.MODIFIED_JULIAN_EPOCH, lastMJD);
                rolloverReference.compareAndSet(null, reference);
            }

            // fix GNSS week rollover
            final int cycleD = WEEK_D * cycleW;
            while (dc.getJ2000Day() < reference.getJ2000Day() - cycleD / 2) {
                dc = new DateComponents(dc, cycleD);
                w += cycleW;
            }

        }

        this.weekNumber    = w;
        this.secondsInWeek = secondsInWeek;
        this.system        = system;

        date = new AbsoluteDate(dc, new TimeComponents(secondsInDay), getTimeScale(system, timeScales));

    }

    /**
     * Build an instance corresponding to a GNSS date.
     * <p>
     * GNSS dates are provided as a week number starting at the GNSS reference epoch and
     * as a number of seconds since week start.
     * </p>
     *
     * @param weekNumber    week number
     * @param secondsInWeek number of seconds since week start
     * @param system        satellite system to consider
     * @param reference     reference date for rollover, the generated date will be less
     *                      than one half cycle from this date
     * @param timeScales    the set of time scales. Used to retrieve the appropriate time
     *                      scale for the given {@code system}.
     * @since 12.0
     */
    public GNSSDate(final int weekNumber, final double secondsInWeek,
                    final SatelliteSystem system, final DateComponents reference,
                    final TimeScales timeScales) {

        final int day = (int) FastMath.floor(secondsInWeek / Constants.JULIAN_DAY);
        final double secondsInDay = secondsInWeek - day * Constants.JULIAN_DAY;

        int w = weekNumber;
        DateComponents dc = new DateComponents(getWeekReferenceDateComponents(system), weekNumber * 7 + day);
        final int cycleW = GNSSDateType.getRollOverWeek(system);
        if (weekNumber < cycleW) {

            // fix GNSS week rollover
            final int cycleD = WEEK_D * cycleW;
            while (dc.getJ2000Day() < reference.getJ2000Day() - cycleD / 2) {
                dc = new DateComponents(dc, cycleD);
                w += cycleW;
            }

        }

        this.weekNumber    = w;
        this.secondsInWeek = secondsInWeek;
        this.system        = system;

        date = new AbsoluteDate(dc, new TimeComponents(secondsInDay), getTimeScale(system, timeScales));

    }

    /** Build an instance from an absolute date.
     *
     * <p>This method uses the {@link DataContext#getDefault() default data context}.
     *
     * @param date absolute date to consider
     * @param system satellite system to consider
     * @see #GNSSDate(AbsoluteDate, SatelliteSystem, TimeScales)
     */
    @DefaultDataContext
    public GNSSDate(final AbsoluteDate date, final SatelliteSystem system) {
        this(date, system, DataContext.getDefault().getTimeScales());
    }

    /**
     * Build an instance from an absolute date.
     *
     * @param date       absolute date to consider
     * @param system     satellite system to consider
     * @param timeScales the set of time scales. Used to retrieve the appropriate time
     *                   scale for the given {@code system}.
     * @since 10.1
     */
    public GNSSDate(final AbsoluteDate date,
                    final SatelliteSystem system,
                    final TimeScales timeScales) {

        this.system = system;
        final AbsoluteDate epoch = getWeekReferenceAbsoluteDate(system, timeScales);
        this.weekNumber  = (int) FastMath.floor(date.durationFrom(epoch) / WEEK_S);
        final AbsoluteDate weekStart = new AbsoluteDate(epoch, WEEK_S * weekNumber);
        this.secondsInWeek = date.durationFrom(weekStart);
        this.date          = date;

    }

    /** Set a reference date for ensuring continuity across GNSS week rollover.
     * <p>
     * Instance created using the {@link #GNSSDate(int, double, SatelliteSystem) GNSSDate(weekNumber, secondsInWeek, system)}
     * constructor and with a week number between 0 and the constellation week cycle (cycleW) after this method has been called will
     * fix the week number to ensure they correspond to dates between {@code reference - cycleW / 2 weeks}
     * and {@code reference + cycleW / 2 weeks}.
     * </p>
     * <p>
     * If this method is never called, a default reference date for rollover will be set using
     * the date of the last known EOP entry retrieved from {@link UT1Scale#getEOPHistory() UT1}
     * time scale.
     * </p>
     * @param reference reference date for GNSS week rollover
     * @see #getRolloverReference()
     * @see #GNSSDate(int, double, SatelliteSystem)
     * @since 9.3.1
     */
    public static void setRolloverReference(final DateComponents reference) {
        rolloverReference.set(reference);
    }

    /** Get the reference date ensuring continuity across GNSS week rollover.
     * @return reference reference date for GNSS week rollover
     * @see #setRolloverReference(DateComponents)
     * @see #GNSSDate(int, double, SatelliteSystem)
     * @since 9.3.1
     */
    public static DateComponents getRolloverReference() {
        return rolloverReference.get();
    }

    /** Get the week number since the GNSS reference epoch.
     * <p>
     * The week number returned here has been fixed for GNSS week rollover, i.e.
     * it may be larger than the corresponding week cycle of the constellation.
     * </p>
     * @return week number since since the GNSS reference epoch
     */
    public int getWeekNumber() {
        return weekNumber;
    }

    /** Get the number of milliseconds since week start.
     * @return number of milliseconds since week start
     */
    public double getMilliInWeek() {
        return getSecondsInWeek() * 1000.0;
    }

    /** Get the number of seconds since week start.
     * @return number of seconds since week start
     * @since 12.0
     */
    public double getSecondsInWeek() {
        return secondsInWeek;
    }

    /** {@inheritDoc} */
    @Override
    public AbsoluteDate getDate() {
        return date;
    }

    /** Get the time scale related to the given satellite system.
     * @param satellite satellite system
     * @param timeScales set of time scales.
     * @return the time scale
     */
    private TimeScale getTimeScale(final SatelliteSystem satellite,
                                   final TimeScales timeScales) {
        switch (satellite) {
            case GPS     : return timeScales.getGPS();
            case GALILEO : return timeScales.getGST();
            case QZSS    : return timeScales.getQZSS();
            case BEIDOU  : return timeScales.getBDT();
            case IRNSS   : return timeScales.getIRNSS();
            case SBAS    : return timeScales.getGPS();
            default      : throw new OrekitException(OrekitMessages.INVALID_SATELLITE_SYSTEM, satellite);
        }
    }

    /** Get the reference epoch of the week number for the given satellite system.
     * <p> Returned parameter is an AbsoluteDate. </p>
     * @param satellite satellite system
     * @param timeScales set of time scales.
     * @return the reference epoch
     */
    private AbsoluteDate getWeekReferenceAbsoluteDate(final SatelliteSystem satellite,
                                                      final TimeScales timeScales) {
        switch (satellite) {
            case GPS     : return timeScales.getGpsEpoch();
            case GALILEO : return timeScales.getGalileoEpoch();
            case QZSS    : return timeScales.getQzssEpoch();
            case BEIDOU  : return timeScales.getBeidouEpoch();
            case IRNSS   : return timeScales.getIrnssEpoch();
            case SBAS    : return timeScales.getGpsEpoch();
            default      : throw new OrekitException(OrekitMessages.INVALID_SATELLITE_SYSTEM, satellite);
        }
    }

    /** Get the reference epoch of the week number for the given satellite system.
     * <p> Returned parameter is a DateComponents. </p>
     * @param satellite satellite system
     * @return the reference epoch
     */
    private DateComponents getWeekReferenceDateComponents(final SatelliteSystem satellite) {
        switch (satellite) {
            case GPS     : return DateComponents.GPS_EPOCH;
            case GALILEO : return DateComponents.GALILEO_EPOCH;
            case QZSS    : return DateComponents.QZSS_EPOCH;
            case BEIDOU  : return DateComponents.BEIDOU_EPOCH;
            case IRNSS   : return DateComponents.IRNSS_EPOCH;
            case SBAS    : return DateComponents.GPS_EPOCH;
            default      : throw new OrekitException(OrekitMessages.INVALID_SATELLITE_SYSTEM, satellite);
        }
    }

    /** Replace the instance with a data transfer object for serialization.
     * @return data transfer object that will be serialized
     */
    @DefaultDataContext
    private Object writeReplace() {
        return new DataTransferObject(weekNumber, secondsInWeek, system);
    }

    /** Internal class used only for serialization. */
    @DefaultDataContext
    private static class DataTransferObject implements Serializable {

        /** Serializable UID. */
        private static final long serialVersionUID = 20221228L;

        /** Week number since the GNSS reference epoch. */
        private final int weekNumber;

        /** Number of seconds since week start. */
        private final double secondsInWeek;

        /** Satellite system to consider. */
        private final SatelliteSystem system;

        /** Simple constructor.
         * @param weekNumber week number since the GNSS reference epoch
         * @param secondsInWeek number of seconds since week start
         * @param system satellite system to consider
         */
        DataTransferObject(final int weekNumber, final double secondsInWeek,
                           final SatelliteSystem system) {
            this.weekNumber    = weekNumber;
            this.secondsInWeek = secondsInWeek;
            this.system        = system;
        }

        /** Replace the deserialized data transfer object with a {@link GNSSDate}.
         * @return replacement {@link GNSSDate}
         */
        private Object readResolve() {
            return new GNSSDate(weekNumber, secondsInWeek, system);
        }

    }

    /** Enumerate for GNSS data. */
    private enum GNSSDateType {

        /** GPS. */
        GPS(SatelliteSystem.GPS, 1024),

        /** Galileo. */
        GALILEO(SatelliteSystem.GALILEO, 4096),

        /** QZSS. */
        QZSS(SatelliteSystem.QZSS, 1024),

        /** BeiDou. */
        BEIDOU(SatelliteSystem.BEIDOU, 8192),

        /** IRNSS. */
        IRNSS(SatelliteSystem.IRNSS, 1024),

        /** SBAS. */
        SBAS(SatelliteSystem.SBAS, 1024);

        /** Map for the number of week in one GNSS rollover cycle. */
        private static final Map<SatelliteSystem, Integer> CYCLE_MAP = new HashMap<SatelliteSystem, Integer>();
        static {
            for (final GNSSDateType type : values()) {
                final int             val       = type.getRollOverCycle();
                final SatelliteSystem satellite = type.getSatelliteSystem();
                CYCLE_MAP.put(satellite, val);
            }
        }

        /** Number of week in one rollover cycle. */
        private final int numberOfWeek;

        /** Satellite system. */
        private final SatelliteSystem satelliteSystem;

        /**
         * Build a new instance.
         *
         * @param system satellite system
         * @param rollover number of week in one rollover cycle
         */
        GNSSDateType(final SatelliteSystem system, final int rollover) {
            this.satelliteSystem = system;
            this.numberOfWeek    = rollover;
        }

        /** Get the number of week in one rollover cycle.
         * @return  the number of week in one rollover cycle
         */
        private int getRollOverCycle() {
            return numberOfWeek;
        }

        /** Get the satellite system.
         * @return the satellite system
         */
        private SatelliteSystem getSatelliteSystem() {
            return satelliteSystem;
        }

        /** Get the number of week in one rollover cycle for the given satellite system.
         *
         * @param satellite satellite system
         * @return the number of week in one rollover cycle for the given satellite system
         */
        private static int getRollOverWeek(final SatelliteSystem satellite) {
            return CYCLE_MAP.get(satellite);
        }

    }
}