1
|
Sengupta A, Badal A, Makeev A, Badano A. Computational models of direct and indirect X-ray breast imaging detectors for in silico trials. Med Phys 2022; 49:6856-6870. [PMID: 35997076 DOI: 10.1002/mp.15935] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 07/27/2022] [Accepted: 07/28/2022] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND To facilitate in silico studies that investigate digital mammography (DM) and breast tomosynthesis (DBT), models replicating the variety in imaging performance of the DM and DBT systems, observed across manufacturers are needed. PURPOSE The main purpose of this work is to develop generic physics models for direct and indirect detector technology used in commercially available systems, with the goal of making them available open source to manufacturers to further tweak and develop the exact in silico replicas of their systems. METHODS We recently reported on an in silico version of the SIEMENS Mammomat Inspiration DM/DBT system using an open-source GPU-accelerated Monte Carlo x-ray imaging simulation code (MC-GPU). We build on the previous version of the MC-GPU codes to mimic the imaging performances of two other Food and Drug Administration (FDA)-approved DM/DBT systems, such as Hologic Selenia Dimensions (HSD) and the General Electric Senographe Pristina (GSP) systems. In this work, we developed a hybrid technique to model the optical spread and signal crosstalk observed in the GSP and HSD systems. MC simulations are used to track each x-ray photon till its first interaction within the x-ray detector. On the other hand, the signal spread in the x-ray detectors is modeled using previously developed analytical equations. This approach allows us to preserve the modeling accuracy offered by MC methods in the patient body, while speeding up secondary carrier transport (either electron-hole pairs or optical photons) using analytical equations in the detector. The analytical optical spread model for the indirect detector includes the depth-dependent spread and collection of optical photons and relies on a pre-computed set of point response functions that describe the optical spread as a function of depth. To understand the capabilities of the computational x-ray detector models, we compared image quality metrics like modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE), simulated with our models against measured data. Please note that the purpose of these comparisons with measured data would be to gauge if the model developed as part of this work could replicate commercially used direct and indirect technology in general and not to achieve perfect fits with measured data. RESULTS We found that the simulated image quality metrics such as MTF, NNPS, and DQE were in reasonable agreement with experimental data. To demonstrate the imaging performance of the three DM/DBT systems, we integrated the detector models with the VICTRE pipeline and simulated DM images of a fatty breast model containing a spiculated mass and a calcium oxalate cluster. In general, we found that the images generated using the indirect model appeared more blurred with a different noise texture and contrast as compared to the systems with direct detectors. CONCLUSIONS We have presented computational models of three commercially available FDA-approved DM/DBT systems, which implement both direct and indirect detector technology. The updated versions of the MC-GPU codes that can be used to replicate three systems are available in open source format through GitHub.
Collapse
Affiliation(s)
- Aunnasha Sengupta
- Division of Imaging, Diagnostics and Software Reliability, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U. S. Food and Drug Administration, Silver Spring, Maryland, USA
| | - Andreu Badal
- Division of Imaging, Diagnostics and Software Reliability, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U. S. Food and Drug Administration, Silver Spring, Maryland, USA
| | - Andrey Makeev
- Division of Imaging, Diagnostics and Software Reliability, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U. S. Food and Drug Administration, Silver Spring, Maryland, USA
| | - Aldo Badano
- Division of Imaging, Diagnostics and Software Reliability, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U. S. Food and Drug Administration, Silver Spring, Maryland, USA
| |
Collapse
|
2
|
Virtual Clinical Trials in 2D and 3D X-ray Breast Imaging and Dosimetry: Comparison of CPU-Based and GPU-Based Monte Carlo Codes. Cancers (Basel) 2022; 14:cancers14041027. [PMID: 35205775 PMCID: PMC8870739 DOI: 10.3390/cancers14041027] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2021] [Revised: 02/12/2022] [Accepted: 02/13/2022] [Indexed: 12/10/2022] Open
Abstract
Simple Summary Virtual clinical trials in X-ray breast imaging may permit substantial reduction of the costs, times, and exposure risk to patient of clinical trials. Monte Carlo simulation techniques are increasingly adopted for VCT in breast imaging and dosimetry studies. This work aims to compare three different platforms for breast VCT studies, to develop real-time virtual DM, DBT and BCT examinations, where the in-silico image acquisition process takes a computational time comparable to that typical of a corresponding real clinical examination. Abstract Computational reproductions of medical imaging tests, a form of virtual clinical trials (VCTs), are increasingly being used, particularly in breast imaging research. The accuracy of the computational platform that is used for the imaging and dosimetry simulation processes is a fundamental requirement. Moreover, for practical usage, the imaging simulation computation time should be compatible with the clinical workflow. We compared three different platforms for in-silico X-ray 3D breast imaging: the Agata (University & INFN Napoli) that was based on the Geant4 toolkit and running on a CPU-based server architecture; the XRMC Monte Carlo (University of Cagliari) that was based on the use of variance reduction techniques, running on a CPU hardware; and the Monte Carlo code gCTD (University of Texas Southwestern Medical Center) running on a single GPU platform with CUDA environment. The tests simulated the irradiation of cylindrical objects as well as anthropomorphic breast phantoms and produced 2D and 3D images and 3D maps of absorbed dose. All the codes showed compatible results in terms of simulated dose maps and imaging values within a maximum discrepancy of 3%. The GPU-based code produced a reduction of the computation time up to factor 104, and so permits real-time VCT studies for X-ray breast imaging.
Collapse
|
3
|
Volken W, Zulliger MA, Koller B, Manser P, Fix MK. Investigation on the resolution of a micro cone beam CT scanner scintillating detector using Monte Carlo methods. Phys Med 2018; 53:17-24. [PMID: 30241750 DOI: 10.1016/j.ejmp.2018.08.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/27/2018] [Revised: 06/13/2018] [Accepted: 08/05/2018] [Indexed: 11/19/2022] Open
Abstract
The impact of several physical quantities on the spatial resolution of an X-ray scintillating pixel detector for a micro cone beam CT (µCBCT) is investigated and discussed. The XtremeCT from SCANCO Medical AG was simulated using the EGSnrc/EGS++ Monte Carlo (MC) framework and extensively benchmarked in a previous work. The resolution of the detector was determined by simulating a titanium knife-edge to obtain the edge spread function (ESF) and the modulation transfer function (MTF). Propagation of the scintillation light through the scintillator and its coupling into the fiber optics system was taken into account. The contribution of particles scattered in the main scanner components to the detector signal is very low and does not affect the spatial resolution of the detector. The resolution obtained from the energy deposition in the scintillator without any blurring due to the propagation of the scintillation light into the fiber optics array was 31 µm. By assuming isotropic light propagation in the scintillator, the resolution degraded to 360 µm. A simple light propagation model taking into account the impact of the scintillator's columnar microstructures was developed and compared with the MANTIS Monte Carlo simulation package. By reducing the width of the model's light propagation kernel by a factor of 2 compared to the isotropic case, the detector resolution can be improved to 83 µm, which corresponds well to the measured resolution of 86 µm. The resolution of the detector is limited mainly by the propagation of the scintillation light through the scintillator layer. It offers the greatest potential to improve the resolution of the µCBCT imaging system.
Collapse
Affiliation(s)
- W Volken
- Division of Medical Radiation Physics and Department of Radiation Oncology, Inselspital, Bern University Hospital, and University of Bern, Switzerland
| | | | - B Koller
- SCANCO Medical AG, Brüttisellen, Switzerland
| | - P Manser
- Division of Medical Radiation Physics and Department of Radiation Oncology, Inselspital, Bern University Hospital, and University of Bern, Switzerland
| | - M K Fix
- Division of Medical Radiation Physics and Department of Radiation Oncology, Inselspital, Bern University Hospital, and University of Bern, Switzerland.
| |
Collapse
|
4
|
Sharma D, Sempau J, Badano A. Technical Note: On the efficiency of variance reduction techniques for Monte Carlo estimates of imaging noise. Med Phys 2017; 45:629-634. [DOI: 10.1002/mp.12638] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2017] [Revised: 08/30/2017] [Accepted: 10/15/2017] [Indexed: 02/06/2023] Open
Affiliation(s)
- Diksha Sharma
- Division of Imaging, Diagnostics, and Software Reliability; OSEL/CDRH; U.S. Food and Drug Administration; Silver Spring MD USA
| | - Josep Sempau
- Institute of Energy Technologies; Polytechnic University of Catalonia; Barcelona Spain
| | - Aldo Badano
- Division of Imaging, Diagnostics, and Software Reliability; OSEL/CDRH; U.S. Food and Drug Administration; Silver Spring MD USA
| |
Collapse
|
5
|
Roncali E, Mosleh-Shirazi MA, Badano A. Modelling the transport of optical photons in scintillation detectors for diagnostic and radiotherapy imaging. Phys Med Biol 2017; 62:R207-R235. [PMID: 28976914 PMCID: PMC5739055 DOI: 10.1088/1361-6560/aa8b31] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Computational modelling of radiation transport can enhance the understanding of the relative importance of individual processes involved in imaging systems. Modelling is a powerful tool for improving detector designs in ways that are impractical or impossible to achieve through experimental measurements. Modelling of light transport in scintillation detectors used in radiology and radiotherapy imaging that rely on the detection of visible light plays an increasingly important role in detector design. Historically, researchers have invested heavily in modelling the transport of ionizing radiation while light transport is often ignored or coarsely modelled. Due to the complexity of existing light transport simulation tools and the breadth of custom codes developed by users, light transport studies are seldom fully exploited and have not reached their full potential. This topical review aims at providing an overview of the methods employed in freely available and other described optical Monte Carlo packages and analytical models and discussing their respective advantages and limitations. In particular, applications of optical transport modelling in nuclear medicine, diagnostic and radiotherapy imaging are described. A discussion on the evolution of these modelling tools into future developments and applications is presented.
Collapse
Affiliation(s)
- Emilie Roncali
- Department of Biomedical Engineering, University of California Davis, Davis, USA
| | - Mohammad Amin Mosleh-Shirazi
- Medical Imaging Research Center, and, Physics Unit, Department of Radiotherapy and Oncology, Namazi Hospital, Shiraz University of Medical Sciences, Shiraz 71936-13311, Iran
| | - Aldo Badano
- Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20852, USA
| |
Collapse
|
6
|
Maurino SL, Badano A, Cunningham IA, Karim KS. Theoretical and Monte Carlo optimization of a stacked three-layer flat-panel x-ray imager for applications in multi-spectral medical imaging. PROCEEDINGS OF SPIE--THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING 2016; 9783:97833Z. [PMID: 28845080 PMCID: PMC5568811 DOI: 10.1117/12.2217085] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
We propose a new design of a stacked three-layer flat-panel x-ray detector for dual-energy (DE) imaging. Each layer consists of its own scintillator of individual thickness and an underlying thin-film-transistor-based flat-panel. Three images are obtained simultaneously in the detector during the same x-ray exposure, thereby eliminating any motion artifacts. The detector operation is two-fold: a conventional radiography image can be obtained by combining all three layers' images, while a DE subtraction image can be obtained from the front and back layers' images, where the middle layer acts as a mid-filter that helps achieve spectral separation. We proceed to optimize the detector parameters for two sample imaging tasks that could particularly benefit from this new detector by obtaining the best possible signal to noise ratio per root entrance exposure using well-established theoretical models adapted to fit our new design. These results are compared to a conventional DE temporal subtraction detector and a single-shot DE subtraction detector with a copper mid-filter, both of which underwent the same theoretical optimization. The findings are then validated using advanced Monte Carlo simulations for all optimized detector setups. Given the performance expected from initial results and the recent decrease in price for digital x-ray detectors, the simplicity of the three-layer stacked imager approach appears promising to usher in a new generation of multi-spectral digital x-ray diagnostics.
Collapse
Affiliation(s)
| | - Aldo Badano
- Center for Devices and Radiological Health, Food and Drug Administration, 10903 New Hampshire Ave, Silver Spring, MD 20993, USA
| | - Ian A Cunningham
- Imaging Research Laboratories, Robarts Research Institute, and Department of Medical Biophysics, The University of Western Ontario, London, ON N6A 5B7, Canada
| | - Karim S Karim
- Electrical and Computer Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
- Centre for Bioengineering and Biotechnology, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| |
Collapse
|
7
|
Zhao C, Konstantinidis AC, Zheng Y, Anaxagoras T, Speller RD, Kanicki J. 50μm pixel pitch wafer-scale CMOS active pixel sensor x-ray detector for digital breast tomosynthesis. Phys Med Biol 2015; 60:8977-9001. [DOI: 10.1088/0031-9155/60/23/8977] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
8
|
Nillius P, Klamra W, Sibczynski P, Sharma D, Danielsson M, Badano A. Light output measurements and computational models of microcolumnar CsI scintillators for x-ray imaging. Med Phys 2015; 42:600-605. [DOI: 10.1118/1.4905096] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2014] [Revised: 11/16/2014] [Accepted: 12/11/2014] [Indexed: 01/12/2023] Open
|
9
|
Dong H, Sharma D, Badano A. Web-based, GPU-accelerated, Monte Carlo simulation and visualization of indirect radiation imaging detector performance. Med Phys 2014; 41:121907. [PMID: 25471967 DOI: 10.1118/1.4901516] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Monte Carlo simulations play a vital role in the understanding of the fundamental limitations, design, and optimization of existing and emerging medical imaging systems. Efforts in this area have resulted in the development of a wide variety of open-source software packages. One such package, hybridmantis, uses a novel hybrid concept to model indirect scintillator detectors by balancing the computational load using dual CPU and graphics processing unit (GPU) processors, obtaining computational efficiency with reasonable accuracy. In this work, the authors describe two open-source visualization interfaces, webmantis and visualmantis to facilitate the setup of computational experiments via hybridmantis. METHODS The visualization tools visualmantis and webmantis enable the user to control simulation properties through a user interface. In the case of webmantis, control via a web browser allows access through mobile devices such as smartphones or tablets. webmantis acts as a server back-end and communicates with an NVIDIA GPU computing cluster that can support multiuser environments where users can execute different experiments in parallel. RESULTS The output consists of point response and pulse-height spectrum, and optical transport statistics generated by hybridmantis. The users can download the output images and statistics through a zip file for future reference. In addition, webmantis provides a visualization window that displays a few selected optical photon path as they get transported through the detector columns and allows the user to trace the history of the optical photons. CONCLUSIONS The visualization tools visualmantis and webmantis provide features such as on the fly generation of pulse-height spectra and response functions for microcolumnar x-ray imagers while allowing users to save simulation parameters and results from prior experiments. The graphical interfaces simplify the simulation setup and allow the user to go directly from specifying input parameters to receiving visual feedback for the model predictions.
Collapse
Affiliation(s)
- Han Dong
- Division of Imaging, Diagnostics, and Software Reliability, Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, Maryland 20993
| | - Diksha Sharma
- Division of Imaging, Diagnostics, and Software Reliability, Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, Maryland 20993
| | - Aldo Badano
- Division of Imaging, Diagnostics, and Software Reliability, Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, Maryland 20993
| |
Collapse
|
10
|
Gutierrez B, Badano A, Samuelson F. Analytic variance estimates of Swank and Fano factors. Med Phys 2014; 41:072102. [PMID: 24989397 DOI: 10.1118/1.4881328] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
PURPOSE Variance estimates for detector energy resolution metrics can be used as stopping criteria in Monte Carlo simulations for the purpose of ensuring a small uncertainty of those metrics and for the design of variance reduction techniques. METHODS The authors derive an estimate for the variance of two energy resolution metrics, the Swank factor and the Fano factor, in terms of statistical moments that can be accumulated without significant computational overhead. The authors examine the accuracy of these two estimators and demonstrate how the estimates of the coefficient of variation of the Swank and Fano factors behave with data from a Monte Carlo simulation of an indirect x-ray imaging detector. RESULTS The authors' analyses suggest that the accuracy of their variance estimators is appropriate for estimating the actual variances of the Swank and Fano factors for a variety of distributions of detector outputs. CONCLUSIONS The variance estimators derived in this work provide a computationally convenient way to estimate the error or coefficient of variation of the Swank and Fano factors during Monte Carlo simulations of radiation imaging systems.
Collapse
Affiliation(s)
| | - Aldo Badano
- US Food and Drug Administration, Silver Spring, Maryland 20993
| | - Frank Samuelson
- US Food and Drug Administration, Silver Spring, Maryland 20993
| |
Collapse
|
11
|
Sharma D, Badano A. Validation of columnar CsI x-ray detector responses obtained with hybridMANTIS, a CPU-GPU Monte Carlo code for coupled x-ray, electron, and optical transport. Med Phys 2013; 40:031907. [PMID: 23464322 DOI: 10.1118/1.4791642] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
PURPOSE hybridMANTIS is a Monte Carlo package for modeling indirect x-ray imagers using columnar geometry based on a hybrid concept that maximizes the utilization of available CPU and graphics processing unit processors in a workstation. METHODS The authors compare hybridMANTIS x-ray response simulations to previously published MANTIS and experimental data for four cesium iodide scintillator screens. These screens have a variety of reflective and absorptive surfaces with different thicknesses. The authors analyze hybridMANTIS results in terms of modulation transfer function and calculate the root mean square difference and Swank factors from simulated and experimental results. RESULTS The comparison suggests that hybridMANTIS better matches the experimental data as compared to MANTIS, especially at high spatial frequencies and for the thicker screens. hybridMANTIS simulations are much faster than MANTIS with speed-ups up to 5260. CONCLUSIONS hybridMANTIS is a useful tool for improved description and optimization of image acquisition stages in medical imaging systems and for modeling the forward problem in iterative reconstruction algorithms.
Collapse
Affiliation(s)
- Diksha Sharma
- Division of Imaging and Applied Mathematics, Center for Devices and Radiological Health, Food and Drug Administration, Silver Spring, Maryland 20993, USA
| | | |
Collapse
|