Development and Initial Demonstration of a Low-Dose Dedicated Fully 3D Breast CT System

Reduction of scatter in breast CT yields improved microcalcification visibility

The inadequate visibility of microcalcifications-small calcium deposits that cue radiologists to early stages of cancer-is a major limitation in current designs of dedicated breast computed tomography (bCT). This limitation has previously been attributed to the constituent components, spatial resolution, and utilized dose. Scattered radiation has been considered an occurrence with low-frequency impacts that can be compensated for in post-processing. We hypothesized, however, that the acquisition of scattered radiation has a far more detrimental impact on clinically relevant features than has previously been understood. Critically, acquisition of scatter leads to the reduced visibility of microcalcifications. This hypothesis was investigated and supported via mathematical derivations and simulation studies. We conducted a series of comparative studies in which four bCT systems were simulated under iso-dose and iso-resolution conditions, characterizing the dependencies of microcalcification contrast on accumulated scatter. Included among the simulated systems is a novel bCT design-narrow beam bCT (NB-bCT)-that captures nearly zero scatter. We find that current bCT systems suffer from significant levels of scatter. As validated in theory, depending on the system and size of microcalcifications, between 25% and over 70% of contrast resolution is lost due to scatter. The results in NB-bCT, however, provide evidence that by removing scatter build-up in projections, the contrast of microcalcifications in a bCT image is preserved, regardless of their size or location in the breast.

Effects of kV, filtration, dose, and object size on soft tissue and iodine contrast in dedicated breast CT

PURPOSE

Clinical use of dedicated breast computed tomography (bCT) requires relatively short scan times necessitating systems with high frame rates. This in turn impacts the x-ray tube operating range. We characterize the effects of tube voltage, beam filtration, dose, and object size on contrast and noise properties related to soft tissue and iodine contrast agents as a way to optimize imaging protocols for soft tissue and iodine contrast at high frame rates.

METHODS

This study design uses the signal-difference-to-noise ratio (SDNR), noise-equivalent quanta (NEQ), and detectability (d´) as measures of imaging performance for a prototype breast CT scanner that utilizes a pulsed x-ray tube (with a 4 ms pulse width) at 43.5 fps acquisition rate. We assess a range of kV, filtration, breast phantom size, and mean glandular dose (MGD). Performance measures are estimated from images of adipose-equivalent breast phantoms machined to have a representative size and shape of small, medium, and large breasts. Water (glandular tissue equivalent) and iodine contrast (5 mg/ml) were used to fill two cylindrical wells in the phantoms.

RESULTS

Air kerma levels required for obtaining an MGD of 6 mGy ranged from 7.1 to 9.1 mGy and are reported across all kV, filtration, and breast phantom sizes. However, at 50 kV, the thick filters (0.3 mm of Cu or Gd) exceeded the maximum available mA of the x-ray generator, and hence, these conditions were excluded from subsequent analysis. There was a strong positive association between measurements of SDNR and d' (R²  > 0.97) within the range of parameters investigated in this work. A significant decrease in soft tissue SDNR was observed for increasing phantom size and increasing kV with a maximum SDNR at 50 kV with 0.2 mm Cu or 0.2 mm Gd filtration. For iodine contrast SDNR, a significant decrease was observed with increasing phantom size, but a decrease in SDNR for increasing kV was only observed for 70 kV (50 and 60 kV were not significantly different). Thicker Gd filtration (0.3 mm Gd) resulted in a significant increase in iodine SDNR and decrease in soft tissue SDNR but requires significantly more tube current to deliver the same MGD.

CONCLUSIONS

The choice of 60 kV with 0.2 mm Gd filtration provides a good trade-off for maximizing both soft tissue and iodine contrast. This scanning technique takes advantage of the ~50 keV Gd k-edge to produce contrast and can be achieved within operating range of the x-ray generator used in this work. Imaging at 60 kV allows for a greater range in dose delivered to the large breast sizes when uniform image quality is desired across all breast sizes. While imaging performance metrics (i.e., detectability index and SDNR) were shown to be strongly correlated, the methodologies presented in this work for the estimation of NEQ (and subsequently d') provides a meaningful description of the spatial resolution and noise characteristics of this prototype bCT system across a range of beam quality, dose, and object sizes.

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

A program devoted to performing the first in vivo synchrotron radiation (SR) breast computed tomography (BCT) is ongoing at the Elettra facility. Using the high spatial coherence of SR, phase-contrast (PhC) imaging techniques can be used. The latest high-resolution BCT acquisitions of breast specimens, obtained with the propagation-based PhC approach, are herein presented as part of the SYRMA-3D collaboration effort toward the clinical exam. Images are acquired with a 60 - μ m   pixel dead-time-free single-photon-counting CdTe detector. The samples are imaged at 32 and 38 keV in a continuous rotating mode, delivering 5 to 20 mGy of mean glandular dose. Contrast-to-noise ratio (CNR) and spatial resolution performances are evaluated for both PhC and phase-retrieved images, showing that by applying the phase-retrieval algorithm a five-time CNR increase can be obtained with a minor loss in spatial resolution across soft tissue interfaces. It is shown that, despite having a poorer CNR, PhC images can provide a sharper visualization of microcalcifications, thus being complementary to phase-retrieved images. Furthermore, the first full-volume scan of a mastectomy sample ( 9 × 9 × 3    cm 3 ) is reported. This investigation into surgical specimens indicates that SR BCT in terms of CNR, spatial resolution, scan duration, and scan volume is feasible.

Characterization of CT Hounsfield Units for 3D acquisition trajectories on a dedicated breast CT system

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is important, especially when using multiple sources, novel detector and non-traditional trajectories. Dedicated clinical breast CT (BCT) systems therefore should be similarly evaluated. In this study, uniform cylindrical phantoms filled with various uniform density fluids were used to characterize differences in HU values between simple circular and complex 3D (saddle) orbits. Based on ACR recommendations, the HU accuracy, center-to-edge variability within a slice, and overall variability within the reconstructed volume were characterized for simple and complex acquisitions possible on a single versatile BCT system. Results illustrate the statistically significantly better performance of the saddle orbit, especially close to the chest and nipple regions of what would clinically be a pendant breast volume. The incomplete cone beam acquisition of a simple circular orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values throughout the reconstructed volumes. This study further establishes the advantages of using 3D acquisition trajectories for breast CT as well as other applications by demonstrating the robustness of HU values throughout large reconstructed volumes.

Evaluation of the BreastSimulator software platform for breast tomography

The aim of this work was the evaluation of the software BreastSimulator, a breast x-ray imaging simulation software, as a tool for the creation of 3D uncompressed breast digital models and for the simulation and the optimization of computed tomography (CT) scanners dedicated to the breast. Eight 3D digital breast phantoms were created with glandular fractions in the range 10%-35%. The models are characterised by different sizes and modelled realistic anatomical features. X-ray CT projections were simulated for a dedicated cone-beam CT scanner and reconstructed with the FDK algorithm. X-ray projection images were simulated for 5 mono-energetic (27, 32, 35, 43 and 51 keV) and 3 poly-energetic x-ray spectra typically employed in current CT scanners dedicated to the breast (49, 60, or 80 kVp). Clinical CT images acquired from two different clinical breast CT scanners were used for comparison purposes. The quantitative evaluation included calculation of the power-law exponent, β, from simulated and real breast tomograms, based on the power spectrum fitted with a function of the spatial frequency, f, of the form S(f)  =  α/f   ^β . The breast models were validated by comparison against clinical breast CT and published data. We found that the calculated β coefficients were close to that of clinical CT data from a dedicated breast CT scanner and reported data in the literature. In evaluating the software package BreastSimulator to generate breast models suitable for use with breast CT imaging, we found that the breast phantoms produced with the software tool can reproduce the anatomical structure of real breasts, as evaluated by calculating the β exponent from the power spectral analysis of simulated images. As such, this research tool might contribute considerably to the further development, testing and optimisation of breast CT imaging techniques.

Characterization of X-ray scattering for various phantoms and clinical breast geometries using breast CT on a dedicated hybrid system

OBJECTIVE

The purpose of this study was to utilize a dedicated breast CT system using a 2D beam stop array to physically evaluate the scatter to primary ratios (SPRs) of different geometric phantoms and prospectively acquired clinical patient data.

METHODS

Including clinically unrealizable compositions of 100% glandular and 100% fat, projection images were acquired using three geometrically different phantoms filled with fluids simulating breast tissue. The beam stop array method was used for measuring scatter in projection space, and creating the scatter corrected primary images. 2D SPRs were calculated. Additionally, a new figure of merit, the 3D normalized scatter contribution (NSC) volumes were calculated.

RESULTS

The 2D SPR values (0.52-1.10) were primarily dependent on phantom geometry; a secondary dependence was due to their uniform density; 2D SPRs were low frequency and smoothly varying in the uniformly filled phantoms. SPRs of clinical patient data followed similar trends as phantoms, but with noticeable deviations and high frequency components due to the heterogeneous distribution of glandular tissue. The maximum measured patient 2D SPRs were all <0.6, even for the largest diameter breast. These results demonstrate modest scatter components with changing object geometries and densities; the 3D NSC volumes with higher frequency components help visualize scatter distribution throughout the reconstructed image volumes. Furthermore, the SPRs in the heterogeneous clinical breast cases were underestimated by the equivalent density, uniformly filled phantoms.

CONCLUSIONS

These results provide guidance on the use of uniformly distributed density and differently shaped phantoms when considering simulations. They also clearly demonstrate that results from patients can vary considerably from 2D SPRs of uniformly simulated phantoms.

Three dimensional dose distribution comparison of simple and complex acquisition trajectories in dedicated breast CT

PURPOSE

A novel breast CT system capable of arbitrary 3D trajectories has been developed to address cone beam sampling insufficiency as well as to image further into the patient's chest wall. The purpose of this study was to characterize any trajectory-related differences in 3D x-ray dose distribution in a pendant target when imaged with different orbits.

METHODS

Two acquisition trajectories were evaluated: circular azimuthal (no-tilt) and sinusoidal (saddle) orbit with ±15° tilts around a pendant breast, using Monte Carlo simulations as well as physical measurements. Simulations were performed with tungsten (W) filtration of a W-anode source; the simulated source flux was normalized to the measured exposure of a W-anode source. A water-filled cylindrical phantom was divided into 1 cm(3) voxels, and the cumulative energy deposited was tracked in each voxel. Energy deposited per voxel was converted to dose, yielding the 3D distributed dose volumes. Additionally, three cylindrical phantoms of different diameters (10, 12.5, and 15 cm) and an anthropomorphic breast phantom, initially filled with water (mimicking pure fibroglandular tissue) and then with a 75% methanol-25% water mixture (mimicking 50-50 fibroglandular-adipose tissues), were used to simulate the pendant breast geometry and scanned on the physical system. Ionization chamber calibrated radiochromic film was used to determine the dose delivered in a 2D plane through the center of the volume for a fully 3D CT scan using the different orbits.

RESULTS

Measured experimental results for the same exposure indicated that the mean dose measured throughout the central slice for different diameters ranged from 3.93 to 5.28 mGy, with the lowest average dose measured on the largest cylinder with water mimicking a homogeneously fibroglandular breast. These results align well with the cylinder phantom Monte Carlo studies which also showed a marginal difference in dose delivered by a saddle trajectory in the central slice. Regardless of phantom material or filled fluid density, dose delivered by the saddle scan was negligibly different than the simple circular, no-tilt scans. The average dose measured in the breast phantom was marginally higher for saddle than the circular no tilt scan at 3.82 and 3.87 mGy, respectively.

CONCLUSIONS

Not only does nontraditional 3D-trajectory CT scanning yield more complete sampling of the breast volume but also has comparable dose deposition throughout the breast and anterior chest volume, as verified by Monte Carlo simulation and physical measurements.

X-ray scatter correction method for dedicated breast computed tomography: improvements and initial patient testing

A previously proposed x-ray scatter correction method for dedicated breast computed tomography was further developed and implemented so as to allow for initial patient testing. The method involves the acquisition of a complete second set of breast CT projections covering 360° with a perforated tungsten plate in the path of the x-ray beam. To make patient testing feasible, a wirelessly controlled electronic positioner for the tungsten plate was designed and added to a breast CT system. Other improvements to the algorithm were implemented, including automated exclusion of non-valid primary estimate points and the use of a different approximation method to estimate the full scatter signal. To evaluate the effectiveness of the algorithm, evaluation of the resulting image quality was performed with a breast phantom and with nine patient images. The improvements in the algorithm resulted in the avoidance of introduction of artifacts, especially at the object borders, which was an issue in the previous implementation in some cases. Both contrast, in terms of signal difference and signal difference-to-noise ratio were improved with the proposed method, as opposed to with the correction algorithm incorporated in the system, which does not recover contrast. Patient image evaluation also showed enhanced contrast, better cupping correction, and more consistent voxel values for the different tissues. The algorithm also reduces artifacts present in reconstructions of non-regularly shaped breasts. With the implemented hardware and software improvements, the proposed method can be reliably used during patient breast CT imaging, resulting in improvement of image quality, no introduction of artifacts, and in some cases reduction of artifacts already present. The impact of the algorithm on actual clinical performance for detection, diagnosis and other clinical tasks in breast imaging remains to be evaluated.

Screening for breast cancer with Breast-CT in a ProHTA simulation

Aims: The potential of dedicated Breast-CT is evaluated by simulating its impact onto the performance of the German breast cancer screening program. Attendance rates, cancer detection and economic implications are quantified. Methods: Based on a prospective health technology assessment approach, we simulated screening in different scenarios. Results: In the simulation, attendance rates increase from 54 to up to 72% due to reduced pain. Breast cancers will be detected earlier while nodal positives and distant recurrences decrease. Assuming no additional cost, cost savings of up to €55 million in one screening period are computed. Conclusion: The simulation indicates that earlier cancer detection, fewer unnecessary biopsies and less pain are potential benefits of Breast-CT resulting in cost savings and higher attendance.

Task-based optimization of dedicated breast CT via Hotelling observer metrics

PURPOSE

The purpose of this work is to develop and demonstrate a set of practical metrics for CT systems optimization. These metrics, based on the Hotelling observer (HO) figure of merit, are task-based. The authors therefore take the specific example of optimizing a dedicated breast CT system, including the reconstruction algorithm, for two relevant tasks, signal detection and Rayleigh discrimination.

METHODS

A dedicated breast CT system is simulated using specifications in the literature from an existing prototype. The authors optimize configuration and image reconstruction algorithm parameters for two tasks: the detection of simulated microcalcifications and the discrimination of two adjacent, high-contrast signals, known as the Rayleigh discrimination task. The effects on task performance of breast diameter, signal location, image grid size, projection view number, and reconstruction filter were all investigated. Two HO metrics were evaluated: the percentage of correct decisions in a two-alternative forced choice experiment (equivalent to area under the ROC curve or AUC), and the HO efficiency, defined as the squared ratio of HO signal-to-noise ratio (SNR) in the reconstructed image to HO SNR in the projection data.

RESULTS

The ease and efficiency of the HO metric computation allows a rapid high-resolution survey of many system parameters. Optimization of a range of system parameters using the HO results in images that subjectively appear optimal for the tasks investigated. Further, the results of assessment through the HO reproduce closely many existing results in the literature regarding the impact of parameter selection on image quality.

CONCLUSIONS

This study demonstrates the utility of a task-based approach to system design, evaluation, and optimization. The methodology presented is equally applicable to determining the impact of a wide range of factors, including patient parameters, system and acquisition design, and the reconstruction algorithm. The results demonstrate the versatility of the proposed HO formalism by not only generating a set of parameters that are optimal for a given task but also by qualitatively reproducing many existing results from the breast CT literature. Meanwhile, the implementation of the proposed methodology is straightforward and entirely simulation-based. This is an attractive feature for many system optimization problems, where the goal is to analyze the individual system components such as the image reconstruction algorithm. Final assessment of the system as a whole should be based also on real data studies.