CoverageSeriesExamples.cpp

Go to the documentation of this file.

// Examples
#include "STExamples.h"

//TerraLib
#include <terralib/common.h>
#include <terralib/dataaccess.h>
#include <terralib/geometry.h>
#include <terralib/datatype.h>
#include <terralib/st.h>

// STL
#include <iostream>

void CoverageSeriesExamples()
{
  try
  {  
    //Load coverage series from a set of geotif files
    std::unique_ptr<te::st::CoverageSeries> cvs = LoadCoverageSeriesFromGeotif();

    //Print coverage series info
    PrintCoverageSeriesInfo(cvs.get());

    //Extract time series from coverage series     
    //Load the centroide of the Angra city -> from a shapefile with the Angra city
    const std::string connInfo("file://" TERRALIB_DATA_DIR "/st/coverage/angra_city.shp");
    std::unique_ptr<te::da::DataSource> ds(te::da::DataSourceFactory::make("OGR", connInfo));
    ds->open();
    
    std::unique_ptr<te::da::DataSet> dSet = ds->getDataSet("angra_city"); 
    std::unique_ptr<te::gm::Geometry> geom; 
    
    if(dSet->moveNext())
    {
      std::size_t geomPos = te::da::GetFirstPropertyPos(dSet.get(), te::dt::GEOMETRY_TYPE);
      geom = dSet->getGeometry(geomPos);
      
      //Centroid operation is not implemented! 
      //We will use the center point of the MBR of the Angra geometry  
      const te::gm::Envelope* env = geom->getMBR(); 
      te::gm::Coord2D coord = env->getCenter();
      te::gm::Point point(coord.x, coord.y);

      std::unique_ptr<te::st::TimeSeries> result = cvs->getTimeSeries(point);

      //Print the result time series
      PrintTimeSeries(result.get());
    }
  }
  catch(const std::exception& e)
  {
    std::cout << std::endl << "An exception has occurred in CoverageExamples: " << e.what() << std::endl;
  }
  catch(...)
  {
    std::cout << std::endl << "An unexpected exception has occurred in CoverageExamples!" << std::endl;
  } 
}