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
_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user
