[GRASS-user] v.in.ogr: .gdb has multiple layers, different projections
Rich Shepard
rshepard at appl-ecosys.com
Mon Aug 17 15:43:00 PDT 2020
ogrinfo -al -so on the .gdb produces this output:
Layer name: WSDACrop_2016
Geometry: Multi Polygon
Feature Count: 198132
Extent: (194040.798600, 26290.144100) - (776078.287500, 411036.825500)
Layer SRS WKT:
PROJCRS["NAD83(HARN) / Washington South",
BASEGEOGCRS["NAD83(HARN)",
DATUM["NAD83 (High Accuracy Reference Network)",
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4152]],
CONVERSION["SPCS83 Washington South zone (meters)",
METHOD["Lambert Conic Conformal (2SP)",
ID["EPSG",9802]],
PARAMETER["Latitude of false origin",45.3333333333333,
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",47.3333333333333,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8823]],
PARAMETER["Latitude of 2nd standard parallel",45.8333333333333,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8824]],
PARAMETER["Easting at false origin",500000,
LENGTHUNIT["metre",1],
ID["EPSG",8826]],
PARAMETER["Northing at false origin",0,
LENGTHUNIT["metre",1],
ID["EPSG",8827]]],
CS[Cartesian,2],
AXIS["easting (X)",east,
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["northing (Y)",north,
ORDER[2],
LENGTHUNIT["metre",1]],
USAGE[
SCOPE["unknown"],
AREA["USA - Washington - SPCS83 - S"],
BBOX[45.54,-124.4,47.61,-116.91]],
ID["EPSG",2856]]
Data axis to CRS axis mapping: 1,2
FID Column = OBJECTID
Geometry Column = SHAPE
CropGroup: String (254.0)
Acres: Integer (0.0)
Irrigation: String (254.0)
InitialSurveyDate: DateTime (0.0)
LastSurveyDate: DateTime (0.0)
DataSource: String (254.0)
Notes: String (254.0)
TRS: String (15.0)
County: String (254.0)
RotCropType: String (254.0)
ExactAcres: Real (0.0)
Organic: String (3.0)
SHAPE_Length: Real (0.0)
SHAPE_Area: Real (0.0)
Layer name: CropData
Geometry: None
Feature Count: 198132
Layer SRS WKT:
(unknown)
FID Column = OBJECTID
CropGroup: String (254.0)
CropType: String (254.0)
Acres: Integer (0.0)
Irrigation: String (254.0)
InitialSurveyDate: DateTime (0.0)
LastSurveyDate: DateTime (0.0)
DataSource: String (254.0)
Notes: String (254.0)
TRS: String (15.0)
County: String (254.0)
RotCropType: String (254.0)
ExactAcres: Real (0.0)
Organic: String (3.0)
SHAPE_Length: Real (0.0)
SHAPE_Area: Real (0.0)
I see two layers (WSDACrop_2016 and CropData) but only one projection: EPSG
2856.
Trying to import it results in this error:
ERROR: Detected different projections of input layers. Input layers must be
imported separately.
What am I not seeing? Or, how should I import it other than v.in.ogr in=<>
out<>?
TIA,
Rich
More information about the grass-user
mailing list