<div dir="ltr"><div><div class="gmail_signature"><div dir="ltr"><div><div><div>Hi,</div></div></div><div>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 <a href="http://kyngchaos.com">kyngchaos.com</a> on OS X 10.11.3</div><div><br></div><div>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.</div><div><br></div><div>> gdalinfo dem</div><div><div>gdalinfo dem</div><div>Driver: AIG/Arc/Info Binary Grid</div><div>Files: dem</div><div>       dem/dblbnd.adf</div><div>       dem/hdr.adf</div><div>       dem/metadata.xml</div><div>       dem/prj.adf</div><div>       dem/sta.adf</div><div>       dem/w001001.adf</div><div>       dem/w001001x.adf</div><div>Size is 317, 263</div><div>Coordinate System is:</div><div>PROJCS["unnamed",</div><div>    GEOGCS["WGS 84",</div><div>        DATUM["WGS_1984",</div><div>            SPHEROID["WGS 84",6378137,298.257223563,</div><div>                AUTHORITY["EPSG","7030"]],</div><div>            TOWGS84[0,0,0,0,0,0,0],</div><div>            AUTHORITY["EPSG","6326"]],</div><div>        PRIMEM["Greenwich",0,</div><div>            AUTHORITY["EPSG","8901"]],</div><div>        UNIT["degree",0.0174532925199433,</div><div>            AUTHORITY["EPSG","9108"]],</div><div>        AUTHORITY["EPSG","4326"]],</div><div>    PROJECTION["Lambert_Conformal_Conic_2SP"],</div><div>    PARAMETER["standard_parallel_1",33],</div><div>    PARAMETER["standard_parallel_2",45],</div><div>    PARAMETER["latitude_of_origin",39.00000381472222],</div><div>    PARAMETER["central_meridian",-107],</div><div>    PARAMETER["false_easting",0],</div><div>    PARAMETER["false_northing",0],</div><div>    UNIT["METERS",1]]</div><div>Origin = (-633999.740890799090266,525999.667091785930097)</div><div>Pixel Size = (4000.000000000000000,-4000.000000000000000)</div><div>Corner Coordinates:</div><div>Upper Left  ( -633999.741,  525999.667) (114d52' 9.29"W, 43d30'32.73"N)</div><div>Lower Left  ( -633999.741, -526000.333) (113d52'56.93"W, 34d 1'21.79"N)</div><div>Upper Right (  634000.259,  525999.667) ( 99d 7'50.68"W, 43d30'32.73"N)</div><div>Lower Right (  634000.259, -526000.333) (100d 7' 3.05"W, 34d 1'21.79"N)</div><div>Center      (   0.2591092,  -0.3329082) (106d59'59.99"W, 39d 0' 0.00"N)</div><div>Band 1 Block=256x16 Type=Float32, ColorInterp=Undefined</div><div>  Min=274.508 Max=3824.624</div><div>  NoData Value=-3.4028234663852886e+38</div></div><div><br></div><div><br></div><div>prj.adf reads:</div><div><div>Projection    LAMBERT</div><div>Units         METERS</div><div>Zunits        NO</div><div>Xshift        0.0</div><div>Yshift        0.0</div><div>Parameters    6370000.0  6370000.0</div><div>  33  0  0.0 /* 1st standard parallel</div><div>  45  0  0.0 /* 2nd standard parallel</div><div>-107  0  0.0 /* central meridian</div><div>  39  0  0.013733 /* latitude of projection's origin</div><div>0.0 /* false easting (meters)</div><div>0.0 /* false northing (meters)</div></div><div><br></div><div><div>The file is here: </div><div><a href="https://www.odrive.com/s/4516b311-380e-4d94-8912-16f39c7d0259-56c21076">https://www.odrive.com/s/4516b311-380e-4d94-8912-16f39c7d0259-56c21076</a><br></div></div><div><br></div></div></div></div>
</div>