1 /* 2 ** libproj -- library of cartographic projections 3 ** 4 ** Copyright (c) 2003, 2006 Gerald I. Evenden 5 */ 6 static const char 7 LIBPROJ_ID[] = "Id"; 8 /* 9 ** Permission is hereby granted, free of charge, to any person obtaining 10 ** a copy of this software and associated documentation files (the 11 ** "Software"), to deal in the Software without restriction, including 12 ** without limitation the rights to use, copy, modify, merge, publish, 13 ** distribute, sublicense, and/or sell copies of the Software, and to 14 ** permit persons to whom the Software is furnished to do so, subject to 15 ** the following conditions: 16 ** 17 ** The above copyright notice and this permission notice shall be 18 ** included in all copies or substantial portions of the Software. 19 ** 20 ** THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, 21 ** EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF 22 ** MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. 23 ** IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY 24 ** CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, 25 ** TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE 26 ** SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. 27 */ 28 #define PROJ_PARMS__ \ 29 double A, B; 30 #define PROJ_LIB__ 31 # include <lib_proj.h> 32 PROJ_HEAD(putp5, "Putnins P5") "\n\tPCyl., Sph."; 33 PROJ_HEAD(putp5p, "Putnins P5'") "\n\tPCyl., Sph."; 34 #define C 1.01346 35 #define D 1.2158542 36 FORWARD(s_forward); /* spheroid */ 37 xy.x = C * lp.lam * (P->A - P->B * sqrt(1. + D * lp.phi * lp.phi)); 38 xy.y = C * lp.phi; 39 return (xy); 40 } 41 INVERSE(s_inverse); /* spheroid */ 42 lp.phi = xy.y / C; 43 lp.lam = xy.x / (C * (P->A - P->B * sqrt(1. + D * lp.phi * lp.phi))); 44 return (lp); 45 } 46 FREEUP; if (P) free(P); } 47 static PROJ * 48 setup(PROJ *P) { 49 P->es = 0.; P->inv = s_inverse; P->fwd = s_forward; 50 return P; 51 } 52 ENTRY0(putp5) P->A = 2.; P->B = 1.; ENDENTRY(setup(P)) 53 ENTRY0(putp5p) P->A = 1.5; P->B = 0.5; ENDENTRY(setup(P)) 54 /* 55 ** Log: proj_putp5.c 56 ** Revision 3.1 2006/01/11 01:38:18 gie 57 ** Initial 58 ** 59 */ 60