Multichannel algorithm for fast 3D reconstruction.

Some recent medical imaging applications such as functional imaging (PET and SPECT) or interventional imaging (CT fluoroscopy) involve increasing amounts of data. In order to reduce the image reconstruction time, we develop a new fast 3D reconstruction algorithm based on a divide and conquer approach. The proposed multichannel algorithm performs an indirect frequential subband decomposition of the image f to be reconstructed (f = sigma fj) through the filtering of the projections Rf. The subband images fj are reconstructed on a downsampled grid without information suppression. In order to reduce the computation time, we do not backproject the null filtered projections and we downsample the number of projections according to the Shannon conditions associated with the subband image. Our algorithm is based on filtering and backprojection operators. Using the same algorithms for these basic operators, our approach is three and a half times faster than a classical FBP algorithm for a 2D image 512 x 512 and six times faster for a 3D image 32 x 512 x 512.

[1]  Pierre Grangeat,et al.  A new computation compression scheme based on a multifrequential approach , 2000, 2000 IEEE Nuclear Science Symposium. Conference Record (Cat. No.00CH37149).

[2]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[3]  Achi Brandt,et al.  A Fast and Accurate Multilevel Inversion of the Radon Transform , 1999, SIAM J. Appl. Math..

[4]  M Kono,et al.  Solitary pulmonary nodules: evaluation of blood flow patterns with dynamic CT. , 1997, Radiology.

[5]  Martin L. Brady,et al.  A Fast Discrete Approximation Algorithm for the Radon Transform , 1998, SIAM J. Comput..

[6]  Pierre Grangeat,et al.  Theoretical framework for a dynamic cone-beam reconstruction algorithm based on a dynamic particle model. , 2002, Physics in medicine and biology.

[7]  Maria Magnusson,et al.  Artefacts and remedies in direct Fourier tomographic reconstruction , 1992, IEEE Conference on Nuclear Science Symposium and Medical Imaging.

[8]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[9]  P. Noll,et al.  MPEG digital audio coding , 1997, IEEE Signal Process. Mag..

[10]  Yoram Bresler,et al.  O(N2log2N) filtered backprojection reconstruction algorithm for tomography , 2000, IEEE Trans. Image Process..

[11]  B. Daly,et al.  Real-time CT fluoroscopy: evolution of an interventional tool. , 1999, Radiology.

[12]  M. Grass,et al.  3D cone-beam CT reconstruction for circular trajectories. , 2000, Physics in medicine and biology.