Pansharpening

Fuse the high spatial detail of a panchromatic image with the spectral content of a lower-resolution multispectral image.

Remote Sensing CoreEnd-User DocumentationTheory + Formula + Parameters

1. What This Feature Does

Fuse the high spatial detail of a panchromatic image with the spectral content of a lower-resolution multispectral image.

Cloud workflow: choose data → configure the scientific method → run → review the resulting layer/report. No programming is required.

2. Recommended Workflow

Select multispectral and higher-resolution panchromatic rasters. ↓ Select multispectral bands and PAN band. ↓ Choose a pansharpening algorithm. ↓ Choose resampling/statistics/detail settings. ↓ Review both sharpness and spectral fidelity.

3. Theory, Methods & Equations

Brovey

Normalizes the panchromatic detail by a synthetic multispectral intensity and injects that detail into each band.

Fi = MSi × PAN / ΣwjMSj
Best used when: high visual sharpness is desired and some spectral distortion is acceptable.

IHS

Transforms three RGB bands to Intensity-Hue-Saturation, replaces/matches intensity using PAN, then converts back.

I ≈ (R+G+B)/3
Best used when: exactly three display bands are being sharpened.

Gram-Schmidt

Creates a synthetic PAN from multispectral bands, orthogonalizes components, replaces the first component with matched PAN, and reconstructs bands.

Best used when: spectral preservation across multiple bands is important.

PCA Pansharpening

Replaces the first principal component with a histogram/statistics-matched PAN, then applies the inverse PCA transform.

C = (1/(n−1))XTX
C vkkvk
Best used when: multiband covariance structure is useful for detail injection.

À trous Wavelet

Decomposes PAN into low-frequency and spatial-detail components without downsampling, then injects detail into multispectral bands.

Detaill = Al−1 − Al
Best used when: multi-scale spatial detail needs controlled injection.

4. Input Data

InputTypeRequirementDescription
Multispectral RasterrasterRequired for selected methodMultiband raster providing the spectral information.
Panchromatic RasterrasterRequired for selected methodHigher-resolution raster providing spatial detail.

5. Parameters Available in the Application

ParameterDefaultChoices / RangeHow to Use It
Pansharpening AlgorithmBrovey, IHS (Intensity-Hue-Saturation), Gram-Schmidt, PCA Pansharpening, Wavelet (à trous)Select the pansharpening algorithm first. Only relevant parameters are shown.
Select Input BandsMultispectral bands to sharpen. Default is All bands. IHS requires exactly three selected bands.
Shown for: brovey, ihs, gram_schmidt, pca, wavelet
Panchromatic Band11 to —
Shown for: brovey, ihs, gram_schmidt, pca, wavelet
Multispectral Resamplingbilinearnearest, bilinear, cubic, lanczos
Shown for: brovey, ihs, gram_schmidt, pca, wavelet
Match PAN StatisticsOnMatch PAN mean and standard deviation to the multispectral intensity before spatial-detail injection.
Shown for: brovey, ihs, gram_schmidt, wavelet
Synthetic PAN WeightsOptional weight per selected multispectral band. Leave empty for equal weights.
Shown for: brovey, gram_schmidt
PCA Fitting Samples2000001000 to 5000000
Shown for: pca
Wavelet Levels21 to 6Number of undecimated à trous B3-spline decomposition levels.
Shown for: wavelet
Wavelet Detail Gain1.00.0 to 5.0
Shown for: wavelet
Random State42
Shown for: pca

6. Output & Interpretation

The result should be interpreted according to the selected method and the physical meaning of the input data. Preserve source units, coordinate reference information, NoData meaning, acquisition date, and preprocessing level when comparing results.

Scientific interpretation: an algorithm can produce a numerically valid result even when the input data are unsuitable. Always check masks, units, sensor characteristics, spatial resolution and reference data.
Important: A sharper image is not automatically more radiometrically accurate. Check spectral indices after pansharpening if quantitative spectral fidelity matters.

7. Best Practices

  • Use analysis-ready inputs and remove invalid/cloud/noise artifacts that are not part of the target phenomenon.
  • Choose parameters from the physical scale of the data, not only from visual appearance.
  • Keep categorical and continuous rasters conceptually separate when selecting interpolation/resampling methods.
  • Compare the result with the source image and independent reference information.
  • Document the settings used when results will be compared across dates, sensors or study areas.