DCBDataReader.cpp

Go to the documentation of this file.

//==============================================================================
//
//  This file is part of GNSSTk, the ARL:UT GNSS Toolkit.
//
//  The GNSSTk is free software; you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published
//  by the Free Software Foundation; either version 3.0 of the License, or
//  any later version.
//
//  The GNSSTk is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU Lesser General Public License for more details.
//
//  You should have received a copy of the GNU Lesser General Public
//  License along with GNSSTk; if not, write to the Free Software Foundation,
//  Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin.
//  Copyright 2004-2022, The Board of Regents of The University of Texas System
//
//==============================================================================
 
//==============================================================================
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin, under contract to an agency or agencies
//  within the U.S. Department of Defense. The U.S. Government retains all
//  rights to use, duplicate, distribute, disclose, or release this software.
//
//  Pursuant to DoD Directive 523024
//
//  DISTRIBUTION STATEMENT A: This software has been approved for public
//                            release, distribution is unlimited.
//
//==============================================================================
 
#include "DCBDataReader.hpp"
 
using namespace std;
 
namespace gnsstk
{
 
      // Method to store load ocean tide harmonics data in this class'
      // data map
   void DCBDataReader::loadData()
   {
 
      try
      {
         allDCB.satDCB.clear();
         allDCB.gpsDCB.clear();
         allDCB.glonassDCB.clear();
 
            // a buffer
         string line;
 
            // read first line
         formattedGetLine(line, true);
 
            // Let's skip 6 lines
         for(int i=0; i<6; i++) formattedGetLine(line, true);
 
 
            // Now, let's read data
         while(1)
         {
            formattedGetLine(line, true);
 
            if(line.length() < 46) continue;
 
            string sysFlag = line.substr(0,1);
 
            int satPRN = StringUtils::asInt(line.substr(1,2));
 
            string station = StringUtils::strip(line.substr(6,4));
 
            const double dcbVal = StringUtils::asDouble(line.substr(26,9));
            //const double dcbRms = StringUtils::asDouble(line.substr(38,9));
 
            if(station.length() < 4)       // this is satellite DCB data
            {
 
               SatID sat;
               if(sysFlag == "G")
               {
                  sat = SatID(satPRN,SatelliteSystem::GPS);
               }
               else if(sysFlag == "R")
               {
                  sat = SatID(satPRN,SatelliteSystem::Glonass);
               }
               else
               {
                  // Unexpected and we do nothing here
 
 
               }
 
               allDCB.satDCB[sat] = dcbVal;
 
            }
            else                           // this is receiver DCB data
            {
               if(sysFlag == "G")
               {
                  allDCB.gpsDCB[station] = dcbVal;
               }
               else if(sysFlag == "R")
               {
                  allDCB.glonassDCB[station] = dcbVal;
               }
               else
               {
                  // Unexpected and we do nothing here
 
               }
            }
 
         }  // End of 'while(1)'
 
      }  // End of try block
      catch (EndOfFile& e)
      {
 
            // We should close this data stream before returning
         (*this).close();
 
         return;
      }
      catch (...)
      {
 
         // We should close this data stream before returning
         (*this).close();
 
         return;
 
      }
 
 
   }  // End of 'DCBDataReader::loadData()'
 
 
 
      // Method to open AND load DCB data file.
   void DCBDataReader::open(const char* fn)
   {
 
      // We need to be sure current data stream is closed
      (*this).close();
 
      // Open data stream
      FFTextStream::open(fn, std::ios::in);
      loadData();
 
      return;
 
   }  // End of method 'DCBDataReader::open()'
 
 
 
      // Method to open AND load DCB data file. It doesn't
      // clear data previously loaded.
   void DCBDataReader::open(const std::string& fn)
   {
 
      // We need to be sure current data stream is closed
      (*this).close();
 
      // Open data stream
      FFTextStream::open(fn.c_str(), std::ios::in);
      loadData();
 
      return;
   }  // End of method 'DCBDataReader::open()'
 
      // return P1-P2 or P1-C1 depend what you have loaded
   double DCBDataReader::getDCB(const SatID& sat)
   {
      return allDCB.satDCB[sat];
   }
 
      // Get DCB data of a satellite
      // return P1-P2 or P1-C1 depend what you have loaded
   double DCBDataReader::getDCB(const int& prn,
      const SatelliteSystem& system)
   {
      SatID sat(prn,system);
      return allDCB.satDCB[sat];
   }
 
      // Get DCB data of aReceiver
      // it return P1-P2
   double DCBDataReader::getDCB(const std::string& station,
                                const SatelliteSystem& system)
   {
 
      if(system == SatelliteSystem::GPS)
      {
         return allDCB.gpsDCB[station];
      }
      else if(system == SatelliteSystem::Glonass)
      {
         return allDCB.glonassDCB[station];
      }
      else
      {
            // Unexpected and return 0
         return 0.0;
      }
 
   }  // End of 'double DCBDataReader::getDCB(const string& station...'
 
 
 
}  // End of namespace gnsstk