ConvertC1ToP1.cpp

Go to the documentation of this file.

#pragma ident "$Id: ConvertC1ToP1.cpp 2552 2011-04-12 06:09:42Z yanweignss $"

//============================================================================
//
//  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., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
//  Wei Yan - Chinese Academy of Sciences . 2009, 2010, 2011
//
//============================================================================


#include "ConvertC1ToP1.hpp"

namespace gpstk
{
      // Index initially assigned to this class
   int ConvertC1ToP1::classIndex = 4800000;

      // Returns an index identifying this object.
   int ConvertC1ToP1::getIndex() const
   { return index; }


      // Returns a string identifying this object.
   std::string ConvertC1ToP1::getClassName() const
   { return "ConvertC1ToP1"; }

      /* Sets name of file containing DCBs data.
       * @param name      Name of the file containing DCB(P1-C1)
       */
   ConvertC1ToP1& ConvertC1ToP1::setDCBFile(const string& fileP1C1)
   {
      dcbP1C1.open(fileP1C1);
      
      return (*this);
   }

      /* Returns a satTypeValueMap object, adding the new data generated
       *  when calling this object.
       *
       * @param time      Epoch corresponding to the data.
       * @param gData     Data object holding the data.
       */
   satTypeValueMap& ConvertC1ToP1::Process( const DayTime& time,
                                                satTypeValueMap& gData )
      throw(ProcessingException)
   {
      
      try
      {
         static const double minLimit(15000000.0);
         static const double maxLimit(30000000.0);

         SatIDSet satRejectedSet;

         // Loop through all the satellites
         satTypeValueMap::iterator it;
         for (it = gData.begin(); it != gData.end(); ++it)
         {
            SatID sat = it->first;
            
            typeValueMap::iterator ittC1 = it->second.find(TypeID::C1);
            typeValueMap::iterator ittP1 = it->second.find(TypeID::P1);

            bool hasC1( ittC1 != it->second.end() );
            bool hasP1( ittP1 != it->second.end() );

            // filter out C1 and P1
            if(hasC1)
            {
               if( it->second[TypeID::C1]<minLimit || 
                   it->second[TypeID::C1]>maxLimit      ) { hasC1 = false;}
            }
            if(hasP1)
            {
               if( it->second[TypeID::P1]<minLimit || 
                   it->second[TypeID::P1]>maxLimit      ) { hasP1 = false;}
            }

               // If no desirable data, then schedule this satellite for removal
            if( !hasC1 && !hasP1)
            {
               satRejectedSet.insert(it->first);
            }

            if( hasC1 && !hasP1 )
            {
               double Bp1c1(0.0);      // in ns
               try
               {
                  Bp1c1 = dcbP1C1.getDCB(sat);
               }
               catch(...)
               {
                  Bp1c1 = 0.0;
               }

               it->second[TypeID::P1] = it->second[TypeID::C1] 
                                       +Bp1c1*(C_GPS_M * 1.0e-9);
            }

         }  // End of 'for (it = gData.begin(); it != gData.end(); ++it)'

         gData.removeSatID(satRejectedSet);

         return gData;
            
      }
      catch(Exception& u)
      {

         // Throw an exception if something unexpected happens
         ProcessingException e( getClassName() + ":"
            + StringUtils::asString( getIndex() ) + ":"
            + u.what() );

         GPSTK_THROW(e);

      }

   }  // End of method 'ConvertC1ToP1::Process()'


}  // End of namespace gpstk

---