ossimAdjustmentExecutive Class Reference

#include <ossimAdjustmentExecutive.h>

Inheritance diagram for ossimAdjustmentExecutive:

Public Member Functions
	ossimAdjustmentExecutive (std::ostream &report)
	constructor More...
	ossimAdjustmentExecutive (ossimObservationSet &obsSet, std::ostream &report)
	ossimObservationSet constructor More...
	~ossimAdjustmentExecutive ()
	destructor More...
bool	initializeSolution (ossimObservationSet &obsSet)
	initialize adjustment solution More...
bool	runSolution ()
	run adjustment solution More...
void	summarizeSolution () const
	summarize solution More...
bool	computeResidualStatistics (NEWMAT::Matrix &res)
	compute mean and RMS error of residuals More...
double	computeSEUW ()
	compute SEUW More...
std::ostream &	printParameterCorrectionSummary (std::ostream &out) const
	Print parameter correction method. More...
std::ostream &	printObservationCorrectionSummary (std::ostream &out) const
	Print observation correction method. More...
std::ostream &	printResidualSummary (std::ostream &out) const
	Print residuals method. More...
bool	isValid () const
void	setMaxIter (const int maxIter)
void	setConvCriteria (const int convCriteria)
Public Member Functions inherited from ossimObject
	ossimObject ()
virtual	~ossimObject ()
virtual ossimObject *	dup () const
virtual ossimString	getShortName () const
virtual ossimString	getLongName () const
virtual ossimString	getDescription () const
virtual 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 bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const
virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)
virtual std::ostream &	print (std::ostream &out) const
	Generic print method. More...
virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const
virtual void	accept (ossimVisitor &visitor)
Public Member Functions inherited from ossimReferenced
	ossimReferenced ()
	ossimReferenced (const ossimReferenced &)
ossimReferenced &	operator= (const ossimReferenced &)
void	ref () const
	increment the reference count by one, indicating that this object has another pointer which is referencing it. More...
void	unref () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...
void	unref_nodelete () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...
int	referenceCount () const

Protected Member Functions
bool	updateParameters ()
	Update adjustable parameters. More...
bool	updateObservations ()
	Update observatin. More...
Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes
bool	theExecValid
ossimObservationSet *	theObsSet
ossimWLSBundleSolution *	theSol
ossimAdjSolutionAttributes *	theSolAttributes
double	theConvCriteria
int	theMaxIter
bool	theMaxIterExceeded
bool	theSolDiverged
bool	theSolConverged
int	theNumObsInSet
int	theNumImages
int	theNumParams
int	theNumMeasurements
int	theRankN
NEWMAT::Matrix	theMeasResiduals
NEWMAT::Matrix	theObjPartials
NEWMAT::Matrix	theParPartials
double	theXrms
double	theYrms
double	theXmean
double	theYmean
std::vector< double >	theSEUW
std::vector< double >	theParInitialValues
std::vector< double >	theParInitialStdDev
std::vector< ossimString >	theParDesc
std::vector< int >	theImgs
std::vector< double >	theObsInitialValues
std::vector< double >	theObsInitialStdDev
std::ostream &	theRep

Detailed Description

Definition at line 30 of file ossimAdjustmentExecutive.h.

Constructor & Destructor Documentation

ossimAdjustmentExecutive() [1/2]

ossimAdjustmentExecutive::ossimAdjustmentExecutive

(

std::ostream &

report

)

constructor

Definition at line 31 of file ossimAdjustmentExecutive.cpp.

   :
      theExecValid(false),
      theSol(0),
      theSolAttributes(0),
      theConvCriteria(5.0),
      theMaxIter(7),      
      theMaxIterExceeded(false),
      theSolDiverged(false),
      theSolConverged(false),
      theNumObsInSet(0),
      theNumImages(0),
      theNumParams(0),
      theNumMeasurements(0),
      theRankN(0),
      theMeasResiduals(),
      theObjPartials(),
      theParPartials(),
      theXrms(0.0),
      theYrms(0.0),
      theXmean(0.0),
      theYmean(0.0),
      theSEUW(0),
      theParInitialValues(0),
      theParInitialStdDev(0),
      theParDesc(0),
      theObsInitialValues(0),
      theObsInitialStdDev(0),
      theRep(report)      
{
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "\nossimAdjustmentExecutive::ossimAdjustmentExecutive 1 DEBUG:" << std::endl;
   }
}

ossimAdjustmentExecutive() [2/2]

ossimAdjustmentExecutive::ossimAdjustmentExecutive	(	ossimObservationSet &	obsSet,
		std::ostream &	report
	)

ossimObservationSet constructor

Parameters

obsSet

tiepoint observation set.

Definition at line 76 of file ossimAdjustmentExecutive.cpp.

   :
      theExecValid(false),
      theSol(0),
      theSolAttributes(0),
      theConvCriteria(5.0),
      theMaxIter(7),      
      theMaxIterExceeded(false),
      theSolDiverged(false),
      theSolConverged(false),
      theNumObsInSet(0),
      theNumImages(0),
      theNumParams(0),
      theNumMeasurements(0),
      theRankN(0),
      theMeasResiduals(),
      theObjPartials(),
      theParPartials(),
      theXrms(0.0),
      theYrms(0.0),
      theXmean(0.0),
      theYmean(0.0),
      theSEUW(0),
      theParInitialValues(0),
      theParInitialStdDev(0),
      theParDesc(0),
      theObsInitialValues(0),
      theObsInitialStdDev(0),
      theRep(report)      
{
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "\nossimAdjustmentExecutive::ossimAdjustmentExecutive 2 DEBUG:" << std::endl;
   }

   theExecValid = initializeSolution(obsSet);
}

~ossimAdjustmentExecutive()

ossimAdjustmentExecutive::~ossimAdjustmentExecutive

(

)

destructor

Definition at line 120 of file ossimAdjustmentExecutive.cpp.

References theSol, and theSolAttributes.

{
   delete theSolAttributes;
   delete theSol;

   if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG)
      << "DEBUG: ~ossimAdjustmentExecutive(): returning..." << std::endl;
}

Member Function Documentation

computeResidualStatistics()

bool ossimAdjustmentExecutive::computeResidualStatistics

(

NEWMAT::Matrix &

res

)

compute mean and RMS error of residuals

Parameters

res	image residual matrix.

Returns

true on success, false on error.

Definition at line 592 of file ossimAdjustmentExecutive.cpp.

References n, theNumMeasurements, theRep, theXmean, theXrms, theYmean, and theYrms.

{
   theXmean = 0.0;
   theYmean = 0.0;
   theXrms  = 0.0;
   theYrms  = 0.0;
   for (int n=1; n<=res.Nrows();++n)
   {
      theXmean += res(n,1);
      theYmean += res(n,2);
      theXrms  += res(n,1)*res(n,1);
      theYrms  += res(n,2)*res(n,2);
   }

   theXmean /= theNumMeasurements;
   theYmean /= theNumMeasurements;
   theXrms   = sqrt(theXrms/theNumMeasurements);
   theYrms   = sqrt(theYrms/theNumMeasurements);

   theRep<<"\n";
   theRep<<" ______________Mean:";
   theRep<<setprecision(1)<<setw(8)<<theXmean;
   theRep<<setprecision(1)<<setw(8)<<theYmean;
   theRep<<"    RMS:";
   theRep<<setprecision(1)<<setw(8)<<theXrms;
   theRep<<setprecision(1)<<setw(8)<<theYrms;
   theRep<<endl;

   return true;
}

computeSEUW()

double ossimAdjustmentExecutive::computeSEUW

(

)

compute SEUW

Returns

standard error of unit weight.

Definition at line 630 of file ossimAdjustmentExecutive.cpp.

References ossimObservationSet::observ(), ossimAdjSolutionAttributes::theImagePtCov, ossimAdjSolutionAttributes::theMeasResiduals, theNumMeasurements, theNumObsInSet, theNumParams, ossimAdjSolutionAttributes::theObjectPtCov, theObsSet, theParInitialStdDev, theRankN, theSolAttributes, and ossimAdjSolutionAttributes::theTotalCorrections.

{
   double vuw = 0.0;

   // Observation contributions
   int cMeas = 1;
   for (int obs=0; obs<theNumObsInSet; ++obs)
   {
      NEWMAT::ColumnVector tc(3);
      NEWMAT::Matrix cov(3,3);
      int idx = theNumParams + 3*obs + 1;
      tc  = theSolAttributes->theTotalCorrections.Rows(idx, idx+2);
      idx -= theNumParams;
      cov = theSolAttributes->theObjectPtCov.SubMatrix(idx,idx+2,1,3);
      for (int i=1; i<=3; ++i)
      {
         vuw += tc(i)*tc(i)/cov(i,i);
      }
      for (ossim_uint32 meas=0; meas<theObsSet->observ(obs)->numMeas(); ++meas)
      {
         NEWMAT::RowVector res(2);
         NEWMAT::Matrix cov(2,2);
         res = theSolAttributes->theMeasResiduals.Row(cMeas);
         int start = (cMeas-1)*2 + 1;
         cov = theSolAttributes->theImagePtCov.SubMatrix(start,start+1,1,2);
         for (int i=1; i<=2; ++i)
         {
            vuw += res(i)*res(i)/cov(i,i);
         }
         cMeas++;
      }
   }

   // Parameter contributions
   NEWMAT::ColumnVector tc(theNumParams);
   tc = theSolAttributes->theTotalCorrections.Rows(1, theNumParams);
   for (int par=0; par<theNumParams; ++par)
   {
      vuw += tc(par+1)*tc(par+1)/(theParInitialStdDev[par]*theParInitialStdDev[par]);
   }

   // DF
   double df = theNumMeasurements*2 - theRankN;
   if (df<=0.0)
      df = 1.0;

   // SEUW
   double seuw = sqrt(vuw/df);

   return seuw;
}

initializeSolution()

bool ossimAdjustmentExecutive::initializeSolution

(

ossimObservationSet &

obsSet

)

initialize adjustment solution

Parameters

obsSet

tiepoint observation set.

Returns

true on success, false on error.

Definition at line 136 of file ossimAdjustmentExecutive.cpp.

{
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "\nossimAdjustmentExecutive::initializeSolution DEBUG:" << std::endl; 
   }
   ossimString ts;
   ossim::getFormattedTime("%a %m.%d.%y %H:%M:%S", false, ts);

   theObsSet = &obsSet;

   // Initial report output
   theRep << "\nossimAdjustmentExecutive Report     ";
   theRep << ts;
   theRep << "\n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~";
   theRep << endl;
   theObsSet->print(theRep);
   theRep << endl;

   theExecValid = false;

   // Adjustment traits
   theNumObsInSet     = theObsSet->numObs();
   if (theNumObsInSet == 0)
      return theExecValid;

   theNumImages       = theObsSet->numImages();
   theNumParams       = theObsSet->numAdjPar();
   theNumMeasurements = theObsSet->numMeas();

   theRankN = theNumParams + theNumObsInSet*3;


   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         <<"\n theNumObsInSet     = "<<theNumObsInSet
         <<"\n theNumImages       = "<<theNumImages
         <<"\n theNumParams       = "<<theNumParams
         <<"\n theNumMeasurements = "<<theNumMeasurements
         <<"\n theRankN           = "<<theRankN
         <<std::endl;
   }


   // Instantiate solution attributes
   theSolAttributes = new ossimAdjSolutionAttributes
      (theNumObsInSet, theNumImages, theNumMeasurements, theRankN);


   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         <<"\n\nInitial Parameter Setup....";

      for (int i=0; i<theNumImages; i++)
      {
         ossimAdjustableParameterInterface* iface =
            theObsSet->getImageGeom(i)->getAdjustableParameterInterface();
         int np = iface->getNumberOfAdjustableParameters();
         for (int cp=0; cp<np; ++cp)
         {
            ossimNotify(ossimNotifyLevel_DEBUG)
               <<"\n "<<cp<<"  "<<iface->getParameterDescription(cp)
               // <<"=  "<<iface->getAdjustableParameter(cp)
               <<"=  "<<iface->getParameterCenter(cp)
               <<", units= "<<iface->getParameterUnit(cp)
               <<", sigma= "<<iface->getParameterSigma(cp);
         }
      }
      ossimNotify(ossimNotifyLevel_DEBUG)<<std::endl;
   }

   // Save parameter initial values and variances
   int start = 1;
   theSolAttributes->theAdjParCov.ReSize(theNumParams,theNumParams);
   for (int i=0; i<theNumImages; i++)
   {
      ossimAdjustableParameterInterface* iface =
         theObsSet->getImageGeom(i)->getAdjustableParameterInterface();
      int np = iface->getNumberOfAdjustableParameters();
      NEWMAT::Matrix parCov(np,np);
      parCov = 0.0;
      for (int cp=0; cp<np; ++cp)
      {
         // Get the a priori value
         // theParInitialValues.push_back(iface->getAdjustableParameter(cp));
         theParInitialValues.push_back(iface->getParameterCenter(cp));
         theImgs.push_back(i);

         theParDesc.push_back(iface->getParameterDescription(cp));
         double sig = iface->getParameterSigma(cp);
         
         // If parameter is locked, tighten down the sigma
         // TODO: Eventually need better handling of this
         if (iface->isParameterLocked(cp))
            sig /= 1000;

         theParInitialStdDev.push_back(sig);
         parCov(cp+1,cp+1) = sig*sig;
      }
      theSolAttributes->
         theAdjParCov.SubMatrix(start,start+np-1,start,start+np-1) = parCov;
      start += np;
   }

   // Ensure initial estimates for observations
   //   TODO: Currently uses mean of single-ray point drops; should
   //         eventually use multi-ray intersection, which is only
   //         available via the ossimSensorModelTuple class.
   for (int obs=0; obs<theNumObsInSet; obs++)
   {
      // If ground position not set, initialize it
      if (theObsSet->observ(obs)->hasNans())
      {
         double latMean = 0.0;
         double lonMean = 0.0;
         double hgtMean = 0.0;
         ossimGpt gpt;
         ossimDpt ipt;
         for (ossim_uint32 meas=0; meas<theObsSet->observ(obs)->numMeas(); ++meas)
         {
            ipt = theObsSet->observ(obs)->getMeasurement(meas);
            theObsSet->observ(obs)->getImageGeom(meas)->localToWorld(ipt, gpt);
            latMean += gpt.latr();
            lonMean += gpt.lonr();
            hgtMean += gpt.height();
         }
         latMean /= theObsSet->observ(obs)->numMeas();
         lonMean /= theObsSet->observ(obs)->numMeas();
         hgtMean /= theObsSet->observ(obs)->numMeas();
         theObsSet->observ(obs)->Gpt().latr(latMean);
         theObsSet->observ(obs)->Gpt().lonr(lonMean);
         theObsSet->observ(obs)->Gpt().height(hgtMean);
      }
   }


   // Save observation intial values
   for (int obs=0; obs<theNumObsInSet; obs++)
   {
      theObsInitialValues.push_back(theObsSet->observ(obs)->Gpt().latr());
      theObsInitialValues.push_back(theObsSet->observ(obs)->Gpt().lonr());
      theObsInitialValues.push_back(theObsSet->observ(obs)->Gpt().height());
   }


   // Load obj/image covariance data
   int cMeas = 1;
   theSolAttributes->theImagePtCov.ReSize(theNumMeasurements*2,2);
   theSolAttributes->theObjectPtCov.ReSize(theNumObsInSet*3,3);

   for (int obs=0; obs<theNumObsInSet; ++obs)
   {
      NEWMAT::Matrix ocov;
      theObsSet->observ(obs)->getObsCov(ocov);
      int idx = obs*3+1;
      theSolAttributes->theObjectPtCov.Rows(idx,idx+2) = ocov;
      for (int i=1; i<=3; ++i)
         theObsInitialStdDev.push_back(sqrt(ocov(i,i)));
      for (ossim_uint32 meas=0; meas<theObsSet->observ(obs)->numMeas(); ++meas)
      {
         NEWMAT::Matrix mcov;
         theObsSet->observ(obs)->getMeasCov((int)meas, mcov);
         int start = (cMeas-1)*2 + 1;
         theSolAttributes->theImagePtCov.Rows(start,start+1) = mcov;
         cMeas++;   
      }
   }


   // Load obj/image xref map
   for (int obs=0; obs<theNumObsInSet; ++obs)
   {
      for (ossim_uint32 meas=0; meas<theObsSet->observ(obs)->numMeas(); ++meas)
      {
         int img = theObsSet->imIndex(meas);
         theSolAttributes->theObjImgXref.insert(pair<int, int>(obs, img));
      }
   }

   // Load image/numpar xref map
   for (int img=0; img<theNumImages; ++img)
   {
      int numpar = theObsSet->adjParCount(img);
      theSolAttributes->theImgNumparXref.insert(pair<int, int>(img, numpar));
   }

   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         <<"\ntheObjImgXref multimap  Obs/Img ....\n";
      for (ObjImgMapIter_t it =theSolAttributes->theObjImgXref.begin();
                           it!=theSolAttributes->theObjImgXref.end(); ++it)
      {
         ossimNotify(ossimNotifyLevel_DEBUG)
            <<"  "<<it->first<<"    "<<it->second<<endl;
      }
      ossimNotify(ossimNotifyLevel_DEBUG)
         <<"\ntheImgNumparXref map  Img/Numpar ....\n";
      for (ImgNumparMapIter_t it =theSolAttributes->theImgNumparXref.begin();
                              it!=theSolAttributes->theImgNumparXref.end(); ++it)
      {
         ossimNotify(ossimNotifyLevel_DEBUG)
            <<"  "<<it->first<<"    "<<it->second<<endl;
      }
   }

 
   theRep << "\n Iteration 0...";

   updateParameters();

   // Perform initial (0th iteration) observation evaluation
   theObsSet->evaluate(theMeasResiduals, theObjPartials, theParPartials);

   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         <<"\n theObjPartials\n"
         <<theObjPartials;
   }

   // Report residual summary
   printResidualSummary(theRep);

   // Residual statistics
   computeResidualStatistics(theMeasResiduals);

   // Load partials
   theSolAttributes->theObjPartials = theObjPartials;
   theSolAttributes->theParPartials = theParPartials;

   // Load residuals
   theSolAttributes->theMeasResiduals = theMeasResiduals;

   // Initial standard error of unit weight
   theSEUW.push_back(computeSEUW());

   // Instantiate solution object
   theSol = new ossimWLSBundleSolution();

   theExecValid = true;

   return theExecValid;
}

isValid()

bool ossimAdjustmentExecutive::isValid ( ) const

inline

Definition at line 104 of file ossimAdjustmentExecutive.h.

{ return theExecValid; }

printObservationCorrectionSummary()

std::ostream & ossimAdjustmentExecutive::printObservationCorrectionSummary

(

std::ostream &

out

)

const

Print observation correction method.

Definition at line 716 of file ossimAdjustmentExecutive.cpp.

References DEG_PER_RAD, ossimObservationSet::observ(), ossimAdjSolutionAttributes::theFullCovMatrix, ossimAdjSolutionAttributes::theLastCorrections, theNumObsInSet, theNumParams, theObsInitialStdDev, theObsInitialValues, theObsSet, theSolAttributes, and ossimAdjSolutionAttributes::theTotalCorrections.

{
   out<<"\nObservation Corrections...";
   out<<"\n  n         observation    a_priori  total_corr   last_corr initial_std    prop_std";
   for (int obs=0; obs<theNumObsInSet; ++obs)
   {
      out<<"\n "<<setprecision(5)<<setw(2)<<obs+1;
      out<<" "<<setw(19)<<theObsSet->observ(obs)->ID();
      double mPerRadp = theObsSet->observ(obs)->Gpt().metersPerDegree().y*DEG_PER_RAD;
      double mPerRadl = theObsSet->observ(obs)->Gpt().metersPerDegree().x*DEG_PER_RAD;
      // ossim_uint32 iobs = obs*3+1;

      for (int k=0; k<3; ++k)
      {
         int idx = theNumParams + obs*3 + k + 1;
         if (k<2)
            out<<setw(12)<<theObsInitialValues[obs*3+k]*DEG_PER_RAD;
         else
            out<<setw(12)<<theObsInitialValues[obs*3+k];
         double factor;
         if (k==0)
            factor = mPerRadp;
         else if (k==1)
            factor = mPerRadl;
         else
            factor = 1.0;
         out<<setw(12)<<theSolAttributes->theTotalCorrections(idx)*factor;
         out<<setw(12)<<theSolAttributes->theLastCorrections(idx)*factor;
         out<<setw(12)<<theObsInitialStdDev[obs*3+k]*factor;
         out<<setw(12)<<sqrt(theSolAttributes->theFullCovMatrix(idx,idx))*factor;
         out<<endl<<"                       ";
      }
   }
   out<<endl;

   return out;
}

printParameterCorrectionSummary()

std::ostream & ossimAdjustmentExecutive::printParameterCorrectionSummary

(

std::ostream &

out

)

const

Print parameter correction method.

Definition at line 690 of file ossimAdjustmentExecutive.cpp.

References ossimAdjSolutionAttributes::theFullCovMatrix, theImgs, ossimAdjSolutionAttributes::theLastCorrections, theNumParams, theParDesc, theParInitialStdDev, theParInitialValues, theSolAttributes, and ossimAdjSolutionAttributes::theTotalCorrections.

{
   out<<"\nParameter Corrections...";
   out<<"\n  n im        parameter    a_priori  total_corr   last_corr initial_std    prop_std";
   for (int pc=1; pc<=theNumParams; ++pc)
   {
      out<<"\n "<<setprecision(5)<<setw(2)<<pc;
      out<<setw(3)<<theImgs[pc-1]+1;
      out<<setw(17)<<theParDesc[pc-1];
      out<<setw(12)<<theParInitialValues[pc-1];
      out<<setw(12)<<theSolAttributes->theTotalCorrections(pc);
      out<<setw(12)<<theSolAttributes->theLastCorrections(pc);
      out<<setw(12)<<theParInitialStdDev[pc-1];
      out<<setw(12)<<sqrt(theSolAttributes->theFullCovMatrix(pc,pc));
   }
   out<<endl;

   return out;
}

printResidualSummary()

std::ostream & ossimAdjustmentExecutive::printResidualSummary

(

std::ostream &

out

)

const

Print residuals method.

Definition at line 760 of file ossimAdjustmentExecutive.cpp.

References ossimObservationSet::imIndex(), ossimObservationSet::observ(), theMeasResiduals, theNumObsInSet, and theObsSet.

{
   out<<"\nMeasurement Residuals...";
   out<<"\n observation   image    samp    line    initial_meas";
   int j = 0;
   for (int obs=0; obs<theNumObsInSet; obs++)
   {
      int numMeasPerObs = theObsSet->observ(obs)->numMeas();
      for (int meas=0; meas<numMeasPerObs; ++meas)
      {
         int imIdx = theObsSet->imIndex(j);
         ++j;
         out<<"\n";
         out<<setw(12)<<theObsSet->observ(obs)->ID();
         out<<setw(8)<<imIdx+1;
         out<<setprecision(1)<<setw(8)<<theMeasResiduals(j,1);
         out<<setprecision(1)<<setw(8)<<theMeasResiduals(j,2);
         out<<"    ";
         theObsSet->observ(obs)->getMeasurement(meas).print(out,1);
      }
      out<<endl;
   }

   return out;
}

runSolution()

bool ossimAdjustmentExecutive::runSolution

(

)

run adjustment solution

Returns

true on success, false on error.

Definition at line 391 of file ossimAdjustmentExecutive.cpp.

{
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         << "\nossimAdjustmentExecutive::runSolution DEBUG:" << std::endl;
   }


   // Iterative loop
   int iter = 0;

   while (iter < theMaxIter &&
          !theSolConverged &&
          !theSolDiverged &&
          theExecValid)
   {
      iter++;

      theRep << "\n Iteration "<<iter<<"...";

      // Execute solution
      theExecValid = theSol->run(theSolAttributes);


      if (theExecValid)
      {
         // Report corrections
         printParameterCorrectionSummary(theRep);
         printObservationCorrectionSummary(theRep);

         // Update adjustable parameters
         updateParameters();

         // Update ground points
         updateObservations();

         // Perform observation evaluation
         theObsSet->evaluate(theMeasResiduals, theObjPartials, theParPartials);

         // Load partials
         theSolAttributes->theObjPartials = theObjPartials;
         theSolAttributes->theParPartials = theParPartials;

         // Load residuals
         theSolAttributes->theMeasResiduals = theMeasResiduals;

         // Report residual summary
         printResidualSummary(theRep);

         // Residual statistics
         computeResidualStatistics(theMeasResiduals);

         // Compute SEUW for current iteration
         theSEUW.push_back(computeSEUW());

         // Check convergence
         double percChange = 
         fabs((theSEUW[iter]-theSEUW[iter-1])/theSEUW[iter-1])*100.0;

         if ((percChange <= theConvCriteria) && iter>1)
            theSolConverged = true;
         else if (iter == theMaxIter)
            theMaxIterExceeded = true;
         else if (iter >= 3)
         {
            if (theSEUW[iter]   > theSEUW[iter-1] &&
                theSEUW[iter-1] > theSEUW[iter-2] &&
                theSEUW[iter-2] > theSEUW[iter-3])
               theSolDiverged = true;
         }
         else
            theSolConverged = false;
      }
   }

   return theExecValid;
}

setConvCriteria()

void ossimAdjustmentExecutive::setConvCriteria ( const int  convCriteria )

inline

Definition at line 108 of file ossimAdjustmentExecutive.h.

{ theConvCriteria = convCriteria; }

setMaxIter()

void ossimAdjustmentExecutive::setMaxIter ( const int  maxIter )

inline

Definition at line 106 of file ossimAdjustmentExecutive.h.

{ theMaxIter = maxIter; }

summarizeSolution()

void ossimAdjustmentExecutive::summarizeSolution

(

)

const

summarize solution

Definition at line 556 of file ossimAdjustmentExecutive.cpp.

References ossim::getFormattedTime(), theConvCriteria, theExecValid, theMaxIter, theMaxIterExceeded, theNumImages, theNumMeasurements, theNumObsInSet, theNumParams, theRep, theSEUW, theSolConverged, and theSolDiverged.

{
   theRep << "\nossimAdjustmentExecutive Summary..."<<std::endl;
   theRep << " Valid Exec:         "<<(theExecValid?"true":"false")<<std::endl;
   theRep << " Nbr Ground Pts:     "<<theNumObsInSet<<std::endl;
   theRep << " Nbr Image Points:   "<<theNumMeasurements<<std::endl;
   theRep << " Nbr Images:         "<<theNumImages<<std::endl;
   theRep << " Nbr Parameters:     "<<theNumParams<<std::endl;
   theRep << " -------------------------"<<std::endl;
   theRep << " Solution Converged: "<<(theSolConverged?"true":"false")<<std::endl;
   theRep << " Solution Diverged:  "<<(theSolDiverged?"true":"false")<<std::endl;
   theRep << " Max Iter Exceeded:  "<<(theMaxIterExceeded?"true":"false")<<std::endl;
   theRep << " Max Iterations:     "<<theMaxIter<<std::endl;
   theRep << " Convergence Crit:   "<<theConvCriteria<<"%"<<std::endl;

   // SEUW history
   theRep << "\n SEUW Trace...";
   theRep << "\n   Iter        SEUW";
   for (ossim_uint32 iter=0; iter<theSEUW.size(); ++iter)
      theRep << "\n"<<setw(7)<<iter<<setw(12)<<setprecision(3)<<theSEUW[iter];

   ossimString ts;
   ossim::getFormattedTime("%a %m.%d.%y %H:%M:%S", false, ts);
   theRep << "\n";
   theRep << "\n" << ts;
   theRep << "\n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~";
   theRep << endl;
}

updateObservations()

bool ossimAdjustmentExecutive::updateObservations ( )

protected

Update observatin.

Definition at line 527 of file ossimAdjustmentExecutive.cpp.

References ossimObservationSet::observ(), ossimAdjSolutionAttributes::theLastCorrections, theNumObsInSet, theNumParams, theObsSet, and theSolAttributes.

{
   int currPar = theNumParams + 1;
   for (int obs=0; obs<theNumObsInSet; ++obs)
   {
      double lat = theObsSet->observ(obs)->Gpt().latr();
      lat += theSolAttributes->theLastCorrections(currPar);
      currPar++;    
      double lon = theObsSet->observ(obs)->Gpt().lonr();
      lon += theSolAttributes->theLastCorrections(currPar);
      currPar++;    
      double hgt = theObsSet->observ(obs)->Gpt().height();
      hgt += theSolAttributes->theLastCorrections(currPar);
      currPar++;    
      theObsSet->observ(obs)->Gpt().latr(lat);
      theObsSet->observ(obs)->Gpt().lonr(lon);
      theObsSet->observ(obs)->Gpt().height(hgt);
   }

   return true;
}

updateParameters()

bool ossimAdjustmentExecutive::updateParameters ( )

protected

Update adjustable parameters.

Definition at line 477 of file ossimAdjustmentExecutive.cpp.

References ossimImageGeometry::getAdjustableParameterInterface(), ossimObservationSet::getImageGeom(), ossimAdjustableParameterInterface::getParameterCenter(), ossimObservationSet::imageFile(), ossimObservationSet::observ(), ossimAdjustableParameterInterface::setParameterCenter(), ossimAdjSolutionAttributes::theImgNumparXref, ossimAdjSolutionAttributes::theLastCorrections, theNumImages, theNumObsInSet, theObsSet, and theSolAttributes.

{
   // Update local geometries
   int currPar = 1;
   for (int img=0; img<theNumImages; ++img)
   {
      ossimAdjustableParameterInterface* iface =
         theObsSet->getImageGeom(img)->getAdjustableParameterInterface();

      for (int par=0; par<theSolAttributes->theImgNumparXref[img]; ++par)
      {
         // double middle = iface->getAdjustableParameter(par);
         // middle += theSolAttributes->theLastCorrections(currPar);
         // iface->setAdjustableParameter(par, middle, true);

         double middle = iface->getParameterCenter(par);
         middle += theSolAttributes->theLastCorrections(currPar);
         iface->setParameterCenter(par, middle, true);

         currPar++;    
      }
   }

   // Copy updated local geometries to observation geometries
   for (int img=0; img<theNumImages; ++img)
   {
      for (int obs=0; obs<theNumObsInSet; ++obs)
      {
         for (ossim_uint32 imgInObs=0; imgInObs<theObsSet->observ(obs)->numImages(); ++imgInObs)
         {
            if (theObsSet->observ(obs)->imageFile(imgInObs) ==
                theObsSet->imageFile(img))
            {
               theObsSet->observ(obs)->setImageGeom
                  (imgInObs, theObsSet->getImageGeom(img));
            }
         }
      }
   }

   return true;
}

Member Data Documentation

theConvCriteria

double ossimAdjustmentExecutive::theConvCriteria

protected

Definition at line 123 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution().

theExecValid

bool ossimAdjustmentExecutive::theExecValid

protected

Definition at line 111 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution().

theImgs

std::vector<int> ossimAdjustmentExecutive::theImgs

protected

Definition at line 152 of file ossimAdjustmentExecutive.h.

Referenced by printParameterCorrectionSummary().

theMaxIter

int ossimAdjustmentExecutive::theMaxIter

protected

Definition at line 124 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution().

theMaxIterExceeded

bool ossimAdjustmentExecutive::theMaxIterExceeded

protected

Definition at line 125 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution().

theMeasResiduals

NEWMAT::Matrix ossimAdjustmentExecutive::theMeasResiduals

protected

Definition at line 137 of file ossimAdjustmentExecutive.h.

Referenced by printResidualSummary().

theNumImages

int ossimAdjustmentExecutive::theNumImages

protected

Definition at line 131 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution(), and updateParameters().

theNumMeasurements

int ossimAdjustmentExecutive::theNumMeasurements

protected

Definition at line 133 of file ossimAdjustmentExecutive.h.

Referenced by computeResidualStatistics(), computeSEUW(), and summarizeSolution().

theNumObsInSet

int ossimAdjustmentExecutive::theNumObsInSet

protected

Definition at line 130 of file ossimAdjustmentExecutive.h.

Referenced by computeSEUW(), printObservationCorrectionSummary(), printResidualSummary(), summarizeSolution(), updateObservations(), and updateParameters().

theNumParams

int ossimAdjustmentExecutive::theNumParams

protected

Definition at line 132 of file ossimAdjustmentExecutive.h.

Referenced by computeSEUW(), printObservationCorrectionSummary(), printParameterCorrectionSummary(), summarizeSolution(), and updateObservations().

theObjPartials

NEWMAT::Matrix ossimAdjustmentExecutive::theObjPartials

protected

Definition at line 138 of file ossimAdjustmentExecutive.h.

theObsInitialStdDev

std::vector<double> ossimAdjustmentExecutive::theObsInitialStdDev

protected

Definition at line 156 of file ossimAdjustmentExecutive.h.

Referenced by printObservationCorrectionSummary().

theObsInitialValues

std::vector<double> ossimAdjustmentExecutive::theObsInitialValues

protected

Definition at line 155 of file ossimAdjustmentExecutive.h.

Referenced by printObservationCorrectionSummary().

theObsSet

ossimObservationSet* ossimAdjustmentExecutive::theObsSet

protected

Definition at line 114 of file ossimAdjustmentExecutive.h.

Referenced by computeSEUW(), printObservationCorrectionSummary(), printResidualSummary(), updateObservations(), and updateParameters().

theParDesc

std::vector<ossimString> ossimAdjustmentExecutive::theParDesc

protected

Definition at line 151 of file ossimAdjustmentExecutive.h.

Referenced by printParameterCorrectionSummary().

theParInitialStdDev

std::vector<double> ossimAdjustmentExecutive::theParInitialStdDev

protected

Definition at line 150 of file ossimAdjustmentExecutive.h.

Referenced by computeSEUW(), and printParameterCorrectionSummary().

theParInitialValues

std::vector<double> ossimAdjustmentExecutive::theParInitialValues

protected

Definition at line 149 of file ossimAdjustmentExecutive.h.

Referenced by printParameterCorrectionSummary().

theParPartials

NEWMAT::Matrix ossimAdjustmentExecutive::theParPartials

protected

Definition at line 139 of file ossimAdjustmentExecutive.h.

theRankN

int ossimAdjustmentExecutive::theRankN

protected

Definition at line 134 of file ossimAdjustmentExecutive.h.

Referenced by computeSEUW().

theRep

std::ostream& ossimAdjustmentExecutive::theRep

protected

Definition at line 158 of file ossimAdjustmentExecutive.h.

Referenced by computeResidualStatistics(), and summarizeSolution().

theSEUW

std::vector<double> ossimAdjustmentExecutive::theSEUW

protected

Definition at line 146 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution().

theSol

ossimWLSBundleSolution* ossimAdjustmentExecutive::theSol

protected

Definition at line 117 of file ossimAdjustmentExecutive.h.

Referenced by ~ossimAdjustmentExecutive().

theSolAttributes

ossimAdjSolutionAttributes* ossimAdjustmentExecutive::theSolAttributes

protected

Definition at line 120 of file ossimAdjustmentExecutive.h.

Referenced by computeSEUW(), printObservationCorrectionSummary(), printParameterCorrectionSummary(), updateObservations(), updateParameters(), and ~ossimAdjustmentExecutive().

theSolConverged

bool ossimAdjustmentExecutive::theSolConverged

protected

Definition at line 127 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution().

theSolDiverged

bool ossimAdjustmentExecutive::theSolDiverged

protected

Definition at line 126 of file ossimAdjustmentExecutive.h.

Referenced by summarizeSolution().

theXmean

double ossimAdjustmentExecutive::theXmean

protected

Definition at line 144 of file ossimAdjustmentExecutive.h.

Referenced by computeResidualStatistics().

theXrms

double ossimAdjustmentExecutive::theXrms

protected

Definition at line 142 of file ossimAdjustmentExecutive.h.

Referenced by computeResidualStatistics().

theYmean

double ossimAdjustmentExecutive::theYmean

protected

Definition at line 145 of file ossimAdjustmentExecutive.h.

Referenced by computeResidualStatistics().

theYrms

double ossimAdjustmentExecutive::theYrms

protected

Definition at line 143 of file ossimAdjustmentExecutive.h.

Referenced by computeResidualStatistics().

---

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

• ossimAdjustmentExecutive.h
• ossimAdjustmentExecutive.cpp