[GRASS-user] LiDAR data CRS issue

[Markus Neteler]

On Sat, Apr 29, 2023 at 5:40 PM Rich Shepard <rshepard at appl-ecosys.com> wrote:
>
> I've defined a location using NAD83(2011) / Oregon North (ft). But, the
> prj.adf file for the data is:
> Projection    LAMBERT
> Spheroid      GRS80
> Units         3.280839895013123
> Zunits        NO
> Xshift        0.0
> Yshift        0.0
> Parameters
>    43  0  0.0 /* 1st standard parallel
>    45 30  0.0 /* 2nd standard parallel
> -120 30  0.0 /* central meridian
>    41 45  0.0 /* latitude of projection's origin
> 400000.0 /* false easting (meters)
> 0.0 /* false northing (meters)
>
> This is EPSG: 7019, but it's not a CRS, is it?

Here you may find valid codes for a specific area:

https://crs-explorer.proj.org/?latlng=44.552901,-123.266602&ignoreWorld=false&allowDeprecated=false&authorities=EPSG&activeTypes=PROJECTED_CRS

A click then brings you to EPSG.org, Find in the upper right a tiny
"WKT" button which shows eg.:

https://epsg.org/crs/wkt/id/6559
-> WKT:

PROJCRS["NAD83(2011) / Oregon North
(ft)",BASEGEOGCRS["NAD83(2011)",DATUM["NAD83 (National Spatial
Reference System 2011)",ELLIPSOID["GRS
1980",6378137,298.257222101,LENGTHUNIT["metre",1,ID["EPSG",9001]],ID["EPSG",7019]],ID["EPSG",1116]],ID["EPSG",6318]],CONVERSION["SPCS83
Oregon North zone (International feet)",METHOD["Lambert Conic
Conformal (2SP)",ID["EPSG",9802]],PARAMETER["Latitude of false
origin",43.6666666666669,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9102]],ID["EPSG",8821]],PARAMETER["Longitude
of false origin",-120.5,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9102]],ID["EPSG",8822]],PARAMETER["Latitude
of 1st standard
parallel",46,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9102]],ID["EPSG",8823]],PARAMETER["Latitude
of 2nd standard
parallel",44.3333333333336,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9102]],ID["EPSG",8824]],PARAMETER["Easting
at false origin",8202099.738,LENGTHUNIT["foot",0.3048,ID["EPSG",9002]],ID["EPSG",8826]],PARAMETER["Northing
at false origin",0,LENGTHUNIT["foot",0.3048,ID["EPSG",9002]],ID["EPSG",8827]],ID["EPSG",15351]],CS[Cartesian,2,ID["EPSG",4495]],AXIS["Easting
(X)",east],AXIS["Northing
(Y)",north],LENGTHUNIT["foot",0.3048,ID["EPSG",9002]],ID["EPSG",6559]]

More human readable output you get with `projinfo`:

projinfo EPSG:6559
PROJ.4 string:
+proj=lcc +lat_0=43.6666666666667 +lon_0=-120.5 +lat_1=46
+lat_2=44.3333333333333 +x_0=2500000.0001424 +y_0=0 +ellps=GRS80
+towgs84=0,0,0,0,0,0,0 +units=ft +no_defs +type=crs

WKT2:2019 string:
PROJCRS["NAD83(2011) / Oregon North (ft)",
    BASEGEOGCRS["NAD83(2011)",
        DATUM["NAD83 (National Spatial Reference System 2011)",
            ELLIPSOID["GRS 1980",6378137,298.257222101,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",6318]],
    CONVERSION["SPCS83 Oregon North zone (International feet)",
        METHOD["Lambert Conic Conformal (2SP)",
            ID["EPSG",9802]],
        PARAMETER["Latitude of false origin",43.6666666666667,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8821]],
        PARAMETER["Longitude of false origin",-120.5,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8822]],
        PARAMETER["Latitude of 1st standard parallel",46,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8823]],
        PARAMETER["Latitude of 2nd standard parallel",44.3333333333333,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8824]],
        PARAMETER["Easting at false origin",8202099.738,
            LENGTHUNIT["foot",0.3048],
            ID["EPSG",8826]],
        PARAMETER["Northing at false origin",0,
            LENGTHUNIT["foot",0.3048],
            ID["EPSG",8827]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["foot",0.3048]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["foot",0.3048]],
    USAGE[
        SCOPE["Engineering survey, topographic mapping."],
        AREA["United States (USA) - Oregon - counties of Baker;
Benton; Clackamas; Clatsop; Columbia; Gilliam; Grant; Hood River;
Jefferson; Lincoln; Linn; Marion; Morrow; Multnomah; Polk; Sherman;
Tillamook; Umatilla; Union; Wallowa; Wasco; Washington; Wheeler;
Yamhill."],
        BBOX[43.95,-124.17,46.26,-116.47]],
    ID["EPSG",6559]]

Note that I just showed an example, it is not necessarily the code you
are hunting for.
Just to illustrate a method.

Cheers,
Markus

-- 
Markus Neteler, PhD
https://www.mundialis.de - free data with free software
https://grass.osgeo.org
https://courses.neteler.org/blog