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 CM 0.77777777777777777777777777778 29 #define CN 0.27777777777777777777777777778 30 #define CX 3.6 31 #define CY 1.285714285714285714285714285 32 #define PROJ_LIB__ 33 #include <lib_proj.h> 34 PROJ_HEAD(wag9, "Wagner IX") "\n\tMod. Azim, no inv."; 35 FORWARD(s_forward); /* sphere */ 36 double cosdel, cosalph; 37 (void) P; /* avoid warning */ 38 39 if ((xy.y = proj_acos(cosdel = cos(lp.phi *= CM) * cos(CN * lp.lam))) == 0.) 40 xy.x = xy.y = 0; 41 else { 42 cosalph = sin(lp.phi) / sqrt(fabs(1. - cosdel * cosdel)); 43 xy.x = CX * xy.y * sqrt(fabs(1. - cosalph * cosalph)); 44 if (lp.lam < 0.) 45 xy.x = -xy.x; 46 xy.y *= CY * cosalph; 47 } 48 return (xy); 49 } 50 FREEUP; if (P) free(P); } 51 ENTRY0(wag9) P->fwd = s_forward; P->inv = 0; P->es = 0.; ENDENTRY(P) 52 /* 53 ** Log: proj_wag9.c 54 ** Revision 3.2 2006/01/19 03:28:27 gie 55 ** remove 'del' because unused 56 ** 57 ** Revision 3.1 2006/01/11 01:38:18 gie 58 ** Initial 59 ** 60 */ 61