PhotometricFunction.h

#ifndef PhotometricFunction_h
#define PhotometricFunction_h
 
/* SPDX-License-Identifier: CC0-1.0 */
 
#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