ossimViewshedTool Class Reference

#include <ossimViewshedTool.h>

Inheritance diagram for ossimViewshedTool:

Classes
class	Radial

Public Member Functions
	ossimViewshedTool ()
	~ossimViewshedTool ()
virtual void	setUsage (ossimArgumentParser &ap)
	Initializes the aurgument parser with expected parameters and options. More...
virtual bool	initialize (ossimArgumentParser &ap)
	Initializes from command line arguments. More...
virtual void	initialize (const ossimKeywordlist &kwl)
	Reads processing params from KWL and prepares for execute. More...
virtual ossimRefPtr< ossimImageData >	getChip (const ossimIrect &img_rect)
	Computes the viewshed for the area represented by argument AOI. More...
virtual bool	execute ()
	Writes product to output file. More...
virtual void	clear ()
	Disconnects and clears the DEM and image layers. More...
virtual ossimString	getClassName () const
void	test ()
	For engineering/debug. More...
Public Member Functions inherited from ossimChipProcTool
	ossimChipProcTool ()
	default constructor More...
virtual	~ossimChipProcTool ()
	virtual destructor More...
virtual void	abort ()
virtual bool	isChipProcessor () const
	Overrides base class implementation to indicate this class supports getChip() calls. More...
virtual ossimListenerManager *	getManager ()
virtual ossimObject *	getObject ()
virtual const ossimObject *	getObject () const
virtual ossimRefPtr< ossimImageData >	getChip (const ossimGrect &gnd_rect)
ossimRefPtr< ossimImageData >	getChip (const ossimDrect &map_bounding_rect, const ossimDpt &gsd)
ossimRefPtr< ossimImageData >	getChip ()
	Get chip method that assumes pre-initialized state. More...
const ossimFilename &	getProductFilename () const
Public Member Functions inherited from ossimTool
	ossimTool ()
virtual	~ossimTool ()
virtual void	loadJSON (const Json::Value &json_request)
	Reads processing params from JSON object provided. More...
virtual void	saveJSON (Json::Value &json) const
	Fetch product as JSON object when applicable Always returns true since using exception on error. More...
virtual const std::string &	getResponse () const
	Fetch product as string (typically JSON) when applicable, otherwise passes back empty string. More...
virtual void	getKwlTemplate (ossimKeywordlist &kwl)
	Assigns a template keywordlist to string for initializing derived classes. More...
void	getAPI (std::string &out) const
	Outputs a JSON representation of the Utility's API. More...
std::string	getAPI () const
void	getBuildDate (std::string &s) const
	Gets build date. More...
void	getRevision (std::string &s) const
	Gets revision. More...
void	getVersion (std::string &s) const
	Gets version. More...
void	setOutputStream (std::ostream *os)
	Redirects any console output to the supplied stream for logging or JNI application. More...
bool	helpRequested () const
	Returns true when the initialization detects a "--help" option, so caller can avoid subsequent execute() call. More...
Public Member Functions inherited from ossimObject
	ossimObject ()
virtual	~ossimObject ()
virtual ossimObject *	dup () const
virtual ossimString	getShortName () const
virtual ossimString	getLongName () const
virtual ossimString	getDescription () const
virtual RTTItypeid	getType () const
virtual bool	canCastTo (ossimObject *obj) const
virtual bool	canCastTo (const RTTItypeid &id) const
virtual bool	canCastTo (const ossimString &parentClassName) const
virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const
virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)
virtual std::ostream &	print (std::ostream &out) const
	Generic print method. More...
virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const
virtual void	accept (ossimVisitor &visitor)
Public Member Functions inherited from ossimReferenced
	ossimReferenced ()
	ossimReferenced (const ossimReferenced &)
ossimReferenced &	operator= (const ossimReferenced &)
void	ref () const
	increment the reference count by one, indicating that this object has another pointer which is referencing it. More...
void	unref () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...
void	unref_nodelete () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...
int	referenceCount () const
Public Member Functions inherited from ossim::JsonInterface
	JsonInterface ()
virtual	~JsonInterface ()
Public Member Functions inherited from ossimProcessInterface
	ossimProcessInterface ()
virtual	~ossimProcessInterface ()
bool	needsAborting () const
bool	isAbortRequested () const
bool	isAborted () const
bool	isExecuting () const
ossimProcessStatus	getProcessStatus () const
virtual void	setProcessStatus (ossimProcessStatus processStatus)
virtual double	getPercentComplete () const
virtual void	setPercentComplete (double percentComplete)
void	enableEvents ()
void	disableEvents ()
void	setCurrentMessage (const ossimString &message)
virtual std::ostream &	print (std::ostream &out) const
Public Member Functions inherited from ossimListenerManager
	ossimListenerManager ()
virtual	~ossimListenerManager ()
virtual void	fireEvent (ossimEvent &event)
virtual bool	addListener (ossimListener *listener)
virtual bool	removeListener (ossimListener *listener)
virtual bool	findListener (ossimListener *listener)

Static Public Attributes
static const char *	DESCRIPTION
	Used by ossimUtilityFactory. More...

Protected Member Functions
virtual void	initProcessingChain ()
	Derived classes initialize their custom chains here. More...
virtual void	initializeProjectionGsd ()
	Initializes the projection gsd. More...
virtual void	initializeAOI ()
	Initializes m_aoiViewRect with the output area of interest as specified in master KWL. More...
void	paintReticle ()
void	initRadials ()
bool	writeHorizonProfile ()
void	computeRadius ()
bool	optimizeFOV ()
bool	computeViewshed ()
Protected Member Functions inherited from ossimChipProcTool
void	processRemainingArgs (ossimArgumentParser &ap)
	Intended to be called after derived class has picked off its own options from the parser, and arguments remain (such as input and output filenames). More...
virtual void	finalizeChain ()
	Called after initProcessingChain() to append common items to the processing chain. More...
virtual void	loadImageFiles ()
	Creates chains for image entries associated with specified keyword. More...
virtual void	loadDemFiles ()
	Loads all DEM files specified in master KWL into the elev manager's database. More...
ossimRefPtr< ossimSingleImageChain >	createInputChain (const ossimFilename &image_file, ossim_uint32 entry_index=0)
	Creates the ossimSingleImageChain from image filename and populates the chain with resampler and product output projection on view-side of resampler's IVT. More...
void	createOutputProjection ()
	Creates the output or view projection. More...
ossimRefPtr< ossimMapProjection >	newIdentityProjection ()
	Sets the single image chain for identity operations view to a ossimImageViewAffineTransform. More...
ossimRefPtr< ossimMapProjection >	newUtmProjection ()
	Convenience method to get a utm projection. More...
bool	getProjectionOrigin (ossimGpt &gpt)
	Reads the KWL for origin latitude and central meridian. More...
ossimRefPtr< ossimImageFileWriter >	newWriter ()
	Creates a new writer. More...
void	propagateGeometryToChains ()
	Loops through all chains and sets the output projection. More...
ossimRefPtr< ossimImageSource >	combineLayers (std::vector< ossimRefPtr< ossimSingleImageChain > > &layers) const
	When multiple input sources are present, this method instantiates a combiner and adds inputs. More...
void	computeAdjustedViewFromGrect ()
	Initializes m_aoiViewRect given m_aoiGroundRect. More...
void	setAoiToInputs ()
	Assigns the AOI to be the bounding rect of the union of all inputs. More...
void	getBandList (ossim_uint32 image_idx, std::vector< ossim_uint32 > &bandList) const
	Gets the band list if BANDS keyword is set. More...
void	setReaderProps (ossimImageHandler *ih) const
	Passes reader properties to single image handler if any. More...
	ossimChipProcTool (const ossimChipProcTool &obj)
	Hidden from use copy constructor. More...
const ossimChipProcTool &	operator= (const ossimChipProcTool &)
	Hidden from use assignment operator. More...
ossimRefPtr< ossimGeoPolygon >	createClipPolygon () const
void	findCenterGpt (ossimGpt &gpt)
	Tries to determine the AOI center point based on KWL entries, else returns NaNs in gpt. More...
ossimRefPtr< ossimImageSource >	mosaicDemSources ()
	Some utilities need to work on DEMs as images. More...
Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()
Protected Member Functions inherited from ossimListenerManager
	ossimListenerManager (const ossimListenerManager &rhs)

Protected Attributes
ossimGpt	m_observerGpt
ossimDpt	m_observerVpt
double	m_obsHgtAbvTer
double	m_visRadius
Radial **	m_radials
bool	m_obsInsideAoi
bool	m_displayAsRadar
ossim_uint32	m_halfWindow
ossimRefPtr< ossimImageData >	m_outBuffer
ossimRefPtr< ossimMemoryImageSource >	m_memSource
ossim_uint8	m_visibleValue
ossim_uint8	m_hiddenValue
ossim_uint8	m_overlayValue
ossim_int32	m_reticleSize
bool	m_simulation
std::shared_ptr< ossimJobMultiThreadQueue >	m_jobMtQueue
ossim_uint32	m_numThreads
double	m_startFov
double	m_stopFov
bool	m_threadBySector
ossimFilename	m_horizonFile
std::map< double, double >	m_horizonMap
double	d_accumT
std::mutex	d_mutex
Protected Attributes inherited from ossimChipProcTool
ossimRefPtr< ossimImageGeometry >	m_geom
ossimIrect	m_aoiViewRect
ossimGrect	m_aoiGroundRect
std::vector< ossimRefPtr< ossimSingleImageChain > >	m_imgLayers
std::vector< ossimFilename >	m_demSources
ossimRefPtr< ossimImageFileWriter >	m_writer
ossimRefPtr< ossimImageChain >	m_procChain
ossimRefPtr< ossimRectangleCutFilter >	m_cutRectFilter
bool	m_projIsIdentity
ossimDpt	m_gsd
bool	m_geoScaled
ossimFilename	m_productFilename
ossimScalarType	m_productScalarType
bool	m_needCutRect
Protected Attributes inherited from ossimTool
ossimKeywordlist	m_kwl
std::ostream *	m_consoleStream
bool	m_helpRequested
std::string	m_response
Protected Attributes inherited from ossimProcessInterface
double	thePercentComplete
ossimProcessStatus	theProcessStatus
ossimString	theMessage
bool	theEventFlag
Protected Attributes inherited from ossimListenerManager
std::list< ossimListener * >	theListenerList
bool	theEnabledFlag
std::list< ossimListener * >	theDelayedAdd
std::list< ossimListener * >	theDelayedRemove
bool	theFireEventFlag

Friends
class	SectorProcessorJob
class	RadialProcessorJob
class	RadialProcessor

Additional Inherited Members
Public Types inherited from ossimProcessInterface
enum	ossimProcessStatus {   PROCESS_STATUS_UNKNOWN = 0, PROCESS_STATUS_EXECUTING = 1, PROCESS_STATUS_ABORTED = 2, PROCESS_STATUS_ABORT_REQUEST = 3,   PROCESS_STATUS_NOT_EXECUTING = 4 }

Detailed Description

Class for computing the viewshed on a DEM given the viewer location and max range of visibility

Definition at line 21 of file ossimViewshedTool.h.

Constructor & Destructor Documentation

ossimViewshedTool()

ossimViewshedTool::ossimViewshedTool

(

)

Definition at line 47 of file ossimViewshedTool.cpp.

References m_observerGpt, and ossimGpt::makeNan().

:   m_obsHgtAbvTer (1.5),
    m_visRadius (0.0),
    m_radials (0),
    m_obsInsideAoi (true),
    m_displayAsRadar (false),
    m_halfWindow (0),
    m_outBuffer (NULL),
    m_visibleValue (1),
    m_hiddenValue (128),
    m_overlayValue (255),
    m_reticleSize(25),
    m_simulation (false),
    m_numThreads(1),
    m_startFov(0),
    m_stopFov(0),
    m_threadBySector(false),
    d_accumT(0)
{
   m_observerGpt.makeNan();
}

~ossimViewshedTool()

ossimViewshedTool::~ossimViewshedTool

(

)

Definition at line 69 of file ossimViewshedTool.cpp.

References m_radials.

{
   if (m_radials)
   {
      for (int i=0; i<4; ++i)
         delete [] m_radials[i];
      delete [] m_radials;
   }
}

Member Function Documentation

clear()

void ossimViewshedTool::clear ( )

virtual

Disconnects and clears the DEM and image layers.

Leaves OSSIM initialized.

Reimplemented from ossimChipProcTool.

Definition at line 305 of file ossimViewshedTool.cpp.

References ossimChipProcTool::clear(), m_horizonMap, m_jobMtQueue, m_observerGpt, m_outBuffer, m_visRadius, and ossimGpt::makeNan().

{
   m_observerGpt.makeNan();
   m_visRadius = 0;
   m_outBuffer = 0;
   m_horizonMap.clear();
   m_jobMtQueue = 0;
   ossimChipProcTool::clear();
}

computeRadius()

void ossimViewshedTool::computeRadius ( )

protected

Definition at line 700 of file ossimViewshedTool.cpp.

References ossimGpt::distanceTo(), ossimIrect::hasNans(), ossimGrect::ll(), ossimGrect::lr(), ossimChipProcTool::m_aoiGroundRect, ossimChipProcTool::m_aoiViewRect, m_observerGpt, m_visRadius, ossimGrect::ul(), and ossimGrect::ur().

Referenced by initProcessingChain().

{
   ostringstream xmsg ("ossimViewshedUtil::computeRadius() -- ");

   // AOI is required for computing R
   if (m_aoiViewRect.hasNans())
   {
      xmsg<<"AOI undefined. Cannot compute visibility radius." << ends;
      throw ossimException(xmsg.str());
   }

   // Compute distance from observer to farthest corner of AOI. This is the radius
   m_visRadius= m_observerGpt.distanceTo(m_aoiGroundRect.ul());
   double d = m_observerGpt.distanceTo(m_aoiGroundRect.ur());
   if (d > m_visRadius)
      m_visRadius = d;
   d = m_observerGpt.distanceTo(m_aoiGroundRect.lr());
   if (d > m_visRadius)
      m_visRadius = d;
   d = m_observerGpt.distanceTo(m_aoiGroundRect.ll());
   if (d > m_visRadius)
      m_visRadius = d;
}

computeViewshed()

bool ossimViewshedTool::computeViewshed ( )

protected

Definition at line 486 of file ossimViewshedTool.cpp.

References ossimImageDataFactory::create(), d_accumT, ossimImageData::fill(), ossimImageChain::getNullPixelValue(), ossim::getNumberOfThreads(), ossimIrect::height(), ossimImageData::initialize(), initRadials(), ossimImageDataFactory::instance(), ossimChipProcTool::m_aoiViewRect, m_halfWindow, m_jobMtQueue, m_memSource, m_numThreads, m_outBuffer, ossimChipProcTool::m_procChain, m_radials, m_threadBySector, ossimProcessInterface::needsAborting(), optimizeFOV(), OSSIM_UINT8, ossimNotify(), ossimNotifyLevel_INFO, paintReticle(), ossimMemoryImageSource::setImage(), ossimImageData::setImageRectangle(), ossim::Thread::sleepInMicroSeconds(), ossimRefPtr< T >::valid(), and ossimIrect::width().

Referenced by execute(), and getChip().

{
   // Allocate the output image buffer:
   m_outBuffer = ossimImageDataFactory::instance()->create(0, OSSIM_UINT8, 1, m_aoiViewRect.width(),
                                                           m_aoiViewRect.height());

   ostringstream xmsg;
   if (!m_outBuffer.valid() || !m_memSource.valid())
   {
      xmsg<<"ossimViewshedUtil:"<<__LINE__<<"  Error encountered allocating output image buffer.";
      throw ossimException(xmsg.str());
   }

   // Initialize the image with all points hidden:
   m_outBuffer->initialize();
   m_outBuffer->setImageRectangle(m_aoiViewRect);
   m_outBuffer->fill(m_procChain->getNullPixelValue());
   m_memSource->setImage(m_outBuffer);

   // Initialize the radials after intersecting the requested FOV with the FOV required to see the
   // full AOI (not applicable if observer inside AOI). Skip radial init if no intersection found:
   if (!optimizeFOV())
      return false;
   initRadials();

   // The viewshed process necessarily first fills the output buffer with the complete result before
   // the writer requests a tile. Control is passed later to the base class execute() for writing.
   d_accumT = 0;

   if (m_numThreads == 0)
      m_numThreads = ossim::getNumberOfThreads();

   if (m_numThreads > 1)
   {
      std::shared_ptr<ossimJobQueue> jobQueue = std::make_shared<ossimJobQueue>();
      for (int sector=0; sector<8; ++sector)
      {
         if (m_radials[sector] == 0)
            continue;

         if (m_threadBySector)
         {
            std::shared_ptr<SectorProcessorJob> job = std::make_shared<SectorProcessorJob>(this, sector, m_halfWindow);
            jobQueue->add(job, false);
         }
         else
         {
            for (ossim_uint32 r=0; r<=m_halfWindow; ++r)
            {
               std::shared_ptr<RadialProcessorJob> job = std::make_shared<RadialProcessorJob>(this, sector, r, m_halfWindow);
               jobQueue->add(job, false);
            }
         }
         if (needsAborting())
            return 0;
      }

      ossimNotify(ossimNotifyLevel_INFO) << "\nSubmitting "<<jobQueue->size()<<" jobs..."<<endl;
      m_jobMtQueue = std::make_shared<ossimJobMultiThreadQueue>(jobQueue, m_numThreads);

      // Wait until all radials have been processed before proceeding:
      ossimNotify(ossimNotifyLevel_INFO) << "Waiting for job threads to finish..."<<endl;
      while (m_jobMtQueue->hasJobsToProcess() || m_jobMtQueue->numberOfBusyThreads())
         ossim::Thread::sleepInMicroSeconds(250);
   }
   else
   {
      // Unthreaded processing:
      ossimNotify(ossimNotifyLevel_INFO) << "\nProcessing radials (non-threaded)..."<<endl;

      // Loop over pixels in layer for each sector:
      for (int sector=0; sector<8; ++sector)
      {
         if (m_radials[sector] == 0)
            continue;

         std::shared_ptr<SectorProcessorJob> spj = std::make_shared<SectorProcessorJob>(this, sector, m_halfWindow);
         spj->start();

         if (needsAborting())
            return false;

      } // end loop over sectors
   }

   ossimNotify(ossimNotifyLevel_INFO) << "Finished processing radials."<<endl;
   paintReticle();

   return true;
}

execute()

bool ossimViewshedTool::execute ( )

virtual

Writes product to output file.

Returns true if successful.

Note

Throws ossimException on error.

Reimplemented from ossimChipProcTool.

Definition at line 472 of file ossimViewshedTool.cpp.

References computeViewshed(), ossimString::empty(), ossimChipProcTool::execute(), ossimTool::m_helpRequested, m_horizonFile, ossimNotify(), ossimNotifyLevel_INFO, and writeHorizonProfile().

{
   if (m_helpRequested)
      return true;

   if (!computeViewshed())
      return false;

   if (!m_horizonFile.empty() && writeHorizonProfile())
      ossimNotify(ossimNotifyLevel_INFO) << "Wrote horizon profile to <"<<m_horizonFile<<">" <<endl;

   return ossimChipProcTool::execute();
}

getChip()

ossimRefPtr< ossimImageData > ossimViewshedTool::getChip ( const ossimIrect &  img_rect )

virtual

Computes the viewshed for the area represented by argument AOI.

Reimplemented from ossimChipProcTool.

Definition at line 449 of file ossimViewshedTool.cpp.

References computeViewshed(), ossimImageChain::getTile(), ossimImageGeometry::localToWorld(), ossimChipProcTool::m_aoiGroundRect, ossimChipProcTool::m_aoiViewRect, ossimChipProcTool::m_geom, ossimChipProcTool::m_procChain, ossimImageGeometry::setImageSize(), ossimIrect::size(), and ossimRefPtr< T >::valid().

{
   cerr<<"ossimViewshedUtil:"<<__LINE__<<endl;//TODO:remove debug
   ostringstream xmsg;
   if (!m_geom.valid())
      return 0;

   cerr<<"ossimViewshedUtil:"<<__LINE__<<endl;//TODO:remove debug
   m_aoiViewRect = bounding_irect;
   m_geom->setImageSize( m_aoiViewRect.size() );
   m_geom->localToWorld(m_aoiViewRect, m_aoiGroundRect);

   cerr<<"ossimViewshedUtil:"<<__LINE__<<endl;//TODO:remove debug
   if (computeViewshed())
   {
      // The memory source has been populated, now do the getTile on the full chain to pick up
      // other filters inserted after the memsource:
      return m_procChain->getTile( m_aoiViewRect, 0 );
   }
   // else:
   return 0;
}

getClassName()

virtual ossimString ossimViewshedTool::getClassName ( ) const

inlinevirtual

Reimplemented from ossimTool.

Definition at line 58 of file ossimViewshedTool.h.

{ return "ossimViewshedUtil"; }

initialize() [1/2]

bool ossimViewshedTool::initialize ( ossimArgumentParser &  ap )

virtual

Initializes from command line arguments.

Returns

FALSE if –help option requested or no params provided, so that derived classes can

Note

Throws ossimException on error.

Reimplemented from ossimChipProcTool.

Definition at line 133 of file ossimViewshedTool.cpp.

References ossimKeywordlist::addPair(), ossimChipProcTool::initialize(), ossimTool::m_helpRequested, ossimTool::m_kwl, ossimArgumentParser::read(), and ossimString::toDouble().

{
   ostringstream xmsg;
   xmsg<<"ossimViewshedUtil::initialize(ossimArgumentParser) -- ";

   int numArgsExpected = 4;

   // Base class first:
   if (!ossimChipProcTool::initialize(ap))
      return false;
   if (m_helpRequested)
      return true;

   string ts1;
   ossimArgumentParser::ossimParameter sp1(ts1);
   string ts2;
   ossimArgumentParser::ossimParameter sp2(ts2);
   string ts3;
   ossimArgumentParser::ossimParameter sp3(ts3);

   if ( ap.read("--fov", sp1, sp2) )
   {
      double startFov = ossimString(ts1).toDouble();
      if (startFov < 0)
         startFov += 360.0;
      ostringstream value;
      value<<startFov<<" "<<ts2;
      m_kwl.addPair( FOV_KW, value.str() );
   }

   if ( ap.read("--hgt-of-eye", sp1) || ap.read("--height-of-eye", sp1) )
      m_kwl.addPair( HEIGHT_OF_EYE_KW, ts1 );

   if ( ap.read("--horizon", sp1) || ap.read("--horizon-file", sp1))
      m_kwl.addPair( HORIZON_FILE_KW, ts1 );

   if ( ap.read("--observer", sp1, sp2) )
   {
      ostringstream value;
      value<<ts1<<" "<<ts2;
      m_kwl.addPair( OBSERVER_KW, value.str() );
      numArgsExpected -= 2;
   }

   if ( ap.read("--radius", sp1) )
      m_kwl.addPair( VISIBILITY_RADIUS_KW, ts1 );

   if ( ap.read("--reticle", sp1) )
      m_kwl.addPair( RETICLE_SIZE_KW, ts1 );

   if ( ap.read("--values", sp1, sp2, sp3) || ap.read("--viewshed-coding", sp1, sp2, sp3))
   {
      ostringstream value;
      value<<ts1<<" "<<ts2<<" "<<ts3;
      m_kwl.addPair( VIEWSHED_CODING_KW, value.str() );
   }

   if ( ap.read("--threads", sp1) )
   {
      m_kwl.addPair( THREADS_KW, ts1 );
   }

   // The remaining options are available only via command line (i.e., no KWL entries defined)
   if ( ap.read("--tbs") )
      m_threadBySector = true;

   if ( ap.read("--simulation") )
      m_simulation = true;

   // Moved to: ossimViewshedTool::initialize(const ossimKeywordlist& kwl)
   // if ( ap.read("--threads", sp1) )
   //    m_numThreads = ossimString(ts1).toUInt32();

   if ( ap.argc() < numArgsExpected )
   {
      xmsg<<"Expecting more arguments.";
      throw(ossimException(xmsg.str()));
   }
   else
   {
      ossimString latstr = ap[1];
      ossimString lonstr = ap[2];
      ostringstream value;
      value<<latstr<<" "<<lonstr;
      m_kwl.addPair( OBSERVER_KW, value.str() );
      ap.remove(1,2);
      processRemainingArgs(ap);
   }

   return true;
}

initialize() [2/2]

void ossimViewshedTool::initialize ( const ossimKeywordlist &  kwl )

virtual

Reads processing params from KWL and prepares for execute.

Returns TRUE if successful.

Note

Throws ossimException on error.

Reimplemented from ossimChipProcTool.

Definition at line 225 of file ossimViewshedTool.cpp.

References ossimKeywordlist::findKey().

{
   ossimString value;

   value = kwl.findKey(THREADS_KW);
   {
      if ( value.size() )
      {
         m_numThreads = ossimString(value).toUInt32();
      }
   }
   
   value = kwl.findKey(FOV_KW);
   if (!value.empty())
   {
      vector <ossimString> coordstr;
      value.split(coordstr, ossimString(" ,"), false);
      if (coordstr.size() == 2)
      {
         m_startFov = coordstr[0].toDouble();
         m_stopFov = coordstr[1].toDouble();
         if (m_startFov < 0)
            m_startFov += 360.0;
      }
   }

   value = kwl.findKey(HEIGHT_OF_EYE_KW);
   if (!value.empty())
      m_obsHgtAbvTer = value.toDouble();

   m_horizonFile = kwl.findKey(HORIZON_FILE_KW);

   value = kwl.findKey(OBSERVER_KW);
   if (!value.empty())
   {
      vector <ossimString> coordstr;
      value.split(coordstr, ossimString(" ,"), false);
      if (coordstr.size() == 2)
      {
         m_observerGpt.lat = coordstr[0].toDouble();
         m_observerGpt.lon = coordstr[1].toDouble();
         m_observerGpt.hgt = 0.0;
      }
   }

   value = kwl.findKey(RETICLE_SIZE_KW);
   if (!value.empty())
      m_reticleSize = value.toInt32();

   value = kwl.findKey(VIEWSHED_CODING_KW);
   if (!value.empty())
   {
      vector <ossimString> coordstr;
      value.split(coordstr, ossimString(" ,"), false);
      if (coordstr.size() == 3)
      {
         m_visibleValue = coordstr[0].toUInt8();
         m_hiddenValue = coordstr[1].toUInt8();
         m_overlayValue = coordstr[2].toUInt8();
      }
   }

   value = kwl.findKey(VISIBILITY_RADIUS_KW);
   if (!value.empty())
   {
      m_visRadius = value.toDouble();
      m_displayAsRadar = true;
   }

   // If running simulation, clear out all pre-loaded elevation databases:
   if (m_simulation)
   {
      ossimElevManager::instance()->clear();
      ossimElevManager::instance()->setUseGeoidIfNullFlag(false);
   }

   // Base class does most work:
   ossimChipProcTool::initialize(kwl);
}

initializeAOI()

void ossimViewshedTool::initializeAOI ( )

protectedvirtual

Initializes m_aoiViewRect with the output area of interest as specified in master KWL.

Initialization will either come from user defined cut options or the source bounding rect with user options taking precidence.

Note

Throws ossimException on error.

Reimplemented from ossimChipProcTool.

Definition at line 336 of file ossimViewshedTool.cpp.

References ossimKeywordlist::hasKey(), ossimGrect::hasNans(), ossimGpt::hasNans(), ossimChipProcTool::initializeAOI(), ossimChipProcTool::m_aoiGroundRect, ossimTool::m_kwl, m_observerGpt, m_visRadius, and ossimGrect::midPoint().

{
   ossimChipProcTool::initializeAOI();
   if (!m_aoiGroundRect.hasNans())
   {
      // AOI established by base class, nothing to do except take this opportunity to set the
      // observer location to AOI center if not already defined:
      if (m_observerGpt.hasNans())
         m_observerGpt = m_aoiGroundRect.midPoint();
      return;
   }

   // Not enough info available to base class to determine AOI, maybe can determine from observer
   // position and radius:
   if ((m_visRadius == 0) && m_kwl.hasKey(AOI_SIZE_METERS_KW))
   {
      ossimString lookup = m_kwl.findKey( AOI_SIZE_METERS_KW );
      lookup.trim();
      m_visRadius = 0.5*(lookup.before(" ").toDouble() + lookup.after(" ").toDouble());
      m_displayAsRadar = true;
   }
   if (m_observerGpt.hasNans())
      findCenterGpt(m_observerGpt);

   if ((m_visRadius != 0) && !m_observerGpt.hasNans())
   {
      ossimMapProjection* proj = dynamic_cast<ossimMapProjection*>(m_geom->getProjection());
      if (!proj)
         return;

      ossimDpt metersPerDegree (m_observerGpt.metersPerDegree());
      double dlat = m_visRadius/metersPerDegree.y;
      double dlon = m_visRadius/metersPerDegree.x;
      ossimGpt ulg (m_observerGpt.lat + dlat, m_observerGpt.lon - dlon);
      ossimGpt lrg (m_observerGpt.lat - dlat, m_observerGpt.lon + dlon);

      m_aoiGroundRect = ossimGrect(ulg, lrg);
      proj->setUlTiePoints(ulg);

      computeAdjustedViewFromGrect();
   }
}

initializeProjectionGsd()

void ossimViewshedTool::initializeProjectionGsd ( )

protectedvirtual

Initializes the projection gsd.

This loops through all chains to find the minimum gsd.

Note

Throws ossimException on error.

Reimplemented from ossimChipProcTool.

Definition at line 315 of file ossimViewshedTool.cpp.

References ossimElevManager::getHeightAboveEllipsoid(), ossimElevManager::getMeanSpacingMeters(), ossimImageGeometry::getProjection(), ossimDpt::hasNans(), ossimGpt::hasNans(), ossimChipProcTool::initializeProjectionGsd(), ossimElevManager::instance(), ossimChipProcTool::m_geom, ossimChipProcTool::m_geoScaled, ossimChipProcTool::m_gsd, m_observerGpt, ossimMapProjection::setMetersPerPixel(), ossimMapProjection::setOrigin(), ossimDpt::x, and ossimDpt::y.

{
   // First try normal base class initialization. If that doesn't work, then probably no DEM
   // or input images were provided, so need to use elev manager resoltion at observer point.

   ossimChipProcTool::initializeProjectionGsd();
   if (!m_gsd.hasNans() || m_observerGpt.hasNans())
      return;

   ossimMapProjection* proj = dynamic_cast<ossimMapProjection*>(m_geom->getProjection());
   if (!proj)
      return;

   ossimElevManager* elevMgr = ossimElevManager::instance();
   elevMgr->getHeightAboveEllipsoid(m_observerGpt);
   m_gsd.x = m_gsd.y = ossimElevManager::instance()->getMeanSpacingMeters();
   if (m_geoScaled)
      proj->setOrigin(m_observerGpt);
   proj->setMetersPerPixel(m_gsd);
}

initProcessingChain()

void ossimViewshedTool::initProcessingChain ( )

protectedvirtual

Derived classes initialize their custom chains here.

Implements ossimChipProcTool.

Definition at line 379 of file ossimViewshedTool.cpp.

References ossimImageChain::add(), ossimIrect::area(), ossimIrect::clipToRect(), ossimChipProcTool::combineLayers(), computeRadius(), ossimConnectableObject::connectMyInputTo(), ossimRefPtr< T >::get(), ossimElevManager::getHeightAboveEllipsoid(), ossimImageGeometry::getProjection(), ossimGpt::hasNans(), ossimIrect::hasNans(), ossimIrect::height(), ossimGpt::hgt, ossimImageDataFactory::instance(), ossimElevManager::instance(), ossimImageGeometry::localToWorld(), ossimChipProcTool::m_aoiGroundRect, ossimChipProcTool::m_aoiViewRect, ossimChipProcTool::m_geom, ossimChipProcTool::m_gsd, m_halfWindow, ossimChipProcTool::m_imgLayers, m_memSource, m_observerGpt, m_observerVpt, m_obsHgtAbvTer, m_outBuffer, ossimChipProcTool::m_procChain, m_visRadius, OSSIM_UINT8, ossimMemoryImageSource::setImage(), ossimMemoryImageSource::setImageGeometry(), ossimImageData::setImageRectangle(), ossimImageGeometry::setImageSize(), ossimIrect::size(), size, ossimRefPtr< T >::valid(), ossimIrect::width(), ossimImageGeometry::worldToLocal(), and ossimDpt::x.

{
   ostringstream xmsg;

   if (m_observerGpt.hasNans())
   {
      xmsg<<"ossimViewshedUtil:"<<__LINE__<<" Observer ground position has not been set."<<ends;
      throw ossimException(xmsg.str());
   }

   ossimElevManager* elevMgr = ossimElevManager::instance();

   // Initialize the height of eye component of observer position:
   m_observerGpt.hgt = elevMgr->getHeightAboveEllipsoid(m_observerGpt);
   m_observerGpt.hgt += m_obsHgtAbvTer;
   m_geom->worldToLocal(m_observerGpt, m_observerVpt);

   ossimRefPtr<ossimMapProjection> mapProj =
         dynamic_cast<ossimMapProjection*>(m_geom->getProjection());

   // If no radius specified, need to compute R large enough to cover the requested AOI:
   if (m_visRadius == 0)
      computeRadius();
   if (m_halfWindow == 0)
      m_halfWindow = ossim::round<ossim_int32, double>(m_visRadius/m_gsd.x);
   ossim_uint32 size = 2*m_halfWindow + 1;

   // If no AOI defined, just use the visibility rectangle:
   ossimIrect visRect (ossimIpt(m_observerVpt), size, size);
   if (m_aoiViewRect.hasNans())
   {
      m_aoiViewRect = visRect;
      m_geom->localToWorld(ossimDrect(m_aoiViewRect), m_aoiGroundRect);
      m_geom->setImageSize(m_aoiViewRect.size());
   }

   // Allocate the output image buffer:
   ossimIrect bufViewRect = visRect.clipToRect(m_aoiViewRect);
   if (bufViewRect.area() == 0)
   {
      xmsg<<"ossimViewshedUtil:"<<__LINE__<<" The requested AOI rect is outside the visibility range." << ends;
      throw ossimException(xmsg.str());
   }
   m_outBuffer = ossimImageDataFactory::instance()->
         create(0, OSSIM_UINT8, 1, m_aoiViewRect.width(), m_aoiViewRect.height());
   if(!m_outBuffer.valid())
   {
      xmsg<<"ossimViewshedUtil:"<<__LINE__<<" Output buffer allocation failed." << ends;
      throw ossimException(xmsg.str());
   }
   m_outBuffer->setImageRectangle(m_aoiViewRect);

   // The processing chain for this class is simply a memory source containing the output buffer:
   m_memSource = new ossimMemoryImageSource;
   m_memSource->setImage(m_outBuffer);
   m_memSource->setImageGeometry(m_geom.get());

   // If input image(s) provided, need to combine them with the product:
   if (m_imgLayers.empty())
   {
      m_procChain->add(m_memSource.get());
   }
   else
   {
      ossimRefPtr<ossimImageSource> combiner = combineLayers(m_imgLayers);
      combiner->connectMyInputTo(m_memSource.get());
      m_procChain->add(combiner.get());
   }
}

initRadials()

void ossimViewshedTool::initRadials ( )

protected

Definition at line 724 of file ossimViewshedTool.cpp.

References ossimViewshedTool::Radial::azimuth, m_halfWindow, m_radials, m_startFov, and m_stopFov.

Referenced by computeViewshed().

{
   // All eaight sectors' radials have the same azimuths except that the abscissa and ordinate are
   // reversed between the N-S and E-W sectors, i.e., the N and S sectors use the y-axis as the
   // abscissa (u) and the x-axis is the ordinate (v). The azimuth (dv/du) is therefore DX/DY for
   // the north and south, while the azimuth is DY/DX for the east and west sectors. Nevertheless,
   // each sectors radials must be maintained separately as they contain the max elevation angle.

   // First determine which sectors are involved given the desired FOV:
   bool* sectorInFov = new bool[8];
   memset(sectorInFov, false, 8);
   bool crossed_north = true;
   if (m_stopFov <= m_startFov) // Crosses 0 azimuth
      crossed_north = false;
   double azimuth = m_startFov;
   for (int i=0; (i < 8) && ((azimuth < m_stopFov) || !crossed_north); ++i)
   {
      if ((azimuth >= 0) && (azimuth < 45.0))
         sectorInFov[0] = true;
      else if (azimuth < 90.0)
         sectorInFov[1] = true;
      else if (azimuth < 135.0)
         sectorInFov[2] = true;
      else if (azimuth < 180.0)
         sectorInFov[3] = true;
      else if (azimuth < 225.0)
         sectorInFov[4] = true;
      else if (azimuth < 270.0)
         sectorInFov[5] = true;
      else if (azimuth < 315.0)
         sectorInFov[6] = true;
      else if (azimuth < 360.0)
         sectorInFov[7] = true;

      azimuth += 45.0;
      if (azimuth > 360.0)
      {
         azimuth -= 360.0;
         crossed_north = true;
      }
   }

   // Compute the azimuth slopes for each radial in the sector.
   m_radials = new Radial* [8];
   double du = m_halfWindow;
   for (int sector=0; sector<8; ++sector)
   {
      if (!sectorInFov[sector])
      {
         m_radials[sector] = 0;
         continue;
      }

      ossim_uint32 ridx = 0;
      m_radials[sector] = new Radial [m_halfWindow+1];
      for (ossim_int32 dv = 0; dv <= (ossim_int32) m_halfWindow; ++dv)
      {
         if (sector & 1) // odd-numbered sector, azimuths computed in reverse order
            m_radials[sector][m_halfWindow-ridx].azimuth = ((double)dv)/du;
         else
            m_radials[sector][ridx].azimuth = ((double)dv)/du;
         ++ridx;
      }
   }

   // Cleanup:
   delete [] sectorInFov;
   sectorInFov = 0;
}

optimizeFOV()

bool ossimViewshedTool::optimizeFOV ( )

protected

Definition at line 577 of file ossimViewshedTool.cpp.

References ossimGpt::azimuthTo(), ossimGpt::lat, ossimGrect::ll(), ossimGpt::lon, ossimGrect::lr(), ossimChipProcTool::m_aoiGroundRect, m_observerGpt, m_startFov, m_stopFov, ossimNotify(), ossimNotifyLevel_INFO, ossimGrect::pointWithin(), ossimGrect::ul(), and ossimGrect::ur().

Referenced by computeViewshed(), and test().

{
   bool intersects = false;

   // If the observer position lies outside of the requested AOI, we can reduce the search arc:
   if (m_aoiGroundRect.pointWithin(m_observerGpt))
      return true;

   // Determine cardinal region (N, NE, E, ...) of observer relative to AOI:
   enum CardinalDirections { N=1, S=2, E=4, W=8, NE=5, NW=9, SE=6, SW=10 };
   int direction = 0;
   if (m_observerGpt.lat > m_aoiGroundRect.ul().lat)
      direction = (int) N;
   else if (m_observerGpt.lat < m_aoiGroundRect.ll().lat)
      direction = (int) S;
   if (m_observerGpt.lon < m_aoiGroundRect.ul().lon)
      direction += (int) W;
   else if (m_observerGpt.lon > m_aoiGroundRect.ur().lon)
      direction += (int) E;

   // Calculate start and stop FOV depending on region:
   double start, stop;
   switch ((CardinalDirections) direction)
   {
   case N:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.ur());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.ul());
      break;
   case NE:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.lr());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.ul());
      break;
   case E:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.lr());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.ur());
      break;
   case SE:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.ll());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.ur());
      break;
   case S:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.ll());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.lr());
      break;
   case SW:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.ul());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.lr());
      break;
   case W:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.ul());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.ll());
      break;
   case NW:
   default:
      start = m_observerGpt.azimuthTo(m_aoiGroundRect.ur());
      stop  = m_observerGpt.azimuthTo(m_aoiGroundRect.ll());
      break;
   }

   // Now need to intersect this arc with the requested FOV:
   if (m_startFov == m_stopFov)
   {
      // There was no requested FOV (i.e, FOV = 360). So use the optimized FOV straight away:
      m_startFov = start;
      m_stopFov  = stop;
      intersects = true;
   }
   else
   {
      // Pick m_startFov as reference, and make sure all others are greater:
      double a1 = m_stopFov;
      double a2 = start;
      double a3 = stop;
      if (m_startFov > m_stopFov)
         a1 += 360;
      if (m_startFov > start)
         a2 += 360;
      if (m_startFov > stop)
         a3 += 360;

      // Map to sort remaining azimuths by increasing angle clockwise:
      map<double, int> angle_map;
      angle_map.insert(pair<double, int>(a1, 1));
      angle_map.insert(pair<double, int>(a2, 2));
      angle_map.insert(pair<double, int>(a3, 3));

      map<double, int>::iterator iter = angle_map.begin();
      if (iter->second == 1)
      {
         ++iter;
         if (iter->second == 3)
            intersects = true;
      }
      else if (iter->second == 2)
      {
         m_startFov = start;
         intersects = true;
         ++iter;
         if (iter->second == 3)
            m_stopFov = stop;
      }
      else
      {
         intersects = true;
         m_stopFov = stop;
      }
   }

   if (!intersects)
   {
      ossimNotify(ossimNotifyLevel_INFO)<<
            "ossimViewshedUtil::optimizeFOV() -- No FOV intersection found. Nothing to do."<<endl;
   }
   else
   {
      ossimNotify(ossimNotifyLevel_INFO)<<"ossimViewshedUtil::optimizeFOV() -- "
            "The start and stop FOV azimuths have been optimized to "<<m_startFov<<" -> "
            <<m_stopFov<<" deg."<<endl;
   }

   return intersects;
}

paintReticle()

void ossimViewshedTool::paintReticle ( )

protected

Definition at line 794 of file ossimViewshedTool.cpp.

References ossimIrect::lr(), ossimChipProcTool::m_aoiViewRect, m_displayAsRadar, ossimChipProcTool::m_geom, m_observerGpt, m_observerVpt, m_outBuffer, m_overlayValue, m_reticleSize, ossimIrect::pointWithin(), ossimImageData::setValue(), ossimIrect::ul(), ossimImageGeometry::worldToLocal(), ossimIpt::x, ossimDpt::x, x, ossimIpt::y, ossimDpt::y, and y.

Referenced by computeViewshed().

{
   if (m_reticleSize == 0)
      return;

   // Highlight the observer position with X reticle:
   ossimDpt obsViewPt;
   m_geom->worldToLocal(m_observerGpt, obsViewPt);
   if (m_aoiViewRect.pointWithin(ossimIpt(obsViewPt)))
   {
      for (int i=-m_reticleSize; i<=m_reticleSize; ++i)
      {
         if (m_aoiViewRect.pointWithin(ossimIpt(m_observerVpt.x + i, m_observerVpt.y)))
            m_outBuffer->setValue(m_observerVpt.x + i, m_observerVpt.y    , m_overlayValue);
         if (m_aoiViewRect.pointWithin(ossimIpt(m_observerVpt.x, m_observerVpt.y + i)))
            m_outBuffer->setValue(m_observerVpt.x    , m_observerVpt.y + i, m_overlayValue);
      }
   }

   // Paint boundary rectangle if no visibility radius painted:
   if (!m_displayAsRadar)
   {
      for (int y=m_aoiViewRect.ul().y; y<=m_aoiViewRect.lr().y; y++)
      {
         m_outBuffer->setValue(m_aoiViewRect.ul().x, y, m_overlayValue);
         m_outBuffer->setValue(m_aoiViewRect.lr().x, y, m_overlayValue);
      }
      for (int x=m_aoiViewRect.ul().x; x<=m_aoiViewRect.lr().x; x++)
      {
         m_outBuffer->setValue(x, m_aoiViewRect.ul().y, m_overlayValue);
         m_outBuffer->setValue(x, m_aoiViewRect.lr().y, m_overlayValue);
      }
   }
}

setUsage()

void ossimViewshedTool::setUsage ( ossimArgumentParser &  ap )

virtual

Initializes the aurgument parser with expected parameters and options.

It does not output anything. To see the usage, the caller will need to do something like:

ap.getApplicationUsage()->write(<ostream>);

Reimplemented from ossimChipProcTool.

Definition at line 79 of file ossimViewshedTool.cpp.

References ossimApplicationUsage::addCommandLineOption(), DESCRIPTION, ossimArgumentParser::getApplicationName(), ossimArgumentParser::getApplicationUsage(), ossimApplicationUsage::setCommandLineUsage(), ossimApplicationUsage::setDescription(), and ossimChipProcTool::setUsage().

{
   // Set the general usage:
   ossimApplicationUsage* au = ap.getApplicationUsage();
   ossimString appName = ap.getApplicationName();
   ossimString usageString = appName;
   usageString += " viewshed [options] <obs_lat> <obs_lon> <output-image>";
   au->setCommandLineUsage(usageString);

   // Set the command line options:
   au->addCommandLineOption(
         "--fov <start> <end>", "Optional arguments specifying the field-of"
         "-view boundary azimuths (in degrees). By default, a 360 deg FOV is"
         " computed. The arc is taken clockwise from start to end, so for a"
         " FOV of 225 deg from W, through N to SE, start=270 and end=135");
   au->addCommandLineOption(
         "--height-of-eye <meters>", "Specifies the observers height-of-eye above the "
         "terrain in meters. Defaults to 1.5 m.");
   au->addCommandLineOption(
         "--horizon <filename>", "Experimental. Outputs the max elevation angles "
         "for all azimuths to <filename>, for horizon profiling.");
   au->addCommandLineOption(
         "--radius <meters>", "Specifies max visibility in meters. Required "
         "unless --size is specified. This option constrains output to a circle, "
         "similar to a radar display");
   au->addCommandLineOption(
         "--reticle <int>", "Specifies the size of the reticle at the observer"
         "location in pixels from the center (i.e., the radius of the reticle). "
         "Defaults to 2. A value of 0 hides the reticle. See --values option for "
         "setting reticle color.");
   au->addCommandLineOption(
         "--simulation", "For engineering/debug purposes ");
   au->addCommandLineOption(
         "--tbs", "\"Thread By Sector\". For engineering/debug purposes ");
   au->addCommandLineOption(
         "--threads <n>", "Number of threads. Defaults to use all available cores. "
         "For engineering/debug purposes ");
   au->addCommandLineOption(
         "--viewshed-coding <int int int>", "Specifies the pixel values (0-255) for the visible,"
         " hidden and overlay pixels, respectively. Defaults to visible=null (0), "
         "hidden=128, and overlay (observer position, reticle, and circumference) is "
         "highlighted with 255.");

   // Base class has its own:
   ossimChipProcTool::setUsage(ap);

   ostringstream description;
   description << DESCRIPTION << "\n\nExamples:\n\n"
         "    "<<appName<<" viewshed --radius 50  28.0 -80.5 output-hlz.tif\n"
         "\nAn alternate command line provides switch for observer lat and lon:\n\n"
         "    "<<appName<<" viewshed --rlz 25 --observer 28.0 -80.5  output-hlz.tif \n";
   au->setDescription(description.str());
}

test()

void ossimViewshedTool::test

(

)

For engineering/debug.

Definition at line 1013 of file ossimViewshedTool.cpp.

References ossimChipProcTool::m_aoiGroundRect, m_observerGpt, m_startFov, m_stopFov, and optimizeFOV().

{
   m_aoiGroundRect = ossimGrect(1.0, 0.0, 0.0, 1.0);

   m_observerGpt = ossimGpt(1.5, 0.5);
   m_startFov = 180;
   m_stopFov = 270;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;

   m_observerGpt = ossimGpt(1.5, 0.5);
   m_startFov = 335;
   m_stopFov = 180;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;

   m_observerGpt = ossimGpt(1.5, 0.5);
   m_startFov = 270;
   m_stopFov = 0;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;

   m_observerGpt = ossimGpt(-0.5, 0.5);
   m_startFov = 270;
   m_stopFov = 10;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;

   m_observerGpt = ossimGpt(-0.5, 0.5);
   m_startFov = 350;
   m_stopFov = 90;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;

   m_observerGpt = ossimGpt(-0.5, 0.5);
   m_startFov = 10;
   m_stopFov = 20;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;

   m_observerGpt = ossimGpt(-0.5, 0.5);
   m_startFov = 270;
   m_stopFov = 90;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;

   m_observerGpt = ossimGpt(-0.5, 0.5);
   m_startFov = 90;
   m_stopFov = 270;
   cout<<"Before: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<endl;
   optimizeFOV();
   cout<<"After: m_startFov="<<m_startFov<<"  m_stopFov="<<m_stopFov<<"\n"<<endl;
}

writeHorizonProfile()

bool ossimViewshedTool::writeHorizonProfile ( )

protected

Definition at line 829 of file ossimViewshedTool.cpp.

References ossim::atand(), ossimString::chars(), m_halfWindow, m_horizonFile, m_horizonMap, and m_radials.

Referenced by execute().

{
   // Store the max elevation angles for horizon profiling:
   double az_deg, arctan;
   for (ossim_uint32 sector=0; sector<8; ++sector)
   {
      if (m_radials[sector] == 0)
         continue;

      for (ossim_uint32 radial = 0; radial <= m_halfWindow; ++radial)
      {
         arctan = ossim::atand(m_radials[sector][radial].azimuth);
         switch (sector)
         {
         case 0: // 0 - 45
            az_deg = arctan;
            break;
         case 1: // 45 - 90
            az_deg = 90 - arctan;
            break;
         case 2: // 90 - 135
            az_deg = 90 + arctan;
            break;
         case 3: // 135 - 180
            az_deg = 180 - arctan;
            break;
         case 4: // 180 - 225
            az_deg = 180 + arctan;
            break;
         case 5: // 225 - 270
            az_deg = 270 - arctan;
            break;
         case 6: // 270 - 315
            az_deg = 270 + arctan;
            break;
         case 7: // 315 - 360
            az_deg = 360 - arctan;
            break;
         default:
            break;
         }

         m_horizonMap.insert(pair<double, double>(az_deg, m_radials[sector][radial].elevation));
      }
   }

   // Open output file and write the map:
   ofstream fstr (m_horizonFile.chars());
   if (!fstr.is_open())
      return false;
   map<double, double>::iterator iter = m_horizonMap.begin();
   while (iter != m_horizonMap.end())
   {
      fstr << iter->first << ", " << iter->second << endl;
      ++iter;
   }

   fstr.close();
   return true;
}

Friends And Related Function Documentation

RadialProcessor

friend class RadialProcessor

friend

Definition at line 25 of file ossimViewshedTool.h.

RadialProcessorJob

friend class RadialProcessorJob

friend

Definition at line 24 of file ossimViewshedTool.h.

SectorProcessorJob

friend class SectorProcessorJob

friend

Definition at line 23 of file ossimViewshedTool.h.

Member Data Documentation

d_accumT

double ossimViewshedTool::d_accumT

protected

Definition at line 112 of file ossimViewshedTool.h.

Referenced by computeViewshed().

d_mutex

std::mutex ossimViewshedTool::d_mutex

protected

Definition at line 113 of file ossimViewshedTool.h.

DESCRIPTION

const char * ossimViewshedTool::DESCRIPTION

static

Initial value:

=
      "Computes bitmap image representing the viewshed from specified location using only "
      "DEM information."

Used by ossimUtilityFactory.

Definition at line 61 of file ossimViewshedTool.h.

Referenced by ossimToolFactory::getCapabilities(), and setUsage().

m_displayAsRadar

bool ossimViewshedTool::m_displayAsRadar

protected

Definition at line 94 of file ossimViewshedTool.h.

Referenced by RadialProcessor::doRadial(), and paintReticle().

m_halfWindow

ossim_uint32 ossimViewshedTool::m_halfWindow

protected

Definition at line 95 of file ossimViewshedTool.h.

Referenced by computeViewshed(), RadialProcessor::doRadial(), initProcessingChain(), initRadials(), and writeHorizonProfile().

m_hiddenValue

ossim_uint8 ossimViewshedTool::m_hiddenValue

protected

Definition at line 99 of file ossimViewshedTool.h.

Referenced by RadialProcessor::doRadial().

m_horizonFile

ossimFilename ossimViewshedTool::m_horizonFile

protected

Definition at line 108 of file ossimViewshedTool.h.

Referenced by execute(), and writeHorizonProfile().

m_horizonMap

std::map<double, double> ossimViewshedTool::m_horizonMap

protected

Definition at line 109 of file ossimViewshedTool.h.

Referenced by clear(), and writeHorizonProfile().

m_jobMtQueue

std::shared_ptr<ossimJobMultiThreadQueue> ossimViewshedTool::m_jobMtQueue

protected

Definition at line 103 of file ossimViewshedTool.h.

Referenced by clear(), and computeViewshed().

m_memSource

ossimRefPtr<ossimMemoryImageSource> ossimViewshedTool::m_memSource

protected

Definition at line 97 of file ossimViewshedTool.h.

Referenced by computeViewshed(), and initProcessingChain().

m_numThreads

ossim_uint32 ossimViewshedTool::m_numThreads

protected

Definition at line 104 of file ossimViewshedTool.h.

Referenced by computeViewshed().

m_observerGpt

ossimGpt ossimViewshedTool::m_observerGpt

protected

Definition at line 88 of file ossimViewshedTool.h.

Referenced by clear(), computeRadius(), RadialProcessor::doRadial(), initializeAOI(), initializeProjectionGsd(), initProcessingChain(), optimizeFOV(), ossimViewshedTool(), paintReticle(), and test().

m_observerVpt

ossimDpt ossimViewshedTool::m_observerVpt

protected

Definition at line 89 of file ossimViewshedTool.h.

Referenced by RadialProcessor::doRadial(), initProcessingChain(), and paintReticle().

m_obsHgtAbvTer

double ossimViewshedTool::m_obsHgtAbvTer

protected

Definition at line 90 of file ossimViewshedTool.h.

Referenced by RadialProcessor::doRadial(), and initProcessingChain().

m_obsInsideAoi

bool ossimViewshedTool::m_obsInsideAoi

protected

Definition at line 93 of file ossimViewshedTool.h.

m_outBuffer

ossimRefPtr<ossimImageData> ossimViewshedTool::m_outBuffer

protected

Definition at line 96 of file ossimViewshedTool.h.

Referenced by clear(), computeViewshed(), RadialProcessor::doRadial(), initProcessingChain(), and paintReticle().

m_overlayValue

ossim_uint8 ossimViewshedTool::m_overlayValue

protected

Definition at line 100 of file ossimViewshedTool.h.

Referenced by RadialProcessor::doRadial(), and paintReticle().

m_radials

Radial** ossimViewshedTool::m_radials

protected

Definition at line 92 of file ossimViewshedTool.h.

Referenced by computeViewshed(), RadialProcessor::doRadial(), initRadials(), writeHorizonProfile(), and ~ossimViewshedTool().

m_reticleSize

ossim_int32 ossimViewshedTool::m_reticleSize

protected

Definition at line 101 of file ossimViewshedTool.h.

Referenced by paintReticle().

m_simulation

bool ossimViewshedTool::m_simulation

protected

Definition at line 102 of file ossimViewshedTool.h.

Referenced by RadialProcessor::doRadial().

m_startFov

double ossimViewshedTool::m_startFov

protected

Definition at line 105 of file ossimViewshedTool.h.

Referenced by initRadials(), optimizeFOV(), and test().

m_stopFov

double ossimViewshedTool::m_stopFov

protected

Definition at line 106 of file ossimViewshedTool.h.

Referenced by initRadials(), optimizeFOV(), and test().

m_threadBySector

bool ossimViewshedTool::m_threadBySector

protected

Definition at line 107 of file ossimViewshedTool.h.

Referenced by computeViewshed().

m_visibleValue

ossim_uint8 ossimViewshedTool::m_visibleValue

protected

Definition at line 98 of file ossimViewshedTool.h.

Referenced by RadialProcessor::doRadial().

m_visRadius

double ossimViewshedTool::m_visRadius

protected

Definition at line 91 of file ossimViewshedTool.h.

Referenced by clear(), computeRadius(), initializeAOI(), and initProcessingChain().

---

The documentation for this class was generated from the following files:

• ossimViewshedTool.h
• ossimViewshedTool.cpp