6 LunarLambertMcEwen::LunarLambertMcEwen(Pvl &pvl) : PhotoModel(pvl) {
9 p_photoM3 = -0.00000146;
11 double c30 = cos(30.0 *
Isis::PI / 180.0);
12 double xl30 = 1.0 + p_photoM1 * 30.0 + p_photoM2 * pow(30.0, 2.0) +
13 p_photoM3 * pow(30.0, 3);
14 p_photoR30 = 2.0 * xl30 * c30 / (1.0 + c30) + (1.0 - xl30) * c30;
17 double LunarLambertMcEwen::PhotoModelAlgorithm(
double phase,
double incidence,
19 static double pht_moonpr;
25 static double old_phase = -9999;
26 static double old_incidence = -9999;
27 static double old_emission= -9999;
29 if (old_phase == phase && old_incidence == incidence && old_emission == emission) {
34 old_incidence = incidence;
35 old_emission = emission;
37 incrad = incidence *
Isis::PI / 180.0;
38 emarad = emission *
Isis::PI / 180.0;
42 double xl = 1.0 + p_photoM1 * phase + p_photoM2 * pow(phase, 2) +
43 p_photoM3 * pow(phase, 3);
44 double r = 2.0 * xl * munot / (mu + munot) + (1.0 - xl) * munot;
50 pht_moonpr = p_photoR30 / r;
const double PI
The mathematical constant PI.
Container for cube-like labels.
Moonpr photometric model Computes normalized albedo for the Moon, normalized to 0 degrees emission an...
Namespace for ISIS/Bullet specific routines.