REAL-TIME and high-quality reconstruction of cone-beam
|
|
- Juliet Willis
- 5 years ago
- Views:
Transcription
1 Real-Time 3D Cone Beam Reconstruction Dzmitry Stsepankou, Klaus Kornmesser, Jürgen Hesser, Reinhard Männer Abstract The paper presents a comparison of filtered backprojection and iterative approaches (modified SIRT and OSEM) where backprojection is performed on FPGA hardware and forward projection is realized by a graphics card. The FPGAarchitecture allows scaling parallel processing of the backprojection with a constant memory bandwidth and it performs its task about a factor of 6 faster than conventional graphics card architectures (using a Xilinx Virtex II 6000 compared to a floating-point implementation on a graphics hardware). For the forward projection a graphics card (with NVIDIA GeForce FX 5900) is used. The system reconstructs a volume of 56 3 voxels and 60 projections (51 ) within 6 seconds for each backprojection step and 1.5 s for the forward projection. As assumed filtered backprojection is the fastest method but shows most artifacts. SIRT and OSEM generate results with suppressed artifacts and require between 4 (SIRT) and 3 (OSEM) minutes for reconstruction. The limited precision of on the FPGA does not lead to any degradation in the reconstructed image quality. Index Terms Cone-beam reconstruction, FBP, iterative reconstruction, real-time reconstruction, FPGA, hardware-acceleration. I. INTRODUCTION REAL-TIME and high-quality reconstruction of cone-beam data is highly relevant for operation control and online registration in radiation therapy; and it becomes more and more a limiting factor for further progress. Additionally, the reconstruction algorithms should be able to cope with the limited data problem in order to be able to reduce the dose for the patient. In this paper we discuss an evaluation of some modern reconstruction techniques implemented on our hardware system that fulfills the aforementioned requirements. One standard approach to reconstruct volumes from x-ray images is the filtered backprojection algorithm (FBP) [1]. It essentially performs first a filtering of the projection data and then the backprojection of them towards the source. This algorithm is often used due to its good compromise between reconstruction quality and speed []. If more reconstruction time is available, there are more sensitive approaches like ART (Algebraic Reconstruction Technique) that both can reduce reconstruction artifacts and the necessary data that is required [3]. The reconstruction time for a 56 3 volume typically requires about seconds for single processor machines [4]. There are several approaches to speed up processing. For example [4] uses a PCI-board with 8 RISC like processors. Others implemented the backprojection algorithm on a special-purpose processor [5] or on FPGAs [6] [7]. Texture mapping on graphics hardware was used by Cabral [8] for FBP and Mueller D. Stsepankou (dstepank@rumms.uni-mannheim.de) is with Institute for Computational Medicine, Universities Mannheim and Heidelberg, B6, 6, D Mannheim, Germany. [9] for several iterative techniques (ART and SART), and most recently [10] [11] OSEM (Ordered Subsets Estimation Maximization) [1] using an NVidia GeForce FX graphics hardware. The main problem of backprojection is the memory bandwidth-limitation, i.e. increasing the internal processor speed does not increase the performance of the whole system. While the performance of the processor doubles every 18 months, the memory bandwidth does not cope with that performance increase and therefore there is a widening gap. All approaches so far have to cope with this limitation. In our approach we have overcome this limitation by a dedicated new reconstruction architecture that is implemented on an FPGA system. The performance increase of FPGAs over years is by a factor of 4 every 18 months - even larger than for PCs and therefore it is assumed that in the near future they will dominate high-performance applications. II. ALGORITHMS In our work we consider the transmission X-ray cone-beam CT case with the flat-panel detector. Y P oisson(b 0 e a iµ i i + r) (1) where b 0 is a blank count, µ i is the absorption coefficient, a i are the interpolation coefficients and r stands for both background and detector noise. After acquisition the negative logarithm usually applied on data. p = log (Y/b 0 ) () To evaluate proposed architecture, we have applied several different reconstruction approaches. The Filtered backprojection is the analytical reconstruction method. The FDK [1] algorithm accumulates the weighted projection data on the voxels of the 3D grid where the weighting of the projection data accounts for the divergent beam for flat-panel: D p uv = p uv (3) D + u + v Here D is a source to detector distance, u and v are the offsets from the detector center. Then the projection data are filtered one dimensionally with the ramp filter [13] along lines perpendicular to the rotation axis. And finally the accumulation is performed: µ xyz = 1 N φ 1 k=0 R p(φ, u(φ, r),v(φ, r)) (4) L(r) The functions u(φ, r) and v(φ, r) map the object coordinates r(x, y, z) into detector coordinates u, v for particular view
2 angle φ. R - source to isocenter distance, L(r) is the distance from the source to the point r(x, y, z) to be reconstructed, projected onto the central ray. All these geometric parameters are conveniently obtained from the CT projection matrices. The iterative approaches make use of projection operation: p i = N 3 1 k=0 µ k w ki, (5) where w k are the weighting coefficients, and the backprojection operation: µ k = N φ 1 i=0 p i l ik (6) Here l i is the combination of factors from equations (3) and (4). The simulated projections are used then to calculate the correction factors for the reconstructed volume, which is incrementally updated and backprojected on each iteration. The algebraic reconstruction approaches assume the following equation system: ˆPµ = p (7) Here ˆP is the projection operator, µ is the volume data vector, p is the projection data vector. The gradient descent method solves this equation by the iterative formula: µ n+1 = µ n λ Jµ n, (8) where Jµ n is the gradient of the following functional: J(µ) = 1 ˆPµ p = 1 (µ), (9) and λ is the iteration step, which could be fixed at the simplest case. Adapting λ for each iteration however leads to better convergence. In this work we have used the method of minimum errors, firstly described by Fridman [14]. Here the iteration step is chosen by the condition of minimum the following functional µ n+1 µ, where µ is the solution of equation (7). Minimizing that functional gives: n λ n = Jµ n (10) Since we have the ill-posed problem, we can incorporate some penalty term into equation. N 3 1 R(µ) = w kh V (µ h µ k ), (11) k=0 h N k where N k is the neighborhood of voxel k. Here we use the MRF Gibbs penalty [15] with the edge-preserving Huber potential function: V (µ) = { µ δ, for µ δ, µ δ/ δ, for µ δ. (1) SIRT Simultaneous Iterative Reconstruction Technique) [16] performs the update for whole volume each step. Let label the forward projection operation (5) by FP(µ) and the backprojection (6) by BP(p). Combining the above equations with SIRT we get the following iterative formula: µ n+1 = µ n λ n (BP(FP(µ n ) p 0 )+βr(µ n )), (13) where p 0 is initial projections, β - is the regularization parameter. As an example of statistical iterative approaches we have implemented the OS-EM [1] algorithm. µ n+1 = µ n BP( )/BP (1) (14) FP(µ n ) This algorithm was originally derived for emission tomography, however some researchers had successfully applied it for transmission as well [17]. Since the poisson distribution of the raw scanner data is damaged by log in the transmission case (see ()), we had implemented the NEC scaling, proposed in [18], to restore the poisson distribution of the sinogram. III. IMPLEMENTATION Most of the reconstruction approaches could be subdivided into following steps: preprocessing (i.e. filtering), backprojection, forward projection (for iterative algorithms) and postprocessing (i.e. volume update calculation in iterative approaches). The most time consuming parts are backprojection and (forward) projection. Therefore, these two steps require hardware acceleration. Since graphics hardware is best suited to the projection step we choose the backprojection for FPGA acceleration. A. Backprojection In backprojection the contribution of each projection pixel to each affected voxel has to be computed leading to an overall complexity of the operation of O(N 4 ), assuming a volume of O(N 3 ) and O(N) projections. Similarly, we require the same amount of memory accesses defining the memory bandwidth. Since the internal clock-frequency of FPGAs is much smaller than that of normal CPUs one has to introduce parallelism to cope for that disadvantage. One typical approach is pipelining, i.e. the innermost loop of the reconstruction algorithm, i.e., how to compute the contribution of a projection to one voxel, is done on a hardware pipeline, that produces one result per clock cycle. In order to further increase the reconstruction speed, several pipelines are used where each pipeline processes the contribution of one projection plane to the considered voxel. However, since each pipeline needs access to projection and voxel data a higher number of pipelines require a proportional higher memory bandwidth, i.e. more memory devices need to be accessed in parallel. The main innovation that decouples memory bandwidth from the degree of parallelization is a suited caching mechanism. We essentially assume that the source-detector-unit rotates about one fixed axis about the imaged object. Next, we assume that we reconstruct slice-wise (see Fig 1). The rotation axis is p 0
3 is reduced by a factor of M. Fig 3 represents the block-diagram of the architecture. The backprojection architecture is voxel-driven and uses bilinear interpolation in the projection data. For this system we have developed a schema that assumes the rotation of the sourcedetector-unit about a given axis. Slabs perpendicular to this axis are backprojected using several projections in parallel. In the current system 1 pipelines are implemented. This means that we process the projection data into bundles of 1, i.e. we first apply backprojection on the first bundle obtaining a volume that is improved by the next bundle etc. Fig. 1. Slice-wise backprojection parallel to the slice normal and therefore the projection of one slice onto the acquired projections covers two or more lines depending on the cone beam angle. The larger the angle, the larger the number of such lines. For an efficient caching strategy we should be able to store all lines locally on the FPGA in order to minimize the cache reloading (see Fig ). We assume a pipelined system with identical pipelines and where each pipeline computes the contribution of a projection to the same voxel of the considered volume slice. Thus, for M pipelines we update M contributions for one voxel in parallel before considering the next voxel in the slice. Fig. 3. Block diagram of FPGA design B. Projection With the appearance of the new generation of consumer graphic cards (such as Nvidia GeForce FX), it becomes possible to accelerate the high quality projection of 16bit images using 3bit floating point calculations. The projection is the common operation for such a hardware, therefore the implementation of it is straightforward. For the line integral 3D texture mapping was selected: the volume is loaded into texture memory as a 3D texture with 16bit greyscale precision. Next, the geometry transformation (i.e. CT projection matrix) is applied where all the slices are interpolated tri-linearly. Finally, the accumulation into a high precision 3bit floating point off-screen buffer generates the projection (see Fig 4). Since we have 8 pipelines in current graphics card (GeForce FX 5900), here we have 8 pixels processed at one time. Fig.. For cache efficiency, we have to be able to store all the projection lines necessary for single volume slice Since for that configuration we can process the contribution of M pipelines to N voxels in the slice within N /M steps we have to read O(N +k N M) data. For the whole volume this amounts to O(N 3 +k N M) data accesses. Since we have only processed M instead of O(N) projections, this process has to be repeated N/M times leading to O( N M MN + N M N 3 )= O(N 3 + N 4 M )=O( N 4 M ). In other words, the memory bandwidth IV. RESULTS In our approach we programmed backprojection into a Xilinx Virtex II 6000 FPGA with 1 pipelines operating in parallel [19] [0]; it is approx. a factor of 6 faster than graphics card implementations: for a 56 3 volume and 60 projections (51 ) it requires 3s for backprojection and 3s for configuring the FPGA and data transfer. The latter is due to the fact that the FPGA is not large enough to implement the logic for parallel I/O during reconstruction.
4 Fig. 6. The results for FBP, OSEM and SIRT for sparce projection data (60 projections 51, volume size 56 3 ). Fig. 4. Accumulating of the 3D texture slices into pixel buffer Fig. 7. The results of the iterative algorithms (1, - OSEM, 3 - SIRT) for real CT datasets (10 projections 51, volume size 56 3 ). For evaluating the performance and quality we used several CT volumes of size 56 3 as basis and generate virtual projections of size 51. We have used the 3D extended Shepp- Logan phantom to demonstrate the reconstruction quality and several clinical CT datasets to show how the proposed system could reconstruct the real data. Table I shows the timings for the various algorithms. FBP shows fastest result, but has the disadvantages when the small number of projections is used. Iterative approaches consumes a longer time, however they could be used with the less projection number. TABLE I TIMINGS FOR VARIOUS ALGORITHMS Iterations Overall time Backprojection Projection FBP - 6s 6s - SIRT 10 4min 6 s/iteration 1.5 s/iteration OSEM 6 3min 1 s/iteration 15 s/iteration Firstly, we have shown that the performance of the FPGAdesign scales linearly with the amount of resources that is quadrupling approximately every two years. Currently, the fastest FPGA has four times the resources of the one we used; thus, using the same design, we could reduce backprojection including reload and transfer to 0.75 s. Secondly, the forward projection is ideally suited to graphics hardware, thus with the most recent graphics card we can cut the forward projection time by a factor of. However we have observed that reading the projections from GPU memory is currently a slow operation and it could be a bottleneck for further speedup. In addition, the quality of the results is not affected by the limited precision of computing on FPGAs since we apply a problem specific scaling before performing backprojection. ACKNOWLEDGMENT This project was supported by BMBF grant 01EZ004. REFERENCES Fig. 5. The original of 3D extended Shepp-Logan phantom and the FBP reconstructions using software and hardware implementations (40 projections 51, volume size 56 3 ). V. CONCLUSIONS AND FUTURE WORK In the paper we have shown how our proposed hardware architecture could be used to accelerate different reconstruction techniques. The main advantage of this system is twofold. [1] L. A. Feldkamp, L. C. Davis, and J. W. Kress, Practical cone-beam algorithm, J. Opt. Sco. Amer, vol. 1, no. A6, pp , [] T. H. Lin, G. Wang, and P. C. Weng, A multiple cone beam reconstruction algorithm for X-Ray microtomography, Springer Series in Optical Science, vol. 67, X-Ray Microscopy III, 199. [3] R. Gordon, R. Bender, and G. T. Herman, Algebraic reconstruction techniques for three-dimensional electron microscopy and X-Ray photography, J. Theor. Biol, vol. 9, pp , [4] [5] T. Schmitt, D. Fimmel, M. Kortke, and R. Merker, Parallel processor arrays for filtered backprojection, Computer Aided Systems Theory- EUROCAST 99, pp , Springer, 000. [6] I. Goddard and M. Trepanier, High-speed cone-beam reconstruction an embedded systems approach, Proc. SPIE Medical Imaging, vol. 4681, pp , 00. [7] S. Coric, M. Leeser, E. Miller, and M. Trepanier, Parallel-beam backprojection: an FPGA implementation optimized for medical imaging, Proceedings of the 00 ACM/SIGDA tenth international symposium on Field-programmable gate arrays, pp. 17 6, 00.
5 [8] B. Cabral, N. Cam, and J. Foran, Accelerated volume rendering and tomographic reconstruction using texture mapping hardware, Proceedings of the 1994 symposium on Volume visualization, pp , [9] K. Mueller, Fast and accurate three-dimensional reconstruction from cone-beam projection data using algebraic methods, PhD dissertation. The Ohio State University, [10] K. Chidlow and T. Moeller, Rapid emission tomography reconstruction, Workshop on Volume graphics (VG03), pp. 15 6, 003. [11] F. Xu and K. Mueller, Towards a unified framework for rapid computed tomography on commodity GPUs, IEEE Medical Imaging Conference, 003. [1] H. M. Hudson and R. S. Larkin, Accelerated image reconstruction using ordered subsets of projection data, IEEE Transactions on Medical Imaging, vol. 13, no. 4, pp , Dec [13] L. T. Chang and G. T. Herman, A scientific study of filter selection for a fan-beam convolution algorithm, SIAM J. Appl. Math., vol. 39, pp , [14] V. M. Fridman, The method of minimum iterations with minimum errors for a system of linear algebraic equations with a symmetrical matrix, USSR Comput. Math. Math. Phys., no., pp , [15] N. Villain, Y. Goussard, J. Idier, and M. Allain, Three-dimensional edgepreserving image enhancement for computed tomography, IEEE Trans. Medical Imaging, vol., no. 10, pp , 003. [16] A. C. Kak and M. Slaney, Principles of Computerized Tomographic Imaging. IEEE Press, [17] J. Nuyts, B. D. Man, P. Dupont, M. Defrise, P. Suetens, and L. Mortelmans, Iterative reconstruction for helical ct: a simulation study, Phys Med Biol, no. 43, pp , [18] J. Nuyts, C. Michel, and P. Dupont, Maximum-likelihood expectationmaximization reconstruction of sinograms with arbitrary noise distribution using NEC-transformations, IEEE Trans. Medical Imaging, vol. 0, pp , 001. [19] K. Kornmesser, B. Schaedler, M. Ebert, J. Hesser, W. Schlegel, and R. Maenner, Fast feldkamp-reconstruction for real-time reconstruction using C-Arm-systems, CARS 00 Paris, France, Jun [0] D. Stsepankou, U. Mueller, K. Kornmesser, J. Hesser, and R. Maenner, FPGA-accelerated volume reconstruction from X-Ray, World Conference on Med. Phys. and Biomed. Eng, 003.
Accelerated C-arm Reconstruction by Out-of-Projection Prediction
Accelerated C-arm Reconstruction by Out-of-Projection Prediction Hannes G. Hofmann, Benjamin Keck, Joachim Hornegger Pattern Recognition Lab, University Erlangen-Nuremberg hannes.hofmann@informatik.uni-erlangen.de
More informationDEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE
Rajesh et al. : Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh V Acharya, Umesh Kumar, Gursharan
More informationHigh Resolution Iterative CT Reconstruction using Graphics Hardware
High Resolution Iterative CT Reconstruction using Graphics Hardware Benjamin Keck, Hannes G. Hofmann, Holger Scherl, Markus Kowarschik, and Joachim Hornegger Abstract Particular applications of computed
More informationOptimization of Cone Beam CT Reconstruction Algorithm Based on CUDA
Sensors & Transducers 2013 by IFSA http://www.sensorsportal.com Optimization of Cone Beam CT Reconstruction Algorithm Based on CUDA 1 Wang LI-Fang, 2 Zhang Shu-Hai 1 School of Electronics and Computer
More informationAn approximate cone beam reconstruction algorithm for gantry-tilted CT
An approximate cone beam reconstruction algorithm for gantry-tilted CT Ming Yan a, Cishen Zhang ab, Hongzhu Liang a a School of Electrical & Electronic Engineering, Nanyang Technological University, Singapore;
More informationA Projection Access Scheme for Iterative Reconstruction Based on the Golden Section
A Projection Access Scheme for Iterative Reconstruction Based on the Golden Section Thomas Köhler Philips Research Laboratories Roentgenstrasse - Hamburg Germany Abstract A new access scheme for projections
More informationFeldkamp-type image reconstruction from equiangular data
Journal of X-Ray Science and Technology 9 (2001) 113 120 113 IOS Press Feldkamp-type image reconstruction from equiangular data Ben Wang a, Hong Liu b, Shiying Zhao c and Ge Wang d a Department of Elec.
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationOn the Efficiency of Iterative Ordered Subset Reconstruction Algorithms for Acceleration on GPUs
On the Efficiency of Iterative Ordered Subset Reconstruction Algorithms for Acceleration on GPUs Fang Xu 1, Klaus Mueller 1, Mel Jones 2, Bettina Keszthelyi 2, John Sedat 2, David Agard 2 1 Center for
More informationAdaptive region of interest method for analytical micro-ct reconstruction
Journal of X-Ray Science and Technology 19 (2011) 23 33 23 DOI 10.3233/XST-2010-0274 IOS Press Adaptive region of interest method for analytical micro-ct reconstruction Wanneng Yang, Xiaochun Xu, Kun Bi,
More informationCombined algorithmic and GPU acceleration for ultra-fast circular conebeam backprojection
Combined algorithmic and GPU acceleration for ultra-fast circular conebeam backprojection Jeffrey Brokish a, Paul Sack a, Yoram Bresler ab a InstaRecon, Inc., Champaign, IL USA 61820 b Dept. of Electrical
More informationClinical Evaluation of GPU-Based Cone Beam Computed Tomography.
Clinical Evaluation of GPU-Based Cone Beam Computed Tomography. Peter B. Noël 1,2, Alan M. Walczak 2, Kenneth R. Hoffmann 1,2, Jinhui Xu 1, Jason J. Corso 1, and Sebastian Schafer 2 1 Department of Computer
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/8289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 2017-0-13
More informationARTICLE IN PRESS. computer methods and programs in biomedicine xxx (2009) xxx xxx. journal homepage:
computer methods and programs in biomedicine xxx (2009) xxx xxx journal homepage: www.intl.elsevierhealth.com/journals/cmpb GPU-based cone beam computed tomography Peter B. Noël a,b,, Alan M. Walczak a,c,
More informationMedical Image Reconstruction Term II 2012 Topic 6: Tomography
Medical Image Reconstruction Term II 2012 Topic 6: Tomography Professor Yasser Mostafa Kadah Tomography The Greek word tomos means a section, a slice, or a cut. Tomography is the process of imaging a cross
More informationTranslational Computed Tomography: A New Data Acquisition Scheme
2nd International Symposium on NDT in Aerospace 2010 - We.1.A.3 Translational Computed Tomography: A New Data Acquisition Scheme Theobald FUCHS 1, Tobias SCHÖN 2, Randolf HANKE 3 1 Fraunhofer Development
More informationImplementation of a backprojection algorithm on CELL
Implementation of a backprojection algorithm on CELL Mario Koerner March 17, 2006 1 Introduction X-ray imaging is one of the most important imaging technologies in medical applications. It allows to look
More informationScaling Calibration in the ATRACT Algorithm
Scaling Calibration in the ATRACT Algorithm Yan Xia 1, Andreas Maier 1, Frank Dennerlein 2, Hannes G. Hofmann 1, Joachim Hornegger 1,3 1 Pattern Recognition Lab (LME), Friedrich-Alexander-University Erlangen-Nuremberg,
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner March 17, 2006 1 2D Fan Beam Reconstruction Two-dimensional objects can be reconstructed from projections that were acquired using parallel
More informationMulti-slice CT Image Reconstruction Jiang Hsieh, Ph.D.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced
More informationA numerical simulator in VC++ on PC for iterative image reconstruction
Journal of X-Ray Science and Technology 11 (2003) 61 70 61 IOS Press A numerical simulator in VC++ on PC for iterative image reconstruction Xiang Li a,, Ming Jiang c and Ge Wang a,b a Department of Biomedical
More informationOn the Use of Graphics Hardware to Accelerate Algebraic Reconstruction Methods
See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/2621187 On the Use of Graphics Hardware to Accelerate Algebraic Reconstruction Methods Article
More informationCUDA and OpenCL Implementations of 3D CT Reconstruction for Biomedical Imaging
CUDA and OpenCL Implementations of 3D CT Reconstruction for Biomedical Imaging Saoni Mukherjee, Nicholas Moore, James Brock and Miriam Leeser September 12, 2012 1 Outline Introduction to CT Scan, 3D reconstruction
More informationRGBA Packing for Fast Cone Beam Reconstruction on the GPU
RGBA Packing for Fast Cone Beam Reconstruction on the GPU Fumihiko Ino a, Seiji Yoshida a, and Kenichi Hagihara a a Graduate School of Information Science and Technology, Osaka University, 1-3 Machikaneyama,
More informationPortability of TV-Regularized Reconstruction Parameters to Varying Data Sets
Portability of TV-Regularized Reconstruction Parameters to Varying Data Sets Mario Amrehn 1, Andreas Maier 1,2, Frank Dennerlein 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
More informationhigh performance medical reconstruction using stream programming paradigms
high performance medical reconstruction using stream programming paradigms This Paper describes the implementation and results of CT reconstruction using Filtered Back Projection on various stream programming
More informationTomographic Reconstruction
Tomographic Reconstruction 3D Image Processing Torsten Möller Reading Gonzales + Woods, Chapter 5.11 2 Overview Physics History Reconstruction basic idea Radon transform Fourier-Slice theorem (Parallel-beam)
More informationHigh-performance tomographic reconstruction using graphics processing units
18 th World IMACS / MODSIM Congress, Cairns, Australia 13-17 July 29 http://mssanz.org.au/modsim9 High-performance tomographic reconstruction using graphics processing units Ya.I. esterets and T.E. Gureyev
More informationA Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT
A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT Daniel Schlifske ab and Henry Medeiros a a Marquette University, 1250 W Wisconsin Ave, Milwaukee,
More informationGPU implementation for rapid iterative image reconstruction algorithm
GPU implementation for rapid iterative image reconstruction algorithm and its applications in nuclear medicine Jakub Pietrzak Krzysztof Kacperski Department of Medical Physics, Maria Skłodowska-Curie Memorial
More informationSpiral ASSR Std p = 1.0. Spiral EPBP Std. 256 slices (0/300) Kachelrieß et al., Med. Phys. 31(6): , 2004
Spiral ASSR Std p = 1.0 Spiral EPBP Std p = 1.0 Kachelrieß et al., Med. Phys. 31(6): 1623-1641, 2004 256 slices (0/300) Advantages of Cone-Beam Spiral CT Image quality nearly independent of pitch Increase
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/887/48289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 207-04-3
More informationComparison of Probing Error in Dimensional Measurement by Means of 3D Computed Tomography with Circular and Helical Sampling
nd International Symposium on NDT in Aerospace - We..A. Comparison of Probing Error in Dimensional Measurement by Means of D Computed Tomography with Circular and Helical Sampling Jochen HILLER, Stefan
More informationAccelerating Cone Beam Reconstruction Using the CUDA-enabled GPU
Accelerating Cone Beam Reconstruction Using the CUDA-enabled GPU Yusuke Okitsu, Fumihiko Ino, and Kenichi Hagihara Graduate School of Information Science and Technology, Osaka University 1-3 Machikaneyama,
More informationWorkshop on Quantitative SPECT and PET Brain Studies January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET
Workshop on Quantitative SPECT and PET Brain Studies 14-16 January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET Físico João Alfredo Borges, Me. Corrections in SPECT and PET SPECT and
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction. Mario Koerner
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner Moscow-Bavarian Joint Advanced Student School 2006 March 19 2006 to March 29 2006 Overview 2D Fan Beam Reconstruction Shortscan Reconstruction
More informationRapid Emission Tomography Reconstruction
Volume Graphics (2003) I. Fujishiro, K. Mueller, A. Kaufman (Editors) Rapid Emission Tomography Reconstruction Ken Chidlow Torsten Möller Department of Computer Science, Simon Fraser University, Burnaby,
More informationAxial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction
Axial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction Jeffrey A. Fessler and Donghwan Kim EECS Department University of Michigan Fully 3D Image Reconstruction Conference July
More informationEfficient 3D Crease Point Extraction from 2D Crease Pixels of Sinogram
Efficient 3D Crease Point Extraction from 2D Crease Pixels of Sinogram Ryo Jinnouchi 1, Yutaka Ohtake 1, Hiromasa Suzuki 1, Yukie Nagai 1 1 The University of Tokyo, Japan, e-mail: {jinnouchi, yu-ohtake,
More informationReconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm
Reconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm ABDESSALEM BENAMMAR, AICHA ALLAG and REDOUANE DRAI Research Center in Industrial Technologies (CRTI), P.O.Box 64, Cheraga
More informationApproximating Algebraic Tomography Methods by Filtered Backprojection: A Local Filter Approach
Fundamenta Informaticae 135 (2014) 1 19 1 DOI 10.3233/FI-2014-1109 IOS Press Approximating Algebraic Tomography Methods by Filtered Backprojection: A Local Filter Approach Linda Plantagie Centrum Wiskunde
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationCIVA Computed Tomography Modeling
CIVA Computed Tomography Modeling R. FERNANDEZ, EXTENDE, France S. LEGOUPIL, M. COSTIN, D. TISSEUR, A. LEVEQUE, CEA-LIST, France page 1 Summary Context From CIVA RT to CIVA CT Reconstruction Methods Applications
More informationImprovement of Efficiency and Flexibility in Multi-slice Helical CT
J. Shanghai Jiaotong Univ. (Sci.), 2008, 13(4): 408 412 DOI: 10.1007/s12204-008-0408-x Improvement of Efficiency and Flexibility in Multi-slice Helical CT SUN Wen-wu 1 ( ), CHEN Si-ping 2 ( ), ZHUANG Tian-ge
More informationSPECT reconstruction
Regional Training Workshop Advanced Image Processing of SPECT Studies Tygerberg Hospital, 19-23 April 2004 SPECT reconstruction Martin Šámal Charles University Prague, Czech Republic samal@cesnet.cz Tomography
More informationStatistical iterative reconstruction using fast optimization transfer algorithm with successively increasing factor in Digital Breast Tomosynthesis
Statistical iterative reconstruction using fast optimization transfer algorithm with successively increasing factor in Digital Breast omosynthesis Shiyu Xu a and Zhenxi Zhang b and Ying Chen a,* a Department
More informationIntroduction to Medical Imaging. Cone-Beam CT. Klaus Mueller. Computer Science Department Stony Brook University
Introduction to Medical Imaging Cone-Beam CT Klaus Mueller Computer Science Department Stony Brook University Introduction Available cone-beam reconstruction methods: exact approximate algebraic Our discussion:
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationEstimating 3D Respiratory Motion from Orbiting Views
Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng, Jeffrey A. Fessler, James M. Balter The University of Michigan Oct. 2005 Funding provided by NIH Grant P01 CA59827 Motivation Free-breathing
More informationUnmatched Projector/Backprojector Pairs in an Iterative Reconstruction Algorithm
548 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 19, NO. 5, MAY 2000 Unmatched Projector/Backprojector Pairs in an Iterative Reconstruction Algorithm Gengsheng L. Zeng*, Member, IEEE, and Grant T. Gullberg,
More informationReconstruction from Projections
Reconstruction from Projections M.C. Villa Uriol Computational Imaging Lab email: cruz.villa@upf.edu web: http://www.cilab.upf.edu Based on SPECT reconstruction Martin Šámal Charles University Prague,
More informationComputed Tomography. Principles, Design, Artifacts, and Recent Advances. Jiang Hsieh THIRD EDITION. SPIE PRESS Bellingham, Washington USA
Computed Tomography Principles, Design, Artifacts, and Recent Advances THIRD EDITION Jiang Hsieh SPIE PRESS Bellingham, Washington USA Table of Contents Preface Nomenclature and Abbreviations xi xv 1 Introduction
More informationA Weighted Least Squares PET Image Reconstruction Method Using Iterative Coordinate Descent Algorithms
A Weighted Least Squares PET Image Reconstruction Method Using Iterative Coordinate Descent Algorithms Hongqing Zhu, Huazhong Shu, Jian Zhou and Limin Luo Department of Biological Science and Medical Engineering,
More informationCentral Slice Theorem
Central Slice Theorem Incident X-rays y f(x,y) R x r x Detected p(, x ) The thick line is described by xcos +ysin =R Properties of Fourier Transform F [ f ( x a)] F [ f ( x)] e j 2 a Spatial Domain Spatial
More informationCIVA CT, an advanced simulation platform for NDT
More Info at Open Access Database www.ndt.net/?id=18774 CIVA CT, an advanced simulation platform for NDT Marius Costin 1, David Tisseur 1, Caroline Vienne 1, Ronan Guillamet 1, Hussein Banjak 1, Navnina
More informationARTICLE IN PRESS. computer methods and programs in biomedicine xxx (2009) xxx xxx. journal homepage:
computer methods and programs in biomedicine xxx (2009) xxx xxx journal homepage: www.intl.elsevierhealth.com/journals/cmpb On the efficiency of iterative ordered subset reconstruction algorithms for acceleration
More informationNon-Stationary CT Image Noise Spectrum Analysis
Non-Stationary CT Image Noise Spectrum Analysis Michael Balda, Björn J. Heismann,, Joachim Hornegger Pattern Recognition Lab, Friedrich-Alexander-Universität Erlangen Siemens Healthcare, Erlangen michael.balda@informatik.uni-erlangen.de
More informationTHE FAN-BEAM scan for rapid data acquisition has
190 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 26, NO. 2, FEBRUARY 2007 Hilbert Transform Based FBP Algorithm for Fan-Beam CT Full Partial Scans Jiangsheng You*, Member, IEEE, Gengsheng L. Zeng, Senior
More informationEvaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT
Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig
More informationDiscrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset
Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset Andreas Maier 1, Patrick Kugler 2, Günter Lauritsch 2, Joachim Hornegger 1 1 Pattern Recognition Lab and SAOT Erlangen,
More informationGeneralized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction
Generalized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction Klaus Erhard a, Michael Grass a, Sebastian Hitziger b, Armin Iske b and Tim Nielsen a a Philips Research Europe Hamburg,
More informationTop-level Design and Pilot Analysis of Low-end CT Scanners Based on Linear Scanning for Developing Countries
Top-level Design and Pilot Analysis of Low-end CT Scanners Based on Linear Scanning for Developing Countries Fenglin Liu 1,3, Hengyong Yu 2, Wenxiang Cong 3, Ge Wang 3,* 5 1 1 Engineering Research Center
More informationImage registration for motion estimation in cardiac CT
Image registration for motion estimation in cardiac CT Bibo Shi a, Gene Katsevich b, Be-Shan Chiang c, Alexander Katsevich d, and Alexander Zamyatin c a School of Elec. Engi. and Comp. Sci, Ohio University,
More informationGPU-based cone-beam CT reconstruction with extended FOV
GPU-based cone-beam CT reconstruction with extended FOV Tamás Huszár 1, Gábor Jakab 12, and Attila Rácz 1 1 Mediso Medical Equipment Developing and Service Ltd. Budapest, 1022 Hungary http://www.mediso.hu,
More informationConvolution-Based Truncation Correction for C-Arm CT using Scattered Radiation
Convolution-Based Truncation Correction for C-Arm CT using Scattered Radiation Bastian Bier 1, Chris Schwemmer 1,2, Andreas Maier 1,3, Hannes G. Hofmann 1, Yan Xia 1, Joachim Hornegger 1,2, Tobias Struffert
More informationHardware/Software 2D-3D backprojection on a SoPC Platform
Hardware/Software 2D-3D backprojection on a SoPC Platform Nicolas GAC, Stéphane MANCINI and Michel DESVIGNES Laboratoire des Images et des Signaux 46 av Felix Viallet 38100 Grenoble, France {nicolas.gac,
More informationGPU-based Fast Cone Beam CT Reconstruction from Undersampled and Noisy Projection Data via Total Variation
GPU-based Fast Cone Beam CT Reconstruction from Undersampled and Noisy Projection Data via Total Variation 5 10 15 20 25 30 35 Xun Jia Department of Radiation Oncology, University of California San Diego,
More informationOptimization of GPU-Accelerated Iterative CT Reconstruction Algorithm for Clinical Use
Washington University in St. Louis Washington University Open Scholarship Engineering and Applied Science Theses & Dissertations Engineering and Applied Science Spring 5-18-2018 Optimization of GPU-Accelerated
More informationSINGLE-PHOTON emission computed tomography
1458 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 59, NO. 4, AUGUST 2012 SPECT Imaging With Resolution Recovery Andrei V. Bronnikov Abstract Single-photon emission computed tomography (SPECT) is a method
More informationForward-Projection Architecture for Fast Iterative Image Reconstruction in X-Ray CT
5508 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 60, NO. 10, OCTOBER 2012 Forward-Projection Architecture for Fast Iterative Image Reconstruction in X-Ray CT Jung Kuk Kim, Student Member, IEEE, Jeffrey
More informationGengsheng Lawrence Zeng. Medical Image Reconstruction. A Conceptual Tutorial
Gengsheng Lawrence Zeng Medical Image Reconstruction A Conceptual Tutorial Gengsheng Lawrence Zeng Medical Image Reconstruction A Conceptual Tutorial With 163 Figures Author Prof. Dr. Gengsheng Lawrence
More informationAccelerated quantitative multi-material beam hardening correction(bhc) in cone-beam CT
Accelerated quantitative multi-material beam hardening correction(bhc) in cone-beam CT Award: Poster No.: C-2161 Certificate of Merit Congress: ECR 2012 Type: Authors: Scientific Exhibit Q. Yang, M. Elter,
More informationSplitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction
Splitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction Hung Nien and Jeffrey A. Fessler Department of Electrical Engineering and Computer Science University of Michigan, Ann
More informationEnhancement Image Quality of CT Using Single Slice Spiral Technique
Enhancement Image Quality of CT Using Single Slice Spiral Technique Doaa. N. Al Sheack 1 and Dr.Mohammed H. Ali Al Hayani 2 1 2 Electronic and Communications Engineering Department College of Engineering,
More informationAn Efficient Technique For Multi-Phase Model Based Iterative Reconstruction
1 An Efficient Technique For Multi-Phase Model Based Iterative Reconstruction Shiyu Xu, Debashish Pal and Jean-Baptiste Thibault Abstract Multi-phase scan is a fundamental CT acquisition technology used
More informationTomography at all Scales. Uccle, 7 April 2014
Tomography at all Scales Uccle, 7 April 2014 Outline The Vision Lab ASTRA: All Scale Tomographic Reconstruction Antwerp Tomography Discrete Tomography In situ CT Superresolution Dynamic imaging The ASTRA
More informationMaterial for Chapter 6: Basic Principles of Tomography M I A Integral Equations in Visual Computing Material
Material for Chapter : Integral Equations in Visual Computing Material Basic Principles of Tomography c 00 Bernhard Burgeth 0 Source: Images Figure : Radon Transform: ttenuation http://en.wikimedia.org/wiki/image:radon_transform.png
More informationAdapted acquisition trajectory and iterative reconstruction for few-views CT inspection
Adapted acquisition trajectory and iterative reconstruction for few-views CT inspection Caroline Vienne 1, Marius Costin 1 More info about this article: http://www.ndt.net/?id=21917 1 CEA, LIST, Département
More informationAdaptive algebraic reconstruction technique
Adaptive algebraic reconstruction technique Wenkai Lua) Department of Automation, Key State Lab of Intelligent Technology and System, Tsinghua University, Beijing 10084, People s Republic of China Fang-Fang
More informationComparison of different iterative reconstruction algorithms for X-ray volumetric inspection
Comparison of different iterative reconstruction algorithms for X-ray volumetric inspection More info about this article: http://www.ndt.net/?id=22973 Georgios Liaptsis 1,2, Alan L. Clarke 1 and Perumal
More informationProjection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT
Projection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT Benedikt Lorch 1, Martin Berger 1,2, Joachim Hornegger 1,2, Andreas Maier 1,2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
More informationRapid CT reconstruction on GPU-enabled HPC clusters
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 2011 http://mssanz.org.au/modsim2011 Rapid CT reconstruction on GPU-enabled HPC clusters D. Thompson a, Ya. I.
More informationAcknowledgments. Nesterov s Method for Accelerated Penalized-Likelihood Statistical Reconstruction for C-arm Cone-Beam CT.
June 5, Nesterov s Method for Accelerated Penalized-Likelihood Statistical Reconstruction for C-arm Cone-Beam CT Adam S. Wang, J. Webster Stayman, Yoshito Otake, Gerhard Kleinszig, Sebastian Vogt, Jeffrey
More informationAN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT
1102 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 40, NO. 4, AUGUST 1993 AN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT Grant. T. Gullberg and Gengsheng L. Zeng, Department of Radiology,
More informationTomographic Algorithm for Industrial Plasmas
Tomographic Algorithm for Industrial Plasmas More info about this article: http://www.ndt.net/?id=22342 1 Sudhir K. Chaudhary, 1 Kavita Rathore, 2 Sudeep Bhattacharjee, 1 Prabhat Munshi 1 Nuclear Engineering
More informationAcknowledgments and financial disclosure
AAPM 2012 Annual Meeting Digital breast tomosynthesis: basic understanding of physics principles James T. Dobbins III, Ph.D., FAAPM Director, Medical Physics Graduate Program Ravin Advanced Imaging Laboratories
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 1, 2017 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More informationX-TRACT: software for simulation and reconstruction of X-ray phase-contrast CT
X-TRACT: software for simulation and reconstruction of X-ray phase-contrast CT T.E.Gureyev, Ya.I.Nesterets, S.C.Mayo, A.W.Stevenson, D.M.Paganin, G.R.Myers and S.W.Wilkins CSIRO Materials Science and Engineering
More informationUSING cone-beam geometry with pinhole collimation,
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 56, NO. 3, JUNE 2009 687 A Backprojection-Based Parameter Estimation Technique for Skew-Slit Collimation Jacob A. Piatt, Student Member, IEEE, and Gengsheng L.
More informationIterative and analytical reconstruction algorithms for varying-focal-length cone-beam
Home Search Collections Journals About Contact us My IOPscience Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam projections This content has been downloaded from IOPscience.
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
More informationCover Page. The handle holds various files of this Leiden University dissertation.
Cover Page The handle http://hdl.handle.net/1887/39638 holds various files of this Leiden University dissertation. Author: Pelt D.M. Title: Filter-based reconstruction methods for tomography Issue Date:
More informationTOMOGRAPHIC reconstruction problems are found in
4750 IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 23, NO. 11, NOVEMBER 2014 Improving Filtered Backprojection Reconstruction by Data-Dependent Filtering Abstract Filtered backprojection, one of the most
More informationParallel iterative cone beam CT image reconstruction on a PC cluster
Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 1 Journal of X-Ray Science and Technology 13 (2005) 1 10 1 IOS Press Parallel iterative cone beam CT image reconstruction on a PC cluster
More informationIN transmission computed tomography (CT), standard scan
1 Distributed Computation of Linear Inverse Problems with Application to Computed Tomography Reconstruction Yushan Gao, Thomas Blumensath, arxiv:179.953v1 [cs.dc] 4 Sep 217 Abstract The inversion of linear
More informationA practical software platform of integrated reconstruction algorithms and forward algorithms for 2D industrial CT
Journal of X-Ray Science and Technology 13 (2005) 9 21 9 IOS Press A practical software platform of integrated reconstruction algorithms and forward algorithms for 2D industrial CT Hui Li a, Zhaotian Zhang
More informationMéthodes d imagerie pour les écoulements et le CND
Méthodes d imagerie pour les écoulements et le CND Journée scientifique FED3G CEA LIST/Lab Imagerie Tomographie et Traitement Samuel Legoupil 15 juin 2012 2D/3D imaging tomography Example Petrochemical
More informationColumn-Action Methods in Image Reconstruction
Column-Action Methods in Image Reconstruction Per Christian Hansen joint work with Tommy Elfving Touraj Nikazad Overview of Talk Part 1: the classical row-action method = ART The advantage of algebraic
More informationCIRCULAR scanning trajectory has been widely used in
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 25, NO. 7, JULY 2006 869 Region of Interest Reconstruction From Truncated Data in Circular Cone-Beam CT Lifeng Yu, Yu Zou, Emil Y. Sidky, Charles A. Pelizzari,
More information