ControlNetStatistics.h

#ifndef _CONTROLNETSTATISTICS_H_
#define _CONTROLNETSTATISTICS_H_
 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include <QMap>
#include <QVector>
 
#include "Progress.h"
#include "PvlGroup.h"
#include "SerialNumberList.h"
#include "Statistics.h"
 
namespace Isis {
  class ControlNet;
  class Progress;
  class PvlGroup;
 
  class ControlNetStatistics {
    public:
      ControlNetStatistics(ControlNet *pCNet, const QString &psSerialNumFile, Progress *pProgress = 0);
 
      ControlNetStatistics(ControlNet *pCNet, Progress *pProgress = 0);
 
      ~ControlNetStatistics();
 
      enum ePointDetails { total, ignore, locked, fixed, constrained, freed };
      static const int numPointDetails = 6;
 
      enum ePointIntStats { totalPoints, validPoints, ignoredPoints, fixedPoints, constrainedPoints, freePoints, editLockedPoints,
                            totalMeasures, validMeasures, ignoredMeasures, editLockedMeasures };
      static const int numPointIntStats = 11;
 
      enum ePointDoubleStats { avgResidual, minResidual, maxResidual, minLineResidual, maxLineResidual, minSampleResidual, maxSampleResidual,
                               avgPixelShift, minPixelShift, maxPixelShift, minLineShift, maxLineShift, minSampleShift, maxSampleShift,
                               minGFit, maxGFit, minEccentricity, maxEccentricity, minPixelZScore, maxPixelZScore};
      static const int numPointDblStats = 20;
 
      enum ImageStats { imgSamples, imgLines, imgTotalPoints, imgIgnoredPoints, imgFixedPoints, imgLockedPoints, imgLocked,
                        imgConstrainedPoints, imgFreePoints, imgConvexHullArea, imgConvexHullRatio };
      static const int numImageStats = 11;
 
      void GenerateImageStats();
 
      void PrintImageStats(const QString &psImageFile);
 
      QVector<double> GetImageStatsBySerialNum(QString psSerialNum) const;
 
      void GeneratePointStats(const QString &psPointFile);
 
      void GenerateControlNetStats(PvlGroup &pStatsGrp);
 
      int NumValidPoints() const {
        return (*mPointIntStats.find(validPoints));
      }
 
      int NumFixedPoints() const {
        return (*mPointIntStats.find(fixedPoints));
      }
 
      int NumConstrainedPoints() const {
        return (*mPointIntStats.find(constrainedPoints));
      }
 
      int NumFreePoints() const {
        return (*mPointIntStats.find(freePoints));
      }
 
      int NumIgnoredPoints() const {
        return (*mPointIntStats.find(ignoredPoints));
      }
 
      int NumEditLockedPoints() const {
        return (*mPointIntStats.find(editLockedPoints));
      }
 
      int NumMeasures() const {
        return (*mPointIntStats.find(totalMeasures));
      }
 
      int NumValidMeasures() const {
        return (*mPointIntStats.find(validMeasures));
      }
 
      int NumIgnoredMeasures() const {
        return (*mPointIntStats.find(ignoredMeasures));
      }
 
      int NumEditLockedMeasures() const {
        return (*mPointIntStats.find(editLockedMeasures));
      }
 
      double GetAverageResidual() const {
        return (*mPointDoubleStats.find(avgResidual));
      }
 
      double GetMinimumResidual() const {
        return (*mPointDoubleStats.find(minResidual));
      }
 
      double GetMaximumResidual() const {
        return (*mPointDoubleStats.find(maxResidual));
      }
 
      double GetMinLineResidual() const {
        return (*mPointDoubleStats.find(minLineResidual));
      }
 
      double GetMinSampleResidual() const {
        return (*mPointDoubleStats.find(minSampleResidual));
      }
 
      double GetMaxLineResidual() const {
        return (*mPointDoubleStats.find(maxLineResidual));
      }
 
      double GetMaxSampleResidual() const {
        return (*mPointDoubleStats.find(maxSampleResidual));
      }
 
      double GetMinLineShift() const {
        return (*mPointDoubleStats.find(minLineShift));
      }
 
      double GetMaxLineShift() const {
        return (*mPointDoubleStats.find(maxLineShift));
      }
 
      double GetMinSampleShift() const {
        return (*mPointDoubleStats.find(minSampleShift));
      }
 
      double GetMaxSampleShift() const {
        return (*mPointDoubleStats.find(maxSampleShift));
      }
 
      double GetMinPixelShift() const {
        return (*mPointDoubleStats.find(minPixelShift));
      }
 
      double GetMaxPixelShift() const {
        return (*mPointDoubleStats.find(maxPixelShift));
      }
 
      double GetAvgPixelShift() const {
        return (*mPointDoubleStats.find(avgPixelShift));
      }
 
    protected:
      SerialNumberList mSerialNumList;           
      ControlNet *mCNet;                         
      Progress *mProgress;                       
 
    private:
      QMap<int, int> mPointIntStats;           
      QMap<int, double> mPointDoubleStats;     
      QMap<QString, QVector<double> > mImageMap; 
      QMap<QString, bool> mSerialNumMap;        
 
      void GetPointIntStats();
 
      void GetPointDoubleStats();
 
      void UpdateMinMaxStats(const Statistics & stats,
                             ePointDoubleStats min,
                             ePointDoubleStats max);
 
      void InitPointDoubleStats();
 
      void InitSerialNumMap();
 
      int numCNetImages;
 
      Statistics mConvexHullStats, mConvexHullRatioStats; 
  };
}
#endif