next up previous contents
Next: Deconvolution Up: The Wavelet Transform Previous: Hierarchical adaptive filtering

   
Comparison using a multiresolution quality criterion

It is sometimes useful, as in image restoration where we want to evaluate the quality of the restoration, to compare images with an objective criterion. Very few quantitative parameters can be extracted for that. The correlation between the original image I(i,j) and the restored one $\tilde{I}(i,j)$ gives a classical criterion. The correlation coefficient is:
$\displaystyle C_{or} = \frac{ \sum_{i = 1}^{N}\sum_{j = 1}^{N} I(i,j)
\tilde{I}...
...N}\sum_{j =
1}^{N} I^2(i,j) \sum_{i = 1}^{N}\sum_{j = 1}^{N}
\tilde{I}^2(i,j)}}$     (14.95)

The correlation is 1 if the images are identical, and less if some differences exist. Another way to compare two pictures is to determine the mean-square error:
$\displaystyle E_{ms}^2 = \frac{1}{N^2} \sum_{i = 1}^{N}\sum_{j = 1}^{N}(I(i,j)-
\tilde{I}(i,j))^2$     (14.96)

Ems2 can be normalized by:
$\displaystyle E_{nms}^2 = \frac{\sum_{i = 1}^{N}\sum_{j = 1}^{N}(I(i,j)- \tilde{I}(i,j))^2}
{\sum_{i = 1}^{N}\sum_{j = 1}^{N}I^2(i,j)}$     (14.97)

The Signal-to-Noise Ratio (SNR) corresponding to the above error is:
$\displaystyle SNR_{dB} = 10 \log_{10} \frac{1}{E_{nms}^2} \mbox{ dB}$     (14.98)

These criteria are not sufficient, they give no information on the resulting resolution. A complete criterion must take into account the resolution. For each dyadic scale, we can compute the correlation coefficient and the quadratic error between the wavelet transforms of the original and the restored images. Hence, we can compare, the quality of the restoration for each resolution.

Figures 14.21 and 14.22 show the comparison of three images with a reference image. Data20 is a simulated noisy image, median and wave are the output images after respectively applying a median filter, and a thresholding in the wavelet space. These curves show that the thresholding in the wavelet space is better than the median at all the scales.


  
Figure 14.21: Correlation.
\begin{figure}
\centerline{
\hbox{
\psfig{figure=correl.ps,height=7.5cm,width=13.5cm,angle=270}
}}
\end{figure}


  
Figure 14.22: Signal to noise ratio.
\begin{figure}
\centerline{
\hbox{
\psfig{figure=snr.ps,height=7.5cm,width=13.5cm,angle=270}
}}
\end{figure}


next up previous contents
Next: Deconvolution Up: The Wavelet Transform Previous: Hierarchical adaptive filtering
http://www.eso.org/midas/midas-support.html
1999-06-15