skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Test of 3D CT reconstructions by EM + TV algorithm from undersampled data

Abstract

Computerized tomography (CT) plays an important role in medical imaging for diagnosis and therapy. However, CT imaging is connected with ionization radiation exposure of patients. Therefore, the dose reduction is an essential issue in CT. In 2011, the Expectation Maximization and Total Variation Based Model for CT Reconstruction (EM+TV) was proposed. This method can reconstruct a better image using less CT projections in comparison with the usual filtered back projection (FBP) technique. Thus, it could significantly reduce the overall dose of radiation in CT. This work reports the results of an independent numerical simulation for cone beam CT geometry with alternative virtual phantoms. As in the original report, the 3D CT images of 128 Multiplication-Sign 128 Multiplication-Sign 128 virtual phantoms were reconstructed. It was not possible to implement phantoms with lager dimensions because of the slowness of code execution even by the CORE i7 CPU.

Authors:
; ;  [1];  [2];  [3];  [4]
  1. Universidade Tecnologica Federal do Parana - UTFPR/FB, 85601-970, Caixa Postal 135, Francisco Beltrao - PR (Brazil)
  2. Universidade Tecnologica Federal do Parana-UTFPR/FB,85601-970,Caixa Postal 135,Francisco Beltrao-PR (Brazil) and Faculdades Pequeno Principe-FPP, Av. Iguacu, 333, Rebou (Brazil)
  3. Universidade Federal de Santa Catarina - UFSC/ARA, 88900-000, Rua Pedro Joao Pereira, 150, Ararangua - SC (Brazil)
  4. Universidade Federal do Parana - UFPR Litoral, 80230-901, Rua Jaguaraiva 512, Caioba, Matinhos - PR (Brazil)
Publication Date:
OSTI Identifier:
22111970
Resource Type:
Journal Article
Resource Relation:
Journal Name: AIP Conference Proceedings; Journal Volume: 1529; Journal Issue: 1; Conference: 35. Brazilian workshop on nuclear physics, Sao Paulo (Brazil), 2-6 Sep 2012; Other Information: (c) 2013 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
62 RADIOLOGY AND NUCLEAR MEDICINE; 61 RADIATION PROTECTION AND DOSIMETRY; ALGORITHMS; CAT SCANNING; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; DIAGNOSIS; GAMMA DOSIMETRY; IMAGE PROCESSING; IMAGES; NUMERICAL ANALYSIS; PATIENTS; PHANTOMS; RADIATION DOSES; RADIOTHERAPY; VARIATIONS

Citation Formats

Evseev, Ivan, Ahmann, Francielle, Silva, Hamilton P. da, Schelin, Hugo R., Yevseyeva, Olga, and Klock, Margio C. L. Test of 3D CT reconstructions by EM + TV algorithm from undersampled data. United States: N. p., 2013. Web. doi:10.1063/1.4804089.
Evseev, Ivan, Ahmann, Francielle, Silva, Hamilton P. da, Schelin, Hugo R., Yevseyeva, Olga, & Klock, Margio C. L. Test of 3D CT reconstructions by EM + TV algorithm from undersampled data. United States. doi:10.1063/1.4804089.
Evseev, Ivan, Ahmann, Francielle, Silva, Hamilton P. da, Schelin, Hugo R., Yevseyeva, Olga, and Klock, Margio C. L. Mon . "Test of 3D CT reconstructions by EM + TV algorithm from undersampled data". United States. doi:10.1063/1.4804089.
@article{osti_22111970,
title = {Test of 3D CT reconstructions by EM + TV algorithm from undersampled data},
author = {Evseev, Ivan and Ahmann, Francielle and Silva, Hamilton P. da and Schelin, Hugo R. and Yevseyeva, Olga and Klock, Margio C. L.},
abstractNote = {Computerized tomography (CT) plays an important role in medical imaging for diagnosis and therapy. However, CT imaging is connected with ionization radiation exposure of patients. Therefore, the dose reduction is an essential issue in CT. In 2011, the Expectation Maximization and Total Variation Based Model for CT Reconstruction (EM+TV) was proposed. This method can reconstruct a better image using less CT projections in comparison with the usual filtered back projection (FBP) technique. Thus, it could significantly reduce the overall dose of radiation in CT. This work reports the results of an independent numerical simulation for cone beam CT geometry with alternative virtual phantoms. As in the original report, the 3D CT images of 128 Multiplication-Sign 128 Multiplication-Sign 128 virtual phantoms were reconstructed. It was not possible to implement phantoms with lager dimensions because of the slowness of code execution even by the CORE i7 CPU.},
doi = {10.1063/1.4804089},
journal = {AIP Conference Proceedings},
number = 1,
volume = 1529,
place = {United States},
year = {Mon May 06 00:00:00 EDT 2013},
month = {Mon May 06 00:00:00 EDT 2013}
}
  • Purpose: Cone-beam CT (CBCT) plays an important role in image guided radiation therapy (IGRT). However, the large radiation dose from serial CBCT scans in most IGRT procedures raises a clinical concern, especially for pediatric patients who are essentially excluded from receiving IGRT for this reason. The goal of this work is to develop a fast GPU-based algorithm to reconstruct CBCT from undersampled and noisy projection data so as to lower the imaging dose. Methods: The CBCT is reconstructed by minimizing an energy functional consisting of a data fidelity term and a total variation regularization term. The authors developed a GPU-friendlymore » version of the forward-backward splitting algorithm to solve this model. A multigrid technique is also employed. Results: It is found that 20-40 x-ray projections are sufficient to reconstruct images with satisfactory quality for IGRT. The reconstruction time ranges from 77 to 130 s on an NVIDIA Tesla C1060 (NVIDIA, Santa Clara, CA) GPU card, depending on the number of projections used, which is estimated about 100 times faster than similar iterative reconstruction approaches. Moreover, phantom studies indicate that the algorithm enables the CBCT to be reconstructed under a scanning protocol with as low as 0.1 mA s/projection. Comparing with currently widely used full-fan head and neck scanning protocol of {approx}360 projections with 0.4 mA s/projection, it is estimated that an overall 36-72 times dose reduction has been achieved in our fast CBCT reconstruction algorithm. Conclusions: This work indicates that the developed GPU-based CBCT reconstruction algorithm is capable of lowering imaging dose considerably. The high computation efficiency in this algorithm makes the iterative CBCT reconstruction approach applicable in real clinical environments.« less
  • The authors have implemented the EM reconstruction algorithm for volume acquisition from current generation retracted-septa PET scanners. Although the software was designed for a GE Advance scanner, it is easily adaptable to other 3D scanners. The reconstruction software was written for an Intel iPSC/860 parallel computer with 128 compute nodes. Running on 32 processors, the algorithm requires approximately 55 minutes per iteration to reconstruct a 128 x 128 x 35 image. No projection data compression schemes or other approximations were used in the implementation. Extensive use of EM system matrix (C{sub ij}) symmetries (including the 8-fold in-plane symmetries, 2-fold axialmore » symmetries, and axial parallel line redundancies) reduces the storage cost by a factor of 188. The parallel algorithm operates on distributed projection data which are decomposed by base-symmetry angles. Symmetry operators copy and index the C{sub ij} chord to the form required for the particular symmetry. The use of asynchronous reads, lookup tables, and optimized image indexing improves computational performance.« less
  • Recent developments made MRI guided radiotherapy feasible. Performing simultaneous imaging during fractions can provide information about changing anatomy by means of deformable image registration for either immediate plan adaptations or accurate dose accumulation on the changing anatomy. In 3D MRI, however, acquisition time is considerable and scales with resolution. Furthermore, intra-scan motion degrades image quality.In this work, we investigate the sensitivity of registration quality on imageresolution: potentially, by employing spatial undersampling, the acquisition timeof MR images for the purpose of deformable image registration can be reducedsignificantly.On a volunteer, 3D-MR imaging data was sampled in a navigator-gated manner, acquiring one axialmore » volume (360×260×100mm{sup 3}) per 3s during exhale phase. A T1-weighted FFE sequence was used with an acquired voxel size of (2.5mm{sup 3}) for a duration of 17min. Deformation vector fields were evaluated for 100 imaging cycles with respect to the initial anatomy using deformable image registration based on optical flow. Subsequently, the imaging data was downsampled by a factor of 2, simulating a fourfold acquisition speed. Displacements of the downsampled volumes were then calculated by the same process.In kidneyliver boundaries and the region around stomach/duodenum, prominent organ drifts could be observed in both the original and the downsampled imaging data. An increasing displacement of approximately 2mm was observed for the kidney, while an area around the stomach showed sudden displacements of 4mm. Comparison of the motile points over time showed high reproducibility between the displacements of high-resolution and downsampled volumes: over a 17min acquisition, the componentwise RMS error was not more than 0.38mm.Based on the synthetic experiments, 3D nonrigid image registration shows little sensitivity to image resolution and the displacement information is preserved even when halving the resolution. This can be employed to greatly reduce image acquisition times for interventional applications in real-time. This work was funded by the SoRTS consortium, which includes the industry partners Elekta, Philips and Technolution.« less
  • Purpose: The present work evaluates an iterative reconstruction approach, namely, the ordered subsets convex (OSC) algorithm with regularization via total variation (TV) minimization in the field of cone-beam optical computed tomography (optical CT). One of the uses of optical CT is gel-based 3D dosimetry for radiation therapy, where it is employed to map dose distributions in radiosensitive gels. Model-based iterative reconstruction may improve optical CT image quality and contribute to a wider use of optical CT in clinical gel dosimetry. Methods: This algorithm was evaluated using experimental data acquired by a cone-beam optical CT system, as well as complementary numericalmore » simulations. A fast GPU implementation of OSC-TV was used to achieve reconstruction times comparable to those of conventional filtered backprojection. Images obtained via OSC-TV were compared with the corresponding filtered backprojections. Spatial resolution and uniformity phantoms were scanned and respective reconstructions were subject to evaluation of the modulation transfer function, image uniformity, and accuracy. The artifacts due to refraction and total signal loss from opaque objects were also studied. Results: The cone-beam optical CT data reconstructions showed that OSC-TV outperforms filtered backprojection in terms of image quality, thanks to a model-based simulation of the photon attenuation process. It was shown to significantly improve the image spatial resolution and reduce image noise. The accuracy of the estimation of linear attenuation coefficients remained similar to that obtained via filtered backprojection. Certain image artifacts due to opaque objects were reduced. Nevertheless, the common artifact due to the gel container walls could not be eliminated. Conclusions: The use of iterative reconstruction improves cone-beam optical CT image quality in many ways. The comparisons between OSC-TV and filtered backprojection presented in this paper demonstrate that OSC-TV can potentially improve the rendering of spatial features and reduce cone-beam optical CT artifacts.« less
  • Purpose: To validate Computed Tomography Fractional Flow Reserve (CT-FFR) measurements with accurate 3D printed coronary phantoms. Methods: DICOM data from four phases in two patients imaged with a standard 320 × 0.5mm coronary CT acquisition (70–80% cardiac cycle) underwent semi-automated segmentation using a research workstation. Both patients had a >50% stenosis from the clinical image interpretation. Each volume was saved as a Stereo Lithographic (STL) file with 250 micron resolution. The 3D geometries were qualitatively assessed; the best of the four phases was 3D printed using a Stratasys Eden260V printer in Tango+, a rubber-like material that roughly emulates mechanical propertiesmore » of human vasculature. We connected the model to a programmable pump and measured the pressure drop using pressure sensors embedded proximal and distal to the arterial stenosis. Next, the STL files used for the 3D printed models were uploaded in the ANSYS meshing tool (ICEM CFD 16.1). A standard meshing process was applied and the meshed geometry was directly imported in the ANSYS Fluent for Computational Flow Dynamics simulations. The CFD simulations were used to calculate the CT-FFR and compared to the bench top FFR measured in the 3D printed phantoms. Results: FFR-CT measurements and phantoms were completed in within an hour after the segmentation. Patient 1 had a 60% stenosis that resulted in a CT-FFR of 0.68. The second case had a 50% stenosis and a CT-FFR of 0.75. The average bench top FFR measurements were 0.72 and 0.80, respectively. Conclusion: This pilot investigation demonstrated the use of a bench-top coronary model for CT-FFR validation. The measurements and the CFD simulations agreed within 6%. Project supported by Support: Toshiba America Medical Systems Corp.and NIH grant R01-EB002873. Project supported by Toshiba America Medical Systems Corp.and partial support from NIH grant R01-EB002873.« less