RinexMetFilterOperators.hpp

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.
//
//==============================================================================
 
#ifndef GNSSTK_RINEXMETFILTEROPERATORS_HPP
#define GNSSTK_RINEXMETFILTEROPERATORS_HPP
 
#include <set>
 
#include "CivilTime.hpp"
#include "FileFilter.hpp"
#include "RinexMetData.hpp"
#include "RinexMetHeader.hpp"
 
namespace gnsstk
{
 
 
   typedef std::unary_function<RinexMetHeader, bool> RinexMetDataUnaryOperator;
   typedef std::binary_function<RinexMetData, RinexMetData, bool> RinexMetDataBinaryOperator;
 
   struct RinexMetDataOperatorLessThanFull : public RinexMetDataBinaryOperator
   {
   public:
 
      RinexMetDataOperatorLessThanFull
      (const std::set<RinexMetHeader::RinexMetType>& rmhset)
            : obsSet(rmhset)
      {}
 
      bool operator()(const RinexMetData& l, const RinexMetData& r) const
      {
            // Compare the times, offsets, then only those elements
            // common to both.  This ignores the flags set to 0.
 
         if (l.time < r.time)
            return true;
         else if (l.time != r.time)
            return false;
 
            // Then check that each observation has the same data
            // for each item in the set of common observations.
 
         RinexMetData::RinexMetMap::const_iterator
            lItr, rItr;
         std::set<RinexMetHeader::RinexMetType>::const_iterator
            obsItr = obsSet.begin();
 
         while (obsItr != obsSet.end())
         {
            rItr = r.data.find(*obsItr);
            if (rItr == r.data.end())
               return false;
 
            lItr = l.data.find(*obsItr);
            if (lItr == l.data.end())
               return false;
 
            if ((*lItr).second < (*rItr).second)
               return true;
            if ((*lItr).second > (*rItr).second)
               return false;
 
            obsItr++;
         }
 
            // the data is either == or > at this point
         return false;
      }
 
   private:
      std::set<RinexMetHeader::RinexMetType> obsSet;
   };
 
   struct RinexMetDataOperatorLessThanSimple : public RinexMetDataBinaryOperator
   {
   public:
 
      bool operator()(const RinexMetData& l, const RinexMetData& r) const
      {
         if (l.time < r.time)
            return true;
         return false;
      }
   };
 
   struct RinexMetDataOperatorEqualsSimple : public RinexMetDataBinaryOperator
   {
   public:
 
      bool operator()(const RinexMetData& l, const RinexMetData& r) const
      {
         if (l.time == r.time)
            return true;
         return false;
      }
   };
 
   struct RinexMetHeaderTouchHeaderMerge : public RinexMetDataUnaryOperator
   {
   public:
 
      RinexMetHeaderTouchHeaderMerge()
            : firstHeader(true)
      {}
 
      bool operator()(const RinexMetHeader& l)
      {
         if (firstHeader)
         {
            theHeader = l;
            firstHeader = false;
         }
         else
         {
            std::set<RinexMetHeader::RinexMetType> thisMetSet,
               tempMetSet;
            std::set<std::string> commentSet;
            obsSet.clear();
 
               // insert the comments to the set and let the set take care of uniqueness
            copy(theHeader.commentList.begin(),
                 theHeader.commentList.end(),
                 inserter(commentSet, commentSet.begin()));
            copy(l.commentList.begin(),
                 l.commentList.end(),
                 inserter(commentSet, commentSet.begin()));
               // then copy the comments back into theHeader
            theHeader.commentList.clear();
            copy(commentSet.begin(), commentSet.end(),
                 inserter(theHeader.commentList,
                          theHeader.commentList.begin()));
 
               // find the set intersection of the obs types
            copy(theHeader.obsTypeList.begin(),
                 theHeader.obsTypeList.end(),
                 inserter(thisMetSet, thisMetSet.begin()));
            copy(l.obsTypeList.begin(),
                 l.obsTypeList.end(),
                 inserter(tempMetSet, tempMetSet.begin()));
            set_intersection(thisMetSet.begin(), thisMetSet.end(),
                             tempMetSet.begin(), tempMetSet.end(),
                             inserter(obsSet, obsSet.begin()));
               // then copy the obsTypes back into theHeader
            theHeader.obsTypeList.clear();
            copy(obsSet.begin(), obsSet.end(),
                 inserter(theHeader.obsTypeList,
                          theHeader.obsTypeList.begin()));
         }
         return true;
      }
 
      bool firstHeader;
      RinexMetHeader theHeader;
      std::set<RinexMetHeader::RinexMetType> obsSet;
   };
 
   struct RinexMetDataFilterTime : public RinexMetDataUnaryOperator
   {
 
   public:
 
      RinexMetDataFilterTime(const CommonTime& startTime,
                             const CommonTime& endTime   )
            : start(startTime), end(endTime)
      {}
 
      bool operator() (const RinexMetData& l) const
      {
         if ( l.time < start || l.time >= end )
            return true;
         return false;
      }
 
   private:
 
      CommonTime start, end;
 
   };
 
 
} // namespace gnsstk
 
#endif