Affine.h

#ifndef Affine_h
#define Affine_h
 
/* SPDX-License-Identifier: CC0-1.0 */
 
#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