[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