Tuesday, August 4, 2009


Hello visitor, its been a while since my last post. Starting from now, I will try to post my blog in English language. This is hard for me because I neither English nor Americans. I am Indonesian. I hope you can understand if my grammar or structure maybe very confusing. Please pardon me. I still have to learn. Ok back to the topic, few days ago, I had downloaded the trial/demo version of Leica Photogrammetry Suite version 9.2 from their FTP (ftp.gi-leicageosystems.com). After exploring some new modules, I get excite with a new image fusion/pan sharpening algorithm that as far I know, not implemented in previous IMAGINE or LPS version. The algorithm was developed by Mr Manfred Ehlers from Onasbruck University Germany, and it is called Ehlers fusion. I have tried to use it to fuse quickbird 16 BIT sample data from digitalglobe (multispectral and panchromatic), and the result is awesome. Comparison between original dataset and sharpening result indicated that this algorithm can maintain spectral integrity from original multispectral dataset without losing spatial details from panchromatic dataset. Looking the result, I have idea to compare this algorithm with Mr Yun Zhang PANSHARP algorithm that now implemented in PCI Geomatica (www.pcigeomatics.com), and Gram-schmidt Spectral Sharpening that now implemented in ITT ENVI (www.ittvis.com). The reason why I chosen PANSHARP and Gram-schmidt is both algorithms are the best pan sharpening algorithms I currently know. These algorithms can preserve either spectral integrity from multispectral bands or spatial details from panchromatic band. The picture below is the comparison result that visualized on ERDAS IMAGINE GLT viewer.

images viewed in GLT viewer



From the picture, we can see that Ehlers fusion giving superior spatial details than Gram-schmidt (GS) and PANSHARP. GS and PANSHARP still leave little artifacts (spectral residuals from original multispectral dataset, especially in bright objects) on the sharpened imagery whereas in the ehlers fusion result, they re totally blended.


I have tried to assess the fusion results using correlation analysis between original multispectral bands and the sharpened bands. Below are the scatter plot results from band 1 (VIS Blue) comparison.





yun zhang


And below are the correlation analysis results.

band x (original bands)band y (sharpened bands)Linear regression equationr
ori_1YZ_10,96x + 10,370,96
ori_2YZ_20,96x + 14,920,96
ori_3YZ_30,96x + 10,90,96
ori_4YZ_40,96x + 14,160,96
ori_1ehlers_10,83x + 84,540,83
ori_2ehlers_20,88x + 41,990,89
ori_3ehlers_30,93x + 9,290,95
ori_4ehlers_40,97x + 11,440,97
ori_1GS_10,90x + 23,340,94
ori_2GS_20,91x + 32,040,94
ori_3GS_30,94x + 15,660,95
ori_4GS_40,99x + 4,570,96


Explanation : 1 = band 1 (blue); 2 = band 2 (green ); 3 = band 3 (red); 4 = band 4 (NIR)


Correlation analysis results showed that Yun Zhang PANSHARP algorithm give stable spectral preservation at all bands. Opposite with the spatial quality analysis result, Ehlers fusion can’t maintain the spectral preservation at all bands than GS or YZ PANSHARP. Although at band 4 (NIR) Ehlers give the best correlation with the original band, the correlation was dropped on the other bands. Although some algorithm give better result than the others, we can concluded that all the three algorithms give satisfied spectral preservation result and high correlation coefficient with the original multispectral bands. These algorithms are better than the more traditional algorithms like IHS fusion, CN Spectral sharpening, Brovey, and PC fusion.

Comment please.