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