Dedicated breast CT: radiation dose for circle-plus-line trajectory

Dedicated cone-beam breast CT: Data acquisition strategies based on projection angle-dependent normalized glandular dose coefficients

BACKGROUND

Dedicated cone-beam breast computed tomography (CBBCT) using short-scan acquisition is being actively investigated to potentially reduce the radiation dose to the breast. This would require determining the optimal x-ray source trajectory for such short-scan acquisition.

PURPOSE

To quantify the projection angle-dependent normalized glandular dose coefficient (D g N C T $Dg{N^{CT}}$ ) in CBBCT, referred to as angularD g N C T $Dg{N^{CT}}$ , so that the x-ray ray source trajectory that minimizes the radiation dose to the breast for short-scan acquisition can be determined.

MATERIALS AND METHODS

A cohort of 75 CBBCT clinical datasets was segmented and used to generate three breast models - (I) patient-specific breast with heterogeneous fibroglandular tissue distribution and real breast shape, (II) patient-specific breast shape with homogeneous tissue distribution and matched fibroglandular weight fraction, and (III) homogeneous semi-ellipsoidal breast with patient-specific breast dimensions and matched fibroglandular weight fraction, which corresponds to the breast model used in current radiation dosimetry protocols. For each clinical dataset, the angularD g N C T $Dg{N^{CT}}$ was obtained at 10 discrete angles, spaced 36° apart, for full-scan, circular, x-ray source trajectory from Monte Carlo simulations. Model III is used for validating the Monte Carlo simulation results. Models II and III are used to determine if breast shape contributes to the observed trends in angularD g N C T $Dg{N^{CT}}$ . A geometry-based theory in conjunction with center-of-mass (C O M $COM$ ) based distribution analysis is used to explain the projection angle-dependent variation in angularD g N C T $Dg{N^{CT}}$ .

RESULTS

The theoretical model predicted that the angularD g N C T $Dg{N^{CT}}$ will follow a sinusoidal pattern and the amplitude of the sinusoid increases when the center-of-mass of fibroglandular tissue (C O M f $CO{M_f}$ ) is farther from the center-of-mass of the breast (C O M b $CO{M_b}$ ). It also predicted that the angularD g N C T $Dg{N^{CT}}$ will be minimized at x-ray source positions complementary to theC O M f $CO{M_f}$ . TheC O M f $CO{M_f}$ was superior to theC O M b $CO{M_b}$ in 80% (60/75) of the breasts. From Monte Carlo simulations and for homogeneous breasts (models II and III), the deviation in breast shape from a semi-ellipsoid had minimal effect on angularD g N C T $Dg{N^{CT}}$ and showed less than 4% variation. From Monte Carlo simulations and for model I, as predicted by our theory, the angularD g N C T $Dg{N^{CT}}$ followed a sinusoidal pattern with maxima and minima at x-ray source positions superior and inferior to the breast, respectively. For model I, the projection angle-dependent variation in angularD g N C T $Dg{N^{CT}}$ was 16.4%.

CONCLUSION

The heterogeneous tissue distribution affected the angularD g N C T $Dg{N^{CT}}$ more than the breast shape. For model I, the angularD g N C T $Dg{N^{CT}}$ was lowest when the x-ray source was inferior to the breast. Hence, for short-scan CBBCT acquisition withC O M b $CO{M_b}$ aligned with axis-of-rotation, an x-ray source trajectory inferior to the breast is preferable and such an acquisition spanning 205° can potentially reduce the mean glandular dose by up to 52%.

Artificial Intelligence in Breast X-Ray Imaging

This topical review is focused on the clinical breast x-ray imaging applications of the rapidly evolving field of artificial intelligence (AI). The range of AI applications is broad. AI can be used for breast cancer risk estimation that could allow for tailoring the screening interval and the protocol that are woman-specific and for triaging the screening exams. It also can serve as a tool to aid in the detection and diagnosis for improved sensitivity and specificity and as a tool to reduce radiologists' reading time. AI can also serve as a potential second 'reader' during screening interpretation. During the last decade, numerous studies have shown the potential of AI-assisted interpretation of mammography and to a lesser extent digital breast tomosynthesis; however, most of these studies are retrospective in nature. There is a need for prospective clinical studies to evaluate these technologies to better understand their real-world efficacy. Further, there are ethical, medicolegal, and liability concerns that need to be considered prior to the routine use of AI in the breast imaging clinic.

Breast dosimetry in alternative X-ray-based imaging modalities used in current clinical practices

In X-ray breast imaging, Digital Mammography (DM) and Digital Breast Tomosynthesis (DBT), are the standard and largely used techniques, both for diagnostic and screening purposes. Other techniques, such as dedicated Breast Computed Tomography (BCT) and Contrast Enhanced Mammography (CEM) have been developed as an alternative or a complementary technique to the established ones. The performance of these imaging techniques is being continuously assessed to improve the image quality and to reduce the radiation dose. These imaging modalities are predominantly used in the diagnostic setting to resolve incomplete or indeterminate findings detected with conventional screening examinations and could potentially be used either as an adjunct or as a primary screening tool in select populations, such as for women with dense breasts. The aim of this review is to describe the radiation dosimetry for these imaging techniques, and to compare the mean glandular dose with standard breast imaging modalities, such as DM and DBT.

Flat-panel conebeam CT in the clinic: history and current state

Research into conebeam CT concepts began as soon as the first clinical single-slice CT scanner was conceived. Early implementations of conebeam CT in the 1980s focused on high-contrast applications where concurrent high resolution ( < 200    μ m ), for visualization of small contrast-filled vessels, bones, or teeth, was an imaging requirement that could not be met by the contemporaneous CT scanners. However, the use of nonlinear imagers, e.g., x-ray image intensifiers, limited the clinical utility of the earliest diagnostic conebeam CT systems. The development of consumer-electronics large-area displays provided a technical foundation that was leveraged in the 1990s to first produce large-area digital x-ray detectors for use in radiography and then compact flat panels suitable for high-resolution and high-frame-rate conebeam CT. In this review, we show the concurrent evolution of digital flat panel (DFP) technology and clinical conebeam CT. We give a brief summary of conebeam CT reconstruction, followed by a brief review of the correction approaches for DFP-specific artifacts. The historical development and current status of flat-panel conebeam CT in four clinical areas-breast, fixed C-arm, image-guided radiation therapy, and extremity/head-is presented. Advances in DFP technology over the past two decades have led to improved visualization of high-contrast, high-resolution clinical tasks, and image quality now approaches the soft-tissue contrast resolution that is the standard in clinical CT. Future technical developments in DFPs will enable an even broader range of clinical applications; research in the arena of flat-panel CT shows no signs of slowing down.

Single-pass metal artifact reduction using a dual-layer flat panel detector

PURPOSE

Metal artifact remains a challenge in cone-beam CT images. Many image domain-based segmentation methods have been proposed for metal artifact reduction (MAR), which require two-pass reconstruction. Such methods first segment metal from a first-pass reconstruction and then forward-project the metal mask to identify them in projections. These methods work well in general but are limited when the metal is outside the scan field-of-view (FOV) or when the metal is moving during the scan. In the former, even reconstructing with a larger FOV does not guarantee a good estimate of metal location in the projections; and in the latter, the metal location in each projection is difficult to identify due to motion. Single-pass methods that detect metal in single-energy projections have also been developed, but often have imperfect metal detection that leads to residual artifacts. In this work, we develop a MAR method using a dual-layer (DL) flat panel detector, which improves performance for single-pass reconstruction.

METHODS

In this work, we directly detect metal objects in projections using dual-energy (DE) imaging that generates material-specific images (e.g., soft tissue and bone), where the metal stands out in bone images when nonuniform soft tissue background is removed. Metal is detected via simple thresholding, and entropy filtration is further applied to remove false-positive detections. A DL detector provides DE images with superior temporal and spatial registration and was used to perform the task. Scatter correction was first performed on DE raw projections to improve the accuracy of material decomposition. One phantom mimicking a liver biopsy setup and a cadaver head were used to evaluate the metal reduction performance of the proposed method and compared with that of a standard two-pass reconstruction, a previously published sinogram-based method using a Markov random field (MRF) model, and a single-pass projection-domain method using single-energy imaging. The phantom has a liver steering setup placed in a hollow chest phantom, with embedded metal and a biopsy needle crossing the phantom boundary. The cadaver head has dental fillings and a metal tag attached to its surface. The identified metal regions in each projection were corrected by interpolation using surrounding pixels, and the images were reconstructed using filtered backprojection.

RESULTS

Our current approach removes metal from the projections, which is robust to FOV truncation during imaging acquisition. In case of FOV truncation, the method outperformed the two-pass reconstruction method. The proposed method using DE renders better accuracy in metal segmentation than the MRF method and single-energy method, which were prone to false-positive errors that cause additional streaks. For the liver steering phantom, the average spatial nonuniformity was reduced from 0.127 in uncorrected images to 0.086 using a standard two-pass reconstruction and to 0.077 using the proposed method. For the cadaver head, the average standard deviation within selected soft tissue regions ( σ s ) was reduced from 209.1 HU in uncorrected images to 69.1 HU using a standard two-pass reconstruction and to 46.8 HU using our proposed method. The proposed method reduced the processing time by 31% as compared with the two-pass method.

CONCLUSIONS

We proposed a MAR method that directly detects metal in the projection domain using DE imaging, which is robust to truncation and superior to that of single-energy imaging. The method requires only a single-pass reconstruction that substantially reduces processing time compared with the standard two-pass metal reduction method.

Radiation dosimetry of a clinical prototype dedicated cone-beam breast CT system with offset detector

PURPOSE

A clinical-prototype, dedicated, cone-beam breast computed tomography (CBBCT) system with offset detector is undergoing clinical evaluation at our institution. This study is to estimate the normalized glandular dose coefficients ( DgN CT ) that provide air kerma-to-mean glandular dose conversion factors using Monte Carlo simulations.

MATERIALS AND METHODS

The clinical prototype CBBCT system uses 49 kV x-ray spectrum with 1.39 mm 1st half-value layer thickness. Monte Carlo simulations (GATE, version 8) were performed with semi-ellipsoidal, homogeneous breasts of various fibroglandular weight fractions ( f g = 0.01 , 0.15 , 0.5 , 1 ) , chest wall diameters ( d = 8 , 10 , 14 , 18 , 20  cm), and chest wall to nipple length ( l = 0.75 d ), aligned with the axis of rotation (AOR) located at 65 cm from the focal spot to determine the DgN CT . Three geometries were considered - 40 × 30 -cm detector with no offset that served as reference and corresponds to a clinical CBBCT system, 30 × 30 -cm detector with 5 cm offset, and a 30 × 30 -cm detector with 10 cm offset.

RESULTS

For 5 cm lateral offset, the DgN CT ranged 0.177 - 0.574  mGy/mGy and reduction in DgN CT with respect to reference geometry was observed only for 18 cm ( 6.4 % ± 0.23 % ) and 20 cm ( 9.6 % ± 0.22 % ) diameter breasts. For the 10 cm lateral offset, the DgN CT ranged 0.221 - 0.581  mGy/mGy and reduction in DgN CT was observed for all breast diameters. The reduction in DgN CT was 1.4 % ± 0.48 % , 7.1 % ± 0.13 % , 17.5 % ± 0.19 % , 25.1 % ± 0.15 % , and 27.7 % ± 0.08 % for 8, 10, 14, 18, and 20 cm diameter breasts, respectively. For a given breast diameter, the reduction in DgN CT with offset-detector geometries was not dependent on f g . Numerical fits of DgN CT d , l , f g were generated for each geometry.

CONCLUSION

The DgN CT and the numerical fit, D g N CT d , l , f g would be of benefit for current CBBCT systems using the reference geometry and for future generations using offset-detector geometry. There exists a potential for radiation dose reduction with offset-detector geometry, provided the same technique factors as the reference geometry are used, and the image quality is clinically acceptable.

Dedicated cone-beam breast CT using laterally-shifted detector geometry: Quantitative analysis of feasibility for clinical translation

BACKGROUND

High-resolution, low-noise detectors with minimal dead-space at chest-wall could improve posterior coverage and microcalcification visibility in the dedicated cone-beam breast CT (CBBCT). However, the smaller field-of-view necessitates laterally-shifted detector geometry to enable optimizing the air-gap for x-ray scatter rejection.

OBJECTIVE

To evaluate laterally-shifted detector geometry for CBBCT with clinical projection datasets that provide for anatomical structures and lesions.

METHODS

CBBCT projection datasets (n = 17 breasts) acquired with a 40×30 cm detector (1024×768-pixels, 0.388-mm pixels) were truncated along the fan-angle to emulate 20.3×30 cm, 22.2×30 cm and 24.1×30 cm detector formats and correspond to 20, 120, 220 pixels overlap in conjugate views, respectively. Feldkamp-Davis-Kress (FDK) algorithm with 3 different weighting schemes were used for reconstruction. Visual analysis for artifacts and quantitative analysis of root-mean-squared-error (RMSE), absolute difference between truncated and 40×30 cm reconstructions (Diff), and its power spectrum (PSDiff) were performed.

RESULTS

Artifacts were observed for 20.3×30 cm, but not for other formats. The 24.1×30 cm provided the best quantitative results with RMSE and Diff (both in units of μ, cm-1) of 4.39×10-3±1.98×10-3 and 4.95×10-4±1.34×10-4, respectively. The PSDiff (>0.3 cycles/mm) was in the order of 10-14μ2mm3 and was spatial-frequency independent.

CONCLUSIONS

Laterally-shifted detector CBBCT with at least 220 pixels overlap in conjugate views (24.1×30 cm detector format) provides quantitatively accurate and artifact-free image reconstruction.

X-ray scatter correction for dedicated cone beam breast CT using a forward-projection model

PURPOSE

The quality of dedicated cone-beam breast CT (CBBCT) imaging is fundamentally limited by x-ray scatter contamination due to the large irradiation volume. In this paper, we propose a scatter correction method for CBBCT using a novel forward-projection model with high correction efficacy and reliability.

METHOD

We first coarsely segment the uncorrected, first-pass, reconstructed CBBCT images into binary-object maps and assign the segmented fibroglandular and adipose tissue with the correct attenuation coefficients based on the mean x-ray energy. The modified CBBCT are treated as the prior images toward scatter correction. Primary signals are first estimated via forward projection on the modified CBBCT. To avoid errors caused by inaccurate segmentation, only sparse samples of estimated primary are selected for scatter estimation. A Fourier-Transform based algorithm, herein referred to as local filtration hereafter, is developed to efficiently estimate the global scatter distribution on the detector. The scatter-corrected images are obtained by removing the estimated scatter distribution from measured projection data.

RESULTS

We evaluate the method performance on six patients with different breast sizes and shapes representing the general population. The results show that the proposed method effectively reduces the image spatial non-uniformity from 8.27 to 1.91% for coronal views and from 6.50 to 3.00% for sagittal views. The contrast-to-deviation ratio is improved by an average factor of 1.41. Comparisons on the image details reveal that the proposed scatter correction successfully preserves fine structures of fibroglandular tissues that are lost in the segmentation process.

CONCLUSION

We propose a highly practical and efficient scatter correction algorithm for CBBCT via a forward-projection model. The method is attractive in clinical CBBCT imaging as it is readily implementable on a clinical system without modifications in current imaging protocols or system hardware.

Towards standardization of x-ray beam filters in digital mammography and digital breast tomosynthesis: Monte Carlo simulations and analytical modelling

In digital breast tomosynthesis and digital mammography, the x-ray beam filter material and thickness vary between systems. Replacing K-edge filters with Al was investigated with the intent to reduce exposure duration and to simplify system design. Tungsten target x-ray spectra were simulated with K-edge filters (50 µm Rh; 50 µm Ag) and Al filters of varying thickness. Monte Carlo simulations were conducted to quantify the x-ray scatter from various filters alone, scatter-to-primary ratio (SPR) with compressed breasts, and to determine the radiation dose to the breast. These data were used to analytically compute the signal-difference-to-noise ratio (SDNR) at unit (1 mGy) mean glandular dose (MGD) for W/Rh and W/Ag spectra. At SDNR matched between K-edge and Al filtered spectra, the reductions in exposure duration and MGD were quantified for three strategies: (i) fixed Al thickness and matched tube potential in kilovolts (kV); (ii) fixed Al thickness and varying the kV to match the half-value layer (HVL) between Al and K-edge filtered spectra; and, (iii) matched kV and varying the Al thickness to match the HVL between Al and K-edge filtered spectra. Monte Carlo simulations indicate that the SPR with and without the breast were not different between Al and K-edge filters. Modelling for fixed Al thickness (700 µm) and kV matched to K-edge filtered spectra, identical SDNR was achieved with 37-57% reduction in exposure duration and with 2-20% reduction in MGD, depending on breast thickness. Modelling for fixed Al thickness (700 µm) and HVL matched by increasing the kV over (0,4) range, identical SDNR was achieved with 62-65% decrease in exposure duration and with 2-24% reduction in MGD, depending on breast thickness. For kV and HVL matched to K-edge filtered spectra by varying Al filter thickness over (700, 880) µm range, identical SDNR was achieved with 23-56% reduction in exposure duration and 2-20% reduction in MGD, depending on breast thickness. These simulations indicate that increased fluence with Al filter of fixed or variable thickness substantially decreases exposure duration while providing for similar image quality with moderate reduction in MGD.

Dosimetry in x-ray-based breast imaging

The estimation of the mean glandular dose to the breast (MGD) for x-ray based imaging modalities forms an essential part of quality control and is needed for risk estimation and for system design and optimisation. This review considers the development of methods for estimating the MGD for mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (DBCT). Almost all of the methodology used employs Monte Carlo calculated conversion factors to relate the measurable quantity, generally the incident air kerma, to the MGD. After a review of the size and composition of the female breast, the various mathematical models used are discussed, with particular emphasis on models for mammography. These range from simple geometrical shapes, to the more recent complex models based on patient DBCT examinations. The possibility of patient-specific dose estimates is considered as well as special diagnostic views and the effect of breast implants. Calculations using the complex models show that the MGD for mammography is overestimated by about 30% when the simple models are used. The design and uses of breast-simulating test phantoms for measuring incident air kerma are outlined and comparisons made between patient and phantom-based dose estimates. The most widely used national and international dosimetry protocols for mammography are based on different simple geometrical models of the breast, and harmonisation of these protocols using more complex breast models is desirable.

Library based x-ray scatter correction for dedicated cone beam breast CT

PURPOSE

The image quality of dedicated cone beam breast CT (CBBCT) is limited by substantial scatter contamination, resulting in cupping artifacts and contrast-loss in reconstructed images. Such effects obscure the visibility of soft-tissue lesions and calcifications, which hinders breast cancer detection and diagnosis. In this work, we propose a library-based software approach to suppress scatter on CBBCT images with high efficiency, accuracy, and reliability.

METHODS

The authors precompute a scatter library on simplified breast models with different sizes using the geant4-based Monte Carlo (MC) toolkit. The breast is approximated as a semiellipsoid with homogeneous glandular/adipose tissue mixture. For scatter correction on real clinical data, the authors estimate the breast size from a first-pass breast CT reconstruction and then select the corresponding scatter distribution from the library. The selected scatter distribution from simplified breast models is spatially translated to match the projection data from the clinical scan and is subtracted from the measured projection for effective scatter correction. The method performance was evaluated using 15 sets of patient data, with a wide range of breast sizes representing about 95% of general population. Spatial nonuniformity (SNU) and contrast to signal deviation ratio (CDR) were used as metrics for evaluation.

RESULTS

Since the time-consuming MC simulation for library generation is precomputed, the authors' method efficiently corrects for scatter with minimal processing time. Furthermore, the authors find that a scatter library on a simple breast model with only one input parameter, i.e., the breast diameter, sufficiently guarantees improvements in SNU and CDR. For the 15 clinical datasets, the authors' method reduces the average SNU from 7.14% to 2.47% in coronal views and from 10.14% to 3.02% in sagittal views. On average, the CDR is improved by a factor of 1.49 in coronal views and 2.12 in sagittal views.

CONCLUSIONS

The library-based scatter correction does not require increase in radiation dose or hardware modifications, and it improves over the existing methods on implementation simplicity and computational efficiency. As demonstrated through patient studies, the authors' approach is effective and stable, and is therefore clinically attractive for CBBCT imaging.

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.

Cone beam CT for determining breast cancer margin: an initial experience and its comparison with mammography and specimen radiograph

PURPOSE

To assess the ability of cone beam CT (CBCT) in determining the breast cancer margin using, to compare the results with mammography and specimen radiography, and to explore the clinical potential of CBCT for breast imaging.

METHODS

Specimens of 46 breast cancer patients were imaged by using a prototype CBCT system. Each patient underwent mammography, CBCT and X-ray of breast surgical specimen within 6 months. Images of mammography, breast surgical specimen radiography and CBCT were evaluated by an experienced radiologist. Indicators, such as: morphology, glitch, density, invasion, structural distortion and calcification, were observed.

RESULT

There was no significant difference of the calcification, glitch and morphology among three methods. However, there was significant difference in indicators of breast tumor invasion among three methods. There was statistical significance in detecting invasions of breast cancer cells in peripheral tissues among three methods.

CONCLUSION

CBCT shows no superiority over mammography and specimen radiography in determining tumor's outline and detecting calcification. On the other hand, CBCT demonstrates its advantage in determining the 3 dimensional position of a lesion which could be a potential clinical application in future practices of breast imaging.

Mean glandular dose coefficients (D(g)N) for x-ray spectra used in contemporary breast imaging systems

To develop tables of normalized glandular dose coefficients D(g)N for a range of anode-filter combinations and tube voltages used in contemporary breast imaging systems. Previously published mono-energetic D(g)N values were used with various spectra to mathematically compute D(g)N coefficients. The tungsten anode spectra from TASMICS were used; molybdenum and rhodium anode-spectra were generated using MCNPX Monte Carlo code. The spectra were filtered with various thicknesses of Al, Rh, Mo or Cu. An initial half value layer (HVL) calculation was made using the anode and filter material. A range of the HVL values was produced with the addition of small thicknesses of polymethyl methacrylate (PMMA) as a surrogate for the breast compression paddle, to produce a range of HVL values at each tube voltage. Using a spectral weighting method, D(g)N coefficients for the generated spectra were calculated for breast glandular densities of 0%, 12.5%, 25%, 37.5%, 50% and 100% for a range of compressed breast thicknesses from 3 to 8 cm. Eleven tables of normalized glandular dose (D(g)N) coefficients were produced for the following anode/filter combinations: W + 50 μm Ag, W + 500 μm Al, W + 700 μm Al, W + 200 μm Cu, W + 300 μm Cu, W + 50 μm Rh, Mo + 400 μm Cu, Mo + 30 μm Mo, Mo + 25 μm Rh, Rh + 400 μm Cu and Rh + 25 μm Rh. Where possible, these results were compared to previously published D(g)N values and were found to be on average less than 2% different than previously reported values.Over 200 pages of D(g)N coefficients were computed for modeled x-ray system spectra that are used in a number of new breast imaging applications. The reported values were found to be in excellent agreement when compared to published values.

Large-angle x-ray scatter in Talbot-Lau interferometry for breast imaging

Monte Carlo simulations were used to investigate large-angle x-ray scatter at design energy of 25 keV during small field of view (9.6 cm × 5 cm) differential phase contrast imaging of the breast using Talbot-Lau interferometry. Homogenous, adipose and fibroglandular breasts of uniform thickness ranging from 2 to 8 cm encompassing the field of view were modeled. Theoretically determined transmission efficiencies of the gratings were used to validate the Monte Carlo simulations, followed by simulations to determine the x-ray scatter reaching the detector. The recorded x-ray scatter was classified into x-ray photons that underwent at least one Compton interaction (incoherent scatter) and Rayleigh interaction alone (coherent scatter) for further analysis. Monte Carlo based estimates of transmission efficiencies showed good correspondence [Formula: see text] with theoretical estimates. Scatter-to-primary ratio increased with increasing breast thickness, ranging from 0.11 to 0.22 for 2-8 cm thick adipose breasts and from 0.12 to 0.28 for 2-8 cm thick fibroglandular breasts. The analyzer grating reduced incoherent scatter by ~18% for 2 cm thick adipose breast and by ~35% for 8 cm thick fibroglandular breast. Coherent scatter was the dominant contributor to the total scatter. Coherent-to-incoherent scatter ratio ranged from 2.2 to 3.1 for 2-8 cm thick adipose breasts and from 2.7 to 3.4 for 2-8 cm thick fibroglandular breasts.

Personalized estimates of radiation dose from dedicated breast CT in a diagnostic population and comparison with diagnostic mammography

This study retrospectively analyzed the mean glandular dose (MGD) to 133 breasts from 132 subjects, all women, who participated in a clinical trial evaluating dedicated breast CT in a diagnostic population. The clinical trial was conducted in adherence to a protocol approved by institutional review boards and the study participants provided written informed consent. Individual estimates of MGD to each breast from dedicated breast CT was obtained by combining x-ray beam characteristics with estimates of breast dimensions and fibroglandular fraction from volumetric breast CT images, and using normalized glandular dose coefficients. For each study participant and for the breast corresponding to that imaged with breast CT, an estimate of the MGD from diagnostic mammography (including supplemental views) was obtained from the DICOM image headers for comparison. This estimate uses normalized glandular dose coefficients corresponding to a breast with 50% fibroglandular weight fraction. The median fibroglandular weight fraction for the study cohort determined from volumetric breast CT images was 15%. Hence, the MGD from diagnostic mammography was corrected to be representative of the study cohort. Individualized estimates of MGD from breast CT ranged from 5.7 to 27.8 mGy. Corresponding to the breasts imaged with breast CT, the MGD from diagnostic mammography ranged from 2.6 to 31.6 mGy. The mean (± inter-breast SD) and the median MGD (mGy) from dedicated breast CT exam were 13.9 ± 4.6 and 12.6, respectively. For the corresponding breasts, the mean (± inter-breast SD) and the median MGD (mGy) from diagnostic mammography were 12.4 ± 6.3 and 11.1, respectively. Statistical analysis indicated that at the 0.05 level, the distributions of MGD from dedicated breast CT and diagnostic mammography were significantly different (Wilcoxon signed ranks test, p = 0.007). While the interquartile range and the range (maximum-minimum) of MGD from dedicated breast CT was lower than diagnostic mammography, the median MGD from dedicated breast CT was approximately 13.5% higher than that from diagnostic mammography. The MGD for breast CT is based on a 1.45 mm skin layer and that for diagnostic mammography is based on a 4 mm skin layer; thus, favoring a lower estimate for MGD from diagnostic mammography. The median MGD from dedicated breast CT corresponds to the median MGD from four to five diagnostic mammography views. In comparison, for the same 133 breasts, the mean and the median number of views per breast during diagnostic mammography were 4.53 and 4, respectively. Paired analysis showed that there was approximately equal likelihood of receiving lower MGD from either breast CT or diagnostic mammography. Future work will investigate methods to reduce and optimize radiation dose from dedicated breast CT.

Scaling-law for the energy dependence of anatomic power spectrum in dedicated breast CT

PURPOSE

To determine the x-ray photon energy dependence of the anatomic power spectrum of the breast when imaged with dedicated breast computed tomography (CT).

METHODS

A theoretical framework for scaling the empirically determined anatomic power spectrum at one x-ray photon energy to that at any given x-ray photon energy when imaged with dedicated breast CT was developed. Theory predicted that when the anatomic power spectrum is fitted with a power curve of the form k f(-β), where k and β are fit coefficients and f is spatial frequency, the exponent β would be independent of x-ray photon energy (E), and the amplitude k scales with the square of the difference in energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues. Twenty mastectomy specimens based numerical phantoms that were previously imaged with a benchtop flat-panel cone-beam CT system were converted to 3D distribution of glandular weight fraction (f(g)) and were used to verify the theoretical findings. The 3D power spectrum was computed in terms of f(g) and after converting to linear attenuation coefficients at monoenergetic x-ray photon energies of 20-80 keV in 5 keV intervals. The 1D power spectra along the axes were extracted and fitted with a power curve of the form k f(-β). The energy dependence of k and β were analyzed.

RESULTS

For the 20 mastectomy specimen based numerical phantoms used in the study, the exponent β was found to be in the range of 2.34-2.42, depending on the axis of measurement. Numerical simulations agreed with the theoretical predictions that for a power-law anatomic spectrum of the form k f(-β), β was independent of E and k(E) = k(1)[μ(g)(E) - μ(a)(E)](2), where k(1) is a constant, and μ(g)(E) and μ(a)(E) represent the energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues, respectively.

CONCLUSIONS

Numerical simulations confirmed the theoretical predictions that in dedicated breast CT, the spatial frequency dependence of the anatomic power spectrum will be independent of x-ray photon energy, and the amplitude of the anatomic power spectrum scales by the square of difference in linear attenuation coefficients of fibroglandular and adipose tissues.

X-ray phase contrast imaging of the breast: analysis of tissue simulating materials

PURPOSE

Phase contrast imaging, particularly of the breast, is being actively investigated. The purpose of this work is to investigate the x-ray phase contrast properties of breast tissues and commonly used breast tissue substitutes or phantom materials with an aim of determining the phantom materials best representative of breast tissues.

METHODS

Elemental compositions of breast tissues including adipose, fibroglandular, and skin were used to determine the refractive index, n = 1 - δ + i β. The real part of the refractive index, specifically the refractive index decrement (δ), over the energy range of 5-50 keV were determined using XOP software (version 2.3, European Synchrotron Radiation Facility, France). Calcium oxalate and calcium hydroxyapatite were considered to represent the material compositions of microcalcifications in vivo. Nineteen tissue substitutes were considered as possible candidates to represent adipose tissue, fibroglandular tissue and skin, and four phantom materials were considered as possible candidates to represent microcalcifications. For each material, either the molecular formula, if available, or the elemental composition based on weight fraction, was used to determine δ. At each x-ray photon energy, the absolute percent difference in δ between the breast tissue and the substitute material was determined, from which three candidates were selected. From these candidate tissue substitutes, the material that minimized the absolute percent difference in linear attenuation coefficient μ, and hence β, was considered to be best representative of that breast tissue.

RESULTS

Over the energy range of 5-50 keV, while the δ of CB3 and fibroglandular tissue-equivalent material were within 1% of that of fibroglandular tissue, the μ of fibroglandular tissue-equivalent material better approximated the fibroglandular tissue. While the δ of BR10 and adipose tissue-equivalent material were within 1% of that of adipose tissue, the tissue-equivalent material better approximated the adipose tissue in terms of μ. Polymethyl methacrylate, a commonly used tissue substitute, exhibited δ greater than fibroglandular tissue by ≈ 12%. The A-150 plastic closely approximated the skin. Several materials exhibited δ between that of adipose and fibroglandular tissue. However, there was an energy-dependent mismatch in terms of equivalent fibroglandular weight fraction between δ and μ for these materials. For microcalcifications, aluminum and calcium carbonate were observed to straddle the δ and μ of calcium oxalate and calcium hydroxyapatite. Aluminum oxide, commonly used to represent microcalcifications in the American College of Radiology recommended phantoms for accreditation exhibited δ greater than calcium hydroxyapatite by ≈ 23%.

CONCLUSIONS

A breast phantom comprising A-150 plastic to represent the skin, commercially available adipose and fibroglandular tissue-equivalent formulations to represent adipose and fibroglandular tissue, respectively, was found to be best suited for x-ray phase-sensitive imaging of the breast. Calcium carbonate or aluminum can be used to represent microcalcifications.

Dedicated breast CT: geometric design considerations to maximize posterior breast coverage

An Institutional Review Board-approved protocol was used to quantify breast tissue inclusion in 52 women, under conditions simulating both craniocaudal (CC) and mediolateral oblique (MLO) views in mammography, dedicated breast CT in the upright subject position, and dedicated breast CT in the prone subject position. Using skin as a surrogate for the underlying breast tissue, the posterior aspect of the breast that is aligned with the chest-wall edge of the breast support in a screen-film mammography system was marked with the study participants positioned for CC and MLO views. The union of skin marks with the study participants positioned for CC and MLO views was considered to represent chest-wall tissue available for imaging with mammography and served as the reference standard. For breast CT, a prone stereotactic breast biopsy unit and a custom-fabricated barrier were used to simulate conditions during prone and upright breast CT, respectively. For the same breast marked on the mammography system, skin marks were made along the breast periphery that was just anterior to the apertures of the prone biopsy unit and the upright barrier. The differences in skin marks between subject positioning simulating breast CT (prone, upright) and mammography were quantified at six anatomic locations. For each location, at least one study participant had a skin mark from breast CT (prone, upright) posterior to mammography. However for all study participants, there was at least one anatomic location where the skin mark from mammography was posterior to that from breast CT (prone, upright) positioning. The maximum amount by which the skin mark from mammography was posterior to breast CT (prone and upright) over all six locations was quantified for each study participant and pair-wise comparison did not exhibit statistically significant difference between prone and upright breast CT (paired t- test, p = 0.4). Quantitatively, for 95% of the study participants the skin mark from mammography was posterior to breast CT (prone or upright) by at the most 9 mm over all six locations. Based on the study observations, geometric design considerations targeting chest-wall coverage with breast CT equivalent to mammography, wherein part of the x-ray beam images through the swale during breast CT are provided. Assuming subjects can extend their chest in to a swale, the optimal swale-depth required to achieve equivalent coverage with breast CT images as mammograms for 95% of the subjects varies in the range of ~30-50 mm for clinical prototypes and was dependent on the system geometry.

Dedicated breast CT: fibroglandular volume measurements in a diagnostic population

PURPOSE

To determine the mean and range of volumetric glandular fraction (VGF) of the breast in a diagnostic population using a high-resolution flat-panel cone-beam dedicated breast CT system. This information is important for Monte Carlo-based estimation of normalized glandular dose coefficients and for investigating the dependence of VGF on breast dimensions, race, and pathology.

METHODS

Image data from a clinical trial investigating the role of dedicated breast CT that enrolled 150 women were retrospectively analyzed to determine the VGF. The study was conducted in adherence to a protocol approved by the institutional human subjects review boards and written informed consent was obtained from all study participants. All participants in the study were assigned BI-RADS(®) 4 or 5 as per the American College of Radiology assessment categories after standard diagnostic work-up and underwent dedicated breast CT exam prior to biopsy. A Gaussian-kernel based fuzzy c-means algorithm was used to partition the breast CT images into adipose and fibroglandular tissue after segmenting the skin. Upon determination of the accuracy of the algorithm with a phantom, it was applied to 137 breast CT volumes from 136 women. VGF was determined for each breast and the mean and range were determined. Pathology results with classification as benign, malignant, and hyperplasia were available for 132 women, and were used to investigate if the distributions of VGF varied with pathology.

RESULTS

The algorithm was accurate to within ±1.9% in determining the volume of an irregular shaped phantom. The study mean (± inter-breast SD) for the VGF was 0.172 ± 0.142 (range: 0.012-0.719). VGF was found to be negatively correlated with age, breast dimensions (chest-wall to nipple length, pectoralis to nipple length, and effective diameter at chest-wall), and total breast volume, and positively correlated with fibroglandular volume. Based on pathology, pairwise statistical analysis (Mann-Whitney test) indicated that at the 0.05 significance level, there was no significant difference in distributions of VGF without adjustment for age between malignant and nonmalignant breasts (p = 0.41). Pairwise comparisons of the distributions of VGF in increasing order of mammographic breast density indicated all comparisons were statistically significant (p < 0.002).

CONCLUSIONS

This study used a different clinical prototype breast CT system than that in previous studies to image subjects from a different geographical region, and used a different algorithm for analysis of image data. The mean VGF estimated from this study is within the range reported in previous studies, indicating that the choice of 50% glandular weight fraction to represent an average breast for Monte Carlo-based estimation of normalized glandular dose coefficients in mammography needs revising. In the study, the distributions of VGF did not differ significantly with pathology.

Characterization of the homogeneous tissue mixture approximation in breast imaging dosimetry

PURPOSE

To compare the estimate of normalized glandular dose in mammography and breast CT imaging obtained using the actual glandular tissue distribution in the breast to that obtained using the homogeneous tissue mixture approximation.

METHODS

Twenty volumetric images of patient breasts were acquired with a dedicated breast CT prototype system and the voxels in the breast CT images were automatically classified into skin, adipose, and glandular tissue. The breasts in the classified images underwent simulated mechanical compression to mimic the conditions present during mammographic acquisition. The compressed thickness for each breast was set to that achieved during each patient's last screening cranio-caudal (CC) acquisition. The volumetric glandular density of each breast was computed using both the compressed and uncompressed classified images, and additional images were created in which all voxels representing adipose and glandular tissue were replaced by a homogeneous mixture of these two tissues in a proportion corresponding to each breast's volumetric glandular density. All four breast images (compressed and uncompressed; heterogeneous and homogeneous tissue) were input into Monte Carlo simulations to estimate the normalized glandular dose during mammography (compressed breasts) and dedicated breast CT (uncompressed breasts). For the mammography simulations the x-ray spectra used was that used during each patient's last screening CC acquisition. For the breast CT simulations, two x-ray spectra were used, corresponding to the x-ray spectra with the lowest and highest energies currently being used in dedicated breast CT prototype systems under clinical investigation. The resulting normalized glandular dose for the heterogeneous and homogeneous versions of each breast for each modality was compared.

RESULTS

For mammography, the normalized glandular dose based on the homogeneous tissue approximation was, on average, 27% higher than that estimated using the true heterogeneous glandular tissue distribution (Wilcoxon Signed Rank Test p = 0.00046). For dedicated breast CT, the overestimation of normalized glandular dose was, on average, 8% (49 kVp spectrum, p = 0.00045) and 4% (80 kVp spectrum, p = 0.000089). Only two cases in mammography and two cases in dedicated breast CT with a tube voltage of 49 kVp resulted in lower dose estimates for the homogeneous tissue approximation compared to the heterogeneous tissue distribution.

CONCLUSIONS

The normalized glandular dose based on the homogeneous tissue mixture approximation results in a significant overestimation of dose to the imaged breast. This overestimation impacts the use of dose estimates in absolute terms, such as for risk estimates, and may impact some comparative studies, such as when modalities or techniques with different x-ray energies are used. The error introduced by the homogeneous tissue mixture approximation in higher energy x-ray modalities, such as dedicated breast CT, although statistically significant, may not be of clinical concern. Further work is required to better characterize this overestimation and potentially develop new metrics or correction factors to better estimate the true glandular dose to breasts undergoing imaging with ionizing radiation.