Deconvolution algorithm based on automatic cutoff frequency selection for EPR imaging

The large line-width associated with electron paramagnetic resonance imaging (EPRI) requires effective algorithms to deconvolve the true spatial profiles of spins from the measured projection data. The commonly used Fourier transform (FT) deconvolution algorithm is easy to implement but suffers from the division-by-zero problem. At a result, a couple of parameters are used to control the deconvolution performance. However, this is inconvenient and the deconvolution results are subject to the experience of the operators. In the present work we examined FT deconvolution for EPRI, and proposed an automatic algorithm to determine the cutoff frequency by calculating the piecewise variance of the division result of the Fourier amplitude spectra. The deconvolution algorithm and the filtered back-projection image reconstruction algorithm were implemented and validated using 3D phantom and in vivo imaging data. It was clearly observed that the image resolution improved after deconvolution with the proposed algorithm. Magn Reson Med 50:444-448, 2003. (C) 2003 Wiley-Liss, Inc.