SinexBase.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 "SinexBase.hpp"
#include "YDSTime.hpp"
#include "FormattedDouble.hpp"
 
using namespace gnsstk::StringUtils;
using namespace std;
 
namespace gnsstk
{
namespace Sinex
{
 
   bool
   isValidObsCode(char c, bool toss)
   {
      if (OBS_CODES.find(c) == string::npos)
      {
         if (toss)
         {
            string  errMsg("Invalid Observation Code: ");
            errMsg += c;
            Exception  err(errMsg);
            GNSSTK_THROW(err);
         }
         return false;
      }
      return true;
   }
 
 
   bool
   isValidConstraintCode(char c, bool toss)
   {
      if (CONSTRAINT_CODES.find(c) == string::npos)
      {
         if (toss)
         {
            string  errMsg("Invalid Constraint Code: ");
            errMsg += c;
            Exception  err(errMsg);
            GNSSTK_THROW(err);
         }
         return false;
      }
      return true;
   }
 
 
   bool
   isValidSolutionType(char c, bool toss)
   {
      if (SOLUTION_TYPES.find(c) == string::npos)
      {
         if (toss)
         {
            string  errMsg("Invalid Solution Type: ");
            errMsg += c;
            Exception  err(errMsg);
            GNSSTK_THROW(err);
         }
         return false;
      }
      return true;
   }
 
 
   bool
   isValidLineStructure(const std::string& line,
                                   size_t minLen,
                                   size_t maxLen,
                                   int divs[],
                                   bool toss)
   {
      size_t  sz = line.size();
 
      if (minLen < MIN_LINE_LEN) minLen = MIN_LINE_LEN;
      if (sz < minLen)
      {
         if (toss)
         {
            ostringstream  ss;
            ss << "Missing data; inadequate line length ("
               << sz << " < " << minLen << ")";
            Exception  err(ss.str() );
            GNSSTK_THROW(err);
         }
         return false;
      }
 
      if (maxLen > MAX_LINE_LEN) maxLen = MAX_LINE_LEN;
      if (sz > maxLen)
      {
         if (toss)
         {
            ostringstream  ss;
            ss << "Excessive line length ("
               << sz << " > " << maxLen << ")";
            Exception  err(ss.str() );
            GNSSTK_THROW(err);
         }
         return false;
      }
      if (divs != NULL)
      {
         int  pos;
         for (size_t i = 0; (pos = divs[i]) >= 0; ++i)
         {
            if ( (pos >= sz) || (line[pos] != FIELD_DIV) )
            {
               if (toss)
               {
                  ostringstream  ss;
                  ss << "Field divider '" << FIELD_DIV
                     << "' expected in column " << pos;
                  Exception  err(ss.str() );
                  GNSSTK_THROW(err);
               }
               return false;
            }
         }
      }
      return true;
   }
 
 
   std::string
   formatStr(const std::string& value, size_t width, bool allowBlank)
   {
      if ( (value.size() == 0) && !allowBlank)
      {
         return string(width, '-');
      }
      else
      {
         ostringstream  ss;
         ss << left << setw(width) << value.substr(0, width);
         return ss.str();
      }
   }
 
 
   std::string
   formatInt(long value, size_t width)
   {
      ostringstream  ss;
      ss << right << setw(width) << value;
      string  s = ss.str();
      if (s.size() > width)
      {
         ostringstream  msg;
         msg << "Cannot represent value " << value;
         msg << " within width " << width;
         Exception  err(msg.str() );
         GNSSTK_THROW(err);
      }
      return s;
   }
 
 
   std::string
   formatUint(unsigned long value, size_t width)
   {
      ostringstream  ss;
      ss << right << setw(width) << value;
      string  s = ss.str();
      if (s.size() > width)
      {
         ostringstream  msg;
         msg << "Cannot represent value " << value;
         msg << " within width " << width;
         Exception  err(msg.str() );
         GNSSTK_THROW(err);
      }
      return s;
   }
 
 
   std::string
   formatFor(double value, size_t width, size_t expLen)
   {
      FormattedDouble fd(value, StringUtils::FFLead::Zero, width-6, expLen,
                         width, 'E');
      return static_cast<std::string>(fd);
   }
 
   std::string
   formatFixed(double value, size_t width, size_t precision)
   {
      ostringstream  ss;
      ss << fixed << right << setw(width) << setprecision(precision) << value;
      string  s = ss.str();
 
      if (s.size() > width)
      {
         if ( (s.size() == width+1) && (value < 1.0) && (value > -1.0) )
         {
            if (value < 0.0)
            {
               s.erase(1, 1);  
            }
            else
            {
               s.erase(0, 1);  
            }
         }
         else
         {
            ostringstream  msg;
            msg << "Cannot represent value " << value;
            msg << " within width " << width;
            msg << " with precision " << precision;
            Exception  err(msg.str() );
            GNSSTK_THROW(err);
         }
      }
      return s;
   }
 
 
   Time::operator gnsstk::CommonTime() const
   {
      long  y = (year > 50) ? (year + 1900) : (year + 2000);
      gnsstk::YDSTime  ydsTime(y, doy, (double)sod);
      return ydsTime.convertToCommonTime();
   }
 
 
   Time::operator std::string() const
   {
      std::ostringstream  ss;
      ss << setfill('0');
      ss << setw(2) << (uint16_t)year << ":";
      ss << setw(3) << doy << ":";
      ss << setw(5) << sod;
      return ss.str();
   }
 
 
   void
   Time::operator=(const gnsstk::CommonTime& other)
   {
      gnsstk::YDSTime  ydsTime(other);
      if (ydsTime.year > 2050)
      {
         Exception  err("SINEX does not support years past 2050.");
         GNSSTK_THROW(err);
      }
      year = ydsTime.year % 100;  // Get last 2 digits of time
      doy  = ydsTime.doy;
      sod  = ydsTime.sod;
   }
 
 
   void
   Time::operator=(const std::string& other)
   {
      if ( (other.size() < 12) || (other[2] != ':') || (other[6] != ':') )
      {
         Exception  err("Invalid time syntax: " + other);
         GNSSTK_THROW(err);
      }
      year = asInt(other.substr(0, 2) );
      doy  = asInt(other.substr(3, 3) );
      sod  = asInt(other.substr(7, 5) );
   }
 
}  // namespace Sinex
 
}  // namespace gnsstk