Abstract
In this work, we propose applying a non-linear strategy to estimate the novel synthetic intensity component through a second-order polynomial regression analysis, such that it could be expressed as a linear combination of multispectral bands in addition to their weighted joint products. Details maps are, therefore, obtained as the residual between the panchromatic (PAN) and the novel synthetic intensity component. Indeed, the high-resolution Multispectral (MS) image is produced by transferring the details into the expanded multispectral bands. The quality improvements achievable by our proposal are evaluated by using two high resolution data sets collected by the WorldView-3 and WorldView-4 sensors. The performance of our proposed technique is clearly shown against the existing state-of-the-art pansharpening methods (especially with linear schemes) and particularly for multispectral channels with non-overlapping spectral wavelengths with respect to the PAN image.