gdal_pansharpen - Perform a pansharpen operation.
(Since GDAL 2.1)
gdal_pansharpen [--help-general] pan_dataset {spectral_dataset[,band=num]}+ out_dataset
[-of format] [-b band]* [-w weight_val]*
[-r {nearest,bilinear,cubic,cubicspline,lanczos,average}]
[-threads {ALL_CPUS|number}] [-bitdepth val] [-nodata val]
[-spat_adjust {union,intersection,none,nonewithoutwarning}]
[-co NAME=VALUE]* [-q]
The gdal_pansharpen.py script performs a pan-sharpening operation.
It can create a 'classic' output dataset (such as GeoTIFF), or a VRT dataset
describing the pan-sharpening operation.
More details can be found in the VRT tutorial.
- -of
format:
- Select the output format. Starting with GDAL 2.3, if not specified, the
format is guessed from the extension (previously was GTiff). Use the short
format name. 'VRT' can also be used.
- -b
band:
- Select band band from the input spectral bands for output. Bands
are numbered from 1 in the order spectral bands are specified. Multiple
-b switches may be used. When no -b switch is used, all input
spectral bands are set for output.
- -w
weight_val:
- Specify a weight for the computation of the pseudo panchromatic value.
There must be as many -w switches as input spectral bands.
- -r
{nearest,bilinear,cubic
(default),cubicspline,lanczos,average}:
- Select a resampling algorithm.
- -threads
{ALL_CPUS,number}:
- Specify number of threads to use to do the resampling and pan-sharpening
itself. Can be an integer number or ALL_CPUS.
- -bitdepth
val:
- Specify the bit depth of the panchromatic and spectral bands (e.g. 12). If
not specified, the NBITS metadata item from the panchromatic band will be
used if it exists.
- -nodata
val:
- Specify nodata value for bands. Used for the resampling and pan-sharpening
computation itself. If not set, deduced from the input bands, provided
they have a consistent setting.
- -spat_adjust
{union
(default),intersection,none,nonewithoutwarning}:
- Select behaviour when bands have not the same extent. See
SpatialExtentAdjustment documentation in VRT tutorial
- -co
'NAME=VALUE':
- Passes a creation option to the output format driver. Multiple -co
options may be listed. See format specific documentation for legal
creation options for each format.
- -q:
- Suppress progress monitor and other non-error output.
- pan_dataset
- Dataset with panchromatic band (first band will be used).
- spectral_dataset[,band=num]
- Dataset with one or several spectral bands. If the band option is not
specified, all bands of the datasets are taken into account. Otherwise,
only the specified (num)th band. The same dataset can be repeated several
times.
- out_dataset
- Output dataset
Bands should be in the same projection.
With spectral bands in a single dataset :
gdal_pansharpen.py panchro.tif multispectral.tif pansharpened_out.tif
With a few spectral bands from a single dataset, reordered :
gdal_pansharpen.py panchro.tif multispectral.tif,band=3 multispectral.tif,band=2 multispectral.tif,band=1 pansharpened_out.tif
With spectral bands in several datasets :
gdal_pansharpen.py panchro.tif band1.tif band2.tif band3.tif pansharpened_out.tif
Specify weights:
gdal_pansharpen.py -w 0.7 -w 0.2 -w 0.1 panchro.tif multispectral.tif pansharpened_out.tif
Specify RGB bands from a RGBNir multispectral dataset while
computing the pseudo panchromatic intensity on the 4 RGBNir bands:
gdal_pansharpen.py -b 1 -b 2 -b 3 panchro.tif rgbnir.tif pansharpened_out.tif
Even Rouault <even.rouault at spatialys.com>