ProcessPolygons.h

#ifndef ProcessPolygons_h
#define ProcessPolygons_h

#include <geos/geom/Coordinate.h>
#include <geos/geom/MultiPolygon.h>
#include <geos/geom/Polygon.h>

#include "Brick.h"
#include "FileName.h"
#include "Process.h"
#include "ProjectionFactory.h"

namespace Isis {

  class ProcessPolygons : public Isis::Process {

    public:
      ProcessPolygons();

      void SetStatCubes(const QString &parameter, const int nsamps,
                         const int nlines, int nbands = 1);

      void SetStatCubes(const QString &avgFileName, const QString
                         &countFileName, Isis::CubeAttributeOutput &atts,
                         const int nsamps, const int nlines, int nbands = 1);
      void SetIntersectAlgorithm(const bool useCenter);

      Isis::Cube *AppendOutputCube(const QString &avgFileName,
                                   const QString &countFileName = "");

      void Rasterize(std::vector<double> &samples,
                     std::vector<double> &lines,
                     std::vector<double> &values);

      void Rasterize(std::vector<double> &samples,
                     std::vector<double> &lines,
                     int &band, double &value);

      void EndProcess();
      void Finalize();

    private:
      void FillPolygon(int Flag);
      void GetPolygonCoords();

      bool m_useCenter;
      std::vector<double> m_sampleVertices, m_lineVertices, m_dns;
      int m_band;
      Brick *m_average;
      Brick *m_count;
      geos::geom::Polygon *m_imagePoly;

  };

};

#endif