Clinical application of low-dose phase contrast breast CT: methods for the optimization of the reconstruction workflow

High resolution propagation-based lung imaging at clinically relevant X-ray dose levels

Absorption-based clinical computed tomography (CT) is the current imaging method of choice in the diagnosis of lung diseases. Many pulmonary diseases are affecting microscopic structures of the lung, such as terminal bronchi, alveolar spaces, sublobular blood vessels or the pulmonary interstitial tissue. As spatial resolution in CT is limited by the clinically acceptable applied X-ray dose, a comprehensive diagnosis of conditions such as interstitial lung disease, idiopathic pulmonary fibrosis or the characterization of small pulmonary nodules is limited and may require additional validation by invasive lung biopsies. Propagation-based imaging (PBI) is a phase sensitive X-ray imaging technique capable of reaching high spatial resolutions at relatively low applied radiation dose levels. In this publication, we present technical refinements of PBI for the characterization of different artificial lung pathologies, mimicking clinically relevant patterns in ventilated fresh porcine lungs in a human-scale chest phantom. The combination of a very large propagation distance of 10.7 m and a photon counting detector with [Formula: see text] pixel size enabled high resolution PBI CT with significantly improved dose efficiency, measured by thermoluminescence detectors. Image quality was directly compared with state-of-the-art clinical CT. PBI with increased propagation distance was found to provide improved image quality at the same or even lower X-ray dose levels than clinical CT. By combining PBI with iodine k-edge subtraction imaging we further demonstrate that, the high quality of the calculated iodine concentration maps might be a potential tool for the analysis of lung perfusion in great detail. Our results indicate PBI to be of great value for accurate diagnosis of lung disease in patients as it allows to depict pathological lesions non-invasively at high resolution in 3D. This will especially benefit patients at high risk of complications from invasive lung biopsies such as in the setting of suspected idiopathic pulmonary fibrosis (IPF).

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.

Improving spatial resolution with an edge-enhancement model for low-dose propagation-based X-ray phase-contrast computed tomography

Propagation-based X-ray phase-contrast computed tomography (PB-PCCT) has been increasingly popular for distinguishing low contrast tissues. Phase retrieval is an important step to quantitatively obtain the phase information before the tomographic reconstructions, while typical phase retrieval methods in PB-PCCT, such as homogenous transport of intensity equation (TIE-Hom), are essentially low-pass filters and thus improve the signal to noise ratio at the expense of the reduced spatial resolution of the reconstructed image. To improve the reconstructed spatial resolution, measured phase contrast projections with high edge enhancement and the phase projections retrieved by TIE-Hom were weighted summed and fed into an iterative tomographic algorithm within the framework of the adaptive steepest descent projections onto convex sets (ASD-POCS), which was employed for suppressing the image noise in low dose reconstructions because of the sparse-view scanning strategy or low exposure time for single phase contrast projection. The merging strategy decreases the accuracy of the linear model of PB-PCCT and would finally lead to the reconstruction failure in iterative reconstructions. Therefore, the additive median root prior is also introduced in the algorithm to partly increase the model accuracy. The reconstructed spatial resolution and noise performance can be flexibly balanced by a pair of antagonistic hyper-parameters. Validations were performed by the established phase-contrast Feldkamp-Davis-Kress, phase-retrieved Feldkamp-Davis-Kress, conventional ASD-POCS and the proposed enhanced ASD-POCS with a numerical phantom dataset and experimental biomaterial dataset. Simulation results show that the proposed algorithm outperforms the conventional ASD-POCS in spatial evaluation assessments such as root mean square error (a ratio of 9.78%), contrast to noise ratio (CNR) (a ratio of 7.46%), and also frequency evaluation assessments such as modulation transfer function (a ratio of 66.48% of MTF_50% (50% MTF value)), noise power spectrum (a ratio of 35.25% of f_50% (50% value of the Nyquist frequency)) and noise equivalent quanta (1-2 orders of magnitude at high frequencies). Experimental results again confirm the superiority of proposed strategy relative to the conventional one in terms of edge sharpness and CNR (an average increase of 67.35%).

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.

Deep Efficient End-to-end Reconstruction (DEER) Network for Few-view Breast CT Image Reconstruction

Breast CT provides image volumes with isotropic resolution in high contrast, enabling detection of small calcification (down to a few hundred microns in size) and subtle density differences. Since breast is sensitive to x-ray radiation, dose reduction of breast CT is an important topic, and for this purpose, few-view scanning is a main approach. In this article, we propose a Deep Efficient End-to-end Reconstruction (DEER) network for few-view breast CT image reconstruction. The major merits of our network include high dose efficiency, excellent image quality, and low model complexity. By the design, the proposed network can learn the reconstruction process with as few as O ( N ) parameters, where N is the side length of an image to be reconstructed, which represents orders of magnitude improvements relative to the state-of-the-art deep-learning-based reconstruction methods that map raw data to tomographic images directly. Also, validated on a cone-beam breast CT dataset prepared by Koning Corporation on a commercial scanner, our method demonstrates a competitive performance over the state-of-the-art reconstruction networks in terms of image quality. The source code of this paper is available at: https://github.com/HuidongXie/DEER.

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.

Propagation-based x-ray phase-contrast tomography of mastectomy samples using synchrotron radiation

PURPOSE

Propagation-based phase-contrast computed tomography (PB-CT) is a method for three-dimensional x-ray imaging that utilizes refraction, as well as absorption, of x rays in the tissues to increase the signal-to-noise ratio (SNR) in the resultant images, in comparison with equivalent conventional absorption-only x-ray tomography (CT). Importantly, the higher SNR is achieved without sacrificing spatial resolution or increasing the radiation dose delivered to the imaged tissues. The present work has been carried out in the context of the current development of a breast CT imaging facility at the Australian Synchrotron.

METHODS

Seven unfixed complete mastectomy samples with and without breast cancer lesions have been imaged using absorption-only CT and PB-CT techniques under controlled experimental conditions. The radiation doses delivered to the mastectomy samples during the scans were comparable to those approved for mammographic screening. Physical characteristics of the reconstructed images, such as spatial resolution and SNR, have been measured and compared with the results of the radiological quality assessment of the complete absorption CT and PB-CT image stacks.

RESULTS

Despite the presence of some image artefacts, the PB-CT images have outperformed comparable absorption CT images collected at the same radiation dose, in terms of both the measured objective image characteristics and the radiological image scores. The outcomes of these experiments are shown to be consistent with predictions of the theory of PB-CT imaging and previous reported experimental studies of this imaging modality.

CONCLUSIONS

The results presented in this paper demonstrate that PB-CT holds a high potential for improving on the quality and diagnostic value of images obtained using existing medical x-ray technologies, such as mammography and digital breast tomosynthesis (DBT). If implemented at suitable synchrotron imaging facilities, PB-CT can be used to complement existing imaging modalities, leading to more accurate breast cancer diagnosis.

Free propagation phase-contrast breast CT provides higher image quality than cone-beam breast-CT at low radiation doses: a feasibility study on human mastectomies

In this study we demonstrate the first direct comparison between synchrotron x-ray propagation-based CT (PB-CT) and cone-beam breast-CT (CB-CT) on human mastectomy specimens (N = 12) including different benign and malignant lesions. The image quality and diagnostic power of the obtained data sets were compared and judged by two independent expert radiologists. Two cases are presented in detail in this paper including a comparison with the corresponding histological evaluation. Results indicate that with PB-CT it is possible to increase the level of contrast-to-noise ratio (CNR) keeping the same level of dose used for the CB-CT or achieve the same level of CNR reached by CB-CT at a lower level of dose. In other words, PB-CT can achieve a higher diagnostic potential compared to the commercial breast-CT system while also delivering a considerably lower mean glandular dose. Therefore, we believe that PB-CT technique, if translated to a clinical setting, could have a significant impact in improving breast cancer diagnosis.

Advanced 3D Imaging of Uterine Leiomyoma's Morphology by Propagation-based Phase-Contrast Microtomography

Uterine leiomyoma is the most common benign smooth muscle tumor in women pelvis, originating from the myometrium. It is caused by a disorder of fibrosis, with a large production and disruption of extracellular matrix (ECM). Medical treatments are still very limited and no preventative therapies have been developed. We supposed that synchrotron-based phase-contrast microtomography (PhC-microCT) may be an appropriate tool to assess the 3D morphology of uterine leiomyoma, without the use of any contrast agent. We used this technique to perform the imaging and the quantitative morphometric analysis of healthy myometrium and pathologic leiomyomas. The quantitative morphometric analysis of collagen bundles was coupled to the Roschger approach. This method, previously only used to evaluate mineralized bone density distribution, was applied here to study the fibrosis mass density distribution in healthy and pathologic biopsies from two patients. This protocol was shown to be powerful in studying uterine leiomyomas, detecting also small signs of the ECM alteration. This is of paramount importance not only for the follow-up of the present study, i.e. the investigation of different compounds and their possible therapeutic benefits, but also because it offers new methodologic possibilities for future studies of the ECM in soft tissues of different body districts.

Advancements towards the implementation of clinical phase-contrast breast computed tomography at Elettra

Breast computed tomography (BCT) is an emerging application of X-ray tomography in radiological practice. A few clinical prototypes are under evaluation in hospitals and new systems are under development aiming at improving spatial and contrast resolution and reducing delivered dose. At the same time, synchrotron-radiation phase-contrast mammography has been demonstrated to offer substantial advantages when compared with conventional mammography. At Elettra, the Italian synchrotron radiation facility, a clinical program of phase-contrast BCT based on the free-space propagation approach is under development. In this paper, full-volume breast samples imaged with a beam energy of 32 keV delivering a mean glandular dose of 5 mGy are presented. The whole acquisition setup mimics a clinical study in order to evaluate its feasibility in terms of acquisition time and image quality. Acquisitions are performed using a high-resolution CdTe photon-counting detector and the projection data are processed via a phase-retrieval algorithm. Tomographic reconstructions are compared with conventional mammographic images acquired prior to surgery and with histologic examinations. Results indicate that BCT with monochromatic beam and free-space propagation phase-contrast imaging provide relevant three-dimensional insights of breast morphology at clinically acceptable doses and scan times.

Synchrotron radiation biomedical imaging and radiotherapy: from the UK to the Antipodes

Although the general public might think of 'X-rays' as they are applied to imaging (radiography) and for the treatment of disease (radiotherapy), the use of synchrotron radiation (SR) X-ray beams in these areas of science was a minor activity 50 years ago. The largest gains in science from SR were seen to be in those areas where signals were weakest in laboratory instruments, such as X-ray diffraction and spectroscopy. As the qualities of SR X-rays were explored and more areas of science adopted SR-based methods, this situation changed. About 30 years ago, the clinical advantages of using SR X-ray beams for radiography, radiotherapy and clinical diagnostics started to be investigated. In the UK, a multi-disciplinary group, consisting of clinicians, medical physicists and other scientists working mainly with the Synchrotron Radiation Source (SRS) in Cheshire, started to investigate techniques for diagnosis and potentially a cure for certain cancers. This preliminary work influenced the design of new facilities being constructed around the world, in particular the Imaging and Medical Beam Line on the Australian Synchrotron in Melbourne. Two authors moved from the UK to Australia to participate in this exciting venture. The following is a personal view of some of the highlights of the early-year SRS work, following through to the current activities on the new facility in Australia. This article is part of the theme issue 'Fifty years of synchrotron science: achievements and opportunities'.

Toward Improving Breast Cancer Imaging: Radiological Assessment of Propagation-Based Phase-Contrast CT Technology

RATIONALE AND OBJECTIVES

This study employs clinical/radiological evaluation in establishing the optimum imaging conditions for breast cancer imaging using the X-ray propagation-based phase-contrast tomography.

MATERIALS AND METHODS

Two series of experiments were conducted and in total 161 synchrotron-based computed tomography (CT) reconstructions of one breast mastectomy specimen were produced at different imaging conditions. Imaging factors include sample-to-detector distance, X-ray energy, CT reconstruction method, phase retrieval algorithm applied to the CT projection images and maximum intensity projection. Observers including breast radiologists and medical imaging experts compared the quality of the reconstructed images with reference images approximating the conventional (absorption) CT. Various radiological image quality attributes in a visual grading analysis design were used for the radiological assessments.

RESULTS

The results show that the application of the longest achievable sample-to-detector distance (9.31 m), the lowest employed X-ray energy (32 keV), the full phase retrieval, and the maximum intensity projection can significantly improve the radiological quality of the image. Several combinations of imaging variables resulted in images with very high-quality scores.

CONCLUSION

The results of the present study will support future experimental and clinical attempts to further optimize this innovative approach to breast cancer imaging.

Advantages of breast cancer visualization and characterization using synchrotron radiation phase-contrast tomography

The aim of this study was to highlight the advantages that propagation-based phase-contrast computed tomography (PB-CT) with synchrotron radiation can provide in breast cancer diagnostics. For the first time, a fresh and intact mastectomy sample from a 60 year old patient was scanned on the IMBL beamline at the Australian Synchrotron in PB-CT mode and reconstructed. The clinical picture was described and characterized by an experienced breast radiologist, who underlined the advantages of providing diagnosis on a PB-CT volume rather than conventional two-dimensional modalities. Subsequently, the image quality was assessed by 11 breast radiologists and medical imaging experts using a radiological scoring system. The results indicate that, with the radiation dose delivered to the sample being equal, the accuracy of a diagnosis made on PB-CT images is significantly higher than one using conventional techniques.

Dose and diagnostic performance comparison between phase-contrast mammography with synchrotron radiation and digital mammography: a clinical study report

Two dosimetric quantities [mean glandular dose (MGD) and entrance surface air kerma (ESAK)] and the diagnostic performance of phase-contrast mammography with synchrotron radiation (MSR) are compared to conventional digital mammography (DM). Seventy-one patients (age range, 41 to 82 years) underwent MSR after a DM examination if questionable or suspicious breast abnormalities were not clarified by ultrasonography. The MGD and the ESAK delivered in both examinations were evaluated and compared. Two on-site radiologists rated the images in consensus according to the Breast Imaging Reporting and Data System assessment categories, which were then correlated with the final diagnoses by means of statistical generalized linear models (GLMs). Receiver operating characteristic curves were also used to assess the diagnostic performance by comparing the area under the curve (AUC). An important MGD and ESAK reduction was observed in MSR due to the monoenergetic beam. In particular, an average 43% reduction was observed for the MGD and a reduction of more than 50% for the ESAK. GLM showed higher diagnostic accuracy, especially in terms of specificity, for MSR, confirmed by AUC analysis ([Formula: see text]). The study design implied that the population was characterized by a high prevalence of disease and that the radiologists, who read the DM images before referring the patient to MSR, could have been influenced in their assessments. Within these limitations, the use of synchrotron radiation with the phase-contrast technique applied to mammography showed an important dose reduction and a higher diagnostic accuracy compared with DM. These results could further encourage research on the translation of x-ray phase-contrast imaging into the clinics.

Review of quantitative multiscale imaging of breast cancer

Breast cancer is the most common cancer among women worldwide and ranks second in terms of overall cancer deaths. One of the difficulties associated with treating breast cancer is that it is a heterogeneous disease with variations in benign and pathologic tissue composition, which contributes to disease development, progression, and treatment response. Many of these phenotypes are uncharacterized and their presence is difficult to detect, in part due to the sparsity of methods to correlate information between the cellular microscale and the whole-breast macroscale. Quantitative multiscale imaging of the breast is an emerging field concerned with the development of imaging technology that can characterize anatomic, functional, and molecular information across different resolutions and fields of view. It involves a diverse collection of imaging modalities, which touch large sections of the breast imaging research community. Prospective studies have shown promising results, but there are several challenges, ranging from basic physics and engineering to data processing and quantification, that must be met to bring the field to maturity. This paper presents some of the challenges that investigators face, reviews currently used multiscale imaging methods for preclinical imaging, and discusses the potential of these methods for clinical breast imaging.

SYRMEP Tomo Project: a graphical user interface for customizing CT reconstruction workflows

When considering the acquisition of experimental synchrotron radiation (SR) X-ray CT data, the reconstruction workflow cannot be limited to the essential computational steps of flat fielding and filtered back projection (FBP). More refined image processing is often required, usually to compensate artifacts and enhance the quality of the reconstructed images. In principle, it would be desirable to optimize the reconstruction workflow at the facility during the experiment (beamtime). However, several practical factors affect the image reconstruction part of the experiment and users are likely to conclude the beamtime with sub-optimal reconstructed images. Through an example of application, this article presents SYRMEP Tomo Project (STP), an open-source software tool conceived to let users design custom CT reconstruction workflows. STP has been designed for post-beamtime (off-line use) and for a new reconstruction of past archived data at user's home institution where simple computing resources are available. Releases of the software can be downloaded at the Elettra Scientific Computing group GitHub repository https://github.com/ElettraSciComp/STP-Gui.

Towards breast tomography with synchrotron radiation at Elettra: first images

The aim of the SYRMA-CT collaboration is to set-up the first clinical trial of phase-contrast breast CT with synchrotron radiation (SR). In order to combine high image quality and low delivered dose a number of innovative elements are merged: a CdTe single photon counting detector, state-of-the-art CT reconstruction and phase retrieval algorithms. To facilitate an accurate exam optimization, a Monte Carlo model was developed for dose calculation using GEANT4. In this study, high isotropic spatial resolution (120 μm)(3) CT scans of objects with dimensions and attenuation similar to a human breast were acquired, delivering mean glandular doses in the range of those delivered in clinical breast CT (5-25 mGy). Due to the spatial coherence of the SR beam and the long distance between sample and detector, the images contain, not only absorption, but also phase information from the samples. The application of a phase-retrieval procedure increases the contrast-to-noise ratio of the tomographic images, while the contrast remains almost constant. After applying the simultaneous algebraic reconstruction technique to low-dose phase-retrieved data sets (about 5 mGy) with a reduced number of projections, the spatial resolution was found to be equal to filtered back projection utilizing a four fold higher dose, while the contrast-to-noise ratio was reduced by 30%. These first results indicate the feasibility of clinical breast CT with SR.

A feasibility study of X-ray phase-contrast mammographic tomography at the Imaging and Medical beamline of the Australian Synchrotron

Results are presented of a recent experiment at the Imaging and Medical beamline of the Australian Synchrotron intended to contribute to the implementation of low-dose high-sensitivity three-dimensional mammographic phase-contrast imaging, initially at synchrotrons and subsequently in hospitals and medical imaging clinics. The effect of such imaging parameters as X-ray energy, source size, detector resolution, sample-to-detector distance, scanning and data processing strategies in the case of propagation-based phase-contrast computed tomography (CT) have been tested, quantified, evaluated and optimized using a plastic phantom simulating relevant breast-tissue characteristics. Analysis of the data collected using a Hamamatsu CMOS Flat Panel Sensor, with a pixel size of 100 µm, revealed the presence of propagation-based phase contrast and demonstrated significant improvement of the quality of phase-contrast CT imaging compared with conventional (absorption-based) CT, at medically acceptable radiation doses.