gs48 problems

[proj4-bot]

Reported by fwillett on 10 May 2012 17:36 UTC
System:
OS - Windows 7 Ultimate, SP1, 64 bit
CPU - Intel Core2 Quad CPU @ 2.83 GHz
Memory - 1.86 GB
Running in Oracle VirtualBox 4.16 r74713

Installation:
proj version - Rel. 4.8.0, 6 March 2012
Microsoft Visual Studio 2010 Professional Version 10.0.40219.1 SP1Rel
Build commands:
"C:\Program Files (x86)\Microsoft Visual Studio 10.0\VC\bin\vcvars32.bat"
nmake /f makefile.vc
nmake /f makefile.vc install-all

Note: This builds a 32 bit version of the programs and libraries.

In the following the OS prompt is "$" and the proj prompt is "%"

gs48 - 1) The origin of the projection is at long/lat (-96,-39) and it should
be at (-96,+39)
Modified pj_mod_ster.c to change the sign on phi0 for the gs48 projection
// P->phi0 = DEG_TO_RAD * -39.;
P->phi0 = DEG_TO_RAD * 39.;

$ ./proj.exe +proj=gs48
%-96 39
-0.00 16805051.65

$ ./proj.exe +proj=gs48 -I
%0 0
168dE 39dS
%1 1
95d59'59.958"W 39d11'17.725"S

   2) Inverse problem at the origin with modified code. 
      NOTE: This problem is evident in the unmodified code in the example above.

$ proj +proj=gs48 -I
%0 0
168dE 39dN
%1 1
95d59'59.958"W 39d11'17.79"N

Migrated-From: https://trac.osgeo.org/proj/ticket/162

[kbevers]

Fixed with 0eccceb and 282d6ab