pj_lask.cpp

projapi是一个关于GIS行业投影转换的程序库

#include "stdafx.h"

#ifndef lint
static const char SCCSID[]="@(#)PJ_lask.c	4.1	94/02/15	GIE	REL";
#endif
#define PROJ_PARMS__
#define PJ_LIB__
#include	"projects.h"

PROJ_HEAD(lask, "Laskowski") "\n\tMisc Sph, no inv.";
#define	a10	 0.975534
#define	a12	-0.119161
#define	a32	-0.0143059
#define	a14	-0.0547009
#define	b01	 1.00384
#define	b21	 0.0802894
#define	b03	 0.0998909
#define	b41	 0.000199025
#define	b23	-0.0285500
#define	b05	-0.0491032
FORWARD(s_forward); /* sphere */
	double l2, p2;

	l2 = lp.lam * lp.lam;
	p2 = lp.phi * lp.phi;
	xy.x = lp.lam * (a10 + p2 * (a12 + l2 * a32 + p2 * a14));
	xy.y = lp.phi * (b01 + l2 * (b21 + p2 * b23 + l2 * b41) +
		p2 * (b03 + p2 * b05));
	return (xy);
}
FREEUP; if (P) pj_dalloc(P); }
ENTRY0(lask) P->fwd = s_forward; P->inv = 0; P->es = 0.; ENDENTRY(P)