[Bug 876241] New: GDAL fails to recognize projection in GeoTIFF DEM files (disables DEM in QLandkarteGT)
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c0 Summary: GDAL fails to recognize projection in GeoTIFF DEM files (disables DEM in QLandkarteGT) Classification: openSUSE Product: openSUSE.org Version: unspecified Platform: x86-64 OS/Version: openSUSE 13.1 Status: NEW Severity: Normal Priority: P5 - None Component: 3rd party software AssignedTo: opensuse-communityscreening@forge.provo.novell.com ReportedBy: pierre.baldensperger@numericable.fr QAContact: opensuse-communityscreening@forge.provo.novell.com Found By: --- Blocker: --- User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:29.0) Gecko/20100101 Firefox/29.0 Until today, I was using QLandkarteGT 1.7.6 and GDAL 1.10.1 from the Application:Geo OBS repository. Since January or so, I had noticed that the DEM functionality in QLandkarteGT had ceased working (only in the OpenSuSE version and not in the Windows version). Today, after discussing the matter with QLandkarteGT author Oliver Eichler through the QLandkarteGT mailing list (see http://article.gmane.org/gmane.comp.gis.qlandkartegt.user/1486), we apparently narrowed this down to a failure to recognize (and conversely generate) projection information in GeoTIFFs with the "Application:Geo" version of GDAL (version 1.10.1, arch x86_64 at least). Note that I updated this package several times since January so it is not a recent issue. Just un-installing that version and recompiling the latest GDAL from source magically fixes the problem. According to Oliver Eichler, it may be due to some binary mismatch between GDAL and one of its support libraries, so a simple "compilation refresh" might resolve the whole issue, but we are not sure. Reproducible: Always Steps to Reproduce: 1. Generate a GeoTIFF DEM file from NASA SRTM data following the instructions on the QLandkarteGT wiki : http://sourceforge.net/apps/mediawiki/qlandkartegt/index.php?title=Use_SRTM_... 2. Run the resulting file through the "gdalinfo" utility : the projection information is missing. 3. Uninstall OBS GDAL and recompile latest GDAL (1.11.0) from source (straightforward : configure->make->install). 4. Regenerate (again) the same GeoTIFF DEM file with the newly compiled GDAL. 5. Run the new GeoTIFF file through the new "gdalinfo" utility : the projection information has reappeared. Actual Results: For GeoTIFF files created according to the instructions of the QLandkarteGT wiki (see "steps" above), the "gdalinfo" should explicitly report parameters of the Mercator projection used, but those are completely absent (no "PROJECTION" information is reported at all). Expected Results: On such GeoTIFF DEM files, the "gdalinfo" output should report something along the lines of : PROJECTION["Mercator_1SP"], PARAMETER["central_meridian",0], PARAMETER["scale_factor",1], PARAMETER["false_easting",0], PARAMETER["false_northing",0], This problem prevents QLandkarteGT (also from Application:Geo repository) from working properly with Digital Elevation Models. Adding DEMs to maps silently fails and no altitude information is ever available. Even DEM files that load perfectly on other platforms (e.g. Windows version of QLGT) fail to load. After compiling the latest QLandkarteGT version 1.7.6 from source against the newly compiled GDAL 1.11.0, the QLGT DEM functionality is restored and works perfectly (but any GeoTIFFS generated with the "faulty" GDAL must first be regenerated). -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c1 Angelos Tzotsos <tzotsos@gmail.com> changed: What |Removed |Added ---------------------------------------------------------------------------- CC| |tzotsos@gmail.com --- Comment #1 from Angelos Tzotsos <tzotsos@gmail.com> 2014-05-03 17:11:15 UTC --- I am in the process of updating GDAL to 1.11 final. There is a dev package available at: https://build.opensuse.org/package/show/home:tzotsos:Application:Geo/libgdal Can you please test this version and check if will help? Best, Angelos -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c2 --- Comment #2 from Pierre Baldensperger <pierre.baldensperger@numericable.fr> 2014-05-03 18:54:15 UTC --- Thanks a lot for the quick reply. OK I just tested the new package from your repository, but the same problem is still there. For some reason it won't recognize the projection in my GeoTIFFs. I will try to attach one such GeoTIFF file, but unfortunately they are rather big so I'm not sure it's going to work. If it works, just run it through "gdalinfo" and its output should include Mercator projection parameters similar to the following : Driver: GTiff/GeoTIFF Files: N51E005.tif Size is 898, 1442 Coordinate System is: PROJCS["unnamed", GEOGCS["unnamed ellipse", DATUM["unknown", SPHEROID["unnamed",6378137,0]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]], PROJECTION["Mercator_1SP"], PARAMETER["central_meridian",0], PARAMETER["scale_factor",1], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] Origin = (556551.070845203124918,6800200.793425155803561) Pixel Size = (124.125864256673921,-124.125864256673921) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 556551.071, 6800200.793) ( 4d59'58.50"E, 52d 0' 1.50"N) Lower Left ( 556551.071, 6621211.297) ( 4d59'58.50"E, 50d59'58.32"N) Upper Right ( 668016.097, 6800200.793) ( 6d 0' 3.21"E, 52d 0' 1.50"N) Lower Right ( 668016.097, 6621211.297) ( 6d 0' 3.21"E, 50d59'58.32"N) Center ( 612283.584, 6710706.045) ( 5d30' 0.85"E, 51d30' 9.80"N) Band 1 Block=898x4 Type=Int16, ColorInterp=Gray "Faulty" GDAL binaries report a "LOCAL_CS" instead of the "PROJCS", with no associated "PROJECTION" fields nor Mercator parameters. -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c3 --- Comment #3 from Pierre Baldensperger <pierre.baldensperger@numericable.fr> 2014-05-03 19:00:41 UTC --- Created an attachment (id=588592) --> (http://bugzilla.novell.com/attachment.cgi?id=588592) test GeoTIFF file containing SRTM DEM data with Mercator projection Obtained from the NASA SRTM "hgt" data (file N51E005.hgt) following the instructions on the QLandkarteGT wiki : http://sourceforge.net/apps/mediawiki/qlandkartegt/index.php?title=Use_SRTM_... More precisely, the command line was : gdalwarp -t_srs \"+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs\" -r cubic N51E005.hgt N51E005.tif -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c4 --- Comment #4 from Pierre Baldensperger <pierre.baldensperger@numericable.fr> 2014-05-03 19:20:15 UTC --- Forgot to mention : the problem seems to originate in libgdal.so. Only replacing the OBS version of this file with a manual compilation fixes the problem (no need to replace the gdal utilities themselves). Here's the configure output of my manual GDAL compilation by the way : Installation directory: /usr/local C compiler: gcc -g -O2 -DHAVE_SSE_AT_COMPILE_TIME C++ compiler: g++ -g -O2 -DHAVE_SSE_AT_COMPILE_TIME LIBTOOL support: yes LIBZ support: external LIBLZMA support: no GRASS support: no CFITSIO support: external PCRaster support: internal LIBPNG support: external DDS support: no GTA support: no LIBTIFF support: external (BigTIFF=yes) LIBGEOTIFF support: internal LIBJPEG support: external 12 bit JPEG: no 12 bit JPEG-in-TIFF: no LIBGIF support: external OGDI support: no HDF4 support: yes HDF5 support: yes NetCDF support: yes Kakadu support: no JasPer support: yes (GeoJP2=no) OpenJPEG support: no ECW support: no MrSID support: no MrSID/MG4 Lidar support: no MSG support: no GRIB support: yes EPSILON support: no WebP support: no cURL support (wms/wcs/...):yes PostgreSQL support: no MySQL support: no Ingres support: no Xerces-C support: no NAS support: no Expat support: yes libxml2 support: yes Google libkml support: no ODBC support: yes PGeo support: yes FGDB support: no MDB support: no PCIDSK support: internal OCI support: no GEORASTER support: no SDE support: no Rasdaman support: no DODS support: no SQLite support: yes PCRE support: yes SpatiaLite support: no DWGdirect support no INFORMIX DataBlade support:no GEOS support: no Poppler support: no Podofo support: no OpenCL support: no Armadillo support: no FreeXL support: no SOSI support: no SWIG Bindings: no Statically link PROJ.4: no enable OGR building: yes enable pthread support: yes enable POSIX iconv support:yes hide internal symbols: no Looking at this, I see that it is using an internal LIBGEOTIFF. Now I wonder whether this could be the problem because I just checked the libgeotiff.so on my system and it is not the latest package from Application:Geo (it is version 1.3.0 from the Packman repo). I'll try to update that and test your build again. I'll post back when I'm done. -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c5 --- Comment #5 from Pierre Baldensperger <pierre.baldensperger@numericable.fr> 2014-05-03 19:35:32 UTC --- Bingo. The problem was with my outdated libgeotiff. Installed version 1.4.0 from the Application:Geo repository and everything works. Version 1.11.0 from your (Angelos's) repository also works but QLandkarteGT segfaults when I use it (probably should be recompiled against it) so I reverted back to version 1.10.1 from the Application:Geo repo, and QLandkarteGT (from the same repo) is fully functional. Bottom line : all my fault, very sorry for the inconvenience. -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c6 Pierre Baldensperger <pierre.baldensperger@numericable.fr> changed: What |Removed |Added ---------------------------------------------------------------------------- Status|NEW |RESOLVED Resolution| |INVALID --- Comment #6 from Pierre Baldensperger <pierre.baldensperger@numericable.fr> 2014-05-03 19:37:41 UTC --- Marking this as RESOLVED / INVALID. -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
https://bugzilla.novell.com/show_bug.cgi?id=876241 https://bugzilla.novell.com/show_bug.cgi?id=876241#c7 --- Comment #7 from Angelos Tzotsos <tzotsos@gmail.com> 2014-05-03 19:56:07 UTC --- Thanks. I will test 1.11 a bit more before I make the switch. -- Configure bugmail: https://bugzilla.novell.com/userprefs.cgi?tab=email ------- You are receiving this mail because: ------- You are on the CC list for the bug.
participants (1)
-
bugzilla_noreply＠novell.com