Article thumbnail

Mapping Iterative Medical Imaging Algorithm on Cell Accelerator

By Meilian Xu and Parimala Thulasiraman


Algebraic reconstruction techniques require about half the number of projections as that of Fourier backprojection methods, which makes these methods safer in terms of required radiation dose. Algebraic reconstruction technique (ART) and its variant OS-SART (ordered subset simultaneous ART) are techniques that provide faster convergence with comparatively good image quality. However, the prohibitively long processing time of these techniques prevents their adoption in commercial CT machines. Parallel computing is one solution to this problem. With the advent of heterogeneous multicore architectures that exploit data parallel applications, medical imaging algorithms such as OS-SART can be studied to produce increased performance. In this paper, we map OS-SART on cell broadband engine (Cell BE). We effectively use the architectural features of Cell BE to provide an efficient mapping. The Cell BE consists of one powerPC processor element (PPE) and eight SIMD coprocessors known as synergetic processor elements (SPEs). The limited memory storage on each of the SPEs makes the mapping challenging. Therefore, we present optimization techniques to efficiently map the algorithm on the Cell BE for improved performance over CPU version. We compare the performance of our proposed algorithm on Cell BE to that of Sun Fire ×4600, a shared memory machine. The Cell BE is five times faster than AMD Opteron dual-core processor. The speedup of the algorithm on Cell BE increases with the increase in the number of SPEs. We also experiment with various parameters, such as number of subsets, number of processing elements, and number of DMA transfers between main memory and local memory, that impact the performance of the algorithm

Topics: Research Article
Publisher: Hindawi Publishing Corporation
OAI identifier:
Provided by: PubMed Central

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.

Suggested articles


  1. (1988). A .C .K a ka n dM .S l a n e y ,Principles of Computerized Tomographic Imaging,
  2. (1996). A Comparison of rotation-based methods for iterative reconstruction algorithms,”
  3. (1989). Algebraic reconstruction
  4. (1970). Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and X-ray photography,”
  5. (2007). B r e n n e ra n dE .J .H a l l ,“ C o m p u t e dt o m o g r a p h y — a n increasing source of radiation exposure,” The New England
  6. (1996). Computed tomography using algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and X-ray photography,”
  7. (1998). Fast and accurate three dimensional reconstruction from cone-beam projecton data using algebraic methods,P h . D . thesis,
  8. (1985). Fast calculation of the exact radiological path for a three-dimensional
  9. (1994). H u d s o na n dR .S .L a r k i n ,“ A c c e l e r a t e di m a g er e c o n -struction using ordered subsets of projection data,”
  10. (2008). HPC for iterative image reconstruction in CT,”
  11. (2005). Introduction to the cell multiprocessor,”
  12. (2008). On the efficiency of iterative ordered subset reconstruction algorithms for acceleration on GPUs,”
  13. (2004). Ordered-subset simultaneous algebraicreconstructiontechniques(OS-SART),”JournalofX-Ray Science and Technology,
  14. (1972). P.Gilbert,“Iterativemethodsforthethree-dimensionalrecons t r u c t i o no fa no b j e c tf r o mp r o j e c t i o n s , ”J o u r n a lo fT h e o r e t i c a l Biology,
  15. (2006). Parallel ART for image reconstruction in CT using processor arrays,”
  16. (1998). Parallel image reconstruction on MIMD computers for threedimensional cone-beam tomography,”
  17. (2005). Parallel statistical image reconstruction for cone-beam x-ray CT on a shared memory computation platform,”
  18. (2003). Performance comparison of projectorbackprojector paris for iterative tomographic reconstruction,”
  19. (2007). Programming the cell broadband engine examples and best practices,”
  20. (2008). R.Banton.,“5criticalfactorstoconsiderwhenchoosingaprocessing solution for your HPC application,”
  21. (2000). Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction
  22. (2010). Rotation based algorithm for parallelizing OS-SART for CT on homogeneous multicore architecture,” in
  23. (1984). Simultaneous algebraic reconstruction technique (SART): a superior implementation
  24. (2001). Web-based parallel ML EM reconstruction for SPECT on SMP clusters,”