Hi Jonathan
Sorry for the late reply.these are the image information for two images.
Thanks
C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_ReadImageInfo -in BANDA100.tif
-outkwl BANDA100.geom
2014 Jan 10 14:57:50 : Application.logger (INFO)
Image general informations:
Number of bands : 1
Start index : [0,0]
Size : [1024,1024]
Origin : [271886,2169578]
Spacing : [2.5,-2.5]
Estimated ground spacing (in meters): [2.49388,2.51071]
Image acquisition informations:
Sensor :
Image identification number:
Image projection : PROJCS["WGS 84 / UTM zone 45N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",87],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32645"]]
Image default RGB composition:
[R, G, B] = [0,1,2]
Ground control points information:
Number of GCPs = 0
GCPs projection =
Output parameters value:
indexx: 0
indexy: 0
sizex: 1024
sizey: 1024
spacingx: 2.5
spacingy: -2.5
originx: 271886.4063
originy: 2169578.75
estimatedgroundspacingx: 2.493876696
estimatedgroundspacingy: 2.510709524
numberbands: 1
sensor:
id:
time:
ullat: 0
ullon: 0
urlat: 0
urlon: 0
lrlat: 0
lrlon: 0
lllat: 0
lllon: 0
town:
country:
rgb.r: 0
rgb.g: 1
rgb.b: 2
projectionref: PROJCS["WGS 84 / UTM zone 45N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",87],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32645"]]
keyword:
gcp.count: 0
gcp.proj:
gcp.ids:
gcp.info:
gcp.imcoord:
gcp.geocoord:
C:\Users\Sandeep\Desktop\OTB\PRODUCT1>otbcli_ReadImageInfo -in BANDF100.tif
-outkwl BANDF100.geom
2014 Jan 10 14:58:37 : Application.logger (INFO)
Image general informations:
Number of bands : 1
Start index : [0,0]
Size : [1024,1024]
Origin : [270621,2169808]
Spacing : [2.5,-2.5]
Estimated ground spacing (in meters): [2.49784,2.51071]
Image acquisition informations:
Sensor :
Image identification number:
Image projection : PROJCS["WGS 84 / UTM zone 45N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",87],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32645"]]
Image default RGB composition:
[R, G, B] = [0,1,2]
Ground control points information:
Number of GCPs = 0
GCPs projection =
Output parameters value:
indexx: 0
indexy: 0
sizex: 1024
sizey: 1024
spacingx: 2.5
spacingy: -2.5
originx: 270621.4063
originy: 2169808.75
estimatedgroundspacingx: 2.497842789
estimatedgroundspacingy: 2.510709524
numberbands: 1
sensor:
id:
time:
ullat: 0
ullon: 0
urlat: 0
urlon: 0
lrlat: 0
lrlon: 0
lllat: 0
lllon: 0
town:
country:
rgb.r: 0
rgb.g: 1
rgb.b: 2
projectionref: PROJCS["WGS 84 / UTM zone 45N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",87],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32645"]]
keyword:
gcp.count: 0
gcp.proj:
gcp.ids:
gcp.info:
gcp.imcoord:
gcp.geocoord:
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.
