[GRASS-user] Finding EPSG code based on prj.adf; importing with v.in.gdal

[Markus Metz]

On Thu, Jul 5, 2018 at 12:44 AM, Rich Shepard <rshepard at appl-ecosys.com>
wrote:
>
>   The ESRI GRID formatted LiDAR raw data file can be imported using
> r.in.gdal using the appropriately defined location. Running 'gdalinfo
> <filename> | less' on the prj.adf file reports:
>
> 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",44.33333333333334],
>     PARAMETER["standard_parallel_2",46],
>     PARAMETER["latitude_of_origin",43.66666666666666],
>     PARAMETER["central_meridian",-120.5],
>     PARAMETER["false_easting",8202099.737532808],
>     PARAMETER["false_northing",0],
>     UNIT["user-defined",0.3048]]
> Origin = (7719274.500000000000000,628402.500000000000000)
> Pixel Size = (3.000000000000000,-3.000000000000000)
>
>   Looking up EPSG:4019 on the web shows me that it's an "Unknown datum
based upon the GRS 1980
> ellipsoid."

EPSG:4019 is not the EPSG code of the SRS, but of one of its components
(there are 5 different EPSG codes within this SRS). The SRS itself does not
have an EPSG code, therefore you can create a new location with g.proj
georef= or r.in.gdal location=

Markus M

>
>   I created a new location, named EPSG4019, selected that code for the
> projection, and ran 'r.in.gdal -o
>
in=$HOME/path_to/data/topography/LDQ-45122C3/2009_OLC_Hood-to-Coast/Bare_Earth/\
> be45122c3/hdr.adf out=dem2009 --overwrite' (Entered on one line; split
here
> for MUA accommodation.) grass75 responds,
>
> WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised
>          by GRASS and no parameters found
> Over-riding projection check
> ERROR: Illegal latitude for North: 628402
>
>   Please show me what I've done incorrectly.
>
> TIA,
>
> Rich
>
>
>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180705/5bbfc904/attachment.html>