Hi when using Grid generator something appears weird in the logs epi.rectsizex: 1024 epi.rectsizey: 1024
*epi.baseline: 0* and stereoframework 2014 Jan 10 16:03:39 : Application.logger (INFO) Minimum disparity : -0 2014 Jan 10 16:03:39 : Application.logger (INFO) Maximum disparity : 0 this value is related to transfomation of horizontal disparities when working in epipolar geometry to physical elevation 2014/1/10 Jonathan Guinet <[email protected]> > Hello, > > your couple origin is quite distant could it be possible to relaunch the > process with extract with the same footprint (use extractROI for > band128B.tif with band128A as ref : using option -mode fit -mode.fit.ref ) > > Stereoframework and stereoreconstruction will give better result using > SRTM data instead of elevation default (DEM data can be easily found using > DownloadSTRMTile application) > > > cheer > > Jonathan > > > 2014/1/10 Sunny Depp <[email protected]> > >> Hi, >> >> If i skip that option and run the command all i get is.... >> >> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_StereoFramework >> -input.il"band128A >> .tif" "band128F.tif" -output.res 2 -output.out "dem.tif" -elev.default >> 100 -pos >> tproc.med 1 -stereorect.fwdgridstep 8 -stereorect.invgridssrate 4 >> -bm.maxhoffset >> 20 -bm.minhoffset -20 -bm.radius 3 >> >> 2014 Jan 10 16:03:38 : Application.logger (INFO) Elevation management: >> setting default height above ellipsoid to 100 meters >> 2014 Jan 10 16:03:38 : Application.logger (INFO) 1 Stereo couples >> will be processed >> 2014 Jan 10 16:03:38 : Application.logger (INFO) Processing stereo >> couple 0 (images index 0 and 1) >> Computing epipolar grids...: 100% >> [**************************************************] (0.2 seconds) >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Minimum disparity : -0 >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Maximum disparity : 0 >> Inverting left deformation field ...: 100% >> [**************************************************] (0.05 seconds) >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Using robust SSD >> Metric for BlockMatching. >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Using reverse >> block-matching to filter incoherent disparity values. >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Default elevation set >> for Min/Max elevation : 100 >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Minimum elevation >> found : 100 >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Maximum elevation >> found : 100 >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Disparity mask >> formula : if((hdisp > -0) and (hdisp < 0) and (mask>0),255,0) >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Output DSM parameters >> are estimated from sensor input data. >> 2014 Jan 10 16:03:39 : Application.logger (INFO) Output DSM projection >> reference: >> GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS >> 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTH >> >> ORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]] >> 2014 Jan 10 16:03:40 : Application.logger (FATAL) The following error >> occurred during application execution : bad allocation >> >> On Friday, January 10, 2014 3:49:34 PM UTC+5:30, Jonathan Guinet wrote: >> >>> Hi, >>> >>> in your test command you explicitly ignore by using extended filaname >>> [1] skipcarto >>> >>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_StereoFramework >>> -input.il"band128A.tif? >>> *&skipcarto=true*" >>> >>> where your geo info seems to be correct , coouèld it be possible to >>> relaunch without this option >>> >>> [1] >>> http://wiki.orfeo-toolbox.org/index.php/ExtendedFileName<http://www.google.com/url?q=http%3A%2F%2Fwiki.orfeo-toolbox.org%2Findex.php%2FExtendedFileName&sa=D&sntz=1&usg=AFQjCNFCenNohDwZwjbyIiu9-5SOhk1ssQ> >>> >>> >>> 2014/1/10 Sunny Depp <[email protected]> >>> >>> Hi Jonathan, >>>> >>>> and if i reduce the size to 128x128 for both images and run the >>>> stereoframework as >>>> >>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_StereoFramework >>>> -input.il"band128A.tif?&skipcarto=true" "band128F.tif?&skipcarto=true" >>>> -output.res >>>> 2 -o >>>> utput.out "dem.tif" -elev.default 100 -postproc.med 1 >>>> -stereorect.fwdgridstep 8 -stereorect.invgridssrate 4 -bm.maxhoffset 20 >>>> -bm.minhoffset -20 -bm.r >>>> adius 3 >>>> >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Elevation >>>> management: setting default height above ellipsoid to 100 meters >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) 1 Stereo couples >>>> will be processed >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Processing stereo >>>> couple 0 (images index 0 and 1) >>>> Computing epipolar grids...: 100% >>>> [**************************************************] >>>> (0.2 seconds) >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Minimum disparity : >>>> -0 >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Maximum disparity : >>>> 0 >>>> Inverting left deformation field ...: 100% >>>> [**************************************************] (0.05 seconds) >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Using robust SSD >>>> Metric for BlockMatching. >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Using reverse >>>> block-matching to filter incoherent disparity values. >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Default elevation >>>> set for Min/Max elevation : 100 >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Minimum elevation >>>> found : 100 >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Maximum elevation >>>> found : 100 >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Disparity mask >>>> formula : if((hdisp > -0) and (hdisp < 0) and (mask>0),255,0) >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Output DSM >>>> parameters are estimated from sensor input data. >>>> 2014 Jan 10 15:38:07 : Application.logger (INFO) Output DSM >>>> projection reference: >>>> GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS >>>> 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]], >>>> AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTH >>>> ORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433, >>>> AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]] >>>> Writing dem.tif...: 0% [ >>>> ] >>>> >>>> >>>> I get this and dem file is writing but i dont see any progress here but >>>> the file size is increasing and going on till 32GB.. i felt somewhat >>>> strange so stopped the execution. can u clarify me why this is happening ? >>>> >>>> >>>> Thanks >>>> >>>> On Wednesday, January 8, 2014 1:51:33 PM UTC+5:30, Jonathan Guinet >>>> wrote: >>>> >>>>> Hi Sandeep >>>>> >>>>> could it be possible to send image information using for example >>>>> otbcli_ReadImageInfo -in BANDA100.tif -outkwl BANDA 100.geom >>>>> >>>>> thanks >>>>> >>>>> >>>>> >>>>> >>>>> 2014/1/8 Sunny Depp <[email protected]> >>>>> >>>>>> >>>>>> Hi Everyone, >>>>>> >>>>>> I'm new to the orfeo toolbox , installed directly from OSGEO4W I'm >>>>>> trying to extract a DTM from cartosat stereo pair. However, I'm running >>>>>> into a bit of trouble down the line. I'm following the OTB Cookbook >>>>>> example >>>>>> for VHR images. >>>>>> >>>>>> i had followed this thread (https://groups.google.com/for >>>>>> um/#!searchin/otb-users/stereo/otb-users/kocntH7LjO8/ZIgLpI7EFJgJ) and >>>>>> successfully generated DEM from the inputs which i downloaded from OTB >>>>>> Data, but for the Cartosat stereo pair which i have, i am doing the >>>>>> stereoframework on this pair i am getting error as: >>>>>> >>>>>> Using ExtractROI i had extracted images to size 1024 from original >>>>>> BANDA and BANDF images and applied stereoframework on these images. >>>>>> >>>>>> >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_StereoFramework - >>>>>> input.il "BANDA100.tif" "BANDF100.tif" -elev.default 50 -output.out >>>>>> "dem.tif" -output.res 20 -postproc.med 1 -stereorect.fwdgridstep 8 >>>>>> -stereorect.invgridssrate 4 >>>>>> >>>>>> 2014 Jan 08 12:48:14 : Application.logger (INFO) Elevation >>>>>> management: setting default height above ellipsoid to 50 meters >>>>>> 2014 Jan 08 12:48:14 : Application.logger (INFO) 1 Stereo couples >>>>>> will be processed >>>>>> 2014 Jan 08 12:48:14 : Application.logger (INFO) Processing stereo >>>>>> couple 0 (images index 0 and 1) >>>>>> Computing epipolar grids...: 100% >>>>>> [**************************************************] >>>>>> (0.2 seconds) >>>>>> 2014 Jan 08 12:48:14 : Application.logger (INFO) Minimum disparity >>>>>> : -0 >>>>>> 2014 Jan 08 12:48:14 : Application.logger (INFO) Maximum disparity >>>>>> : 0 >>>>>> Inverting left deformation field ...: 100% >>>>>> [**************************************************] (5e+001 seconds) >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Using robust SSD >>>>>> Metric for BlockMatching. >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Using reverse >>>>>> block-matching to filter incoherent disparity values. >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Default elevation >>>>>> set for Min/Max elevation : 50 >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Minimum elevation >>>>>> found : 50 >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Maximum elevation >>>>>> found : 50 >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Disparity mask >>>>>> formula : if((hdisp > -0) and (hdisp < 0) and (mask>0),255,0) >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Output DSM >>>>>> parameters are estimated from sensor input data. >>>>>> 2014 Jan 08 12:49:08 : Application.logger (INFO) Output DSM >>>>>> projection reference: >>>>>> GEOGCS["WGS84",DATUM["WGS_1984",SPHEROID["WGS84",6378137, >>>>>> 298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]], >>>>>> PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0. >>>>>> 0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]] >>>>>> 2014 Jan 08 12:49:10 : Application.logger (FATAL) The following >>>>>> error occurred during application execution : bad allocation >>>>>> >>>>>> >>>>>> So i had followed this http://orfeo-toolbox.org/ >>>>>> CookBook/CookBooksu41.html#tailCookBooksu41.html step by step. still >>>>>> i am getting error at otbcli_DisparityMapToElevationMap step >>>>>> >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_StereoRectifica >>>>>> tionGridGenerator -io.inleft "BANDA100.tif?&geom=BANDA100.geom" >>>>>> -io.inright "BANDF100.tif? >>>>>> &geom=BANDF100.geom"-epi.elevation.dem srtm_directory\ >>>>>> -epi.elevation.geoid egm96.grd -epi.step 20 -io.outleft SRGG_A.tif >>>>>> -io.outright SRGG_F.tif >>>>>> >>>>>> 2014 Jan 08 01:00:47 : Application.logger (INFO) Elevation >>>>>> management: setting default height above ellipsoid to 0 meters >>>>>> 2014 Jan 08 01:00:48 : Application.logger (INFO) Elevation >>>>>> management: using geoid file (egm96.grd) >>>>>> 2014 Jan 08 01:00:48 : Application.logger (INFO) Elevation >>>>>> management: using DEM directory (srtm_directory\) >>>>>> Computing epipolar grids ...: 100% >>>>>> [**************************************************] >>>>>> (0.3 seconds) >>>>>> Writing SRGG_A.tif...: 100% >>>>>> [**************************************************] >>>>>> (0.2 seconds) >>>>>> Writing SRGG_F.tif...: 100% >>>>>> [**************************************************] >>>>>> (0.2 seconds) >>>>>> Output parameters value: >>>>>> epi.rectsizex: 1024 >>>>>> epi.rectsizey: 1024 >>>>>> epi.baseline: 0 >>>>>> >>>>>> >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_GridBasedImageResampling >>>>>> -io.in BANDA100.tif -io.out GBIR_A.tif -grid.in SRGG_A.tif >>>>>> -out.sizex 1024 -out. >>>>>> sizey 1024 >>>>>> >>>>>> Writing GBIR_A.tif...: 100% >>>>>> [**************************************************] >>>>>> (0.9 seconds) >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_GridBasedImageResampling >>>>>> -io.in BANDF100.tif -io.out GBIR_F.tif -grid.in SRGG_F.tif >>>>>> -out.sizex 1024 -out. >>>>>> sizey 1024 >>>>>> >>>>>> Writing GBIR_F.tif...: 100% >>>>>> [**************************************************] >>>>>> (0.9 seconds) >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_BandMath -il >>>>>> GBIR_A.tif -out BM_A.tif -exp "if(im1b1<=0,0,255)" >>>>>> >>>>>> 2014 Jan 08 02:03:52 : Application.logger (INFO) Image #1 has 1 >>>>>> components >>>>>> >>>>>> Writing BM_A.tif...: 100% >>>>>> [**************************************************] >>>>>> (0.3 seconds) >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_BandMath -il >>>>>> GBIR_F.tif -out BM_F.tif -exp "if(im1b1<=0,0,255)" >>>>>> >>>>>> 2014 Jan 08 02:04:06 : Application.logger (INFO) Image #1 has 1 >>>>>> components >>>>>> >>>>>> Writing BM_F.tif...: 100% >>>>>> [**************************************************] >>>>>> (0.3 seconds) >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_BlockMatching >>>>>> -io.inleft GBIR_A.tif -io.inright GBIR_F.tif -io.out >>>>>> disparity_map_ncc.tif >>>>>> -bm.minhd -45 -b >>>>>> m.maxhd 5 -bm.minvd 1 -bm.maxvd 1 -mask.inleft BM_A.tif -mask.inright >>>>>> BM_F.tif -io.outmetric 1 -bm.metric ncc -bm.subpixel dichotomy >>>>>> -bm.medianfilter. >>>>>> radius 5 -bm.medianfilter.incoherence 2.0 >>>>>> >>>>>> 2014 Jan 08 02:10:45 : Application.logger (INFO) Masking criterion >>>>>> on left image: if(inmask > 0,255,0) >>>>>> 2014 Jan 08 02:10:45 : Application.logger (INFO) Masking criterion >>>>>> on right image: if(inmask > 0,255,0) >>>>>> 2014 Jan 08 02:10:45 : Application.logger (INFO) Applies a median >>>>>> filter to smooth the disparity map >>>>>> >>>>>> NCC block matching: 100% >>>>>> [**************************************************] >>>>>> (3e+001 seconds) >>>>>> Sub-pixel refinement: 100% >>>>>> [**************************************************] >>>>>> (9e+001 seconds) >>>>>> Writing disparity_map_ncc.tif...: 100% >>>>>> [**************************************************] >>>>>> (1e+002 seconds) >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_BandMath -il >>>>>> disparity_map_ncc.tif -out thres_hdisparity.tif uint8 -exp >>>>>> "if(im1b3>0.9,255,0)" >>>>>> >>>>>> 2014 Jan 08 02:15:03 : Application.logger (INFO) Image #1 has 3 >>>>>> components >>>>>> >>>>>> Writing thres_hdisparity.tif...: 100% >>>>>> [**************************************************] >>>>>> (0.4 seconds) >>>>>> >>>>>> >>>>>> C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_DisparityMapToElevationMap >>>>>> -io.in disparity_map_ncc.tif -io.left BANDA100.tif -io.right >>>>>> BANDF100.tif -io.lgrid SRGG_A.tif -io.rgrid SRGG_F.tif -io.mask >>>>>> thres_hdisparity.tif -io.out disparity_map_ssd_to_elevation.tif >>>>>> -hmin 10 -hmax 400 -elev.default 50 >>>>>> >>>>>> 2014 Jan 08 02:26:29 : Application.logger (INFO) Elevation >>>>>> management: setting default height above ellipsoid to 50 meters >>>>>> 2014 Jan 08 02:26:29 : Application.logger (INFO) Elevation map >>>>>> origin : [271886,2.16958e+006] >>>>>> >>>>>> 2014 Jan 08 02:26:29 : Application.logger (INFO) Elevation map >>>>>> size : [19966642,51874882] >>>>>> >>>>>> 2014 Jan 08 02:26:30 : Application.logger (FATAL) The following >>>>>> error occurred during application execution : bad allocation >>>>>> >>>>>> >>>>>> How do i interpret this error? and does this OTB work on Cartosat >>>>>> Stereo images. I am struck in this. >>>>>> >>>>>> >>>>>> Thanks in advance, >>>>>> Sandeep. >>>>>> >>>>>> -- >>>>>> -- >>>>>> 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. >>>>>> >>>>> >>>>> -- >>>> -- >>>> 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. >>>> >>> >>> -- >> -- >> 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. >> > > -- -- 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.
