On jeudi 26 octobre 2017 13:52:21 CEST Rich Shepard wrote: > I am trying to import/open data files from www.streamnet.org but their > GDB format is quite different from what I've seen before now and > grass-7.3.svn will not open them. As an aside, I was thrown for a while > looking for OpenFileGDB in the v.in.ogr supported format list until I > noticed it's been renamed to ESRIFileGDB. > > The results of running 'ogrinfo -al so' on the directory, > BullTrout_201201_v931.gdb/, produces this result: > > Had to open data source read-only. > INFO: Open of `BullTrout_201201_v931.gdb/' > using driver `OpenFileGDB' successful. > > Layer name: FishD_BullTrout_January2012 > Geometry: Measured Multi Line String > Feature Count: 4620 > Extent: (1158336.501590, 288917.146174) - (4143404.913877, 3148258.748272) > Layer SRS WKT: > PROJCS["NAD_1983_Lambert_Conformal_Conic", > GEOGCS["GCS_North_American_1983", > DATUM["North_American_Datum_1983", > SPHEROID["GRS_1980",6378137.0,298.257222101]], > PRIMEM["Greenwich",0.0], > UNIT["Degree",0.0174532925199433]], > PROJECTION["Lambert_Conformal_Conic_2SP"], > PARAMETER["False_Easting",2999999.999988], > PARAMETER["False_Northing",0.0], > PARAMETER["Central_Meridian",-117.0], > PARAMETER["Standard_Parallel_1",42.33333333333334], > PARAMETER["Standard_Parallel_2",48.66666666666666], > PARAMETER["Latitude_Of_Origin",41.0], > UNIT["Foot_US",0.3048006096012192]] > FID Column = OBJECTID > Geometry Column NOT NULL = Shape > RecordID: Integer (0.0) > LocationID: String (13.0) > BegFt: Integer (0.0) > EndFt: Integer (0.0) > StreamName: String (100.0) > StreamAndTribNames: String (160.0) > SpecieID: Integer(Int16) (0.0) > Species: String (50.0) > SciName: String (50.0) > RunID: Integer(Int16) (0.0) > Run: String (50.0) > SubRunID: Integer(Int16) (0.0) > SubRun: String (50.0) > LifeHistoryID: Integer(Int16) (0.0) > LifeHistoryType: String (80.0) > UseTypeID: Integer(Int16) (0.0) > UseType: String (50.0) > RefID: Integer (0.0) > Citation: String (3600.0) > URL: String (1000.0) > CompAgy: String (12.0) > BasisID: Integer(Int16) (0.0) > Basis: String (80.0) > UpdDate: DateTime (0.0) > EndExtentID: Integer(Int16) (0.0) > Shape_Length: Real (0.0) > > There's no ESPG code there. So I wrote a proj4 file > > +datum=NAD83 +ellps=GRS80 +lat_0=41.0 +lat_1=42.33333333333334 > +lat_2=48.66666666666666 +proj=lcc +units=foot_us +x_0=2999999.999988 > +y_0=0.0
You forgot +lon_0, and a bit counter-intuitively , proj.4 requires x_0 to be expressed in meters (*) So try instead +proj=lcc +lat_1=42.33333333333334 +lat_2=48.66666666666666 +lat_0=41 +lon_0=-117 +x_0=914401.8287999999 +y_0=0 +datum=NAD83 +units=us-ft +no_defs (hint: was given by "gdalsrsinfo the.gdb") Even (*) See http://proj4.org/parameters.html#false-easting-northing -- Spatialys - Geospatial professional services http://www.spatialys.com
_______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
