Baiyu Chen1, Erich Kobler2, Matthew J Muckley1, Aaron D Sodickson3, Thomas O'Donnell4, Thomas Flohr5, Bernhard Schmidt5, Daniel K Sodickson1, Ricardo Otazo1,6,7. 1. Department of Radiology, NYU School of Medicine, New York, NY, 10016, USA. 2. Institute of Computer Graphics and Vision, Graz University of Technology, Graz, 8010, Austria. 3. Department of Radiology, Brigham and Women's Hospital, Boston, MA, 02115, USA. 4. Siemens Healthineers, Malvern, PA, 19355, USA. 5. Siemens Healthineers, Forchheim, 91301, Germany. 6. Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA. 7. Department of Radiology, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA.
Abstract
PURPOSE: SparseCT, an undersampling scheme for compressed sensing (CS) computed tomography (CT), has been proposed to reduce radiation dose by acquiring undersampled projection data from clinical CT scanners (Koesters et al. in, SparseCT: Interrupted-Beam Acquisition and Sparse Reconstruction for Radiation Dose Reduction; 2017). SparseCT partially blocks the x-ray beam with a multislit collimator (MSC) to perform a multidimensional undersampling along the view and detector row dimensions. SparseCT undersamples the projection data within each view and moves the MSC along the z-direction during gantry rotation to change the undersampling pattern. It enables reconstruction of images from undersampled data using CS algorithms. The purpose of this work is to design the spacing and width of the MSC slits and the MSC motion patterns based on beam separation, undersampling efficiency, and image quality. The development and testing of a SparseCT prototype with the designed MSC will be described in a following paper. METHODS: We chose a few initial MSC designs based on the guidance from two metrics: beam separation and undersampling efficiency. Both beam separation and undersampling efficiency were measured from numerically simulated photon distribution with MSC taken into consideration. Beam separation measures the separation between x-ray beams from consecutive slits, taking into account penumbra effects on both sides of each slit. Undersampling efficiency measures the dose-weighted similarity between penumbra undersampling and binary undersampling, in other words, the effective contribution of the incident dose to the signal to noise ratio of the projection data. We then compared the initially chosen MSC designs in terms of their reconstruction image quality. SparseCT projections were simulated from fully sampled patient projection data according to the MSC design and motion pattern, reconstructed iteratively using a sparsity-enforcing penalized weighted least squares cost function with ordered subsets/momentum algorithm, and compared visually and quantitatively. RESULTS: Simulated photon distributions indicate that the size of the penumbra is dominated by the size of the focal spot. Therefore, a wider MSC slit and a smaller focal spot lead to increased beam separation and undersampling efficiency. For fourfold undersampling with a 1.2 mm focal spot, a minimum MSC slit width of three detector rows (projected to the detector surface) is needed for beam separation; for threefold undersampling, a minimum slit width of four detector rows is needed. Simulations of SparseCT projection and reconstruction indicate that the motion pattern of the MSC does not have a visible impact on image quality. An MSC slit width of three or four detector rows yields similar image quality. CONCLUSION: The MSC is the key component of the SparseCT method. Simulations of MSC designs incorporating x-ray beam penumbra effects showed that for threefold and fourfold dose reductions, an MSC slit width of four detector rows provided reasonable beam separation, undersampling efficiency, and image quality.
PURPOSE: SparseCT, an undersampling scheme for compressed sensing (CS) computed tomography (CT), has been proposed to reduce radiation dose by acquiring undersampled projection data from clinical CT scanners (Koesters et al. in, SparseCT: Interrupted-Beam Acquisition and Sparse Reconstruction for Radiation Dose Reduction; 2017). SparseCT partially blocks the x-ray beam with a multislit collimator (MSC) to perform a multidimensional undersampling along the view and detector row dimensions. SparseCT undersamples the projection data within each view and moves the MSC along the z-direction during gantry rotation to change the undersampling pattern. It enables reconstruction of images from undersampled data using CS algorithms. The purpose of this work is to design the spacing and width of the MSC slits and the MSC motion patterns based on beam separation, undersampling efficiency, and image quality. The development and testing of a SparseCT prototype with the designed MSC will be described in a following paper. METHODS: We chose a few initial MSC designs based on the guidance from two metrics: beam separation and undersampling efficiency. Both beam separation and undersampling efficiency were measured from numerically simulated photon distribution with MSC taken into consideration. Beam separation measures the separation between x-ray beams from consecutive slits, taking into account penumbra effects on both sides of each slit. Undersampling efficiency measures the dose-weighted similarity between penumbra undersampling and binary undersampling, in other words, the effective contribution of the incident dose to the signal to noise ratio of the projection data. We then compared the initially chosen MSC designs in terms of their reconstruction image quality. SparseCT projections were simulated from fully sampled patient projection data according to the MSC design and motion pattern, reconstructed iteratively using a sparsity-enforcing penalized weighted least squares cost function with ordered subsets/momentum algorithm, and compared visually and quantitatively. RESULTS: Simulated photon distributions indicate that the size of the penumbra is dominated by the size of the focal spot. Therefore, a wider MSC slit and a smaller focal spot lead to increased beam separation and undersampling efficiency. For fourfold undersampling with a 1.2 mm focal spot, a minimum MSC slit width of three detector rows (projected to the detector surface) is needed for beam separation; for threefold undersampling, a minimum slit width of four detector rows is needed. Simulations of SparseCT projection and reconstruction indicate that the motion pattern of the MSC does not have a visible impact on image quality. An MSC slit width of three or four detector rows yields similar image quality. CONCLUSION: The MSC is the key component of the SparseCT method. Simulations of MSC designs incorporating x-ray beam penumbra effects showed that for threefold and fourfold dose reductions, an MSC slit width of four detector rows provided reasonable beam separation, undersampling efficiency, and image quality.
Authors: Lifeng Yu; Michael R Bruesewitz; Kristen B Thomas; Joel G Fletcher; James M Kofler; Cynthia H McCollough Journal: Radiographics Date: 2011 May-Jun Impact factor: 5.333
Authors: Jodie A Christner; Vanessa A Zavaletta; Christian D Eusemann; Alisa I Walz-Flannigan; Cynthia H McCollough Journal: AJR Am J Roentgenol Date: 2010-01 Impact factor: 3.959
Authors: Matthew J Muckley; Baiyu Chen; Thomas Vahle; Thomas O'Donnell; Florian Knoll; Aaron D Sodickson; Daniel K Sodickson; Ricardo Otazo Journal: Phys Med Biol Date: 2019-08-07 Impact factor: 3.609