<div dir="ltr"><div><div>Even,<br><br>On Sun, Apr 15, 2018 at 8:17 PM, Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br>><br>> Markus,<br>><br>> From a quick test generating a GeoTIFF in EPSG:25832, I can't reproduce.<br>><br>> I suspect you may link at runtime against different proj versions.<br><br></div>You are right. I compiled GDAL with<br><br>libproj.so -> libproj.so.13.0.1<br></div><div>there is no other libproj.so on my system<br></div><div><br></div>but gdal links against libproj.so.12<br><div><div><br>> How did you<br>> build GDAL ? I'd strongly recommend you use --with-static-proj4 (whose name is<br>> a bit misleading since it can actually link against a libproj.so), and then<br>> check with "ldd libgdal.so" that you link against a single libproj.<br><br></div><div>OK, I recompiled gdal-2.2.4 --with-static-proj4 and got<br><br>ldd libgdal.so | grep proj<br>    libproj.so.13 => /usr/local/lib64/libproj.so.13 (0x00007fe2ff181000)<br>    libproj.so.12 => /lib64/libproj.so.12 (0x00007fe2f31c7000)<br>???<br></div><div><br></div><div>gdalinfo now runs fine and produces expected results.<br><br></div><div>I'm still concerned about the output of ldd libgdal.so<br><br></div><div>Markus M<br></div><div>><br>> Even<br>><br>> > I observed a segfault in gdalinfo (2.2.3, 2.2.4, 2.3.0dev-71e2ada881) when<br>> > compiled against proj-5.0.0/proj-5.0.1<br>> ><br>> > The dataset is a GeoTIFF with EPSG:25832, and I get<br>> ><br>> > Corner Coordinates:<br>> > ERROR 1: illegal axis orientation combination<br>> > Upper Left  (  375000.000, 5631000.000)<br>> > ERROR 1: illegal axis orientation combination<br>> > Lower Left  (  375000.000, 5628000.000)<br>> > ERROR 1: illegal axis orientation combination<br>> > Upper Right (  379000.000, 5631000.000)<br>> > ERROR 1: illegal axis orientation combination<br>> > Lower Right (  379000.000, 5628000.000)<br>> > ERROR 1: illegal axis orientation combination<br>> > Center      (  377000.000, 5629500.000)<br>> > Segmentation fault (core dumped)<br>> ><br>> > The error message is printed by GDAL and comes from PROJ. I think the<br>> > segfault is a symptom of the error, therefore I would like to find out the<br>> > reason of the error.<br>> ><br>> > Using gdal-2.2+ with proj-4.9.3, all is fine.<br>> ><br>> > Testing proj-5.0.1 with<br>> ><br>> > cs2cs +proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m<br>> > +no_defs +to +proj=latlong +ellps=GRS80<br>> ><br>> > proj -I +proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m<br>> > +no_defs<br>> ><br>> > cct +proj=pipeline +step +inv +proj=utm +zone=32 +ellps=GRS80<br>> > +towgs84=0,0,0,0,0,0,0 +units=m +no_defs<br>> ><br>> > all is fine, too.<br>> ><br>> > Is this a bug in (the old API still available with) proj-5.0.x or a bug in<br>> > gdal-2.2+?<br>> > Any hints where and how to dig deeper?<br>> ><br>> > Markus M<br>><br>><br>> --<br>> Spatialys - Geospatial professional services<br>> <a href="http://www.spatialys.com">http://www.spatialys.com</a><br><br></div></div></div>