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 (MASMICS_M-T , RASMICS_M-T , and TASMICS_M-T ; subscript "M-T" indicating mammography and tomosynthesis). When compared against reference models (MASMIP_M , RASMIP_M , and TASMIP_M ; 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. TASMICS_bCT -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 (MASMICS_M-T , RASMICS_M-T , and TASMICS_M-T ) and the W anode spectra for bCT (TASMICS_bCT ) 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.

Simulated lesion, human observer performance comparison between thin-section dedicated breast CT images versus computed thick-section simulated projection images of the breast

The objective of this study was to compare the lesion detection performance of human observers between thin-section computed tomography images of the breast, with thick-section (>40 mm) simulated projection images of the breast. Three radiologists and six physicists each executed a two alterative force choice (2AFC) study involving simulated spherical lesions placed mathematically into breast images produced on a prototype dedicated breast CT scanner. The breast image data sets from 88 patients were used to create 352 pairs of image data. Spherical lesions with diameters of 1, 2, 3, 5, and 11 mm were simulated and adaptively positioned into 3D breast CT image data sets; the native thin section (0.33 mm) images were averaged to produce images with different slice thicknesses; average section thicknesses of 0.33, 0.71, 1.5 and 2.9 mm were representative of breast CT; the average 43 mm slice thickness served to simulate simulated projection images of the breast.The percent correct of the human observer's responses were evaluated in the 2AFC experiments. Radiologists lesion detection performance was significantly (p < 0.05) better in the case of thin-section images, compared to thick section images similar to mammography, for all but the 1 mm lesion diameter lesions. For example, the average of three radiologist's performance for 3 mm diameter lesions was 92% correct for thin section breast CT images while it was 67% for the simulated projection images. A gradual reduction in observer performance was observed as the section thickness increased beyond about 1 mm. While a performance difference based on breast density was seen in both breast CT and the projection image results, the average radiologist performance using breast CT images in dense breasts outperformed the performance using simulated projection images in fatty breasts for all lesion diameters except 11 mm. The average radiologist performance outperformed that of the average physicist observer, however trends in performance were similar. Human observers demonstrate significantly better mass-lesion detection performance on thin-section CT images of the breast, compared to thick-section simulated projection images of the breast.

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.

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°.

TU-E-217BCD-03: Characterization of Anatomical Noise in Mammography, Tomosynthesis and Breast CT

PURPOSE

The role of breast density in cancer detection has been well characterized, and newer modalities such as breast tomosynthesis and breast CT (bCT) were developed to improve cancer detection in the dense breast by reducing anatomical complexity. Anatomical noise was characterized on a small cohort of patient images and compared across digital mammography, tomosynthesis, and bCT images.

METHODS AND MATERIALS

An IRB-approved and HIPPA-compliant clinical study was performed on women undergoing breast biopsy, and mammography, tomosynthesis, and bCT were performed on both breasts immediately prior to biopsy. A total of 23 women participated in this study, and the unaffected breast (no lesion) was evaluated. A total of 1000 regions of interest were sampled on each image data set, and the 2D noise power spectrum (NPS) was evaluated. This was radially averaged to produce a 1D NPS, and the NPS was fit to a power law: ln{NPS(f)} = alpha+betaxln(f), over an anatomically-relevant range of spatial frequencies. The slope, beta, was averaged across patients and compared between modalities and projections.

RESULTS

The value of beta was determined for bCT 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, beta was 3.06 (0.361) and 3.10 (0.315) for the CC and 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 beta for breast CT were significantly different than those for tomosynthesis and mammography (p<0.001, all 12 comparisons).

CONCLUSIONS

The results of this investigation demonstrate that the anatomical complexity of the breast, as characterized by the parameter beta, is statistically similar between mammography and tomosynthesis, a somewhat surprising finding. The breast CT image data, however, demonstrate a statistically-significant reduction in beta across all projections. Funded in part by Hologic Corporation and by a grant from the National Institute of Biomedical Imaging and Bioengineering, EB002138.

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.