ATP::AtpGenerator Class Reference

Base class for OSSIM-based ATP generators. More...

#include <AtpGenerator.h>

Inheritance diagram for ATP::AtpGenerator:

Public Types
enum	Algorithm { ALGO_UNASSIGNED =0, CROSSCORR, DESCRIPTOR, NASA }

Public Member Functions
	AtpGenerator (Algorithm algo)
virtual	~AtpGenerator ()
void	setRefImage (std::shared_ptr< ossim::Image > ref_image)
void	setCmpImage (std::shared_ptr< ossim::Image > cmp_image)
virtual void	initialize ()
	Needs to be called after ref and cmp images are set. More...
void	setBandWeights (const std::vector< double > weights)
	When the input images are multiband, the bands must be combined into a single-band image. More...
virtual bool	generateTiePointList (ossim::TiePointList &tpList)
	This is the main workhorse method. More...
ossimRefPtr< ossimImageViewProjectionTransform >	getRefIVT ()
ossimRefPtr< ossimImageViewProjectionTransform >	getCmpIVT ()
ossimRefPtr< ossimImageChain >	getRefChain ()
ossimRefPtr< ossimImageChain >	getCmpChain ()
std::string	getRefImageID ()
std::string	getCmpImageID ()
std::string	getRefFilename ()
std::string	getCmpFilename ()
ossimRefPtr< AtpTileSource >	getAtpTileSource ()

Static Public Member Functions
static void	writeTiePointList (ostream &out, const AtpList &tpList)
	For engineering use. More...

Public Attributes
ossimRefPtr< AtpAnnotatedImage >	m_annotatedRefImage
ossimRefPtr< AtpAnnotatedImage >	m_annotatedCmpImage

Protected Member Functions
	AtpGenerator ()
virtual ossimRefPtr< ossimImageChain >	constructChain (std::shared_ptr< ossim::Image > image, ossimRefPtr< ossimImageViewProjectionTransform > &ivt)
	Constructs the processing chain for the input image according to the needs of the generator. More...
bool	getValidVertices (ossimRefPtr< ossimImageChain > chain, ossimRefPtr< ossimImageViewProjectionTransform > &ivt, std::vector< ossimDpt > &validVertices)
	Establishes valid image vertices in view space for later computing overlap. More...
void	layoutSearchTileRects (ossimPolygon &overlapPoly)
	Finds optimum layout of patches within the intersect area for feature search. More...
ossimImageHandler *	getImageHandler (ossimRefPtr< ossimImageChain > &chain)

Protected Attributes
Algorithm	m_algorithm
std::shared_ptr< ossim::Image >	m_refImage
std::shared_ptr< ossim::Image >	m_cmpImage
ossimRefPtr< AtpTileSource >	m_atpTileSource
ossimRefPtr< ossimImageChain >	m_refChain
ossimRefPtr< ossimImageChain >	m_cmpChain
ossimRefPtr< ossimImageViewProjectionTransform >	m_refIVT
ossimRefPtr< ossimImageViewProjectionTransform >	m_cmpIVT
ossimRefPtr< ossimImageGeometry >	m_viewGeom
ossimIrect	m_aoiView
std::vector< double >	m_bandWeights
double	m_refEllipHgt
std::vector< ossimIrect >	m_searchTileRects

Static Protected Attributes
static std::shared_ptr< AutoTiePoint >	s_referenceATP

Friends
class	AtpTileSource

Detailed Description

Base class for OSSIM-based ATP generators.

TODO: Currently, ATP generation only works with pair-wise matching. Eventually it should support N-way matching with an arbitrary number of input images. The multiple image case is presently handled by the AutoTiePointService.

Definition at line 33 of file AtpGenerator.h.

Member Enumeration Documentation

Algorithm

enum ATP::AtpGenerator::Algorithm

Enumerator
ALGO_UNASSIGNED
CROSSCORR
DESCRIPTOR
NASA

Definition at line 38 of file AtpGenerator.h.

{ ALGO_UNASSIGNED=0, CROSSCORR, DESCRIPTOR, NASA };

Constructor & Destructor Documentation

AtpGenerator() [1/2]

ATP::AtpGenerator::AtpGenerator

(

Algorithm

algo

)

Definition at line 39 of file AtpGenerator.cpp.

:   m_algorithm (algo),
    m_refEllipHgt(0)
{
}

~AtpGenerator()

ATP::AtpGenerator::~AtpGenerator ( )

virtual

Definition at line 46 of file AtpGenerator.cpp.

{
}

AtpGenerator() [2/2]

ATP::AtpGenerator::AtpGenerator ( )

inlineprotected

Definition at line 95 of file AtpGenerator.h.

: m_algorithm (ALGO_UNASSIGNED), m_refEllipHgt(0) {}

Member Function Documentation

constructChain()

ossimRefPtr< ossimImageChain > ATP::AtpGenerator::constructChain ( std::shared_ptr< ossim::Image >  image, ossimRefPtr< ossimImageViewProjectionTransform > &  ivt  )

protectedvirtual

Constructs the processing chain for the input image according to the needs of the generator.

As a convenience, the image-to-common view IVT is initialized

Definition at line 211 of file AtpGenerator.cpp.

References ossimImageChain::add(), ossim::JsonParam::asBool(), ossim::JsonParam::asFloat(), ossim::JsonParam::asString(), ossimImageHandler::buildHistogram(), CINFO, CROSSCORR, ossim::JsonConfig::diagnosticLevel(), ossimString::empty(), ossimRefPtr< T >::get(), ossimImageGeometry::getAsMapProjection(), ossimImageGeometry::getBoundingRect(), ossimImageHandler::getImageGeometry(), ossimImageHandler::getImageHistogram(), ossimImageHandler::getImageID(), ossimImageGeometry::getMetersPerPixel(), ossimImageSource::getNumberOfOutputBands(), ossimImageChain::getOutputScalarType(), ossim::JsonConfig::getParameter(), ossimImageGeometry::getTiePoint(), ossimImageViewProjectionTransform::imageToView(), ATP::AtpConfig::instance(), ossimImageHandlerRegistry::instance(), ossim::isnan(), ossimHistogramRemapper::LINEAR_AUTO_MIN_MAX, m_algorithm, m_viewGeom, ossimImageHandlerRegistry::open(), OSSIM_UINT8, ossimImageHandler::setCurrentEntry(), ossimHistogramRemapper::setHistogram(), ossimImageHandler::setImageID(), ossimImageRenderer::setImageViewTransform(), ossimMapProjection::setMetersPerPixel(), ossimUtmProjection::setOrigin(), ossimBandSelector::setOutputBandList(), ossimScalarRemapper::setOutputScalarType(), ossimImageRenderer::setProperty(), ossimHistogramRemapper::setStretchMode(), ossimMapProjection::setUlTiePoints(), ossimImageViewProjectionTransform::viewToImage(), ossimDpt::x, and ossimDpt::y.

Referenced by initialize().

{
   static const char *MODULE = "AtpGenerator::constructChain()  ";
   AtpConfig &config = AtpConfig::instance();

   ossimImageHandlerRegistry *factory = ossimImageHandlerRegistry::instance();
   ossimRefPtr<ossimImageHandler> handler = factory->open(image->getFilename());
   if (!handler)
      return NULL;
   handler->setCurrentEntry(image->getEntryIndex());

   // Check if an imageID was provided in job JSON. otherwise, check if image handler can get it:
   ossimString imageId = image->getImageId();
   if (imageId.empty())
   {
      imageId = handler->getImageID();
      if (imageId.empty())
         imageId = image->getFilename();
   }
   handler->setImageID(imageId);
   image->setImageId(imageId);

   if (config.diagnosticLevel(5))
   {
      ossimIrect vuRect;
      handler->getImageGeometry()->getBoundingRect(vuRect);
      CINFO << MODULE << "m_aoiView: " << vuRect << endl;;
   }

   ossimRefPtr<ossimImageChain> chain = new ossimImageChain;
   chain->add(handler.get());

   // Add histogram (experimental)
   bool doHistogramStretch = config.getParameter("doHistogramStretch").asBool();
   if (doHistogramStretch)
   {
      ossimRefPtr<ossimHistogramRemapper> histo_remapper = new ossimHistogramRemapper();
      chain->add(histo_remapper.get());
      histo_remapper->setStretchMode(ossimHistogramRemapper::LINEAR_AUTO_MIN_MAX);
      ossimRefPtr<ossimMultiResLevelHistogram> histogram = handler->getImageHistogram();
      if (!histogram)
      {
         handler->buildHistogram();
         histogram = handler->getImageHistogram();
      }
      histo_remapper->setHistogram(histogram);

      // Need to select a band if multispectral:
      unsigned int numBands = handler->getNumberOfOutputBands();
      if (numBands > 1)
      {
         unsigned int band = image->getActiveBand();
         if ((band > numBands) || (band <= 0))
         {
            CINFO << MODULE << "Specified band (" << band << ") is outside allowed range (1 to "
                  << numBands
                  << "). Using default band 1 which may not be ideal." << endl;
            band = 1;
         }
         band--; // shift to zero-based
         ossimRefPtr<ossimBandSelector> band_selector = new ossimBandSelector();
         chain->add(band_selector.get());
         vector<ossim_uint32> bandList;
         bandList.push_back(band);
         band_selector->setOutputBandList(bandList);
      }
   }

   // Can only work in 8-bit radiometry:
   if (chain->getOutputScalarType() != OSSIM_UINT8)
   {
      ossimRefPtr<ossimScalarRemapper> remapper = new ossimScalarRemapper();
      chain->add(remapper.get());
      remapper->setOutputScalarType(OSSIM_UINT8);
   }

   // Finally, add a cache to the input chain:
   ossimRefPtr<ossimCacheTileSource> cache = new ossimCacheTileSource();
   chain->add(cache.get());

   // Check for specified GSD in config parameters (experimental):
   double viewGsd = config.getParameter("viewGSD").asFloat();

   // Set up a common view geometry, even if derived class works exclusively in image space.
   ossimRefPtr<ossimImageGeometry> image_geom = handler->getImageGeometry();
   if (!m_viewGeom)
   {
      ossimDpt gsd;
      if (!ossim::isnan(viewGsd) && (viewGsd != 0.0))
         gsd = ossimDpt(viewGsd, viewGsd);
      else
      {
         gsd = image_geom->getMetersPerPixel();
         if (gsd.x > gsd.y)
            gsd.y = gsd.x;
         else
            gsd.x = gsd.y;
      }
      ossimGpt proj_origin;
      image_geom->getTiePoint(proj_origin, false);
      ossimUtmProjection *proj = new ossimUtmProjection();
      proj->setOrigin(proj_origin);
      proj->setMetersPerPixel(gsd);
      proj->setUlTiePoints(proj_origin);
      m_viewGeom = new ossimImageGeometry(NULL, proj);
   }
   else
   {
      // View geometry was already set up by first image of the pair (REF). Only need to check the
      // GSD of the CMP image and adjust the view GSD to match if larger:
      ossimDpt cmpGsd(image_geom->getMetersPerPixel());
      if (cmpGsd.x > cmpGsd.y)
         cmpGsd.y = cmpGsd.x;
      else
         cmpGsd.x = cmpGsd.y;
      ossimDpt refGsd(m_viewGeom->getMetersPerPixel());
      if (cmpGsd.x > refGsd.x)
      {
         refGsd.x = cmpGsd.x;
         refGsd.y = cmpGsd.y;
         m_viewGeom->getAsMapProjection()->setMetersPerPixel(refGsd);
      }
   }

   ivt = new ossimImageViewProjectionTransform(image_geom.get(), m_viewGeom.get());

   if (m_algorithm == CROSSCORR)
   {
      // Need to add a renderer since correlations are done in view-space:
      ossimImageRenderer* renderer = new ossimImageRenderer;
      chain->add(renderer);
      renderer->setImageViewTransform(ivt.get());

      // Set the appropriate resampling filter for ATP. TODO: Need experimentation.
      ossimString filterType = AtpConfig::instance().getParameter("resamplingMethod").asString();
      if (filterType.empty())
         filterType = "bilinear";
      ossimRefPtr<ossimProperty> p = new ossimStringProperty("filter_type", filterType);
      renderer->setProperty(p);

      // Add cache after resampler:
      ossimRefPtr<ossimCacheTileSource> cache = new ossimCacheTileSource();
      chain->add(cache.get());

      if (AtpConfig::instance().diagnosticLevel(3))
      {
         ossimDpt imagePt(1,1);
         ossimDpt viewPt, rtIpt;
         ivt->imageToView(imagePt, viewPt);
         ivt->viewToImage(viewPt, rtIpt);
         CINFO<<MODULE<<"\n  imagePt: "<<imagePt<<
              "\n  viewPt: "<<viewPt<<"\n  rtIpt: "<<rtIpt<<endl;
      }

   }
   return chain;
}

generateTiePointList()

bool ATP::AtpGenerator::generateTiePointList ( ossim::TiePointList &  tpList )

virtual

This is the main workhorse method.

It launches the auto-tie-point generation process given the REF and CMP images. TODO: Currently, ATP generation only works with pair-wise matching. Eventually it should support N-way matching with an arbitrary number of input images.

Definition at line 450 of file AtpGenerator.cpp.

References CINFO, ossim::JsonConfig::diagnosticLevel(), initialize(), ATP::AtpConfig::instance(), m_annotatedCmpImage, m_annotatedRefImage, m_atpTileSource, m_cmpChain, m_refChain, and m_searchTileRects.

{
   const char* MODULE = "AtpGenerator::generateTiePointList()  ";
   initialize();

   if (!m_refChain || !m_cmpChain || !m_atpTileSource)
      return false;

   AtpConfig& config = AtpConfig::instance();

   // Need to search for features in the REF tile first, then consider only those locations
   // for correlating. The search tiles were previously established in the base class initialization
   ossimRefPtr<ossimImageData> ref_tile = 0;
   ossim_uint32 tileId = 0;
   for (auto &searchTileRect : m_searchTileRects)
   {
      ref_tile = m_atpTileSource->getTile(searchTileRect);

      AtpList &tileATPs = m_atpTileSource->getTiePoints();
      if (tileATPs.empty() && config.diagnosticLevel(1))
      {
         CINFO << "\n" << MODULE << "No TPs found in tile." << endl;
         ++tileId;
         continue;
      }

      // Add remaining tile ATPs to master ATP list:
      if (!tileATPs.empty())
         atpList.insert(atpList.end(), tileATPs.begin(), tileATPs.end());

      ++tileId;
   }

   if (config.diagnosticLevel(1))
      CINFO<<"\n"<<MODULE<<"Total number of TPs found in pair = "<<atpList.size()<<endl;

   if (config.diagnosticLevel(3))
   {
      m_annotatedRefImage->write();
      if (config.diagnosticLevel(4))
         m_annotatedCmpImage->write();
   }
   return true;
}

getAtpTileSource()

ossimRefPtr<AtpTileSource> ATP::AtpGenerator::getAtpTileSource ( )

inline

Definition at line 90 of file AtpGenerator.h.

References m_atpTileSource.

{ return m_atpTileSource; } // For engineering use

getCmpChain()

ossimRefPtr<ossimImageChain> ATP::AtpGenerator::getCmpChain ( )

inline

Definition at line 83 of file AtpGenerator.h.

References m_cmpChain.

{ return m_cmpChain; }

getCmpFilename()

string ATP::AtpGenerator::getCmpFilename

(

)

Definition at line 87 of file AtpGenerator.cpp.

References ossimImageHandler::getFilename(), getImageHandler(), m_cmpChain, and ossimString::string().

{
   ossimImageHandler* handler = getImageHandler(m_cmpChain);
   if (!handler)
      return "";
   return handler->getFilename().string();
}

getCmpImageID()

string ATP::AtpGenerator::getCmpImageID

(

)

Definition at line 71 of file AtpGenerator.cpp.

References getImageHandler(), ossimImageHandler::getImageID(), m_cmpChain, and ossimString::string().

{
   ossimImageHandler* handler = getImageHandler(m_cmpChain);
   if (!handler)
      return "";
   return handler->getImageID().string();
}

getCmpIVT()

ossimRefPtr<ossimImageViewProjectionTransform> ATP::AtpGenerator::getCmpIVT ( )

inline

Definition at line 81 of file AtpGenerator.h.

References m_cmpIVT.

{ return m_cmpIVT; }

getImageHandler()

ossimImageHandler * ATP::AtpGenerator::getImageHandler ( ossimRefPtr< ossimImageChain > &  chain )

protected

Definition at line 95 of file AtpGenerator.cpp.

References ossimImageChain::accept(), and ossimCollectionVisitor::getObjectAs().

Referenced by getCmpFilename(), getCmpImageID(), getRefFilename(), and getRefImageID().

{
   if (!chain)
      return 0;

   ossimTypeNameVisitor visitor ("ossimImageHandler");
   chain->accept(visitor);
   ossimImageHandler* handler = (ossimImageHandler*) visitor.getObjectAs<ossimImageHandler>();

   return handler;
}

getRefChain()

ossimRefPtr<ossimImageChain> ATP::AtpGenerator::getRefChain ( )

inline

Definition at line 82 of file AtpGenerator.h.

References m_refChain.

{ return m_refChain; }

getRefFilename()

string ATP::AtpGenerator::getRefFilename

(

)

Definition at line 79 of file AtpGenerator.cpp.

References ossimImageHandler::getFilename(), getImageHandler(), m_refChain, and ossimString::string().

{
   ossimImageHandler* handler = getImageHandler(m_refChain);
   if (!handler)
      return "";
   return handler->getFilename().string();
}

getRefImageID()

string ATP::AtpGenerator::getRefImageID

(

)

Definition at line 63 of file AtpGenerator.cpp.

References getImageHandler(), ossimImageHandler::getImageID(), m_refChain, and ossimString::string().

{
   ossimImageHandler* handler = getImageHandler(m_refChain);
   if (!handler)
      return "";
   return handler->getImageID().string();
}

getRefIVT()

ossimRefPtr<ossimImageViewProjectionTransform> ATP::AtpGenerator::getRefIVT ( )

inline

Definition at line 80 of file AtpGenerator.h.

References m_refIVT.

{ return m_refIVT; }

getValidVertices()

bool ATP::AtpGenerator::getValidVertices ( ossimRefPtr< ossimImageChain >  chain, ossimRefPtr< ossimImageViewProjectionTransform > &  ivt, std::vector< ossimDpt > &  validVertices  )

protected

Establishes valid image vertices in view space for later computing overlap.

Definition at line 370 of file AtpGenerator.cpp.

References ossimImageHandler::createDefaultValidVerticesFilename(), ossimImageHandler::getBoundingRect(), ossimImageHandler::getImageGeometry(), ossimImageChain::getLastSource(), ossimImageGeometry::getProjection(), ossimImageHandler::getValidImageVertices(), ossimImageViewProjectionTransform::imageToView(), and ossimImageHandler::openValidVertices().

Referenced by initialize().

{
   validVertices.clear();

   // Need to establish the valid image vertices for establishing overlap polygon. Compute in
   // image space then transform the vertices to view space:
   if (!chain)
      return false;
   ossimImageHandler *handler = dynamic_cast<ossimImageHandler *>(chain->getLastSource());
   if (!handler)
      return false;

   // Check if there is already a vertices file alongside the image:
   if (handler->openValidVertices())
   {
      vector<ossimIpt> iverts;
      handler->getValidImageVertices(iverts);
      for (const auto &vert : iverts)
         validVertices.emplace_back(vert);
   }

   if (validVertices.empty())
   {
      // If the projection is a sensor model, then assume that the four image-space coordinates are
      // valid image coordinates.
      ossimRefPtr<ossimImageGeometry> geom = handler->getImageGeometry();
      if (geom)
      {
         if (dynamic_cast<ossimSensorModel*>(geom->getProjection()))
         {
            ossimIrect rect (handler->getBoundingRect());
            validVertices.emplace_back(rect.ul());
            validVertices.emplace_back(rect.ur());
            validVertices.emplace_back(rect.lr());
            validVertices.emplace_back(rect.ll());
         }
      }

      if (validVertices.empty())
      {
         ossimFilename verticesFile(handler->createDefaultValidVerticesFilename());
         ossimRefPtr<ossimVertexExtractor> ve = new ossimVertexExtractor(handler);
         ve->setOutputName(verticesFile);
         ve->setAreaOfInterest(handler->getBoundingRect(0));
         ve->execute();
         const vector<ossimIpt> &iverts = ve->getVertices();
         for (const auto &vert : iverts)
            validVertices.emplace_back(vert);
      }
   }

   // The image vertices are necessarily in image space coordinates. Need to transform to a common
   // "view" space before determining the overlap polygon:
   ossimDpt vpt;
   for (auto &vertex : validVertices)
   {
      ivt->imageToView(vertex, vpt);
      vertex = vpt;
   }

   return true;
}

initialize()

void ATP::AtpGenerator::initialize ( )

virtual

Needs to be called after ref and cmp images are set.

Definition at line 108 of file AtpGenerator.cpp.

References ALGO_UNASSIGNED, ossim::JsonParam::asString(), ossim::JsonParam::asUint(), CFATAL, CINFO, constructChain(), CROSSCORR, DESCRIPTOR, ossim::JsonConfig::diagnosticLevel(), ossimPolyArea2d::getBoundingRect(), ossimIrect::getCenter(), ossimElevManager::getHeightAboveEllipsoid(), ossimImageChain::getImageGeometry(), ossim::JsonConfig::getParameter(), getValidVertices(), ossimPolyArea2d::getVisiblePolygons(), ossimGpt::height(), ATP::AtpConfig::instance(), ossimElevManager::instance(), layoutSearchTileRects(), ossimImageGeometry::localToWorld(), ossimIrect::lr(), m_algorithm, m_annotatedCmpImage, m_annotatedRefImage, m_aoiView, m_atpTileSource, m_cmpChain, m_cmpImage, m_cmpIVT, m_refChain, m_refEllipHgt, m_refImage, m_refIVT, m_viewGeom, NASA, s_referenceATP, ossimIrect::set_lr(), ossimIrect::set_ul(), ossimElevManager::setDefaultHeightAboveEllipsoid(), ossimImageGeometry::setImageSize(), ossimElevManager::setUseGeoidIfNullFlag(), ossimIrect::size(), ossimIrect::ul(), ossimIpt::x, and ossimIpt::y.

Referenced by generateTiePointList().

{
   static const char* MODULE=" AtpGenerator::initialize()  ";
   AtpConfig &config = AtpConfig::instance();

   // Initialize the reference ATP containing the static data members used by all ATPs:
   s_referenceATP.reset(new AutoTiePoint());
   m_viewGeom = 0;

   // Establish image processing chains:
   vector<ossimDpt> refVertices, cmpVertices;
   m_refChain = constructChain(m_refImage, m_refIVT);
   m_cmpChain = constructChain(m_cmpImage, m_cmpIVT);
   if (!m_refChain || !m_cmpChain)
   {
      CFATAL<<MODULE<<"Null input chain(s)."<<endl;
      return;
   }

   if (config.diagnosticLevel(1))
   {
      CINFO<<"\n"<<MODULE<<"Initializing image pair with:"
            <<"\n  REF_ID: "<<m_refImage->getImageId()
            <<"\n  CMP_ID: "<<m_cmpImage->getImageId()<<endl;
   }

   // Establish the actual overlap polygon and set AOI as its bounding rect:
   getValidVertices(m_refChain, m_refIVT, refVertices);
   getValidVertices(m_cmpChain, m_cmpIVT, cmpVertices);
   ossimPolyArea2d refPoly (refVertices);
   ossimPolyArea2d cmpPoly (cmpVertices);
   ossimPolyArea2d overlapPoly (refPoly & cmpPoly);
   vector<ossimPolygon> intersectPolys;
   overlapPoly.getVisiblePolygons(intersectPolys);
   if (intersectPolys.empty())
   {
      CINFO<<MODULE<<"No intersect polygon found."<<endl;
      return;
   }
   ossimPolygon intersectPoly (intersectPolys[0]);

   ossimDrect viewDrect;
   overlapPoly.getBoundingRect(viewDrect);
   m_aoiView = viewDrect;
   m_viewGeom->setImageSize(m_aoiView.size());

   if (config.diagnosticLevel(3))
   {
      ossimString basename (config.getParameter("annotatedImageBaseName").asString());
      m_annotatedRefImage = new AtpAnnotatedImage(m_refChain, m_refIVT, m_aoiView);
      m_annotatedRefImage->setImageName(basename+"Ref.tif");
      m_annotatedCmpImage = new AtpAnnotatedImage(m_cmpChain, m_cmpIVT, m_aoiView);
      m_annotatedCmpImage->setImageName(basename+"Cmp.tif");
   }

   // Establish the feature search tile locations:
   layoutSearchTileRects(intersectPoly);

   // Establish the nominal height above ellipsoid for MSL at this location:
   ossimRefPtr<ossimImageGeometry> refGeom = m_refChain->getImageGeometry();
   if (!refGeom)
   {
      CFATAL<<MODULE<<"Null ref image geometry."<<endl;
      return;
   }
   ossimDpt viewCenter;
   ossimGpt geoCenter;
   m_aoiView.getCenter(viewCenter);
   refGeom->localToWorld(viewCenter, geoCenter);
   m_refEllipHgt = ossimElevManager::instance()->getHeightAboveEllipsoid(geoCenter);
   geoCenter.height(m_refEllipHgt); // Sea level
   ossimElevManager::instance()->setDefaultHeightAboveEllipsoid(m_refEllipHgt);
   ossimElevManager::instance()->setUseGeoidIfNullFlag(true);

   // Now need the algorithm-specific tile source that does the pixel pulling:
   switch (m_algorithm)
   {
   case CROSSCORR:
      m_atpTileSource = new ossimCorrelationSource(this);
      break;
   case DESCRIPTOR:
      m_atpTileSource = new ossimDescriptorSource(this);
      break;
   case NASA:
   case ALGO_UNASSIGNED:
   default:
      CFATAL<<MODULE<<"Unhandled algorithm type."<<endl;
      return;
   }

   // Adjust AOI for half-width of correlation window -- this really should only apply to crosscorr,
   // but descriptors also have a sampling window (just not readily known):
   int patch_center = (config.getParameter("corrWindowSize").asUint() + 1) / 2;
   ossimIpt first_pos(m_aoiView.ul().x + patch_center, m_aoiView.ul().y + patch_center);
   ossimIpt last_pos(m_aoiView.lr().x - patch_center, m_aoiView.lr().y - patch_center);
   m_aoiView.set_ul(first_pos);
   m_aoiView.set_lr(last_pos);

   m_atpTileSource->initialize();
}

layoutSearchTileRects()

void ATP::AtpGenerator::layoutSearchTileRects ( ossimPolygon &  overlapPoly )

protected

Finds optimum layout of patches within the intersect area for feature search.

Definition at line 495 of file AtpGenerator.cpp.

References ossimPolygon::area(), ossim::JsonParam::asUint(), CINFO, CWARN, ossim::JsonConfig::diagnosticLevel(), ossimPolygon::getBoundingRect(), ossimPolygon::getNumberOfVertices(), ossim::JsonConfig::getParameter(), ossimIrect::height(), ATP::AtpConfig::instance(), ossimPolygon::isRectWithin(), ossimIrect::lr(), m_annotatedRefImage, m_aoiView, m_searchTileRects, ossimPolygon::removeVertex(), ossimIrect::ul(), ossimIrect::width(), ossimIpt::x, x, ossimIpt::y, and y.

Referenced by initialize().

{
   const char* MODULE = "AtpGenerator::layoutSearchTiles()  ";
   AtpConfig& config = AtpConfig::instance();

   m_searchTileRects.clear();

   // Note this seems to be specific to correlation-based matching, but should be generalized:
   unsigned int featureSideSize = config.getParameter("corrWindowSize").asUint();

   // ### Critical value set here such that we can fit at least 1/2 the number of desired TPs:
   unsigned int numTpsPerTile = config.getParameter("numFeaturesPerTile").asUint();
   double min_area = numTpsPerTile*featureSideSize*featureSideSize/2.0;

   // Found issue with last vertex = first vertex. Remove last if so:
   int lastVertex = intersectPoly.getNumberOfVertices() - 1;
   if (lastVertex < 2)
   {
      CINFO<<MODULE<<"Intersect polygon num vertices is less than 3! This should never happen."<<endl;
      return;
   }
   if (intersectPoly[0] == intersectPoly[lastVertex])
   {
      intersectPoly.removeVertex(lastVertex);
      if (lastVertex < 3)
      {
         CWARN<<MODULE<<"Intersect polygon num vertices is less than 3! This should never happen."<<endl;
         return;
      }
   }

   double intersectArea = fabs(intersectPoly.area());
   if (intersectArea < min_area)
   {
      CINFO<<MODULE<<"Intersect area too small"<<endl;
      return;
   }

   if (config.diagnosticLevel(3) && m_annotatedRefImage)
      m_annotatedRefImage->annotateOverlap(intersectPoly);

   // Prepare to layout evenly-spaced search tiles across overlap bounding rect:
   ossimIrect aoiRect;
   intersectPoly.getBoundingRect(aoiRect);
   int side_length = (int) config.getParameter("featureSearchTileSize").asUint();
   int dx, dy; // tile width and height
   int aoiWidth = aoiRect.width();
   int aoiHeight = aoiRect.height();

   // Determine optimal X and Y-direction space between tiles (actually between corresponding
   // UL corners)
   int gridSize = (int) config.getParameter("searchGridSize").asUint();
   if (gridSize == 0)
   {
      // Check for overlap requirement to account for sampling kernel width:
      int marginSize = (int) config.getParameter("featureSearchTileMargin").asUint();

      // No gridsize indicates contiguous tiles (full coverage, usually for DEM generation ATPs:
      dx = side_length-marginSize;
      dy = side_length-marginSize;
   }
   else
   {
      for (int nx=gridSize; nx>1; nx--)
      {
         dx = (int) floor(aoiWidth/nx);
         if (dx >= side_length)
            break;
      }
      for (int ny=gridSize; ny>1; ny--)
      {
         dy = (int) floor(aoiHeight/ny);
         if (dy >= side_length)
            break;
      }
   }

   // Create tile rects now that the intervals are known:
   for (int y=aoiRect.ul().y; y<=aoiRect.lr().y-side_length; y+=dy)
   {
      for (int x=aoiRect.ul().x; x<=aoiRect.lr().x-side_length; x+=dx)
      {
         m_searchTileRects.push_back(ossimIrect (x, y, x+side_length-1, y+side_length-1));
      }
   }

   // Keep only those rects that are inside the overlap polygon. The layout was done over the
   // overlap's bounding rect, so some may spill outside the overlap:
   vector<ossimIrect>::iterator b = m_searchTileRects.begin();
   while (b != m_searchTileRects.end())
   {
      if(!intersectPoly.isRectWithin(*b))
         b = m_searchTileRects.erase(b);
      else
         b++;
   }

   // Don't bother with less than N boxes, just make the entire intersection a search patch:
   if (gridSize && (m_searchTileRects.size() < 4))
   {
      m_searchTileRects.clear();
      m_searchTileRects.push_back(m_aoiView);
   }

   if (config.diagnosticLevel(2))
   {
      for (int tileIdx=0; tileIdx<m_searchTileRects.size(); ++tileIdx)
         CINFO<<MODULE<<"m_searchTileRects["<<tileIdx<<"] = "<<m_searchTileRects[tileIdx]<<endl;

      if (config.diagnosticLevel(3) && m_annotatedRefImage)
         m_annotatedRefImage->annotateFeatureSearchTiles(m_searchTileRects);
   }
}

setBandWeights()

void ATP::AtpGenerator::setBandWeights ( const std::vector< double >  weights )

inline

When the input images are multiband, the bands must be combined into a single-band image.

Parameters

weights

Vector of normalized weights (0-1) applied to the input bands for combining.

Definition at line 54 of file AtpGenerator.h.

References m_bandWeights.

{ m_bandWeights = weights; }

setCmpImage()

void ATP::AtpGenerator::setCmpImage

(

std::shared_ptr< ossim::Image >

cmp_image

)

Definition at line 57 of file AtpGenerator.cpp.

References m_cmpImage.

{
   m_cmpImage = cmp_image;
}

setRefImage()

void ATP::AtpGenerator::setRefImage

(

std::shared_ptr< ossim::Image >

ref_image

)

Definition at line 51 of file AtpGenerator.cpp.

References m_refImage.

{
   m_refImage = ref_image;
}

writeTiePointList()

void ATP::AtpGenerator::writeTiePointList ( ostream &  out, const AtpList &  tpList  )

static

For engineering use.

Renders relevant tiepoint data to disk file(s). Optionally accepts an LUT file for raster-based remapping.

Parameters

outBaseName	This is the base name of possible multiple outputs, as needed by derived generator.
lutFile	Optional LUT file formatted as specified in ossimIndexToRgbLutFilter.

Returns

True if successful. For engineering use. Convenience method for outputting tiepoint list to output stream.

Definition at line 435 of file AtpGenerator.cpp.

{
   if (tpList.empty())
   {
      out<<"\nTiepoint list is empty.\n"<<endl;
   }
   else
   {
      for (size_t i=0; i<tpList.size(); ++i)
      {
         out<<"\ntiepoint["<<i<<"]: "<<tpList[i]<<endl;
      }
   }
}

Friends And Related Function Documentation

AtpTileSource

friend class AtpTileSource

friend

Definition at line 35 of file AtpGenerator.h.

Member Data Documentation

m_algorithm

Algorithm ATP::AtpGenerator::m_algorithm

protected

Definition at line 116 of file AtpGenerator.h.

Referenced by constructChain(), and initialize().

m_annotatedCmpImage

ossimRefPtr<AtpAnnotatedImage> ATP::AtpGenerator::m_annotatedCmpImage

Definition at line 92 of file AtpGenerator.h.

Referenced by generateTiePointList(), and initialize().

m_annotatedRefImage

ossimRefPtr<AtpAnnotatedImage> ATP::AtpGenerator::m_annotatedRefImage

Definition at line 91 of file AtpGenerator.h.

Referenced by generateTiePointList(), initialize(), and layoutSearchTileRects().

m_aoiView

ossimIrect ATP::AtpGenerator::m_aoiView

protected

Definition at line 125 of file AtpGenerator.h.

Referenced by initialize(), and layoutSearchTileRects().

m_atpTileSource

ossimRefPtr<AtpTileSource> ATP::AtpGenerator::m_atpTileSource

protected

Definition at line 119 of file AtpGenerator.h.

Referenced by generateTiePointList(), getAtpTileSource(), and initialize().

m_bandWeights

std::vector<double> ATP::AtpGenerator::m_bandWeights

protected

Definition at line 126 of file AtpGenerator.h.

Referenced by setBandWeights().

m_cmpChain

ossimRefPtr<ossimImageChain> ATP::AtpGenerator::m_cmpChain

protected

Definition at line 121 of file AtpGenerator.h.

Referenced by generateTiePointList(), getCmpChain(), getCmpFilename(), getCmpImageID(), and initialize().

m_cmpImage

std::shared_ptr<ossim::Image> ATP::AtpGenerator::m_cmpImage

protected

Definition at line 118 of file AtpGenerator.h.

Referenced by initialize(), and setCmpImage().

m_cmpIVT

ossimRefPtr<ossimImageViewProjectionTransform> ATP::AtpGenerator::m_cmpIVT

protected

Definition at line 123 of file AtpGenerator.h.

Referenced by getCmpIVT(), and initialize().

m_refChain

ossimRefPtr<ossimImageChain> ATP::AtpGenerator::m_refChain

protected

Definition at line 120 of file AtpGenerator.h.

Referenced by generateTiePointList(), getRefChain(), getRefFilename(), getRefImageID(), and initialize().

m_refEllipHgt

double ATP::AtpGenerator::m_refEllipHgt

protected

Definition at line 127 of file AtpGenerator.h.

Referenced by initialize().

m_refImage

std::shared_ptr<ossim::Image> ATP::AtpGenerator::m_refImage

protected

Definition at line 117 of file AtpGenerator.h.

Referenced by initialize(), and setRefImage().

m_refIVT

ossimRefPtr<ossimImageViewProjectionTransform> ATP::AtpGenerator::m_refIVT

protected

Definition at line 122 of file AtpGenerator.h.

Referenced by getRefIVT(), and initialize().

m_searchTileRects

std::vector<ossimIrect> ATP::AtpGenerator::m_searchTileRects

protected

Definition at line 128 of file AtpGenerator.h.

Referenced by generateTiePointList(), and layoutSearchTileRects().

m_viewGeom

ossimRefPtr<ossimImageGeometry> ATP::AtpGenerator::m_viewGeom

protected

Definition at line 124 of file AtpGenerator.h.

Referenced by constructChain(), and initialize().

s_referenceATP

std::shared_ptr< AutoTiePoint > ATP::AtpGenerator::s_referenceATP

staticprotected

Definition at line 129 of file AtpGenerator.h.

Referenced by initialize().

---

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

• AtpGenerator.h
• AtpGenerator.cpp