[GRASS-user] rshepard at appl-ecosys.com sent you a file

[Markus Neteler]

On Sun, Jun 10, 2018 at 6:41 PM, Rich Shepard <rshepard at appl-ecosys.com> wrote:
> On Sun, 10 Jun 2018, Micha Silver wrote:
>
>> Now it works

.. also here, git the file.

>> I created a new LOCATION based on "a georeferenced file" and pointed to
>> the hdr.adf file.
>> The location projection parameters:
>
>
> Micha,
>
>   Now to find why it's not working the same way here with 7.5svn release
> 72793.

When using gdalinfo, it does not recognize the SRS (on my system at least).
[mneteler at oboe Bare_Earth ]$ gdalinfo bh45122c3
Driver: AIG/Arc/Info Binary Grid
Files: bh45122c3
       bh45122c3.ovr
       bh45122c3.aux.xml
       bh45122c3/sta.adf
       bh45122c3/w001001x.adf
       bh45122c3/prj.adf
       bh45122c3/metadata.xml
       bh45122c3/hdr.adf
       bh45122c3/w001001.adf
       bh45122c3/dblbnd.adf
Size is 11052, 15426
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["Unknown datum based upon the GRS 1980 ellipsoid",
        DATUM["Not_specified_based_on_GRS_1980_ellipsoid",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6019"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4019"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",43],
    PARAMETER["standard_parallel_2",45.5],
    PARAMETER["latitude_of_origin",41.75],
    PARAMETER["central_meridian",-120.5],
    PARAMETER["false_easting",1312335.958005249],
    PARAMETER["false_northing",0],
    UNIT["user-defined",0.3048]]
Origin = (829504.500000000000000,1326862.500000000000000)
Pixel Size = (3.000000000000000,-3.000000000000000)
Corner Coordinates:
Upper Left  (  829504.500, 1326862.500) (122d22'44.53"W, 45d22'29.77"N)
Lower Left  (  829504.500, 1280584.500) (122d22'29.72"W, 45d14'52.95"N)
Upper Right (  862660.500, 1326862.500) (122d15' 0.15"W, 45d22'37.00"N)
Lower Right (  862660.500, 1280584.500) (122d14'46.37"W, 45d15' 0.16"N)
Center      (  846082.500, 1303723.500) (122d18'45.19"W, 45d18'45.04"N)
Band 1 Block=256x16 Type=Float32, ColorInterp=Undefined
  Min=206.380 Max=1157.470
  Minimum=206.380, Maximum=1157.470, Mean=526.722, StdDev=308.445
  NoData Value=-3.4028234663852886e+38
  Overviews: 5526x7713, 2763x3857, 1382x1929, 691x965, 346x483, 173x242
  Metadata:
    STATISTICS_MAXIMUM=1157.4699707031
    STATISTICS_MEAN=526.72212274996
    STATISTICS_MINIMUM=206.38000488281
    STATISTICS_STDDEV=308.44451906832

However, I discovered the SRS in the included XML file:

grep SRS bh45122c3.aux.xml
<SRS>PROJCS["NAD_1983_2011_Oregon_Statewide_Lambert_Feet_Intl",GEOGCS["GCS_NAD_1983_2011",DATUM["D_NAD_1983_2011",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["false_easting",1312335.958005249],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",-120.5],PARAMETER["standard_parallel_1",43.0],PARAMETER["standard_parallel_2",45.5],PARAMETER["latitude_of_origin",41.75],UNIT["Foot",0.3048]]</SRS>

Apparently the SRS lookup mechanism fails in GDAL? GRASS GIS mostly
relies on it.

I have GDAL 2.2.4, released 2018/03/19, does anyone have a more recent
version and could try?

Markus