ControlNetStatistics.h

Go to the documentation of this file.

#ifndef _CONTROLNETSTATISTICS_H_
#define _CONTROLNETSTATISTICS_H_

#include <map>
#include <iostream>
#include <vector>
#include "Progress.h"
#include "PvlGroup.h"
#include "SerialNumberList.h"
#include "Statistics.h"

namespace Isis {
  class ControlNet;
  class ControlCubeGraphNode;
  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);

      std::vector<double> GetImageStatsBySerialNum(QString psSerialNum) const;

      void GeneratePointStats(const QString &psPointFile);

      void GenerateControlNetStats(PvlGroup &pStatsGrp);

      int NumValidPoints() const {
        return (*mPointIntStats.find(validPoints)).second;
      }

      int NumFixedPoints() const {
        return (*mPointIntStats.find(fixedPoints)).second;
      }

      int NumConstrainedPoints() const {
        return (*mPointIntStats.find(constrainedPoints)).second;
      }

      int NumFreePoints() const {
        return (*mPointIntStats.find(freePoints)).second;
      }

      int NumIgnoredPoints() const {
        return (*mPointIntStats.find(ignoredPoints)).second;
      }

      int NumEditLockedPoints() const {
        return (*mPointIntStats.find(editLockedPoints)).second;
      }

      int NumMeasures() const {
        return (*mPointIntStats.find(totalMeasures)).second;
      }

      int NumValidMeasures() const {
        return (*mPointIntStats.find(validMeasures)).second;
      }

      int NumIgnoredMeasures() const {
        return (*mPointIntStats.find(ignoredMeasures)).second;
      }

      int NumEditLockedMeasures() const {
        return (*mPointIntStats.find(editLockedMeasures)).second;
      }

      double GetAverageResidual() const {
        return (*mPointDoubleStats.find(avgResidual)).second;
      }

      double GetMinimumResidual() const {
        return (*mPointDoubleStats.find(minResidual)).second;
      }

      double GetMaximumResidual() const {
        return (*mPointDoubleStats.find(maxResidual)).second;
      }

      double GetMinLineResidual() const {
        return (*mPointDoubleStats.find(minLineResidual)).second;
      }

      double GetMinSampleResidual() const {
        return (*mPointDoubleStats.find(minSampleResidual)).second;
      }

      double GetMaxLineResidual() const {
        return (*mPointDoubleStats.find(maxLineResidual)).second;
      }

      double GetMaxSampleResidual() const {
        return (*mPointDoubleStats.find(maxSampleResidual)).second;
      }

      double GetMinLineShift() const {
        return (*mPointDoubleStats.find(minLineShift)).second;
      }

      double GetMaxLineShift() const {
        return (*mPointDoubleStats.find(maxLineShift)).second;
      }

      double GetMinSampleShift() const {
        return (*mPointDoubleStats.find(minSampleShift)).second;
      }

      double GetMaxSampleShift() const {
        return (*mPointDoubleStats.find(maxSampleShift)).second;
      }

      double GetMinPixelShift() const {
        return (*mPointDoubleStats.find(minPixelShift)).second;
      }

      double GetMaxPixelShift() const {
        return (*mPointDoubleStats.find(maxPixelShift)).second;
      }

      double GetAvgPixelShift() const {
        return (*mPointDoubleStats.find(avgPixelShift)).second;
      }

    protected:
      SerialNumberList mSerialNumList;           
      ControlNet *mCNet;                         
      Progress *mProgress;                       
      QList<ControlCubeGraphNode *> mCubeGraphNodes;

    private:
      std::map<int, int> mPointIntStats;           
      std::map<int, double> mPointDoubleStats;     
      std::map<QString, std::vector<double> > mImageMap; 
      std::map<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