A novel hardware duo of beam modulation and shielding to reduce scatter acquisition and dose in cone-beam breast CT

PURPOSE

In cone-beam breast CT, scattered photons form a large portion of the acquired signal, adversely impacting image quality throughout the frequency response of the imaging system. Prior simulation studies provided proof of concept for utilization of a hardware solution to prevent scatter acquisition. Here, we report the design, implementation, and characterization of an auxiliary apparatus of fluence modulation and scatter shielding that does indeed lead to projections with a reduced level of scatter.

METHODS

An apparatus was designed for permanent installation within an existing cone-beam CT system. The apparatus is composed of two primary assemblies: a "Fluence Modulator" (FM) and a "Scatter Shield" (SS). The design of the assemblies enables them to operate in synchrony during image acquisition, converting the sourced x-rays into a moving narrow beam. During a projection, this narrow beam sweeps the entire fan angle coverage of the imaging system. As the two assemblies are contingent on one another, their joint implementation is described in the singular as apparatus FM-SS. The FM and the SS assemblies are each comprised a metal housing, a sensory system, and a robotic system. A controller unit handles their relative movements. A series of comparative studies were conducted to evaluate the performance of a cone-beam CT system in two "modes" of operation: with and without FM-SS installed, and to compare the results of physical implementation with those previously simulated. The dynamic range requirements of the utilized detector in the cone-beam CT imaging system were first characterized, independent of the mode of operation. We then characterized and compared the spatial resolution of the imaging system with, and without, FM-SS. A physical breast phantom, representative of an average size breast, was developed and imaged. Actual differences in signal level obtained with, versus without, FM-SS were then compared to the expected level gains based on previously reported simulations. Following these initial assessments, the scatter acquisition in each projection in both modes of operation was investigated. Finally, as an initial study of the impact of FM-SS on radiation dose in an average size breast, a series of Monte Carlo simulations were coupled with physical measurements of air kerma, with and without FM-SS.

RESULTS

With implementation of FM-SS, the detector's required dynamic range was reduced by a factor of 5.5. Substantial reduction in the acquisition of the scattered rays, by a factor of 5.1 was achieved. With the implementation of FM-SS, deposited dose was reduced by 27% in the studied breast.

CONCLUSIONS

The disclosed implementation of FM-SS, within a cone-beam breast CT system, results in reduction of scatter-components in acquired projections, reduction of dose deposit to the breast, and relaxation of requirements for the detector's dynamic range. Controlling or correcting for patient motion occurring during image acquisition remains an open problem to be solved prior to practical clinical usage of FM-SS cone-beam breast CT.

Shading artifact correction in breast CT using an interleaved deep learning segmentation and maximum-likelihood polynomial fitting approach

PURPOSE

The purpose of this work was twofold: (a) To provide a robust and accurate method for image segmentation of dedicated breast CT (bCT) volume data sets, and (b) to improve Hounsfield unit (HU) accuracy in bCT by means of a postprocessing method that uses the segmented images to correct for the low-frequency shading artifacts in reconstructed images.

METHODS

A sequential and iterative application of image segmentation and low-order polynomial fitting to bCT volume data sets was used in the interleaved correction (IC) method. Image segmentation was performed through a deep convolutional neural network (CNN) with a modified U-Net architecture. A total of 45 621 coronal bCT images from 111 patient volume data sets were segmented (using a previously published segmentation algorithm) and used for neural network training, validation, and testing. All patient data sets were selected from scans performed on four different prototype breast CT systems. The adipose voxels for each patient volume data set, segmented using the proposed CNN, were then fit to a three-dimensional low-order polynomial. The polynomial fit was subsequently used to correct for the shading artifacts introduced by scatter and beam hardening in a method termed "flat fielding." An interleaved utilization of image segmentation and flat fielding was repeated until a convergence criterion was satisfied. Mathematical and physical phantom studies were conducted to evaluate the dependence of the proposed algorithm on breast size and the distribution of fibroglandular tissue. In addition, a subset of patient scans (not used in the CNN training, testing or validation) were used to investigate the accuracy of the IC method across different scanner designs and beam qualities.

RESULTS

The IC method resulted in an accurate classification of different tissue types with an average Dice similarity coefficient > 95%, precision > 97%, recall > 95%, and F1-score > 96% across all tissue types. The flat fielding correction of bCT images resulted in a significant reduction in either cupping or capping artifacts in both mathematical and physical phantom studies as measured by the integral nonuniformity metric with an average reduction of 71% for cupping and 30% for capping across different phantom sizes, and the Uniformity Index with an average reduction of 53% for cupping and 34% for capping.

CONCLUSION

The validation studies demonstrated that the IC method improves Hounsfield Units (HU) accuracy and effectively corrects for shading artifacts caused by scatter contamination and beam hardening. The postprocessing approach described herein is relevant to the broad scope of bCT devices and does not require any modification in hardware or existing scan protocols. The trained CNN parameters and network architecture are available for interested users.

Creation and characterization of an ultrasound and CT phantom for noninvasive ultrasound thermometry calibration

Ultrasound thermometry provides noninvasive 2-D temperature monitoring, and in this paper, we have investigated the use of computed tomography (CT) radiodensity to characterize tissues to improve the accuracy of ultrasound thermometry. Agarose-based tissue-mimicking phantoms were created with glyceryl trioleate (a fat-mimicking material) concentration of 0%, 10%, 20%, 30%, 40%, and 50%. The speed of sound (SOS) of the phantoms was measured over a temperature range of 22.1-41.1 °C. CT images of the phantoms were acquired by a clinical dedicated breast CT scanner, followed by calculation of the Hounsfield units (HU). The phantom was heated with a therapeutic acoustic pulse (1.54 MHz), while RF data were acquired with a 10-MHz linear-array transducer. Two-dimensional speckle tracking was used to calculate the thermal strain offline. The tissue-dependent thermal strain parameter required for ultrasound thermometry was analyzed and correlated with CT radiodensity, followed by the validation of the temperature prediction. Results showed that the change in SOS with the temperature increase was opposite in sign between the 0%-10% and 20%-50% trioleate phantoms. The inverse of the tissue-dependent thermal strain parameter of the phantoms was correlated with the CT radiodensity (R(2) = 0.99). A blinded ultrasound thermometry study on phantoms with a trioleate range of 5%-35% demonstrated the capability to estimate the tissue-dependent thermal strain parameter and estimate temperature with error less than ~1 °C. In conclusion, CT radiodensity may provide a method for improving ultrasound thermometry in heterogeneous tissues.