NonParametricXMLMieModel.h

Go to the documentation of this file.

#ifndef _atm_NonParametricXMLMieModel_h_
#define _atm_NonParametricXMLMieModel_h_

#include <atm/VMieModel.h>
#include <utl/ShadowPtr.h>
#include <utl/TabulatedFunction.h>

#include <vector>

namespace atm {

  class NonParametricXMLMieModel : public atm::VMieModel {

  public:
    virtual ~NonParametricXMLMieModel() { }

    atm::ScatteringResult EvaluateMieScattering(const utl::Point&  xA,
                                                const utl::Point&  xB,
                                                const double angle,
                                                const double distance,
                                                const std::vector<double>& wLength) const;

    atm::ScatteringResult EvaluateMieScattering(const utl::Point&  xA,
                                                const utl::Point&  xB,
                                                const double angle,
                                                const double distance,
                                                const atm::AttenuationResult& mieAttenuation) const;

    atm::AttenuationResult EvaluateMieAttenuation(const utl::Point& xInit,
                                                  const utl::Point& xFinal,
                                                  const std::vector<double>& wLength) const;

    double EvaluateMieScattering(const utl::Point& xA,
                                 const utl::Point& xB,
                                 const double angle,
                                 const double distance,
                                 const double wLength) const;

    double EvaluateMieScattering(const utl::Point& xA,
                                 const utl::Point& xB,
                                 const double angle,
                                 const double distance,
                                 const double wLength,
                                 const double mieAttenuation) const;

    double EvaluateMieAttenuation(const utl::Point& xInit,
                                  const utl::Point& xFinal,
                                  const double wLength) const;

    double GetVerticalAerosolOpticalDepth(const unsigned int eyeId,
                                          const double altitude) const;

    double GetAttenuationLength(const utl::Point& p,
                                const double wLength) const;

    double EvaluateScatteringAngle(const utl::Point& p,
                                   const double angle,
                                   const double wLength) const;

    void Init();

    bool HasData() const { return true; }

    void SetOpticalDepth(const utl::TabulatedFunction& vaod, const utl::TabulatedFunction& alpha);

    typedef std::vector<double> PhaseFunction;

    void SetPhaseFunction(const PhaseFunction& phaseFunc);
  private:

    utl::ShadowPtr<utl::TabulatedFunction> fVAODVsHeight;

    utl::ShadowPtr<utl::TabulatedFunction> fAlphaVsHeight;

    PhaseFunction fPF;

    double fGamma;
    double fGammaError;

    REGISTER_MODEL(atm::VMieModel, "NonParametricXML", NonParametricXMLMieModel);
  };

}

#endif