Purpose:
Monte Carlo methods were used to generate lightly filtered high resolution x-ray spectra spanning from 20 kV to 640 kV.
Methods:
X-ray spectra were simulated for a conventional tungsten ...anode. The Monte Carlo N-Particle eXtended radiation transport code (MCNPX 2.6.0) was used to produce 35 spectra over the tube potential range from 20 kV to 640 kV, and cubic spline interpolation procedures were used to create piecewise polynomials characterizing the photon fluence per energy bin as a function of x-ray tube potential. Using these basis spectra and the cubic spline interpolation, 621 spectra were generated at 1 kV intervals from 20 to 640 kV. The tungsten anode spectral model using interpolating cubic splines (TASMICS) produces minimally filtered (0.8 mm Be) x-ray spectra with 1 keV energy resolution. The TASMICS spectra were compared mathematically with other, previously reported spectra.
Results:
Using pairedt-test analyses, no statistically significant difference (i.e., p > 0.05) was observed between compared spectra over energy bins above 1% of peak bremsstrahlung fluence. For all energy bins, the correlation of determination (R
2) demonstrated good correlation for all spectral comparisons. The mean overall difference (MOD) and mean absolute difference (MAD) were computed over energy bins (above 1% of peak bremsstrahlung fluence) and over all the kV permutations compared. MOD and MAD comparisons with previously reported spectra were 2.7% and 9.7%, respectively (TASMIP), 0.1% and 12.0%, respectively R. Birch and M. Marshall, “Computation of bremsstrahlung x-ray spectra and comparison with spectra measured with a Ge(Li) detector,” Phys. Med. Biol.
24, 505–517 (1979), 0.4% and 8.1%, respectively (Poludniowski), and 0.4% and 8.1%, respectively (AAPM TG 195). The effective energy of TASMICS spectra with 2.5 mm of added Al filtration ranged from 17 keV (at 20 kV) to 138 keV (at 640 kV); with 0.2 mm of added Cu filtration the effective energy was 9 keV at 20 kV and 169 keV at 640 kV.
Conclusions:
Ranging from 20 kV to 640 kV, 621 x-ray spectra were produced and are available at 1 kV tube potential intervals. The spectra are tabulated at 1 keV intervals. TASMICS spectra were shown to be largely equivalent to published spectral models and are available in spreadsheet format for interested users by emailing the corresponding author (JMB).
Purpose:
Current dosimetry methods in mammography assume that the breast is comprised of a homogeneous mixture of glandular and adipose tissues. Three‐dimensional (3D) dedicated breast CT (bCT) data ...sets were used previously to assess the complex anatomical structure within the breast, characterizing the statistical distribution of glandular tissue in the breast. The purpose of this work was to investigate the effect of bCT‐derived heterogeneous glandular distributions on dosimetry in mammography.
Methods:
bCT‐derived breast diameters, volumes, and 3D fibroglandular distributions were used to design realistic compressed breast models comprised of heterogeneous distributions of glandular tissue. The bCT‐derived glandular distributions were fit to biGaussian functions and used as probability density maps to assign the density distributions within compressed breast models. The MCNPX 2.6.0 Monte Carlo code was used to estimate monoenergetic normalized mean glandular dose “DgN(E)” values in mammography geometry. The DgN(E) values were then weighted by typical mammography x‐ray spectra to determine polyenergetic DgN (pDgN) coefficients for heterogeneous (pDgNhetero) and homogeneous (pDgNhomo) cases. The dependence of estimated pDgN values on phantom size, volumetric glandular fraction (VGF), x‐ray technique factors, and location of the heterogeneous glandular distributions was investigated.
Results:
The pDgNhetero coefficients were on average 35.3% (SD, 4.1) and 24.2% (SD, 3.0) lower than the pDgNhomo coefficients for the Mo–Mo and W–Rh x‐ray spectra, respectively, across all phantom sizes and VGFs when the glandular distributions were centered within the breast phantom in the coronal plane. At constant breast size, increasing VGF from 7.3% to 19.1% lead to a reduction in pDgNhetero relative to pDgNhomo of 23.6%–27.4% for a W–Rh spectrum. Displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width in the superior and inferior directions, resulted in a 37.3% and a −26.6% change in the pDgNhetero coefficient, respectively, relative to the centered distribution for the Mo–Mo spectrum. Lateral displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width, resulted in a 1.5% change in the pDgNhetero coefficient relative to the centered distribution for the W–Rh spectrum.
Conclusions:
Introducing bCT‐derived heterogeneous glandular distributions into mammography phantom design resulted in decreased glandular dose relative to the widely used homogeneous assumption. A homogeneous distribution overestimates the amount of glandular tissue near the entrant surface of the breast, where dose deposition is exponentially higher. While these findings are based on clinically measured distributions of glandular tissue using a large cohort of women, future work is required to improve the classification of glandular distributions based on breast size and overall glandular fraction.
The purpose of this study is to determine whether Hounsfield numbers of soft tissues on unenhanced abdominal CT of the same patient vary on repeat scans done on two different manufacturers' MDCT ...scanners.
A database search was performed to identify patients older than 18 years who underwent unenhanced CT of the abdomen and pelvis performed both on a Volume CT (GE Healthcare) and a Definition AS Plus (Siemens Healthcare) 64-MDCT scanner within 12 months of each other. After excluding those patients for whom Hounsfield unit measurements would be affected by mitigating factors, 48 patients (mean age, 58.8 years) were identified. Hounsfield unit measurements were obtained in nine different soft-tissue anatomic locations on each scan, and the location of these sites was kept identical on each scan pair. Data were analyzed to evaluate Hounsfield unit differences between these scanners.
In general, there was a low consistency in the Hounsfield unit measurements for each of these sites on scans obtained by the two scanners, with the subcutaneous fat in the left posterolateral flank showing the lowest correlation (intraclass correlation coefficient, 0.198). There were differences in the Hounsfield unit measurements obtained in all anatomic sites on scans obtained by both scanners. Mean Hounsfield unit measurements obtained on the Definition AS Plus scanner were lower than those obtained on the Volume CT scanner, with the intriguing exception of the anterior midline subcutaneous fat Hounsfield unit measurements, which were higher on the Definition AS Plus scanner. All differences were statistically significant (p < 0.05).
Hounsfield unit measurements for unenhanced abdominal soft tissues of the same patient vary between scanners of two common MDCT manufacturers.
Purpose:
Dosimetry in computed tomography (CT) is increasingly based on Monte Carlo studies that define the dose in the patient (in mGy) as a function of air kerma (free in air) at isocenter (mGy). ...The accuracy of Monte Carlo studies depends in part on the accuracy of the characterization of the bow tie filter for a given CT scanner model. A simple method for characterizing the bow tie filter attenuation profile in CT scanners would therefore be very useful. The theory behind such a method is proposed.
Methods:
A measurement protocol is discussed mathematically and demonstrated using computer simulation. The proposed method requires the placement of a radiation monitor at the periphery of the CT field, and the time domain signal (kerma rate versus time) is measured with good temporal resolution (
∼
200
Hz
or better) and with all other objects (e.g., patient couch) retracted from the field of view. Knowledge of the source to isocenter distance (or alternately, the isocenter to probe distance) is required. The stationary detector records the kerma rate versus time signal as the gantry rotates through several revolutions. From this temporal data, signal processing techniques are used to extract in-phase peaks, as well as out-of-phase kerma rate levels. From these data, the distance from isocenter to the probe can be determined (or, alternatively, the source to isocenter distance), and the angle-dependent bow tie filter attenuation can be computed. By measuring the angle-dependent bow tie filter attenuation at several kVp settings, the bow tie composition versus fan angle can be computed using basis decomposition techniques.
Results:
The simulations illustrated that with 2% added noise in the kerma rate versus time signal, the attenuation properties of a hypothetical two component (aluminum and polymethyl methacrylate) bow tie filter could be determined
(
r
2
>
0.99
)
. Although the computed basis material thicknesses were not exactly equal to the actual thicknesses, their combined attenuation factors matched that of the actual filter across kVp’s to within an average of 0.057%.
Conclusions:
It is concluded that the proposed method may provide a simple noninvasive approach to characterizing the performance of bow tie filters in CT systems; however, experimental validation is necessary.
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.
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’ (R2 > 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.
The studies by Sahbaee et al were well performed, and the results are provocative, but the authors of this editorial suggest that there are limitations to all modeling studies and that the results ...should be considered as only the first chapter in a much longer story about the role of contrast agents on radiation dose at CT.