[GRASS-user] LiDAR data CRS issue

[Rich Shepard]

On Sat, 29 Apr 2023, Markus Neteler wrote:

> 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

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

Markus,

I found that ESPG code yesterday and created a location for it. The problem
is that the prj.adf file does not use that EPSG code:

ERROR: Projection of dataset does not appear to match current location.

        Location PROJ_INFO is:
        name: NAD83(2011) / Oregon North (ft)
        datum: nad83_2011
        ellps: grs80
        proj: lcc
        lat_0: 43.6666666666667
        lon_0: -120.5
        lat_1: 46
        lat_2: 44.3333333333333
        x_0: 2500000.0001424
        y_0: 0
        no_defs: defined

        Dataset PROJ_INFO is:
        name: unnamed
        ellps: grs80
        proj: lcc
        lat_0: 41.75
        lon_0: -120.5
        lat_1: 43
        lat_2: 45.5
        x_0: 400000
        y_0: 0
        no_defs: defined

        Difference in: datum

I've uploaded the Bare_Earth subdirectory to
<http://www.fileconvoy.com/dfl.php?id=g0f185de302ff7e1b100049146825ac005b4609a539>
where is will stay for 5 days.

There's also a Shapefiles/ directory with several maps. The
OLC_METRO_TAF_7_5_OGIC_2011.prj file contains:

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]

Should I use that one instead? I see the DEM and vector files using two
different CRS'.

Regards,

Rich