PhotometricFunction.h

#ifndef PhotometricFunction_h
#define PhotometricFunction_h

#include "IString.h"
#include "Camera.h"
#include "DbProfile.h"
#include "SpecialPixel.h"

#include <iostream>
#include <sstream>
#include <iomanip>

namespace Isis {

  template<typename T> inline T MIN( const T &A, const T &B ) {
    if (A < B) {
      return (A);
    }
    else {
      return (B);
    }
  }

  template<typename T> inline T MAX( const T &A, const T &B ) {
    if (A > B) {
      return (A);
    }
    else {
      return (B);
    }
  }

  class PvlObject;
  class Camera;

  class PhotometricFunction {

    public:
      PhotometricFunction(PvlObject &pvl, Cube &cube, bool useCamera);
      virtual ~PhotometricFunction();

      void setCamera(Camera *cam);
      static QString algorithmName( const PvlObject &pvl );
      virtual double compute( const double &line, const double &sample, int band = 1, bool useDem = false);
      virtual double photometry( double i, double e, double g, int band = 1 ) const = 0;
      virtual void report( PvlContainer &pvl ) = 0;
      virtual void setMinimumIncidenceAngle( double angle );
      virtual void setMaximumIncidenceAngle( double angle );
      virtual void setMinimumEmissionAngle( double angle );
      virtual void setMaximumEmissionAngle( double angle );
      virtual void setMinimumPhaseAngle( double angle );
      virtual void setMaximumPhaseAngle( double angle );
      virtual double minimumIncidenceAngle();
      virtual double maximumIncidenceAngle();
      virtual double minimumEmissionAngle();
      virtual double maximumEmissionAngle();
      virtual double minimumPhaseAngle();
      virtual double maximumPhaseAngle();


    protected:

      Camera *m_camera;               //<! Camera used for calculating photmetric angles
      double m_iRef;                  //<! Incidence refernce angle
      double m_eRef;                  //<! Emission reference angle
      double m_gRef;                  //<! Phase reference angle
      double m_minimumIncidenceAngle; //<! The minimum incidence angle to perform computations
      double m_maximumIncidenceAngle; //<! The maximum incidence angle to perform computations
      double m_minimumEmissionAngle;  //<! The minimum emission angle to perform computations
      double m_maximumEmissionAngle;  //<! The maximum emission angle to perform computations
      double m_minimumPhaseAngle;     //<! The minimum phase angle to perform computations
      double m_maximumPhaseAngle;     //<! The maximum phase angle to perform computations
      DbProfile m_normProf;           //<! Parameters for the normalization model

      template<typename T>
      T ConfKey( const DbProfile &conf, const QString &keyname, const T &defval, int index = 0 ) const {
        if (!conf.exists(keyname)) {
          return (defval);
        }
        if (conf.count(keyname) < index) {
          return (defval);
        }
         return conf.value(keyname, index);
      }
  };
};
#endif