Selection.h

#ifndef Selection_h
#define Selection_h
 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include "Chip.h"
#include "Statistics.h"
 
typedef struct Ellipse
{
  double A[3];
  double semiMajor;
  double semiMinor;
  double majorAxis[2];
  double minorAxis[2];
  double cen[2];
  double area;
}Ellipse;
 
namespace Isis {
  class Selection {
  public:
    Selection();
    virtual ~Selection();
 
    virtual int select(Chip *inputChip,Chip *selectionChip)=0;  //pure virtual function to be redefined in each child class
 
    //Reduction methods--methods used to trim off outliers in the selection based on aprior knowledge of the expected shape
    virtual int elipticalReduction(Chip *selectionChip, double percent_selected, double play, int patience_limit);    //a general, but slow elliptical trimming is provided, it is virtual so it can redefined to take advantage of differing levels of apriori knowlege of the characteristics of the ellipse
 
    //Observation Methods--methods used to reduce a selection to a single sub-pixel observation
    int centerOfMass(Chip *selectionChip, double *sample, double *line);
    int centerOfMassWeighted(Chip *inputChip, Chip *selectionChip, double *sample, double *line);
 
    //Sub-tasks for elliptical trimming and ellipse fitting
    bool ellipseFrom5Pts(Ellipse *ell,double pts[5][2]);
    double elipsePercentSelected(Chip *selectionChip,Ellipse *ell);
    std::vector<double> minimumBoundingElipse( std::vector< std::vector<int> > pts,Ellipse *ell);  
    bool ellipseFromCubic(Ellipse *ell, double cubic[6]);            
    bool ellipseAxesAreaFromMatrix(Ellipse *ell);
    bool bestFitEllipse(Ellipse *ell, std::vector < std::vector<int> > *pts,double play,unsigned int max_iter);
    bool pointInEllipse(Ellipse *ell, double pt[2],double play);
    bool ellipseInChip(Ellipse *ell, Chip *chip);
    bool ellipseFromCenterAxesAngle(Ellipse *ell, double centerSample, double centerLine, double semiMajor, double semiMinor, double theta);
 
    void selectionEdge(Chip *selectionChip, std::vector < std::vector <int> > *pts);    
  };
};
 
#endif