00001 #include "gis.h"
00002 #include "pi.h"
00003
00004 extern double sin(), sqrt();
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00053 double
00054 G_meridional_radius_of_curvature (lon, a, e2)
00055 double lon, a, e2;
00056 {
00057 double x;
00058 double s;
00059
00060 s = sin (Radians(lon));
00061 x = 1 - e2 * s*s;
00062
00063 return a * (1 - e2) / (x * sqrt(x));
00064 }
00065
00066
00067
00084 double
00085 G_transverse_radius_of_curvature (lon, a, e2)
00086 double lon, a, e2;
00087 {
00088 double x;
00089 double s;
00090
00091 s = sin (Radians(lon));
00092 x = 1 - e2 * s*s;
00093
00094 return a / sqrt(x);
00095 }
00096
00097
00114 double
00115 G_radius_of_conformal_tangent_sphere (lon, a, e2)
00116 double lon, a, e2;
00117 {
00118 double x;
00119 double s;
00120
00121 s = sin (Radians(lon));
00122 x = 1 - e2 * s*s;
00123
00124 return a * sqrt (1 - e2) / x;
00125 }