RadialProcessor Class Reference

This class provides a common entry point for both SectorProcessorJob and RadialProcessorJob for processing a single radial. More...

#include <ossimViewshedTool.h>

Static Public Member Functions
static void	doRadial (ossimViewshedTool *vs, ossim_uint32 s, ossim_uint32 r)

Private Member Functions
	RadialProcessor ()

Static Private Attributes
static std::mutex	m_bufMutex

Detailed Description

This class provides a common entry point for both SectorProcessorJob and RadialProcessorJob for processing a single radial.

Eventually, SectorProcessorJob can likely go away (invoked with the "--tbs" command-line option, and doRadial() method can be moved into RadialProcessorJob class. In the meantime, both Sector/thread and Radial/thread schemes are supported to continue evaluating performance.

Definition at line 167 of file ossimViewshedTool.h.

Constructor & Destructor Documentation

RadialProcessor()

RadialProcessor::RadialProcessor ( )

inlineprivate

Definition at line 174 of file ossimViewshedTool.h.

{};

Member Function Documentation

doRadial()

void RadialProcessor::doRadial ( ossimViewshedTool *  vs, ossim_uint32  s, ossim_uint32  r  )

static

Definition at line 906 of file ossimViewshedTool.cpp.

References ossimViewshedTool::Radial::azimuth, ossimViewshedTool::Radial::elevation, ossimGpt::hasNans(), ossimGpt::hgt, ossimViewshedTool::Radial::insideAoi, ossim::isnan(), ossimImageGeometry::localToWorld(), ossimChipProcTool::m_aoiViewRect, ossimViewshedTool::m_displayAsRadar, ossimChipProcTool::m_geom, ossimViewshedTool::m_halfWindow, ossimViewshedTool::m_hiddenValue, ossimViewshedTool::m_observerGpt, ossimViewshedTool::m_observerVpt, ossimViewshedTool::m_obsHgtAbvTer, ossimViewshedTool::m_outBuffer, ossimViewshedTool::m_overlayValue, ossimViewshedTool::m_radials, ossimViewshedTool::m_simulation, ossimViewshedTool::m_visibleValue, ossimIrect::pointWithin(), ossimImageData::setValue(), ossimIpt::x, ossimDpt::x, ossimIpt::y, and ossimDpt::y.

Referenced by SectorProcessorJob::run(), and RadialProcessorJob::run().

{
   double u, v;
   ossimDpt pt_i, vpt_i;
   ossimGpt gpt_i;
   double elev_i, elev;
   double r2_max = vsUtil->m_halfWindow*vsUtil->m_halfWindow;

   // Establish shorthand access to radial:
   ossimViewshedTool::Radial& radial = vsUtil->m_radials[sector_idx][radial_idx];

   // Walk along the radial using the appropriate coordinate abscissa for that sector and
   // compute ordinate using the radials azimuth:
   for (u=1.0; u <= (double) vsUtil->m_halfWindow; u += 1.0)
   {
      // Compute ordinate from abscissa and slope of this radial:
      v = radial.azimuth*(u);
      switch (sector_idx)
      {
      case 0: // N-NE, (u, v) = (-y, x)
         pt_i.y = -u;
         pt_i.x = v;
         break;
      case 1: // NE-E, (u, v) = (x, -y)
         pt_i.x = u;
         pt_i.y = -v;
         break;
      case 2: // E-SE, (u, v) = (x, y)
         pt_i.x = u;
         pt_i.y = v;
         break;
      case 3: // SE-S, (u, v) = (y, x)
         pt_i.y = u;
         pt_i.x = v;
         break;
      case 4: // S-SW, (u, v) = (y, -x)
         pt_i.y = u;
         pt_i.x = -v;
         break;
      case 5: // SW-W, (u, v) = (-x, y)
         pt_i.x = -u;
         pt_i.y = v;
         break;
      case 6: // W-NW, (u, v) = (-x, -y)
         pt_i.x = -u;
         pt_i.y = -v;
         break;
      case 7: // NW-N, (u, v) = (-y, -x)
         pt_i.y = -u;
         pt_i.x = -v;
         break;
      default:
         break;
      }

      // Shift to actual view coordinates:
      vpt_i = pt_i + vsUtil->m_observerVpt;
      ossimIpt ipt (vpt_i);

      // Check if alread accounted for at this location:
      //if (!vsUtil->m_outBuffer->isNull(vpt_i))
      //   continue;

      // Check if we are exiting the AOI (no more processing required for this radial):
      bool pointInsideAoi = vsUtil->m_aoiViewRect.pointWithin(ipt);
      if (radial.insideAoi && !pointInsideAoi)
         break;

      // Alternatively, check if we were OUTSIDE and now moving INSIDE:
      if (!radial.insideAoi && pointInsideAoi)
         radial.insideAoi = true;

      // Check if we passed beyong the visibilty radius, and exit loop if so:
      if (vsUtil->m_displayAsRadar && ((u*u + v*v) >= r2_max))
      {
         vsUtil->m_outBuffer->setValue(ipt.x, ipt.y, vsUtil->m_overlayValue);
         break;
      }

      // Fetch the pixel value as the elevation value and compute elevation angle from
      // the observer pt as dz/dx
      vsUtil->m_geom->localToWorld(vpt_i, gpt_i);
      if (vsUtil->m_simulation && ossim::isnan(gpt_i.hgt))
         gpt_i.hgt = vsUtil->m_observerGpt.hgt-vsUtil->m_obsHgtAbvTer; // ground level

      if (!gpt_i.hasNans())
      {
         // Compare elev angle to max angle latched so far along this radial:
         elev_i = (gpt_i.hgt - vsUtil->m_observerGpt.hgt) / u;
         elev = radial.elevation;
         if (elev_i > elev)
         {
            // point is visible, latch this line-of-sight as the new max elevation angle for this
            // radial, and mark the output pixel as visible:
            radial.elevation = elev_i;
            vsUtil->m_outBuffer->setValue(ipt.x, ipt.y, vsUtil->m_visibleValue);
         }
         else
         {
            vsUtil->m_outBuffer->setValue(ipt.x, ipt.y, vsUtil->m_hiddenValue);
         }
      }
   } // end loop over radial's abscissas
}

Member Data Documentation

m_bufMutex

std::mutex RadialProcessor::m_bufMutex

staticprivate

Definition at line 173 of file ossimViewshedTool.h.

---

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

• ossimViewshedTool.h
• ossimViewshedTool.cpp