Affine.h

Go to the documentation of this file.

#ifndef Affine_h
#define Affine_h

#include <vector>
#include "tnt/tnt_array2d.h"

namespace Isis {
  class Affine {
    public:
      typedef TNT::Array2D<double> AMatrix; 
      Affine();
      Affine(const AMatrix &a);
      ~Affine();
      void Solve(const double x[], const double y[],
                 const double xp[], const double yp[], int n);
      static AMatrix getIdentity();
      void Identity();
      void Translate(double tx, double ty);
      void Rotate(double rot);
      void Scale(double scaleFactor);

      void Compute(double x, double y);

      double xp() const {
        return p_xp;
      };

      double yp() const {
        return p_yp;
      };

      void ComputeInverse(double xp, double yp);

      double x() const {
        return p_x;
      };

      double y() const {
        return p_y;
      };

      std::vector<double> Coefficients(int var);
      std::vector<double> InverseCoefficients(int var);
      
      AMatrix Forward() const {
        return (p_matrix.copy());
      }
      
      AMatrix Inverse() const {
        return (p_invmat.copy());
      }

    private:
      AMatrix p_matrix;   
      AMatrix p_invmat;   

      double p_x;         
      double p_y;         
      double p_xp;        
      double p_yp;        

      void checkDims(const AMatrix &am) const;
      AMatrix invert(const AMatrix &a) const;
  };
};

#endif