EOPC04FilesLoader.java

  1. /* Copyright 2002-2020 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.frames;

  19. import java.io.BufferedReader;
  20. import java.io.IOException;
  21. import java.io.InputStream;
  22. import java.io.InputStreamReader;
  23. import java.nio.charset.StandardCharsets;
  24. import java.util.ArrayList;
  25. import java.util.Collection;
  26. import java.util.List;
  27. import java.util.SortedSet;
  28. import java.util.function.Supplier;
  29. import java.util.regex.Matcher;
  30. import java.util.regex.Pattern;

  32. import org.orekit.data.DataProvidersManager;
  33. import org.orekit.errors.OrekitException;
  34. import org.orekit.errors.OrekitMessages;
  35. import org.orekit.time.AbsoluteDate;
  36. import org.orekit.time.DateComponents;
  37. import org.orekit.time.TimeScale;
  38. import org.orekit.utils.Constants;
  39. import org.orekit.utils.IERSConventions;
  40. import org.orekit.utils.IERSConventions.NutationCorrectionConverter;

  42. /** Loader for EOP xx C04 files.
  43.  * <p>EOP xx C04 files contain {@link EOPEntry
  44.  * Earth Orientation Parameters} consistent with ITRF20xx for one year periods, with various
  45.  * xx (05, 08, 14) depending on the data source.</p>
  46.  * <p>The EOP xx C04 files retrieved from the ftp site
  47.  * <a href="ftp://ftp.iers.org/products/eop/long-term/">ftp://ftp.iers.org/products/eop/long-term/</a>
  48.  * are recognized thanks to their base names, which must match one of the patterns
  49.  * {@code eopc04_##_IAU2000.##} or {@code eopc04_##.##} (or the same ending with <code>.gz</code> for
  50.  * gzip-compressed files) where # stands for a digit character.</p>
  51.  * <p>
  52.  * Beware that files retrieved from the web site
  53.  * <a href="http://hpiers.obspm.fr/eoppc/eop/">http://hpiers.obspm.fr/eoppc/eop/</a> do
  54.  * <em>not</em> follow the same naming convention. They lack the _05, _08 or _14
  55.  * marker (and for EOP 14 C04 even the directory does not have this marker anymore...),
  56.  * so all the files in this site have the same name and can easily be confused. This is the reason
  57.  * why the default regular expression in {@link FramesFactory} does not recognize them and
  58.  * enforces the year marker presence.
  59.  * </p>
  60.  * <p>Between 2002 and 2007, another series of Earth Orientation Parameters was
  61.  * in use: EOPC04 (without the _##). These parameters were consistent with the
  62.  * previous ITRS realization: ITRF2000.</p>
  63.  * <p>Since 2008, several series of Earth Orientation Parameters have been available:
  64.  * EOP 05 C04 consistent with ITRF 2005, EOP 08 C04 consistent with ITRF 2008, and
  65.  * EOP 14 C04 consistent with ITRF 2014. As of mid 2017, only the EOP 08 C04 and
  66.  * EOP 14 C04 are updated for current years.</p>
  67.  * <p>Files are available at IERS FTP site: <a
  68.  * href="ftp://ftp.iers.org/products/eop/long-term/">ftp://ftp.iers.org/products/eop/long-term/</a>.</p>
  69.  * <p>
  70.  * This class is immutable and hence thread-safe
  71.  * </p>
  72.  * @author Luc Maisonobe
  73.  */
  74. class EOPC04FilesLoader extends AbstractEopLoader implements EOPHistoryLoader {

  76.     /** Pattern to match the columns header. */
  77.     private static final Pattern COLUMNS_HEADER_PATTERN;

  79.     /** Pattern for data lines. */
  80.     private static final Pattern DATA_LINE_PATTERN;

  82.     /** Year field. */
  83.     private static final int YEAR_FIELD;

  85.     /** Month field. */
  86.     private static final int MONTH_FIELD;

  88.     /** Day field. */
  89.     private static final int DAY_FIELD;

  91.     /** MJD field. */
  92.     private static final int MJD_FIELD;

  94.     /** X component of pole motion field. */
  95.     private static final int POLE_X_FIELD;

  97.     /** Y component of pole motion field. */
  98.     private static final int POLE_Y_FIELD;

  100.     /** UT1-UTC field. */
  101.     private static final int UT1_UTC_FIELD;

  103.     /** LoD field. */
  104.     private static final int LOD_FIELD;

  106.     /** Correction for nutation first field (either dX or dPsi). */
  107.     private static final int NUT_0_FIELD;

  109.     /** Correction for nutation second field (either dY or dEps). */
  110.     private static final int NUT_1_FIELD;

  112.     static {
  113.         // Header have either the following form:
  114.         //       Date      MJD      x          y        UT1-UTC       LOD         dPsi      dEps       x Err     y Err   UT1-UTC Err  LOD Err    dPsi Err   dEpsilon Err
  115.         //                          "          "           s           s            "         "        "          "          s           s            "         "
  116.         //      (0h UTC)
  117.         // or the following form:
  118.         //       Date      MJD      x          y        UT1-UTC       LOD         dX        dY        x Err     y Err   UT1-UTC Err  LOD Err     dX Err       dY Err
  119.         //                          "          "           s           s          "         "           "          "          s         s            "           "
  120.         //      (0h UTC)
  121.         //
  122.         COLUMNS_HEADER_PATTERN = Pattern.compile("^ *Date +MJD +x +y +UT1-UTC +LOD +((?:dPsi +dEps)|(?:dX +dY)) .*");

  124.         // The data lines in the EOP C04 yearly data files have the following fixed form:
  125.         // year month day MJD ...12 floating values fields in decimal format...
  126.         // 2000   1   1  51544   0.043242   0.377915   0.3554777   ...
  127.         // 2000   1   2  51545   0.043515   0.377753   0.3546065   ...
  128.         // 2000   1   3  51546   0.043623   0.377452   0.3538444   ...
  129.         // the corresponding fortran format is:
  130.         // 3(I4),I7,2(F11.6),2(F12.7),2(F12.6),2(F11.6),2(F12.7),2F12.6
  131.         DATA_LINE_PATTERN = Pattern.compile("^\\d+ +\\d+ +\\d+ +\\d+(?: +-?\\d+\\.\\d+){12}$");

  133.         YEAR_FIELD    = 0;
  134.         MONTH_FIELD   = 1;
  135.         DAY_FIELD     = 2;
  136.         MJD_FIELD     = 3;
  137.         POLE_X_FIELD  = 4;
  138.         POLE_Y_FIELD  = 5;
  139.         UT1_UTC_FIELD = 6;
  140.         LOD_FIELD     = 7;
  141.         NUT_0_FIELD   = 8;
  142.         NUT_1_FIELD   = 9;

  144.     }

  146.     /** Build a loader for IERS EOP C04 files.
  147.      * @param supportedNames regular expression for supported files names
  148.      * @param manager provides access to the EOP C04 files.
  149.      * @param utcSupplier UTC time scale.
  150.      */
  151.     EOPC04FilesLoader(final String supportedNames,
  152.                       final DataProvidersManager manager,
  153.                       final Supplier<TimeScale> utcSupplier) {
  154.         super(supportedNames, manager, utcSupplier);
  155.     }

  157.     /** {@inheritDoc} */
  158.     public void fillHistory(final IERSConventions.NutationCorrectionConverter converter,
  159.                             final SortedSet<EOPEntry> history) {
  160.         final ItrfVersionProvider itrfVersionProvider = new ITRFVersionLoader(
  161.                 ITRFVersionLoader.SUPPORTED_NAMES,
  162.                 getDataProvidersManager());
  163.         final Parser parser = new Parser(converter, itrfVersionProvider, getUtc());
  164.         final EopParserLoader loader = new EopParserLoader(parser);
  165.         this.feed(loader);
  166.         history.addAll(loader.getEop());
  167.     }

  169.     /** Internal class performing the parsing. */
  170.     static class Parser extends AbstractEopParser {

  172.         /**
  173.          * Simple constructor.
  174.          *
  175.          * @param converter           converter to use
  176.          * @param itrfVersionProvider to use for determining the ITRF version of the EOP.
  177.          * @param utc                 time scale for parsing dates.
  178.          */
  179.         Parser(final NutationCorrectionConverter converter,
  180.                final ItrfVersionProvider itrfVersionProvider,
  181.                final TimeScale utc) {
  182.             super(converter, itrfVersionProvider, utc);
  183.         }

  185.         /** {@inheritDoc} */
  186.         public Collection<EOPEntry> parse(final InputStream input, final String name)
  187.             throws IOException, OrekitException {

  189.             final List<EOPEntry> history = new ArrayList<>();
  190.             ITRFVersionLoader.ITRFVersionConfiguration configuration = null;

  192.             // set up a reader for line-oriented bulletin B files
  193.             final BufferedReader reader = new BufferedReader(new InputStreamReader(input, StandardCharsets.UTF_8));

  195.             // reset parse info to start new file (do not clear history!)
  196.             int lineNumber              = 0;
  197.             boolean inHeader            = true;
  198.             boolean isNonRotatingOrigin = false;

  200.             // read all file
  201.             for (String line = reader.readLine(); line != null; line = reader.readLine()) {
  202.                 ++lineNumber;
  203.                 boolean parsed = false;

  205.                 if (inHeader) {
  206.                     final Matcher matcher = COLUMNS_HEADER_PATTERN.matcher(line);
  207.                     if (matcher.matches()) {
  208.                         if (matcher.group(1).startsWith("dX")) {
  209.                             isNonRotatingOrigin = true;
  210.                         }
  211.                     }
  212.                 }

  214.                 if (DATA_LINE_PATTERN.matcher(line).matches()) {
  215.                     inHeader = false;
  216.                     // this is a data line, build an entry from the extracted fields
  217.                     final String[] fields = line.split(" +");
  218.                     final DateComponents dc = new DateComponents(Integer.parseInt(fields[YEAR_FIELD]),
  219.                                                                  Integer.parseInt(fields[MONTH_FIELD]),
  220.                                                                  Integer.parseInt(fields[DAY_FIELD]));
  221.                     final int    mjd   = Integer.parseInt(fields[MJD_FIELD]);
  222.                     if (dc.getMJD() != mjd) {
  223.                         throw new OrekitException(OrekitMessages.INCONSISTENT_DATES_IN_IERS_FILE,
  224.                                                   name, dc.getYear(), dc.getMonth(), dc.getDay(), mjd);
  225.                     }
  226.                     final AbsoluteDate date = new AbsoluteDate(dc, getUtc());

  228.                     // the first six fields are consistent with the expected format
  229.                     final double x     = Double.parseDouble(fields[POLE_X_FIELD]) * Constants.ARC_SECONDS_TO_RADIANS;
  230.                     final double y     = Double.parseDouble(fields[POLE_Y_FIELD]) * Constants.ARC_SECONDS_TO_RADIANS;
  231.                     final double dtu1  = Double.parseDouble(fields[UT1_UTC_FIELD]);
  232.                     final double lod   = Double.parseDouble(fields[LOD_FIELD]);
  233.                     final double[] equinox;
  234.                     final double[] nro;
  235.                     if (isNonRotatingOrigin) {
  236.                         nro = new double[] {
  237.                             Double.parseDouble(fields[NUT_0_FIELD]) * Constants.ARC_SECONDS_TO_RADIANS,
  238.                             Double.parseDouble(fields[NUT_1_FIELD]) * Constants.ARC_SECONDS_TO_RADIANS
  239.                         };
  240.                         equinox = getConverter().toEquinox(date, nro[0], nro[1]);
  241.                     } else {
  242.                         equinox = new double[] {
  243.                             Double.parseDouble(fields[NUT_0_FIELD]) * Constants.ARC_SECONDS_TO_RADIANS,
  244.                             Double.parseDouble(fields[NUT_1_FIELD]) * Constants.ARC_SECONDS_TO_RADIANS
  245.                         };
  246.                         nro = getConverter().toNonRotating(date, equinox[0], equinox[1]);
  247.                     }
  248.                     if (configuration == null || !configuration.isValid(mjd)) {
  249.                         // get a configuration for current name and date range
  250.                         configuration = getItrfVersionProvider().getConfiguration(name, mjd);
  251.                     }
  252.                     history.add(new EOPEntry(mjd, dtu1, lod, x, y, equinox[0], equinox[1], nro[0], nro[1],
  253.                                              configuration.getVersion(), date));
  254.                     parsed = true;

  256.                 }
  257.                 if (!(inHeader || parsed)) {
  258.                     throw new OrekitException(OrekitMessages.UNABLE_TO_PARSE_LINE_IN_FILE,
  259.                             lineNumber, name, line);
  260.                 }
  261.             }

  263.             // check if we have read something
  264.             if (inHeader) {
  265.                 throw new OrekitException(OrekitMessages.NOT_A_SUPPORTED_IERS_DATA_FILE, name);
  266.             }

  268.             return history;
  269.         }

  271.     }

  273. }
Created with JaCoCo 0.8.5.201910111838