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] 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.
