450 likes | 693 Vues
Spatial and Spectral Evaluation of Image Fusion Methods. Sascha Klonus Manfred Ehlers Institute for Geoinformatics and Remote Sensing University of Osnabrück. Content. Introduction Image Fusion Test Site Fusion Results Color Distortions Evaluation Methods and Results Ehlers Fusion
E N D
Spatial and Spectral Evaluation of Image Fusion Methods Sascha Klonus Manfred Ehlers Institute for Geoinformatics and Remote Sensing University of Osnabrück
Content • Introduction Image Fusion • Test Site • Fusion Results • Color Distortions • Evaluation Methods and Results • Ehlers Fusion • Conclusions and Future Work
Data Fusion: Why is it Necessary? • Remote sensors have different spatial resolution for panchromatic and multispectral imagery • The ratios vary between 1:2 and 1:5 • For multisensor fusion the ratios can exceed 1:30(e.g. Ikonos/Landsat)
Objectives of Image Fusion • Sharpen images • Improve geometric corrections • Provide stereo-viewing capabilities • Enhance certain features • Complement data sets • Detect changes • Substitute missing information • Replace defective data Pohl & van Genderen (1998)
Meaning of Pan-Sharpening Spatial + Spectral panchromatic & high geometric resolution multi-/hyperspectral image & low geometric resolution multi-/hyperspectral & high geometric resolution
Fusion Methods • Color Transformations • Modified IHS Transformation • Statistical Methods • Principal Component Merge • Numerical Methods • Brovey • CN Spectral Sharpening • Gram-Schmidt Spectral Sharpening • Wavelet based Fusion • Combined Methods • Ehlers Fusion
Original Data Quickbird Panchromatic image (2004-09-04) Quickbird Multispectral image (2004-09-04) Formosat Multispectral image (2004-01-30) Ikonos Multispectral image (2005-08-03)
Single Sensor Fusion: Quickbird Fused with CN Spectral Sharpening Fused with Ehlers Fused with Brovey Fused with Gram-Schmidt Quickbird Multispectral image Fused with Wavelet Fused with modified IHS Fused with PC
Multisensor Fusion: Ikonos Fused with Ehlers Fused with CN Spectral Sharpening Fused with Brovey Fused with Gram-Schmidt Fused with PC Fused with Wavelet Fused with modified IHS Ikonos Multispectral image
Multisensor Fusion: Formosat Fused with Ehlers Fused with PC Fused with modified IHS Fused with Gram-Schmidt Fused with Brovey Fused with CN Spectral Sharpening Fused with Wavelet Formosat Multispectral image
Fusion Problem: Color Distortion • Panchromatic band has a different spectral sensitivity • Multisensoral differences (e.g. Ikonos and SPOT merge) • Multitemporal (seasonal) changes between pan and ms image data • Inconsistent panchromatic information is fused into the multispectral bands
Spectral Comparison Methods (1) • RMSE s = standard deviation org = Original image fused = Fused image x = Mean • Correlation coefficients • Visual (Structure and Colour Preservation)
Spectral Comparison Methods (2) Per Pixel Deviation Fused image (Formosat 2m) Degraded to ground resolution of original image(Formosat = 8m) Result: Vector containing the deviation per pixel Degrade Original multispectral image (Formosat 8m)
Spatial Comparison Methods (1) - Edge Detection - -
Spatial Comparison Methods (2) Highpass Filtering Correlation
FFT Filter Based Data Fusion (Ehlers Fusion) FFT FFT LPF HPF FourierSpectrum PanHP FFT-1 FourierSpectrum ILP ILP+PanHP H S R‘ G‘ B‘ I H S R G B IHS-1 Basis: IHS Transform and Filtering in the Fourier Domain Multispectral Image Panchromatic Image
Panchromatic image and its spectrum Original panchromatic image Panchromatic Spectrum
Filtersetting effects Intensity Frequency fn Cut-off Frequency Filtered Panchromatic Spectrum
Effects in the spatial domain Filtered panchromatic image Fused image
Filtersetting effects Intensity Frequency fn Cut-off Frequency Filtered Panchromatic Spectrum
Effects in the spatial domain Filtered panchromatic image Fused image
Filtersetting effects Intensity Frequency fn Cut-off Frequency Filtered Panchromatic Spectrum
Effects in the spatial domain Filtered panchromatic image Fused image
Results • Ehlers Fusion shows the best overall results in all images • It works also if the panchromatic Information does not match the spectral sensitivity of the merged bands (multitemporal and multisensoral fusion) • Its performance is superior to standard fusion techniques (IHS, Brovey Transform, PC Merge) • Wavelet preserves the spectral characteristics at the cost of spatial improvement • Ehlers Fusion is integrated in a commercial image processing system (Erdas Imagine 9.1)
Future Work • Fusion of radar- and optical Data • Development of one method to evaluate the spatial and spectral quality of an fused image • Comparison with the algorithm of Zhang (PCI Geomatica) • Research on automation for filter design
Thanks for your Attention Questions???
Multispectral image and its spectrum Original multispectral intensity Multispectral intensity spectrum
Filtersetting effects Intensity Frequency fn Cut-off Frequency Filtered multispectral spectrum
Filtersetting effects Intensity Frequency fn Cut-off Frequency Filtered multispectral spectrum
Filtersetting effects Intensity Frequency fn Cut-off Frequency Filtered multispectral spectrum