ossimSrtmFactory.cpp

Go to the documentation of this file.

//----------------------------------------------------------------------------
//
// License:  See top level LICENSE.txt file.
//
// Author:  David Burken
//
// Description:
// 
// Shuttle Radar Topography Mission (SRTM) factory to return an
// ossimSrtmElevSource given a ground point.
//
//----------------------------------------------------------------------------
// $Id: ossimSrtmFactory.cpp 16123 2009-12-17 22:07:31Z dburken $

#include <cstdlib> /* abs() */
#include <iostream>
#include <iomanip>
#include <sstream>
#include <ossim/elevation/ossimSrtmFactory.h>
#include <ossim/elevation/ossimSrtmHandler.h>
#include <ossim/base/ossimTrace.h>
#include <ossim/base/ossimGpt.h>
#include <ossim/elevation/ossimElevManager.h>
#include <ossim/base/ossimNotifyContext.h>
#include <ossim/base/ossimStreamFactoryRegistry.h>
#include <ossim/base/ossimDirectory.h>
#include <ossim/support_data/ossimSrtmFilename.h>
#include <ossim/support_data/ossimSrtmSupportData.h>

static ossimTrace traceDebug ("ossimSrtmFactory:debug");

RTTI_DEF1(ossimSrtmFactory, "ossimSrtmFactory", ossimElevSourceFactory)

ossimSrtmFactory::ossimSrtmFactory()
   : ossimElevSourceFactory()
{}

ossimSrtmFactory::ossimSrtmFactory(const ossimFilename& dir)
   : ossimElevSourceFactory()
{
   theDirectory = dir;
}
      
ossimSrtmFactory::~ossimSrtmFactory()
{}

ossimElevSource* ossimSrtmFactory::getNewElevSource(const ossimGpt& gpt) const
{
   
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "DEBUG ossimSrtmFactory::getNewElevSource: Entered..."
         << std::endl;
   }

   ossimRefPtr<ossimElevCellHandler> srtmPtr;

   if (theDirectory == ossimFilename::NIL)
   {
      ossimNotify(ossimNotifyLevel_FATAL)
         << "ossimSrtmFactory::getNewElevSource: "
         << "SRTM directory has not been set!"
         << "\nReturning null elevation source..."
         << std::endl;
      
      return srtmPtr.release();
   }

   //---
   // Build up a srtm file name.
   //
   // Standard for name is upper case 'N' and 'W' lower case "hgt" like:
   // N27W081.hgt
   //---
   ossimFilename srtmFileBasename;

   int ilat =  static_cast<int>(floor(gpt.latd()));
   if (ilat < 0)
   {
      srtmFileBasename = "S";
   }
   else
   {
      srtmFileBasename = "N";
   }

   ilat = abs(ilat);
   std::ostringstream  os1;
   
   os1 << std::setfill('0') << std::setw(2) <<ilat;
   
   srtmFileBasename += os1.str().c_str();

   int ilon = static_cast<int>(floor(gpt.lond()));
   
   if (ilon < 0)
   {
      srtmFileBasename += "W";
   }
   else
   {
      srtmFileBasename += "E";
   }

   ilon = abs(ilon);
   std::ostringstream  os2;
   os2 << std::setfill('0') << std::setw(3) << ilon;
   
   srtmFileBasename += os2.str().c_str();
   srtmFileBasename.setExtension(".hgt");

   ossimFilename srtmFile = theDirectory.dirCat(srtmFileBasename);

   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "DEBUG ossimSrtmFactory::getNewElevSource:"
         << "\nSearching for file:  " << srtmFile
         << std::endl;
   }
   // ossimRefPtr<ossimIFStream> is = ossimStreamFactoryRegistry::instance()->
   std::shared_ptr<ossim::istream> is = ossim::StreamFactoryRegistry::instance()->
                                           createIstream(srtmFile);

   // Look for the file mix case, then all lower case, then all upper case.
   if ( is )
   {
      if(is->fail())
      {
         // Try down casing the whole thing.
         srtmFileBasename = srtmFileBasename.downcase();
         srtmFile = theDirectory.dirCat(srtmFileBasename);
         
         is = ossim::StreamFactoryRegistry::instance()->
               createIstream(srtmFile);
         if ( is )
         {
            if(is->fail())
            {
               // OK, try upcasing the whole thing.
               srtmFileBasename = srtmFileBasename.upcase();
               srtmFile = theDirectory.dirCat(srtmFileBasename);
               is = ossim::StreamFactoryRegistry::instance()->
                     createIstream(srtmFile);
            }
         }
      }
   }

   if ( is && (!is->fail()) )
   {
      is.reset();
      srtmPtr = new ossimSrtmHandler();
      if(srtmPtr->open(srtmFile)&&srtmPtr->pointHasCoverage(gpt) )
      {
         return srtmPtr.release();
      }
      else
      {
         srtmPtr = 0;
      }
   }
   return srtmPtr.release();
}