In-line phase-contrast breast tomosynthesis: a phantom feasibility study at a synchrotron radiation facility

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.

A voxel-by-voxel method for mixing two filaments during a 3D printing process for soft-tissue replication in an anthropomorphic breast phantom

Objective. In this study, a novel voxel-by-voxel mixing method is presented, according to which two filaments of different material are combined during the three dimensional (3D) printing process.Approach. In our approach, two types of filaments were used for the replication of soft-tissues, a polylactic acid (PLA) filament and a polypropylene (PP) filament. A custom-made software was used, while a series of breast patient CT scan images were directly associated to the 3D printing process. Each phantom´s layer was printed twice, once with the PLA filament and a second time with the PP filament. For each material, the filament extrusion rate was controlled voxel-by-voxel and was based on the Hounsfield units (HU) of the imported CT images. The phantom was scanned at clinical CT, breast tomosynthesis and micro CT facilities, as the major processing was performed on data from the CT. A side by side comparison between patient´s and phantom´s CT slices by means of profile and histogram comparison was accomplished. Further, in case of profile comparison, the Pearson´s coefficients were calculated.Main results. The visual assessment of the distribution of the glandular tissue in the CT slices of the printed breast anatomy showed high degree of radiological similarity to the corresponding patient´s glandular distribution. The profile plots´ comparison showed that the HU of the replicated and original patient soft tissues match adequately. In overall, the Pearson´s coefficients were above 0.91, suggesting a close match of the CT images of the phantom with those of the patient. The overall HU were close in terms of HU ranges. The HU mean, median and standard deviation of the original and the phantom CT slices were -149, -167, ±65 and -121, -130, ±91, respectively.Significance. The results suggest that the proposed methodology is appropriate for manufacturing of anthropomorphic soft tissue phantoms for x-ray imaging and dosimetry purposes, since it may offer an accurate replication of these tissues.

Synchrotron X-ray Radiation (SXR) in Medical Imaging: Current Status and Future Prospects

Synchrotron X-ray radiation (SXR) has been widely studied to explore the structure of matter. Recently, there has been an intense focus on the medical application of SXR in imaging. This review is intended to explore the latest applications of SXR in medical imaging and to shed light on the advantages and drawbacks of this modality. The article highlights the latest developments in other fields that can greatly enhance the capability and applicability of SXR. The potentials of using machine and deep learning (DL)-based methods to generate synthetic images to use in regular clinics along with the use of photon counting X-ray detectors for spectral medical imaging with SXR are also discussed.

BriXS, a new X-ray inverse Compton source for medical applications

MariX is a research infrastructure conceived for multi-disciplinary studies, based on a cutting-edge system of combined electron accelerators at the forefront of the world-wide scenario of X-ray sources. The generation of X-rays over a large photon energy range will be enabled by two unique X-ray sources: a Free Electron Laser and an inverse Compton source, called BriXS (Bright compact X-ray Source). The X-ray beam provided by BriXS is expected to have an average energy tunable in the range 20-180 keV and intensities between 10¹¹ and 10¹³ photon/s within a relative bandwidth ΔE/E=1-10%. These characteristics, together with a very small source size (~20 μm) and a good transverse coherence, will enable a wide range of applications in the bio-medical field. An additional unique feature of BriXS will be the possibility to make a quick switch of the X-ray energy between two values for dual-energy and K-edge subtraction imaging. In this paper, the expected characteristics of BriXS will be presented, with a particular focus on the features of interest to its possible medical applications.

Elastographic Tomosynthesis From X-Ray Strain Imaging of Breast Cancer

Noncancerous breast tissue and cancerous breast tissue have different elastic properties. In particular, cancerous breast tumors are stiff when compared to the noncancerous surrounding tissue. This difference in elasticity can be used as a means for detection through the method of elastographic tomosynthesis by means of physical modulation. This paper deals with a method to visualize elasticity of soft tissues, particularly breast tissues, via x-ray tomosynthesis. X-ray tomosynthesis is now used to visualize breast tissues with better resolution than the conventional single-shot mammography. The advantage of X-ray tomosynthesis over X-ray CT is that fewer projections are needed than CT to perform the reconstruction, thus radiation exposure and cost are both reduced. Two phantoms were used for the testing of this method, a physical phantom and an in silico phantom. The standard root mean square error in the tomosynthesis for the physical phantom was 2.093 and the error in the in silico phantom was negligible. The elastographs were created through the use of displacement and strain graphing. A Gaussian Mixture Model with an expectation–maximization clustering algorithm was applied in three dimensions with an error of 16.667%. The results of this paper have been substantial when using phantom data. There are no equivalent comparisons yet in 3D x-ray elastographic tomosynthesis. Tomosynthesis with and without physical modulation in the 3D elastograph can identify feature groupings used for biopsy. The studies have potential to be applied to human test data used as a guide for biopsy to improve accuracy of diagnosis results. Further research on this topic could prove to yield new techniques for human patient diagnosis purposes.

Impact of a single distance phase retrieval algorithm on spatial resolution in X-ray inline phase sensitive imaging

A single-projection based phase retrieval method based on the phase attenuation duality principle (PAD) was used to compare the spatial resolution of the acquired phase sensitive and PAD processed phase retrieved images. An inline phase sensitive prototype was used to acquire the phase sensitive images. The prototype incorporates a micro-focus x-ray source and a flat panel detector with a 50 μm pixel pitch. A phantom composed of a 2 cm thick 50-50 adipose-glandular mimicking slab sandwiched with a 0.82 cm thick slanted PMMA sharp edge was used. Phase sensitive image of the phantom was acquired at 120 kV, 3.35 mAs with a 16 μm tube focal spot size under a geometric magnification (M) of 2.5. The PAD based method was applied to the acquired phase sensitive image for the retrieval of phase values. With necessary data processing, modulation transfer function (MTF) curves were determined for the estimation and comparison of the spatial resolution. The PAD processed phase retrieved values of the phantom were in good agreement with the theoretically calculated values. Phase sensitive images showed higher spatial resolution at all spatial frequencies compared to the phase retrieved images. It was noted that the high-frequency signal components in the retrieved image were suppressed that resulted in lower MTF values. When compared to the phase sensitive image, the cutoff resolution (10% MTF) for phase retrieved image dropped 32% from 15.6 lp/mm (32μm) to 10.6 lp/mm (47μm). The resolution offered by this phase sensitive prototype is radiographically enough to detect breast cancer.

Detectability comparison of simulated tumors in digital breast tomosynthesis using high-energy X-ray inline phase sensitive and commercial imaging systems

This study compared the detectability of simulated tumors using a high-energy X-ray inline phase sensitive digital breast tomosynthesis (DBT) prototype and a commercial attenuation-based DBT system. Each system imaged a 5-cm thick modular breast phantom with 50-50 adipose-glandular percentage density containing contrast-detail (CD) test objects to simulate different tumor sizes. A commercial DBT system acquired 15 projection views over 15 degrees (15d-15p) was used to acquire the attenuation-based projection views and to reconstruct the conventional DBT slices. Attenuation-based projection views were acquired at 32 kV, 46 mAs with a mean glandular dose (D_g) of 1.6 mGy. For acquiring phase sensitive projection views, the prototype utilized two acquisition geometries: 11 projection views were acquired over 15 degrees (15d-11p), and 17 projection views were acquired over 16 degrees (16d-17p) at 120 kV, 5.27 mAs with 1.51 mGy under the magnification (M) of 2. A phase retrieval algorithm based on the phase-attenuation duality (PAD) was applied to each projection view, and a modified Feldkamp-Davis-Kress (FDK) algorithm was used to reconstruct the phase sensitive DBT slices. Simulated tumor margins were rated as more conspicuous and better visualized for both phase sensitive acquisition geometries versus conventional DBT imaging. The CD curves confirmed the improvement in both contrast and spatial resolutions with the phase sensitive DBT imaging. The superiority of the phase sensitive DBT imaging was further endorsed by higher contrast to noise ratio (CNR) and figure-of-merit (FOM) values. The CNR improvements provided by the phase sensitive DBT prototype were sufficient to offset the noise reduction provided by the attenuation-based DBT imaging.

Digital breast tomosynthesis: Dose and image quality assessment

The aim of this work was to evaluate how different acquisition geometries and reconstruction parameters affect the performance of four digital breast tomosynthesis (DBT) systems (Senographe Essential - GE, Mammomat Inspiration - Siemens, Selenia Dimensions - Hologic and Amulet Innovality - Fujifilm) on the basis of a physical characterization. Average Glandular Dose (AGD) and image quality parameters such as in-plane/in-depth resolution, signal difference to noise ratio (SDNR) and artefact spread function (ASF) were examined. Measured AGD values resulted below EUREF limits for 2D imaging. A large variability was recorded among the investigated systems: the mean dose ratio DBT/2D ranged between 1.1 and 1.9. In-plane resolution was in the range: 2.2mm^-1-3.8mm^-1 in chest wall-nipple direction. A worse resolution was found for all devices in tube travel direction. In-depth resolution improved with increasing scan angle but was also affected by the choice of reconstruction and post-processing algorithms. The highest z-resolution was provided by Siemens (50°, FWHM=2.3mm) followed by GE (25°, FWHM=2.8mm), while the Fujifilm HR showed the lowest one, despite its wide scan angle (40°, FWHM=4.1mm). The ASF was dependent on scan angle: smaller range systems showed wider ASF curves; however a clear relationship was not found between scan angle and ASF, due to the different post processing and reconstruction algorithms. SDNR analysis, performed on Fujifilm system, demonstrated that pixel binning improves detectability for a fixed dose/projection. In conclusion, we provide a performance comparison among four DBT systems under a clinical acquisition mode.