#include <ossimApplanixUtmModel.h>

Inheritance diagram for ossimApplanixUtmModel:

Public Member Functions
	ossimApplanixUtmModel ()

	ossimApplanixUtmModel (const ossimApplanixUtmModel &src)

virtual ossimObject *	dup () const

virtual void	imagingRay (const ossimDpt &image_point, ossimEcefRay &image_ray) const

void	lineSampleToWorld (const ossimDpt &image_point, ossimGpt &gpt) const

virtual void	lineSampleHeightToWorld (const ossimDpt &image_point, const double &heightEllipsoid, ossimGpt &worldPoint) const

virtual void	worldToLineSample (const ossimGpt &world_point, ossimDpt &image_point) const

virtual void	updateModel ()

void	setPrincipalPoint (ossimDpt principalPoint)

void	setPixelSize (const ossimDpt &pixelSize)

void	setImageRect (const ossimDrect &rect)

void	setFocalLength (double focalLength)

void	setPlatformPosition (const ossimGpt &gpt)

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

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

virtual void	initAdjustableParameters ()

virtual bool	useForward () const

virtual bool	setupOptimizer (const ossimString &init_file)
	ground to image faster (you don't need DEM) //TBC More...

Public Member Functions inherited from ossimSensorModel
	ossimSensorModel ()

	ossimSensorModel (const ossimSensorModel &copy_this)

	ossimSensorModel (const ossimKeywordlist &geom_kwl)

const ossimSensorModel &	operator= (const ossimSensorModel &rhs)
	assignment operator More...

virtual ossimObject *	getBaseObject ()

virtual const ossimObject *	getBaseObject () const

virtual ossimGpt	origin () const

virtual ossimDpt	getMetersPerPixel () const

virtual const double &	getNominalPosError () const
	Returns the estimated Absolute horizontal position error (CE90) of the sensor model. More...

virtual const double &	getRelativePosError () const
	Returns the estimated RELATIVE horizontal position error (CE90) of the sensor model. More...

virtual void	setNominalPosError (const double &ce90)
	Assigns the absolute image position error uncertainty (abs CE90) More...

virtual void	setRelativePosError (const double &ce90)
	Assigns the relative image position error uncertainty (rel CE90) More...

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

void	setRefImgPt (const ossimDpt &pt)
	Sets the center line sampe of the image. More...

void	setRefGndPt (const ossimGpt &pt)
	Sets the center latitude, longitude, height of the image. More...

void	setImageRect (const ossimDrect &imageRect)

void	setGroundRect (const ossimGpt &ul, const ossimGpt &ur, const ossimGpt &lr, const ossimGpt &ll)

ossimDpt	imageSize () const

void	setImageSize (const ossimDpt &size)

virtual void	adjustableParametersChanged ()

virtual bool	insideImage (const ossimDpt &p) const

virtual bool	operator== (const ossimProjection &proj) const

const ossimString &	getImageID () const
	Access methods: More...

const ossimDrect &	getImageClipRect () const

virtual ossim_uint32	degreesOfFreedom () const

virtual bool	needsInitialState () const
	needsInitialState() More...

virtual double	optimizeFit (const ossimTieGptSet &tieSet, double *targetVariance=0)

virtual ossimDpt	getForwardDeriv (int parmIdx, const ossimGpt &gpos, double hdelta=1e-11)

virtual ossimGpt	getInverseDeriv (int parmIdx, const ossimDpt &ipos, double hdelta=1e-11)

virtual ossimSensorModel::CovMatStatus	getObsCovMat (const ossimDpt &ipos, NEWMAT::SymmetricMatrix &Cov, const ossim_float64 defPointingSigma=0.5) const
	Gives 2X2 covariance matrix of observations. More...

virtual bool	isAffectedByElevation () const
	Implementation of pure virtual ossimProjection::isAffectedByElevation method. More...

void	computeGsd ()
	This method computes the ground sample distance(gsd) and sets class attributes theGSD and theMeanGSD by doing a lineSampleHeightToWorld on four points and calculating the distance from them. More...

virtual bool	getImageGeometry (const ossimString &, const ossimString &, ossimKeywordlist &) const
	Extracts geometry info from a non-ossim key,value pair to an ossim keyword list. More...

Public Member Functions inherited from ossimProjection
	ossimProjection ()

virtual	~ossimProjection ()

virtual ossimDpt	forward (const ossimGpt &wp) const

virtual ossimGpt	inverse (const ossimDpt &pp) const

virtual void	getRoundTripError (const ossimDpt &imagePoint, ossimDpt &errorResult) const

virtual void	getRoundTripError (const ossimGpt &groundPoint, ossimDpt &errorResult) const

virtual void	getGroundClipPoints (ossimGeoPolygon &gpts) const

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

virtual bool	operator!= (const ossimProjection &projection) 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

Public Member Functions inherited from ossimErrorStatusInterface
	ossimErrorStatusInterface ()

virtual	~ossimErrorStatusInterface ()

virtual ossimErrorCode	getErrorStatus () const

virtual ossimString	getErrorStatusString () const

virtual void	setErrorStatus (ossimErrorCode error_status) const

virtual void	setErrorStatus () const

virtual void	clearErrorStatus () const

bool	hasError () const

Public Member Functions inherited from ossimOptimizableProjection
	ossimOptimizableProjection ()

	ossimOptimizableProjection (const ossimOptimizableProjection &source)

virtual	~ossimOptimizableProjection ()

virtual ossimOptimizableProjection &	operator= (const ossimOptimizableProjection &source)

Public Member Functions inherited from ossimAdjustableParameterInterface
	ossimAdjustableParameterInterface ()

	ossimAdjustableParameterInterface (const ossimAdjustableParameterInterface &rhs)

virtual	~ossimAdjustableParameterInterface ()

void	newAdjustment (ossim_uint32 numberOfParameters=0)

void	setAdjustmentDescription (const ossimString &description)

ossimString	getAdjustmentDescription () const

ossimString	getAdjustmentDescription (ossim_uint32 adjustmentIdx) const
	Returns adjustmentDescription of specific adjustmentInfo. More...

void	setCurrentAdjustment (ossim_uint32 adjustmentIndex, bool notify=false)

bool	setCurrentAdjustment (const ossimString &description, bool notify=false)
	Sets the current adjustment to the adjustment with a matching description. More...

void	eraseAdjustment (bool notify)

void	eraseAdjustment (ossim_uint32 idx, bool notify)

void	resetAdjustableParameters (bool notify=false)

void	copyAdjustment (ossim_uint32 idx, bool notify)

void	copyAdjustment (bool notify=false)

void	keepAdjustment (ossim_uint32 idx, bool createCopy)

virtual void	keepAdjustment (bool createCopy=true)

const ossimAdjustableParameterInterface &	operator= (const ossimAdjustableParameterInterface &rhs)

void	removeAllAdjustments ()

ossim_uint32	getNumberOfAdjustableParameters () const

double	getAdjustableParameter (ossim_uint32 idx) const

virtual void	setAdjustableParameter (ossim_uint32 idx, double value, bool notify=false)

virtual void	setAdjustableParameter (ossim_uint32 idx, double value, double sigma, bool notify=false)

double	getParameterSigma (ossim_uint32 idx) const

void	setParameterSigma (ossim_uint32 idx, double value, bool notify=false)

ossimUnitType	getParameterUnit (ossim_uint32 idx) const

void	setParameterUnit (ossim_uint32 idx, ossimUnitType unit)

void	setParameterUnit (ossim_uint32 idx, const ossimString &unit)

void	setParameterCenter (ossim_uint32 idx, double center, bool notify=false)

double	getParameterCenter (ossim_uint32 idx) const

double	computeParameterOffset (ossim_uint32 idx) const

void	setParameterOffset (ossim_uint32 idx, ossim_float64 value, bool notify=false)

ossimString	getParameterDescription (ossim_uint32 idx) const

void	setParameterDescription (ossim_uint32 idx, const ossimString &descrption)

ossim_int32	findParameterIdxGivenDescription (ossim_uint32 adjustmentIdx, const ossimString &name) const

ossim_int32	findParameterIdxContainingDescription (ossim_uint32 adjustmentIdx, const ossimString &name) const

bool	isParameterLocked (ossim_uint32 idx) const

void	setParameterLockFlag (ossim_uint32 idxParam, bool flag)

bool	getParameterLockFlag (ossim_uint32 idx) const

void	lockAllParametersCurrentAdjustment ()

void	unlockAllParametersCurrentAdjustment ()

void	lockAllParameters (ossim_uint32 idxAdjustment)

void	unlockAllParameters (ossim_uint32 idxAdjustment)

void	resizeAdjustableParameterArray (ossim_uint32 numberOfParameters)

void	setAdjustment (const ossimAdjustmentInfo &adj, bool notify=false)

void	setAdjustment (ossim_uint32 idx, const ossimAdjustmentInfo &adj, bool notify=false)

void	addAdjustment (const ossimAdjustmentInfo &adj, bool notify)

void	getAdjustment (ossimAdjustmentInfo &adj) const

void	getAdjustment (ossim_uint32 idx, ossimAdjustmentInfo &adj) const

ossim_uint32	getNumberOfAdjustments () const

ossim_uint32	getCurrentAdjustmentIdx () const

void	setDirtyFlag (bool flag=true)

void	setAllDirtyFlag (bool flag=true)

bool	hasDirtyAdjustments () const

virtual void	saveCurrentAdjustmentOnly (ossimKeywordlist &kwl, const ossimString &prefix=ossimString(""))
	Saves the current active adjustment to the KWL file. More...

bool	saveAdjustments (ossimKeywordlist &kwl, const ossimString &prefix=ossimString("")) const
	Save all adjustments to the KWL file. More...

bool	loadAdjustments (const ossimKeywordlist &kwl, const ossimString &prefix=ossimString(""))

std::ostream &	print (std::ostream &out) const
	Dumps the currently active adjustment to ostream. More...

Protected Attributes
NEWMAT::Matrix	theCompositeMatrix
	uses file path to init model More...

NEWMAT::Matrix	theCompositeMatrixInverse

double	theOmega

double	thePhi

double	theKappa

double	theBoreSightTx

double	theBoreSightTy

double	theBoreSightTz

ossimDpt	thePrincipalPoint

ossimDpt	thePixelSize

double	theFocalLength

ossimEcefPoint	theEcefPlatformPosition

ossimGpt	thePlatformPosition

ossim_int32	theUtmZone

ossim_int8	theUtmHemisphere

ossimDpt3d	theUtmPlatformPosition

ossimEcefVector	theShiftValues

ossimRefPtr< ossimMeanRadialLensDistortion >	theLensDistortion

ossimEcefPoint	theAdjEcefPlatformPosition

Protected Attributes inherited from ossimSensorModel
ossimIpt	theImageSize

ossimDpt	theSubImageOffset

ossimString	theImageID

ossimString	theSensorID

ossimDpt	theGSD

ossim_float64	theMeanGSD

ossimGpt	theRefGndPt

ossimDpt	theRefImgPt

ossimPolygon	theBoundGndPolygon

ossimDrect	theImageClipRect

ossim_float64	theRelPosError

ossim_float64	theNominalPosError

ossimDpt	theParWRTx
	Partials for current point. More...

ossimDpt	theParWRTy

ossimDpt	theParWRTz

ossimDpt	theObs
	Observations & residuals for current point. More...

ossimDpt	theResid

ossimRefPtr< ossimProjection >	theSeedFunction
	Used as an initial guess for iterative solutions and a guess for points outside the support bounds. More...

bool	theExtrapolateImageFlag

bool	theExtrapolateGroundFlag

Protected Attributes inherited from ossimErrorStatusInterface
ossimErrorCode	theErrorStatus

Additional Inherited Members
Public Types inherited from ossimSensorModel
enum	CovMatStatus { COV_INVALID = 0, COV_PARTIAL = 1, COV_FULL = 2 }

enum	DeriveMode { OBS_INIT =-99, EVALUATE =-98, P_WRT_X = -1, P_WRT_Y = -2, P_WRT_Z = -3 }

Static Public Member Functions inherited from ossimSensorModel
static void	writeGeomTemplate (ostream &os)

Protected Member Functions inherited from ossimSensorModel
virtual	~ossimSensorModel ()

virtual ossimDpt	extrapolate (const ossimGpt &gp) const

virtual ossimGpt	extrapolate (const ossimDpt &ip, const double &height=ossim::nan()) const

void	buildNormalEquation (const ossimTieGptSet &tieSet, NEWMAT::SymmetricMatrix &A, NEWMAT::ColumnVector &residue, NEWMAT::ColumnVector &projResidue, double pstep_scale)

NEWMAT::ColumnVector	getResidue (const ossimTieGptSet &tieSet)

NEWMAT::ColumnVector	solveLeastSquares (NEWMAT::SymmetricMatrix &A, NEWMAT::ColumnVector &r) const

NEWMAT::Matrix	invert (const NEWMAT::Matrix &m) const
	stable invert stolen from ossimRpcSolver More...

Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Member Functions inherited from ossimAdjustableParameterInterface
bool	paramChanged (ossim_uint32 param_idx) const
	Returns true if specified parameter has been modified since last setAllChangeFlag(false) call. More...

void	setAllChangeFlags (bool areChanged)
	Sets all the change flags to the boolean indicated to indicate parameters are changed (TRUE) or not (FALSE). More...

void	initChangeFlags ()
	Initializes the change flags to TRUE. More...

Detailed Description

Definition at line 16 of file ossimApplanixUtmModel.h.

Constructor & Destructor Documentation

◆ ossimApplanixUtmModel() [1/2]

ossimApplanixUtmModel::ossimApplanixUtmModel ( )

Definition at line 32 of file ossimApplanixUtmModel.cpp.

References ossimMatrix4x4::createIdentity(), initAdjustableParameters(), theCompositeMatrix, theCompositeMatrixInverse, theEcefPlatformPosition, theFocalLength, ossimSensorModel::theGSD, theLensDistortion, ossimSensorModel::theMeanGSD, thePixelSize, ossimDpt::x, and ossimDpt::y.

Referenced by dup().

    :theOmega(0.0),
     thePhi(0.0),
     theKappa(0.0),
     theBoreSightTx(0.0),
     theBoreSightTy(0.0),
     theBoreSightTz(0.0)
 {
    theCompositeMatrix          = ossimMatrix4x4::createIdentity();
    theCompositeMatrixInverse   = ossimMatrix4x4::createIdentity();
    theFocalLength              = 55.0;
    thePixelSize = ossimDpt(.009, .009);
    theEcefPlatformPosition = ossimGpt(0.0,0.0, 1000.0);
    theGSD.x = 0.1524;
    theGSD.y = 0.1524;
    theMeanGSD = 0.1524;
    theLensDistortion = new ossimMeanRadialLensDistortion;
    initAdjustableParameters();
 
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << "ossimApplanixUtmModel::ossimApplanixUtmModel DEBUG:" << endl;
 #ifdef OSSIM_ID_ENABLED
       ossimNotify(ossimNotifyLevel_DEBUG)<< "OSSIM_ID:  " << OSSIM_ID << endl;
 #endif
    }
 }

◆ ossimApplanixUtmModel() [2/2]

ossimApplanixUtmModel::ossimApplanixUtmModel ( const ossimApplanixUtmModel & src )

Definition at line 61 of file ossimApplanixUtmModel.cpp.

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

    :ossimSensorModel(src),
     theCompositeMatrix(src.theCompositeMatrix),
     theCompositeMatrixInverse(src.theCompositeMatrixInverse),
     theOmega(src.theOmega),
     thePhi(src.thePhi),
     theKappa(src.theKappa),
     theBoreSightTx(src.theBoreSightTx),
     theBoreSightTy(src.theBoreSightTy),
     theBoreSightTz(src.theBoreSightTz),
     thePrincipalPoint(src.thePrincipalPoint),
     thePixelSize(src.thePixelSize),
     theFocalLength(src.theFocalLength),
     theEcefPlatformPosition(src.theEcefPlatformPosition),
     thePlatformPosition(src.thePlatformPosition),
     theUtmZone(src.theUtmZone),
     theUtmHemisphere(src.theUtmHemisphere),
     theUtmPlatformPosition(src.theUtmPlatformPosition)
 {
    if(src.theLensDistortion.valid())
    {
       theLensDistortion = new ossimMeanRadialLensDistortion(*(src.theLensDistortion.get()));
    }
    else
    {
       theLensDistortion = new ossimMeanRadialLensDistortion();
    }
 }

Member Function Documentation

◆ dup()

ossimObject * ossimApplanixUtmModel::dup ( ) const

virtual

Implements ossimProjection.

Definition at line 90 of file ossimApplanixUtmModel.cpp.

References ossimApplanixUtmModel().

 {
    return new ossimApplanixUtmModel(*this);
 }

◆ imagingRay()

void ossimApplanixUtmModel::imagingRay	(	const ossimDpt &	image_point,
		ossimEcefRay &	image_ray
	)		const

virtual

METHOD: imagingRay(image_point, &ossimEcefRay) Given an image point, returns a ray originating at some arbitrarily high point (ideally at the sensor position) and pointing towards the target. Implemented here but should be overriden for more efficient solution.

Reimplemented from ossimSensorModel.

Definition at line 95 of file ossimApplanixUtmModel.cpp.

References ossimEcefVector::magnitude(), ossimEcefRay::setDirection(), ossimEcefRay::setOrigin(), theAdjEcefPlatformPosition, theCompositeMatrix, theFocalLength, theLensDistortion, thePixelSize, thePrincipalPoint, ossimSensorModel::theRefImgPt, ossimMeanRadialLensDistortion::undistort(), ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

Referenced by lineSampleHeightToWorld().

 {
 //    if(traceDebug())
 //    {
 //       ossimNotify(ossimNotifyLevel_DEBUG) << "ossimApplanixUtmModel::imagingRay: ..... entered" << std::endl;
 //    }
     ossimDpt f1 ((image_point) - theRefImgPt);
    f1.x *= thePixelSize.x;
    f1.y *= -thePixelSize.y;
    ossimDpt film (f1 - thePrincipalPoint);
 //    if(traceDebug())
 //    {
 //       ossimNotify(ossimNotifyLevel_DEBUG) << "pixel size   = " << thePixelSize << std::endl;
 //       ossimNotify(ossimNotifyLevel_DEBUG) << "principal pt = " << thePrincipalPoint << std::endl;
 //       ossimNotify(ossimNotifyLevel_DEBUG) << "film pt      = " << film << std::endl;
 //    }
    if (theLensDistortion.valid())
    {
       ossimDpt filmOut;
       theLensDistortion->undistort(film, filmOut);
       film = filmOut;
    }
    
    ossimColumnVector3d cam_ray_dir (film.x,
                                     film.y,
                                     -theFocalLength);
    ossimEcefVector     ecf_ray_dir (theCompositeMatrix*cam_ray_dir);
    ecf_ray_dir = ecf_ray_dir*(1.0/ecf_ray_dir.magnitude());
    
    image_ray.setOrigin(theAdjEcefPlatformPosition);
    image_ray.setDirection(ecf_ray_dir);
    
 //    if(traceDebug())
 //    {
 //       ossimNotify(ossimNotifyLevel_DEBUG) << "ossimApplanixUtmModel::imagingRay: ..... leaving" << std::endl;
 //    }
 }

◆ initAdjustableParameters()

void ossimApplanixUtmModel::initAdjustableParameters ( )

virtual

Reimplemented from ossimAdjustableParameterInterface.

Definition at line 320 of file ossimApplanixUtmModel.cpp.

References ossimAdjustableParameterInterface::resizeAdjustableParameterArray(), ossimAdjustableParameterInterface::setAdjustableParameter(), ossimAdjustableParameterInterface::setParameterDescription(), ossimAdjustableParameterInterface::setParameterSigma(), and ossimAdjustableParameterInterface::setParameterUnit().

Referenced by ossimApplanixUtmModel().

 {
    resizeAdjustableParameterArray(6);
    
    setAdjustableParameter(0, 0.0);
    setParameterDescription(0, "x_offset");
    setParameterUnit(0, "pixels");
 
    setAdjustableParameter(1, 0.0);
    setParameterDescription(1, "y_offset");
    setParameterUnit(1, "pixels");
 
    setAdjustableParameter(2, 0.0);
    setParameterDescription(2, "orientation x");
    setParameterUnit(2, "degrees");
 
    setAdjustableParameter(3, 0.0);
    setParameterDescription(3, "orientation y");
    setParameterUnit(3, "degrees");
 
    setAdjustableParameter(4, 0.0);
    setParameterDescription(4, "orientation z");
    setParameterUnit(4, "degrees");
 
    setAdjustableParameter(5, 0.0);
    setParameterDescription(5, "Altitude delta");
    setParameterUnit(5, "meters");
 
    
    
    setParameterSigma(0, 20.0);
    setParameterSigma(1, 20.0);
    setParameterSigma(2, .1);
    setParameterSigma(3, .1);
    setParameterSigma(4, .1);
    setParameterSigma(5, 50);
 }

◆ lineSampleHeightToWorld()

void ossimApplanixUtmModel::lineSampleHeightToWorld	(	const ossimDpt &	lineSampPt,
		const double &	heightEllipsoid,
		ossimGpt &	worldPt
	)		const

virtual

METHOD: lineSampleHeightToWorld This is the pure virtual that performs the actual work of projecting the image point to the given elevation above Ellipsoid.

Implements ossimSensorModel.

Definition at line 173 of file ossimApplanixUtmModel.cpp.

References imagingRay(), and ossimEcefRay::intersectAboveEarthEllipsoid().

 {
 //   if (!insideImage(image_point))
 //   {
       //worldPoint.makeNan();
 //       worldPoint = extrapolate(image_point, heightEllipsoid);
 //   }
 //   else
    {
       //***
       // First establish imaging ray from image point:
       //***
       ossimEcefRay ray;
       imagingRay(image_point, ray);
       ossimEcefPoint Pecf (ray.intersectAboveEarthEllipsoid(heightEllipsoid));
       worldPoint = ossimGpt(Pecf);
    }
 }

◆ lineSampleToWorld()

void ossimApplanixUtmModel::lineSampleToWorld	(	const ossimDpt &	image_point,
		ossimGpt &	world_point
	)		const

virtual

Implementation of base-class pure virtual projection methods. These methods may be overriden by derived classes if those have more efficient schemes. The implementations here are iterative (relatively slow). Both depend on calls to the pure virtual lineSampleHeightToWorld() method.

Reimplemented from ossimSensorModel.

Definition at line 134 of file ossimApplanixUtmModel.cpp.

Referenced by updateModel().

 {
    if (traceDebug())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimApplanixEcefModel::lineSampleToWorld:entering..." << std::endl;
    
    if(image_point.hasNans())
    {
       gpt.makeNan();
       return;
    }
    //***
    // Extrapolate if image point is outside image:
    //***
   // if (!insideImage(image_point))
   // {
    //   gpt.makeNan();
 //       gpt = extrapolate(image_point);
    //   return;
    //}
 
    //***
    // Determine imaging ray and invoke elevation source object's services to
    // intersect ray with terrain model:
    //***
    ossimEcefRay ray;
    imagingRay(image_point, ray);
    ossimElevManager::instance()->intersectRay(ray, gpt);
 
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "image_point = " << image_point << std::endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "ray = " << ray << std::endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "gpt = " << gpt << std::endl;
    }
 
    if (traceDebug())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimSensorModel::lineSampleToWorld: returning..." << std::endl;
    return;
 }

◆ loadState()

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

virtual

Method to the load (recreate) the state of the object from a keyword list. Return true if ok or false on error.

Reimplemented from ossimSensorModel.

Definition at line 443 of file ossimApplanixUtmModel.cpp.

Referenced by setupOptimizer().

 {
    if(traceDebug())
    {
       std::cout << "ossimApplanixUtmModel::loadState: ......... entered" << std::endl;
    }
    theImageClipRect = ossimDrect(0,0,4076,4091);
    theRefImgPt      = ossimDpt(2046.0, 2038.5);
    
    ossimSensorModel::loadState(kwl, prefix);
 
    if(getNumberOfAdjustableParameters() < 1)
    {
       initAdjustableParameters();
    }
 
    const char* eo_file           = kwl.find(prefix, "eo_file");
    const char* eo_id             = kwl.find(prefix, "eo_id");
    const char* omega             = kwl.find(prefix, "omega");
    const char* phi               = kwl.find(prefix, "phi");
    const char* kappa             = kwl.find(prefix, "kappa");
    const char* bore_sight_tx     = kwl.find(prefix, "bore_sight_tx");
    const char* bore_sight_ty     = kwl.find(prefix, "bore_sight_ty");
    const char* bore_sight_tz     = kwl.find(prefix, "bore_sight_tz");
    
    const char* principal_point   = kwl.find(prefix, "principal_point");
    const char* pixel_size        = kwl.find(prefix, "pixel_size");
    const char* focal_length      = kwl.find(prefix, "focal_length");
    const char* latlonh_platform_position = kwl.find(prefix, "latlonh_platform_position");
    const char* utm_platform_position = kwl.find(prefix, "utm_platform_position");
    const char* compute_gsd_flag  = kwl.find(prefix, "compute_gsd_flag");
    const char* utm_zone          = kwl.find(prefix, "utm_zone");
    const char* utm_hemisphere    = kwl.find(prefix, "utm_hemisphere");
    const char* camera_file       = kwl.find(prefix, "camera_file");
    const char* shift_values      = kwl.find(prefix, "shift_values");
    
    theCompositeMatrix          = ossimMatrix3x3::createIdentity();
    theCompositeMatrixInverse   = ossimMatrix3x3::createIdentity();
    theOmega                    = 0.0;
    thePhi                      = 0.0;
    theKappa                    = 0.0;
    theBoreSightTx              = 0.0;
    theBoreSightTy              = 0.0;
    theBoreSightTz              = 0.0;
    theFocalLength              = 55.0;
    thePixelSize = ossimDpt(.009, .009);
    theEcefPlatformPosition = ossimGpt(0.0,0.0, 1000.0);
 
    bool loadedFromEoFile = false;
 
    if(eo_id)
    {
       theImageID = eo_id;
    }
    // loading from standard eo file with given record id
    //
    if(eo_file)
    {
       ossimApplanixEOFile eoFile;
       
       if(eoFile.parseFile(ossimFilename(eo_file)))
       {
          ossimRefPtr<ossimApplanixEORecord> record = eoFile.getRecordGivenId(theImageID);
 
          if(record.valid())
          {
             loadedFromEoFile = true;
             theBoreSightTx = eoFile.getBoreSightTx()/60.0;
             theBoreSightTy = eoFile.getBoreSightTy()/60.0;
             theBoreSightTz = eoFile.getBoreSightTz()/60.0;
             theShiftValues = ossimEcefVector(eoFile.getShiftValuesX(),
                                              eoFile.getShiftValuesY(),
                                              eoFile.getShiftValuesZ());
             ossim_int32 easting  = eoFile.getFieldIdxLike("EASTING");
             ossim_int32 northing = eoFile.getFieldIdxLike("NORTHING");
             ossim_int32 height   = eoFile.getFieldIdxLike("HEIGHT");
             ossim_int32 omega    = eoFile.getFieldIdxLike("OMEGA");
             ossim_int32 phi      = eoFile.getFieldIdxLike("PHI");
             ossim_int32 kappa    = eoFile.getFieldIdxLike("KAPPA");
 
             if((omega>=0)&&
                (phi>=0)&&
                (kappa>=0)&&
                (height>=0)&&
                (easting>=0)&&
                (northing>=0))
             {
                theOmega = (*record)[omega].toDouble();
                thePhi   = (*record)[phi].toDouble();
                theKappa = (*record)[kappa].toDouble();
                double h = (*record)[height].toDouble();
                ossimString heightType = kwl.find(prefix, "height_type");
                if(eoFile.isUtmFrame())
                {
                   theUtmZone = eoFile.getUtmZone();
                   theUtmHemisphere = eoFile.getUtmHemisphere()=="North"?'N':'S';
                   ossimUtmProjection utmProj;
                   utmProj.setZone(theUtmZone);
                   utmProj.setHemisphere((char)theUtmHemisphere);
                   theUtmPlatformPosition.x = (*record)[easting].toDouble();
                   theUtmPlatformPosition.y = (*record)[northing].toDouble();
                   theUtmPlatformPosition.z = h;
                   thePlatformPosition = utmProj.inverse(ossimDpt(theUtmPlatformPosition.x,
                                                                  theUtmPlatformPosition.y));
                   thePlatformPosition.height(h);
 
                   if(eoFile.isHeightAboveMSL())
                   {
                      double offset = ossimGeoidManager::instance()->offsetFromEllipsoid(thePlatformPosition);
                      if(!ossim::isnan(offset))
                      {
                         thePlatformPosition.height(h + offset);
                         theUtmPlatformPosition.z = h + offset;
                      }
                   }
                }
                else
                {
                   return false;
                }
             }
             theEcefPlatformPosition = thePlatformPosition;
          }
          else
          {
             return false;
          }
       }
    }
    if(!loadedFromEoFile)
    {
       if(shift_values)
       {
          std::vector<ossimString> splitString;
          ossimString tempString(shift_values);
          tempString = tempString.trim();
          tempString.split(splitString, " " );
          if(splitString.size() == 3)
          {
             theShiftValues = ossimEcefVector(splitString[0].toDouble(),
                                              splitString[1].toDouble(),
                                              splitString[2].toDouble());
          }
       }
       if(omega&&phi&&kappa)
       {
          theOmega = ossimString(omega).toDouble();
          thePhi   = ossimString(phi).toDouble();
          theKappa = ossimString(kappa).toDouble();
       }
       if(bore_sight_tx&&bore_sight_ty&&bore_sight_tz)
       {
          theBoreSightTx = ossimString(bore_sight_tx).toDouble()/60.0;
          theBoreSightTy = ossimString(bore_sight_ty).toDouble()/60.0;
          theBoreSightTz = ossimString(bore_sight_tz).toDouble()/60.0;
       }
       double lat=0.0, lon=0.0, h=0.0;
       if(utm_zone)
       {
          theUtmZone = ossimString(utm_zone).toInt32();
       }
       if(utm_hemisphere)
       {
          ossimString hem = utm_hemisphere;
          hem = hem.trim();
          hem = hem.upcase();
          theUtmHemisphere = *(hem.begin());
       }
       if(utm_platform_position)
       {
          ossimUtmProjection utmProj;
          std::vector<ossimString> splitString;
          ossimString tempString(utm_platform_position);
          tempString = tempString.trim();
          ossimString datumString;
          utmProj.setZone(theUtmZone);
          utmProj.setHemisphere((char)theUtmHemisphere);
          tempString.split(splitString, " ");
          if(splitString.size() > 2)
          {
             theUtmPlatformPosition.x = splitString[0].toDouble();
             theUtmPlatformPosition.y = splitString[1].toDouble();
             theUtmPlatformPosition.z = splitString[2].toDouble();
          }
          if(splitString.size() > 3)
          {
             datumString = splitString[3];
          }
          const ossimDatum* datum = ossimDatumFactory::instance()->create(datumString);
          if(datum)
          {
             utmProj.setDatum(datum);
          }
          
          thePlatformPosition = utmProj.inverse(ossimDpt(theUtmPlatformPosition.x,
                                                         theUtmPlatformPosition.y));
          thePlatformPosition.height(theUtmPlatformPosition.z);
          
          ossimString heightType = kwl.find(prefix, "height_type");
          if(heightType == "msl")
          {
             double offset = ossimGeoidManager::instance()->offsetFromEllipsoid(thePlatformPosition);
             if(ossim::isnan(offset) == false)
             {
                thePlatformPosition.height(thePlatformPosition.height() + offset);
                theUtmPlatformPosition.z = thePlatformPosition.height();
             }
          }
          theEcefPlatformPosition = thePlatformPosition;
       }
       else if(latlonh_platform_position)
       {
          std::vector<ossimString> splitString;
          ossimString tempString(latlonh_platform_position);
          std::string datumString;
          tempString = tempString.trim();
          tempString.split(splitString, " ");
          if(splitString.size() > 2)
          {
             lat = splitString[0].toDouble();
             lon = splitString[1].toDouble();
             h = splitString[2].toDouble();
          }
          if(splitString.size() > 3)
          {
             datumString = splitString[3].string();
          }
          thePlatformPosition.latd(lat);
          thePlatformPosition.lond(lon);
          thePlatformPosition.height(h);
          const ossimDatum* datum = ossimDatumFactory::instance()->create(datumString);
          if(datum)
          {
             thePlatformPosition.datum(datum);
          }
          ossimString heightType = kwl.find(prefix, "height_type");
          if(heightType == "msl")
          {
             double offset = ossimGeoidManager::instance()->offsetFromEllipsoid(thePlatformPosition);
             if(ossim::isnan(offset) == false)
             {
                thePlatformPosition.height(thePlatformPosition.height() + offset);
             }
          }
          theEcefPlatformPosition = thePlatformPosition;
       }
    }
    if(principal_point)
    {
       std::vector<ossimString> splitString;
       ossimString tempString(principal_point);
       tempString = tempString.trim();
       tempString.split(splitString, " ");
       if(splitString.size() == 2)
       {
          thePrincipalPoint.x = splitString[0].toDouble();
          thePrincipalPoint.y = splitString[1].toDouble();
       }
    }
    if(pixel_size)
    {
       std::vector<ossimString> splitString;
       ossimString tempString(principal_point);
       tempString = tempString.trim();
       tempString.split(splitString, " ");
       if(splitString.size() == 2)
       {
          thePixelSize.x = splitString[0].toDouble();
          thePixelSize.y = splitString[1].toDouble();
       }
    }
    if(focal_length)
    {
       theFocalLength = ossimString(focal_length).toDouble();
    }
 
    if(camera_file)
    {
       ossimKeywordlist cameraKwl;
       ossimKeywordlist lensKwl;
       cameraKwl.add(camera_file);
       const char* sensor = cameraKwl.find("sensor");
       const char* image_size      = cameraKwl.find(prefix, "image_size");
       principal_point = cameraKwl.find("principal_point");
       focal_length    = cameraKwl.find("focal_length");
       pixel_size      = cameraKwl.find(prefix, "pixel_size");
       focal_length    = cameraKwl.find(prefix, "focal_length");
       const char* distortion_units = cameraKwl.find(prefix, "distortion_units");
       ossimUnitConversionTool tool;
       ossimUnitType unitType = OSSIM_MILLIMETERS;
 
       if(distortion_units)
       {
          unitType = (ossimUnitType)ossimUnitTypeLut::instance()->getEntryNumber(distortion_units);
 
          if(unitType == OSSIM_UNIT_UNKNOWN)
          {
             unitType = OSSIM_MILLIMETERS;
          }
       }
       if(image_size)
       {
          std::vector<ossimString> splitString;
          ossimString tempString(image_size);
          tempString = tempString.trim();
          tempString.split(splitString, " ");
          double w=1, h=1;
          if(splitString.size() == 2)
          {
             w = splitString[0].toDouble();
             h = splitString[1].toDouble();
          }
          theImageClipRect = ossimDrect(0,0,w-1,h-1);
          theRefImgPt      = ossimDpt(w/2.0, h/2.0);
          theImageSize     = ossimDpt(w, h);
       }
       if(sensor)
       {
          theSensorID = sensor;
       }
       if(principal_point)
       {
          std::vector<ossimString> splitString;
          ossimString tempString(principal_point);
          tempString = tempString.trim();
          tempString.split(splitString, " ");
          if(splitString.size() == 2)
          {
             thePrincipalPoint.x = splitString[0].toDouble();
             thePrincipalPoint.y = splitString[1].toDouble();
          }
       }
       if(pixel_size)
       {
          std::vector<ossimString> splitString;
          ossimString tempString(pixel_size);
          tempString = tempString.trim();
          tempString.split(splitString, " ");
          if(splitString.size() == 1)
          {
             thePixelSize.x = splitString[0].toDouble();
             thePixelSize.y = thePixelSize.x;
          }
          else if(splitString.size() == 2)
          {
             thePixelSize.x = splitString[0].toDouble();
             thePixelSize.y = splitString[1].toDouble();
          }
       }
       if(focal_length)
       {
          theFocalLength = ossimString(focal_length).toDouble();
       }
 
       cameraKwl.trimAllValues();
       
       
       ossimString regExpression =  ossimString("^(") + "d[0-9]+)";
       vector<ossimString> keys;
       cameraKwl.getSubstringKeyList( keys, regExpression );
       long numberOfDistortions = (long)keys.size();
       int offset = (int)ossimString("d").size();
       ossim_uint32 idx = 0;
       std::vector<int> numberList(numberOfDistortions);
       for(idx = 0; idx < numberList.size();++idx)
       {
          ossimString numberStr(keys[idx].begin() + offset,
                                keys[idx].end());
          numberList[idx] = numberStr.toInt();
       }
       std::sort(numberList.begin(), numberList.end());
       double distance=0.0, distortion=0.0;
      
       for(idx = 0; idx < numberList.size(); ++idx)
       {
          ossimString value = cameraKwl.find(ossimString("d")+ossimString::toString(numberList[idx]));
          if(!value.empty())
          {
             std::istringstream inStr(value.c_str());
             inStr >> distance;
             ossim::skipws(inStr);
             inStr >> distortion;
 #if 0
             std::vector<ossimString> splitString;
             ossimString tempString(value);
             tempString = tempString.trim();
             tempString.split(splitString, " ");
             std::cout << splitString.size() << std::endl;
             if(splitString.size() >= 2)
             {
                distance = splitString[0].toDouble();
                distortion = splitString[1].toDouble();
             }
 #endif
             
             tool.setValue(distortion, unitType);
             lensKwl.add(ossimString("distance") + ossimString::toString(idx),
                         distance,
                         true);
             lensKwl.add(ossimString("distortion") + ossimString::toString(idx),
                         tool.getMillimeters(),
                         true);
          }
          lensKwl.add("convergence_threshold",
                      .00001,
                      true);
          if(pixel_size)
          {
             lensKwl.add("dxdy",
                         ossimString::toString(thePixelSize.x) + " " +ossimString::toString(thePixelSize.y),
                         true);
          }
          else
          {
             lensKwl.add("dxdy",
                         ".009 .009",
                         true);
          }
       }
       if(theLensDistortion.valid())
       {
          theLensDistortion->loadState(lensKwl,"");
       }
    }
    else
    {
       
       if(principal_point)
       {
          std::vector<ossimString> splitString;
          ossimString tempString(principal_point);
          tempString = tempString.trim();
          tempString.split(splitString, " ");
          if(splitString.size() >= 2)
          {
             thePrincipalPoint.x = splitString[0].toDouble();
             thePrincipalPoint.y = splitString[1].toDouble();
          }
       }
       else
       {
          if(traceDebug())
          {
             ossimNotify(ossimNotifyLevel_DEBUG) << "No principal_point found" << std::endl;
             return false;
          }
       }
       if(pixel_size)
       {
          std::vector<ossimString> splitString;
          ossimString tempString(pixel_size);
          tempString = tempString.trim();
          tempString.split(splitString, " ");
          if(splitString.size() == 1)
          {
             thePixelSize.x = splitString[0].toDouble();
             thePixelSize.y = thePixelSize.x;
          }
          else if(splitString.size() == 2)
          {
             thePixelSize.x = splitString[0].toDouble();
             thePixelSize.y = splitString[1].toDouble();
          }
       }
       else
       {
          if(traceDebug())
          {
             ossimNotify(ossimNotifyLevel_DEBUG) << "No pixel size found" << std::endl;
             return false;
          }
       }
       if(focal_length)
       {
          theFocalLength = ossimString(focal_length).toDouble();
       }
       else
       {
          if(traceDebug())
          {
             ossimNotify(ossimNotifyLevel_DEBUG) << "No focal length found" << std::endl;
             return false;
          }
       }
       
       if(theLensDistortion.valid())
       {
          ossimString lensPrefix = ossimString(prefix)+"distortion.";
          theLensDistortion->loadState(kwl,
                                       lensPrefix.c_str());
       }
    }
    theRefGndPt = thePlatformPosition;
    
    updateModel();
 
    lineSampleToWorld(theRefImgPt, theRefGndPt);
    if(compute_gsd_flag)
    {
       if(ossimString(compute_gsd_flag).toBool())
       {
          ossimGpt right;
          ossimGpt top;
          lineSampleToWorld(theRefImgPt + ossimDpt(1.0, 0.0),
                            right);
          lineSampleToWorld(theRefImgPt + ossimDpt(1.0, 0.0),
                            top);
 
          ossimEcefVector horizontal = ossimEcefPoint(theRefGndPt)-ossimEcefPoint(right);
          ossimEcefVector vertical   = ossimEcefPoint(theRefGndPt)-ossimEcefPoint(top);
 
          theGSD.x = horizontal.length();
          theGSD.y = vertical.length();
          theMeanGSD = (theGSD.x+theGSD.y)*.5;
       }
    }
    if(traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "theOmega:              " << theOmega << std::endl
                                           << "thePhi:                " << thePhi   << std::endl
                                           << "theKappa:              " << theKappa << std::endl;
       std::cout << "platform position:     " << thePlatformPosition << std::endl;
       std::cout << "platform position ECF: " << theEcefPlatformPosition << std::endl;
       std::cout << "ossimApplanixModel::loadState: ......... leaving" << std::endl;
    }
 
    return true;
 }

◆ saveState()

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

virtual

METHODS: saveState, loadState Fulfills ossimObject base-class pure virtuals. Loads and saves geometry KWL files.

Reimplemented from ossimSensorModel.

Definition at line 387 of file ossimApplanixUtmModel.cpp.

References ossimKeywordlist::add(), ossimString::c_str(), ossimDatum::code(), ossimGpt::datum(), ossimGpt::height(), ossimGpt::latd(), ossimGpt::lond(), ossimMeanRadialLensDistortion::saveState(), ossimSensorModel::saveState(), theBoreSightTx, theBoreSightTy, theBoreSightTz, theEcefPlatformPosition, theFocalLength, theKappa, theLensDistortion, theOmega, thePhi, thePixelSize, thePlatformPosition, thePrincipalPoint, theShiftValues, theUtmHemisphere, theUtmPlatformPosition, theUtmZone, ossimString::toString(), ossimRefPtr< T >::valid(), ossimEcefPoint::x(), ossimEcefVector::x(), ossimDpt3d::x, ossimDpt::x, ossimEcefPoint::y(), ossimEcefVector::y(), ossimDpt3d::y, ossimDpt::y, ossimEcefPoint::z(), ossimEcefVector::z(), and ossimDpt3d::z.

 {
    ossimSensorModel::saveState(kwl, prefix);
    
    kwl.add(prefix, "type", "ossimApplanixUtmModel", true);
 
    kwl.add(prefix, "omega", theOmega, true);
    kwl.add(prefix, "phi", thePhi, true);
    kwl.add(prefix, "kappa", theKappa, true);
    kwl.add(prefix, "bore_sight_tx", theBoreSightTx*60, true);
    kwl.add(prefix, "bore_sight_ty", theBoreSightTy*60, true);
    kwl.add(prefix, "bore_sight_tz", theBoreSightTz*60, true);
    kwl.add(prefix, "principal_point", ossimString::toString(thePrincipalPoint.x) + " " + ossimString::toString(thePrincipalPoint.y));
    kwl.add(prefix, "pixel_size",      ossimString::toString(thePixelSize.x) + " " + ossimString::toString(thePixelSize.y));
    kwl.add(prefix, "focal_length", theFocalLength);
    kwl.add(prefix, "ecef_platform_position",
            ossimString::toString(theEcefPlatformPosition.x()) + " " +
            ossimString::toString(theEcefPlatformPosition.y()) + " " +
            ossimString::toString(theEcefPlatformPosition.z()));
    kwl.add(prefix, "latlonh_platform_position",
            ossimString::toString(thePlatformPosition.latd()) + " " +
            ossimString::toString(thePlatformPosition.lond()) + " " +
            ossimString::toString(thePlatformPosition.height()) + " " +
            thePlatformPosition.datum()->code());
    kwl.add(prefix,
            "utm_platform_position",
            ossimString::toString(theUtmPlatformPosition.x) + " " +
            ossimString::toString(theUtmPlatformPosition.y) + " " +
            ossimString::toString(theUtmPlatformPosition.z) + " " +
            thePlatformPosition.datum()->code(),
            true);
    kwl.add(prefix,
            "utm_zone",
            theUtmZone,
            true);
    kwl.add(prefix,
            "utm_hemisphere",
            theUtmHemisphere,
            true);
    kwl.add(prefix,
            "shift_values",
            ossimString::toString(theShiftValues.x()) + " " + 
            ossimString::toString(theShiftValues.y()) + " " +
            ossimString::toString(theShiftValues.z()),
            true);
    if(theLensDistortion.valid())
    {
       ossimString lensPrefix = ossimString(prefix)+"distortion.";
       theLensDistortion->saveState(kwl,
                                    lensPrefix.c_str());
    }
    
    return true;
 }

◆ setFocalLength()

void ossimApplanixUtmModel::setFocalLength ( double focalLength )

Definition at line 374 of file ossimApplanixUtmModel.cpp.

References theFocalLength.

 {
    theFocalLength = focalLength;
 }

◆ setImageRect()

void ossimApplanixUtmModel::setImageRect ( const ossimDrect & rect )

Definition at line 368 of file ossimApplanixUtmModel.cpp.

References ossimDrect::midPoint(), ossimSensorModel::theImageClipRect, and ossimSensorModel::theRefImgPt.

 {
    theImageClipRect = rect;
    theRefImgPt = rect.midPoint();
 }

◆ setPixelSize()

void ossimApplanixUtmModel::setPixelSize ( const ossimDpt & pixelSize )

Definition at line 363 of file ossimApplanixUtmModel.cpp.

References thePixelSize.

 {
    thePixelSize = pixelSize;
 }

◆ setPlatformPosition()

void ossimApplanixUtmModel::setPlatformPosition ( const ossimGpt & gpt )

Definition at line 379 of file ossimApplanixUtmModel.cpp.

References theEcefPlatformPosition, ossimSensorModel::theRefGndPt, and updateModel().

 {
    theRefGndPt             = gpt;
    theEcefPlatformPosition = gpt;
    updateModel();
    
 }

◆ setPrincipalPoint()

void ossimApplanixUtmModel::setPrincipalPoint ( ossimDpt principalPoint )

Definition at line 358 of file ossimApplanixUtmModel.cpp.

References thePrincipalPoint.

 {
    thePrincipalPoint = principalPoint;
 }

◆ setupOptimizer()

bool ossimApplanixUtmModel::setupOptimizer ( const ossimString & init_file )

virtual

ground to image faster (you don't need DEM) //TBC

Reimplemented from ossimOptimizableProjection.

Definition at line 973 of file ossimApplanixUtmModel.cpp.

References ossimKeywordlist::addFile(), ossimString::c_str(), and loadState().

 {
    ossimKeywordlist kwl;
    kwl.addFile(init_file.c_str());
 
    return loadState(kwl);
 }

◆ updateModel()

void ossimApplanixUtmModel::updateModel ( )

virtual

VIRTUAL METHOD: updateModel() Following a change to the adjustable parameter set, this virtual is called to permit instances to compute derived quantities after parameter change.

Reimplemented from ossimSensorModel.

Definition at line 234 of file ossimApplanixUtmModel.cpp.

References ossimAdjustableParameterInterface::computeParameterOffset(), ossimMatrix3x3::createRotationXMatrix(), ossimMatrix3x3::createRotationYMatrix(), ossimMatrix3x3::createRotationZMatrix(), ossimEcefVector::cross(), ossimUtmProjection::forward(), ossimGpt::height(), ossimUtmProjection::inverse(), ossimGpt::latd(), lineSampleToWorld(), ossimDrect::ll(), ossimGpt::lond(), ossimDrect::lr(), ossimLsrSpace::lsrToEcefRotMatrix(), ossimGpt::metersPerDegree(), ossimEcefVector::normalize(), OSSIM_LEFT_HANDED, ossimPolygon::resize(), ossimMapProjection::setMetersPerPixel(), ossimUtmProjection::setZone(), theAdjEcefPlatformPosition, ossimSensorModel::theBoundGndPolygon, theCompositeMatrix, theCompositeMatrixInverse, ossimSensorModel::theImageClipRect, theKappa, ossimSensorModel::theMeanGSD, theOmega, thePhi, thePlatformPosition, theUtmHemisphere, theUtmZone, ossimDrect::ul(), ossimDrect::ur(), and ossimDpt::x.

Referenced by setPlatformPosition().

 {
 
    ossimGpt wgs84Pt;
    double metersPerDegree = wgs84Pt.metersPerDegree().x;
    double degreePerMeter = 1.0/metersPerDegree;
    double latShift = -computeParameterOffset(1)*theMeanGSD*degreePerMeter;
    double lonShift = computeParameterOffset(0)*theMeanGSD*degreePerMeter;
 
    ossimGpt gpt = thePlatformPosition;
 //   gpt.height(0.0);
    double height = gpt.height();
    gpt.height(height + computeParameterOffset(5));
    gpt.latd(gpt.latd() + latShift);
    gpt.lond(gpt.lond() + lonShift);
    
    theAdjEcefPlatformPosition = gpt;
    ossimLsrSpace lsrSpace(theAdjEcefPlatformPosition);
    // ORIENT TO A UTM AXIS
    //
     NEWMAT::ColumnVector v(3);
    
     v[0] = 0.0;
     v[1] = 0.0;
     v[2] = 1.0;
     NEWMAT::ColumnVector v2 = lsrSpace.lsrToEcefRotMatrix()*v;
     ossimEcefVector zVector(v2[0], v2[1], v2[2]);
     zVector.normalize();
    
    // now lets create a UTM axis by creating a derivative at the center
    // by shift over a few pixels and subtracting
    //
    ossimUtmProjection utmProj;
    
    utmProj.setZone(theUtmZone);
    utmProj.setZone(theUtmHemisphere);
    utmProj.setMetersPerPixel(ossimDpt(1.0,1.0));
    ossimDpt midPt  = utmProj.forward(theAdjEcefPlatformPosition);
    ossimDpt rPt    = midPt + ossimDpt(10, 0.0);
    ossimDpt uPt    = midPt + ossimDpt(0.0, 10.0);
    ossimGpt wMidPt = utmProj.inverse(midPt);
    ossimGpt wRPt   = utmProj.inverse(rPt);
    ossimGpt wUPt   = utmProj.inverse(uPt);
    
    ossimEcefPoint ecefMidPt = wMidPt;
    ossimEcefPoint ecefRPt   = wRPt;
    ossimEcefPoint ecefUPt   = wUPt;
    
    ossimEcefVector east  = ecefRPt-ecefMidPt;
    ossimEcefVector north = ecefUPt-ecefMidPt;
    east.normalize();
    north.normalize();
    
    // now use the lsr space constructors to construct an orthogonal set of axes
    //
    lsrSpace = ossimLsrSpace(thePlatformPosition,
                             0,
                             north,
                             east.cross(north));
 //   lsrSpace = ossimLsrSpace(thePlatformPosition);
    // DONE ORIENT TO UTM AXIS
    
    NEWMAT::Matrix platformLsrMatrix = lsrSpace.lsrToEcefRotMatrix();
    NEWMAT::Matrix orientationMatrix = (ossimMatrix3x3::createRotationXMatrix(theOmega+computeParameterOffset(2), OSSIM_LEFT_HANDED)*
                                        ossimMatrix3x3::createRotationYMatrix(thePhi+computeParameterOffset(3), OSSIM_LEFT_HANDED)*
                                        ossimMatrix3x3::createRotationZMatrix(theKappa+computeParameterOffset(4), OSSIM_LEFT_HANDED));
    
    theCompositeMatrix         = platformLsrMatrix*orientationMatrix;
    theCompositeMatrixInverse  = theCompositeMatrix.i();
 
 //   theAdjEcefPlatformPosition = theEcefPlatformPosition; 
 
    theBoundGndPolygon.resize(4);
    // ossim_float64 w = theImageClipRect.width();//*2.0;
    // ossim_float64 h = theImageClipRect.height();//*2.0;
    
    lineSampleToWorld(theImageClipRect.ul(),gpt);//+ossimDpt(-w, -h), gpt);
    theBoundGndPolygon[0] = gpt;
    lineSampleToWorld(theImageClipRect.ur(),gpt);//+ossimDpt(w, -h), gpt);
    theBoundGndPolygon[1] = gpt;
    lineSampleToWorld(theImageClipRect.lr(),gpt);//+ossimDpt(w, h), gpt);
    theBoundGndPolygon[2] = gpt;
    lineSampleToWorld(theImageClipRect.ll(),gpt);//+ossimDpt(-w, h), gpt);
    theBoundGndPolygon[3] = gpt;
 }

◆ useForward()

virtual bool ossimApplanixUtmModel::useForward ( ) const

inlinevirtual

ossimOptimizableProjection

Implements ossimOptimizableProjection.

Definition at line 55 of file ossimApplanixUtmModel.h.

55 {return true;}

◆ worldToLineSample()

void ossimApplanixUtmModel::worldToLineSample	(	const ossimGpt &	worldPoint,
		ossimDpt &	lineSampPt
	)		const

virtual

METHOD: worldToLineSample() Performs the forward projection from ground point to line, sample.

Reimplemented from ossimSensorModel.

Definition at line 194 of file ossimApplanixUtmModel.cpp.

References ossimEcefVector::data(), ossimLensDistortion::distort(), ossimPolygon::getNumberOfVertices(), ossimPolygon::hasNans(), ossimDpt::makeNan(), ossimPolygon::pointWithin(), theAdjEcefPlatformPosition, ossimSensorModel::theBoundGndPolygon, theCompositeMatrixInverse, theFocalLength, theLensDistortion, thePixelSize, thePrincipalPoint, ossimSensorModel::theRefImgPt, ossimRefPtr< T >::valid(), ossimDpt::x, and ossimDpt::y.

 {
 #if 0
    if((theBoundGndPolygon.getNumberOfVertices() > 0)&&
       (!theBoundGndPolygon.hasNans()))
    {
       if (!(theBoundGndPolygon.pointWithin(world_point)))
       {
          image_point.makeNan();
 //          image_point = extrapolate(world_point);
          return;
       }         
    }
 #endif
    ossimEcefPoint g_ecf(world_point);
    ossimEcefVector ecfRayDir(g_ecf - theAdjEcefPlatformPosition);
    ossimColumnVector3d camRayDir (theCompositeMatrixInverse*ecfRayDir.data());
    
       
    double scale = -theFocalLength/camRayDir[2];
    ossimDpt film (scale*camRayDir[0], scale*camRayDir[1]);
       
    if (theLensDistortion.valid())
    {
       ossimDpt filmOut;
       theLensDistortion->distort(film, filmOut);
       film = filmOut;
     }
    
      ossimDpt f1(film + thePrincipalPoint);
     ossimDpt p1(f1.x/thePixelSize.x,
                 -f1.y/thePixelSize.y);
 
     ossimDpt p0 (p1.x + theRefImgPt.x,
                  p1.y + theRefImgPt.y);
     
     image_point = p0;
 }

Member Data Documentation

◆ theAdjEcefPlatformPosition

ossimEcefPoint ossimApplanixUtmModel::theAdjEcefPlatformPosition

protected

Definition at line 78 of file ossimApplanixUtmModel.h.

Referenced by imagingRay(), updateModel(), and worldToLineSample().

◆ theBoreSightTx

double ossimApplanixUtmModel::theBoreSightTx

protected

Definition at line 64 of file ossimApplanixUtmModel.h.

Referenced by saveState().

◆ theBoreSightTy

double ossimApplanixUtmModel::theBoreSightTy

protected

Definition at line 65 of file ossimApplanixUtmModel.h.

Referenced by saveState().

◆ theBoreSightTz

double ossimApplanixUtmModel::theBoreSightTz

protected

Definition at line 66 of file ossimApplanixUtmModel.h.

Referenced by saveState().

◆ theCompositeMatrix

NEWMAT::Matrix ossimApplanixUtmModel::theCompositeMatrix

protected

uses file path to init model

Definition at line 59 of file ossimApplanixUtmModel.h.

Referenced by imagingRay(), ossimApplanixUtmModel(), and updateModel().

◆ theCompositeMatrixInverse

NEWMAT::Matrix ossimApplanixUtmModel::theCompositeMatrixInverse

protected

Definition at line 60 of file ossimApplanixUtmModel.h.

Referenced by ossimApplanixUtmModel(), updateModel(), and worldToLineSample().

◆ theEcefPlatformPosition

ossimEcefPoint ossimApplanixUtmModel::theEcefPlatformPosition

protected

Definition at line 70 of file ossimApplanixUtmModel.h.

Referenced by ossimApplanixUtmModel(), saveState(), and setPlatformPosition().

◆ theFocalLength

double ossimApplanixUtmModel::theFocalLength

protected

Definition at line 69 of file ossimApplanixUtmModel.h.

Referenced by imagingRay(), ossimApplanixUtmModel(), saveState(), setFocalLength(), and worldToLineSample().

◆ theKappa

double ossimApplanixUtmModel::theKappa

protected

Definition at line 63 of file ossimApplanixUtmModel.h.

Referenced by saveState(), and updateModel().

◆ theLensDistortion

ossimRefPtr<ossimMeanRadialLensDistortion> ossimApplanixUtmModel::theLensDistortion

protected

Definition at line 77 of file ossimApplanixUtmModel.h.

Referenced by imagingRay(), ossimApplanixUtmModel(), saveState(), and worldToLineSample().

◆ theOmega

double ossimApplanixUtmModel::theOmega

protected

Definition at line 61 of file ossimApplanixUtmModel.h.

Referenced by saveState(), and updateModel().

◆ thePhi

double ossimApplanixUtmModel::thePhi

protected

Definition at line 62 of file ossimApplanixUtmModel.h.

Referenced by saveState(), and updateModel().

◆ thePixelSize

ossimDpt ossimApplanixUtmModel::thePixelSize

protected

Definition at line 68 of file ossimApplanixUtmModel.h.

Referenced by imagingRay(), ossimApplanixUtmModel(), saveState(), setPixelSize(), and worldToLineSample().

◆ thePlatformPosition

ossimGpt ossimApplanixUtmModel::thePlatformPosition

protected

Definition at line 71 of file ossimApplanixUtmModel.h.

Referenced by saveState(), and updateModel().

◆ thePrincipalPoint

ossimDpt ossimApplanixUtmModel::thePrincipalPoint

protected

Definition at line 67 of file ossimApplanixUtmModel.h.

Referenced by imagingRay(), saveState(), setPrincipalPoint(), and worldToLineSample().

◆ theShiftValues

ossimEcefVector ossimApplanixUtmModel::theShiftValues

protected

Definition at line 75 of file ossimApplanixUtmModel.h.

Referenced by saveState().

◆ theUtmHemisphere

ossim_int8 ossimApplanixUtmModel::theUtmHemisphere

protected

Definition at line 73 of file ossimApplanixUtmModel.h.

Referenced by saveState(), and updateModel().

◆ theUtmPlatformPosition

ossimDpt3d ossimApplanixUtmModel::theUtmPlatformPosition

protected

Definition at line 74 of file ossimApplanixUtmModel.h.

Referenced by saveState().

◆ theUtmZone

ossim_int32 ossimApplanixUtmModel::theUtmZone

protected

Definition at line 72 of file ossimApplanixUtmModel.h.

Referenced by saveState(), and updateModel().

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

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ossimApplanixUtmModel() [1/2]

◆ ossimApplanixUtmModel() [2/2]

Member Function Documentation

◆ dup()

◆ imagingRay()

◆ initAdjustableParameters()

◆ lineSampleHeightToWorld()

◆ lineSampleToWorld()

◆ loadState()

◆ saveState()

◆ setFocalLength()

◆ setImageRect()

◆ setPixelSize()

◆ setPlatformPosition()

◆ setPrincipalPoint()

◆ setupOptimizer()

◆ updateModel()

◆ useForward()

◆ worldToLineSample()

Member Data Documentation

◆ theAdjEcefPlatformPosition

◆ theBoreSightTx

◆ theBoreSightTy

◆ theBoreSightTz

◆ theCompositeMatrix

◆ theCompositeMatrixInverse

◆ theEcefPlatformPosition

◆ theFocalLength

◆ theKappa

◆ theLensDistortion

◆ theOmega

◆ thePhi

◆ thePixelSize

◆ thePlatformPosition

◆ thePrincipalPoint

◆ theShiftValues

◆ theUtmHemisphere

◆ theUtmPlatformPosition

◆ theUtmZone