NormModel.cpp

#include "Pvl.h"
#include "NormModel.h"
#include "PhotoModel.h"
#include "AtmosModel.h"
#include "IException.h"

namespace Isis {
  NormModel::NormModel(Pvl &pvl, PhotoModel &pmodel) {
    p_normAlgorithmName = "Unknown";
    p_normPM = &pmodel;
    p_normAM = NULL;
    p_normWavelength = 1.0;
  }

  NormModel::NormModel(Pvl &pvl, PhotoModel &pmodel, AtmosModel &amodel) {
    p_normAlgorithmName = "Unknown";
    p_normPM = &pmodel;
    p_normAM = &amodel;
    p_normWavelength = 1.0;
  }

  void NormModel::SetNormWavelength(double wavelength) {
    p_normWavelength = wavelength;
  }

  void NormModel::CalcNrmAlbedo(double pha, double inc, double ema, double dn,
                                double &albedo, double &mult, double &base) {

    // Check validity of photometric angles
    //if (pha > 180.0 || inc > 90.0 || ema > 90.0 || pha < 0.0 ||
    //    inc < 0.0 || ema < 0.0) {
    //  std::string msg = "Invalid photometric angles";
    //  throw iException::Message(iException::Programmer,msg,_FILEINFO_);
    //}

    // Perform normalization
    NormModelAlgorithm(pha, inc, ema, dn, albedo, mult, base);
  }

  void NormModel::CalcNrmAlbedo(double pha, double inc, double ema,
                                double deminc, double demema, double dn, double &albedo,
                                double &mult, double &base) {

    // Check validity of photometric angles
    //if (pha > 180.0 || inc > 90.0 || ema > 90.0 || pha < 0.0 ||
    //    inc < 0.0 || ema < 0.0) {
    //  std::string msg = "Invalid photometric angles";
    //  throw iException::Message(iException::Programmer,msg,_FILEINFO_);
    //}

    // Perform normalization
    NormModelAlgorithm(pha, inc, ema, deminc, demema, dn, albedo, mult, base);
  }
}