Three-dimensional reconstruction from cone-beam data using an efficient Fourier technique combined with a special interpolation filter.
We here present LINCON(FAST) which is an exact method for 3D reconstruction from cone-beam projection data. The new method is compared to the LINCON method which is known to be fast and to give good image quality. Both methods have O(N3 log N) complexity and are based on Grangeat's result which states that the derivative of the Radon transform of the object function can be obtained from cone-beam projections. One disadvantage with LINCON is that the rather computationally intensive chirp z-transform is frequently used. In LINCON(FAST), FFT and interpolation in the Fourier domain are used instead, which are less computationally demanding. The computation tools involved in LINCON(FAST) are solely FFT, 1D eight-point interpolation, multiplicative weighting and tri-linear interpolation. We estimate that LINCON(FAST) will be 2-2.5 times faster than LINCON. The interpolation filter belongs to a special class of filters developed by us. It turns out that the filter must be very carefully designed to keep a good image quality. Visual inspection of experimental results shows that the image quality is almost the same for LINCON and the new method LINCON(FAST). However, it should be remembered that LINCON(FAST) can never give better image quality than LINCON, since LINCON(FAST) is designed to approximate LINCON as well as possible.
Seger
M M
MM
Department of Electrical Engineering, Linköping University, Sweden. maria@isy.liu.se
D002138:Calibration; D005583:Fourier Analysis; D007091:Image Processing, Computer-Assisted / Q000295:instrumentation / Q000379:methods*; D008962:Models, Theoretical; D019047:Phantoms, Imaging; D011886:Radon; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D014055:Tomography, Emission-Computed
D011886:Radon