ossimImageGeometry Class Reference

Container class that holds both 2D transform and 3D projection information for an image Only one instance of this is alive per image. More...

#include <ossimImageGeometry.h>

Inheritance diagram for ossimImageGeometry:

Public Types
typedef std::vector< ossim_float64 >	DeltaParamList

Public Member Functions
	ossimImageGeometry ()
	Default constructor defaults to unity transform with no projection. More...
virtual	~ossimImageGeometry ()
	ossimImageGeometry (const ossimImageGeometry &copy_this)
	Copy constructor – performs a deep copy. More...
	ossimImageGeometry (ossim2dTo2dTransform *transform, ossimProjection *projection)
	Constructs with projection and transform objects available for referencing. More...
bool	open (const ossimFilename &image)
	Shortcut way of getting to an image's geometry when access to pixels is not needed. More...
void	rnToRn (const ossimDpt &inRnPt, ossim_uint32 inResolutionLevel, ossim_uint32 outResolutionLevel, ossimDpt &outRnPt) const
	rnToRn is a utility method that takes a rn resolution image point and maps it to the another rn resolution image point. More...
void	rnToFull (const ossimDpt &rnPt, ossim_uint32 resolutionLevel, ossimDpt &fullPt) const
	rnToFull is a utility method that takes a rn resolution image point and maps it to the full image point. More...
void	fullToRn (const ossimDpt &fullPt, ossim_uint32 resolutionLevel, ossimDpt &rnPt) const
	fullToRn is a utility method that takes a full image point and maps it to a rn resolution image point. More...
void	rnToWorld (const ossimDpt &rnPt, ossim_uint32 resolutionLevel, ossimGpt &wpt) const
	rnToWorld is a utility method that takes a rn resolution image point and maps it to the world point. More...
void	worldToRn (const ossimGpt &wpt, ossim_uint32 resolutionLevel, ossimDpt &rnPt) const
	worldToRn is a utility method that takes a world point allows one to transform all the way back to an rn point. More...
bool	localToWorld (const ossimDpt &local_pt, ossimGpt &world_pt) const
	Exposes the 3D projection from image to world coordinates. More...
bool	localToWorld (const ossimDrect &local_rect, ossimGrect &world_rect) const
bool	localToWorld (const ossimDpt &local_pt, const double &h_ellipsoid, ossimGpt &world_pt) const
	Exposes the 3D projection from image to world coordinates given a constant height above ellipsoid. More...
bool	worldToLocal (const ossimGpt &world_pt, ossimDpt &local_pt) const
	Exposes the 3D world-to-local image coordinate reverse projection. More...
bool	worldToLocal (const ossimGrect &world_rect, ossimDrect &local_rect) const
void	setTransform (ossim2dTo2dTransform *transform)
	Sets the transform to be used for local-to-full-image coordinate transformation. More...
void	setProjection (ossimProjection *projection)
	Sets the projection to be used for local-to-world coordinate transformation. More...
const ossim2dTo2dTransform *	getTransform () const
	Access methods for transform (may be NULL pointer). More...
ossim2dTo2dTransform *	getTransform ()
const ossimProjection *	getProjection () const
	Access methods for projection (may be NULL pointer). More...
ossimProjection *	getProjection ()
bool	isMapProjected () const
	Returns true if underlying projection is derived from ossimMapProjection. More...
const ossimMapProjection *	getAsMapProjection () const
ossimMapProjection *	getAsMapProjection ()
bool	hasProjection () const
	Returns TRUE if valid projection defined. More...
bool	hasTransform () const
	Returns TRUE if valid transform defined. More...
bool	isAffectedByElevation () const
	Returns TRUE if this geometry is sensitive to elevation. More...
ossimDpt	getMetersPerPixel () const
	Returns the GSD associated with this image in the active projection. More...
void	getMetersPerPixel (ossimDpt &gsd) const
	Get the ground sample distance(GSD) associated with this image in the active projection. More...
ossimDpt	getDegreesPerPixel () const
	Returns the resolution of this image in degrees/pixel. More...
void	getDegreesPerPixel (ossimDpt &dpp) const
	Get the resolution of this image in degrees/pixel. More...
bool	getCornerGpts (ossimGpt &ul, ossimGpt &ur, ossimGpt &lr, ossimGpt &ll) const
	Assigns the ossimGpts with the ground coordinates of the four corresponding image corner points. More...
void	getTiePoint (ossimGpt &tie, bool edge) const
	Get the latitude, longitude of the tie point. More...
void	getTiePoint (ossimDpt &tie, bool edge) const
	Get the easting, northing of the tie point. More...
std::ostream &	print (std::ostream &out) const
	Prints contents to output stream. More...
bool	operator== (const ossimImageGeometry &compare_this) const
	Returns TRUE if argument geometry has exactly the same instances for transform and m_projection. More...
const ossimImageGeometry &	operator= (const ossimImageGeometry &copy_this)
	Replaces any existing transform and projection with those in the copy_this object. More...
ossimDpt	decimationFactor (ossim_uint32 r_index) const
	Returns the decimation factor from R0 for the resolution level specified. More...
void	decimationFactor (ossim_uint32 r_index, ossimDpt &result) const
	Method to get the decimation factor for a given resolution level. More...
void	decimationFactors (std::vector< ossimDpt > &decimations) const
	Gets array of all decimation levels. More...
void	setDiscreteDecimation (const std::vector< ossimDpt > &decimation_list)
	Sets the decimation scheme to a discrete list of decimation factors. More...
ossim_uint32	getNumberOfDecimations () const
void	setImageSize (const ossimIpt &size)
const ossimIpt &	getImageSize () const
bool	getCrossesDateline () const
void	getImageEdgePoints (std::vector< ossimDpt > &result, ossim_uint32 partitions=25) const
void	calculatePolyBounds (ossimPolyArea2d &result, ossim_uint32 partitions=25) const
	This is the first stage implementation. More...
void	getBoundingRect (ossimIrect &bounding_rect) const
	Get the bounding rect of (0, 0) to (imageSize.x-1, imageSize.y-1). More...
void	getBoundingRect (ossimDrect &bounding_rect) const
void	getBoundingGroundRect (ossimGrect &bounding_grect) const
virtual ossimObject *	dup () const
	Creates a new instance of ossimImageGeometry with the same transform and projection. More...
virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)
	Attempts to initialize a transform and a projection given the KWL. More...
virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const
	Saves the transform (if any) and projection (if any) states to the KWL. More...
void	setTargetRrds (ossim_uint32 rrds)
	Set m_targetRrds data member. More...
ossim_uint32	getTargetRrds () const
void	applyScale (const ossimDpt &scale, bool recenterTiePoint)
	Changes the GSD and image size to reflect the scale provided. More...
virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const
virtual ossimAdjustableParameterInterface *	getAdjustableParameterInterface ()
virtual const ossimAdjustableParameterInterface *	getAdjustableParameterInterface () const
bool	computeImageToGroundPartialsWRTAdjParam (ossimDpt &result, const ossimGpt &gpt, ossim_uint32 idx, ossim_float64 paramDelta=1.0)
bool	computeImageToGroundPartialsWRTAdjParams (NEWMAT::Matrix &result, const ossimGpt &gpt, ossim_float64 paramDelta=1.0)
bool	computeImageToGroundPartialsWRTAdjParams (NEWMAT::Matrix &result, const ossimGpt &gpt, const DeltaParamList &deltas)
bool	computeGroundToImagePartials (NEWMAT::Matrix &result, const ossimGpt &gpt, const ossimDpt3d &deltaLlh)
bool	computeGroundToImagePartials (NEWMAT::Matrix &result, const ossimGpt &gpt)
ossim_float64	upIsUpAngle (const ossimDpt &) const
ossim_float64	upIsUpAngle () const
ossim_float64	northUpAngle () const
Public Member Functions inherited from ossimObject
	ossimObject ()
virtual	~ossimObject ()
virtual ossimString	getShortName () const
virtual ossimString	getLongName () const
virtual ossimString	getDescription () const
virtual ossimString	getClassName () 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 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

Protected Member Functions
void	undecimatePoint (const ossimDpt &rnPt, ossim_uint32 resolutionLevel, ossimDpt &outPt) const
	Method to back out decimation of a point. More...
void	decimatePoint (const ossimDpt &inPt, ossim_uint32 resolutionLevel, ossimDpt &rnPt) const
	Method to apply decimation of a point. More...
Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes
ossimRefPtr< ossim2dTo2dTransform >	m_transform
	Maintains local_image-to-full_image transformation. More...
ossimRefPtr< ossimProjection >	m_projection
	Maintains full_image-to-world_space transformation. More...
std::vector< ossimDpt >	m_decimationFactors
	List of decimation factors for R-levels. More...
ossimIpt	m_imageSize
ossim_uint32	m_targetRrds
	Target rrds for localToWorld and worldToLocal methods. More...

Detailed Description

Container class that holds both 2D transform and 3D projection information for an image Only one instance of this is alive per image.

This is the object that is returned from a call to ossimImageSource::getImageGeometry()

All images in OSSIM have at least a 2D transformation from local (relative to start of file) to full-image coordinates. Often this transform is an identity transform (does nothing). In some cases however, the first pixel in the file does not correspond to the upper-left pixel of the original, full image. Since the image's sensor model may only apply to full-image pixel coordinates, it is necessary to reference the latter when making computing ground point location. Also the local image could be a reduced resolution level of the full image. So the image could be a decimated sub image of the full image.

Additionally, images typically are accompanied by some form of metadata that defines the mapping from 2D image coordinates to 3D world coordinates. This mapping may be in the form of a map projection for orthorectified images, or a perspective projection such as a sensor model.

This object class maintains both 2D transform and 3D projection information for the associated image. This object will typically be created and "owned" by the image handler object. Therefore, only one copy per image will exist at a time.

Notes: 1) Shifts are relative to "full image space". So if you have a sub image from r2 the shift given to the transform should be relative to "full image space". 2) A decimation of 1.0 is the full image. This may or may not be r0 as r0 can be decimated.

Definition at line 65 of file ossimImageGeometry.h.

Member Typedef Documentation

DeltaParamList

typedef std::vector<ossim_float64> ossimImageGeometry::DeltaParamList

Definition at line 68 of file ossimImageGeometry.h.

Constructor & Destructor Documentation

ossimImageGeometry() [1/3]

ossimImageGeometry::ossimImageGeometry

(

)

Default constructor defaults to unity transform with no projection.

Definition at line 31 of file ossimImageGeometry.cpp.

References m_imageSize, and ossimIpt::makeNan().

:   
m_transform (0),
m_projection(0),
m_decimationFactors(0),
m_imageSize(),
m_targetRrds(0)
{
   m_imageSize.makeNan();
}

~ossimImageGeometry()

ossimImageGeometry::~ossimImageGeometry ( )

virtual

Definition at line 79 of file ossimImageGeometry.cpp.

{
   // Nothing to do
}

ossimImageGeometry() [2/3]

ossimImageGeometry::ossimImageGeometry

(

const ossimImageGeometry &

copy_this

)

Copy constructor – performs a deep copy.

This is needed when an imageSource in the chain needs to modify the image geometry at that point in the chain. The image geometry to the left of that imageSource must remain unchanged. This constructor permits that imageSource to maintain its own geometry based on the input geometry. All chained objects to the right would see this alternate geometry. See ossimScaleFilter for an example of this case.

Definition at line 49 of file ossimImageGeometry.cpp.

:
ossimObject(copy_this),
m_transform(copy_this.m_transform.valid()?(ossim2dTo2dTransform*)copy_this.m_transform->dup():((ossim2dTo2dTransform*)0)),
m_projection(copy_this.m_projection.valid()?(ossimProjection*)copy_this.m_projection->dup():(ossimProjection*)0),
m_decimationFactors(copy_this.m_decimationFactors),
m_imageSize(copy_this.m_imageSize),
m_targetRrds(copy_this.m_targetRrds)
{
}

ossimImageGeometry() [3/3]

ossimImageGeometry::ossimImageGeometry	(	ossim2dTo2dTransform *	transform,
		ossimProjection *	proj
	)

Constructs with projection and transform objects available for referencing.

Either pointer can be NULL – the associated mapping would be identity.

Definition at line 64 of file ossimImageGeometry.cpp.

References m_imageSize, and ossimIpt::makeNan().

:
ossimObject(),
m_transform(transform),
m_projection(proj),
m_decimationFactors(0),
m_imageSize(),
m_targetRrds(0)   
{
   m_imageSize.makeNan();
}

Member Function Documentation

applyScale()

void ossimImageGeometry::applyScale	(	const ossimDpt &	scale,
		bool	recenterTiePoint
	)

Changes the GSD and image size to reflect the scale provided.

Parameters

scale	The scale to be applied in x and y directions
recenterTiePoint	If true the will adjust the tie point by shifting the original tie to the upper left corner, applying scale, then shifting back by half of either the new theDeltaLat/lon or theMetersPerPixel depending on if underlying projection isGeographic.

Definition at line 1285 of file ossimImageGeometry.cpp.

References ossimMapProjection::applyScale(), ossimRefPtr< T >::get(), m_imageSize, m_projection, ossimIpt::x, ossimDpt::x, ossimIpt::y, and ossimDpt::y.

{
   if ((scale.x != 0.0) && (scale.y != 0.0))
   {
      ossimMapProjection* map_Proj = dynamic_cast<ossimMapProjection*>(m_projection.get());
      if ( map_Proj )
      {
         m_imageSize.x = ossim::round<ossim_int32>(m_imageSize.x / scale.x);
         m_imageSize.y = ossim::round<ossim_int32>(m_imageSize.y / scale.y);
         map_Proj->applyScale(scale, recenterTiePoint);
      }
   }
}

calculatePolyBounds()

void ossimImageGeometry::calculatePolyBounds	(	ossimPolyArea2d &	result,
		ossim_uint32	partitions = 25
	)		const

This is the first stage implementation.

It will determine if it crosses the dateline and do a special multi polygon for the ossimPolyArea2d result. For each edge of the image we walk "partitions" number of points. If we cross a dateline we calculate the crossing lat by using the parametric form of the equation. For example, if we are crossing from positive to negative that means we are going through the 180 degree lon location and we will need to solve for the paramtetric parameter t and plug back into the parametric equation to solve for the latitude. Basic form: start + (end-start)*t = 180. If we are coming from negative to positive then we solve basic form: start + (end-start)*t = -180. Where deltaPoint is (end-start)

              ossim_float64 t = ((180-start.x)/deltaPoint.x);
              ossim_float64 lat = (start.y+deltaPoint.y*t);

If we do not cross the dateline then if the image is affected by elevation (i.e. a sensor model) then we use the partitions parameter to calculate that number of partitions to sample along each edge.

if The geometry is not affected by elevation then we just use the corner points and ignore the partitions parameter

Parameters

poly	holds the resulting polygon. This could be a MultiPolygon depending on the geometry
partitions	These are the number of steps you want when walking the border.

Definition at line 1146 of file ossimImageGeometry.cpp.

References ossimPolyArea2d::add(), ossimPolyArea2d::clear(), ossimPolygon::clear(), ossimPolygon::clipToRect(), ossim_hdf5::crossesDateline(), getBoundingRect(), getCrossesDateline(), getImageEdgePoints(), ossimDrect::hasNans(), ossimGpt::isLatNan(), ossimGpt::isLonNan(), ossimPolyArea2d::isValid(), localToWorld(), ossimGpt::lond(), ossimDrect::midPoint(), ossimPolyArea2d::setToBufferedShape(), and ossim::sgn().

{
   std::vector<ossimDpt> points;
   std::vector<ossimGpt> gPoints;
   std::vector<ossimPolygon> polyList;
   ossimDrect imageRect;
   ossimGpt cg;
   ossimDrect worldRect(ossimDpt(-180,-90),
                        ossimDpt(-180,90),
                        ossimDpt(180,90),
                        ossimDpt(180,-90));
   getBoundingRect(imageRect);
   // bool affectedByElevation = isAffectedByElevation();
   bool crossesDateline     = getCrossesDateline();
   result.clear();
   if(imageRect.hasNans())
   {
      // error out
      return;
   }
   localToWorld(imageRect.midPoint(), cg);
   if(cg.isLatNan() || cg.isLonNan())
   {
      return;
   }
   ossim_int32 sgn = static_cast<ossim_int32>(
                                ossim::sgn(cg.lond()));

   getImageEdgePoints(points, partitions);
   if(points.empty())
   {
      return;
   }

   if(crossesDateline)
   {
      for(std::vector<ossimDpt>::const_iterator iter=points.begin(); 
          iter != points.end();
          ++iter)
      {
         ossimGpt testGpt;
         localToWorld(*iter, testGpt); 

         if(!testGpt.isLatNan()&&!testGpt.isLonNan())
         {
            bool needToShift =( (std::fabs(testGpt.lond()-cg.lond()) > 180.0)||
                                (std::fabs(testGpt.lond()>180.0)) ); 
            if(needToShift)
            {
               testGpt.lond(testGpt.lond()+sgn*360);
            }
            gPoints.push_back(testGpt);        
         }
      }
      // now clip the ground list to the full ground rect
      //
      ossimPolygon tempPoly(gPoints);
      tempPoly.clipToRect(polyList, worldRect);
      for(std::vector<ossimPolygon>::const_iterator iter = polyList.begin();
         iter!=polyList.end();++iter)
      {
         result.add(ossimPolyArea2d(*iter));
      }
      // now shift gpoints to the other side
      //
      for(std::vector<ossimGpt>::iterator iter = gPoints.begin();
          iter != gPoints.end();
          ++iter)
      {
        (*iter).lond((*iter).lond()+(-sgn*360));
      }
      tempPoly = gPoints;
      polyList.clear();
      tempPoly.clipToRect(polyList, worldRect);
      for(std::vector<ossimPolygon>::const_iterator iter = polyList.begin();
         iter!=polyList.end();++iter)
      {
         result.add(ossimPolyArea2d(*iter));
      }
   }
   else
   {
      for(std::vector<ossimDpt>::const_iterator iter=points.begin(); 
          iter != points.end();++iter)
      {
         ossimGpt testGpt;
         localToWorld(*iter, testGpt); 

         if(!testGpt.isLatNan()&&!testGpt.isLonNan())
         {
            gPoints.push_back(testGpt);        
         }
      }
      result.add(ossimPolygon(gPoints));
   }
   if(!result.isValid()) result.setToBufferedShape();
}

computeGroundToImagePartials() [1/2]

bool ossimImageGeometry::computeGroundToImagePartials	(	NEWMAT::Matrix &	result,
		const ossimGpt &	gpt,
		const ossimDpt3d &	deltaLlh
	)

Definition at line 1400 of file ossimImageGeometry.cpp.

References ossimGpt::datum(), DEG_PER_RAD, getProjection(), ossimGpt::height(), ossim::isnan(), ossimGpt::latd(), ossimGpt::lond(), worldToLocal(), ossimDpt3d::x, ossimDpt::x, ossimDpt3d::y, ossimDpt::y, and ossimDpt3d::z.

Referenced by computeGroundToImagePartials(), and ossimPointObservation::getObjSpacePartials().

{
   if(!getProjection()) return false;
   ossimDpt p1;
   ossimDpt p2;
   
   ossimDpt deltaWithRespectToLon;
   ossimDpt deltaWithRespectToLat;
   ossimDpt deltaWithRespectToH;
   ossim_float64 h = ossim::isnan(gpt.height())?0.0:gpt.height();

   // do the change in lon first for the dx, dy
   //
   worldToLocal(ossimGpt(gpt.latd(), gpt.lond()+deltaLlh.x, h, gpt.datum()), p1);
   worldToLocal(ossimGpt(gpt.latd(), gpt.lond()-deltaLlh.x, h, gpt.datum()), p2);
   
   double den = 0.5/deltaLlh.x; // this is the same as dividing by 2*delta
   deltaWithRespectToLon = (p2-p1)*den;
   
    
   // do the change in lat for the dx, dy
   //
   worldToLocal(ossimGpt(gpt.latd()+deltaLlh.y, gpt.lond(), h, gpt.datum()), p1);
   worldToLocal(ossimGpt(gpt.latd()-deltaLlh.y, gpt.lond(), h, gpt.datum()), p2);
   
   den = 0.5/deltaLlh.y; // this is the same as dividing by 2*delta
   deltaWithRespectToLat = (p2-p1)*den;

   
   // do the change in height first for the dx, dy
   //
   worldToLocal(ossimGpt(gpt.latd(), gpt.lond(), h+deltaLlh.z, gpt.datum()), p1);
   worldToLocal(ossimGpt(gpt.latd(), gpt.lond(), h-deltaLlh.z, gpt.datum()), p2);
   
   den = 0.5/deltaLlh.z; // this is the same as dividing by 2*delta
   deltaWithRespectToH = (p2-p1)*den;
   
   
   result = NEWMAT::Matrix(3,2);
   
   // set the matrix
   //
   result[1][0] = deltaWithRespectToLon.x*DEG_PER_RAD; 
   result[1][1] = deltaWithRespectToLon.y*DEG_PER_RAD; 
   result[0][0] = deltaWithRespectToLat.x*DEG_PER_RAD; 
   result[0][1] = deltaWithRespectToLat.y*DEG_PER_RAD; 
   result[2][0] = deltaWithRespectToH.x; 
   result[2][1] = deltaWithRespectToH.y; 
   
   
   return true; 
}

computeGroundToImagePartials() [2/2]

bool ossimImageGeometry::computeGroundToImagePartials	(	NEWMAT::Matrix &	result,
		const ossimGpt &	gpt
	)

Definition at line 1455 of file ossimImageGeometry.cpp.

References computeGroundToImagePartials(), FLT_EPSILON, getMetersPerPixel(), ossimDpt::length(), and ossimGpt::metersPerDegree().

{
   ossimDpt mpp = getMetersPerPixel();
   ossimGpt originPoint;
   ossim_float64 len = mpp.length();
   
   if(len > FLT_EPSILON)
   {
      ossim_float64 delta = originPoint.metersPerDegree().length();
   
      delta = len/delta;
      
      return computeGroundToImagePartials(result, 
                                          gpt, 
                                          ossimDpt3d(delta,
                                                     delta,
                                                     len));
   }
   
   return false;
}

computeImageToGroundPartialsWRTAdjParam()

bool ossimImageGeometry::computeImageToGroundPartialsWRTAdjParam	(	ossimDpt &	result,
		const ossimGpt &	gpt,
		ossim_uint32	idx,
		ossim_float64	paramDelta = 1.0
	)

Definition at line 1310 of file ossimImageGeometry.cpp.

References getAdjustableParameterInterface(), ossimAdjustableParameterInterface::getNumberOfAdjustableParameters(), ossimAdjustableParameterInterface::getParameterCenter(), ossimDpt::hasNans(), ossimAdjustableParameterInterface::setParameterCenter(), and worldToLocal().

Referenced by computeImageToGroundPartialsWRTAdjParams().

{
   double den = 0.5/paramDelta; // this is the same as dividing by 2*delta
   
   result = ossimDpt(0.0,0.0);
   ossimAdjustableParameterInterface* adjustableParamInterface = getAdjustableParameterInterface();
   
   if(!adjustableParamInterface) return false;
   if(idx >= adjustableParamInterface->getNumberOfAdjustableParameters()) return false;
   
   ossimDpt p1, p2;
   // double middle = adjustableParamInterface->getAdjustableParameter(idx);
   double middle = adjustableParamInterface->getParameterCenter(idx);
   
   //set parm to high value
   // adjustableParamInterface->setAdjustableParameter(idx, middle + paramDelta, true);
   adjustableParamInterface->setParameterCenter(idx, middle + paramDelta, true);
   worldToLocal(gpt, p1);
   
   //set parm to low value and gte difference
   // adjustableParamInterface->setAdjustableParameter(idx, middle - paramDelta, true);
   adjustableParamInterface->setParameterCenter(idx, middle - paramDelta, true);
   worldToLocal(gpt, p2);
   
   //get partial derivative
   result = (p2-p1)*den;
   
   //reset param
   // adjustableParamInterface->setAdjustableParameter(idx, middle, true);
   adjustableParamInterface->setParameterCenter(idx, middle, true);
   
   return !result.hasNans();
}

computeImageToGroundPartialsWRTAdjParams() [1/2]

bool ossimImageGeometry::computeImageToGroundPartialsWRTAdjParams	(	NEWMAT::Matrix &	result,
		const ossimGpt &	gpt,
		ossim_float64	paramDelta = 1.0
	)

Definition at line 1347 of file ossimImageGeometry.cpp.

References computeImageToGroundPartialsWRTAdjParam(), getAdjustableParameterInterface(), ossimAdjustableParameterInterface::getNumberOfAdjustableParameters(), ossimDpt::x, and ossimDpt::y.

Referenced by ossimPointObservation::getParameterPartials().

{
   ossimAdjustableParameterInterface* adjustableParamInterface = getAdjustableParameterInterface();
   
   ossim_uint32 nAdjustables = adjustableParamInterface->getNumberOfAdjustableParameters();
   
   ossim_uint32 idx = 0;
   
   result = NEWMAT::Matrix(nAdjustables, 2);
   for(;idx < nAdjustables; ++idx)
   {
      ossimDpt paramResResult;
      computeImageToGroundPartialsWRTAdjParam(paramResResult,
                                                        gpt,
                                                        idx,
                                                        paramDelta);
      result[idx][0] = paramResResult.x;
      result[idx][1] = paramResResult.y;
   }
   
   return true;
}

computeImageToGroundPartialsWRTAdjParams() [2/2]

bool ossimImageGeometry::computeImageToGroundPartialsWRTAdjParams	(	NEWMAT::Matrix &	result,
		const ossimGpt &	gpt,
		const DeltaParamList &	deltas
	)

Definition at line 1372 of file ossimImageGeometry.cpp.

References computeImageToGroundPartialsWRTAdjParam(), getAdjustableParameterInterface(), ossimAdjustableParameterInterface::getNumberOfAdjustableParameters(), ossimDpt::x, and ossimDpt::y.

{
   
   ossimAdjustableParameterInterface* adjustableParamInterface = getAdjustableParameterInterface();
   
   ossim_uint32 nAdjustables = adjustableParamInterface->getNumberOfAdjustableParameters();
   
   if(nAdjustables != deltas.size()) return false;
   ossim_uint32 idx = 0;
   
   result = NEWMAT::Matrix(nAdjustables, 2);
   for(;idx < nAdjustables; ++idx)
   {
      ossimDpt paramResResult;
      computeImageToGroundPartialsWRTAdjParam(paramResResult,
                                                        gpt,
                                                        idx,
                                                        deltas[idx]);
      // ROWxCOL
      result[idx][0] = paramResResult.x;
      result[idx][1] = paramResResult.y;
   }
   
   return true;
}

decimatePoint()

void ossimImageGeometry::decimatePoint ( const ossimDpt &  inPt, ossim_uint32  resolutionLevel, ossimDpt &  rnPt  ) const

protected

Method to apply decimation of a point.

Parameters

inPt	Is a point with no decimation.
resolutionLevel	Is the resolution of the point rnPt.
rnPt	Is the result of the transform

Definition at line 948 of file ossimImageGeometry.cpp.

References decimationFactor(), ossimDpt::hasNans(), ossimDpt::makeNan(), ossimDpt::x, and ossimDpt::y.

Referenced by fullToRn(), and rnToRn().

{
   
   // Apply the decimation.
   ossimDpt decimation = decimationFactor(resolutionLevel);

   //---
   // If no nans and one or both of the factors is not 1.0 decimation should
   // be applied.
   //---
   if ( (decimation.x != 1.0) || (decimation.y != 1.0) ) 
   {
      if ( !decimation.hasNans() && !inPt.hasNans() )
      {
         rnPt.x = inPt.x * decimation.x;
         rnPt.y = inPt.y * decimation.y;
      }
      else
      {
        rnPt.makeNan(); 
      }
   }
   else
   {
      rnPt = inPt; // No decimation.
   }
}

decimationFactor() [1/2]

ossimDpt ossimImageGeometry::decimationFactor

(

ossim_uint32

r_index

)

const

Returns the decimation factor from R0 for the resolution level specified.

For r_index=0, the decimation factor is by definition 1.0. For the non-discrete case, r_index=1 returns a decimation of 0.5. If the vector of discrete decimation factors (m_decimationFactors) is empty, the factor will be computed as f=1/2^n

Definition at line 489 of file ossimImageGeometry.cpp.

References m_decimationFactors, and size.

Referenced by decimatePoint(), ossimMrSidReader::getDecimationFactor(), ossimMrSidReader::getTile(), ossimImageViewProjectionTransform::imageToView(), operator==(), undecimatePoint(), and ossimImageViewProjectionTransform::viewToImage().

{
   ossim_uint32 size = (ossim_uint32)m_decimationFactors.size();
   if (size)
   {
      if (r_index < size)
      {
         return m_decimationFactors[r_index];
      }
      
      // Return the last defined decimation if the index requested exceeds list size:
      return m_decimationFactors[size-1];
   }

   // Compute the decimation factor:
   ossim_float64 factor = 1.0/(ossim_float64)(1 << r_index);
   
   return ossimDpt(factor, factor);
}

decimationFactor() [2/2]

void ossimImageGeometry::decimationFactor	(	ossim_uint32	r_index,
		ossimDpt &	result
	)		const

Method to get the decimation factor for a given resolution level.

If the array of decimations is not initialized by owner, the default is: r_index=0 is by definition 1.0. r_index=n Where n is some level the factor will be computed as f=1/2^n.

Parameters

resLevel	Reduced resolution set for requested decimation.
result	ossimDpt to initialize with requested decimation.

Definition at line 509 of file ossimImageGeometry.cpp.

References m_decimationFactors, ossimDpt::x, and ossimDpt::y.

{
   const ossim_uint32 SIZE = (ossim_uint32)m_decimationFactors.size();
   if (SIZE)
   {
      if (r_index < SIZE)
      {
         result = m_decimationFactors[r_index];
      }
      else
      {
         //---
         // Return the last defined decimation if the index requested
         // exceeds list size:
         //---
         result = m_decimationFactors[SIZE-1];
      }
   }
   else
   {
      // Compute the decimation factor:
      ossim_float64 factor = 1.0/(ossim_float64)(1 << r_index);
      result.x = factor;
      result.y = factor;
   }
}

decimationFactors()

void ossimImageGeometry::decimationFactors

(

std::vector< ossimDpt > &

decimations

)

const

Gets array of all decimation levels.

Parameters

decimations

Array to initialiaze.

Definition at line 537 of file ossimImageGeometry.cpp.

References m_decimationFactors.

Referenced by ossimMrSidReader::getDecimationFactors().

{
   decimations = m_decimationFactors;
}

dup()

virtual ossimObject* ossimImageGeometry::dup ( ) const

inlinevirtual

Creates a new instance of ossimImageGeometry with the same transform and projection.

Overrides base-class version requiring loadState() and saveState() (not implemented yet)

Reimplemented from ossimObject.

Definition at line 378 of file ossimImageGeometry.h.

Referenced by ossimEquationUtil::initInputSources().

{ return new ossimImageGeometry(*this); }

fullToRn()

void ossimImageGeometry::fullToRn	(	const ossimDpt &	fullPt,
		ossim_uint32	resolutionLevel,
		ossimDpt &	rnPt
	)		const

fullToRn is a utility method that takes a full image point and maps it to a rn resolution image point.

Parameters

fullPt	Is a point in full image space.
resolutionLevel	Is the resolution of the point rnPt. A value of 0 is the local image.
fullPt	Is the result of the transform

Definition at line 135 of file ossimImageGeometry.cpp.

References decimatePoint(), ossim2dTo2dTransform::inverse(), m_transform, and ossimRefPtr< T >::valid().

Referenced by ossimImageViewProjectionTransform::imageToView(), ossimImageViewProjectionTransform::viewToImage(), and worldToLocal().

{
   // Apply shift/rotation.
   ossimDpt localPt;
   if (m_transform.valid())
   {
      m_transform->inverse(fullPt, localPt);
   }
   else
   {
      localPt = fullPt; // No transform (shift/rotation)
   }

   // Apply the decimation.
   decimatePoint(localPt, resolutionLevel, rnPt);
}

getAdjustableParameterInterface() [1/2]

ossimAdjustableParameterInterface * ossimImageGeometry::getAdjustableParameterInterface ( )

virtual

Definition at line 1299 of file ossimImageGeometry.cpp.

References getProjection(), and PTR_CAST.

Referenced by ossimPointObservation::addMeasurement(), computeImageToGroundPartialsWRTAdjParam(), computeImageToGroundPartialsWRTAdjParams(), ossimAutRegUtil::saveAdjustment(), and ossimAdjustmentExecutive::updateParameters().

{
   return PTR_CAST(ossimAdjustableParameterInterface, getProjection());
}

getAdjustableParameterInterface() [2/2]

const ossimAdjustableParameterInterface * ossimImageGeometry::getAdjustableParameterInterface ( ) const

virtual

Definition at line 1304 of file ossimImageGeometry.cpp.

References getProjection().

{
   return dynamic_cast<const ossimAdjustableParameterInterface*>(getProjection());

}

getAsMapProjection() [1/2]

const ossimMapProjection* ossimImageGeometry::getAsMapProjection ( ) const

inline

Returns

const ossimMapProjection* or NULL if projection not set or not derived from ossimMapProjection.

Definition at line 179 of file ossimImageGeometry.h.

Referenced by ossimGpkgWriter::addLevels(), ossimGmlSupportData::configureGmljp2V2(), ATP::AtpGenerator::constructChain(), ossimChipProcTool::getChip(), ossimGmlSupportData::getGeoBounds(), ossimGpkgWriter::getGsd(), ossimPdfWriter::getLgiDictDatum(), ossimPdfWriter::getLgiDictDescription(), ossimPdfWriter::getLgiDictProjectionType(), ossimGmlSupportData::getMapOrigin(), ossimGpkgWriter::getNewOutputProjection(), ossimGmlSupportData::initialize(), ossimPointCloudTool::loadPC(), ossimPointCloudTool::setGSD(), and ossimGpkgWriter::writeEntry().

   { return dynamic_cast<const ossimMapProjection*>( m_projection.get() ); }

getAsMapProjection() [2/2]

ossimMapProjection* ossimImageGeometry::getAsMapProjection ( )

inline

Returns

ossimMapProjection* or NULL if projection not set or not derived from ossimMapProjection.

Definition at line 186 of file ossimImageGeometry.h.

   { return dynamic_cast<ossimMapProjection*>( m_projection.get() ); }

getBoundingGroundRect()

void ossimImageGeometry::getBoundingGroundRect

(

ossimGrect &

bounding_grect

)

const

Definition at line 1275 of file ossimImageGeometry.cpp.

References getCornerGpts(), ossimGpt::lat, ossimGpt::lon, ossimGrect::lr(), and ossimGrect::ul().

Referenced by getTiePoint(), and ossimChipProcTool::setAoiToInputs().

{
   ossimGpt ul, ur, lr, ll;
   getCornerGpts(ul, ur, lr, ll);
   bounding_grect.ul().lat = ul.lat > ur.lat ? ul.lat : ur.lat;
   bounding_grect.ul().lon = ul.lon < ll.lon ? ul.lon : ll.lon;
   bounding_grect.lr().lat = ll.lat < lr.lat ? ll.lat : lr.lat;
   bounding_grect.lr().lon = ur.lon > lr.lon ? ur.lon : lr.lon;
}

getBoundingRect() [1/2]

void ossimImageGeometry::getBoundingRect

(

ossimIrect &

bounding_rect

)

const

Get the bounding rect of (0, 0) to (imageSize.x-1, imageSize.y-1).

Relies on image size being initialized.

Parameters

bounding_rect

Initialized by this method. Will do a ossimIrect::makeNan() if the image size is not initialized.

Definition at line 1245 of file ossimImageGeometry.cpp.

References ossimIpt::hasNans(), m_imageSize, ossimIrect::makeNan(), ossimIrect::set_lrx(), ossimIrect::set_lry(), ossimIrect::set_ulx(), ossimIrect::set_uly(), ossimIpt::x, and ossimIpt::y.

Referenced by calculatePolyBounds(), ATP::AtpGenerator::constructChain(), ossimPotraceTool::convertToBitmap(), ossimInfo::getCenterGround(), getCrossesDateline(), ossimInfo::getImageBounds(), getImageEdgePoints(), ossimInfo::getImg2grd(), ossimRangeDomeTileSource::getTile(), ossimImageViewProjectionTransform::getViewSegments(), ossimImageRenderer::initializeBoundingRects(), northUpAngle(), ossimPotraceTool::transformLineStrings(), and upIsUpAngle().

{
   if (m_imageSize.hasNans())
   {
      bounding_rect.makeNan();
   }
   else
   {
      bounding_rect.set_ulx(0);
      bounding_rect.set_uly(0);
      bounding_rect.set_lrx(m_imageSize.x-1);
      bounding_rect.set_lry(m_imageSize.y-1);
   }
}

getBoundingRect() [2/2]

void ossimImageGeometry::getBoundingRect

(

ossimDrect &

bounding_rect

)

const

Definition at line 1260 of file ossimImageGeometry.cpp.

References ossimIpt::hasNans(), m_imageSize, ossimDrect::makeNan(), ossimDrect::set_lrx(), ossimDrect::set_lry(), ossimDrect::set_ulx(), ossimDrect::set_uly(), ossimIpt::x, and ossimIpt::y.

{
   if (m_imageSize.hasNans())
   {
      bounding_rect.makeNan();
   }
   else
   {
      bounding_rect.set_ulx(0);
      bounding_rect.set_uly(0);
      bounding_rect.set_lrx(m_imageSize.x-1);
      bounding_rect.set_lry(m_imageSize.y-1);
   }
}

getCornerGpts()

bool ossimImageGeometry::getCornerGpts	(	ossimGpt &	ul,
		ossimGpt &	ur,
		ossimGpt &	lr,
		ossimGpt &	ll
	)		const

Assigns the ossimGpts with the ground coordinates of the four corresponding image corner points.

Returns the ossimGpts for the four image corner points.

Returns true if points are valid.

Definition at line 714 of file ossimImageGeometry.cpp.

References localToWorld(), m_imageSize, status, ossimIpt::x, and ossimIpt::y.

Referenced by getBoundingGroundRect(), ossimTiledElevationDatabase::getBoundingRect(), ossimImageElevationHandler::open(), and ossimNitfWriterBase::writeGeometry().

{
   ossimDpt iul (0,0);
   ossimDpt iur (m_imageSize.x-1, 0);
   ossimDpt ilr (m_imageSize.x-1, m_imageSize.y-1);
   ossimDpt ill (0, m_imageSize.y-1);

   bool status = true;

   status &= localToWorld(iul, gul);
   status &= localToWorld(iur, gur);
   status &= localToWorld(ilr, glr);
   status &= localToWorld(ill, gll);

   return status;
}

getCrossesDateline()

bool ossimImageGeometry::getCrossesDateline

(

)

const

Definition at line 1027 of file ossimImageGeometry.cpp.

References ossimGpt::changeDatum(), ossimGpt::datum(), getBoundingRect(), ossimGpt::isLonNan(), ossimIrect::ll(), localToWorld(), ossimGpt::lond(), ossimIrect::lr(), ossimIrect::midPoint(), ossimIrect::ul(), and ossimIrect::ur().

Referenced by calculatePolyBounds(), ossimImageRenderer::initializeBoundingRects(), ossimImageViewProjectionTransform::initializeDatelineCrossing(), and ossimImageViewProjectionTransform::initializeViewSize().

{
   bool result = false;
   ossimIrect rect;
   ossimGpt ul;
   ossimGpt ur;
   ossimGpt lr;
   ossimGpt ll;
   ossimGpt center;
   ossimGpt wgs84;

   getBoundingRect(rect);

   localToWorld(rect.ul(), ul);
   localToWorld(rect.ur(), ur);
   localToWorld(rect.lr(), lr);
   localToWorld(rect.ll(), ll);
   localToWorld(rect.midPoint(), center);

   if(ul.isLonNan()||ur.isLonNan()||lr.isLonNan()||ll.isLonNan()||center.isLonNan())
   {
      return result;
   }
   else
   {
      ul.changeDatum(wgs84.datum());
      ur.changeDatum(wgs84.datum());
      lr.changeDatum(wgs84.datum());
      ll.changeDatum(wgs84.datum());
      center.changeDatum(wgs84.datum());
      result = ( (fabs(center.lond()-ul.lond()) > 180.0) ||
                 (fabs(center.lond()-ur.lond()) > 180.0) ||
                 (fabs(center.lond()-lr.lond()) > 180.0) ||
                 (fabs(center.lond()-ll.lond()) > 180.0) ||
                 (fabs(ul.lond()) > 180.0) ||
                 (fabs(ur.lond()) > 180.0) ||
                 (fabs(lr.lond()) > 180.0) ||
                 (fabs(ll.lond()) > 180.0)

                 );
   }

   return result; 
}

getDegreesPerPixel() [1/2]

ossimDpt ossimImageGeometry::getDegreesPerPixel

(

)

const

Returns the resolution of this image in degrees/pixel.

Note that this only makes sense if there is a projection associated with the image. Returns NaNs if no projection defined.

Definition at line 389 of file ossimImageGeometry.cpp.

Referenced by ossimChipProcTool::computeAdjustedViewFromGrect(), ossimGmlSupportData::getGeoOrigin(), ossimGpkgWriter::getGsd(), ossimLasReader::getImageGeometry(), and ossimPdfWriter::getLgiDictCtm().

{
   ossimDpt dpp;
   getDegreesPerPixel(dpp);
   return dpp;
}

getDegreesPerPixel() [2/2]

void ossimImageGeometry::getDegreesPerPixel

(

ossimDpt &

dpp

)

const

Get the resolution of this image in degrees/pixel.

Degrees/pixel is taken from projection if there is no transform set; else, taken from three localToWorld calls (one pixel apart) at the image center. If the projection or the image size(if needed) is not set the point will be set to NaNs.

Parameters

dpp	Point to intialize with degrees/pixel.

Note

This only makes sense if there is a projection associated with the image.

Result should be checked for NaNs after call.

Definition at line 401 of file ossimImageGeometry.cpp.

References ossimRefPtr< T >::get(), ossimMapProjection::getDecimalDegreesPerPixel(), ossimIpt::hasNans(), ossimGpt::height(), ossimDpt::lat, ossimGpt::lat, localToWorld(), ossimDpt::lon, ossimGpt::lon, m_imageSize, m_projection, m_transform, ossimDpt::makeNan(), ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

{
   const ossimMapProjection *map_proj = dynamic_cast<const ossimMapProjection *>(m_projection.get());
   if (map_proj && !m_transform.valid())
   {
      // No transform present, so simply query the projection for resolution:
      dpp = map_proj->getDecimalDegreesPerPixel();
   }
   else if (m_projection.valid() && (m_imageSize.hasNans() == false))
   {
      // A transform is involved, so need to use localToWorld call below:
      ossimDpt pL0 (m_imageSize/2);
      ossimDpt pLx (pL0.x+1, pL0.y);
      ossimDpt pLy (pL0.x, pL0.y+1);
      ossimGpt g0, gx, gy;

      localToWorld(pL0, g0);
      localToWorld(pLx, g0.height(), gx);
      localToWorld(pLy, g0.height(), gy);

      // Compute horizontal distance for one pixel:
      double dlatx = std::fabs(g0.lat - gx.lat);
      double dlaty = std::fabs(g0.lat - gy.lat);
      double dlonx = std::fabs(g0.lon - gx.lon);
      double dlony = std::fabs(g0.lon - gy.lon);
      dpp.lat = sqrt(dlatx*dlatx + dlaty*dlaty);
      dpp.lon = sqrt(dlonx*dlonx + dlony*dlony);
   }
   else
   {
      // This object is not fully initialized:
      dpp.makeNan();
   }
}

getImageEdgePoints()

void ossimImageGeometry::getImageEdgePoints	(	std::vector< ossimDpt > &	result,
		ossim_uint32	partitions = 25
	)		const

Definition at line 1072 of file ossimImageGeometry.cpp.

References ossimDrect::expand(), getBoundingRect(), ossimDrect::hasNans(), ossimDrect::ll(), ossimDrect::lr(), ossimDrect::ul(), ossimDrect::ur(), ossimDpt::x, and ossimDpt::y.

Referenced by calculatePolyBounds(), and ossimImageViewProjectionTransform::getViewSegments().

{

   ossimDrect imageRect;
   getBoundingRect(imageRect);
   if(imageRect.hasNans())
   {
      // error out
      return;
   }

   // Make edge to edge.
   imageRect.expand( ossimDpt(0.5, 0.5) );
   
   result.clear();
   // First get the image points we will be transforming
   if(partitions > 2)
   {
      ossimDpt uli = imageRect.ul();
      ossimDpt uri = imageRect.ur();
      ossimDpt lri = imageRect.lr();
      ossimDpt lli = imageRect.ll();

      ossim_float32 stepSize = partitions;
      ossimDpt deltaUpper = (uri-uli)*(1.0/stepSize);
      ossimDpt deltaRight = (lri-uri)*(1.0/stepSize);
      ossimDpt deltaLower = (lli-lri)*(1.0/stepSize);
      ossimDpt deltaLeft  = (uli-lli)*(1.0/stepSize);

      ossimDpt p;
      ossim_int32 idx = 0;
      ossimDpt initialPoint= uli;
      for(idx = 0; idx < stepSize;++idx)
      {
         result.push_back(initialPoint);
         initialPoint.x+=deltaUpper.x;
         initialPoint.y+=deltaUpper.y;
      }
      initialPoint= uri;
      for(idx = 0; idx < stepSize;++idx)
      {
         result.push_back(initialPoint);
         initialPoint.x+=deltaRight.x;
         initialPoint.y+=deltaRight.y;
      }

      initialPoint= lri;
      for(idx = 0; idx < stepSize;++idx)
      {
         result.push_back(initialPoint);
         initialPoint.x+=deltaLower.x;
         initialPoint.y+=deltaLower.y;
      }

      initialPoint= lli;
      for(idx = 0; idx < stepSize;++idx)
      {
         result.push_back(initialPoint);
         initialPoint.x+=deltaLeft.x;
         initialPoint.y+=deltaLeft.y;
      }
   }
   else // If not enough partitions then we will just use the corners
   {
      result.resize(4);

      result[0] = imageRect.ul();
      result[1] = imageRect.ur();
      result[2] = imageRect.lr();
      result[3] = imageRect.ll();
   }   

}

getImageSize()

const ossimIpt& ossimImageGeometry::getImageSize ( ) const

inline

Definition at line 330 of file ossimImageGeometry.h.

Referenced by ossimChipProcTool::execute(), ossimRangeDomeTileSource::getNumberOfLines(), ossimPointCloudImageHandler::getNumberOfLines(), ossimRangeDomeTileSource::getNumberOfSamples(), ossimPointCloudImageHandler::getNumberOfSamples(), ossimImageViewProjectionTransform::initializeViewSize(), ossimImageHandler::initImageParameters(), and ossimGeomFileWriter::writeFile().

   {
      return m_imageSize;
   }

getMetersPerPixel() [1/2]

ossimDpt ossimImageGeometry::getMetersPerPixel

(

)

const

Returns the GSD associated with this image in the active projection.

Returns GSD in meters in x- and y-directions.

Note that this only makes sense if there is a projection associated with the image. Returns NaNs if no projection defined.

This may not be the same as the member projection's stated resolution since a transform may be involved.

Definition at line 343 of file ossimImageGeometry.cpp.

Referenced by computeGroundToImagePartials(), ATP::AtpGenerator::constructChain(), ossimGpkgWriter::getGsd(), ossimLasReader::getImageGeometry(), ossimInfo::getImageGeometryInfo(), ossimGmlSupportData::getMapOrigin(), ossimGeoAnnotationEllipseObject::getWidthHeightInPixels(), ossimSFIMFusion::initialize(), ossimImageToPlaneNormalFilter::initialize(), ossimTilingRect::initializeBase(), ossimImageRenderer::initializeBoundingRects(), ossimEquationUtil::initInputSources(), ossimImageElevationHandler::open(), saveState(), and ossimGeneralRasterElevHandler::setFilename().

{
   ossimDpt gsd;
   getMetersPerPixel(gsd);
   return gsd;
}

getMetersPerPixel() [2/2]

void ossimImageGeometry::getMetersPerPixel

(

ossimDpt &

gsd

)

const

Get the ground sample distance(GSD) associated with this image in the active projection.

Computes GSD in meters in x- and y-directions.

GSD is taken from projection if there is no transform set; else, taken from three localToWorld calls (one pixel apart) at the image center. If the projection or the image size(if needed) is not set the point will be set to NaNs.

Parameters

gsd	Point to intialize with GSD.

Note

This only makes sense if there is a projection associated with the image.

Result should be checked for NaNs after call.

This may not be the same as the member projection's stated resolution since a transform may be involved.

Definition at line 354 of file ossimImageGeometry.cpp.

References ossimGpt::distanceTo(), ossimProjection::getMetersPerPixel(), ossimIpt::hasNans(), ossimGpt::height(), localToWorld(), m_imageSize, m_projection, m_transform, ossimDpt::makeNan(), ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

{
   if (m_projection.valid() && !m_transform.valid())
   {
      // No transform present, so simply query the projection for GSD:
      gsd = m_projection->getMetersPerPixel();
   }
   else if (m_projection.valid() && (m_imageSize.hasNans() == false))
   {
      // A transform is involved, so need to use localToWorld call below:
      ossimDpt pL0 (m_imageSize/2);
      ossimDpt pLx (pL0.x+1, pL0.y);
      ossimDpt pLy (pL0.x, pL0.y+1);
      ossimGpt g0, gx, gy;

      localToWorld(pL0, g0);
      localToWorld(pLx, g0.height(), gx);
      localToWorld(pLy, g0.height(), gy);

      // Compute horizontal distance for one pixel:
      gsd.x = g0.distanceTo(gx);
      gsd.y = g0.distanceTo(gy);
   }
   else
   {
      // This object is not fully initialized:
      gsd.makeNan();
   }
}

getNumberOfDecimations()

ossim_uint32 ossimImageGeometry::getNumberOfDecimations ( ) const

inline

Returns

The number of decimation factors

Definition at line 322 of file ossimImageGeometry.h.

Referenced by ossimMrSidReader::getDecimationFactors(), and ossimImageHandler::initImageParameters().

   {
      return (ossim_uint32)m_decimationFactors.size();
   }

getProjection() [1/2]

const ossimProjection* ossimImageGeometry::getProjection ( ) const

inline

Access methods for projection (may be NULL pointer).

Definition at line 165 of file ossimImageGeometry.h.

Referenced by ossimCoarseGridModel::buildGrid(), ossimGeoPdfReader::computeBoundingRect(), computeGroundToImagePartials(), ossimElevRemapper::elevRemap(), ossimHdf5Tool::execute(), ossimSubImageTool::execute(), ossimKakaduJpipImageGeometryFactory::extendGeometry(), ossimImageGeometryFactory::extendGeometry(), getAdjustableParameterInterface(), ossimPdfWriter::getEpsgCode(), ossimERSTileSource::getImageGeometry(), ossimKakaduJpipHandler::getImageGeometry(), ossimGeneralRasterTileSource::getImageGeometry(), ossimUsgsDemTileSource::getImageGeometry(), ossimOpjJp2Reader::getImageGeometry(), ossimJpegTileSource::getImageGeometry(), ossimMG4LidarReader::getImageGeometry(), ossimGeoPdfReader::getImageGeometry(), ossimMrSidReader::getImageGeometry(), ossimH5ImageHandler::getImageGeometry(), ossimLasReader::getImageGeometry(), ossimKakaduJp2Reader::getImageGeometry(), ossimGpkgReader::getImageGeometry(), ossimImageViewProjectionTransform::getImageToViewBounds(), ossimImageViewProjectionTransform::getInputMetersPerPixel(), ossimPdfWriter::getLgiDictCtm(), ossimChipperUtil::getMapProjection(), ossimImageViewProjectionTransform::getOutputMetersPerPixel(), ATP::AtpGenerator::getValidVertices(), ossimGeoAnnotationEllipseObject::getWidthHeightInPixels(), ossimChipperUtil::hasSensorModelInput(), ossimImageViewProjectionTransform::imageToView(), ossimLasReader::init(), ossimElevRemapper::initialize(), ossimImageProjectionModel::initialize(), ossimViewshedTool::initializeProjectionGsd(), ossimChipProcTool::initializeProjectionGsd(), ossimTiledElevationDatabase::initializeReferenceProjection(), ossimImageViewProjectionTransform::initializeViewSize(), ossimViewshedTool::initProcessingChain(), ossimMapCompositionSource::inputMapProjection(), ossimGdalOgrVectorAnnotation::loadExternalGeometryFile(), ossimGdalOgrVectorAnnotation::loadLineString(), ossimGdalOgrVectorAnnotation::loadMultiLineString(), ossimGdalOgrVectorAnnotation::loadMultiPoint(), ossimGdalOgrVectorAnnotation::loadMultiPolygon(), ossimGdalOgrVectorAnnotation::loadPoint(), ossimGdalOgrVectorAnnotation::loadPolygon(), ossimChipProcTool::newIdentityProjection(), ossimGeoidImage::open(), ossimVpfAnnotationSource::ossimVpfAnnotationSource(), ossimTilingPoly::parseShpFile(), ossimPointCloudImageHandler::setGSD(), ossimRpcSolver::solveCoefficients(), ossimOrthoImageMosaic::updateGeometry(), ossimRLevelFilter::updateGeometry(), ossimScaleFilter::updateGeometry(), ossimGdalOgrVectorAnnotation::verifyViewParams(), ossimImageViewProjectionTransform::viewToImage(), ossimERSFileWriter::writeFile(), ossimEnviHeaderFileWriter::writeFile(), ossimFgdcFileWriter::writeFile(), ossimNitfWriterBase::writeGeometry(), ossimOpjCompressor::writeGeotiffBox(), ossimKakaduCompressor::writeGeotiffBox(), ossimGdalWriter::writeProjectionInfo(), and ossimWriter::writeTiffTags().

{ return m_projection.get(); }

getProjection() [2/2]

ossimProjection* ossimImageGeometry::getProjection ( )

inline

Definition at line 166 of file ossimImageGeometry.h.

{ return m_projection.get(); }

getTargetRrds()

ossim_uint32 ossimImageGeometry::getTargetRrds

(

)

const

Returns

The target zero based reduced resolution data set used for localToWorld and worldToLocal.

See also

setTargetRrds

Definition at line 679 of file ossimImageGeometry.cpp.

References m_targetRrds.

{
   return m_targetRrds;
}

getTiePoint() [1/2]

void ossimImageGeometry::getTiePoint	(	ossimGpt &	tie,
		bool	edge
	)		const

Get the latitude, longitude of the tie point.

This is the bounding upper left point of the image which is not necessarily the image (0, 0) point.

Parameters

tie	ossimGpt to assign.
edge	If true the tie point is shifted up and to the right by half the gsd.

Note

Requires projection and image size to be initialized or the latitude and longitue will be set to nan.

The height parameter of the ossimGpt is not touched by this method.

Definition at line 732 of file ossimImageGeometry.cpp.

References ossimProjection::forward(), ossimRefPtr< T >::get(), getBoundingGroundRect(), ossimMapProjection::getDecimalDegreesPerPixel(), ossimProjection::getMetersPerPixel(), ossimIpt::hasNans(), ossimDpt::hasNans(), ossimProjection::inverse(), ossimMapProjection::isGeographic(), ossimGpt::lat, localToWorld(), ossimGpt::lon, m_imageSize, m_projection, ossimGpt::makeNan(), ossim::nan(), ossimGrect::ul(), ossimRefPtr< T >::valid(), ossimIpt::x, ossimDpt::x, ossimIpt::y, and ossimDpt::y.

Referenced by ATP::AtpGenerator::constructChain(), getTiePoint(), and ossimEquationUtil::initInputSources().

{
   ossimGrect grect;
   getBoundingGroundRect(grect);
   if ( m_projection.valid() && (m_imageSize.hasNans() == false) )
   {
      // Use the easting/northing version of this method if underlying projection is meters:
      const ossimMapProjection* map_proj = 
         dynamic_cast<const ossimMapProjection*>(m_projection.get());
      if (map_proj)
      {
         if(!map_proj->isGeographic())
         {
            ossimDpt enTie;
            getTiePoint(enTie, edge);
            if (!enTie.hasNans())
               tie = m_projection->inverse(enTie);
            else
               tie.makeNan();
         }
         else
         {
            tie = grect.ul();
            if(edge)
            {
               ossimDpt half_pixel_shift =  map_proj->getDecimalDegreesPerPixel() * 0.5;
               tie.lat += half_pixel_shift.y;
               tie.lon -= half_pixel_shift.x;
            }
         }
      }
      else
      {

         // must be a sensor model so we will set to the upper left bounds of the image
         tie = grect.ul();
//         std::cout << "GRECT ======= " << grect << "\n";
//         std::cout << m_projection->getClassName() << std::endl;
//         ossimDpt pt = m_projection->forward( tie );
//         std::cout << "FORWARD:" <<pt << std::endl;
      }
   }


#if 0
  // std::cout << "GRECT ======= " << grect << "\n";
   if ( m_projection.valid() && (m_imageSize.hasNans() == false) )
   {
      // Use the easting/northing version of this method if underlying projection is meters:
      const ossimMapProjection* map_proj = 
         dynamic_cast<const ossimMapProjection*>(m_projection.get());
      if (map_proj && !map_proj->isGeographic())
      {
         ossimDpt enTie;
         getTiePoint(enTie, edge);
         if (!enTie.hasNans())
            tie = m_projection->inverse(enTie);
         else
            tie.makeNan();
         return; // return here only because it bugs Dave
      }

      // Use projection to ground to establish UL extreme of image:
      ossimDpt iul (0,0);
      ossimDpt iur (m_imageSize.x-1, 0);
      ossimDpt ilr (m_imageSize.x-1, m_imageSize.y-1);
      ossimDpt ill (0, m_imageSize.y-1);
      ossimDpt iRight(1, 0);
      ossimDpt iDown(0, 1);
      
      ossimGpt gul;
      ossimGpt gur; 
      ossimGpt glr;
      ossimGpt gll;
      ossimGpt gRight;
      ossimGpt gDown;
      
      localToWorld(iul, gul);
      localToWorld(iur, gur);
      localToWorld(ilr, glr);
      localToWorld(ill, gll);
      localToWorld(iRight, gRight);
      localToWorld(iDown, gDown);
      
      std::cout << "gul: " << gul << "\n"
                << "gur: " << gur << "\n"
                << "glr: " << glr << "\n"
                << "gll: " << gll << "\n"
                << "gRight: " << gRight << "\n"
                << "gDown:  " << gDown << "\n";
      // Determine the direction of the image:
      if ( gul.lat > gDown.lat ) // oriented north up
      {
         if ( gul.lat >= gRight.lat ) // straight up or leaning right
         {
            std::cout << "ONE\n";
            tie.lat = gul.lat;
            tie.lon = gll.lon;
         }
         else // leaning left
         {
            std::cout << "TWO\n";
            tie.lat = gur.lat;
            tie.lon = gul.lon;
         }
      }
      else // south or on side
      {
         if ( gRight.lat >= gul.lat ) // straight down or leaning right
         {
            std::cout << "THREE\n";
            tie.lat = glr.lat;
            tie.lon = gur.lon;
         }
         else // leaning left
         {
            std::cout << "FOUR\n";
            tie.lat = gll.lat;
            tie.lon = glr.lon;
         }   
      }
      std::cout << "TIE BEFORE ==== " << tie << "\n";

      if ( edge )
      {
         ossimDpt pt = m_projection->forward( tie );
         ossimDpt half_pixel_shift = m_projection->getMetersPerPixel() * 0.5;
         pt.y += half_pixel_shift.y;
         pt.x -= half_pixel_shift.x;
         tie = m_projection->inverse( pt );
      }
      
      std::cout << "TIE ==== " << tie << "\n";
   } // if ( (m_imageSize.hasNans() == false) && m_projection.valid() )
   else
   { 
      tie.lat = ossim::nan();
      tie.lon = ossim::nan();
   }
#endif
}

getTiePoint() [2/2]

void ossimImageGeometry::getTiePoint	(	ossimDpt &	tie,
		bool	edge
	)		const

Get the easting, northing of the tie point.

This is the bounding upper left point of the image which is not necessarily the image (0, 0) point.

Parameters

tie	ossimDpt to assign.
edge	If true the tie point is shifted up and to the right by half the gsd.

Note

Requires projection and image size to be initialized or the easting and northing will be set to nan.

Definition at line 878 of file ossimImageGeometry.cpp.

References ossimProjection::forward(), ossimRefPtr< T >::get(), ossimMapProjection::getMetersPerPixel(), getTiePoint(), ossimIpt::hasNans(), ossimDpt::hasNans(), ossimGpt::hasNans(), ossimGpt::hgt, ossimMapProjection::isGeographic(), ossimMapProjection::lineSampleToEastingNorthing(), m_imageSize, m_projection, ossimDpt::makeNan(), rnToFull(), ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

{
   if (!m_projection.valid() || m_imageSize.hasNans())
   {
      tie.makeNan();
      return; // return here only because it bugs Dave
   }

   // Use the geographic version of this method if underlying projection is NOT map in meters:
   const ossimMapProjection* map_proj = dynamic_cast<const ossimMapProjection*>(m_projection.get());
   if (!map_proj || map_proj->isGeographic())
   {
      // Use the geographic version of this method to establish UL:
      ossimGpt gTie;
      gTie.hgt = 0.0;
      getTiePoint(gTie, edge);
      if (!gTie.hasNans())
         tie = m_projection->forward( gTie );
      else
         tie.makeNan();
      return; // return here only because it bugs Dave
   }

   // The underlying projection is a proper map projection in meters. Use easting northing 
   // directly to avoid shifting the UL tiepoint because of skewed edge in geographic. Note:
   // assume the image is North up:
   ossimDpt iul (0,0);
   ossimDpt ful;
   rnToFull(iul, 0, ful);
   map_proj->lineSampleToEastingNorthing(ful, tie);
   if (edge && !tie.hasNans())
   {
      // Shift from pixel-is-point to pixel-is-area
      ossimDpt half_pixel_shift = map_proj->getMetersPerPixel() * 0.5;
      tie.y += half_pixel_shift.y;
      tie.x -= half_pixel_shift.x;
   }
}

getTransform() [1/2]

const ossim2dTo2dTransform* ossimImageGeometry::getTransform ( ) const

inline

Access methods for transform (may be NULL pointer).

Definition at line 161 of file ossimImageGeometry.h.

Referenced by ossimKakaduJpipImageGeometryFactory::extendGeometry(), ossimImageGeometryFactory::extendGeometry(), ossimQuickbirdNitfTileSource::getImageGeometry(), ossimImageViewProjectionTransform::imageToView(), and ossimImageViewProjectionTransform::viewToImage().

{ return m_transform.get(); }

getTransform() [2/2]

ossim2dTo2dTransform* ossimImageGeometry::getTransform ( )

inline

Definition at line 162 of file ossimImageGeometry.h.

{ return m_transform.get(); }

hasProjection()

bool ossimImageGeometry::hasProjection ( ) const

inline

Returns TRUE if valid projection defined.

Definition at line 190 of file ossimImageGeometry.h.

Referenced by ossimImageGeometryFactory::extendGeometry(), ossimImageViewProjectionTransform::getImageToViewBounds(), ossimImageViewProjectionTransform::getInputMetersPerPixel(), ossimImageViewProjectionTransform::getOutputMetersPerPixel(), ossimGeoAnnotationEllipseObject::getWidthHeightInPixels(), ossimGeoAnnotationSource::loadState(), and ossimNitfWriterBase::writeGeometry().

{ return m_projection.valid(); }

hasTransform()

bool ossimImageGeometry::hasTransform ( ) const

inline

Returns TRUE if valid transform defined.

Definition at line 193 of file ossimImageGeometry.h.

Referenced by ossimKakaduJpipImageGeometryFactory::extendGeometry(), and ossimImageGeometryFactory::extendGeometry().

{ return m_transform.valid(); }

isAffectedByElevation()

bool ossimImageGeometry::isAffectedByElevation

(

)

const

Returns TRUE if this geometry is sensitive to elevation.

Definition at line 332 of file ossimImageGeometry.cpp.

References ossimProjection::isAffectedByElevation(), m_projection, and ossimRefPtr< T >::valid().

Referenced by ossimInfo::getUpIsUpAngle(), ossimImageViewProjectionTransform::getViewSegments(), ossimEquationUtil::initInputSources(), and worldToLocal().

{
    if (m_projection.valid())
        return m_projection->isAffectedByElevation();
    return false;
}

isEqualTo()

bool ossimImageGeometry::isEqualTo ( const ossimObject &  obj, ossimCompareType  compareType = OSSIM_COMPARE_FULL  ) const

virtual

< Maintains full_image-to-world_space transformation

Reimplemented from ossimObject.

Definition at line 978 of file ossimImageGeometry.cpp.

References ossimRefPtr< T >::get(), ossimObject::isEqualTo(), ossimProjection::isEqualTo(), ossimIpt::isEqualTo(), m_decimationFactors, m_imageSize, m_projection, m_targetRrds, m_transform, OSSIM_COMPARE_IMMEDIATE, and ossimRefPtr< T >::valid().

Referenced by ossimGeoPolyCutter::setView().

{
   bool result = ossimObject::isEqualTo(obj, compareType);
   const ossimImageGeometry* rhs = dynamic_cast<const ossimImageGeometry*> (&obj);
   if(rhs&&result) // we know the types are the same
   {
      result = ((m_decimationFactors.size() == rhs->m_decimationFactors.size())&&
                m_imageSize.isEqualTo(rhs->m_imageSize)&& 
                (m_targetRrds == rhs->m_targetRrds)); 
      
      if(result)
      {
         ossim_uint32 decimationIdx = 0;
         for(decimationIdx = 0; result&&(decimationIdx < m_decimationFactors.size());++decimationIdx)
         {
            result = m_decimationFactors[decimationIdx].isEqualTo(rhs->m_decimationFactors[decimationIdx]);
         }
      }
      
      if(result && compareType==OSSIM_COMPARE_IMMEDIATE)
      {
         result = ((m_transform.get()  == rhs->m_transform.get())&& 
                   (m_projection.get() == rhs->m_projection.get()));  
         
      }
      else
      {
         if(m_transform.valid()&&rhs->m_transform.valid())
         {
            result = m_transform->isEqualTo(*rhs->m_transform.get());
         }
         else if(reinterpret_cast<ossim_uint64>(m_transform.get()) | reinterpret_cast<ossim_uint64>(rhs->m_transform.get())) // one is null
         {
            result = false;
         }
         if(m_projection.valid()&&rhs->m_projection.valid())
         {
            result = m_projection->isEqualTo(*rhs->m_projection.get());
         }
         else if(reinterpret_cast<ossim_uint64>(m_projection.get()) | reinterpret_cast<ossim_uint64>(rhs->m_projection.get())) // one is null
         {
            result = false;
         }
         
      }
   }
   return result;
}

isMapProjected()

bool ossimImageGeometry::isMapProjected

(

)

const

Returns true if underlying projection is derived from ossimMapProjection.

Returns

true if map projected; false, if not.

Definition at line 1681 of file ossimImageGeometry.cpp.

References ossimRefPtr< T >::get(), and m_projection.

{
   return dynamic_cast<const ossimMapProjection*>( m_projection.get() ) != 0;
}

loadState()

bool ossimImageGeometry::loadState ( const ossimKeywordlist &  kwl, const char *  prefix = 0  )

virtual

Attempts to initialize a transform and a projection given the KWL.

Returns TRUE if no error encountered.

Reimplemented from ossimObject.

Definition at line 547 of file ossimImageGeometry.cpp.

References ossimString::c_str(), ossimKeywordlist::find(), ossim2dTo2dTransformRegistry::instance(), ossimProjectionFactoryRegistry::instance(), ossimObject::loadState(), m_decimationFactors, m_imageSize, m_projection, m_targetRrds, m_transform, setProjection(), ossimString::size(), STATIC_TYPE_NAME, ossimIpt::toPoint(), ossimString::toUInt32(), ossim::toVector(), ossimKeywordNames::TYPE_KW, and ossimRefPtr< T >::valid().

Referenced by ossimImageGeometryFactory::createGeometry(), ossimImageHandler::getExternalImageGeometry(), ossimDoqqTileSource::getImageGeometry(), ossimERSTileSource::getImageGeometry(), ossimRpfCacheTileSource::getImageGeometry(), ossimImageCacheBase::getImageGeometry(), ossimGdalOgrVectorAnnotation::loadExternalGeometryFile(), ossimGeoAnnotationSource::loadState(), ossimMemoryImageSource::loadState(), ossimGeoPolyCutter::loadState(), ossimImageViewProjectionTransform::loadState(), ossimOgrVectorTileSource::loadState(), ossimEsriShapeFileFilter::loadState(), and ossimGeneralRasterElevHandler::setFilename().

{
   const char* lookup = kwl.find(prefix, ossimKeywordNames::TYPE_KW);
   if (lookup)
   {
      if ( ossimString(lookup) == STATIC_TYPE_NAME(ossimImageGeometry) )
      {
         ossimObject::loadState(kwl, prefix);

         // m_transform
         ossimString transformPrefix = ossimString(prefix) + "transform.";
         ossimRefPtr<ossim2dTo2dTransform> transform = ossim2dTo2dTransformRegistry::instance()->
            createNativeObjectFromRegistry(kwl, transformPrefix.c_str());
         if( transform.valid() )
         {
            m_transform = transform;
         }

         // m_projection:
         ossimString projectionPrefix = ossimString(prefix) + "projection.";
         ossimRefPtr<ossimProjection> projection = ossimProjectionFactoryRegistry::instance()->
            createProjection(kwl, projectionPrefix.c_str());
         if( projection.valid() )
         {
            m_projection = projection;
         }

         // m_decimationFactors:
         ossimString decimations = kwl.find(prefix, "decimations");
         if( decimations.size() )
         {
            m_decimationFactors.clear();
            ossim::toVector(m_decimationFactors, decimations);
         }

         // m_imageSize:
         ossimString imageSize = kwl.find(prefix, "image_size");
         if( imageSize.size() )
         {
            m_imageSize.toPoint(imageSize);
         }

         // m_targetRrds:
         ossimString targetRrds = kwl.find(prefix, "target_rrds");
         if ( targetRrds.size() )
         {
            m_targetRrds = ossimString(targetRrds).toUInt32();
         }
      }
      else
      {
         // Now look for projection spec (for backwards compatibility):
         ossimProjection* projection = ossimProjectionFactoryRegistry::instance()->
            createProjection(kwl, prefix);
         if (projection)
         {
            setProjection(projection);
         }
      }
   }
   else
   {
      //---
      // Old geometry file with no type keyword:
      //---
      ossimProjection* projection = ossimProjectionFactoryRegistry::instance()->
         createProjection(kwl, prefix);
      if (projection)
      {
         setProjection(projection);
      } 
   }
   
   return true;
}

localToWorld() [1/3]

bool ossimImageGeometry::localToWorld	(	const ossimDpt &	local_pt,
		ossimGpt &	world_pt
	)		const

Exposes the 3D projection from image to world coordinates.

The caller should verify that a valid projection exists before calling this method. Returns TRUE if a valid ground point is available in the ground_pt argument. This method depends on the existence of elevation information. If no DEM is available, the results will be incorrect or inaccurate.

Return a NAN point of no projection is defined:

Definition at line 178 of file ossimImageGeometry.cpp.

References ossimProjection::lineSampleToWorld(), m_projection, m_targetRrds, ossimGpt::makeNan(), rnToFull(), and ossimRefPtr< T >::valid().

Referenced by ossimCoarseGridModel::buildGrid(), calculatePolyBounds(), ATP::AtpTileSource::computeParallax(), RadialProcessor::doRadial(), ossimEsriShapeFileFilter::drawAnnotations(), ossimInfo::getCenterGround(), ossimHlzTool::getChip(), ossimShorelineTool::getChip(), ossimViewshedTool::getChip(), getCornerGpts(), getCrossesDateline(), getDegreesPerPixel(), ossimGdalOgrVectorAnnotation::getFeatures(), ossimGmlSupportData::getGeoBounds(), ossimGmlSupportData::getGeoOrigin(), ossimInfo::getImageBounds(), ossimInfo::getImageGeometryInfo(), ossimImageViewProjectionTransform::getImageToViewBounds(), ossimPdfWriter::getLgiDictCtm(), ossimGmlSupportData::getMapOrigin(), getMetersPerPixel(), getTiePoint(), ossimPointCloudUtilityFilter::getTile(), ossimImageViewProjectionTransform::getViewSegments(), ossimIvtGeomXform::imageToGround(), ossimImageViewProjectionTransform::imageToView(), ATP::AtpGenerator::initialize(), ossimImageViewProjectionTransform::initializeViewSize(), ossimViewshedTool::initProcessingChain(), ossimGeoPolyCutter::invertPolygon(), localToWorld(), main(), ossimTiledElevationDatabase::mapRegion(), northUpAngle(), rnToWorld(), ossimGeneralRasterElevHandler::setFilename(), ossimPotraceTool::transformLineStrings(), ossimPotraceTool::transformPolygons(), upIsUpAngle(), and ossimImageViewProjectionTransform::viewToImage().

{
   if (!m_projection.valid())
   {
      world_pt.makeNan();
      return false;
   }

   // First transform local pixel to full-image pixel:
   ossimDpt full_image_pt;
   rnToFull(local_pt, m_targetRrds, full_image_pt);

   // Perform projection to world coordinates:
   m_projection->lineSampleToWorld(full_image_pt, world_pt);


    // Put longitude between -180 and +180 and latitude between -90 and +90 if not so. 
   // world_pt.wrap();
    
   return true;
}

localToWorld() [2/3]

bool ossimImageGeometry::localToWorld	(	const ossimDrect &	local_rect,
		ossimGrect &	world_rect
	)		const

Definition at line 201 of file ossimImageGeometry.cpp.

References ossimDrect::ll(), localToWorld(), ossimDrect::lr(), ossimDrect::ul(), and ossimDrect::ur().

{
   ossimGpt gp1, gp2, gp3, gp4;
   if (  localToWorld(local_rect.ul(), gp1) && localToWorld(local_rect.ur(), gp2) &&
         localToWorld(local_rect.lr(), gp3) && localToWorld(local_rect.ll(), gp4))
   {
      world_rect = ossimGrect(gp1, gp2, gp3, gp4);
      return true;
   }
   return false;
}

localToWorld() [3/3]

bool ossimImageGeometry::localToWorld	(	const ossimDpt &	local_pt,
		const double &	h_ellipsoid,
		ossimGpt &	world_pt
	)		const

Exposes the 3D projection from image to world coordinates given a constant height above ellipsoid.

The caller should verify that a valid projection exists before calling this method. Returns TRUE if a valid ground point is available in the ground_pt argument.

Return a NAN point of no projection is defined:

Definition at line 218 of file ossimImageGeometry.cpp.

References ossimProjection::lineSampleHeightToWorld(), m_projection, m_targetRrds, ossimGpt::makeNan(), rnToFull(), ossimRefPtr< T >::valid(), and ossimGpt::wrap().

{
   if (!m_projection.valid())
   {
      world_pt.makeNan();
      return false;
   }

   // First transform local pixel to full-image pixel:
   ossimDpt full_image_pt;
   rnToFull(local_pt, m_targetRrds, full_image_pt);

   // Perform projection to world coordinates:
   m_projection->lineSampleHeightToWorld(full_image_pt, h_ellipsoid, world_pt);

   // Put longitude between -180 and +180 and latitude between -90 and +90 if not so. 
   world_pt.wrap();
    
   return true;
}

northUpAngle()

ossim_float64 ossimImageGeometry::northUpAngle

(

)

const

Returns

Returns the angle to "north up" in decimal degrees, ossim::nan on error.

Definition at line 1645 of file ossimImageGeometry.cpp.

References getBoundingRect(), ossimGpt::hasNans(), ossimDrect::hasNans(), localToWorld(), ossimDrect::lr(), m_projection, ossimDrect::midPoint(), ossim::nan(), ossimDrect::ul(), ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

Referenced by ossimInfo::getNorthUpAngle().

{
   ossim_float64 result = ossim::nan();

   if ( m_projection.valid() )
   {
      ossimDrect bounds;
      getBoundingRect( bounds );
      
      if( !bounds.hasNans() )
      {
         ossimDpt centerIpt = bounds.midPoint();

         ossimDpt midBottomIpt( centerIpt.x, bounds.lr().y );
         ossimDpt midTopIpt( centerIpt.x, bounds.ul().y );

         ossimGpt midBottomGpt;
         ossimGpt midTopGpt;
         
         localToWorld( midBottomIpt, midBottomGpt );
         localToWorld( midTopIpt, midTopGpt );

         if ( !midBottomGpt.hasNans() && !midTopGpt.hasNans() )
         {
            result = midBottomGpt.azimuthTo( midTopGpt );

            if(result < 0) result += 360.0;
         }
      }
      
   } // Matches: if ( m_projection.valid() )

   return result;
   
} // End: ossimImageGeometry::northUpAngle()

open()

bool ossimImageGeometry::open

(

const ossimFilename &

image

)

Shortcut way of getting to an image's geometry when access to pixels is not needed.

Parameters

image

Image to open

Returns

true if open was successful and a valid geometry was initialized.

Definition at line 84 of file ossimImageGeometry.cpp.

References ossimImageHandler::getImageGeometry(), ossimImageHandlerRegistry::instance(), ossimImageHandlerRegistry::open(), and ossimRefPtr< T >::valid().

Referenced by ossimChipProcTool::newIdentityProjection(), and ossimChipProcTool::setAoiToInputs().

{
   ossimRefPtr<ossimImageHandler> handler = ossimImageHandlerRegistry::instance()->open(image);
   if (!handler.valid())
      return false;

   ossimRefPtr<ossimImageGeometry> geom = handler->getImageGeometry();
   if (!geom.valid())
      return false;

   *this = *geom;
   return true;
}

operator=()

const ossimImageGeometry & ossimImageGeometry::operator=

(

const ossimImageGeometry &

copy_this

)

Replaces any existing transform and projection with those in the copy_this object.

Definition at line 687 of file ossimImageGeometry.cpp.

References ossimObject::dup(), ossimProjection::dup(), m_decimationFactors, m_imageSize, m_projection, m_targetRrds, m_transform, and ossimRefPtr< T >::valid().

{
   if (this != &copy_this)
   {
      // Get deep copy of 2D transform if one exists:
      if (copy_this.m_transform.valid())
      {
         m_transform = (ossim2dTo2dTransform*) copy_this.m_transform->dup();
      }
      
      // Now establish a deep copy of the projection, if any:
      if (copy_this.m_projection.valid())
      {
         m_projection = (ossimProjection*) copy_this.m_projection->dup();
      }
      
      // the Gsd should already be solved from the source we are copying from
      m_imageSize         = copy_this.m_imageSize;
      m_decimationFactors = copy_this.m_decimationFactors;
      m_targetRrds        = copy_this.m_targetRrds;
   }
   return *this;
}

operator==()

bool ossimImageGeometry::operator==

(

const ossimImageGeometry &

other

)

const

Returns TRUE if argument geometry has exactly the same instances for transform and m_projection.

This should be expanded to permit different instances of same transforms

This may need to be expanded to permit different instances of same transforms

Definition at line 477 of file ossimImageGeometry.cpp.

References decimationFactor(), m_projection, and m_transform.

{
    return ((m_transform == other.m_transform) && (m_projection == other.m_projection) &&
            (decimationFactor(0) == other.decimationFactor(0)) );
}

print()

std::ostream & ossimImageGeometry::print ( std::ostream &  out ) const

virtual

Prints contents to output stream.

Reimplemented from ossimObject.

Definition at line 439 of file ossimImageGeometry.cpp.

References m_decimationFactors, m_imageSize, m_projection, m_targetRrds, m_transform, ossim2dTo2dTransform::print(), ossimProjection::print(), and ossimRefPtr< T >::valid().

Referenced by ossimImageViewProjectionTransform::print().

{
   out << "type: ossimImageGeometry" << std::endl;
   if(m_transform.valid())
   {
      out << "  m_transform: ";
      m_transform->print(out);
   }
   else
   {
      out << "  No transform defined. Using identity transform.\n";
   }
   
   if(m_projection.valid())
   {
      out << "  m_projection: ";
      m_projection->print(out);
   }
   else
   {
      out << "  No projection defined. ";
   }

   for ( std::vector<ossimDpt>::size_type i = 0; i < m_decimationFactors.size(); ++i )
   {
      cout << "m_decimationFactors[" << i << "]: " << m_decimationFactors[i] << "\n";
   }

   out << "m_imageSize: " << m_imageSize
       << "\nm_targetRrds: " << m_targetRrds << "\n";

   return out;
}

rnToFull()

void ossimImageGeometry::rnToFull	(	const ossimDpt &	rnPt,
		ossim_uint32	resolutionLevel,
		ossimDpt &	fullPt
	)		const

rnToFull is a utility method that takes a rn resolution image point and maps it to the full image point.

Parameters

rnPt	Is a point in resolution n.
resolutionLevel	Is the resolution of the point rnPt. a value of 0 is the local image
fullPt	Is the result of the transform

Definition at line 116 of file ossimImageGeometry.cpp.

References ossim2dTo2dTransform::forward(), ossimDpt::hasNans(), m_transform, undecimatePoint(), and ossimRefPtr< T >::valid().

Referenced by getTiePoint(), ossimImageViewProjectionTransform::imageToView(), localToWorld(), and ossimImageViewProjectionTransform::viewToImage().

{
   // Back out the decimation.
   ossimDpt localPt;
   undecimatePoint(rnPt, resolutionLevel, localPt);

   // Remove any shift/rotation.
   if ( m_transform.valid() && !localPt.hasNans() )
   {
      m_transform->forward(localPt, fullPt);
   }
   else
   {
      fullPt = localPt; // No transform (shift/rotation)
   }
}

rnToRn()

void ossimImageGeometry::rnToRn	(	const ossimDpt &	inRnPt,
		ossim_uint32	inResolutionLevel,
		ossim_uint32	outResolutionLevel,
		ossimDpt &	outRnPt
	)		const

rnToRn is a utility method that takes a rn resolution image point and maps it to the another rn resolution image point.

Parameters

inRnPt	Is a point in resolution n.
inResolutionLevel	Is the resolution of the point inRnPt.
outResolutionLevel	Is the resolution of the point outRnPt.
outRnPt	Is the result of the transform.

Definition at line 98 of file ossimImageGeometry.cpp.

References decimatePoint(), and undecimatePoint().

Referenced by rnToWorld(), and worldToRn().

{
   if (inResolutionLevel != outResolutionLevel)
   {
      // Back out the decimation.
      ossimDpt pt;
      undecimatePoint(inRnPt, inResolutionLevel, pt);

      // Decimate to new level.
      decimatePoint(pt, outResolutionLevel, outRnPt);
   }
   else
   {
      outRnPt = inRnPt; // No transform required.
   }
}

rnToWorld()

void ossimImageGeometry::rnToWorld	(	const ossimDpt &	rnPt,
		ossim_uint32	resolutionLevel,
		ossimGpt &	wpt
	)		const

rnToWorld is a utility method that takes a rn resolution image point and maps it to the world point.

Parameters

rnPt	Is a point in resolution n.
resolutionLevel	Is the resolution of the point rnPt. a value of 0 is the local image
wpt	Is the result of the transform

Definition at line 154 of file ossimImageGeometry.cpp.

References localToWorld(), m_targetRrds, and rnToRn().

Referenced by ossimElevRemapper::elevRemap(), and ossimPointCloudImageHandler::getTile().

{
   ossimDpt localPt;
   rnToRn(rnPt, resolutionLevel, m_targetRrds, localPt);
   localToWorld(localPt, wpt);
}

saveState()

bool ossimImageGeometry::saveState ( ossimKeywordlist &  kwl, const char *  prefix = 0  ) const

virtual

Saves the transform (if any) and projection (if any) states to the KWL.

Reimplemented from ossimObject.

Definition at line 627 of file ossimImageGeometry.cpp.

References ossimKeywordlist::add(), ossimString::c_str(), getMetersPerPixel(), m_decimationFactors, m_imageSize, m_projection, m_targetRrds, m_transform, ossim2dTo2dTransform::saveState(), ossimObject::saveState(), ossimProjection::saveState(), ossimIpt::toString(), ossimDpt::toString(), ossim::toStringList(), and ossimRefPtr< T >::valid().

Referenced by ossimHdf5Tool::execute(), ossimSubImageTool::execute(), ossimInfo::getImageGeometryInfo(), ossimImageHandler::saveImageGeometry(), ossimElevRemapper::saveState(), ossimGdalOgrVectorAnnotation::saveState(), ossimMemoryImageSource::saveState(), ossimImageViewProjectionTransform::saveState(), ossimOgrVectorTileSource::saveState(), ossimEsriShapeFileFilter::saveState(), and ossimGeoPdfReader::saveState().

{
   bool good_save = true;

   // m_transform:
   if (m_transform.valid())
   {
      ossimString transformPrefix = ossimString(prefix) + "transform.";
      good_save = m_transform->saveState(kwl, transformPrefix.c_str());
   }

   // m_projection:
   if (m_projection.valid())
   {
      ossimString projectionPrefix = ossimString(prefix) + "projection.";
      good_save &= m_projection->saveState(kwl, projectionPrefix.c_str());
   }

   // m_gsd:
   ossimDpt gsd;
   getMetersPerPixel(gsd);
   kwl.add(prefix, "gsd", gsd.toString(), true);

   // m_decimationFactors:
   if(m_decimationFactors.size())
   {
      ossimString resultPoints;
      ossim::toStringList(resultPoints, m_decimationFactors);
      kwl.add(prefix, "decimations", resultPoints, true);
   }
   else
   {
      kwl.add(prefix, "decimations", "", true);
   }

   // m_imageSize:
   kwl.add(prefix, "image_size", m_imageSize.toString(), true);

   // m_targetRrds;
   kwl.add(prefix, "target_rrds", m_targetRrds, true);

   // Base class:
   good_save &= ossimObject::saveState(kwl, prefix);
   
   return good_save;
}

setDiscreteDecimation()

void ossimImageGeometry::setDiscreteDecimation ( const std::vector< ossimDpt > &  decimation_list )

inline

Sets the decimation scheme to a discrete list of decimation factors.

Definition at line 318 of file ossimImageGeometry.h.

Referenced by ossimMrSidReader::getDecimationFactors(), and ossimImageHandler::initImageParameters().

      { m_decimationFactors = decimation_list; }

setImageSize()

void ossimImageGeometry::setImageSize ( const ossimIpt &  size )

inline

Definition at line 326 of file ossimImageGeometry.h.

References size.

Referenced by ossimSubImageTool::execute(), ossimHlzTool::getChip(), ossimShorelineTool::getChip(), ossimViewshedTool::getChip(), ossimChipProcTool::getChip(), ossimChipperUtil::getChip(), ossimLasReader::getImageGeometry(), ATP::AtpGenerator::initialize(), ossimRangeDomeTileSource::initialize(), ossimTilingRect::initializeBase(), ossimIgen::initializeChain(), ossimImageViewProjectionTransform::initializeViewSize(), ossimImageHandler::initImageParameters(), ossimViewshedTool::initProcessingChain(), and ossimPointCloudImageHandler::setGSD().

   {
      m_imageSize = size;
   }

setProjection()

void ossimImageGeometry::setProjection

(

ossimProjection *

projection

)

Sets the projection to be used for local-to-world coordinate transformation.

Definition at line 324 of file ossimImageGeometry.cpp.

References m_projection.

Referenced by ossimCoarseGridModel::buildGrid(), ossimKakaduJpipImageGeometryFactory::createGeometry(), ossimImageGeometryFactory::extendGeometry(), ossimMrSidReader::getExternalImageGeometry(), ossimGdalTileSource::getExternalImageGeometryFromXml(), ossimQuickbirdTiffTileSource::getImageGeometry(), ossimQbTileFilesHandler::getImageGeometry(), ossimRpfTocEntry::getImageGeometry(), ossimDtedTileSource::getImageGeometry(), ossimGeneralRasterTileSource::getImageGeometry(), ossimPointCloudImageHandler::getImageGeometry(), ossimLasReader::getImageGeometry(), ossimplugins::ossimRadarSat2TiffReader::getInternalImageGeometry(), ossimplugins::ossimTerraSarTiffReader::getInternalImageGeometry(), ossimJpegTileSource::getInternalImageGeometry(), ossimH5ImageHandler::getInternalImageGeometry(), ossimEquationUtil::initInputSources(), ossimGdalOgrVectorAnnotation::loadExternalImageGeometryFromXml(), loadState(), ossimGeoAnnotationSource::setView(), ossimGdalOgrVectorAnnotation::setView(), ossimImageViewProjectionTransform::setView(), ossimEsriShapeFileFilter::setView(), and ossimRpcSolver::solveCoefficients().

{ 
   m_projection = projection; 
}

setTargetRrds()

void ossimImageGeometry::setTargetRrds

(

ossim_uint32

rrds

)

Set m_targetRrds data member.

This is used by methods worldToRn and localToWorld that do not take a rrds argument and methods rnToWorld and worldToRn. If the target rrds is set to 2, then the resulting point from worldToLocal is relative to reduced resolution level 2. Vice versa for localToWorld it is assumed the local point is relative to the target reduced resolution data set.

Parameters

rrds	Target (zero based) reduced resolution data set.

Definition at line 674 of file ossimImageGeometry.cpp.

References m_targetRrds.

{
   m_targetRrds = rrds;
}

setTransform()

void ossimImageGeometry::setTransform

(

ossim2dTo2dTransform *

transform

)

Sets the transform to be used for local-to-full-image coordinate transformation.

Definition at line 316 of file ossimImageGeometry.cpp.

References m_transform.

Referenced by ossimKakaduJpipImageGeometryFactory::createGeometry(), ossimKakaduJpipImageGeometryFactory::extendGeometry(), ossimImageGeometryFactory::extendGeometry(), ossimQuickbirdNitfTileSource::getImageGeometry(), and ossimQbTileFilesHandler::getImageGeometry().

{ 
   m_transform = transform; 
}

undecimatePoint()

void ossimImageGeometry::undecimatePoint ( const ossimDpt &  rnPt, ossim_uint32  resolutionLevel, ossimDpt &  outPt  ) const

protected

Method to back out decimation of a point.

Parameters

rnPt	Is a point in resolution n.
resolutionLevel	Is the resolution of the point rnPt.
outPt	Is the result of the transform a non-decimated point.

Definition at line 918 of file ossimImageGeometry.cpp.

References decimationFactor(), ossimDpt::hasNans(), ossimDpt::makeNan(), ossimDpt::x, and ossimDpt::y.

Referenced by rnToFull(), and rnToRn().

{
   // Back out the decimation.
   ossimDpt decimation = decimationFactor(resolutionLevel);
   
   //---
   // If no nans and one or both of the factors is not 1.0 decimation should
   // be applied.
   //---
   if ( (decimation.x != 1.0) || (decimation.y != 1.0) ) 
   {
      if ( !decimation.hasNans() && !rnPt.hasNans() )
      {
         outPt.x = rnPt.x * (1.0/decimation.x);
         outPt.y = rnPt.y * (1.0/decimation.y);
      }
      else
      {
         outPt.makeNan();
      }
   }
   else
   {
      outPt = rnPt; // No decimation.
   }
}

upIsUpAngle() [1/2]

ossim_float64 ossimImageGeometry::upIsUpAngle

(

const ossimDpt &

pt

)

const

Returns

Returns the angle to "up is up" in decimal degrees, 0.0 if image is not affected by elevation, ossim::nan on error.

Definition at line 1479 of file ossimImageGeometry.cpp.

References ossim::almostEqual(), ossim::atan2d(), FLT_EPSILON, getBoundingRect(), ossimDpt::hasNans(), ossimDrect::hasNans(), ossimProjection::isAffectedByElevation(), ossim::isnan(), ossimDpt::length(), localToWorld(), m_projection, ossimDrect::midPoint(), ossim::nan(), ossimRefPtr< T >::valid(), worldToLocal(), ossimDpt::x, and ossimDpt::y.

{
   ossim_float64 result = ossim::nan();

   if ( m_projection.valid() )
   {
      if ( m_projection->isAffectedByElevation() )
      {
         const int NUMBER_OF_SAMPLES = 9;
         
         // In meters.  This is about a height of a 6 to 7 story building.
         const double ELEVATION_DISPLACEMENT = 20;
         
         ossimDrect bounds;
         getBoundingRect( bounds );
         
         if( !bounds.hasNans() )
         {
            // only need an average
            // so will keep the test close to the pixel location
            // we are sampling
            // let's get rid of the scale
            ossim_float64 widthPercent  = 1;//bounds.width()*.1;
            ossim_float64 heightPercent = 1;//bounds.height()*.1;
            
            //---
            // Sanity check to make sure that taking 10 percent out on the image
            // gets us to at least 1 pixel away.
            //---
           // if(widthPercent < 1.0) widthPercent = 1.0;
           // if(heightPercent < 1.0) heightPercent = 1.0;
            
            // set up some work variables to help calculate the average partial
            //
            std::vector<ossimDpt> ipts(NUMBER_OF_SAMPLES);
            std::vector<ossimGpt> gpts(NUMBER_OF_SAMPLES);
            std::vector<ossimDpt> iptsDisplacement(NUMBER_OF_SAMPLES);
            std::vector<ossimDpt> partials(NUMBER_OF_SAMPLES);
            ossimDpt averageDelta(0.0,0.0);
            ossimDpt centerIpt = pt;
            if(centerIpt.hasNans())
            {
               centerIpt = bounds.midPoint();
            }
            //---
            // Lets take an average displacement about the center point (3x3 grid)
            // we will go 10 percent out of the width and height of the image and
            // look at the displacement at those points and average them we will
            // use the average displacement to compute the up is up rotation.
            //---
            
            // top row:

            // 45 degree left quadrant
            // RP - FIX POINTS IN FIRST ROW
            ipts[0] = centerIpt + ossimDpt(-widthPercent, -heightPercent);
            // 45 degree middle top 
            ipts[1] = centerIpt + ossimDpt(0.0,heightPercent);
            // 45 degree right quadrant
            ipts[2] = centerIpt + ossimDpt(widthPercent, -heightPercent);
            
            // middle row
            ipts[3] = centerIpt + ossimDpt(-widthPercent, 0.0); // left middle
            ipts[4] = centerIpt;
            ipts[5] = centerIpt + ossimDpt(widthPercent, 0.0);
            
            // bottom row
            ipts[6] = centerIpt + ossimDpt(-widthPercent, heightPercent);
            ipts[7] = centerIpt + ossimDpt(0.0, heightPercent);
            ipts[8] = centerIpt + ossimDpt(widthPercent, heightPercent);
            
            ossim_uint32 idx = 0;
            for(idx = 0; idx < ipts.size(); ++idx)
            {
               double h = 0.0;
               localToWorld(ipts[idx], gpts[idx]);
               h = gpts[idx].height();
               if(ossim::isnan(h)) h = 0.0;
               gpts[idx].height(h + ELEVATION_DISPLACEMENT);
               worldToLocal(gpts[idx], iptsDisplacement[idx]);
               averageDelta = averageDelta + (iptsDisplacement[idx] - ipts[idx]);
            }
            
            ossim_float64 averageLength = averageDelta.length();
            if(averageLength < 1) return false;
            
            if(!ossim::almostEqual((double)0.0, (double)averageLength))
            {
               averageDelta = averageDelta/averageLength;
            }
            // RP - the rest of this calculation seems to be leading to less accurate results.  
            // We already have the average image space displacement from raising the height for 9 points in the image 
            // Just rotate based on this displacement, 
            // scaling based on the image size distorts the direction of the up vector, leading to inaccurate results for non-square images
            //
            //ossimDpt averageLocation =  (ossimDpt(averageDelta.x*bounds.width(), 
            //                                      averageDelta.y*bounds.height())+centerIpt) ;
            //ossimGpt averageLocationGpt;
            //ossimGpt centerGpt;
            
            //localToWorld(averageLocation, averageLocationGpt);
            //localToWorld(ipts[0], centerGpt);
            
            ossimDpt deltaPt = averageDelta; //averageLocation - centerIpt; 
            ossimDpt deltaUnitPt = deltaPt;
            double len = deltaPt.length();
            if(len > FLT_EPSILON)
            {
               deltaUnitPt  = deltaUnitPt/len;
            }
            
            // Image space model positive y is down.  Let's reflect to positve y up.
            deltaUnitPt.y *= -1.0; // Reflect y to be right handed.
            
            result = ossim::atan2d(deltaUnitPt.x, deltaUnitPt.y);
            
            //---
            // We are essentially simulating camera rotation so negate the rotation
            // value.
            //---
            result *= -1;
            
            if(result < 0) result += 360.0;
            
         }  // Matches: if( bounds.hasNans() == false )
      }
      else
      {
         result = 0;
      }
      
   } // Matches: if ( m_projection.valid() && m_projection->isAffectedByElevation() )

   return result;

}

upIsUpAngle() [2/2]

ossim_float64 ossimImageGeometry::upIsUpAngle

(

)

const

Definition at line 1617 of file ossimImageGeometry.cpp.

References getBoundingRect(), ossimDrect::hasNans(), ossimProjection::isAffectedByElevation(), m_projection, ossimDrect::midPoint(), ossim::nan(), and ossimRefPtr< T >::valid().

{
   ossim_float64 result = ossim::nan();

   if ( m_projection.valid() )
   {
      if ( m_projection->isAffectedByElevation() )
      {
         
         ossimDrect bounds;
         getBoundingRect( bounds );
         
         if( !bounds.hasNans() )
         {
            return upIsUpAngle(bounds.midPoint());
         }
      }
      else
      {
         result = 0;
      }
      
   } // Matches: if ( m_projection.valid() && m_projection->isAffectedByElevation() )

   return result;
   
} // End: ossimImageGeometry::upIsUpAngle()

worldToLocal() [1/2]

bool ossimImageGeometry::worldToLocal	(	const ossimGpt &	world_pt,
		ossimDpt &	local_pt
	)		const

Exposes the 3D world-to-local image coordinate reverse projection.

The caller should verify that a valid projection exists before calling this method. Returns TRUE if a valid image point is available in the local_pt argument.

Definition at line 247 of file ossimImageGeometry.cpp.

References fullToRn(), ossimElevManager::getHeightAboveEllipsoid(), ossimGpt::hgt, ossimElevManager::instance(), isAffectedByElevation(), ossimGpt::isHgtNan(), m_projection, m_targetRrds, ossimDpt::makeNan(), ossimRefPtr< T >::valid(), and ossimProjection::worldToLineSample().

Referenced by ossimChipperUtil::addCrossHairAnnotation(), ossimChipProcTool::computeAdjustedViewFromGrect(), ossimGeoPdfReader::computeBoundingRect(), computeGroundToImagePartials(), computeImageToGroundPartialsWRTAdjParam(), ATP::AtpTileSource::computeParallax(), ossimChipProcTool::getChip(), ossimImageElevationHandler::getHeightAboveMSL(), ossimGeneralRasterElevHandler::getHeightAboveMSLFileTemplate(), ossimGeneralRasterElevHandler::getHeightAboveMSLMemoryTemplate(), ossimImageViewProjectionTransform::getImageToViewBounds(), ossimChipperUtil::getIrect(), ossimPointObservation::getResiduals(), ossimPointCloudImageHandler::getValidImageVertices(), ossimImageViewProjectionTransform::getViewSegments(), ossimIvtGeomXform::groundToImage(), ossimImageViewProjectionTransform::imageToView(), ossimRangeDomeTileSource::initialize(), ossimGdalOgrVectorAnnotation::initializeBoundingRec(), ossimImageViewProjectionTransform::initializeViewSize(), ossimViewshedTool::initProcessingChain(), ossimTiledElevationDatabase::mapRegion(), ossimGeoidImage::offsetFromEllipsoidTemplate(), ossimKmlSuperOverlayReader::open(), ossimViewshedTool::paintReticle(), ossimImageElevationHandler::pointHasCoverage(), ossimPointCloudImageHandler::setGSD(), ossimGeoAnnotationPolyLineObject::transform(), ossimGeoAnnotationLineObject::transform(), ossimGeoAnnotationMultiPolyObject::transform(), ossimGeoAnnotationMultiPolyLineObject::transform(), ossimGeoAnnotationMultiEllipseObject::transform(), ossimGeoAnnotationBitmap::transform(), ossimGeoAnnotationGdBitmapFont::transform(), ossimGeoAnnotationEllipseObject::transform(), ossimGeoAnnotationPolyObject::transform(), ossimGeoAnnotationFontObject::transform(), ossimGeoPolyCutter::transformVertices(), upIsUpAngle(), ossimImageViewProjectionTransform::viewToImage(), worldToLocal(), and worldToRn().

{
   bool result = true;
   
   if ( m_projection.valid() )
   {
      //const ossimEquDistCylProjection* eqProj =
      //   dynamic_cast<const ossimEquDistCylProjection*>( m_projection.get() );
      
      ossimDpt full_image_pt;

      //***** GCP
      // I am having major problems with the call and is messing up on Image that are edge to edge -180 to 180.
      // It appears to wrap and think that the image is onl 1 pixel wide.  I am commenting out for now
      // until a better solution can be done for points that wrap.  We need a general implementation that will work
      // with any projector
      //

      //if ( eqProj && (m_imageSize.hasNans() == false) )
     // {
         // Call specialized method to handle wrapping...
     //    eqProj->worldToLineSample( world_pt, m_imageSize, full_image_pt );
     // }
     // else if( isAffectedByElevation() )
      if( isAffectedByElevation() )
      {
         ossimGpt copyPt( world_pt );
         if(world_pt.isHgtNan())
         {
            copyPt.hgt = ossimElevManager::instance()->getHeightAboveEllipsoid(copyPt);
         }     

         // Perform projection from world coordinates to full-image space:
         m_projection->worldToLineSample(copyPt, full_image_pt);
      }
      else
      {
         // Perform projection from world coordinates to full-image space:
         m_projection->worldToLineSample(world_pt, full_image_pt);
      }
      
      // Transform to local space:
      fullToRn(full_image_pt, m_targetRrds, local_pt);
   }
   else // No projection set:
   {
      local_pt.makeNan();
      result = false;  
   }
   
   return result;
   
} // End: ossimImageGeometry::worldToLocal(const ossimGpt&, ossimDpt&)

worldToLocal() [2/2]

bool ossimImageGeometry::worldToLocal	(	const ossimGrect &	world_rect,
		ossimDrect &	local_rect
	)		const

Definition at line 301 of file ossimImageGeometry.cpp.

References ossimGrect::ll(), ossimGrect::lr(), ossimGrect::ul(), ossimGrect::ur(), and worldToLocal().

{
   ossimDpt dp1, dp2, dp3, dp4;
   if (  worldToLocal(world_rect.ul(), dp1) && worldToLocal(world_rect.ur(), dp2) &&
         worldToLocal(world_rect.lr(), dp3) && worldToLocal(world_rect.ll(), dp4))
   {
      local_rect = ossimDrect(dp1, dp2, dp3, dp4);
      return true;
   }
   return false;
}

worldToRn()

void ossimImageGeometry::worldToRn	(	const ossimGpt &	wpt,
		ossim_uint32	resolutionLevel,
		ossimDpt &	rnPt
	)		const

worldToRn is a utility method that takes a world point allows one to transform all the way back to an rn point.

Parameters

wpt	Ground point.
resolutionLevel	Is the resolution of the point rnPt. a value of 0 is the local image
rnPt	Is the resoltion point.

Definition at line 163 of file ossimImageGeometry.cpp.

References m_targetRrds, rnToRn(), and worldToLocal().

Referenced by ossimPointCloudImageHandler::getTile().

{   
   ossimDpt localPt;
   worldToLocal(wpt, localPt);
   rnToRn(localPt, m_targetRrds, resolutionLevel, rnPt);
}

Member Data Documentation

m_decimationFactors

std::vector<ossimDpt> ossimImageGeometry::m_decimationFactors

protected

List of decimation factors for R-levels.

Definition at line 481 of file ossimImageGeometry.h.

Referenced by decimationFactor(), decimationFactors(), isEqualTo(), loadState(), operator=(), print(), and saveState().

m_imageSize

ossimIpt ossimImageGeometry::m_imageSize

protected

Definition at line 482 of file ossimImageGeometry.h.

Referenced by applyScale(), getBoundingRect(), getCornerGpts(), getDegreesPerPixel(), getMetersPerPixel(), getTiePoint(), isEqualTo(), loadState(), operator=(), ossimImageGeometry(), print(), and saveState().

m_projection

ossimRefPtr<ossimProjection> ossimImageGeometry::m_projection

protected

Maintains full_image-to-world_space transformation.

Definition at line 480 of file ossimImageGeometry.h.

Referenced by applyScale(), getDegreesPerPixel(), getMetersPerPixel(), getTiePoint(), isAffectedByElevation(), isEqualTo(), isMapProjected(), loadState(), localToWorld(), northUpAngle(), operator=(), operator==(), print(), saveState(), setProjection(), upIsUpAngle(), and worldToLocal().

m_targetRrds

ossim_uint32 ossimImageGeometry::m_targetRrds

protected

Target rrds for localToWorld and worldToLocal methods.

Definition at line 485 of file ossimImageGeometry.h.

Referenced by getTargetRrds(), isEqualTo(), loadState(), localToWorld(), operator=(), print(), rnToWorld(), saveState(), setTargetRrds(), worldToLocal(), and worldToRn().

m_transform

ossimRefPtr<ossim2dTo2dTransform> ossimImageGeometry::m_transform

protected

Maintains local_image-to-full_image transformation.

Definition at line 479 of file ossimImageGeometry.h.

Referenced by fullToRn(), getDegreesPerPixel(), getMetersPerPixel(), isEqualTo(), loadState(), operator=(), operator==(), print(), rnToFull(), saveState(), and setTransform().

---

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

• ossimImageGeometry.h
• ossimImageGeometry.cpp