[Gdal-dev] NAD83 Fipszone 903 to WGS84 Error

[Eric Lindahl]

Hi,

I attempted to fix the PRJ.ADF import that didn't recognize the FIPSZONE
key word.


691,692c694,704
         int nZone = (int) OSR_GDV( papszPrj, "zone", 0.0 );
         nZone = ESRIToUSGSZone( nZone );
---
         double def = -10000.0;
         double dZone = OSR_GDV( papszPrj, "zone", def );
         int nZone = 0;
         if( dZone == def ){
           dZone = OSR_GDV( papszPrj, "Fipszone", def );
           nZone = (dZone==def)?0:(int)dZone;
         }else{
           nZone = (int)dZone;
           nZone = ESRIToUSGSZone( nZone );
         }


After I did this, my gdalinfo command output this description with
the *wrong* corner coordinates:
Driver: AIG/Arc/Info Binary Grid
Size is 16771, 13439
Coordinate System is:
LOCAL_CS["State Plane Zone 2500 / NAD83",
    UNIT["Foot (International)",0.3048]]
Origin = (1599823.500000,411154.500000)
Pixel Size = (3.00000000,-3.00000000)
Corner Coordinates:
Upper Left  ( 1599823.500,  411154.500) 
Lower Left  ( 1599823.500,  370837.500) 
Upper Right ( 1650136.500,  411154.500) 
Lower Right ( 1650136.500,  370837.500) 
Center      ( 1624980.000,  390996.000) 
Band 1 Block=256x4 Type=Int16, ColorInterp=Undefined
  Min=-10.000 Max=384.000 
  NoData Value=-32768
Driver: AIG/Arc/Info Binary Grid
Size is 16771, 13439
Coordinate System is:
PROJCS["NAD83 / Florida North",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree (supplier to define representation)",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",30.75],
    PARAMETER["standard_parallel_2",29.58333333333333],
    PARAMETER["latitude_of_origin",29],
    PARAMETER["central_meridian",-84.5],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",0],
    UNIT["Foot (International)",0.3048],
    AUTHORITY["EPSG","26960"]]
Origin = (1599823.500000,411154.500000)
Pixel Size = (3.00000000,-3.00000000)
Corner Coordinates:
Upper Left  ( 1599823.500,  411154.500) ( 81d20'17.01"W, 30d 5'32.78"N)
Lower Left  ( 1599823.500,  370837.500) ( 81d20'29.73"W, 29d58'53.84"N)
Upper Right ( 1650136.500,  411154.500) ( 81d10'44.52"W, 30d 5'18.62"N)
Lower Right ( 1650136.500,  370837.500) ( 81d10'57.87"W, 29d58'39.70"N)
Center      ( 1624980.000,  390996.000) ( 81d15'37.28"W, 30d 2'6.32"N)
Band 1 Block=256x4 Type=Int16, ColorInterp=Undefined
  Min=-10.000 Max=384.000 
  NoData Value=-32768


The PROJ4 format for the source ends up being:
+proj=lcc +lat_1=30.75 +lat_2=29.58333333333333 +lat_0=29 +lon_0=-84.5 
+x_0=182880 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=ft +no_defs


But if the coordinates are calculated by hand:
Converting Upper Left coordinates to meters
1599823.5 * 1200 / 3937 = 487627.178 meters
411154.5 * 1200 / 3937 = 125320.14224 meters

Using these numbers as input at
http://www.ngs.noaa.gov/cgi-bin/spc_getgp.prl
===========================================================
                 North(Meters)     East(Meters)     Datum     Zone
 INPUT =  125320.14224      487627.178       NAD83     0903
 ===========================================================
  LATITUDE         LONGITUDE         AREA    CONVERGENCE   SCALE FACTOR
 DD MM SS.sssss   DDD MM SS.sssss           DD MM SS.ss      
 --------------        ---------------        ----     -----------          ----
--------
 30 07 31.37149   085 39 58.31866   FL N    -0 35  9.76        0.9999487
 
 USING TOPOZONE AND LOOKING AT THE ADF FILE DISPLAYED IN GEOMATICA,
I WOULD SAY THAT THIS IS YOUR EXACT UPPER LEFT HAND CORNER

 North 30 07 31.37149   West 085 39 58.31866
    Or in decimal ......
 30.1254°N     85.6662°W    
 
And then converting Lower Right coordinates to meters
1650136.500 * 1200 / 3937 = 502962.6111 meters
370837.500 * 1200 / 3937 =  113031.496 meters
 
===========================================================
               North(Meters)     East(Meters)     Datum     Zone
 INPUT =  113031.496        502962.6111      NAD83     0903
===========================================================
 LATITUDE             LONGITUDE         AREA    CONVERGENCE   SCALE FACTOR
 DD MM SS.sssss   DDD MM SS.sssss           DD MM SS.ss      
 --------------         ---------------        ----      -----------          --
----------
 30 00 57.03675   085 30 21.35891   FL N    -0 30 19.83          0.9999519