ProcessRubberSheet.h

Go to the documentation of this file.

#ifndef ProcessRubberSheet_h
#define ProcessRubberSheet_h

#include "Process.h"
#include "Buffer.h"
#include "Transform.h"
#include "Interpolator.h"
#include "Portal.h"
#include "TileManager.h"

namespace Isis {
  class ProcessRubberSheet : public Isis::Process {
    public:
      ProcessRubberSheet(int startSize = 128, int endSize = 8) {
        p_bandChangeFunct = NULL;
        p_forceSamp = Isis::Null;
        p_forceLine = Isis::Null;
        p_startQuadSize = startSize;
        p_endQuadSize = endSize;
      };

      ~ProcessRubberSheet() {};

      // Line Processing method for one input and output cube
      void StartProcess(Isis::Transform &trans, Isis::Interpolator &interp);

      // Register a function to be called when the band number changes
      void BandChange(void (*funct)(const int band));

      void ForceTile(double Samp, double Line) {
        p_forceSamp = Samp;
        p_forceLine = Line;
      }

      void SetTiling(int start, int end) {
        p_startQuadSize = start;
        p_endQuadSize = end;
      }

    private:

      class Quad {
        public:
          int slineTile; 
          int ssampTile; 
          int sline;     
          int ssamp;     
          int eline;     
          int esamp;     
      };

      void ProcessQuad(std::vector<Quad *> &quadTree, Isis::Transform &trans,
                       std::vector< std::vector<double> > &lineMap,
                       std::vector< std::vector<double> > &sampMap);

      void SplitQuad(std::vector<Quad *> &quadTree);
      void SlowQuad(std::vector<Quad *> &quadTree, Isis::Transform &trans,
                    std::vector< std::vector<double> > &lineMap, 
                    std::vector< std::vector<double> > &sampMap);
      double Det4x4(double m[4][4]);
      double Det3x3(double m[3][3]);

      // SlowGeom method is never used but saved for posterity
      void SlowGeom(Isis::TileManager &otile, Isis::Portal &iportal,
                    Isis::Transform &trans, Isis::Interpolator &interp);
      void QuadTree(Isis::TileManager &otile, Isis::Portal &iportal,
                    Isis::Transform &trans, Isis::Interpolator &interp,
                    bool useLastTileMap);

      bool TestLine(Isis::Transform &trans, int ssamp, int esamp, int sline, 
                    int eline, int increment);

      void (*p_bandChangeFunct)(const int band);

      std::vector< std::vector<double> > p_sampMap; 
      std::vector< std::vector<double> > p_lineMap; 

      double p_forceSamp; 
      double p_forceLine; 

      int p_startQuadSize; 
      int p_endQuadSize;   
  };
};

#endif