Menu
▾
▴
Re: [OSSIM] Weird behaviour when reading tiff files (geographic referencing)
|
From: Roberta R. <rob...@un...> - 2018-09-22 21:54:13
|
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 > > > |
