ComputeMelbourneWubbena.cpp

Go to the documentation of this file.

#pragma ident "$Id: ComputeMelbourneWubbena.cpp 1325 2008-07-29 14:33:43Z architest $"

//============================================================================
//
//  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
//
//  Dagoberto Salazar - gAGE ( http://www.gage.es ). 2007, 2008
//
//============================================================================


#include "ComputeMelbourneWubbena.hpp"


namespace gpstk
{

      // Index initially assigned to this class
   int ComputeMelbourneWubbena::classIndex = 2200000;


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


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


      // Default constructor
   ComputeMelbourneWubbena::ComputeMelbourneWubbena()
      : type3(TypeID::L1), type4(TypeID::L2), DEN1(L1_FREQ + L2_FREQ),
        DEN2(L1_FREQ - L2_FREQ)
   {
      type1 = TypeID::P1;
      type2 = TypeID::P2;
      resultType = TypeID::MWubbena;
      setIndex();
   }


      /* Returns a satTypeValueMap object, adding the new data generated
       *  when calling this object.
       *
       * @param gData     Data object holding the data.
       */
   satTypeValueMap& ComputeMelbourneWubbena::Process(satTypeValueMap& gData)
      throw(ProcessingException)
   {

      try
      {

         double value1(0.0);
         double value2(0.0);
         double value3(0.0);
         double value4(0.0);

         SatIDSet satRejectedSet;

            // Loop through all the satellites
         satTypeValueMap::iterator it;
         for (it = gData.begin(); it != gData.end(); ++it) 
         {
            try
            {
                  // Try to extract the values
               value1 = (*it).second(type1);
               value2 = (*it).second(type2);
               value3 = (*it).second(type3);
               value4 = (*it).second(type4);
            }
            catch(...)
            {
                  // If some value is missing, then schedule this satellite
                  // for removal
               satRejectedSet.insert( (*it).first );
               continue;
            }

               // If everything is OK, then get the new value inside
               // the structure
            (*it).second[resultType] = getCombination( value1,
                                                       value2,
                                                       value3,
                                                       value4 );
         }

            // Remove satellites with missing data
         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 'ComputeMelbourneWubbena::Process()'



      // Compute the combination of observables.
   double ComputeMelbourneWubbena::getCombination( const double& p1,
                                                   const double& p2,
                                                   const double& l1,
                                                   const double& l2 )
   {

      return (   ( L1_FREQ*l1 - L2_FREQ*l2 ) / ( DEN2 )
               - ( L1_FREQ*p1 + L2_FREQ*p2 ) / ( DEN1 ) );

   }  // End of method 'ComputeMelbourneWubbena::getCombination()'



} // End of namespace gpstk

---