*/
#define PROJ_PARMS__ \
void *en, *en2; \
double phi1; \
double phit; \
double n, kRF; \
int czech;
#define PROJ_LIB__
#include <lib_proj.h>
PROJ_HEAD(kocc, "Krovak Oblique Conformal Conic")
"\n\tConic, Sph&Ell\n\tlat_1= lat_t=";
FORWARD(e_forward); /* ellipsoid */
double rho, theta;
lp = proj_translate(proj_gauss(lp, P->en), P->en2);
rho = P->kRF / pow(tan(.5 * lp.phi + FORTPI), P->n);
theta = P->n * lp.lam;
if (P->czech) { /* Czech grid mode */
xy.x = rho * cos(theta);
xy.y = - rho * sin(theta);
} else { /* proper math mode */
xy.x = rho * sin(theta);
xy.y = - rho * cos(theta);
}
return (xy);
}
INVERSE(e_inverse); /* ellipsoid */
double x, y, rho, theta;
if (P->czech) {
double cosc0, sinc0; \
double R2; \
void *en;
#define PROJ_LIB__
#include <lib_proj.h>
PROJ_HEAD(sterea, "Oblique Stereographic Alternative")
"\n\tAzimuthal, Sph&Ell";
# define DEL_TOL 1.e-14
# define MAX_ITER 10
FORWARD(e_forward); /* ellipsoid */
double cosc, sinc, cosl, k;
lp = proj_gauss(lp, P->en);
sinc = sin(lp.phi);
cosc = cos(lp.phi);
cosl = cos(lp.lam);
k = P->k0 * P->R2 / (1. + P->sinc0 * sinc + P->cosc0 * cosc * cosl);
xy.x = k * cosc * sin(lp.lam);
xy.y = k * (P->cosc0 * sinc - P->sinc0 * cosc * cosl);
return (xy);
}
INVERSE(e_inverse); /* ellipsoid */
double rho, c, sinc, cosc;
xy.x /= P->k0;
xy.y /= P->k0;
if((rho = hypot(xy.x, xy.y))) {
c = 2. * atan2(rho, P->R2);
