ossimSmacCallibrationSystem.h

Go to the documentation of this file.

//*******************************************************************
//
// License:  See top level LICENSE.txt file.
//
// Author:  Garrett Potts
//
// Description:
//
// Proecedure for compensation of aerial camera lens distortion
// as computed by the Simultaneous Multiframe Analytical Calibration
// (SMAC) System:
// 
//*******************************************************************
//  $Id$
#ifndef ossimSmacCallibrationSystem_HEADER
#define ossimSmacCallibrationSystem_HEADER
#include <ossim/projection/ossimLensDistortion.h>

class OSSIM_DLL ossimSmacCallibrationSystem : public ossimLensDistortion
{
public:
   ossimSmacCallibrationSystem(const ossimDpt& callibratedCenter = ossimDpt(0.0, 0.0))
   :ossimLensDistortion(callibratedCenter)
   {
      _symmetricRadialDistortionCoefficients[0] = 0;
      _symmetricRadialDistortionCoefficients[1] = 0;
      _symmetricRadialDistortionCoefficients[2] = 0;
      _symmetricRadialDistortionCoefficients[3] = 0;
      _symmetricRadialDistortionCoefficients[4] = 0;
      _decenteringDistortionCoefficients[0]     = 0;
      _decenteringDistortionCoefficients[1]     = 0;
      _decenteringDistortionCoefficients[2]     = 0;
      _decenteringDistortionCoefficients[3]     = 0;
   }
   ossimSmacCallibrationSystem(double k0, double k1, double k2, double k3, double k4,
                               double p1, double p2, double p3, double p4,
                               const ossimDpt& callibratedCenter = ossimDpt(0.0,0.0))
   :ossimLensDistortion(callibratedCenter)
   {
      setSymmetricRadialDistortionCoefficients(k0, k1, k2, k3, k4);
      _decenteringDistortionCoefficients[0]     = p1;
      _decenteringDistortionCoefficients[1]     = p2;
      _decenteringDistortionCoefficients[2]     = p3;
      _decenteringDistortionCoefficients[3]     = p4;
   }
   ossimSmacCallibrationSystem(const ossimSmacCallibrationSystem& src)
   {
      std::copy(src._symmetricRadialDistortionCoefficients,
                src._symmetricRadialDistortionCoefficients+5,
                _symmetricRadialDistortionCoefficients);
      std::copy(src._decenteringDistortionCoefficients,
                src._decenteringDistortionCoefficients+4,
                _decenteringDistortionCoefficients);
   }
   virtual ossimObject* dup()const
   {
      return new ossimSmacCallibrationSystem(*this);
   }
   const double* symmetricRadialDistortionCoefficients()const
   {
      return _symmetricRadialDistortionCoefficients;
   }
   void setSymmetricRadialDistortionCoefficients(double k0, double k1, double k2, double k3, double k4)
   {
      _symmetricRadialDistortionCoefficients[0] = k0;
      _symmetricRadialDistortionCoefficients[1] = k1;
      _symmetricRadialDistortionCoefficients[2] = k2;
      _symmetricRadialDistortionCoefficients[3] = k3;
      _symmetricRadialDistortionCoefficients[4] = k4;
   }
   const double* decenteringDistortionCoefficients()const
   {
      return _decenteringDistortionCoefficients;
   }
   void setDecenteringDistortionCoefficients(double p1, double p2, double p3, double p4)
   {
      _decenteringDistortionCoefficients[0]     = p1;
      _decenteringDistortionCoefficients[1]     = p2;
      _decenteringDistortionCoefficients[2]     = p3;
      _decenteringDistortionCoefficients[3]     = p4;
   }
   const double* radialDistortionCoefficients()const
   {
      return _symmetricRadialDistortionCoefficients;
   }
   virtual void undistort(const ossimDpt& input, ossimDpt& output)const;
   
   static void test();
protected:
   double _symmetricRadialDistortionCoefficients[5];
   
   double _decenteringDistortionCoefficients[4];
   
};

#endif