Re: [OSSIM] Weird behaviour when reading tiff files (geographic referencing)

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 422-6466

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... <mailto:rob...@un...>
> 
> 
> Il giorno sab 22 set 2018 alle ore 22:53 Garrett Potts <gar...@gm... <mailto: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... <mailto: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 <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 <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... <mailto:rob...@un...>
>> 
>> 
>> Il giorno sab 22 set 2018 alle ore 15:12 Garrett Potts <gar...@gm... <mailto: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... <mailto: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... <mailto:rob...@un...>
>>> 
>>> 
>>> Il giorno ven 21 set 2018 alle ore 21:00 David Burken <db...@gm... <mailto: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... <mailto: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 <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... <mailto:rob...@un...>
>>>>>> 
>>>>>> 
>>>>>> 
>>>>>> _______________________________________________
>>>>>> www.ossim.org <http://www.ossim.org/>
>>>>>> Ossim-developer mailing list
>>>>>> Oss...@li... <mailto:Oss...@li...>
>>>>>> https://lists.sourceforge.net/lists/listinfo/ossim-developer <https://lists.sourceforge.net/lists/listinfo/ossim-developer>
>>>>> 
>>>>> _______________________________________________
>>>>> www.ossim.org <http://www.ossim.org/>
>>>>> Ossim-developer mailing list
>>>>> Oss...@li... <mailto:Oss...@li...>
>>>>> https://lists.sourceforge.net/lists/listinfo/ossim-developer <https://lists.sourceforge.net/lists/listinfo/ossim-developer>
>>>> 
>>> 
>>> _______________________________________________
>>> www.ossim.org <http://www.ossim.org/>
>>> Ossim-developer mailing list
>>> Oss...@li... <mailto:Oss...@li...>
>>> https://lists.sourceforge.net/lists/listinfo/ossim-developer <https://lists.sourceforge.net/lists/listinfo/ossim-developer>
>> 
>> _______________________________________________
>> www.ossim.org <http://www.ossim.org/>
>> Ossim-developer mailing list
>> Oss...@li... <mailto:Oss...@li...>
>> https://lists.sourceforge.net/lists/listinfo/ossim-developer <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