Monochromatic breast computed tomography with synchrotron radiation: phase-contrast and phase-retrieved image comparison and full-volume reconstruction

Optimization of a customized simultaneous algebraic reconstruction technique algorithm for phase-contrast breast computed tomography

Objective. To introduce the optimization of a customized GPU-based simultaneous algebraic reconstruction technique (cSART) in the field of phase-contrast breast computed tomography (bCT). The presented algorithm features a 3D bilateral regularization filter that can be tuned to yield optimal performance for clinical image visualization and tissues segmentation. Approach. Acquisitions of a dedicated test object and a breast specimen were performed at Elettra, the Italian synchrotron radiation (SR) facility (Trieste, Italy) using a large area CdTe single-photon counting detector. Tomographic images were obtained at 5 mGy of mean glandular dose, with a 32 keV monochromatic x-ray beam in the free-space propagation mode. Three independent algorithms parameters were optimized by using contrast-to-noise ratio (CNR), spatial resolution, and noise texture metrics. The results obtained with the cSART algorithm were compared with conventional SART and filtered back projection (FBP) reconstructions. Image segmentation was performed both with gray scale-based and supervised machine-learning approaches. Main results. Compared to conventional FBP reconstructions, results indicate that the proposed algorithm can yield images with a higher CNR (by 35% or more), retaining a high spatial resolution while preserving their textural properties. Alternatively, at the cost of an increased image ‘patchiness’, the cSART can be tuned to achieve a high-quality tissue segmentation, suggesting the possibility of performing an accurate glandularity estimation potentially of use in the realization of realistic 3D breast models starting from low radiation dose images. Significance. The study indicates that dedicated iterative reconstruction techniques could provide significant advantages in phase-contrast bCT imaging. The proposed algorithm offers great flexibility in terms of image reconstruction optimization, either toward diagnostic evaluation or image segmentation.

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.

Motion artifacts assessment and correction using optical tracking in synchrotron radiation breast CT

Purpose

The SYRMA‐3D collaboration is setting up a breast computed tomography (bCT) clinical program at the Elettra synchrotron radiation facility in Trieste, Italy. Unlike the few dedicated scanners available at hospitals, synchrotron radiation bCT requires the patient's rotation, which in turn implies a long scan duration (from tens of seconds to few minutes). At the same time, it allows the achievement of high spatial resolution. These features make synchrotron radiation bCT prone to motion artifacts. This article aims at assessing and compensating for motion artifacts through an optical tracking approach.

Methods

In this study, patients’ movements due to breathing have been first assessed on seven volunteers and then simulated during the CT scans of a breast phantom and a surgical specimen, by adding a periodic oscillatory motion (constant speed, 1 mm amplitude, 12 cycles/minute). CT scans were carried out at 28 keV with a mean glandular dose of 5 mGy. Motion artifacts were evaluated and a correction algorithm based on the optical tracking of fiducial marks was introduced. A quantitative analysis based on the structural similarity (SSIM) index and the normalized mean square error (nMSE) was performed on the reconstructed CT images.

Results

CT images reconstructed through the optical tracking procedure were found to be as good as the motionless reference image. Moreover, the analysis of SSIM and nMSE demonstrated that an uncorrected motion of the order of the system's point spread function (around 0.1 mm in the present case) can be tolerated.

Conclusions

Results suggest that a motion correction procedure based on an optical tracking system would be beneficial in synchrotron radiation bCT.

Quantitative X-ray phase contrast computed tomography with grating interferometry : Biomedical applications of quantitative X-ray grating-based phase contrast computed tomography

The ability of biomedical imaging data to be of quantitative nature is getting increasingly important with the ongoing developments in data science. In contrast to conventional attenuation-based X-ray imaging, grating-based phase contrast computed tomography (GBPC-CT) is a phase contrast micro-CT imaging technique that can provide high soft tissue contrast at high spatial resolution. While there is a variety of different phase contrast imaging techniques, GBPC-CT can be applied with laboratory X-ray sources and enables quantitative determination of electron density and effective atomic number. In this review article, we present quantitative GBPC-CT with the focus on biomedical applications.

Development and preclinical evaluation of a patient-specific high energy x-ray phase sensitive breast tomosynthesis system

BACKGROUND

This article reports the first x-ray phase sensitive breast tomosynthesis (PBT) system that is aimed for direct translation to clinical practice for the diagnosis of breast cancer.

PURPOSE

To report the preclinical evaluation and comparison of the newly built PBT system with a conventional digital breast tomosynthesis (DBT) system.

METHODS AND MATERIALS

The PBT system is developed based on a comprehensive inline phase contrast theoretical model. The system consists of a polyenergetic microfocus x-ray source and a flat panel detector mounted on an arm that is attached to a rotating gantry. It acquires nine projections over a 15° angular span in a stop-and-shoot manner. A dedicated phase retrieval algorithm is integrated with a filtered back-projection method that reconstructs tomographic slices. The American College of Radiology (ACR) accreditation phantom, a contrast detail (CD) phantom and mastectomy tissue samples were imaged at the same glandular dose levels by both the PBT and a standard of care DBT system for image quality characterizations and comparisons.

RESULTS

The PBT imaging scores with the ACR phantom are in good to excellent range and meet the quality assurance criteria set by the Mammography Quality Standard Act. The CD phantom image comparison and associated statistical analyses from two-alternative forced-choice reader studies confirm the improvement offered by the PBT system in terms of contrast resolution, spatial resolution, and conspicuity. The artifact spread function (ASF) analyses revealed a sizable lateral spread of metal artifacts in PBT slices as compared to DBT slices. Signal-to-noise ratio values for various inserts of the ACR and CD phantoms further validated the superiority of the PBT system. Mastectomy sample images acquired by the PBT system showed a superior depiction of microcalcifications vs the DBT system.

CONCLUSION

The PBT imaging technology can be clinically employed for improving the accuracy of breast cancer screening and diagnosis.

Recent advances in X-ray imaging of breast tissue: From two- to three-dimensional imaging

Breast cancer is a globally widespread disease whose detection has already been significantly improved by the introduction of screening programs. Nevertheless, mammography suffers from low soft tissue contrast and the superposition of diagnostically relevant anatomical structures as well as from low values for sensitivity and specificity especially for dense breast tissue. In recent years, two techniques for X-ray breast imaging have been developed that bring advances for the early detection of breast cancer. Grating-based phase-contrast mammography is a new imaging technique that is able to provide three image modalities simultaneously (absorption-contrast, phase-contrast and dark-field signal). Thus, an enhanced detection and delineation of cancerous structures in the phase-contrast image and an improved visualization and characterization of microcalcifications in the dark-field image is possible. Furthermore, latest studies about this approach show that dose-compatible imaging with polychromatic X-ray sources is feasible. In order to additionally overcome the limitations of projection-based imaging, efforts were also made towards the development of breast computed tomography (BCT), which recently led to the first clinical installation of an absorption-based BCT system. Further research combining the benefits of both imaging technologies is currently in progress. This review article summarizes the latest advances in phase-contrast imaging for the female breast (projection-based and three-dimensional view) with special focus on possible clinical implementations in the future.

Experimental optimization of the energy for breast-CT with synchrotron radiation

Breast Computed Tomography (bCT) is a three-dimensional imaging technique that is raising interest among radiologists as a viable alternative to mammographic planar imaging. In X-rays imaging it would be desirable to maximize the capability of discriminating different tissues, described by the Contrast to Noise Ratio (CNR), while minimizing the dose (i.e. the radiological risk). Both dose and CNR are functions of the X-ray energy. This work aims at experimentally investigating the optimal energy that, at fixed dose, maximizes the CNR between glandular and adipose tissues. Acquisitions of both tissue-equivalent phantoms and actual breast specimens have been performed with the bCT system implemented within the Syrma-3D collaboration at the Syrmep beamline of the Elettra synchrotron (Trieste). The experimental data have been also compared with analytical simulations and the results are in agreement. The CNR is maximized at energies around 26–28 keV. These results are in line with the outcomes of a previously presented simulation study which determined an optimal energy of 28 keV for a large set of breast phantoms with different diameters and glandular fractions. Finally, a study on photon starvation has been carried out to investigate how far the dose can be reduced still having suitable images for diagnostics.

Radiochromic film dosimetry in synchrotron radiation breast computed tomography: a phantom study

This study relates to the INFN project SYRMA-3D for in vivo phase-contrast breast computed tomography using the SYRMEP synchrotron radiation beamline at the ELETTRA facility in Trieste, Italy. This peculiar imaging technique uses a novel dosimetric approach with respect to the standard clinical procedure. In this study, optimization of the acquisition procedure was evaluated in terms of dose delivered to the breast. An offline dose monitoring method was also investigated using radiochromic film dosimetry. Various irradiation geometries have been investigated for scanning the prone patient's pendant breast, simulated by a 14 cm-diameter polymethylmethacrylate cylindrical phantom containing pieces of calibrated radiochromic film type XR-QA2. Films were inserted mid-plane in the phantom, as well as wrapped around its external surface, and irradiated at 38 keV, with an air kerma value that would produce an estimated mean glandular dose of 5 mGy for a 14 cm-diameter 50% glandular breast. Axial scans were performed over a full rotation or over 180°. The results point out that a scheme adopting a stepped rotation irradiation represents the best geometry to optimize the dose distribution to the breast. The feasibility of using a piece of calibrated radiochromic film wrapped around a suitable holder around the breast to monitor the scan dose offline is demonstrated.

Image quality comparison between a phase-contrast synchrotron radiation breast CT and a clinical breast CT: a phantom based study

In this study we compared the image quality of a synchrotron radiation (SR) breast computed tomography (BCT) system with a clinical BCT in terms of contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), noise power spectrum (NPS), spatial resolution and detail visibility. A breast phantom consisting of several slabs of breast-adipose equivalent material with different embedded targets (i.e., masses, fibers and calcifications) was used. Phantom images were acquired using a dedicated BCT system installed at the Radboud University Medical Center (Nijmegen, The Netherlands) and the SR BCT system at the SYRMEP beamline of Elettra SR facility (Trieste, Italy) based on a photon-counting detector. Images with the SR setup were acquired mimicking the clinical BCT conditions (i.e., energy of 30 keV and radiation dose of 6.5 mGy). Images were reconstructed with an isotropic cubic voxel of 273 µm for the clinical BCT, while for the SR setup two phase-retrieval (PhR) kernels (referred to as “smooth” and “sharp”) were alternatively applied to each projection before tomographic reconstruction, with voxel size of 57 × 57 × 50 µm³. The CNR for the clinical BCT system can be up to 2-times higher than SR system, while the SNR can be 3-times lower than SR system, when the smooth PhR is used. The peak frequency of the NPS for the SR BCT is 2 to 4-times higher (0.9 mm⁻¹ and 1.4 mm⁻¹ with smooth and sharp PhR, respectively) than the clinical BCT (0.4 mm⁻¹). The spatial resolution (MTF_10%) was estimated to be 1.3 lp/mm for the clinical BCT, and 5.0 lp/mm and 6.7 lp/mm for the SR BCT with the smooth and sharp PhR, respectively. The smallest fiber visible in the SR BCT has a diameter of 0.15 mm, while for the clinical BCT is 0.41 mm. Calcification clusters with diameter of 0.13 mm are visible in the SR BCT, while the smallest diameter for the clinical BCT is 0.29 mm. As expected, the image quality of the SR BCT outperforms the clinical BCT system, providing images with higher spatial resolution and SNR, and with finer granularity. Nevertheless, this study assesses the image quality gap quantitatively, giving indications on the benefits associated with SR BCT and providing a benchmarking basis for its clinical implementation. In addition, SR-based studies can provide a gold-standard in terms of achievable image quality, constituting an upper-limit to the potential clinical development of a given technique.

Optimization of the energy for Breast monochromatic absorption X-ray Computed Tomography

The limits of mammography have led to an increasing interest on possible alternatives such as the breast Computed Tomography (bCT). The common goal of all X-ray imaging techniques is to achieve the optimal contrast resolution, measured through the Contrast to Noise Ratio (CNR), while minimizing the radiological risks, quantified by the dose. Both dose and CNR depend on the energy and the intensity of the X-rays employed for the specific imaging technique. Some attempts to determine an optimal energy for bCT have suggested the range 22 keV-34 keV, some others instead suggested the range 50 keV-60 keV depending on the parameters considered in the study. Recent experimental works, based on the use of monochromatic radiation and breast specimens, show that energies around 32 keV give better image quality respect to setups based on higher energies. In this paper we report a systematic study aiming at defining the range of energies that maximizes the CNR at fixed dose in bCT. The study evaluates several compositions and diameters of the breast and includes various reconstruction algorithms as well as different dose levels. The results show that a good compromise between CNR and dose is obtained using energies around 28 keV.

Quantitative characterization of breast tissues with dedicated CT imaging

A quantitative characterization of the soft tissues composing the human breast is achieved by means of a monochromatic CT phase-contrast imaging system, through accurate measurements of their attenuation coefficients within the energy range of interest for breast CT clinical examinations. Quantitative measurements of linear attenuation coefficients are performed on tomographic reconstructions of surgical samples, using monochromatic x-ray beams from a synchrotron source and a free space propagation setup. An online calibration is performed on the obtained reconstructions, in order to reassess the validity of the standard calibration procedure of the CT scanner. Three types of healthy tissues (adipose, glandular, and skin) and malignant tumors, when present, are considered from each sample. The measured attenuation coefficients are in very good agreement with the outcomes of similar studies available in the literature, although they span an energy range that was mostly neglected in the previous studies. No globally significant differences are observed between healthy and malignant dense tissues, although the number of considered samples does not appear sufficient to address the issue of a quantitative differentiation of tumors. The study assesses the viability of the proposed methodology for the measurement of linear attenuation coefficients, and provides a denser sampling of attenuation data in the energy range useful to breast CT.

Post-reconstruction 3D single-distance phase retrieval for multi-stage phase-contrast tomography with photon-counting detectors

In the case of single-distance propagation-based phase-contrast X-ray computed tomography with synchrotron radiation, the conventional reconstruction pipeline includes an independent 2D phase retrieval filtering of each acquired projection prior to the actual reconstruction. In order to compensate for the limited height of the X-ray beam or the small sensitive area of most modern X-ray photon-counting detectors, it is quite common to image large objects with a multi-stage approach, i.e. several acquisitions at different vertical positions of the sample. In this context, the conventional reconstruction pipeline may introduce artifacts at the margins of each vertical stage. This article presents a modified computational protocol where a post-reconstruction 3D volume phase retrieval is applied. By comparing the conventional 2D and the proposed 3D reconstructions of a large mastectomy specimen (9 cm in diameter and 3 cm in height), it is here shown that the 3D approach compensates for the multi-stage artifacts, it avoids refined projection stitching, and the image quality in terms of spatial resolution, contrast and contrast-to-noise ratio is preserved.