[GRASS-user] Finding EPSG code based on prj.adf; importing with v.in.gdal
Rich Shepard
rshepard at appl-ecosys.com
Wed Jul 4 15:44:33 PDT 2018
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."
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
More information about the grass-user
mailing list