ossimLeastSquaresBilin.h

Go to the documentation of this file.

//*******************************************************************
//
// License:  See top level LICENSE.txt file.
//
// Author:  Garrett Potts (gpotts@imagelinks.com
//
// Description: Source code produced by Dave Knopp
//
//*******************************************************************
//  $Id: ossimLeastSquaresBilin.h 9968 2006-11-29 14:01:53Z gpotts $
#ifndef    ossimLeastSquaresBilin_INCLUDE
#define    ossimLeastSquaresBilin_INCLUDE
#include <ossim/base/ossimConstants.h>
#include <ossim/matrix/newmat.h>

class OSSIMDLLEXPORT ossimLeastSquaresBilin
{
public:

   ossimLeastSquaresBilin(const ossimLeastSquaresBilin &);
   ossimLeastSquaresBilin();
   
   ossimLeastSquaresBilin & operator = (const ossimLeastSquaresBilin &);
   
   virtual ~ossimLeastSquaresBilin();
   

   virtual void clear();

   virtual void addSample(double x, double yy, double zmea);
   
   virtual bool getLSParms(double& pa, double& pb_x, double& pc_y, double& pd_xy)const;

   virtual void setLSParams(double pa, //constant
                            double pb_x, // linear coefficient of x,
                            double pc_y, // linear coefficient of y
                            double pd_xy); // coefficient of x*y
                           
   virtual inline double lsFitValue(double xx, double yy) const
      {
         return (bl_a + bl_b*xx + bl_c*yy + bl_d*xx*yy);
      }
   
   bool solveLS();
   
private:
   double bl_a;  

   double bl_b;

   double bl_c;

   double bl_d;   

   NEWMAT::Matrix*  AtA;

   NEWMAT::Matrix*  Atb;
};

#endif //  LeastSquaresBilinilin_INCL_