ossimObservationSet Class Reference

#include <ossimObservationSet.h>

Inheritance diagram for ossimObservationSet:

Public Member Functions
	ossimObservationSet ()
	~ossimObservationSet ()
ossim_uint32	numImages () const
ossim_uint32	numObs () const
ossim_uint32	numAdjPar () const
ossim_uint32	numMeas () const
ossimRefPtr< ossimPointObservation >	observ (const int &i) const
ossimImageGeometry *	getImageGeom (const int index)
void	setImageGeom (const int index, ossimImageGeometry *geom)
int	imIndex (const int &i) const
int	adjParCount (const int &i) const
ossimFilename	imageFile (const int index) const
bool	evaluate (NEWMAT::Matrix &measResiduals, NEWMAT::Matrix &objPartials, NEWMAT::Matrix &parPartials)
	standard evaluation [1] measResiduals: [x,y] residuals (numMeas X 2) [2] objPartials: object point partial derivatives (numMeas*3 X 2) [3] parPartials: image parameter partial derivatives (numParams X 2) More...
bool	addObservation (ossimRefPtr< ossimPointObservation > obs)
	operations More...
std::ostream &	print (std::ostream &os) const
	text output : header + tab separated tie points More...
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 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 Attributes
int	theNumAdjPar
int	theNumMeas
int	theNumPartials
std::vector< ossimRefPtr< ossimPointObservation > >	theObs
std::vector< int >	theImageIndex
std::vector< ossimFilename >	theImageFiles
std::vector< int >	theNumAdjParams
std::vector< ossimRefPtr< ossimImageHandler > >	theImageHandlers

Additional Inherited Members
Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Detailed Description

Definition at line 22 of file ossimObservationSet.h.

Constructor & Destructor Documentation

ossimObservationSet()

ossimObservationSet::ossimObservationSet

(

)

Definition at line 22 of file ossimObservationSet.cpp.

                                         :
theNumAdjPar(0),
theNumMeas(0),
theNumPartials(0)
{
   if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG)
      << "DEBUG: ossimObservationSet(): returning..." << std::endl;
}

~ossimObservationSet()

ossimObservationSet::~ossimObservationSet

(

)

Definition at line 32 of file ossimObservationSet.cpp.

References theImageHandlers.

{
   for (ossim_uint32 i=0; i<theImageHandlers.size(); ++i)
      theImageHandlers[i] = 0;
   // for (int i=0; i<theObs.size(); ++i)
   //    theObs[i] = 0;

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

Member Function Documentation

addObservation()

bool ossimObservationSet::addObservation

(

ossimRefPtr< ossimPointObservation >

obs

)

operations

Definition at line 45 of file ossimObservationSet.cpp.

{
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)
         <<"\n ossimObservationSet::addObservation: "<<obs->ID()<<endl;
   }

   theObs.push_back(obs);
   theNumMeas += obs->numMeas();


   // Update image list
   for (ossim_uint32 i=0; i<obs->numImages(); ++i)
   {
      if (traceDebug())
      {
         ossimNotify(ossimNotifyLevel_DEBUG)<<"  i="<<i<<endl;
      }

      bool found = false;

      // Check for image already in list
      for (ossim_uint32 j=0; j<theImageFiles.size(); ++j)
      {
         if (theImageFiles[j] == obs->imageFile(i))
         {
            found = true;
            theImageIndex.push_back(j);
            int nAdjPar = obs->numPars(i);
            theNumPartials += nAdjPar;
         }
      }

      // If not found yet, add to image list
      if (!found)
      {
         theImageFiles.push_back(obs->imageFile(i));

         // Geometry
         ossimRefPtr<ossimImageHandler> ih;
         ih = ossimImageHandlerRegistry::instance()->open(obs->imageFile(i));        
         theImageHandlers.push_back(ih);

         theImageIndex.push_back( (int)(theImageHandlers.size()-1) );

         // Handle ossimAdjustableParameterInterface
         ossimAdjustableParameterInterface* adjParIface = 
            obs->getImageGeom(i)->getAdjustableParameterInterface();
         if (!adjParIface)
         {
            return false;
         }
         else
         {
            if (traceDebug())
            {
               ossimNotify(ossimNotifyLevel_DEBUG)<<"  Found new image... "<<i<<endl;
            }
            int nAdjPar = adjParIface->getNumberOfAdjustableParameters();
            theNumAdjPar   += nAdjPar;
            theNumPartials += nAdjPar;
            theNumAdjParams.push_back(nAdjPar);
         }
      }
   }
   
   if (traceDebug())
   {
      ossimNotify(ossimNotifyLevel_DEBUG)<<"    theNumAdjPar   = "<<theNumAdjPar<<endl;
      ossimNotify(ossimNotifyLevel_DEBUG)<<"    theNumPartials = "<<theNumPartials<<endl;
   }

   return true;
}

adjParCount()

int ossimObservationSet::adjParCount ( const int &  i ) const

inline

Definition at line 45 of file ossimObservationSet.h.

{ return theNumAdjParams[i]; }

evaluate()

bool ossimObservationSet::evaluate	(	NEWMAT::Matrix &	measResiduals,
		NEWMAT::Matrix &	objPartials,
		NEWMAT::Matrix &	parPartials
	)

standard evaluation [1] measResiduals: [x,y] residuals (numMeas X 2) [2] objPartials: object point partial derivatives (numMeas*3 X 2) [3] parPartials: image parameter partial derivatives (numParams X 2)

Definition at line 134 of file ossimObservationSet.cpp.

References numMeas(), numObs(), theNumPartials, and theObs.

{

   // Dimension output matrices
   measResiduals = NEWMAT::Matrix(numMeas(), 2);
   objPartials   = NEWMAT::Matrix(numMeas()*3, 2);
   parPartials   = NEWMAT::Matrix(theNumPartials, 2);

   int img = 1;
   int cParIndex = 1;
   int cObjIndex = 1;
   for (ossim_uint32 cObs=0; cObs<numObs(); ++cObs)
   {
      int numMeasPerObs = theObs[cObs]->numMeas();
      if (traceDebug())
      {
         ossimNotify(ossimNotifyLevel_DEBUG)<<"\n cObs= "<<cObs;
      }

      for (int cImg=0; cImg<numMeasPerObs; ++cImg)
      {
         NEWMAT::Matrix cResid(1, 2);
         theObs[cObs]->getResiduals(cImg, cResid);
         if (traceDebug())
         {
            ossimNotify(ossimNotifyLevel_DEBUG)
               <<"\n   cImg, img, cObjIndex, cParIndex, cResid: "
               <<cImg<<" "<<img<<" "<<cObjIndex<<" "<<cParIndex<<" "<<cResid;
         }
         measResiduals.Row(img) = cResid;
         img++;

         NEWMAT::Matrix cObjPar(3, 2);
         theObs[cObs]->getObjSpacePartials(cImg, cObjPar);
         objPartials.SubMatrix(cObjIndex,cObjIndex+2,1,2) << cObjPar;
         cObjIndex += 3;

         int numPar = theObs[cObs]->numPars(cImg);
         NEWMAT::Matrix cParamPar(numPar, 2);
         theObs[cObs]->getParameterPartials(cImg, cParamPar);
         parPartials.SubMatrix(cParIndex,cParIndex+numPar-1,1,2) << cParamPar;
         cParIndex += numPar;
      }
   }

   return true;
}

getImageGeom()

ossimImageGeometry * ossimObservationSet::getImageGeom

(

const int

index

)

Definition at line 122 of file ossimObservationSet.cpp.

References theImageHandlers.

Referenced by ossimAutRegUtil::saveAdjustment(), and ossimAdjustmentExecutive::updateParameters().

{
   return theImageHandlers[index]->getImageGeometry().get();
}

imageFile()

ossimFilename ossimObservationSet::imageFile ( const int  index ) const

inline

Definition at line 48 of file ossimObservationSet.h.

Referenced by ossimAutRegUtil::saveAdjustment(), and ossimAdjustmentExecutive::updateParameters().

{ return theImageFiles[index]; }

imIndex()

int ossimObservationSet::imIndex ( const int &  i ) const

inline

Definition at line 42 of file ossimObservationSet.h.

Referenced by ossimAdjustmentExecutive::printResidualSummary().

{ return theImageIndex[i]; }

numAdjPar()

ossim_uint32 ossimObservationSet::numAdjPar ( ) const

inline

Definition at line 31 of file ossimObservationSet.h.

{ return (ossim_uint32)theNumAdjPar; }

numImages()

ossim_uint32 ossimObservationSet::numImages ( ) const

inline

Definition at line 29 of file ossimObservationSet.h.

Referenced by print(), and ossimAutRegUtil::saveAdjustment().

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

numMeas()

ossim_uint32 ossimObservationSet::numMeas ( ) const

inline

Definition at line 32 of file ossimObservationSet.h.

Referenced by evaluate().

{ return (ossim_uint32)theNumMeas; }

numObs()

ossim_uint32 ossimObservationSet::numObs ( ) const

inline

Definition at line 30 of file ossimObservationSet.h.

Referenced by evaluate(), and print().

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

observ()

ossimRefPtr<ossimPointObservation> ossimObservationSet::observ ( const int &  i ) const

inline

Definition at line 35 of file ossimObservationSet.h.

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

{ return theObs[i]; }

print()

std::ostream & ossimObservationSet::print ( std::ostream &  os ) const

virtual

text output : header + tab separated tie points

Reimplemented from ossimObject.

Definition at line 185 of file ossimObservationSet.cpp.

References ossim::isnan(), numImages(), numObs(), size, theImageFiles, theImageIndex, theNumAdjParams, and theObs.

{
   int idx = 0;

   os<<"\nImages:\n";
   ossim_uint32 spaces = 0;
   for (ossim_uint32 i=0; i<numImages(); ++i)
   {
      if (theImageFiles[i].size() > spaces)
         spaces = (ossim_uint32)theImageFiles[i].size();
   }
   for (ossim_uint32 i=0; i<numImages(); ++i)
   {
      int nAdj = theNumAdjParams[i];
      os << " "<<i+1<<": " << setw(spaces) << theImageFiles[i] << "   nPar: "<< setw(2) << nAdj <<endl;
   }

   os<<"\nObservations:";
   for (ossim_uint32 cObs=0; cObs<numObs(); ++cObs)
   {
      os << "\n " << theObs[cObs]->ID();
      if (ossim::isnan(theObs[cObs]->getScore()) == false)
      {
         os << std::setiosflags(std::ios::fixed) << std::setprecision(2);
         os << "\t" << theObs[cObs]->getScore();
      }
      else
      {
         os << "\t" << "nan";
      }
      os << "\t";
      theObs[cObs]->getGroundPoint().print(os);

      int numMeasPerObs = theObs[cObs]->numMeas();
      for (int cImg=0; cImg<numMeasPerObs; ++cImg)
      {
         os << "\n\t\t" << theImageIndex[idx]+1;
         os << "\t";
         theObs[cObs]->getMeasurement(cImg).print(os,1);
         idx++;
      }
   }
   return os;
}

setImageGeom()

void ossimObservationSet::setImageGeom	(	const int	index,
		ossimImageGeometry *	geom
	)

Definition at line 128 of file ossimObservationSet.cpp.

References theImageHandlers.

{
   theImageHandlers[index]->setImageGeometry(geom);
}

Member Data Documentation

theImageFiles

std::vector<ossimFilename> ossimObservationSet::theImageFiles

protected

Definition at line 82 of file ossimObservationSet.h.

Referenced by print().

theImageHandlers

std::vector< ossimRefPtr<ossimImageHandler> > ossimObservationSet::theImageHandlers

protected

Definition at line 87 of file ossimObservationSet.h.

Referenced by getImageGeom(), setImageGeom(), and ~ossimObservationSet().

theImageIndex

std::vector<int> ossimObservationSet::theImageIndex

protected

Definition at line 79 of file ossimObservationSet.h.

Referenced by print().

theNumAdjPar

int ossimObservationSet::theNumAdjPar

protected

Definition at line 71 of file ossimObservationSet.h.

theNumAdjParams

std::vector<int> ossimObservationSet::theNumAdjParams

protected

Definition at line 85 of file ossimObservationSet.h.

Referenced by print().

theNumMeas

int ossimObservationSet::theNumMeas

protected

Definition at line 72 of file ossimObservationSet.h.

theNumPartials

int ossimObservationSet::theNumPartials

protected

Definition at line 73 of file ossimObservationSet.h.

Referenced by evaluate().

theObs

std::vector< ossimRefPtr<ossimPointObservation> > ossimObservationSet::theObs

protected

Definition at line 76 of file ossimObservationSet.h.

Referenced by evaluate(), and print().

---

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

• ossimObservationSet.h
• ossimObservationSet.cpp