Menu
▾
▴
Re: [OSSIM] Weird behaviour when reading tiff files (geographic referencing)
|
From: Roberta R. <rob...@un...> - 2018-09-23 09:11:46
|
Hi Garrett, now everything works fine! Thank you again for all the support!! Cheers, Roberta -------------------------------------------------------------------------------- Roberta Ravanelli, *PhD * Geodesy and Geomatics Division University of Rome "La Sapienza" Via Eudossiana, 18 - 00184 Rome Italy E-mail rob...@un... Il giorno dom 23 set 2018 alle ore 02:46 Garrett Potts < gar...@gm...> ha scritto: > Hello: > > Shoot! I had it still checked out locally and never pushed! Sorry! > Try it now > > > Take care > > Garrett > > > On Sep 22, 2018, at 5:53 PM, Roberta Ravanelli < > rob...@un...> wrote: > > Thank you again, but now Ossim does not even build: > > > */home/roberta/library_sources/ossim/ossim/src/support_data/ossimGeoTiff.cpp: > In member function ‘bool ossimGeoTiff::addImageGeometry(ossimKeywordlist&, > const char*) const’:* > */home/roberta/library_sources/ossim/ossim/src/support_data/ossimGeoTiff.cpp:1711:26: > error: assignment of member ‘ossimGeoTiff::theAngularUnits’ in read-only > object* > * theAngularUnits = ANGULAR_DEGREE;* > [ 59%] Building CXX object > ossim/src/CMakeFiles/ossim.dir/support_data/ossimRpfFrameFileIndexSubsection.cpp.o > [ 59%] Building CXX object > ossim/src/CMakeFiles/ossim.dir/support_data/ossimFfL5.cpp.o > ossim/src/CMakeFiles/ossim.dir/build.make:20438: set di istruzioni per > l'obiettivo > "ossim/src/CMakeFiles/ossim.dir/support_data/ossimGeoTiff.cpp.o" non > riuscito > make[2]: *** > [ossim/src/CMakeFiles/ossim.dir/support_data/ossimGeoTiff.cpp.o] Errore 1 > make[2]: *** Attesa per i processi non terminati.... > CMakeFiles/Makefile2:144: set di istruzioni per l'obiettivo > "ossim/src/CMakeFiles/ossim.dir/all" non riuscito > make[1]: *** [ossim/src/CMakeFiles/ossim.dir/all] Errore 2 > Makefile:127: set di istruzioni per l'obiettivo "all" non riuscito > make: *** [all] Errore 2 > > > Bests, > > Roberta > > > -------------------------------------------------------------------------------- > Roberta Ravanelli, *PhD * > Geodesy and Geomatics Division > University of Rome "La Sapienza" > Via Eudossiana, 18 - 00184 Rome Italy > E-mail rob...@un... > > > > Il giorno sab 22 set 2018 alle ore 22:53 Garrett Potts < > gar...@gm...> ha scritto: > >> Hello: >> >> I added a check to see if the angular units is not set then if the model >> type is geographic then default it to DEGREES. I just checked in the code. >> >> >> Let me know if that works out >> >> >> Take care >> >> Garrett >> >> >> On Sep 22, 2018, at 1:33 PM, Roberta Ravanelli < >> rob...@un...> wrote: >> >> Hi Garrett, >> >> the problem occurs with all the SRTM tiles that I have in my pc. >> >> I try to explain better my situation: we are trying to resume the >> development of the OSSIM DATE plugin that Martina has implented during her >> GSoCs. >> The bug happens at very beginning of the workflow, when we try to >> retrieve the heights of the grid points (line 354 of >> https://github.com/martidi/opencv_dsm/blob/imageStack/apps/ossim-opencv.cpp >> file). >> >> In particular, the program remains stuck within the while loop of the >> *ossimImageElevationDatabase::createCell* method (line 127 of >> https://github.com/ossimlabs/ossim/blob/dev/src/elevation/ossimImageElevationDatabase.cpp: >> for an unknown reason Ossim cannot read the metadata of the SRTM tiles. For >> instance, you can see my shell output below: >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_37_04.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_39_05.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_40_04.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_59_08.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_22_03.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_42_05.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_39_03.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_38_04.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_38_03.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_40_05.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_39_04.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> >> file scansionato >> /home/roberta/library_sources/ossim/data/elevation/srtm/v4_1/srtm_37_04.tif >> >> WARNING ossimGeoTiff::addImageGeometry: >> >> Not coded yet for unit type: 0 >> >> Upper Left and Low Right ( nan, nan, nan, WGE ), ( nan, nan, nan, WGE ) >> >> and so on forever... >> >> >> >> Here you can see all the calls that have brought the program to the >> warning messages: >> >> Thread 1 "ossim-DATE-plug" hit Breakpoint 1, 0x00007ffff6988720 in >> ossimGeoTiff::addImageGeometry(ossimKeywordlist&, char const*) const@plt >> () >> >> from /home/roberta/library_sources/ossim/build/lib/libossim.so.1 >> >> (gdb) bt >> >> #0 0x00007ffff6988720 in >> ossimGeoTiff::addImageGeometry(ossimKeywordlist&, char const*) const@plt >> () >> >> from /home/roberta/library_sources/ossim/build/lib/libossim.so.1 >> >> #1 0x00007ffff704bba7 in ossimTiffProjectionFactory::createProjection ( >> >> this=0x6d3860, handler=0x6dcac0) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/projection/ossimTiffProjectionFactory.cpp:99 >> >> #2 0x00007ffff70f1308 in >> ossimProjectionFactoryRegistry::createProjection ( >> >> this=0x7ffff7b6cbc0 >> <ossimProjectionFactoryRegistry::instance()::inst>, >> >> handler=0x6dcac0) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/projection/ossimProjectionFactoryRegistry.cpp:81 >> >> #3 0x00007ffff6c6fe50 in ossimImageGeometryFactory::createProjection ( >> >> this=0x6d3780, handler=0x6dcac0) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/imaging/ossimImageGeometryFactory.cpp:134 >> >> #4 0x00007ffff6c6fb37 in ossimImageGeometryFactory::extendGeometry ( >> >> this=0x6d3780, handler=0x6dcac0) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/imaging/ossimImageGeometryFactory.cpp:87 >> >> #5 0x00007ffff6d4a65b in ossimImageGeometryRegistry::extendGeometry ( >> >> this=0x6d3720, handler=0x6dcac0) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/imaging/ossimImageGeometryR---Type >> <return> to continue, or q <return> to quit--- >> >> egistry.cpp:34 >> >> #6 0x00007ffff6c0b332 in ossimImageHandler::getImageGeometry >> (this=0x6dcac0) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/imaging/ossimImageHandler.cpp:712 >> >> #7 0x00007ffff6bb0e9c in ossimImageElevationHandler::open >> (this=0x6dc490, >> >> file=...) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/elevation/ossimImageElevationHandler.cpp:76 >> >> #8 0x00007ffff6ba3de8 in ossimImageElevationDatabase::createCell ( >> >> this=0x6d44e0, gpt=...) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/elevation/ossimImageElevationDatabase.cpp:136 >> >> #9 0x00007ffff6ba4395 in >> ossimImageElevationDatabase::getOrCreateCellHandler ( >> >> this=0x6d44e0, gpt=...) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/elevation/ossimImageElevationDatabase.cpp:258 >> >> #10 0x00007ffff6ba3b62 in ossimImageElevationDatabase::getHeightAboveMSL ( >> >> this=0x6d44e0, gpt=...) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/elevation/ossimImageElevationDatabase.cpp:95 >> >> #11 0x00007ffff6bb52c2 in ossimElevManager::getHeightAboveMSL ( >> >> this=0x7ffff79e3240 <ossimElevManager::instance()::inst>, gpt=...) >> >> at >> /home/roberta/library_sources/ossim/ossim/src/elevation/ossimElevManager.---Type >> <return> to continue, or q <return> to quit--- >> >> cpp:146 >> >> #12 0x0000000000419b94 in main (argc=10, argv=0x7fffffffdc18) >> >> at >> /home/roberta/library_sources/ossim/opencv_dsm/apps/ossim-opencv.cpp:354 >> >> >> >> I do not understand what is happening! Maybe now the problems depends on >> my pc? But this seems strange to me, since before the fix the program was >> working (reading a wrong georeferencing, but at least it arrived to the >> end). >> >> Thank you for all the help! >> >> Roberta >> >> -------------------------------------------------------------------------------- >> Roberta Ravanelli, *PhD * >> Geodesy and Geomatics Division >> University of Rome "La Sapienza" >> Via Eudossiana, 18 - 00184 Rome Italy >> E-mail rob...@un... >> >> >> >> Il giorno sab 22 set 2018 alle ore 15:12 Garrett Potts < >> gar...@gm...> ha scritto: >> >>> Hello Roberta: >>> >>> >>> Seems that the error message means it couldn’t determine the angular >>> units. Could I get a sample or a link from where you downloaded the image? >>> >>> >>> Take care >>> >>> >>> Garrett >>> >>> >>> On Sep 21, 2018, at 3:57 PM, Roberta Ravanelli < >>> rob...@un...> wrote: >>> >>> Thank you very much for the very fast fix, but now I get the following >>> warning when I call the method >>> ossimElevManager::instance()->getHeightAboveMSL(world_point); >>> >>> WARNING ossimGeoTiff::addImageGeometry: >>> Not coded yet for unit type: 0 >>> >>> Any hint? >>> >>> Thank you again! >>> >>> Roberta >>> >>> >>> -------------------------------------------------------------------------------- >>> Roberta Ravanelli, *PhD * >>> Geodesy and Geomatics Division >>> University of Rome "La Sapienza" >>> Via Eudossiana, 18 - 00184 Rome Italy >>> E-mail rob...@un... >>> >>> >>> >>> Il giorno ven 21 set 2018 alle ore 21:00 David Burken <db...@gm...> >>> ha scritto: >>> >>>> Thanks Garrett, look good now... >>>> >>>> >>>> On 9/21/18 1:49 PM, Garrett Potts wrote: >>>> >>>> Hello All: >>>> >>>> >>>> The fix has been checked into the dev baseline. >>>> >>>> >>>> Take care >>>> >>>> Garrett >>>> >>>> >>>> On Sep 21, 2018, at 9:03 AM, David Burken <db...@gm...> wrote: >>>> >>>> Hi Roberta, >>>> >>>> Good catch... >>>> >>>> Garrett, Roberta is correct, you're not picking up the raster type, >>>> i.e. point, area and you're off by 1/2 pixel. >>>> >>>> $ ossim-info -d srtm_47_04.tif >>>> tiff.image0.raster_type: pixel_is_area >>>> tiff.image0.model_pixel_scale: """0.000833333333333333 >>>> 0.000833333333333333 0 """ >>>> tiff.image0.model_tie_point: """0 0 0 49.999583817611 45.0004168845861 >>>> 0 """ >>>> >>>> // Bad: >>>> < image0.geometry.tie_point_lat: 45.0004168845861 >>>> < image0.geometry.tie_point_lon: 49.999583817611 >>>> --- >>>> // Good: >>>> > image0.geometry.tie_point_lat: 45.0000002179194 >>>> > image0.geometry.tie_point_lon: 50.0000004842777 >>>> >>>> You can download some samples at: >>>> >>>> http://srtm.csi.cgiar.org/SELECTION/listImages.asp >>>> >>>> Take care, >>>> Dave >>>> >>>> >>>> >>>> On 9/20/18 4:34 PM, Roberta Ravanelli wrote: >>>> >>>> Hi all, >>>> >>>> I would like to update my Ossim installation from an old version >>>> installed in an old pc (Ubuntu 12.04, libgeotiff 2, opencv 2) to a newer >>>> version in a more recent pc (Ubuntu 16.04, libgeotiff 5, opencv 4). >>>> >>>> I've actually managed to compile and run Ossim in the newer pc, but >>>> I've noticed a different behaviour between the two versions when reading >>>> tiff files. >>>> Specifically, if the tiff has a geographic georeferencing (lon/lat), >>>> there is a wrong shift of a half pixel in the coordinate of the upper left >>>> coordinate (the tiepoint, as it is called in the ossimImageGeometry object). >>>> >>>> Below an example on a SRTM v4_1 tile (srtm_39_03.tif): >>>> >>>> *OLDER VERSION (right: UL as in the .tfw file)* >>>> type: ossimImageGeometry >>>> No transform defined. Using identity transform. >>>> m_projection: >>>> // ossimMapProjection::print >>>> type: ossimEquDistCylProjection >>>> major_axis: 6378137.000000000000000 >>>> minor_axis: 6356752.314199999906123 >>>> origin_latitude: 0.000000000000000 >>>> central_meridian: 0.000000000000000 >>>> origin: ( 0.000000000000000, 0.000000000000000, 0.000, WGE ) >>>> datum: WGE >>>> meters_per_pixel_x: 92.6625433148566 >>>> meters_per_pixel_y: 92.6625433148566 >>>> false_easting_northing: (0,0) >>>> false_easting_northing_units: meters >>>> pcs_code: 4326 >>>> tie_point_xy:* (10.0000002905424,50.0000003631855)* >>>> tie_point_units: degrees >>>> pixel_scale_xy: (0.000833333333333333,0.000833333333333333) >>>> pixel_scale_units: degrees >>>> ossimErrorStatusInterface::print >>>> theErrorStatus: 0 >>>> theErrorStatus string: OSSIM_OK >>>> m_decimationFactors[0]: ( 1.000000000000000, 1.000000000000000 ) >>>> m_imageSize: ( 6001, 6001 ) >>>> m_targetRrds: 0 >>>> >>>> >>>> *Newer version (WRONG: UL as in the .hdr file)* >>>> Type: ossimImageGeometry >>>> No transform defined. Using identity transform. >>>> m_projection: >>>> // ossimMapProjection::print >>>> type: ossimEquDistCylProjection >>>> major_axis: 6378137.000000000000000 >>>> minor_axis: 6356752.314199999906123 >>>> origin_latitude: 0.000000000000000 >>>> central_meridian: 0.000000000000000 >>>> origin: ( 0.000000000000000, 0.000000000000000, 0.000, WGE ) >>>> datum: WGE >>>> meters_per_pixel_x: 92.6625433148566 >>>> meters_per_pixel_y: 92.6625433148566 >>>> false_easting_northing: (0,0) >>>> false_easting_northing_units: meters >>>> pcs_code: 4326 >>>> tie_point_xy: *(9.99958362387572,50.0004170298522)* >>>> tie_point_units: degrees >>>> pixel_scale_xy: (0.000833333333333333,0.000833333333333333) >>>> pixel_scale_units: degrees >>>> ossimErrorStatusInterface::print >>>> theErrorStatus: 0 >>>> theErrorStatus string: OSSIM_OK >>>> m_decimationFactors[0]: ( 1.000000000000000, 1.000000000000000 ) >>>> m_imageSize: ( 6001, 6001 ) >>>> m_targetRrds: 0 >>>> >>>> >>>> I would like to know if such a different behaviour may depend on the >>>> changes performed in Ossim (for instance, I've noticed that the >>>> ossimGeoTiff class is changed a lot beetween the two verions) or maybe >>>> it is caused by the changes in the dependencies (maybe libgeotiff?) >>>> >>>> Can you give me a hint, since I have to abandon the oldest pc? >>>> >>>> Thank you very much! >>>> >>>> Roberta >>>> >>>> >>>> -------------------------------------------------------------------------------- >>>> Roberta Ravanelli, *PhD * >>>> Geodesy and Geomatics Division >>>> University of Rome "La Sapienza" >>>> Via Eudossiana, 18 - 00184 Rome Italy >>>> E-mail rob...@un... >>>> >>>> >>>> >>>> _______________________________________________www.ossim.org >>>> Ossim-developer mailing lis...@li...https://lists.sourceforge.net/lists/listinfo/ossim-developer >>>> >>>> >>>> _______________________________________________ >>>> www.ossim.org >>>> Ossim-developer mailing list >>>> Oss...@li... >>>> https://lists.sourceforge.net/lists/listinfo/ossim-developer >>>> >>>> >>>> >>>> _______________________________________________ >>> www.ossim.org >>> Ossim-developer mailing list >>> Oss...@li... >>> https://lists.sourceforge.net/lists/listinfo/ossim-developer >>> >>> >>> _______________________________________________ >> www.ossim.org >> Ossim-developer mailing list >> Oss...@li... >> https://lists.sourceforge.net/lists/listinfo/ossim-developer >> >> >> _______________________________________________ > www.ossim.org > Ossim-developer mailing list > Oss...@li... > https://lists.sourceforge.net/lists/listinfo/ossim-developer > > > |
