Virtual Clinical Trials in 2D and 3D X-ray Breast Imaging and Dosimetry: Comparison of CPU-Based and GPU-Based Monte Carlo Codes

Virtual dosimetry study with three cone-beam breast computed tomography scanners using a fast GPU-based Monte Carlo code

Objective. To compare the dosimetric performance of three cone-beam breast computed tomography (BCT) scanners, using real-time Monte Carlo-based dose estimates obtained with the virtual clinical trials (VCT)-BREAST graphical processing unit (GPU)-accelerated platform dedicated to VCT in breast imaging. Approach. A GPU-based Monte Carlo (MC) code was developed for replicatingin silicothe geometric, x-ray spectra and detector setups adopted, respectively, in two research scanners and one commercial BCT scanner, adopting 80 kV, 60 kV and 49 kV tube voltage, respectively. Our cohort of virtual breasts included 16 anthropomorphic voxelized breast phantoms from a publicly available dataset. For each virtual patient, we simulated exams on the three scanners, up to a nominal simulated mean glandular dose of 5 mGy (primary photons launched, in the order of 1011-1012per scan). Simulated 3D dose maps (recorded for skin, adipose and glandular tissues) were compared for the same phantom, on the three scanners. MC simulations were implemented on a single NVIDIA GeForce RTX 3090 graphics card.Main results.Using the spread of the dose distribution as a figure of merit, we showed that, in the investigated phantoms, the glandular dose is more uniform within less dense breasts, and it is more uniformly distributed for scans at 80 kV and 60 kV, than at 49 kV. A realistic virtual study of each breast phantom was completed in about 3.0 h with less than 1% statistical uncertainty, with 109primary photons processed in 3.6 s computing time.Significance. We reported the first dosimetric study of the VCT-BREAST platform, a fast MC simulation tool for real-time virtual dosimetry and imaging trials in BCT, investigating the dose delivery performance of three clinical BCT scanners. This tool can be adopted to investigate also the effects on the 3D dose distribution produced by changes in the geometrical and spectrum characteristics of a cone-beam BCT scanner.

Simulation of X-ray projections on GPU: Benchmarking gVirtualXray with clinically realistic phantoms

BACKGROUND AND OBJECTIVES

This study provides a quantitative comparison of images created using gVirtualXray (gVXR) to both Monte Carlo (MC) and real images of clinically realistic phantoms. gVirtualXray is an open-source framework that relies on the Beer-Lambert law to simulate X-ray images in realtime on a graphics processor unit (GPU) using triangular meshes.

METHODS

Images are generated with gVirtualXray and compared with a corresponding ground truth image of an anthropomorphic phantom: (i) an X-ray projection generated using a Monte Carlo simulation code, (ii) real digitally reconstructed radiographs (DRRs), (iii) computed tomography (CT) slices, and (iv) a real radiograph acquired with a clinical X-ray imaging system. When real images are involved, the simulations are used in an image registration framework so that the two images are aligned.

RESULTS

The mean absolute percentage error (MAPE) between the images simulated with gVirtualXray and MC is 3.12%, the zero-mean normalised cross-correlation (ZNCC) is 99.96% and the structural similarity index (SSIM) is 0.99. The run-time is 10 days for MC and 23 ms with gVirtualXray. Images simulated using surface models segmented from a CT scan of the Lungman chest phantom were similar to (i) DRRs computed from the CT volume and (ii) an actual digital radiograph. CT slices reconstructed from images simulated with gVirtualXray were comparable to the corresponding slices of the original CT volume.

CONCLUSIONS

When scattering can be ignored, accurate images that would take days using MC can be generated in milliseconds with gVirtualXray. This speed of execution enables the use of repetitive simulations with varying parameters, e.g. to generate training data for a deep-learning algorithm, and to minimise the objective function of an optimisation problem in image registration. The use of surface models enables the combination of X-ray simulation with real-time soft-tissue deformation and character animation, which can be deployed in virtual reality applications.

Air Kerma Calculation in Diagnostic Medical Imaging Devices Using Group Method of Data Handling Network

The air kerma, which is the amount of energy given off by a radioactive substance, is essential for medical specialists who use radiation to diagnose cancer problems. The amount of energy that a photon has when it hits something can be described as the air kerma (the amount of energy that was deposited in the air when the photon passed through it). Radiation beam intensity is represented by this value. Hospital X-ray equipment has to account for the heel effect, which means that the borders of the picture obtain a lesser radiation dosage than the center, and that air kerma is not symmetrical. The voltage of the X-ray machine can also affect the uniformity of the radiation. This work presents a model-based approach to predict air kerma at various locations inside the radiation field of medical imaging instruments, making use of just a small number of measurements. Group Method of Data Handling (GMDH) neural networks are suggested for this purpose. Firstly, a medical X-ray tube was modeled using Monte Carlo N Particle (MCNP) code simulation algorithm. X-ray tubes and detectors make up medical X-ray CT imaging systems. An X-ray tube's electron filament, thin wire, and metal target produce a picture of the electrons' target. A small rectangular electron source modeled electron filaments. An electron source target was a thin, 19,290 kg/m³ tungsten cube in a tubular hoover chamber. The electron source-object axis of the simulation object is 20° from the vertical. For most medical X-ray imaging applications, the kerma of the air was calculated at a variety of discrete locations within the conical X-ray beam, providing an accurate data set for network training. Various locations were taken into account in the aforementioned voltages inside the radiation field as the input of the GMDH network. For diagnostic radiology applications, the trained GMDH model could determine the air kerma at any location in the X-ray field of view and for a wide range of X-ray tube voltages with a Mean Relative Error (MRE) of less than 0.25%. This study yielded the following results: (1) The heel effect is included when calculating air kerma. (2) Computing the air kerma using an artificial neural network trained with minimal data. (3) An artificial neural network quickly and reliably calculated air kerma. (4) Figuring out the air kerma for the operating voltage of medical tubes. The high accuracy of the trained neural network in determining air kerma guarantees the usability of the presented method in operational conditions.