Implementing Efficient iterative 3D Deconvolution for Microscopy

Abstract: Both Gauss-Seidel Iterative 3D deconvolution and Richardson-Lucy like algorithms are used due to their stability and high quality results in high noise microscopic medical image processing. An approach to determine the difference between these two algorithms is presented in this paper. It is shown that the convergence rate and the quality of these two algorithms are influenced by the size of the point spread function (PSF). Larger PSF sizes causes faster convergence but this effect falls off for larger sizes . It is furthermore shown that the relaxation factor and the number of iterations are influencing the convergence rate of the two algorithms. It has been found that increasing relaxation factor and number of iterations improve convergence and can reduce the error of the deblurred image. It also found that overrelaxation converges faster than underrelaxation for small number of iterations. However, it can be achieved smaller final error with under-relaxation. The choice of underrelaxation factor and overrelaxation factor value are highly problem specific and different from one type of images. In addition, when it comes to 3D iterative deconvolution, the influence of boundary conditions for these two algorithms is discussed. Implementation aspects are discussed and it is concluded that cache memory is vital for achieving a fast implementation of iterative 3D deconvolution. A mix of the two algorithms have been developed and compared with the previously mentioned Gauss-Seidel and the Richardson-Lucy-like algorithms. The experiments indicate that, if the value of the relaxation parameter is optimized, then the Richardson-Lucy-like algorithm has the best performance for 3D iterative deconvolution.