Hi all.
I am trying to extract DSM using KOMPSAT-3 stereo pair on StereoFramework.
KOMPSAT-3 spec. as below :
- Product Lv. : 1R
- Ground Resolution : 0.7m(panchromatic), 2.8(Multispectral, RGB & NIR
Channels)
- Projection : UTM
- with rpc parameters(_rpc.txt) and metadata(.xml)
my desktop env. is Windows 7 64bits, 8GB ram.
1. first, I'm try on OTB inside QGIS (2.6.1)
1) installed OSGeo4W64
2) run QGIS and set option for running orfeotoolbox
- Geoid file : C:/Orfeo/egm96.grd
- OTB applications folder : C:\OSGeo4W64\bin
- OTB command line tools folder : C:\OSGeo4W64\bin
- SRTM tiles folder :
- Active : check
3) set parameters on stereoframework dialog.
- input images list : select kompsat stereo layers, I use NIR channel
images
- default elevation : 200, I found this parameter rpc txt, height offset
- output resolution : 10
- nodata value : -32768 (it set as default)
- method to fuse measure in each DSM cell : max
- parameters estimation modes : fit
- upper left x, upper left y, size x, size y, pixel size x, pixel size y
these parameters is not set, because, as i know, estimation mode is "fit",
these is not necessary
- output cartographic map projection : select wgs
- zone number, northern hemisphere and epsg code is not set
and other parameters are not set except output DSM.
it was failed with otbApplicationLauncherCommandLine.exe error message
including "debug" and "close" buttons. I couldn't read logs, because log
information windows is closed when I close qgis message window.
2. next, I was trying to perform on Monteverdi2, StereoFramework function
with same parameters.
- input images list : select kompsat stereo layers, I use NIR channel
images
- default elevation : 200, I found this parameter rpc txt, height offset
- output resolution : 10
- nodata value : -32768 (it set as default)
- method to fuse measure in each DSM cell : max
- parameters estimation modes : fit to sensor image
- output cartographic map projection : wgs
- block-matching metric : sum of squared distances divided by mean of
block
- minimum altitude offset : 0 meters
- maximum altitude offset : 1000 meters (refer to srtm elevation model)
- others are not set
it was failed.
3. last, StereoFramework was performed on windows command line in OSGeo4W
Shell as :
StereoFramework -input.il
F:\2015_spain\K3_20130205130436_03856_00841310_L1R_N.tif
F:\2015_spain\K3_20130205130557_03856_00841310_L1R_N.tif -elev.geoid
F:\2015_spain\egm96.grd -elev.default 200.0 -output.res 10 -output.nodata
-32768 -output.fusionmethod max -output.out F:\2015_spain\dem_oncmd.tif
-output.mode fit -map wgs -bm.metric ssdmean -bm.minhoffset 0.0
-bm.maxhofffset 1000.0 -ram 1024
2015 Feb 05 14:17:14 : Application.logger (INFO) Elevation management:
setting default height above ellipsoid to 0 meters
2015 Feb 05 14:17:14 : Application.logger (INFO) 1 Stereo couples will
be processed
2015 Feb 05 14:17:15 : Application.logger (INFO) Processing stereo
couple 0 (images index 0 and 1)
Computing epipolar grids...: 0% [
]
Computing epipolar grids...: 2% [*
]
Computing epipolar grids...: 4% [**
]
Computing epipolar grids...: 6% [***
]
Computing epipolar grids...: 8% [****
]
Computing epipolar grids...: 10% [*****
]
Computing epipolar grids...: 12% [******
]
Computing epipolar grids...: 14% [*******
]
Computing epipolar grids...: 16% [********
]
Computing epipolar grids...: 18% [*********
]
Computing epipolar grids...: 20% [**********
]
Computing epipolar grids...: 22% [***********
]
Computing epipolar grids...: 24% [************
]
Computing epipolar grids...: 26% [*************
]
Computing epipolar grids...: 28% [**************
]
Computing epipolar grids...: 30% [***************
]
Computing epipolar grids...: 32% [****************
]
Computing epipolar grids...: 34% [*****************
]
Computing epipolar grids...: 36% [******************
]
Computing epipolar grids...: 38% [*******************
]
Computing epipolar grids...: 40% [********************
]
Computing epipolar grids...: 42% [*********************
]
Computing epipolar grids...: 44% [**********************
]
Computing epipolar grids...: 46% [***********************
]
Computing epipolar grids...: 48% [************************
]
Computing epipolar grids...: 50% [*************************
]
Computing epipolar grids...: 52% [**************************
]
Computing epipolar grids...: 54% [***************************
]
Computing epipolar grids...: 56% [****************************
]
Computing epipolar grids...: 58% [*****************************
]
Computing epipolar grids...: 60% [******************************
]
Computing epipolar grids...: 62% [*******************************
]
Computing epipolar grids...: 64% [********************************
]
Computing epipolar grids...: 66% [*********************************
]
Computing epipolar grids...: 68% [**********************************
]
Computing epipolar grids...: 70% [***********************************
]
Computing epipolar grids...: 72% [************************************
]
Computing epipolar grids...: 74% [*************************************
]
Computing epipolar grids...: 76% [**************************************
]
Computing epipolar grids...: 78% [***************************************
]
Computing epipolar grids...: 80% [****************************************
]
Computing epipolar grids...: 82% [*****************************************
]
Computing epipolar grids...: 84%
[****************************************** ]
Computing epipolar grids...: 86%
[******************************************* ]
Computing epipolar grids...: 88%
[******************************************** ]
Computing epipolar grids...: 90%
[********************************************* ]
Computing epipolar grids...: 92%
[********************************************** ]
Computing epipolar grids...: 94%
[*********************************************** ]
Computing epipolar grids...: 96%
[************************************************ ]
Computing epipolar grids...: 98%
[************************************************* ]
Computing epipolar grids...: 100%
[**************************************************] (0.7 seconds)
2015 Feb 05 14:17:15 : Application.logger (INFO) Minimum disparity : -256
2015 Feb 05 14:17:15 : Application.logger (INFO) Maximum disparity : -0
Inverting left displacement field ...: 0% [
]
Inverting left displacement field ...: 9% [****
]
Inverting left displacement field ...: 19% [*********
]
Inverting left displacement field ...: 29% [**************
]
Inverting left displacement field ...: 39% [*******************
]
Inverting left displacement field ...: 49% [************************
]
Inverting left displacement field ...: 59% [*****************************
]
Inverting left displacement field ...: 69%
[********************************** ]
Inverting left displacement field ...: 79%
[*************************************** ]
Inverting left displacement field ...: 89%
[******************************************** ]
Inverting left displacement field ...: 99%
[************************************************* ]
Inverting left displacement field ...: 100%
[**************************************************] (8e+001 seconds)
2015 Feb 05 14:18:34 : Application.logger (INFO) Using robust SSD Metric
for BlockMatching.
2015 Feb 05 14:18:34 : Application.logger (INFO) Using reverse
block-matching to filter incoherent disparity values.
2015 Feb 05 14:18:34 : Application.logger (INFO) Default elevation set
for Min/Max elevation : 0
2015 Feb 05 14:18:34 : Application.logger (INFO) Minimum elevation found
: 0
2015 Feb 05 14:18:34 : Application.logger (INFO) Maximum elevation found
: 0
2015 Feb 05 14:18:34 : Application.logger (INFO) Disparity mask formula
: if((hdisp > -256) and (hdisp < -0) and (mask>0),255,0)
2015 Feb 05 14:18:34 : Application.logger (INFO) Output DSM parameters
are estimated from sensor input data.
2015 Feb 05 14:18:34 : 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,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]
Writing F:\2015_spain\output_dem.tif...: 0% [
]
it was failed with otbApplicationLauncherCommandLine.exe error message
including "debug" and "close" buttons.
as I know, kompsat-3 relative model is same IKONOS rpc model without any
converting rpc file.
so, I generate rpc having .geom file extention using readImageInfo on cmd
line, and I try again. it was failed!!
how to generate DEM using KOMPSAT-3 stereo pair?
--
--
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/d/optout.
