Logo Search packages:      
Sourcecode: magics++ version File versions  Download package

ObsDecoder.cc

Go to the documentation of this file.
/******************************** LICENSE ********************************

 Copyright 2007 European Centre for Medium-Range Weather Forecasts (ECMWF)

 Licensed under the Apache License, Version 2.0 (the "License");
 you may not use this file except in compliance with the License.
 You may obtain a copy of the License at 

    http://www.apache.org/licenses/LICENSE-2.0

 Unless required by applicable law or agreed to in writing, software
 distributed under the License is distributed on an "AS IS" BASIS,
 WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 See the License for the specific language governing permissions and
 limitations under the License.

 ******************************** LICENSE ********************************/

/*! \file ObsDecoder.cc
    \brief Implementation of the Template class ObsDecoder.
    
    Magics Team - ECMWF 2005
    
    Started: Wed 23-Mar-2005
    
    Changes:
    
*/

#include "ObsDecoder.h"
#include "MvObs.h"
#include "TextVisitor.h"
#include "CustomisedPoint.h"
#include "expat.h"
#include "Factory.h"

using namespace magics;

map<int, bool> multilevels_;

class  BufrSubType : public map<string, string> 
{
public:
      BufrSubType() {}
      ~BufrSubType() {}
protected :
      void print(ostream& out) const {
            for (const_iterator subtype = begin(); subtype != end(); ++subtype) 
                  out << "subtype[" << subtype->first << "=" << subtype->second << "]";
      } 
      //! Overloaded << operator to call print().
      friend ostream& operator<<(ostream& s,const BufrSubType& p)
            { p.print(s); return s; }
};


class BufrFamily : public map<string, BufrSubType> 
{
public:
      BufrFamily(const string&);
      ~BufrFamily() {};
      
      void currentType(const map<string, string>& def) { type_ = def.find("value")->second; }

      void currentSubtype(const map<string, string>& def)
      {
            map<string, string>::const_iterator value = def.find("value");
            if ( value == def.end() )
            {
                  Log::warning() << " No code defined for subtype--> Contact graphics team.\n";
                  subtype_ = "???";
            }
            else 
              subtype_ = def.find("value")->second;

            map<string, string>::const_iterator tmpl = def.find("template");
            template_ = ( tmpl != def.end() ) ? tmpl->second : "";
            subtypeToBeSet_ = true; 
      }

      void addSubtype(const string& value)
      {
            if (template_ == "") template_ = value;
            (*this)[type_].insert(make_pair(subtype_, template_)); 
            subtypeToBeSet_ = false;
      }

      bool subtype() { return subtypeToBeSet_; }

      void subtypeToReset() { subtypeToBeSet_ = false; }
      
      string getType(const string& type, const string& subtype) const
      {
            const_iterator t = find(type);
            if ( t != end() )
            {
                  BufrSubType::const_iterator s = t->second.find(subtype);
                  if (s != t->second.end() ) return s->second;
            }
            return "";
      }

protected :
      void print(ostream& out) const
      {
            for (const_iterator type = begin(); type != end(); ++type) 
                  out << type->first << "=[" << type->second << "]";
      } 

      string type_;
      string subtype_;
      string template_;
      bool   subtypeToBeSet_;
      string centre_;

      //! Overloaded << operator to call print().
      friend ostream& operator<<(ostream& s,const BufrFamily& p)
            { p.print(s); return s; }
};


class BufrFamilyTable : public map<string, BufrFamily*>
{
public :
      BufrFamilyTable() {}
      ~BufrFamilyTable() {}
      static const BufrFamily& get(const string& centre)
      {
            const_iterator table = table_.find(centre);
            if ( table != table_.end() ) return *(table->second);
            BufrFamily*  idents = new BufrFamily(centre);
            table_[centre] = idents;
            return *idents; 
      }

      static BufrFamilyTable table_;
};


static void XMLCALL
startElement(void *userData, const char *name, const char **atts)
{
      BufrIdentifiers* ident  = (BufrIdentifiers*) userData; 
      
        while (*atts) {
            if (string(*atts) == "descriptor") {
                  ident->insert(make_pair(name,*(atts+1)));
            }
            atts+=2;
        }
}

static void XMLCALL
endElement(void *, const char *)
{
}


static void XMLCALL
startFamilyElement(void *userData, const char *name, const char **atts)
{
      BufrFamily* ident  = (BufrFamily*) userData; 
      string token(name);
      ident->subtypeToReset();
      map<string, string> def;
      if (token != "type" && token != "subtype") return;
      while (*atts) {         
          def.insert(make_pair(*atts,*(atts+1)));
          atts+=2;
      }
      if ( token == "type") { ident->currentType(def); }
      if ( token == "subtype") { ident->currentSubtype(def); }
}

static void XMLCALL
endFamilyElement(void *, const char *)
{
}

static void XMLCALL character (void *userData,
                            const char *s,
                            int len)
{
      //int *depthPtr = (int*)userData;
      BufrFamily* object  = (BufrFamily*) userData;
      if (  std::string(s, len) == "\n" ) return;
      if (object->subtype()) {
            object->addSubtype(string(s, len));
      }
}


BufrIdentifiers::BufrIdentifiers(int centre) : centre_(centre) 
{
      ostringstream file;
      file << getEnvVariable("MAGPLUS_HOME") <<  MAGPLUS_PATH_TO_SHARE_ << "bufr_" << centre << ".xml";
      char buf[BUFSIZ];
      XML_Parser parser = XML_ParserCreate(NULL);
      int done;
      XML_SetUserData(parser, this);
      XML_SetElementHandler(parser, startElement, endElement);

      FILE* in  = fopen(file.str().c_str(), "r");
      if (!in) return;

      do
      {
            size_t len = fread(buf, 1, sizeof(buf), in);
            done = len < sizeof(buf);
            if (XML_Parse(parser, buf, len, done) == XML_STATUS_ERROR)
            {
                  ostringstream s;
                  s << "XmlException : " << XML_ErrorString(XML_GetErrorCode(parser))  << " at line  " <<  XML_GetCurrentLineNumber(parser)  << ends;
                  cerr <<  s.str() << "\n";
            }
      } while (!done);
      XML_ParserFree(parser);
      fclose(in);
}




int BufrIdentifiers::ident(const string& token) const
{
      const_iterator ident = find(token);
      if ( ident != end() ) return atoi(ident->second.c_str());
      if (centre_ == 98) return -1;
      else return BufrIdentTable::get(98).ident(token);
}


BufrIdentTable BufrIdentTable::table_;
BufrFamilyTable BufrFamilyTable::table_;

BufrFamily::BufrFamily(const string& centre) : centre_(centre)
{
      ostringstream file;
      file << getEnvVariable("MAGPLUS_HOME") <<  MAGPLUS_PATH_TO_SHARE_ << "bufr_" << centre << ".xml";
      char buf[BUFSIZ];
      XML_Parser parser = XML_ParserCreate(NULL);
      int done;
      XML_SetUserData(parser, this);
      XML_SetElementHandler(parser, startFamilyElement, endFamilyElement);
      XML_SetCharacterDataHandler(parser, character);

      FILE* in  = fopen(file.str().c_str(), "r");
      if (!in) return;

      do
      {
            size_t len = fread(buf, 1, sizeof(buf), in);
            done = len < sizeof(buf);
            if (XML_Parse(parser, buf, len, done) == XML_STATUS_ERROR)
            {
                  ostringstream s;
                  s << "XmlException : " << XML_ErrorString(XML_GetErrorCode(parser))  << " at line  " <<  XML_GetCurrentLineNumber(parser)  << ends;
                  cerr <<  s.str() << "\n";
            }
      } while (!done);
      XML_ParserFree(parser);
      fclose(in);
}

const BufrIdentifiers& BufrIdentTable::get(int centre)
{
      const_iterator table = table_.find(centre);
      if ( table != table_.end() ) return *(table->second);
      BufrIdentifiers*  idents = new BufrIdentifiers(centre);
      table_[centre] = idents;
      return *idents; 
}


class BufrAccessor 
{
public:
      BufrAccessor() {};
      BufrAccessor(const string& descriptor) : descriptor_(descriptor) {}
      virtual ~BufrAccessor() {};
      virtual void operator()(const ObsDecoder&, MvObs&, string&) const {}
      virtual void operator()(const ObsDecoder&, MvObs& obs, double& val) const
      {
            const BufrIdentifiers& table =  BufrIdentTable::get(obs.originatingCentre());
            val = obs.value(table.ident(descriptor_));
            Log::dev()<< "BufrAccessor-Descriptor--->" << descriptor_ << " Value--->" << val << endl;
      }
      virtual void print() {}
      const string& keyword() { return descriptor_; }
protected:
      string descriptor_;
};

class BufrMultiValueAccessor : public BufrAccessor 
{
public:
      BufrMultiValueAccessor(const string& descriptor, int index) : BufrAccessor(descriptor), index_(index) {}
      virtual ~BufrMultiValueAccessor() {};
      virtual void operator()(const ObsDecoder&, MvObs& obs, double& val) const {
            const BufrIdentifiers& table =  BufrIdentTable::get(obs.originatingCentre());
            val = obs.valueByOccurrence(index_, table.ident(descriptor_));
            Log::dev()<< "BufrMultiValueAccessor-Descriptor--->" << descriptor_ << " INDEX--->" << index_ << " Value--->" << val << endl;
      }
      virtual void print() {}
protected:
      int index_;
};

class BufrLowCloudAccessor : public BufrMultiValueAccessor 
{
public:
      BufrLowCloudAccessor() : BufrMultiValueAccessor("low_cloud", 1) {}
      virtual ~BufrLowCloudAccessor() {}
};

class BufrMediumCloudAccessor : public BufrMultiValueAccessor 
{
public:
      BufrMediumCloudAccessor() : BufrMultiValueAccessor("medium_cloud", 2) {}
      virtual ~BufrMediumCloudAccessor() {}
};

class BufrHighCloudAccessor : public BufrMultiValueAccessor 
{
public:
      BufrHighCloudAccessor() : BufrMultiValueAccessor("high_cloud", 3) {}
      virtual ~BufrHighCloudAccessor() {}
};

class BufrLowCloudNebulosityAccessor : public BufrMultiValueAccessor 
{
public:
      BufrLowCloudNebulosityAccessor() : BufrMultiValueAccessor("low_cloud_nebulosity", 1) {}
      virtual ~BufrLowCloudNebulosityAccessor() {}
};

class BufrMediumCloudNebulosityAccessor : public BufrMultiValueAccessor 
{
public:
      BufrMediumCloudNebulosityAccessor() : BufrMultiValueAccessor("medium_cloud_nebulosity", 2) {}
      virtual ~BufrMediumCloudNebulosityAccessor() {}
};

class BufrHighCloudNebulosityAccessor : public BufrMultiValueAccessor 
{
public:
      BufrHighCloudNebulosityAccessor() : BufrMultiValueAccessor("high_cloud_nebulosity", 3) {}
      virtual ~BufrHighCloudNebulosityAccessor() {}
};
class BufrLowCloudHeightAccessor : public BufrMultiValueAccessor 
{
public:
      BufrLowCloudHeightAccessor() : BufrMultiValueAccessor("low_cloud_height", 1) {}
      virtual ~BufrLowCloudHeightAccessor() {}
};

class BufrMediumCloudHeightAccessor : public BufrMultiValueAccessor 
{
public:
      BufrMediumCloudHeightAccessor() : BufrMultiValueAccessor("medium_cloud_height", 2) {}
      virtual ~BufrMediumCloudHeightAccessor() {}
};

class BufrHighCloudHeightAccessor : public BufrMultiValueAccessor 
{
public:
      BufrHighCloudHeightAccessor() : BufrMultiValueAccessor("high_cloud_height", 3) {}
      virtual ~BufrHighCloudHeightAccessor() {}
};


class BufrTypeAccessor : public BufrAccessor 
{
public:
      BufrTypeAccessor() { descriptor_ = "type"; };
      virtual ~BufrTypeAccessor() {};
      void operator()(const ObsDecoder&, MvObs& obs, string& val) const {
            val = obs.messageType();
          string type = tostring(obs.messageType());
          string subtype = tostring(obs.messageSubtype());
          string centre = tostring(obs.originatingCentre());
            const BufrFamily& family =  BufrFamilyTable::get(centre);               
            val = family.getType(type, subtype);            
            if (val == "") {
                  Log::warning() << "Cannot find template type for [" << centre << ", " << type << ", " << subtype << "][centre, type, subtype]" << endl; 
                  Log::warning() << "Contact graphics team.\n";
            }
      }
      
};

class BufrThicknessAccessor : public BufrAccessor
{
public:
      BufrThicknessAccessor() { descriptor_ = "thickness"; keyword_ = "geopotential"; }

      virtual ~BufrThicknessAccessor() {}
      void operator()(const ObsDecoder& decoder, MvObs& obs, double& val) const {
            const BufrIdentifiers& table =  BufrIdentTable::get(obs.originatingCentre()); 
            int type  = obs.messageType();
            
            map<int, bool>::const_iterator multilevel = multilevels_.find(type);
            if ( multilevel == multilevels_.end() ) {
                  Log::warning() << "Unknown observation type [" << val << "]: Contact graphics group\n";
                  val =  kBufrMissingValue;
            }
            if ( type == 0 || type == 1) {
                  // surface data 
                  val = 0;                
            }
            else {
                  //Log::dev()<< " look for --->" << table.ident(altitude_) << " at " << decoder.getLevel();
                  // upper-air data 
                  if  ( type == 5 || type == 4 ) 
                        val = 0;                
                  else // Multi-level data {
                        val = abs(obs.valueByPressureLevel(decoder.getLevel(), table.ident(keyword_)) - obs.valueByPressureLevel(decoder.getLevel2(), table.ident(keyword_))) ;     
                  //Log::dev()<< " : get --->" << val << endl;
            }
      }
protected:
      string keyword_;
};


class BufrMultiLevelAccessor : public BufrAccessor
{
public:
      BufrMultiLevelAccessor(const string& descriptor, const string& surface, const string& altitude) :
            surface_(surface), altitude_(altitude) { descriptor_ = descriptor; }
      BufrMultiLevelAccessor(const string& descriptor) :
                  surface_(descriptor), altitude_(descriptor) { descriptor_ = descriptor; }
      virtual ~BufrMultiLevelAccessor() {}
      void operator()(const ObsDecoder& decoder, MvObs& obs, double& val) const {
            const BufrIdentifiers& table =  BufrIdentTable::get(obs.originatingCentre()); 
            int type  = obs.messageType();
            
            map<int, bool>::const_iterator multilevel = multilevels_.find(type);
            if ( multilevel == multilevels_.end() ) {
                  Log::warning() << "Unknown observation type [" << val << "]: Contact graphics group\n";
                  val =  kBufrMissingValue;
            }
            if ( type == 0 || type == 1) {
                  // surface data 
                  val = obs.value(table.ident(surface_));               
            }
            else {
                  //Log::dev()<< " look for --->" << table.ident(altitude_) << " at " << decoder.getLevel();
                  // upper-air data 
                  if  ( type == 5 || type == 4 ) 
                        val = obs.value(table.ident(altitude_));              
                  else // Multi-level data
                        val = obs.valueByPressureLevel(decoder.getLevel(), table.ident(altitude_));   
                  //Log::dev()<< " : get --->" << val << endl;
            }
      }
protected:
      string surface_;
      string altitude_;
};

class BufrGeopotentialAccessor : public BufrMultiLevelAccessor 
{
public:
      BufrGeopotentialAccessor() : BufrMultiLevelAccessor("geopotential") {}
      virtual ~BufrGeopotentialAccessor() {}
};

class BufrTemperatureAccessor : public BufrMultiLevelAccessor 
{
public:
      BufrTemperatureAccessor() : BufrMultiLevelAccessor("temperature", "temperature_2meters",  "temperature") {}
      virtual ~BufrTemperatureAccessor() {}
};
class BufrDewPointAccessor : public BufrMultiLevelAccessor 
{
public:
      BufrDewPointAccessor() : BufrMultiLevelAccessor("dewpoint", "dewpoint_2meters",  "dewpoint") {}
      virtual ~BufrDewPointAccessor() {}
};

class BufrWindSpeedAccessor : public BufrMultiLevelAccessor 
{
public:
      BufrWindSpeedAccessor() : BufrMultiLevelAccessor("wind_speed", "wind_speed_10meters",  "wind_speed") {}
      virtual ~BufrWindSpeedAccessor() {}
};

class BufrWindDirectionAccessor : public BufrMultiLevelAccessor 
{
public:
      BufrWindDirectionAccessor() : BufrMultiLevelAccessor("wind_direction", "wind_direction_10meters",  "wind_direction") {}
      virtual ~BufrWindDirectionAccessor() {}
};

SimpleObjectMaker<BufrTypeAccessor, BufrAccessor> type_accessor("type");
SimpleObjectMaker<BufrTemperatureAccessor, BufrAccessor> temperature_accessor("temperature");
SimpleObjectMaker<BufrDewPointAccessor, BufrAccessor> dewpoint_accessor("dewpoint");
SimpleObjectMaker<BufrGeopotentialAccessor, BufrAccessor> geopotential_accessor("geopotential");
SimpleObjectMaker<BufrThicknessAccessor, BufrAccessor> thickness_accessor("thickness");

SimpleObjectMaker<BufrWindSpeedAccessor, BufrAccessor> wind_speed_accessor("wind_speed");
SimpleObjectMaker<BufrWindDirectionAccessor, BufrAccessor> wind_direction_accessor("wind_direction");

SimpleObjectMaker<BufrLowCloudAccessor, BufrAccessor> low_cloud_accessor("low_cloud");
SimpleObjectMaker<BufrMediumCloudAccessor, BufrAccessor> medium_cloud_accessor("medium_cloud");
SimpleObjectMaker<BufrHighCloudAccessor, BufrAccessor> high_cloud_accessor("high_cloud");

SimpleObjectMaker<BufrLowCloudHeightAccessor, BufrAccessor> low_cloud_nebulosity_accessor("low_cloud_height");
SimpleObjectMaker<BufrMediumCloudHeightAccessor, BufrAccessor> medium_cloud_nebulosity_accessor("medium_cloud_height");
SimpleObjectMaker<BufrHighCloudHeightAccessor, BufrAccessor> high_cloud_nebulosity_accessor("high_cloud_height");

SimpleObjectMaker<BufrLowCloudNebulosityAccessor, BufrAccessor> low_cloud_height_accessor("low_cloud_nebulosity");
SimpleObjectMaker<BufrMediumCloudNebulosityAccessor, BufrAccessor> medium_cloud_height_accessor("medium_cloud_nebulosity");
SimpleObjectMaker<BufrHighCloudNebulosityAccessor, BufrAccessor> high_cloud_heightaccessor("high_cloud_nebulosity");



/*!
 
\class ObsDecoder

*/
ObsDecoder::ObsDecoder()
{
      if ( multilevels_.empty() )
      {
            multilevels_[0] = false;
            multilevels_[1] = false;
            multilevels_[2] = true;
            multilevels_[3] = true;
            multilevels_[4] = false;
            multilevels_[5] = false;
      }
}

ObsDecoder::~ObsDecoder() 
{
}

void ObsDecoder::decode() 
{
      // Read observation file
      MvObsSet obsSet(file_name_.c_str());

      // Test 
        MvObsSetIterator obsIterator(obsSet);
      MvObs obs = obsIterator();

      ostringstream title;
      title << "Observation: " << obs.obsTime() << " [ type = " << obs.messageType() << " , subtype = " << obs.messageSubtype() << "]";
      title_ = title.str();
        
      while (obs)
      {
            MvLocation location = obs.location();
#ifdef OBS_DEBUG_
            cout << obsIterator.msgNumber() << " " << location << " " << time << " " << obs.WmoIdentNumber()
                 << "\n\ttype=" << obs.messageType() << ", subtype=" << obs.messageSubtype()
                 << "\n\tlat=" << obs.value(5001) << ", lon=" << obs.value(6001)
                 << "\n\ttempe=" << obs.value(12004)  << endl;
#endif
            push_back(GeoPoint(location.x(), location.y(), obs.value(12004)));
            obs = obsIterator();
      }
}


bool ObsDecoder::findInTypes(const string& val) 
{
      if ( types_.empty() ) 
      // No restriction specified by the user --> we accept all the types.
            return true;
      
      for (stringarray::const_iterator type = types_.begin(); type != types_.end(); ++type) {
            return magCompare(*type, val);      
      }
      return false;
}


bool ObsDecoder::checkLevel(double level) 
{
    if ( level == kBufrMissingValue ) return true; // Surface data ....
      level = level/100; // levels are expressed in Pascal...
    return (level_ - tolerance_ < level && level < level_ + tolerance_ );
}


void ObsDecoder::getInfo(const std::set<string>& tokens, multimap<string, string>& values)
{
      MvObsSet obsSet(file_name_.c_str());

      for (std::set<string>::const_iterator token = tokens.begin(); token != tokens.end(); ++token)
      {
            std::set<string> noduplicate;
            string value;
            try {
                  BufrAccessor* accessor = SimpleObjectMaker<BufrAccessor>::create(*token);
                  MvObsSetIterator   obsIterator(obsSet);
                  MvObs obs;
                  accessor->print();
                  obs = obsIterator();        
                  while (obs)
                  {
                        (*accessor)(*this, obs, value);
                        if ( value != "") noduplicate.insert(value);
                        obs = obsIterator();
                  }
            }
            catch ( NoFactoryException) {}

            for (std::set<string>::const_iterator no = noduplicate.begin(); no != noduplicate.end(); ++no)
            {
#ifdef OBS_DEBUG_
                  cout << " ObsDecoderToken: "<< *token << " -> " << *no << "\n";
#endif
                  // here we add it to the list, only if it is the obs_types_list!
                  if ( findInTypes(*no) ) values.insert(make_pair(*token, *no));
            }
      }
}

/*!
   For Each points ... create an CustomisedPoints using the tokens..
   then add it to the list! 
*/
00628 void ObsDecoder::customisedPoints(const Transformation& transformation, const std::set<string>& tokens, CustomisedPointsList& values)
{
      MvObsSet obsSet(file_name_.c_str());

      // Test 
      MvObsSetIterator obsIterator(obsSet);
      MvObs obs = obsIterator();

      ostringstream title;
      title << "Observation: " << obs.obsTime() << " [ type = " << obs.messageType() << " , subtype = " << obs.messageSubtype() << "]";
      title_ = title.str();
#ifdef OBS_DEBUG_
      cout << " ObsTitle: "<< title_ << "\n";
#endif 
      VectorOfPointers<vector<BufrAccessor*> > accessors;
      for (std::set<string>::const_iterator token = tokens.begin(); token != tokens.end(); ++token)
      {
#ifdef MAGICS_EXCEPTION
            try {
                  BufrAccessor* accessor = SimpleObjectMaker<BufrAccessor>::create(*token);
                  accessors.push_back(accessor);
            }
            catch (NoFactoryException&) {
                  Log::dev() << " No Specialised accessor for " << *token << ": take the default\n";
                  BufrAccessor* accessor = new BufrAccessor(*token);
                  accessors.push_back(accessor);
            }
#else 
            BufrAccessor* accessor = SimpleObjectMaker<BufrAccessor>::create(*token);
            if (accessor)
            {
                  accessors.push_back(accessor);
            }
            else
            {
                  Log::dev() << " No Specialised accessor for " << *token << ": take the default\n";
                  BufrAccessor* accessor = new BufrAccessor(*token);
                  accessors.push_back(accessor);
            }
#endif
      }

      // Create the type accessor!
      BufrAccessor* type_accessor = SimpleObjectMaker<BufrAccessor>::create("type");
      string type;
      BufrAccessor* level_accessor = new BufrAccessor("pressure");
      double level;
    
      while (obs)
      {
            (*type_accessor)(*this, obs, type);
            
            if ( findInTypes(type) )
            {
                  (*level_accessor)(*this, obs, level);
                  if (multilevels_[obs.messageType()] == false || checkLevel(level) )
                  {                       
                     MvLocation loc = obs.location();
                     double value;        
                     GeoPoint geo(loc.x(), loc.y());
                     if ( transformation.in(geo) ) {
                           CustomisedPoint* point = new CustomisedPoint(loc.x(), loc.y(),  obs.findSomeIdent());        
                        
                           for (VectorOfPointers<vector<BufrAccessor*> >::const_iterator accessor = accessors.begin(); accessor != accessors.end(); ++accessor)
                           {
                                 (**accessor)(*this, obs, value);     
                                 if (value != kBufrMissingValue)  (*point)[(*accessor)->keyword()] = value;
                           }
                           if ( !(*point).empty())
                           {
                                 TMetTime date = obs.msgTime();
                                 (*point)["time"] = date.GetHour()*100 + date.GetMin();
                                    values.push_back(point);
                           }
                           else 
                                 delete point;
                     }
                  }
            }
            obs = obsIterator();
      }
} 

PointsHandler<GeoPoint>& ObsDecoder::points()
{
      decode();

      pointsHandlers_.push_back(new PointsHandler<GeoPoint>(*this));
      return *(pointsHandlers_.back());
}

/*!
 Class information are given to the output-stream.
*/          
00722 void ObsDecoder::print(ostream& out)  const
{
      out << "ObsDecoder[";
      out << "]";
}

void ObsDecoder::visit(TitleNode& /*title*/)
{
      //title.add(new TitleEntry(title_));
}

Generated by  Doxygen 1.6.0   Back to index