Sorry to double post on this thread but I still couldn't find a solution for my problem concerning the use of a Pleiades mosaic produced with Monteverdi from Pleiades tiles with otbcli_HomologousPointsExtraction.
As I told in the previous post, results from otbcli_ReadImageInfo are a bit strange: R1C1 tile: Estimated ground spacing (in meters): [0.604252,0.544926] ullat: 36.26757812 ullon: 27.97968864 urlat: 36.27081299 urlon: 28.18044281 lrlat: 36.0558815 lrlon: 28.17984772 lllat: 36.05049896 lllon: 27.98065186 mosaic: Estimated ground spacing (in meters): [342.499,49.2681] ullat: 0 ullon: 0 urlat: 0 urlon: 0 lrlat: 0 lrlon: 0 lllat: 0 lllon: 0 Other informations (Size, Origin, ...) are the same. I wonder wether I did a mistake or something. Thank you, Arthur Le lundi 21 juillet 2014 16:11:08 UTC+2, arthurd a écrit : > > Hello, > > I bring back this discussion, as I have an issue with the mosaic module in > Monteverdi (I think), still with Pleiades TIF data. > > otbcli_HomologousPointsExtraction doesn't find any point when I use the 2 > mosaics generated with Monteverdi. On the other hand, when I try the same > command with the R1C1 tiles directly (which represent 4/5 of the data, R1C2 > being the last 1/5), it finds a lot of points. > > otbcli_HomologousPointsExtraction -in1 IMG1_mos.tif -in2 IMG2_OTB.tif > -algorithm sift -mode geobins -mode.geobins.binstep 1024 > -mode.geobins.binsize 512 -precision 21 -elev.dem srtm/ -elev.geoid > egm96.grd -out ptsHomol_test.txt -mfilter 1 -backmatching 1 > > The positions of the 2 regions considered at each step are really strange > sometimes (very distant, like 900px from where it should find an homologous > point). When I use otbcli_ReadImageInfo on a mosaic, I obtain an Estimated > ground spacing (in meters) of [80.8625,53.4112], where it should be near > 50cm. Maybe that's the source of my troubles. > > To open the Pleiades tiles, I used the "Real image" Data type in the "Open > dataset" window. I also tried with "Pleiades image", but I can't find a way > to export the mosaic... With the "Real image" way, I can right click on the > mosaic to export it in a file, but with the "Pleiades image" way, the > option is not available. So maybe the metadata would be right, but I can't > produce it. > > Any help would be greatly appreciated! > > Thank you, > > Arthur > > Le mardi 28 mai 2013 11:43:27 UTC+2, arthurd a écrit : >> >> Dear Julien, >> >> You're right, the mosaic module is working, I made a confusion between 2 >> images when I checked on ENVI, my bad... >> >> I tried to make gdalbuildvrt work after removing the TFW files (and also >> the -separate option, as it seems that it's for multi-band datasets), >> but I still got an error: >> >> 0...10...20..Warning 6: gdalbuildvrt does not support ungeoreferenced >> image. Skipping IMG_PHR1A_P_201304200754474_SEN_617691101-001_R1C1.TIF >> >> Anyway, I think that the mosaic module will be fine for me, thanks again >> for the help! >> >> Arthur >> >> Le 28/05/13 11:08, Julien Michel a �crit : >> > Dear Arthur, >> > >> > Two comments : >> > - The mosaic module should work with TIF Pleiades ... Did you add all >> > tiles to the module ? >> > - GdalBuildVRT will not be able to process sensor geometry image, and >> it >> > is confused by the TFW files (I suggest you remove them or backup them >> > somewhere else, these TFW file are very confusing for most software). >> > However, the VRT file format is simple xml and should be able to deal >> > with this case ( http://www.gdal.org/gdal_vrttut.html ). >> > >> > Regards, >> > >> > Julien >> > >> > Le 28/05/2013 11:00, arthurd a �crit : >> >> Thanks a lot to both of you for the help. >> >> >> >> I tried to generate a VRT file with gdal, but I get a warning for each >> >> one of my tiles: >> >> >> >> 0...10...20..Warning 6: gdalbuildvrt does not support rotated geo >> >> transforms. Skipping >> >> IMG_PHR1A_P_201304200754474_SEN_617691101-001_R1C1.TIF >> >> >> >> And no result in the end. gdalinfo's output of the first tile is >> showing: >> >> >> >> ... >> >> GeoTransform = >> >> 40.5539979235211, 4.61860309e-06, -7.271299999999999e-10 >> >> 12.61045455273367, 1.518228e-08, -4.53298586e-06 >> >> ... >> >> >> >> Then I tried to read directly a TIF file with otbcli_ReadImageInfo and >> >> it works, but I need to read the 4 tiles as a whole. Is it possible? >> >> >> >> I also gave a look at Monteverdi (v1.14.0) and found the mosaic >> >> option. After opening the 4 tiles separately, (I tried with both >> >> "Pleiades image" and "real image" options), I successfuly imported >> >> them into the mosaic module and it generated a mosaic output. But when >> >> I display this result, all I can see is the first tile. The following >> >> message appears on the Terminal when I try to display the mosaic: >> >> >> >> libGL error: failed to load driver: swrast >> >> libGL error: Try again with LIBGL_DEBUG=verbose for more details. >> >> >> >> Arthur >> >> >> >> Le mardi 28 mai 2013 09:15:39 UTC+2, Julien Michel a �crit : >> >> >> >> Dear Arthur, >> >> >> >> The -in option from otbcli_ReadImageInfo awaits an image file, not >> >> an XML one. >> >> >> >> otbcli_ReadImageInfo -in >> >> DIM_PHR1A_P_201304200754474_SEN_617691101-001.TIF -outkwl test >> >> >> >> should work fine, and it will find the XML file by itself. >> >> >> >> Regards, >> >> >> >> Julien >> >> >> >> Le 27/05/2013 10:38, arthurd a �crit : >> >>> Hi everyone, >> >>> >> >>> I read some informations in the CookBook, about how to open >> >>> jpeg2000 Pleiades data, but I didn't find the same info for the >> >>> tif format. >> >>> >> >>> I have a Pleiade image with several tif tiles which I tried to >> >>> open, using otbcli_ReadImageInfo, indicating the DIMAP XML file. >> >>> Here is the output: >> >>> >> >>> otbcli_ReadImageInfo -in >> >>> DIM_PHR1A_P_201304200754474_SEN_617691101-001.XML -outkwl test >> >>> ReadImageInfo -in >> >>> DIM_PHR1A_P_201304200754474_SEN_617691101-001.XML -outkwl test >> >>> >> >>> ERROR 4: Failed to find <Raster_Dimensions> in document. >> >>> ImageFileReader (0x9972e0) >> >>> RTTI typeinfo: otb::ImageFileReader<otb::VectorImage<float, >> 2u> > >> >>> Reference Count: 1 >> >>> Modified Time: 7597 >> >>> Debug: Off >> >>> Observers: >> >>> none >> >>> Number Of Required Inputs: 0 >> >>> Number Of Required Outputs: 1 >> >>> Number Of Threads: 4 >> >>> ReleaseDataFlag: Off >> >>> ReleaseDataBeforeUpdateFlag: Off >> >>> No Inputs >> >>> Output 0: (0x807110) >> >>> AbortGenerateData: Off >> >>> Progress: 0 >> >>> Multithreader: >> >>> RTTI typeinfo: itk::MultiThreader >> >>> Reference Count: 1 >> >>> Modified Time: 7585 >> >>> Debug: Off >> >>> Observers: >> >>> none >> >>> Thread Count: 4 >> >>> Global Maximum Number Of Threads: 128 >> >>> Global Default Number Of Threads: 4 >> >>> ImageIO: (null) >> >>> UserSpecifiedImageIO flag: 0 >> >>> m_FileName: DIM_PHR1A_P_201304200754474_SEN_617691101-001.XML >> >>> m_UseStreaming: 1 >> >>> ImageIO: (null) >> >>> UserSpecifiedImageIO flag: 0 >> >>> m_FileName: DIM_PHR1A_P_201304200754474_SEN_617691101-001.XML >> >>> 2013 May 27 10:15:02 : Application.logger (CRITICAL) Invalid >> >>> image filename DIM_PHR1A_P_201304200754474_SEN_617691101-001.XML. >> >>> 2013 May 27 10:15:02 : Application.logger (FATAL) The >> >>> following error occurred during application execution : >> >>> >> /build/buildd/otb-3.16.0/Code/ApplicationEngine/otbWrapperInputImageParameter.cxx:159: >> >> >> >>> itk::ERROR: InputImageParameter(0x9a7af0): No input image or >> >>> filename detected... >> >>> >> >>> I checked the DIMAP file: <Raster_Dimensions> is present. I >> >>> googled the error and found some results with gdal problems >> >>> (http://trac.osgeo.org/gdal/ticket/4826 >> >>> <http://trac.osgeo.org/gdal/ticket/4826>). So I guess it might >> be >> >>> a reason. I work on Ubuntu 12.04.2 with gdal 1.9. >> >>> >> >>> Arthur >> >>> -- >> >>> -- >> >>> Check the OTB FAQ at >> >>> http://www.orfeo-toolbox.org/FAQ.html >> >>> <http://www.orfeo-toolbox.org/FAQ.html> >> >>> >> >>> You received this message because you are subscribed to the >> Google >> >>> Groups "otb-users" group. >> >>> To post to this group, send email to [email protected] >> >>> <javascript:> >> >>> To unsubscribe from this group, send email to >> >>> [email protected] <javascript:> >> >>> For more options, visit this group at >> >>> http://groups.google.com/group/otb-users?hl=en >> >>> <http://groups.google.com/group/otb-users?hl=en> >> >>> --- >> >>> You received this message because you are subscribed to the >> >>> Google Groups "otb-users" group. >> >>> To unsubscribe from this group and stop receiving emails from it, >> >>> send an email to [email protected] <javascript:>. >> >>> For more options, visit https://groups.google.com/groups/opt_out >> >>> <https://groups.google.com/groups/opt_out>. >> >>> >> >>> >> >> >> >> >> >> -- >> >> Julien MICHEL >> >> CNES - DCT/SI/AP - BPI 1219 >> >> 18, avenue Edouard Belin >> >> 31401 Toulouse Cedex 09 - France >> >> Tel: +33 561 282 894 - Fax: +33 561 283 109 >> >> >> >> -- >> >> -- >> >> Check the OTB FAQ at >> >> http://www.orfeo-toolbox.org/FAQ.html >> >> >> >> You received this message because you are subscribed to the Google >> >> Groups "otb-users" group. >> >> To post to this group, send email to [email protected] >> >> To unsubscribe from this group, send email to >> >> [email protected] >> >> For more options, visit this group at >> >> http://groups.google.com/group/otb-users?hl=en >> >> --- >> >> You received this message because you are subscribed to the Google >> >> Groups "otb-users" group. >> >> To unsubscribe from this group and stop receiving emails from it, send >> >> an email to [email protected]. >> >> For more options, visit https://groups.google.com/groups/opt_out. >> >> >> >> >> > >> > >> > -- >> > Julien MICHEL >> > CNES - DCT/SI/AP - BPI 1219 >> > 18, avenue Edouard Belin >> > 31401 Toulouse Cedex 09 - France >> > Tel: +33 561 282 894 - Fax: +33 561 283 109 >> > >> > -- >> > -- >> > Check the OTB FAQ at >> > http://www.orfeo-toolbox.org/FAQ.html >> > >> > You received this message because you are subscribed to the Google >> > Groups "otb-users" group. >> > To post to this group, send email to [email protected] >> > To unsubscribe from this group, send email to >> > [email protected] >> > For more options, visit this group at >> > http://groups.google.com/group/otb-users?hl=en >> > --- >> > You received this message because you are subscribed to a topic in the >> > Google Groups "otb-users" group. >> > To unsubscribe from this topic, visit >> > >> https://groups.google.com/d/topic/otb-users/7vkOHINal3k/unsubscribe?hl=en. >> >> > To unsubscribe from this group and all its topics, send an email to >> > [email protected]. >> > For more options, visit https://groups.google.com/groups/opt_out. >> > >> > >> > -- -- Check the OTB FAQ at http://www.orfeo-toolbox.org/FAQ.html You received this message because you are subscribed to the Google Groups "otb-users" group. To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/otb-users?hl=en --- You received this message because you are subscribed to the Google Groups "otb-users" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
