[gdal-dev] gdal not recognizing ellpsoid properly in esri grid format

[Dominik Schneider]

Hi,
I was sent an esri grid that uses a spherical ellpsoid with +a and +b
parameters but this not being read properly. It is being assigned
+ellps=WGS84 when read with raster::raster() in R 3.2.3 but also when using
gdalinfo on the command line.  I'm using GDAL 1.11.3 from kyngchaos.com on
OS X 10.11.3

Not sure if this is a known issue or subsequently fixed in later versions
so I wanted to report it here. I didn't see a ticket that addresses this in
particular.

> gdalinfo dem
gdalinfo dem
Driver: AIG/Arc/Info Binary Grid
Files: dem
       dem/dblbnd.adf
       dem/hdr.adf
       dem/metadata.xml
       dem/prj.adf
       dem/sta.adf
       dem/w001001.adf
       dem/w001001x.adf
Size is 317, 263
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9108"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",33],
    PARAMETER["standard_parallel_2",45],
    PARAMETER["latitude_of_origin",39.00000381472222],
    PARAMETER["central_meridian",-107],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["METERS",1]]
Origin = (-633999.740890799090266,525999.667091785930097)
Pixel Size = (4000.000000000000000,-4000.000000000000000)
Corner Coordinates:
Upper Left  ( -633999.741,  525999.667) (114d52' 9.29"W, 43d30'32.73"N)
Lower Left  ( -633999.741, -526000.333) (113d52'56.93"W, 34d 1'21.79"N)
Upper Right (  634000.259,  525999.667) ( 99d 7'50.68"W, 43d30'32.73"N)
Lower Right (  634000.259, -526000.333) (100d 7' 3.05"W, 34d 1'21.79"N)
Center      (   0.2591092,  -0.3329082) (106d59'59.99"W, 39d 0' 0.00"N)
Band 1 Block=256x16 Type=Float32, ColorInterp=Undefined
  Min=274.508 Max=3824.624
  NoData Value=-3.4028234663852886e+38


prj.adf reads:
Projection    LAMBERT
Units         METERS
Zunits        NO
Xshift        0.0
Yshift        0.0
Parameters    6370000.0  6370000.0
  33  0  0.0 /* 1st standard parallel
  45  0  0.0 /* 2nd standard parallel
-107  0  0.0 /* central meridian
  39  0  0.013733 /* latitude of projection's origin
0.0 /* false easting (meters)
0.0 /* false northing (meters)

The file is here:
https://www.odrive.com/s/4516b311-380e-4d94-8912-16f39c7d0259-56c21076
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20160215/050d8cd7/attachment.html>