Go to the documentation of this file.
72 double xx,yy,xys,c1,c2,c3;
86 xys = xx * xx + yy * yy;
87 c1 = -
fabs(yy) * (1.0 + xys);
88 c2 = c1 - 2.0 * yy * yy + xx * xx;
89 c3 = -2.0 * c1 + 1.0 + 2.0 * yy * yy + xys * xys;
90 d = yy * yy / c3 + (2.0 * c2 * c2 * c2 / c3 / c3 / c3 - 9.0 * c1 * c2 / c3 /c3)
92 a1 = (c1 - c2 * c2 / 3.0 / c3) / c3;
93 m1 = 2.0 * sqrt( -a1 / 3.0);
94 con = ((3.0 * d) / a1) / m1;
102 th1 = acos(con) / 3.0;
104 *
lat = (-m1 * cos(th1 +
PI / 3.0) - c2 / 3.0 / c3) *
PI;
106 *
lat = -(-m1 * cos(th1 +
PI / 3.0) - c2 / 3.0 / c3) *
PI;
114 (xx * xx - yy * yy) + xys * xys)) / 2.0 / xx);
void gctp_print_title(const char *proj_name)
void gctp_print_cenlon(double A)
double adjust_lon(double x)
long vandginv(double x, double y, double *lon, double *lat)
void gctp_print_offsetp(double A, double B)
long vandginvint(double r, double center_long, double false_east, double false_north)
void gctp_print_radius(double radius)