diff options
| -rw-r--r-- | src/Panorama.H | 3 | ||||
| -rw-r--r-- | src/Panorama.cxx | 19 |
2 files changed, 18 insertions, 4 deletions
diff --git a/src/Panorama.H b/src/Panorama.H index aea48b8..bfa70d5 100644 --- a/src/Panorama.H +++ b/src/Panorama.H @@ -35,6 +35,7 @@ class Panorama { double distance(double phi, double lam); double alpha(const Hill *m); double nick(const Hill *m); + double refraction(const Hill *m); double comp_center_angle(double alph_a, double alph_b, double d1, double d2); double comp_scale(double alph_a, double alph_b, double d1, double d2); int get_matrix(double m[]); @@ -73,7 +74,7 @@ class Panorama { double get_view_long(); double get_view_height(); double get_earth_radius(double latitude); - double get_real_distance(Hill *m); + double get_real_distance(const Hill *m); int comp_params(Hills *h); ProjectionLSQ::Projection_t get_projection(); void set_projection(ProjectionLSQ::Projection_t p); diff --git a/src/Panorama.cxx b/src/Panorama.cxx index 8bab3e4..5f59103 100644 --- a/src/Panorama.cxx +++ b/src/Panorama.cxx @@ -432,13 +432,26 @@ Panorama::alpha(const Hill *m) { } double +Panorama::refraction(const Hill *m) { + double a, b, c, alpha = 6.5, T0 = 0.0; + + a = 2.9e-4 * exp (-view_height / 10000.0) / (1.0 + 2.9 * T0 / 760.0); + b = 2.9 * alpha / (760.0 * (1.0 + 2.9 * T0 / 760.0)); + c = a * (b - 1.0 / 10.0); + + return c * get_real_distance(m) / (2000.0 * (1.0 + a)); +} +double Panorama::nick(const Hill *m) { - double b, c; + double b, c, theta = refraction(m); b = m->height + get_earth_radius(m->phi); c = view_height + get_earth_radius(view_phi); - return atan((cos(m->dist) * b - c) / (sin(m->dist) * b)); +if (get_real_distance(m) < 100000) +fprintf(stderr, "=== %s, %g, %g\n", m->name, get_real_distance(m), theta / deg2rad); + + return atan((cos(m->dist) * b - c) / (sin(m->dist) * b)) - theta; } // return local distance to center of WGS84 ellipsoid @@ -454,7 +467,7 @@ Panorama::get_earth_radius(double phi) { } double -Panorama::get_real_distance(Hill *m) { +Panorama::get_real_distance(const Hill *m) { double a, b, c, gam; a = view_height + get_earth_radius(view_phi); |