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] <javascript:>> > >> >> 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/forum/#!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_StereoRectificationGridGenerator >> >> -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]<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 >> --- >> 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. >> > > -- -- 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.
