Re: [gdal-dev] The "Bertin 1953" projection

2018-11-13 Thread Nikos Alexandris
* Even Rouault [2018-11-07 12:10:19 +0100]: ``` echo 12 55 0 0 | proj "+proj=bertin1953 +wktext" Rel. 6.0.0, March 1st, 2019 : projection initialization failure cause: unknown projection id program abnormally terminated ``` 'proj' syntax is without quoting. And you only need +wktext for GDAL,

Re: [gdal-dev] ogr2ogr clip operation very slow compared to QGIS

2018-11-13 Thread Nyall Dawson
On Wed, 14 Nov 2018 at 01:03, Even Rouault wrote: > Looking at QGIS code, it uses GEOS prepared geometry mechanism with the clip > geometry, so as to speed up intersection operations on it. Whereas ogr2ogr - > clipsrc uses the 'naive' intersection function. I guess this is the reason for > the sp

Re: [gdal-dev] Non-standard field types in OGR

2018-11-13 Thread Even Rouault
Hi Matthias, > > I'm looking for a way to load a geopackage layer with a type affinity > set to a custom type and to detect this type. This would likely be illegal for core GeoPackage spec since data types are standardized per http://www.geopackage.org/spec/#r5 . But could be made an extension

[gdal-dev] Non-standard field types in OGR

2018-11-13 Thread Matthias Kuhn
Hi, I'm looking for a way to load a geopackage layer with a type affinity set to a custom type and to detect this type. I.e. what I have is a table with a field called myJson, the type of which is JSON. Looking at the debug output of QGIS, I currently get this message which seems to be generated

Re: [gdal-dev] Confusing error behavior of OGR_L_GetSpatialRef

2018-11-13 Thread Even Rouault
On mardi 13 novembre 2018 14:13:27 CET Sean Gillies wrote: > Hi, > > While debugging a problem with my Fiona project I realized that if > GDAL_DATA is not known, OGR_L_GetSpatialRef() (master branch and also > 2.3.2) can return a not-NULL spatial reference while also putting a > CE_Failure error o

[gdal-dev] Confusing error behavior of OGR_L_GetSpatialRef

2018-11-13 Thread Sean Gillies
Hi, While debugging a problem with my Fiona project I realized that if GDAL_DATA is not known, OGR_L_GetSpatialRef() (master branch and also 2.3.2) can return a not-NULL spatial reference while also putting a CE_Failure error on the error stack. Is this intended or is the error code perhaps one le

Re: [gdal-dev] ogr2ogr clip operation very slow compared to QGIS

2018-11-13 Thread Even Rouault
Peter, > > I'm trying to clip some GML data via a fairly complex ShapeFile boundary > using standard ogr2ogr and it is taking roughly 5 minutes before completing > successfully. Does your input dataset (input_file.gz) has many features ? (let's say ~1000) > This same operation takes less than a

[gdal-dev] ogr2ogr clip operation very slow compared to QGIS

2018-11-13 Thread Peter Marlow
Hi, I'm trying to clip some GML data via a fairly complex ShapeFile boundary using standard ogr2ogr and it is taking roughly 5 minutes before completing successfully. This same operation takes less than a second using the Clip operation in QGIS with exactly the same data. I was lead to believe