date:20230429

Re: [GRASS-user] LiDAR data CRS issue

2023-04-29 Thread 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.9,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.6,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.7 +lon_0=-120.5 +lat_1=46
+lat_2=44.3 +x_0=250.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.7
   lon_0: -120.5
   lat_1: 46
   lat_2: 44.3
   x_0: 250.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: 40
   y_0: 0
   no_defs: defined

   Difference in: datum

I've uploaded the Bare_Earth subdirectory to

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

___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] LiDAR data CRS issue

2023-04-29 Thread Markus Neteler

On Sat, Apr 29, 2023 at 5:40 PM Rich Shepard  wrote:
>
> I've defined a location using NAD83(2011) / Oregon North (ft). But, the
> prj.adf file for the data is:
> ProjectionLAMBERT
> Spheroid  GRS80
> Units 3.280839895013123
> ZunitsNO
> Xshift0.0
> Yshift0.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
> 40.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=false=false=EPSG=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.9,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.6,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.7 +lon_0=-120.5 +lat_1=46
+lat_2=44.3 +x_0=250.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.7,
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.3,
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
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] LiDAR data CRS issue

2023-04-29 Thread Rich Shepard


I've defined a location using NAD83(2011) / Oregon North (ft). But, the
prj.adf file for the data is:
ProjectionLAMBERT
Spheroid  GRS80
Units 3.280839895013123
ZunitsNO
Xshift0.0
Yshift0.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
40.0 /* false easting (meters)
0.0 /* false northing (meters)

This is EPSG: 7019, but it's not a CRS, is it?

How do I define a location that will accept these data?

Rich
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] LiDAR data CRS issue

Re: [GRASS-user] LiDAR data CRS issue

[GRASS-user] LiDAR data CRS issue

3 matches

Site Navigation

Mail list logo

Footer information