Rinex3ObsFilterOperators.hpp

Go to the documentation of this file.

#pragma ident "$Id: Rinex3ObsFilterOperators.hpp 3319 2012-09-19 16:58:10Z prestonherrmann $"

#ifndef GPSTK_RINEX3OBSFILTEROPERATORS_HPP
#define GPSTK_RINEX3OBSFILTEROPERATORS_HPP

//============================================================================
//
//  This file is part of GPSTk, the GPS Toolkit.
//
//  The GPSTk 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 2.1 of the License, or
//  any later version.
//
//  The GPSTk 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 GPSTk; if not, write to the Free Software Foundation,
//  Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//  
//  Copyright 2004, The University of Texas at Austin
//
//============================================================================

//============================================================================
//
//This software 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 <set>
#include <algorithm>

#include "FileFilter.hpp"
#include "Rinex3ObsData.hpp"
#include "Rinex3ObsHeader.hpp"
#include "ObsID.hpp"

namespace gpstk
{

   struct Rinex3ObsDataOperatorLessThanFull : 
      public std::binary_function<Rinex3ObsData, Rinex3ObsData, bool>
   {
   public:
      Rinex3ObsDataOperatorLessThanFull
      (const std::vector<ObsID>& rohset)
         : obsSet(rohset)
      {}

      bool operator()(const Rinex3ObsData& l, const Rinex3ObsData& r) const
         {
               // compare the times, offsets, then only those elements
               // that are common to both.  this ignores the flags
               // that are set to 0
            if (l.time < r.time)
               return true;
            else if (l.time == r.time)
            {
               if (l.epochFlag < r.epochFlag)
                  return true;
               else if (l.epochFlag == r.epochFlag)
               {
                  if (l.clockOffset < r.clockOffset)
                     return true;
                  else if (l.clockOffset > r.clockOffset)
                     return false;
               }
               else
                  return false;
            }
            else
               return false;
         
               // for the obs, first check that they're the same size
               // i.e. - contain the same number of PRNs
            if (l.obs.size() < r.obs.size())
               return true;

            if (l.obs.size() > r.obs.size())
               return false;

               // then check that each PRN has the same data for each of the 
               // shared fields
            Rinex3ObsData::DataMap::const_iterator lItr = l.obs.begin(), rItr;
         
            SatID sat;

            while (lItr != l.obs.end())
            {
               sat = (*lItr).first;
               rItr = r.obs.find(sat);
               if (rItr == r.obs.end())
                  return false;
         
               std::vector<Rinex3ObsData::RinexDatum> lObs = lItr->second,
                                                      rObs = rItr->second;

         
               for(int i = 0; i < obsSet.size(); ++i)
               {
                  Rinex3ObsData::RinexDatum lData, rData;
                  lData = lObs[i];
                  rData = rObs[i];

                  if (lData.data < rData.data)
                     return true;
               
                  if ( lData.lli != 0 && rData.lli != 0 )
                     if (lData.lli < rData.lli)
                        return true;
               
                  if ( lData.ssi != 0 && rData.ssi != 0 )
                     if (lData.ssi < rData.ssi)
                        return true;
               
               }

               lItr++;
            }

               // the data is either == or > at this point
            return false;
         }

   private:
      std::vector<ObsID> obsSet;
   };

   struct Rinex3ObsDataOperatorLessThanSimple : 
      public std::binary_function<Rinex3ObsData, Rinex3ObsData, bool>
   {
   public:
      bool operator()(const Rinex3ObsData& l, const Rinex3ObsData& r) const
         {
            if (l.time < r.time)
               return true;
            return false;
         }
   };

   struct Rinex3ObsDataOperatorEqualsSimple : 
      public std::binary_function<Rinex3ObsData, Rinex3ObsData, bool>
   {
   public:
      bool operator()(const Rinex3ObsData& l, const Rinex3ObsData& r) const
         {
            if (l.time == r.time)
               return true;
            return false;
         }
   };

   struct Rinex3ObsHeaderTouchHeaderMerge :
      public std::unary_function<Rinex3ObsHeader, bool>
   {
   public:
      Rinex3ObsHeaderTouchHeaderMerge()
         : firstHeader(true)
      {}

      bool operator()(const Rinex3ObsHeader& l)
         {
            if (firstHeader)
            {
               theHeader = l;
               firstHeader = false;
            }
            else
            {
               std::vector<ObsID> thisObsSet, tempObsSet;
               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(thisObsSet, thisObsSet.begin()));
               copy(l.obsTypeList.begin(),
                    l.obsTypeList.end(),
                    inserter(tempObsSet, tempObsSet.begin()));
               set_intersection(thisObsSet.begin(), thisObsSet.end(),
                                tempObsSet.begin(), tempObsSet.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;
      Rinex3ObsHeader theHeader;
      std::vector<ObsID> obsSet;
   };


} // namespace gpstk

#endif // GPSTK_RINEX3OBSFILTEROPERATORS_HPP

---