Perform P+XS pansharpening
This application performs P+XS pansharpening. The default mode use Pan and XS sensor models to estimate the transformation to superimpose XS over Pan before the fusion (“default mode”). The application provides also a PHR mode for Pleiades images which does not use sensor models as Pan and XS products are already coregistered but only estimate an affine transformation to superimpose XS over the Pan.Note that this option is automatically activated in case Pleiades images are detected as input.
This application has several output images and supports “multi-writing”. Instead of computing and writing each image independently, the streamed image blocks are written in a synchronous way for each output. The output images will be computed strip by strip, using the available RAM to compute the strip size, and a user defined streaming mode can be specified using the streaming extended filenames (type, mode and value). Note that multi-writing can be disabled using the multi-write extended filename option: &multiwrite=false, in this case the output images will be written one by one. Note that multi-writing is not supported for MPI writers.
Input PAN Image
-inp image Mandatory
Input panchromatic image.
Input XS Image
-inxs image Mandatory
Input XS image.
-out image [dtype] Mandatory
This group of parameters allows managing elevation values. Supported formats are SRTM, DTED or any geotiff. DownloadSRTMTiles application could be a useful tool to list/download tiles related to a product.
This parameter allows selecting a directory containing Digital Elevation Model files. Note that this directory should contain only DEM files. Unexpected behaviour might occurs if other images are found in this directory.
-elev.geoid filename [dtype]
Use a geoid grid to get the height above the ellipsoid in case there is no DEM available, no coverage for some points or pixels with no_data in the DEM tiles. A version of the geoid can be found on the OTB website(https://gitlab.orfeo-toolbox.org/orfeotoolbox/otb-data/blob/master/Input/DEM/egm96.grd).
-elev.default float Default value: 0
This parameter allows setting the default height above ellipsoid when there is no DEM available, no coverage for some points or pixels with no_data in the DEM tiles, and no geoid file has been set. This is also used by some application as an average elevation value.
-mode [default|phr] Default value: default
- Default mode
Default superimposition mode : uses any projection reference or sensor model found in the images
- Pleiades mode
Pleiades superimposition mode, designed for the case of a P+XS bundle in SENSOR geometry. It uses a simple transform on the XS image : a scaling and a residual translation.
-method [rcs|lmvm|bayes] Default value: rcs
Selection of the pan-sharpening method.
Simple RCS Pan sharpening operation.
Local Mean and Variance Matching (LMVM) Pan sharpening.
-method.lmvm.radiusx int Default value: 3
Set the x radius of the sliding window.
-method.lmvm.radiusy int Default value: 3
Set the y radius of the sliding window.
-method.bayes.lambda float Default value: 0.9999
Set the weighting value.
-method.bayes.s float Default value: 1
Set the S coefficient.
Spacing of the deformation field
-lms float Default value: 4
Spacing of the deformation field. Default is 10 times the PAN image spacing.
-interpolator [bco|nn|linear] Default value: bco
This group of parameters allows defining how the input image will be interpolated during resampling.
- Bicubic interpolation
Bicubic interpolation leads to very good image quality but is slow.
- Nearest Neighbor interpolation
Nearest neighbor interpolation leads to poor image quality, but it is very fast.
- Linear interpolation
Linear interpolation leads to average image quality but is quite fast
Radius for bicubic interpolation
-interpolator.bco.radius int Default value: 2
This parameter allows controlling the size of the bicubic interpolation filter. If the target pixel size is higher than the input pixel size, increasing this parameter will reduce aliasing artifacts.
-fv float Default value: 0
Fill value for area outside the reprojected image
Available RAM (MB)
-ram int Default value: 256
Available memory for processing (in MB).
From the command-line:
otbcli_BundleToPerfectSensor -inp QB_Toulouse_Ortho_PAN.tif -inxs QB_Toulouse_Ortho_XS.tif -out BundleToPerfectSensor.png uchar
import otbApplication app = otbApplication.Registry.CreateApplication("BundleToPerfectSensor") app.SetParameterString("inp", "QB_Toulouse_Ortho_PAN.tif") app.SetParameterString("inxs", "QB_Toulouse_Ortho_XS.tif") app.SetParameterString("out", "BundleToPerfectSensor.png") app.SetParameterOutputImagePixelType("out", 1) app.ExecuteAndWriteOutput()