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.
