ossimFormosatModel.h

Go to the documentation of this file.

//----------------------------------------------------------------------------
//
// "Copyright Centre National d'Etudes Spatiales"
//
// License:  LGPL
//
// See LICENSE.txt file in the top level directory for more details.
//
//----------------------------------------------------------------------------
// $Id$



#ifndef ossimFormosatModel_HEADER
#define ossimFormosatModel_HEADER

#include <iostream>
using namespace std;

#include <ossim/plugin/ossimPluginConstants.h>
#include <ossim/projection/ossimSensorModel.h>
#include <ossim/base/ossimIpt.h>
#include <ossim/base/ossimFilename.h>
#include <ossim/base/ossimGpt.h>
#include <ossim/base/ossimDpt.h>
#include <ossim/base/ossimEcefRay.h>
#include <ossim/base/ossimEcefPoint.h>
#include <ossim/base/ossimMatrix3x3.h>

class ossimFormosatDimapSupportData;

namespace ossimplugins
{
class OSSIM_PLUGINS_DLL ossimFormosatModel : public ossimSensorModel
{
public:
   ossimFormosatModel();
   ossimFormosatModel(ossimFormosatDimapSupportData* sd);
   ossimFormosatModel(const ossimFilename& init_file);
   ossimFormosatModel(const ossimKeywordlist& geom_kwl);
   ossimFormosatModel(const ossimFormosatModel& rhs);

   virtual ~ossimFormosatModel();

   enum AdjustParamIndex
   {
      ROLL_OFFSET = 0,
      PITCH_OFFSET,
      YAW_OFFSET,
      ROLL_RATE,
      PITCH_RATE,
      YAW_RATE,
      FOCAL_LEN_OFFSET,
      NUM_ADJUSTABLE_PARAMS // not an index
   };   

   virtual ossimObject* dup() const;

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

   virtual bool saveState(ossimKeywordlist& kwl,
                          const char* prefix=NULL) const;
   
   virtual bool loadState(const ossimKeywordlist& kwl,
                          const char* prefix=NULL);

   static void writeGeomTemplate(ostream& os);

   virtual void lineSampleHeightToWorld(const ossimDpt& image_point,
                                        const ossim_float64& heightEllipsoid,
                                        ossimGpt& worldPoint) const;
   
   virtual void imagingRay(const ossimDpt& image_point,
                           ossimEcefRay&   image_ray) const;

   virtual void updateModel();

   inline virtual bool useForward()const {return false;} 
   virtual bool setupOptimizer(const ossimString& init_file); 

   bool initFromMetadata(ossimFormosatDimapSupportData* sd);
protected:

   void initAdjustableParameters();
   
   void loadGeometry(FILE*);
   void loadSupportData();
   //void computeSatToOrbRotation(ossim_float64 t)const;
   void computeSatToOrbRotation(NEWMAT::Matrix& result, ossim_float64 t)const;

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

   ossimRefPtr<ossimFormosatDimapSupportData> theSupportData;

   //---
   // Image constant parameters:
   //---
   ossimFilename  theMetaDataFile;
   ossim_float64  theIllumAzimuth;  
   ossim_float64  theIllumElevation;
   ossim_float64  theSatAzimuth;
   ossim_float64  thePositionError;
   ossim_float64  theRefImagingTime;

   ossim_float64  theRefImagingTimeLine;
   
   ossim_float64  theLineSamplingPeriod;
   ossimDpt       theSpotSubImageOffset;
//   mutable NEWMAT::Matrix theSatToOrbRotation;
//   mutable NEWMAT::Matrix theOrbToEcfRotation;

   //---
   // Adjustable parameters:
   //---
   ossim_float64  theRollOffset;      // degrees
   ossim_float64  thePitchOffset;     // degrees
   ossim_float64  theYawOffset;       // degrees
   ossim_float64  theRollRate;        // degrees/sec
   ossim_float64  thePitchRate;       // degrees/sec
   ossim_float64  theYawRate;         // degrees/sec
   ossim_float64  theFocalLenOffset;  // percent deviation from nominal


TYPE_DATA
};
}
#endif /* #ifndef ossimFormosatModel_HEADER */