[GRASS-user] Unusual [Open|ESRI]FileGDB: how to import

[Even Rouault]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20171026/543c3c3a/attachment-0001.html>