Isis 3 Programmer Reference
Isis::HapkeAtm2 Class Reference

#include <HapkeAtm2.h>

Inheritance diagram for Isis::HapkeAtm2:
Inheritance graph
Collaboration diagram for Isis::HapkeAtm2:
Collaboration graph

Public Member Functions

 HapkeAtm2 (Pvl &pvl, PhotoModel &pmodel)
 
void CalcAtmEffect (double pha, double inc, double ema, double *pstd, double *trans, double *trans0, double *sbar, double *transs)
 Calculate the atmospheric scattering effect using photometric angle information. More...
 
virtual void SetStandardConditions (bool standard)
 Used to calculate atmosphere at standard conditions. More...
 
void GenerateAhTable ()
 This method computes the values of the atmospheric Ah table and sets the properties of the atmospheric Ah spline. More...
 
void GenerateHahgTables ()
 This method computes the values of the atmospheric Hahg and Hahg0 tables and sets the properties of the atmospheric Hahg and Hahg0 splines. More...
 
void GenerateHahgTablesShadow ()
 This method is a modified version of the GenerateHahgTables method and is used solely for shadow modeling. More...
 
void SetAtmosAtmSwitch (const int atmswitch)
 Set the switch that controls the function that will be integrated. More...
 
void SetAtmosBha (const double bha)
 Set the Anisotropic Atmospheric function parameter. More...
 
void SetAtmosHga (const double hga)
 Set the Hapke atmospheric function parameter. More...
 
void SetAtmosInc (const double inc)
 Set the incidence angle. More...
 
void SetAtmosNulneg (const string nulneg)
 Set the Atmospheric function parameter. More...
 
void SetAtmosPhi (const double phi)
 Set the azimuth angle. More...
 
void SetAtmosTau (const double tau)
 Set the Atmospheric function parameter. More...
 
void SetAtmosTauref (const double tauref)
 Set the Atmospheric function parameter. More...
 
void SetAtmosWha (const double wha)
 Set the Atmospheric function parameter. More...
 
void SetAtmosHnorm (const double hnorm)
 Set the Atmospheric function parameter. More...
 
void SetAtmosIord (const string offset)
 Set additive offset in fit. More...
 
void SetAtmosEstTau (const string esttau)
 Estimate the optical depth tau using shadows. More...
 
string AlgorithmName () const
 Return atmospheric algorithm name. More...
 
bool AtmosAdditiveOffset () const
 Allow additive offset in fit? More...
 
double AtmosHnorm () const
 Return atmospheric Hnorm value. More...
 
double AtmosBha () const
 Return atmospheric Bha value. More...
 
double AtmosTau () const
 Return atmospheric Tau value. More...
 
double AtmosWha () const
 Return atmospheric Wha value. More...
 
double AtmosHga () const
 Return atmospheric Hga value. More...
 
double AtmosTauref () const
 Return atmospheric Tauref value. More...
 
bool AtmosNulneg () const
 Return atmospheric Nulneg value. More...
 
double AtmosAb () const
 Return atmospheric Ab value. More...
 
double AtmosHahgsb () const
 Return atmospheric Hahgsb value. More...
 
int AtmosNinc () const
 Return atmospheric Ninc value. More...
 
double AtmosMunot () const
 Return atmospheric Munot value. More...
 
vector< double > AtmosIncTable ()
 Return atmospheric IncTable value. More...
 
vector< double > AtmosAhTable ()
 Return atmospheric AhTable value. More...
 
vector< double > AtmosHahgtTable ()
 Return atmospheric HahgtTable value. More...
 
vector< double > AtmosHahgt0Table ()
 Return atmospheric Hahgt0Table value. More...
 
NumericalApproximation AtmosAhSpline ()
 If GenerateAhTable() has been called this returns a clamped cubic spline of the data set (p_atmosIncTable, p_atmosAhTable) with first derivatives of the endpoints equal to 1.0e30. More...
 
NumericalApproximation AtmosHahgtSpline ()
 If GenerateHahgTables() has been called this returns a clamped cubic spline of the data set (p_atmosIncTable, p_atmosHahgtTable) with first derivatives of the endpoints equal to 1.0e30. More...
 
NumericalApproximation AtmosHahgt0Spline ()
 If GenerateHahgTables() has been called this returns a clamped cubic spline of the data set (p_atmosIncTable, p_atmosHahgt0Table) with first derivatives of the endpoints equal to 1.0e30. More...
 

Static Public Member Functions

static double G11Prime (double tau)
 Perform Chandra and Van de Hulst's series approximation for the g'11 function needed in second order scattering theory. More...
 
static double Ei (double x)
 This routine computes the exponential integral, Ei(x). More...
 
static double En (unsigned int n, double x)
 This routine evaluates the generalized exponential integral, En(x). More...
 

Protected Member Functions

virtual void AtmosModelAlgorithm (double phase, double incidence, double emission)
 Henyey-Greenstein atmos scattering in the 1st approximation. More...
 
void SetAtmosNulneg (bool nulneg)
 
void SetAtmosIord (bool offset)
 
void SetAtmosEstTau (bool esttau)
 
void SetAlgorithmName (string name)
 
void SetOldTau (double tau)
 
void SetOldWha (double wha)
 
PhotoModelGetPhotoModel () const
 
bool StandardConditions () const
 
bool TauOrWhaChanged () const
 Checks whether tau or wha have changed. More...
 
double Eulgam () const
 

Protected Attributes

int p_atmosAtmSwitch
 
int p_atmosNinc
 
double p_atmosBha
 
double p_atmosBhasave
 
double p_atmosHgasave
 
double p_atmosTauref
 
double p_atmosTausave
 
double p_atmosWhasave
 
double p_pstd
 Pure atmospheric-scattering term. More...
 
double p_trans
 Transmission of surface reflected light through the atmosphere overall. More...
 
double p_trans0
 Transmission of surface reflected light through the atmosphere with no scatterings in the atmosphere. More...
 
double p_transs
 Transmission of light that must be subtracted from the flat surface model to get the shadow model. More...
 
double p_sbar
 Illumination of the ground by the sky. More...
 
double p_atmosHga
 
double p_atmosTau
 
double p_atmosWha
 
double p_atmosAb
 
double p_atmosHnorm
 Atmospheric shell thickness normalized to planet radius. More...
 
bool p_atmosAddOffset
 Allow additive offset in fit. More...
 
bool p_atmosEstTau
 Estimate optical depth tau using shadows. More...
 
vector< double > p_atmosIncTable
 
vector< double > p_atmosAhTable
 
double p_atmosHahgsb
 
vector< double > p_atmosHahgtTable
 
vector< double > p_atmosHahgt0Table
 
double p_atmosInc
 
double p_atmosPhi
 
double p_atmosMunot
 
double p_atmosSini
 
double p_atmosCosphi
 
double p_atmosEulgam
 
NumericalApproximation p_atmosAhSpline
 Spline object for the atmospheric Ah Table. Properties are set in GenerateAhTable(). More...
 
NumericalApproximation p_atmosHahgtSpline
 Spline object for the atmospheric Hahg Table. Properties are set in GenerateHahgTables(). More...
 
NumericalApproximation p_atmosHahgt0Spline
 Spline object for the atmospheric Hahg0 Table. Properties are set in GenerateHahgTables(). More...
 

Private Attributes

double p_e1
 
double p_e1_2
 
double p_e2
 
double p_e3
 
double p_e4
 
double p_e5
 
double p_em
 
double p_e
 
double p_g11p
 
double p_g12p
 
double p_g13p
 
double p_g14p
 
double p_x0
 
double p_y0
 
double p_wha2
 
double p_g12
 
double p_g13
 
double p_g14
 
double p_f1
 
double p_f2
 
double p_f3
 
double p_f4
 
double p_f1m
 
double p_f2m
 
double p_f3m
 
double p_f4m
 
double p_delta
 
double p_fixcon
 
double p_alpha0
 
double p_alpha1
 
double p_alpha2
 
double p_beta0
 
double p_beta1
 
double p_beta2
 
double p_gammax
 
double p_gammay
 

Detailed Description

Author
2000-07-07 Randy Kirk
History:

2000-07-07 Randy Kirk - USGS, Flagstaff - Original code

2007-02-20 Janet Barrett - Imported from Isis2.

2007-08-15 Steven Lambright - Refactored code

2008-03-07 Janet Barrett - Moved code to set standard conditions to the AtmosModel class

2008-06-18 Stuart Sides - Fixed doc error

2008-11-05 Jeannie Walldren - Modified references to NumericalMethods class. Added documentation from Isis2.

2011-08-18 Sharmila Prasad Moved common HNORM to base AtmosModel

2011-12-19 Janet Barrett - Added code to estimate the shadow brightness value (transs). Also got rid of unnecessary check for identical photometric angle values between successive calls. This check should only be made in the photometric models.

Definition at line 56 of file HapkeAtm2.h.

Member Function Documentation

◆ AlgorithmName()

string Isis::AtmosModel::AlgorithmName ( ) const
inlineinherited

Return atmospheric algorithm name.

Definition at line 116 of file AtmosModel.h.

◆ AtmosAb()

double Isis::AtmosModel::AtmosAb ( ) const
inlineinherited

Return atmospheric Ab value.

Definition at line 155 of file AtmosModel.h.

Referenced by Isis::AlbedoAtm::NormModelAlgorithm().

◆ AtmosAdditiveOffset()

bool Isis::AtmosModel::AtmosAdditiveOffset ( ) const
inlineinherited

Allow additive offset in fit?

Definition at line 121 of file AtmosModel.h.

◆ AtmosAhSpline()

NumericalApproximation Isis::AtmosModel::AtmosAhSpline ( )
inlineinherited

If GenerateAhTable() has been called this returns a clamped cubic spline of the data set (p_atmosIncTable, p_atmosAhTable) with first derivatives of the endpoints equal to 1.0e30.

Otherwise, it is a natural cubic spline with an empty data set.

Returns
NumericalApproximation Cubic spline
History:
2008-11-05 Jeannie Walldren - Original version

Definition at line 199 of file AtmosModel.h.

Referenced by Isis::AlbedoAtm::AlbedoAtm(), and Isis::AlbedoAtm::NormModelAlgorithm().

◆ AtmosAhTable()

vector<double> Isis::AtmosModel::AtmosAhTable ( )
inlineinherited

Return atmospheric AhTable value.

Definition at line 176 of file AtmosModel.h.

◆ AtmosBha()

double Isis::AtmosModel::AtmosBha ( ) const
inlineinherited

Return atmospheric Bha value.

Definition at line 131 of file AtmosModel.h.

◆ AtmosHahgsb()

double Isis::AtmosModel::AtmosHahgsb ( ) const
inlineinherited

Return atmospheric Hahgsb value.

Definition at line 159 of file AtmosModel.h.

◆ AtmosHahgt0Spline()

NumericalApproximation Isis::AtmosModel::AtmosHahgt0Spline ( )
inlineinherited

If GenerateHahgTables() has been called this returns a clamped cubic spline of the data set (p_atmosIncTable, p_atmosHahgt0Table) with first derivatives of the endpoints equal to 1.0e30.

Otherwise, it is a natural cubic spline with an empty data set.

Returns
NumericalApproximation Cubic spline
History:
2008-11-05 Jeannie Walldren - Original version

Definition at line 227 of file AtmosModel.h.

◆ AtmosHahgt0Table()

vector<double> Isis::AtmosModel::AtmosHahgt0Table ( )
inlineinherited

Return atmospheric Hahgt0Table value.

Definition at line 184 of file AtmosModel.h.

◆ AtmosHahgtSpline()

NumericalApproximation Isis::AtmosModel::AtmosHahgtSpline ( )
inlineinherited

If GenerateHahgTables() has been called this returns a clamped cubic spline of the data set (p_atmosIncTable, p_atmosHahgtTable) with first derivatives of the endpoints equal to 1.0e30.

Otherwise, it is a natural cubic spline with an empty data set.

Returns
NumericalApproximation Cubic spline
History:
2008-11-05 Jeannie Walldren - Original version

Definition at line 213 of file AtmosModel.h.

◆ AtmosHahgtTable()

vector<double> Isis::AtmosModel::AtmosHahgtTable ( )
inlineinherited

Return atmospheric HahgtTable value.

Definition at line 180 of file AtmosModel.h.

◆ AtmosHga()

double Isis::AtmosModel::AtmosHga ( ) const
inlineinherited

Return atmospheric Hga value.

Definition at line 143 of file AtmosModel.h.

Referenced by Isis::NumericalAtmosApprox::InrFunc2Bint().

◆ AtmosHnorm()

double Isis::AtmosModel::AtmosHnorm ( ) const
inlineinherited

Return atmospheric Hnorm value.

Definition at line 126 of file AtmosModel.h.

◆ AtmosIncTable()

vector<double> Isis::AtmosModel::AtmosIncTable ( )
inlineinherited

Return atmospheric IncTable value.

Definition at line 172 of file AtmosModel.h.

◆ AtmosModelAlgorithm()

void Isis::HapkeAtm2::AtmosModelAlgorithm ( double  phase,
double  incidence,
double  p_emission 
)
protectedvirtual

Henyey-Greenstein atmos scattering in the 1st approximation.

Isotropic atmospheric scattering in the second approximation, with corrections to the singly-scattered terms (in the spirit of Hapke's photometric function for surfaces) for a strongly anisotropic single-particle phase function. The particular phase function implemented is a single-term Henyey-Greenstein. The model for scattering for a general, non-Lambertian surface with an atmosphere looks like this:

P = Pstd + trans*(rho*Ah*munot)/(1.d0-rho*Ab*sbar) + trans0*rho*(psurf-Ah*munot)

where P is the overall photometric function (the model of the data), PSTD is the pure atmospheric-scattering term, PSURF is the surface photometric function, AH*MUNOT is a Lambertian approximation to this with hemispheric albedo AH, TRANS and TRANS0 quantify transmission of surface reflected light through the atmosphere overall and with no scatterings in the atmosphere, and finally SBAR quantifies the illumination of the ground by the sky. RHO is the ratio of the sur- face albedo to the albedo assumed in the functional form of PSURF.

Parameters
phaseValue of the phase angle.
incidenceValue of the incidence angle.
emissionValue of the emission angle.
History:

2000-07-07 Randy Kirk - USGS, Flagstaff - Original code

2000-12-18 K Teal Thompson Port to Unix/ISIS;

2007-02-20 Janet Barrett - Imported from Isis2 pht_atm_functions to Isis3.

2008-11-05 Jeannie Walldren - Modified references to NumericalMethods class. Replaced Isis::PI with PI since this is in Isis namespace.

2011-12-19 Janet Barrett - Added code to estimate the shadow brightness value (transs). Also got rid of unnecessary check for identical photometric angle values between successive calls. This check should only be made in the photometric models.

Implements Isis::AtmosModel.

Definition at line 56 of file HapkeAtm2.cpp.

References _FILEINFO_, Isis::AtmosModel::AtmosWha(), Isis::AtmosModel::Ei(), Isis::AtmosModel::En(), Isis::NumericalApproximation::Evaluate(), Isis::NumericalApproximation::Extrapolate, Isis::AtmosModel::G11Prime(), Isis::AtmosModel::GenerateHahgTables(), Isis::AtmosModel::GenerateHahgTablesShadow(), Isis::NumericalAtmosApprox::OuterFunction, Isis::AtmosModel::p_atmosEstTau, Isis::AtmosModel::p_atmosHahgt0Spline, Isis::AtmosModel::p_atmosHahgtSpline, Isis::AtmosModel::p_atmosHnorm, Isis::AtmosModel::p_pstd, Isis::AtmosModel::p_sbar, Isis::AtmosModel::p_trans, Isis::AtmosModel::p_trans0, Isis::AtmosModel::p_transs, Isis::PI, Isis::NumericalApproximation::Reset(), Isis::NumericalAtmosApprox::RombergsMethod(), Isis::AtmosModel::TauOrWhaChanged(), and Isis::IException::Unknown.

◆ AtmosMunot()

double Isis::AtmosModel::AtmosMunot ( ) const
inlineinherited

Return atmospheric Munot value.

Definition at line 167 of file AtmosModel.h.

◆ AtmosNinc()

int Isis::AtmosModel::AtmosNinc ( ) const
inlineinherited

Return atmospheric Ninc value.

Definition at line 163 of file AtmosModel.h.

◆ AtmosNulneg()

bool Isis::AtmosModel::AtmosNulneg ( ) const
inlineinherited

Return atmospheric Nulneg value.

Definition at line 151 of file AtmosModel.h.

◆ AtmosTau()

double Isis::AtmosModel::AtmosTau ( ) const
inlineinherited

Return atmospheric Tau value.

Definition at line 135 of file AtmosModel.h.

Referenced by Isis::NumericalAtmosApprox::InrFunc2Bint().

◆ AtmosTauref()

double Isis::AtmosModel::AtmosTauref ( ) const
inlineinherited

Return atmospheric Tauref value.

Definition at line 147 of file AtmosModel.h.

◆ AtmosWha()

double Isis::AtmosModel::AtmosWha ( ) const
inlineinherited

Return atmospheric Wha value.

Definition at line 139 of file AtmosModel.h.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), and AtmosModelAlgorithm().

◆ CalcAtmEffect()

void Isis::AtmosModel::CalcAtmEffect ( double  pha,
double  inc,
double  ema,
double *  pstd,
double *  trans,
double *  trans0,
double *  sbar,
double *  transs 
)
inherited

Calculate the atmospheric scattering effect using photometric angle information.

Parameters
phaphase angle
incincidence angle
emaemission angle
pstdpure atmospheric-scattering term
transtransmission of surface reflected light through the atmosphere overall
trans0transmission of surface reflected light through the atmosphere with no scatterings in the atmosphere
sbarillumination of the ground by the sky
Exceptions
Isis::IException::User"Invalid photometric angles"

Definition at line 469 of file AtmosModel.cpp.

Referenced by Isis::AlbedoAtm::AlbedoAtm(), and Isis::AlbedoAtm::NormModelAlgorithm().

◆ Ei()

double Isis::AtmosModel::Ei ( double  x)
staticinherited

This routine computes the exponential integral, Ei(x).

This is defined as

\[ Ei(x) = - \int_{-x}^{\infty} \frac{e^{-t}}{t}\mathrm{d}t = \int_{-\infty}^{x} \frac{e^{t}}{t}\mathrm{d}t \]

for x > 0.

For small x, this method uses the power series,

\[ Ei(x) = \gamma + \ln x + \frac{x}{1*1!} + \frac{x^2}{2*2!} + ... \]

where \( \gamma = 0.5772156649... \) is Euler's constant, and for large x, the asymptotic series is used,

\[ Ei(x) \sim \frac{e^x}{x} (1 + \frac{1!}{x} + \frac{2!}{x^2} + ...) \]

See also
mathworld.wolfram.com/ExponentialIntegral.html
Parameters
xValue at which the exponential integral will be computed, x > 0.0
Returns
double The exponential integral for the given x.
Exceptions
Isis::IException::Programmer"Invalid arguments. Definition requires x > 0.0."
Isis::IException::Math"Power series failed to converge"
History:

1999-08-11 K Teal Thompson - Original version named pht_r8ei in Isis2.

2007-02-20 Janet Barrett - Imported to Isis3 in NumericalMethods class named r8ei().

2008-11-05 Jeannie Walldren - Renamed and modified input parameters. Added documentation.

Definition at line 223 of file AtmosModel.cpp.

References _FILEINFO_.

Referenced by Isis::Isotropic2::AtmosModelAlgorithm(), AtmosModelAlgorithm(), and Isis::Anisotropic2::AtmosModelAlgorithm().

◆ En()

double Isis::AtmosModel::En ( unsigned int  n,
double  x 
)
staticinherited

This routine evaluates the generalized exponential integral, En(x).

This is defined as

\[ E_n(x) = \int_{1}^{\infty} \frac{e^{-xt}}{t^n}\mathrm{d}t \mathrm{ for } x \in \mathbb{R }, n \in \mathbb{Z } \mathrm{such that } x \ge 0.0, n \ge 0 \]

Notice that if x = 0, then the equation is not defined for n = 0 or n = 1, as these will result in \( E_n(x) = \infty \).

This method uses the following for cases:

  • For n > 1 and x = 0.0,

    \[ E_n(0) = \frac{1}{n-1} \]

  • For n = 0 and x > 0.0,

    \[ E_0(x) = \frac{e^{-x}}{x} \]

  • For \( x >\sim 1.0 \), Lentz's continued fraction algorithm is used,

    \[ E_n(x) = e^{-x}(\frac{1}{x + n - \frac{1*n}{x + n + 2 - \frac{2(n+1)}{x + n + 4 - ...}}}) \]

  • For \( 0.0 < x <\sim 1.0 \), series representation is used,

    \[ E_n(x) = \frac{(-x)^{n-1}}{(n-1)!}[- \ln x + \psi (n)] - \sum_{m=0, m \neq n-1}^{\infty}\frac{(-x)^m}{(m-n+1)m!} \]

    where \( \psi (1) = \gamma \) and \( \psi (n) = - \gamma + \sum_{m=1}^{n-1} \frac{1}{m}\) and \( \gamma = 0.5772156649... \) is Euler's constant.

The parameter n is of type unsigned int since n cannot be negative.

The routine allows fast evaluation of En(x) to any accuracy, \(\epsilon\), within the reach of your machine's word length for floating-point numbers. The only modification required for increased accuracy is to supply Euler's constant with enough significant digits.

See also
mathworld.wolfram.com/En-Function.html
Parameters
nInteger value at which the exponential integral will be evaluated (n >= 0)
xThe exponential integral En(x) will be evaluated
Returns
double Value of the exponential integral for the give n and x.
Exceptions
Isis::IException::Programmer"Invalid arguments. Definition requires (x > 0.0 and n >=0 ) or (x >= 0.0 and n > 1)." @throws Isis::IException::Math "Continued fraction failed to converge" @throws Isis::IException::Math "Series representation failed to converge"
History:

1999-08-10 K Teal Thompson - Original version in named pht_r8expint Isis2.

2007-02-20 Janet Barrett - Imported to Isis3 in NumericalMethods class named r8expint().

2008-11-05 Jeannie Walldren - Renamed and modified input parameters. Added documentation.

Definition at line 364 of file AtmosModel.cpp.

References _FILEINFO_.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), Isis::Isotropic2::AtmosModelAlgorithm(), AtmosModelAlgorithm(), Isis::Anisotropic1::AtmosModelAlgorithm(), and Isis::Anisotropic2::AtmosModelAlgorithm().

◆ G11Prime()

double Isis::AtmosModel::G11Prime ( double  tau)
staticinherited

Perform Chandra and Van de Hulst's series approximation for the g'11 function needed in second order scattering theory.

Parameters
taunormal optical depth of atmosphere
Returns
double Value of the g'11 function evaluated at the given tau
History:

1998-12-21 Randy Kirk, USGS - Flagstaff - Original code in Isis2 pht_am_functions

1999-03-12 K Teal Thompson - Port to Unix/ISIS; declare vars; cadd implicit none.

2007-02-20 Janet Barrett - Imported from Isis2 to Isis3 in NumericalMethods class.

2008-11-05 Jeannie Walldren - Moved this method from NumericalMethods class.

Definition at line 147 of file AtmosModel.cpp.

Referenced by Isis::Isotropic2::AtmosModelAlgorithm(), AtmosModelAlgorithm(), and Isis::Anisotropic2::AtmosModelAlgorithm().

◆ GenerateAhTable()

void Isis::AtmosModel::GenerateAhTable ( )
inherited

This method computes the values of the atmospheric Ah table and sets the properties of the atmospheric Ah spline.

It obtains the hemispheric albedo by integrating the photometric function times mu over mu = 0 to 1 and then over phi = 0 to 2*pi to calculate the hemispheric reflectance Ah needed for the photometric model with atmosphere. The Trapezoid rule is applied to the table of Ah to obtain bihemispheric albedo Ab. The parameter p_atmosAtmSwitch is set to 0 to integrate Ah.

History:

1998-12-21 Randy Kirk - USGS, Flagstaff - Original specs

1999-01-15 K Teal Thompson - Original code

2000-12-29 Randy Kirk - Modified /hide_inc/ so phi gets passed, etc. moved factors to outside integration

2006-05-30 Randy Kirk - Added code to report the directional hemispheric albedo that relates to thermal balance for simulating Themis images

2007-02-20 Janet Barrett - Imported to Isis3 from Isis2 pht_get_ah_table

2008-11-05 Jeannie Walldren - Modified references to NumericalMethods and replaced r8qromb() with NumericalAtmosApprox::RombergsMethod(). Changed name from PhtGetAhTable().

Definition at line 531 of file AtmosModel.cpp.

References Isis::PI, Isis::NumericalAtmosApprox::RombergsMethod(), and Isis::IString::UpCase().

Referenced by Isis::AlbedoAtm::AlbedoAtm().

◆ GenerateHahgTables()

void Isis::AtmosModel::GenerateHahgTables ( )
inherited

This method computes the values of the atmospheric Hahg and Hahg0 tables and sets the properties of the atmospheric Hahg and Hahg0 splines.

It integrates functions involving the single particle phase function (assumed to be Hapke Henyey-Greenstein) over mu = 0 to 1 and then over phi = 0 to 2*pi to calculate the corrections needed for the anisotropic photometric model with atmosphere. The Trapezoid rule is applied to the table of Ah to obtain bihemispheric albedo Ab. The parameter p_atmosAtmSwitch is set to 1, 2, 3 to evaluate the 3 required integrals.

History:

1998-12-21 Randy Kirk - USGS, Flagstaff - Original specs

1999-01-15 K Teal Thompson - Original code

2006-07-07 Randy Kirk - Modify get_ah_table to get other integrals

2000-12-29 Randy Kirk - Modified /hide_inc/ so phi gets passed, etc. moved factors to outside integration

2007-02-20 Janet Barrett - Imported to Isis3 from Isis2 pht_get_ah_table

2008-11-05 Jeannie Walldren - Modified references to NumericalMethods and replaced r8qromb() with NumericalAtmosApprox::RombergsMethod(). Changed name from GetHahgTables().

Definition at line 631 of file AtmosModel.cpp.

References Isis::PI, Isis::NumericalApproximation::Reset(), and Isis::NumericalAtmosApprox::RombergsMethod().

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), and AtmosModelAlgorithm().

◆ GenerateHahgTablesShadow()

void Isis::AtmosModel::GenerateHahgTablesShadow ( )
inherited

This method is a modified version of the GenerateHahgTables method and is used solely for shadow modeling.

Unlike the GenerateHahgTables method, this method does not tabulate the first or third integrals. It only evaluates the middle integral that corrects the sbar variable (illumination of the ground by the sky).

History:

1998-12-21 Randy Kirk - USGS, Flagstaff - Original specs

1999-01-15 K Teal Thompson - Original code

2006-07-07 Randy Kirk - Modify get_ah_table to get other integrals

2000-09-15 Randy Kirk - Delete first, last integrals and table of second integral

2000-12-29 Randy Kirk - Modified /hide_inc/ so phi gets passed, etc. moved factors to outside integration

2011-12-16 Janet Barrett - Imported to Isis3 from Isis2 get_hahgsb_shad routine in the not_shadow_tau set of routines

Definition at line 720 of file AtmosModel.cpp.

References Isis::PI, Isis::NumericalApproximation::Reset(), and Isis::NumericalAtmosApprox::RombergsMethod().

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), and AtmosModelAlgorithm().

◆ SetAtmosAtmSwitch()

void Isis::AtmosModel::SetAtmosAtmSwitch ( const int  atmswitch)
inherited

Set the switch that controls the function that will be integrated.

This method is only used for testing the methods in this class. This parameter is limited to the values 0, 1, 2, and 3.

Parameters
atmswitchInternal atmospheric function parameter, there is no default
Exceptions
Isis::IException::User"Invalid value of atmospheric atmswitch"

Definition at line 773 of file AtmosModel.cpp.

References _FILEINFO_.

◆ SetAtmosBha()

void Isis::AtmosModel::SetAtmosBha ( const double  bha)
inherited

Set the Anisotropic Atmospheric function parameter.

This is the coefficient of the single particle Legendre phase function. This parameter is limited to values that are >=-1 and <=1.

Parameters
bhaAnisotropic atmospheric function parameter, default is 0.85
Exceptions
Isis::IException::User"Invalid value of atmospheric bha"

Definition at line 793 of file AtmosModel.cpp.

References _FILEINFO_.

◆ SetAtmosEstTau()

void Isis::AtmosModel::SetAtmosEstTau ( const string  esttau)
inherited

Estimate the optical depth tau using shadows.

Author
Janet Barrett (12/16/2011)
Parameters
esttautrue/false

Definition at line 996 of file AtmosModel.cpp.

References _FILEINFO_, and Isis::IString::UpCase().

◆ SetAtmosHga()

void Isis::AtmosModel::SetAtmosHga ( const double  hga)
inherited

Set the Hapke atmospheric function parameter.

This is the coefficient of the single particle Henyey-Greenstein phase function. This parameter is limited to values that are >-1 and <1.

Parameters
hgaHapke atmospheric function parameter, default is 0.68
Exceptions
Isis::IException::User"Invalid value of atmospheric hga"

Definition at line 814 of file AtmosModel.cpp.

References _FILEINFO_.

◆ SetAtmosHnorm()

void Isis::AtmosModel::SetAtmosHnorm ( const double  hnorm)
inherited

Set the Atmospheric function parameter.

This is the atmospheric shell thickness normalized to the planet radius and is used to modify angles to get more accurate path lengths near the terminator (ratio of scale height to the planetary radius). This parameter is limited to values that are >=0.

Parameters
hnormAtmospheric function parameter, default is 0.003

Definition at line 962 of file AtmosModel.cpp.

References _FILEINFO_, and Isis::toString().

◆ SetAtmosInc()

void Isis::AtmosModel::SetAtmosInc ( const double  inc)
inherited

Set the incidence angle.

This method is only used for testing the methods in this class. This parameter is limited to values >=0 and <=90.

Parameters
incInternal atmospheric function parameter, there is no default
Exceptions
Isis::IException::User"Invalid value of atmospheric inc"

Definition at line 833 of file AtmosModel.cpp.

References _FILEINFO_, and Isis::PI.

◆ SetAtmosIord()

void Isis::AtmosModel::SetAtmosIord ( const string  offset)
inherited

Set additive offset in fit.

Author
Sharmila Prasad (8/18/2011)
Parameters
offsettrue/false

Definition at line 977 of file AtmosModel.cpp.

References _FILEINFO_, and Isis::IString::UpCase().

◆ SetAtmosNulneg()

void Isis::AtmosModel::SetAtmosNulneg ( const string  nulneg)
inherited

Set the Atmospheric function parameter.

This determines if negative values after removal of atmospheric effects will be set to NULL. This parameter is limited to values of YES or NO.

Parameters
nulnegAtmospheric function parameter, default is NO
Exceptions
Isis::IException::User"Invalid value of atmospheric nulneg"

Definition at line 855 of file AtmosModel.cpp.

References _FILEINFO_, and Isis::IString::UpCase().

◆ SetAtmosPhi()

void Isis::AtmosModel::SetAtmosPhi ( const double  phi)
inherited

Set the azimuth angle.

This method is only used for testing the methods in this class. This parameter is limited to values >=0 and <=360.

Parameters
phiInternal atmospheric function parameter, there is no default
Exceptions
Isis::IException::User"Invalid value of atmospheric phi"

Definition at line 878 of file AtmosModel.cpp.

References _FILEINFO_, and Isis::PI.

◆ SetAtmosTau()

void Isis::AtmosModel::SetAtmosTau ( const double  tau)
inherited

Set the Atmospheric function parameter.

This specifies the normal optical depth of the atmosphere. This parameter is limited to values that are >=0.

Parameters
tauAtmospheric function parameter, default is 0.28
Exceptions
Isis::IException::User"Invalid value of atmospheric tau"

Definition at line 897 of file AtmosModel.cpp.

References _FILEINFO_.

◆ SetAtmosTauref()

void Isis::AtmosModel::SetAtmosTauref ( const double  tauref)
inherited

Set the Atmospheric function parameter.

This specifies the reference optical depth of the atmosphere to which the image will be normalized. This parameter is limited to values that are >=0.

Parameters
taurefAtmospheric function parameter, default is 0.0
Exceptions
Isis::IException::User"Invalid value of atmospheric tauref"

Definition at line 916 of file AtmosModel.cpp.

References _FILEINFO_.

◆ SetAtmosWha()

void Isis::AtmosModel::SetAtmosWha ( const double  wha)
inherited

Set the Atmospheric function parameter.

This is the single- scattering albedo of atmospheric particles. This parameter is limited to values that are >0 and <=1.

Parameters
whaAtmospheric function parameter, default is 0.95
Exceptions
Isis::IException::User"Invalid value of atmospheric wha"

Definition at line 935 of file AtmosModel.cpp.

References _FILEINFO_.

◆ SetStandardConditions()

void Isis::AtmosModel::SetStandardConditions ( bool  standard)
virtualinherited

Used to calculate atmosphere at standard conditions.

Definition at line 492 of file AtmosModel.cpp.

Referenced by Isis::AlbedoAtm::AlbedoAtm().

◆ TauOrWhaChanged()

bool Isis::AtmosModel::TauOrWhaChanged ( ) const
protectedinherited

Checks whether tau or wha have changed.

Returns
bool True if either have changed.

Definition at line 948 of file AtmosModel.cpp.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), Isis::Isotropic2::AtmosModelAlgorithm(), AtmosModelAlgorithm(), Isis::Anisotropic1::AtmosModelAlgorithm(), and Isis::Anisotropic2::AtmosModelAlgorithm().

Member Data Documentation

◆ p_atmosAddOffset

bool Isis::AtmosModel::p_atmosAddOffset
protectedinherited

Allow additive offset in fit.

Definition at line 284 of file AtmosModel.h.

◆ p_atmosAhSpline

NumericalApproximation Isis::AtmosModel::p_atmosAhSpline
protectedinherited

Spline object for the atmospheric Ah Table. Properties are set in GenerateAhTable().

Definition at line 299 of file AtmosModel.h.

◆ p_atmosEstTau

bool Isis::AtmosModel::p_atmosEstTau
protectedinherited

Estimate optical depth tau using shadows.

Definition at line 285 of file AtmosModel.h.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), and AtmosModelAlgorithm().

◆ p_atmosHahgt0Spline

NumericalApproximation Isis::AtmosModel::p_atmosHahgt0Spline
protectedinherited

Spline object for the atmospheric Hahg0 Table. Properties are set in GenerateHahgTables().

Definition at line 303 of file AtmosModel.h.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), and AtmosModelAlgorithm().

◆ p_atmosHahgtSpline

NumericalApproximation Isis::AtmosModel::p_atmosHahgtSpline
protectedinherited

Spline object for the atmospheric Hahg Table. Properties are set in GenerateHahgTables().

Definition at line 301 of file AtmosModel.h.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), and AtmosModelAlgorithm().

◆ p_atmosHnorm

double Isis::AtmosModel::p_atmosHnorm
protectedinherited

◆ p_pstd

double Isis::AtmosModel::p_pstd
protectedinherited

◆ p_sbar

double Isis::AtmosModel::p_sbar
protectedinherited

◆ p_trans

double Isis::AtmosModel::p_trans
protectedinherited

Transmission of surface reflected light through the atmosphere overall.

Definition at line 275 of file AtmosModel.h.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), Isis::Isotropic2::AtmosModelAlgorithm(), AtmosModelAlgorithm(), Isis::Anisotropic1::AtmosModelAlgorithm(), and Isis::Anisotropic2::AtmosModelAlgorithm().

◆ p_trans0

double Isis::AtmosModel::p_trans0
protectedinherited

Transmission of surface reflected light through the atmosphere with no scatterings in the atmosphere.

Definition at line 276 of file AtmosModel.h.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), Isis::Isotropic2::AtmosModelAlgorithm(), AtmosModelAlgorithm(), Isis::Anisotropic1::AtmosModelAlgorithm(), and Isis::Anisotropic2::AtmosModelAlgorithm().

◆ p_transs

double Isis::AtmosModel::p_transs
protectedinherited

Transmission of light that must be subtracted from the flat surface model to get the shadow model.

Definition at line 277 of file AtmosModel.h.

Referenced by Isis::HapkeAtm1::AtmosModelAlgorithm(), Isis::Isotropic2::AtmosModelAlgorithm(), AtmosModelAlgorithm(), Isis::Anisotropic1::AtmosModelAlgorithm(), and Isis::Anisotropic2::AtmosModelAlgorithm().


The documentation for this class was generated from the following files: