Parallel algorithm based on a frequential decomposition for dynamic 3D computed tomography

We present a tomographic reconstruction algorithm based on a frequential decomposition of the data. We show that the frequential components of the attenuation function to be identified can be reconstructed from the frequential decomposition of the data. Moreover, downsampling techniques added to the identification of null components and coupled to compression techniques, speed up the reconstruction time up to six compared to the classical FBP. We identify the optimal number of frequential components. We show reconstructions from real data. A parallel implementation of our new algorithm is then proposed and evaluated on two small PC clusters.

[1]  Johan Waldén,et al.  Analysis of the direct Fourier method for computer tomography , 2000, IEEE Transactions on Medical Imaging.

[2]  Reinhard Männer,et al.  Fast Feldkamp-Reconstruction for Real-Time Reconstruction using C-Arm-Systems , 2002 .

[3]  F. Natterer,et al.  Sampling in Fan Beam Tomography , 1993, SIAM J. Appl. Math..

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

[5]  Bertil Gustafsson,et al.  On the direct Fourier method for computer tomography , 2000, IEEE Transactions on Medical Imaging.

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

[7]  Ronald R. Coifman,et al.  Fast Slant Stack: a Notion of Radon Transform for Data in a Cartesian Grid Which Is Rapidly Computible, Algebraically Exact, Geometrically Faithful and Invertible , 2001 .

[8]  Thomas Rodet Algorithmes rapides de reconstruction en tomographie par compression des calculs : application à la tomofluoroscopie 3D , 2002 .

[9]  G. Herman,et al.  Image reconstruction from linograms: implementation and evaluation. , 1988, IEEE transactions on medical imaging.

[10]  Henrik Turbell,et al.  Cone-Beam Reconstruction Using Filtered Backprojection , 2001 .

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

[12]  G. Herman,et al.  Linograms in Image Reconstruction from Projections , 1987, IEEE Transactions on Medical Imaging.

[13]  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.

[14]  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).

[15]  A L Oberg,et al.  Benefits and safety of CT fluoroscopy in interventional radiologic procedures. , 2001, Radiology.

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

[17]  Pierre Grangeat,et al.  Multichannel algorithm for fast 3D reconstruction. , 2002, Physics in medicine and biology.

[18]  W H Wright,et al.  CT-guided pericardial drainage catheter placement with subsequent pericardial sclerosis. , 1995, Journal of computer assisted tomography.

[19]  P.-E. Danielsson,et al.  Backprojection in O(N/sup 2/logN) time [X-ray CT image reconstruction] , 1997, 1997 IEEE Nuclear Science Symposium Conference Record.

[20]  L. Desbat,et al.  Efficient Sampling in 3D Tomography , 1996 .

[21]  Iain Goddard,et al.  High-speed cone-beam reconstruction: an embedded systems approach , 2002, SPIE Medical Imaging.

[22]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[23]  K. Katada,et al.  Guidance with real-time CT fluoroscopy: early clinical experience. , 1996, Radiology.

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

[25]  Thierry Ducourant,et al.  New CsI/a-Si 17" x 17" x-ray flat-panel detector provides superior detectivity and immediate direct digital output for general radiography systems , 1998, Medical Imaging.