[Proj] [gdal-dev] gdalinfo error and subsequent segfault with proj-5.0.x

[Even Rouault]

Markus,

>From a quick test generating a GeoTIFF in EPSG:25832, I can't reproduce.

I suspect you may link at runtime against different proj versions. How did you 
build GDAL ? I'd strongly recommend you use --with-static-proj4 (whose name is 
a bit misleading since it can actually link against a libproj.so), and then 
check with "ldd libgdal.so" that you link against a single libproj.

Even

> I observed a segfault in gdalinfo (2.2.3, 2.2.4, 2.3.0dev-71e2ada881) when
> compiled against proj-5.0.0/proj-5.0.1
> 
> The dataset is a GeoTIFF with EPSG:25832, and I get
> 
> Corner Coordinates:
> ERROR 1: illegal axis orientation combination
> Upper Left  (  375000.000, 5631000.000)
> ERROR 1: illegal axis orientation combination
> Lower Left  (  375000.000, 5628000.000)
> ERROR 1: illegal axis orientation combination
> Upper Right (  379000.000, 5631000.000)
> ERROR 1: illegal axis orientation combination
> Lower Right (  379000.000, 5628000.000)
> ERROR 1: illegal axis orientation combination
> Center      (  377000.000, 5629500.000)
> Segmentation fault (core dumped)
> 
> The error message is printed by GDAL and comes from PROJ. I think the
> segfault is a symptom of the error, therefore I would like to find out the
> reason of the error.
> 
> Using gdal-2.2+ with proj-4.9.3, all is fine.
> 
> Testing proj-5.0.1 with
> 
> cs2cs +proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m
> +no_defs +to +proj=latlong +ellps=GRS80
> 
> proj -I +proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m
> +no_defs
> 
> cct +proj=pipeline +step +inv +proj=utm +zone=32 +ellps=GRS80
> +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
> 
> all is fine, too.
> 
> Is this a bug in (the old API still available with) proj-5.0.x or a bug in
> gdal-2.2+?
> Any hints where and how to dig deeper?
> 
> Markus M


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com