ossimDtedFactory.cpp

Go to the documentation of this file.

//**************************************************************************
// FILE: ossimDtedFactory.cc
//
// License:  See top level LICENSE.txt file.
//
// Author:  David Burken
//
// Description:
//
// Class definition for ossimDtedFactory.
//
// This is the base class interface for elevation source factories.  Contains
// pure virtual methods that all elevation source factories must implement.
//
//**************************************************************************
// $Id: ossimDtedFactory.cpp 16104 2009-12-17 18:09:59Z gpotts $

#include <cstdlib> /* abs() */
#include <sstream>
#include <iomanip>
#include <ossim/elevation/ossimDtedFactory.h>
#include <ossim/elevation/ossimDtedHandler.h>
#include <ossim/base/ossimTrace.h>
#include <ossim/base/ossimGpt.h>
#include <ossim/elevation/ossimElevManager.h>
#include <ossim/base/ossimNotifyContext.h>

static ossimTrace traceDebug ("ossimDtedFactory:debug");

RTTI_DEF1(ossimDtedFactory, "ossimDtedFactory", ossimElevSourceFactory)

ossimDtedFactory::ossimDtedFactory()
   : ossimElevSourceFactory()
{}

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

ossimDtedFactory::~ossimDtedFactory()
{}

ossimElevSource* ossimDtedFactory::getNewElevSource(const ossimGpt& gpt) const
{
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "DEBUG ossimDtedFactory::getNewElevSource: Entered..."
         << std::endl;
   }
   
   ossimRefPtr<ossimElevSource> dted_ptr;
   
   if (theDirectory == ossimFilename::NIL)
   {
      ossimNotify(ossimNotifyLevel_FATAL)
         << "FATAL ossimDtedFactory::getNewElevSource: "
         << "DTED directory has not been set!"
         << "\nReturning null elevation source..."
         << std::endl;
      return dted_ptr.release();
   }
   
   // Build up a dted file name.
   ossimFilename dted_file_base;
   
   int ilon = static_cast<int>(floor(gpt.lond()));
   
   if (ilon < 0)
   {
      dted_file_base = "w";
   }
   else
   {
      dted_file_base = "e";
   }
   
   ilon = abs(ilon);
   std::ostringstream  s1;
   s1 << std::setfill('0') << std::setw(3)<< ilon;
   
   dted_file_base += s1.str().c_str();//ossimString::toString(ilon);
   dted_file_base += "/";
   
   int ilat =  static_cast<int>(floor(gpt.latd()));
   if (ilat < 0)
   {
      dted_file_base += "s";
   }
   else
   {
      dted_file_base += "n";
   }

   ilat = abs(ilat);
   std::ostringstream  s2;

   s2<< std::setfill('0') << std::setw(2)<< ilat;
   
   dted_file_base += s2.str().c_str();

   // Look for a dted file with a level 2 extension first.
   ossimFilename dted_name = theDirectory.dirCat(dted_file_base);
   dted_name += ".dt2";

   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "DEBUG ossimDtedFactory::getNewElevSource:"
         << "\nSearching for file:  " << dted_name
         << std::endl;
   }
   
   if (dted_name.exists())
   {
      dted_ptr = new ossimDtedHandler(dted_name);
      if ( (!(dted_ptr->getErrorStatus())) &&
           dted_ptr->pointHasCoverage(gpt) )
      {
         return dted_ptr.release();
      }
      else
      {
         dted_ptr = 0;
      }
   }

   // Look for a dted file with a level 1 extension next.
   dted_name = theDirectory.dirCat(dted_file_base);
   dted_name += ".dt1";

   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "DEBUG ossimDtedFactory::getNewElevSource:"
         << "\nSearching for file:  " << dted_name
         << std::endl;
   }
   
   if (dted_name.exists() )
   {
      dted_ptr = new ossimDtedHandler(dted_name);
      if ( (!(dted_ptr->getErrorStatus())) &&
           dted_ptr->pointHasCoverage(gpt) )
      {
         return dted_ptr.release();
      }
      else
      {
         dted_ptr = 0;
      }
   }

   // Look for a dted file with a level 0 extension next.
   dted_name = theDirectory.dirCat(dted_file_base);
   dted_name += ".dt0";

   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "DEBUG ossimDtedFactory::getNewElevSource:"
         << "\nSearching for file:  " << dted_name
         << std::endl;
   }
   
   if (dted_name.exists() )
   {
      dted_ptr = new ossimDtedHandler(dted_name);
      if ( (!(dted_ptr->getErrorStatus())) &&
           dted_ptr->pointHasCoverage(gpt) )
      {
         return dted_ptr.release();
      }
      else
      {
         dted_ptr = 0;
      }
   }

   return dted_ptr.release();
}

void ossimDtedFactory::createIndex()
{
}