Generation and analysis of clinically relevant breast imaging x-ray spectra.

PURPOSE: The purpose of this work was to develop and make available x-ray spectra for some of the most widely used digital mammography (DM), breast tomosynthesis (BT), and breast CT (bCT) systems in North America. METHODS: The Monte Carlo code MCNP6 was used to simulate minimally filtered (only beryllium) x-ray spectra at 8 tube potentials from 20 to 49 kV for DM/BT, and 9 tube potentials from 35 to 70 kV for bCT. Vendor-specific anode compositions, effective anode angles, focal spot sizes, source-to-detector distances, and beryllium filtration were simulated. For each 0.5 keV energy bin in all simulated spectra, the fluence was interpolated using cubic splines across the range of simulated tube potentials to produce spectra in 1 kV increments from 20 to 49 kV for DM/BT and from 35 to 70 kV for bCT. The HVL of simulated spectra with conventional filtration (at 35 kV for DM/BT and 49 kV for bCT) was used to assess spectral differences resulting from variations in: (a) focal spot size (0.1 and 0.3 mm IEC), (b) solid angle at the detector (i.e., small and large FOV size), and (c) geometrical specifications for vendors that employ the same anode composition. RESULTS: Averaged across all DM/BT vendors, variations in focal spot and FOV size resulted in HVL differences of 2.2% and 0.9%, respectively. Comparing anode compositions separately, the HVL differences for Mo (GE, Siemens) and W (Hologic, Philips, and Siemens) spectra were 0.3% and 0.6%, respectively. Both the commercial Koning and prototype "Doheny" (UC Davis) bCT systems utilize W anodes with a 0.3 mm focal spot. Averaged across both bCT systems, variations in FOV size resulted in a 2.2% difference in HVL. In addition, the Koning spectrum was slightly harder than Doheny with a 4.2% difference in HVL. Therefore to reduce redundancy, a generic DM/BT system and a generic bCT system were used to generate the new spectra reported herein. The spectral models for application to DM/BT were dubbed the Molybdenum, Rhodium, and Tungsten Anode Spectral Models using Interpolating Cubic Splines (MASMICSM-T , RASMICSM-T , and TASMICSM-T ; subscript "M-T" indicating mammography and tomosynthesis). When compared against reference models (MASMIPM , RASMIPM , and TASMIPM ; subscript "M" indicating mammography), the new spectral models were in close agreement with mean differences of 1.3%, -1.3%, and -3.3%, respectively, across tube potential comparisons of 20, 30, and 40 kV with conventional filtration. TASMICSbCT -generated bCT spectra were also in close agreement with the reference TASMIP model with a mean difference of -0.8%, across tube potential comparisons of 35, 49, and 70 kV with 1.5 mm Al filtration. CONCLUSIONS: The Mo, Rh, and W anode spectra for application in DM and BT (MASMICSM-T , RASMICSM-T , and TASMICSM-T ) and the W anode spectra for bCT (TASMICSbCT ) as described in this study should be useful for individuals interested in modeling the performance of modern breast x-ray imaging systems including dual-energy mammography which extends to 49 kV. These new spectra are tabulated in spreadsheet form and are made available to any interested party.

Differentiation of ductal carcinoma in-situ from benign micro-calcifications by dedicated breast computed tomography.

PURPOSE: Compare conspicuity of ductal carcinoma in-situ (DCIS) to benign calcifications on unenhanced (bCT), contrast-enhanced dedicated breast CT (CEbCT) and mammography (DM). METHODS AND MATERIALS: The institutional review board approved this HIPAA-compliant study. 42 women with Breast Imaging Reporting and Data System 4 or 5 category micro-calcifications had breast CT before biopsy. Three subjects with invasive disease at surgery were excluded. Two breast radiologists independently compared lesion conspicuity scores (CS) for CEbCT, to bCT and DM. Enhancement was measured in Hounsfield units (HU). Mean CS ± standard deviations are shown. Receiver operating characteristic analysis (ROC) measured radiologists' discrimination performance by comparing CS to enhancement alone. Statistical measurements were made using ANOVA F-test, Wilcoxon rank-sum test and robust linear regression analyses. RESULTS: 39 lesions (17 DCIS, 22 benign) were analyzed. DCIS (8.5 ± 0.9, n=17) was more conspicuous than benign micro-calcifications (3.6 ± 2.9, n=22; p<0.0001) on CEbCT. DCIS was equally conspicuous on CEbCT and DM (8.5 ± 0.9, 8.7 ± 0.8, n=17; p=0.85) and more conspicuous when compared to bCT (5.3 ± 2.6, n=17; p<0.001). All DCIS enhanced; mean enhancement (90HU ± 53HU, n=17) was higher compared to benign lesions (33 ± 30HU, n=22) (p<0.0001). ROC analysis of the radiologists' CS showed high discrimination performance (AUC=0.94) compared to enhancement alone (AUC=0.85) (p<0.026). CONCLUSION: DCIS is more conspicuous than benign micro-calcifications on CEbCT. DCIS visualization on CEbCT is equal to mammography but improved compared to bCT. Radiologists' discrimination performance using CEBCT is significantly higher than enhancement values alone. CEbCT may have an advantage over mammography by reducing false positive examinations when calcifications are analyzed.

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.

Radiation dosimetry in digital breast tomosynthesis: report of AAPM Tomosynthesis Subcommittee Task Group 223.

The radiation dose involved in any medical imaging modality that uses ionizing radiation needs to be well understood by the medical physics and clinical community. This is especially true of screening modalities. Digital breast tomosynthesis (DBT) has recently been introduced into the clinic and is being used for screening for breast cancer in the general population. Therefore, it is important that the medical physics community have the required information to be able to understand, estimate, and communicate the radiation dose levels involved in breast tomosynthesis imaging. For this purpose, the American Association of Physicists in Medicine Task Group 223 on Dosimetry in Tomosynthesis Imaging has prepared this report that discusses dosimetry in breast imaging in general, and describes a methodology and provides the data necessary to estimate mean breast glandular dose from a tomosynthesis acquisition. In an effort to maximize familiarity with the procedures and data provided in this Report, the methodology to perform the dose estimation in DBT is based as much as possible on that used in mammography dose estimation.

Comprehensive assessment of the slice sensitivity profiles in breast tomosynthesis and breast CT.

PURPOSE: This study experimentally evaluated the slice sensitivity profile (SSP) and its relationship between acquisition angle, object size, and cone angle. The sensitivity profile metric was used to characterize a breast tomosynthesis system's resolution in the z-axis. The SSP was also measured on a prototype breast computed tomography (bCT) system. METHODS: The SSP was measured using brass disks placed within adipose tissue-equivalent breast phantoms. The digital tomosynthesis system (Selenia Dimensions, Hologic Corporation, Bedford, MA) acquires projection images over a 15° angular range and the bCT scanner acquires projection images over a 360° angular range. Angular ranges between 15° and 360° were studied by using a subset of the projection images acquired on the bCT scanner. The SSP was determined by measuring a background-corrected mean gray scale value as a function of the z-position (axis normal to the plane of the detector). RESULTS: The results show that SSP improves when the angular acquisition range is increased and the SSP approaches a delta function for angles greater than 180°. Smaller objects have a narrower SSP and the SSP is not significantly dependent on the cone angle. For a 2.5, 5, 10 mm disk, the full width at half maximum of the SSP was 35, 61, 115 mm, respectively, on the tomosynthesis system (at 15°) and was 0.5 mm for all disk diameters on the bCT scanner (at 360°). CONCLUSIONS: The SSP is dependent on object size and angular acquisition range. These dependencies are overcome once the angular acquisition range is increased beyond 180°.

Anatomical complexity in breast parenchyma and its implications for optimal breast imaging strategies.

PURPOSE: The purpose of this investigation was to assess the anatomical noise in breast images using a mathematically derived parameter ß as a surrogate for detection performance, across the same patient cohort but in different imaging modalities including mammography, tomosynthesis, and breast CT. METHODS: Women who were scheduled for breast biopsy were approached for participation in this IRB and HIPPA-compliant investigation. A total of 23 women had all views of each modality and represent the cohort studied in this investigation. Image data sets across all modalities were analyzed using 1000 regions of interest per image data set, and the anatomical noise power spectrum, NPS(a)(f), was computed and averaged for each breast image data set. After windowing the total noise power spectrum NPS(t)(f) to a specific frequency range corresponding to anatomical noise, the power-law slope (ß) of the NPS(a)(f) was computed where NPS(a)(f) = α f(-) (ß). RESULTS: The value of ß was determined for breast CT data sets, and they were 1.75 (0.424), 1.83 (0.352), and 1.79 (0.397), for the coronal, sagittal, and axial views, respectively. For tomosynthesis, ß was 3.06 (0.361) and 3.10 (0.315) for the craniocaudal (CC) and medial lateral oblique (MLO) views, respectively. For mammography, these values were 3.17 (0.226) and 3.30 (0.236), for the CC and MLO views, respectively. The values of ß for breast CT were significantly different than those for tomosynthesis and mammography (p < 0.001, all 12 comparisons). CONCLUSIONS: Based on the parameter ß which is thought to describe anatomical noise in breast images, breast CT was shown to have a statistically significant lower ß than mammography or tomosynthesis. It has been suggested in the literature that a lower ß may correspond to increased cancer detection performance; however, this has yet to be demonstrated unequivocally.

Experimental validation of a method characterizing bow tie filters in CT scanners using a real-time dose probe.

PURPOSE: Beam-shaping or "bow tie" (BT) filters are used to spatially modulate the x-ray beam in a CT scanner, but the conventional method of step-and-shoot measurement to characterize a beam's profile is tedious and time-consuming. The theory for characterization of bow tie relative attenuation (COBRA) method, which relies on a real-time dosimeter to address the issues of conventional measurement techniques, was previously demonstrated using computer simulations. In this study the feasibility of the COBRA theory is further validated experimentally through the employment of a prototype real-time radiation meter and a known BT filter. METHODS: The COBRA method consisted of four basic steps: (1) The probe was placed at the edge of a scanner's field of view; (2) a real-time signal train was collected as the scanner's gantry rotated with the x-ray beam on; (3) the signal train, without a BT filter, was modeled using peak values measured in the signal train of step 2; and (4) the relative attenuation of the BT filter was estimated from filtered and unfiltered data sets. The prototype probe was first verified to have an isotropic and linear response to incident x-rays. The COBRA method was then tested on a dedicated breast CT scanner with a custom-designed BT filter and compared to the conventional step-and-shoot characterization of the BT filter. Using basis decomposition of dual energy signal data, the thickness of the filter was estimated and compared to the BT filter's manufacturing specifications. The COBRA method was also demonstrated with a clinical whole body CT scanner using the body BT filter. The relative attenuation was calculated at four discrete x-ray tube potentials and used to estimate the thickness of the BT filter. RESULTS: The prototype probe was found to have a linear and isotropic response to x-rays. The relative attenuation produced from the COBRA method fell within the error of the relative attenuation measured with the step-and-shoot method. The BT filter thickness estimates resulting from the dual energy scans on the breast CT system were equivalent to the manufacturing specifications. The clinical CT evaluation produced data conceptually similar to previous computer simulations and plausible relative attenuation profiles were observed. CONCLUSIONS: The COBRA method is a fast and accurate method for BT filter characterization, which requires a simple experimental setup in a clinical environment. Because of the ease of data acquisition, multienergy scans can be acquired which allow characterization of the BT filter thickness.