1
|
Tanveer MS, Wiedeman C, Li M, Shi Y, De Man B, Maltz JS, Wang G. Deep-silicon photon-counting x-ray projection denoising through reinforcement learning. JOURNAL OF X-RAY SCIENCE AND TECHNOLOGY 2024; 32:173-205. [PMID: 38217633 DOI: 10.3233/xst-230278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2024]
Abstract
BACKGROUND In recent years, deep reinforcement learning (RL) has been applied to various medical tasks and produced encouraging results. OBJECTIVE In this paper, we demonstrate the feasibility of deep RL for denoising simulated deep-silicon photon-counting CT (PCCT) data in both full and interior scan modes. PCCT offers higher spatial and spectral resolution than conventional CT, requiring advanced denoising methods to suppress noise increase. METHODS In this work, we apply a dueling double deep Q network (DDDQN) to denoise PCCT data for maximum contrast-to-noise ratio (CNR) and a multi-agent approach to handle data non-stationarity. RESULTS Using our method, we obtained significant image quality improvement for single-channel scans and consistent improvement for all three channels of multichannel scans. For the single-channel interior scans, the PSNR (dB) and SSIM increased from 33.4078 and 0.9165 to 37.4167 and 0.9790 respectively. For the multichannel interior scans, the channel-wise PSNR (dB) increased from 31.2348, 30.7114, and 30.4667 to 31.6182, 30.9783, and 30.8427 respectively. Similarly, the SSIM improved from 0.9415, 0.9445, and 0.9336 to 0.9504, 0.9493, and 0.0326 respectively. CONCLUSIONS Our results show that the RL approach improves image quality effectively, efficiently, and consistently across multiple spectral channels and has great potential in clinical applications.
Collapse
Affiliation(s)
- Md Sayed Tanveer
- Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Christopher Wiedeman
- Department of Electrical, Computer, and Systems Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Mengzhou Li
- Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Yongyi Shi
- Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Bruno De Man
- GE HealthCare, One Research Circle, Niskayuna, NY, USA
| | | | - Ge Wang
- Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| |
Collapse
|
2
|
Foldyna B, Basmagi S, Zangeneh FA, Wagner M, Doktorov K, Matveeva A, Denecke T, Gohmann RF, Lücke C, Gutberlet M, Lehmkuhl L. CT-derived coronary artery calcium density is affected by regional lesion distribution and image reconstruction parameters. Clin Imaging 2023; 103:109980. [PMID: 37677856 DOI: 10.1016/j.clinimag.2023.109980] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Revised: 08/17/2023] [Accepted: 08/28/2023] [Indexed: 09/09/2023]
Abstract
PURPOSE The prognostic relevance of coronary artery calcium (CAC) density, assessed from cardiac CT scans, is established. However, the influence of CAC distribution, volume, image reconstruction, and clinical factors on CAC density warrants further examination. METHODS In this study, 120 patients underwent non-contrast ECG-gated cardiac CT scans using a prospectively defined CAC scoring protocol with 1-, 3-, and 5-mm thick image reconstructions, both with and without a 20% image overlap. We segmented CAC in all reconstructions and assessed the relationship between CAC density, volume, and number of detected calcifications/patient. RESULTS Overall, 75/120 (63%) patients (66% men, mean age 63 ± 11 years) presented CAC across 342 segments. CAC density, CAC volume, and the number of detected calcifications decreased with increasing slice thickness (p < 0.001 for all); these effects were slightly reduced by image overlap (p < 0.001 for all). Higher CAC density correlated with greater CAC volume (ρ = 0.62; p < 0.001) and more calcified segments per person (ρ = 0.32; p = 0.006). Higher CAC density was also associated with lower patient weight (beta: -0.6, 95%CI: -1.1--0.1, p = 0.022) and increased high-density lipoprotein (HDL) levels (beta: 0.7, 95%CI: 0.0-1.4, p = 0.046). In a multivariable analysis adjusted for clinical covariates, lower CAC density was associated with broader CAC distribution (i.e., a higher number of calcified segments at a given CAC volume; beta-coefficient: -58.9; 95%CI: -84.7 to -33.1; p < 0.001). CONCLUSION CAC density is significantly impacted by regional CAC distribution and image reconstruction, potentially confounding its prognostic value. Accounting for these factors may improve patient risk assessment, management, and cardiovascular health outcomes.
Collapse
Affiliation(s)
- Borek Foldyna
- Cardiovascular Imaging Research Center, Massachusetts General Hospital - Harvard Medical School, 165 Cambridge Street, Suite 400, 02114 Boston, USA; Clinic for Radiology, Heart Center Bad Neustadt a.d. Saale, Bad Neustadt a.d. Saale, Germany.
| | - Said Basmagi
- Clinic for Radiology, Heart Center Bad Neustadt a.d. Saale, Bad Neustadt a.d. Saale, Germany
| | | | - Matthias Wagner
- Clinic for Radiology, Heart Center Bad Neustadt a.d. Saale, Bad Neustadt a.d. Saale, Germany
| | - Kalin Doktorov
- Clinic for Radiology, Heart Center Bad Neustadt a.d. Saale, Bad Neustadt a.d. Saale, Germany
| | - Anna Matveeva
- Clinic for Radiology, Heart Center Bad Neustadt a.d. Saale, Bad Neustadt a.d. Saale, Germany
| | - Timm Denecke
- Clinic for Diagnostic and Interventional Radiology, University of Leipzig, Leipzig, Germany
| | - Robin F Gohmann
- Department of Diagnostic and Interventional Radiology, Heart Center Leipzig at University of Leipzig, Leipzig, Germany
| | - Christian Lücke
- Department of Diagnostic and Interventional Radiology, Heart Center Leipzig at University of Leipzig, Leipzig, Germany
| | - Matthias Gutberlet
- Department of Diagnostic and Interventional Radiology, Heart Center Leipzig at University of Leipzig, Leipzig, Germany
| | - Lukas Lehmkuhl
- Clinic for Radiology, Heart Center Bad Neustadt a.d. Saale, Bad Neustadt a.d. Saale, Germany
| |
Collapse
|
3
|
Inkinen SI, Juntunen MAK, Ketola J, Korhonen K, Sepponen P, Kotiaho A, Pohjanen VM, Nieminen M. Virtual monochromatic imaging reduces beam hardening artefacts in cardiac interior photon counting computed tomography: a phantom study with cadaveric specimens. Biomed Phys Eng Express 2021; 8. [PMID: 34911047 DOI: 10.1088/2057-1976/ac4397] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2021] [Accepted: 12/15/2021] [Indexed: 11/11/2022]
Abstract
In interior cardiac computed tomography (CT) imaging, the x-ray beam is collimated to a limited field-of-view covering the heart volume, which decreases the radiation exposure to surrounding tissues. Spectral CT enables the creation of virtual monochromatic images (VMIs) through a computational material decomposition process. This study investigates the utility of VMIs for beam hardening (BH) reduction in interior cardiac CT, and further, the suitability of VMIs for coronary artery calcium (CAC) scoring and volume assessment is studied using spectral photon counting detector CT (PCD-CT).Ex vivocoronary artery samples (N = 18) were inserted in an epoxy rod phantom. The rod was scanned in the conventional CT geometry, and subsequently, the rod was positioned in a torso phantom and re-measured in the interior PCD-CT geometry. The total energy (TE) 10-100 keV reconstructions from PCD-CT were used as a reference. The low energy 10-60 keV and high energy 60-100 keV data were used to perform projection domain material decomposition to polymethyl methacrylate and calcium hydroxylapatite basis. The truncated basis-material sinograms were extended using the adaptive detruncation method. VMIs from 30-180 keV range were computed from the detruncated virtual monochromatic sinograms using filtered back projection. Detrending was applied as a post-processing method prior to CAC scoring. The results showed that BH artefacts from the exterior structures can be suppressed with high (≥100 keV) VMIs. With appropriate selection of the monoenergy (46 keV), the underestimation trend of CAC scores and volumes shown in Bland-Altman (BA) plots for TE interior PCD-CT was mitigated, as the BA slope values were -0.02 for the 46 keV VMI compared to -0.21 the conventional TE image. To conclude, spectral PCD-CT imaging using VMIs could be applied to reduce BH artefacts interior CT geometry, and further, optimal selection of VMI may improve the accuracy of CAC scoring assessment in interior PCD-CT.
Collapse
Affiliation(s)
- Satu I Inkinen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland
| | - Mikael A K Juntunen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland.,Oulu University Hospital, Department of Diagnostic Radiology, Oulu, Finland
| | - Juuso Ketola
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland.,The South Savo Social and Health Care Authority, Mikkeli Central Hospital, Mikkeli, Finland
| | - Kristiina Korhonen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland
| | - Pasi Sepponen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland
| | - Antti Kotiaho
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland.,Oulu University Hospital, Department of Diagnostic Radiology, Oulu, Finland
| | - Vesa-Matti Pohjanen
- Cancer and Translational Medicine Research Unit, Medical Research Center, University of Oulu and Oulu University Hospital, Oulu, Finland
| | - Miika Nieminen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland.,Oulu University Hospital, Department of Diagnostic Radiology, Oulu, Finland
| |
Collapse
|
4
|
Ketola JHJ, Heino H, Juntunen MAK, Nieminen MT, Siltanen S, Inkinen SI. Generative adversarial networks improve interior computed tomography angiography reconstruction. Biomed Phys Eng Express 2021; 7. [PMID: 34673559 DOI: 10.1088/2057-1976/ac31cb] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Accepted: 10/21/2021] [Indexed: 11/12/2022]
Abstract
In interior computed tomography (CT), the x-ray beam is collimated to a limited field-of-view (FOV) (e.g. the volume of the heart) to decrease exposure to adjacent organs, but the resulting image has a severe truncation artifact when reconstructed with traditional filtered back-projection (FBP) type algorithms. In some examinations, such as cardiac or dentomaxillofacial imaging, interior CT could be used to achieve further dose reductions. In this work, we describe a deep learning (DL) method to obtain artifact-free images from interior CT angiography. Our method employs the Pix2Pix generative adversarial network (GAN) in a two-stage process: (1) An extended sinogram is computed from a truncated sinogram with one GAN model, and (2) the FBP reconstruction obtained from that extended sinogram is used as an input to another GAN model that improves the quality of the interior reconstruction. Our double GAN (DGAN) model was trained with 10 000 truncated sinograms simulated from real computed tomography angiography slice images. Truncated sinograms (input) were used with original slice images (target) in training to yield an improved reconstruction (output). DGAN performance was compared with the adaptive de-truncation method, total variation regularization, and two reference DL methods: FBPConvNet, and U-Net-based sinogram extension (ES-UNet). Our DGAN method and ES-UNet yielded the best root-mean-squared error (RMSE) (0.03 ± 0.01), and structural similarity index (SSIM) (0.92 ± 0.02) values, and reference DL methods also yielded good results. Furthermore, we performed an extended FOV analysis by increasing the reconstruction area by 10% and 20%. In both cases, the DGAN approach yielded best results at RMSE (0.03 ± 0.01 and 0.04 ± 0.01 for the 10% and 20% cases, respectively), peak signal-to-noise ratio (PSNR) (30.5 ± 2.6 dB and 28.6 ± 2.6 dB), and SSIM (0.90 ± 0.02 and 0.87 ± 0.02). In conclusion, our method was able to not only reconstruct the interior region with improved image quality, but also extend the reconstructed FOV by 20%.
Collapse
Affiliation(s)
- Juuso H J Ketola
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, FI-90014, Finland.,The South Savo Social and Health Care Authority, Mikkeli Central Hospital, FI-50100, Finland
| | - Helinä Heino
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, FI-90014, Finland
| | - Mikael A K Juntunen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, FI-90014, Finland.,Department of Diagnostic Radiology, Oulu University Hospital, FI-90029, Finland
| | - Miika T Nieminen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, FI-90014, Finland.,Department of Diagnostic Radiology, Oulu University Hospital, FI-90029, Finland.,Medical Research Center Oulu, University of Oulu and Oulu University Hospital, FI-90014, Finland
| | - Samuli Siltanen
- Department of Mathematics and Statistics, University of Helsinki, Helsinki, FI-00014, Finland
| | - Satu I Inkinen
- Research Unit of Medical Imaging, Physics and Technology, University of Oulu, FI-90014, Finland
| |
Collapse
|
5
|
Juntunen MAK, Kotiaho AO, Nieminen MT, Inkinen SI. Optimizing iterative reconstruction for quantification of calcium hydroxyapatite with photon counting flat-detector computed tomography: a cardiac phantom study. J Med Imaging (Bellingham) 2021; 8:052102. [PMID: 33718518 PMCID: PMC7946398 DOI: 10.1117/1.jmi.8.5.052102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2020] [Accepted: 01/28/2021] [Indexed: 11/28/2022] Open
Abstract
Purpose: Coronary artery calcium (CAC) scoring with computed tomography (CT) has been proposed as a screening tool for coronary artery disease, but concerns remain regarding the radiation dose of CT CAC scoring. Photon counting detectors and iterative reconstruction (IR) are promising approaches for patient dose reduction, yet the preservation of CAC scores with IR has been questioned. The purpose of this study was to investigate the applicability of IR for quantification of CAC using a photon counting flat-detector. Approach: We imaged a cardiac rod phantom with calcium hydroxyapatite (CaHA) inserts with different noise levels using an experimental photon counting flat-detector CT setup to simulate the clinical CAC scoring protocol. We applied filtered back projection (FBP) and two IR algorithms with different regularization strengths. We compared the air kerma values, image quality parameters [noise magnitude, noise power spectrum, modulation transfer function (MTF), and contrast-to-noise ratio], and CaHA quantification accuracy between FBP and IR. Results: IR regularization strength influenced CAC scores significantly ( p < 0.05 ). The CAC volumes and scores between FBP and IRs were the most similar when the IR regularization strength was chosen to match the MTF of the FBP reconstruction. Conclusion: When the regularization strength is selected to produce comparable spatial resolution with FBP, IR can yield comparable CAC scores and volumes with FBP. Nonetheless, at the lowest radiation dose setting, FBP produced more accurate CAC volumes and scores compared to IR, and no improved CAC scoring accuracy at low dose was demonstrated with the utilized IR methods.
Collapse
Affiliation(s)
- Mikael A. K. Juntunen
- University of Oulu, Research Unit of Medical Imaging, Physics, and Technology, Oulu, Finland
- Oulu University Hospital, Department of Diagnostic Radiology, Oulu, Finland
| | - Antti O. Kotiaho
- Oulu University Hospital, Department of Diagnostic Radiology, Oulu, Finland
| | - Miika T. Nieminen
- University of Oulu, Research Unit of Medical Imaging, Physics, and Technology, Oulu, Finland
- Oulu University Hospital, Department of Diagnostic Radiology, Oulu, Finland
- Medical Research Center, University of Oulu, Oulu University Hospital, Oulu, Finland
| | - Satu I. Inkinen
- University of Oulu, Research Unit of Medical Imaging, Physics, and Technology, Oulu, Finland
| |
Collapse
|