Grid and slot scan scatter reduction in mammography: comparison by using Monte Carlo techniques

The impact of X-ray scatter correction software on abdomen radiography in terms of image quality and radiation dose

INTRODUCTION

The conventional anti-scatter grid is widely used in X-ray radiography to reduce scattered X-rays, but it increases patient dose. Scatter-correction software offers a dose-reducing alternative by correcting for scattered X-rays without a physical grid. Grids and software correction are necessary to reduce scatter radiation and improve image quality especially for the large body parts. The scatter correction can be beneficial in situations where the use of grid is challenging. The implementation of grids and advanced software correction techniques is imperative to ensure that radiographic images maintain high levels of clarity, contrast, and resolution, and ultimately facilitating more accurate diagnoses. This study compares image quality and radiation dose for abdomen exams using scatter correction software and physical grids.

METHODS

An anthropomorphic phantom (abdomen) underwent imaging with varying fat and lean tissue layers and body mass index (BMI) configurations. Imaging parameters included 70 kVp tube voltage, 110 cm SID, and Automatic Exposure Control (AEC) both lateral and central chambers. AP abdomen X-ray projections were acquired with and without an anti-scatter grid, and scatter correction software was applied. Image quality was assessed using contrast to noise ratio (CNR) and signal to noise ratio (SNR) metrics. The tube current mAs was considered an exposure factor that affected radiation dose and was used to compare the VG software and physical grid. Radiation dose was measured using Dose Area Products (DAP). The effective dose was estimated using Monte Carlo simulation-PCXMC software. Paired t-tests were used to investigate the image quality difference between the Gridless and VG software, Gridless and PG, and VG software and PG approaches. For the DAP and effective dose, paired t-test was used to investigate the difference between VG software and PG.

RESULTS

Images acquired with a grid had the highest mean CNR (71.3 ± 32) compared to Gridless (50 ± 33.8) and scatter correction software (59.3 ± 37.9). The mean SNR of the grid images was (82.7.3 ± 38.9), which is 18% higher than the scatter correction software images (70.4 ± 36.7) and 29% higher than in the Gridless images (62.9.3 ± 34). The mean DAP value was reduced by 81% when the scatter correction software was used compared to the grid (mean: 65.4 μGy.m2 and 338.2 μGy.m2, respectively) with a significant difference (p = 0.001). Scatter correction software resulted in a lower effective dose compared to physical grid use, (mean difference± SD = -0.3 ± 0.18 mSv) with a significant difference (P = 0.02).

CONCLUSION

Scatter correction software reduced the radiation dose required but images employing a grid yielded higher CNR and SNR. However, the radiation dose reduction might affect the image quality to a level that impacts the diagnostic information available. Thus, further research needs to be conducted to optimise the use of the scatter correction software.

IMPLICATION FOR PRACTICE

Objectively, X-ray scatter correction software might be promising in conditions where a grid cannot be applied.

Improving input contrast estimation to an x-ray imaging system

Objective.An appropriate parameter to study the performance of an x-ray imaging system is contrast transfer, or the system's ability to capture contrast in the radiation beam and bring it to the image. However, determining the input contrast to the system is difficult, as it is heavily affected by secondary radiation, which in turn depends on a multitude of factors. This work presents a method to improve the calculation of input contrast to the imaging system when PMMA phantoms are used.Approach.An expression to obtain input contrast from primary radiation attenuation and scatter-to-primary ratio is shown, and the approximation upon which it is based is discussed. Primary and secondary radiation emerging from the phantoms are calculated for monoenergetic pencil beams impinging on planar PMMA phantoms of different thicknesses. Monte Carlo simulations of two types of anti-scatter grids are also incorporated into the calculations.Main results.The primary and secondary components of radiant energy and grid transmission factors are presented for monoenergetic beams with energies from 10 to 150 keV. These results are then used to calculate input contrast for polyenergetic beams when using a commercial image quality phantom combined with different thicknesses of PMMA and anti-scatter grids.Significance.The information of the object contrast carried by the beam constitutes the input to the imaging system. An accurate determination of this input contrast can be carried out in a wide variety of situations from the study of a reduced number of cases, as those presented in this work for monoenergetic beams, PMMA phantoms of different thicknesses and anti-scatter grids. The relationship between the input contrast and the contrast due to primary radiation used in this work provides a good approximation for the different combinations of inserts, phantoms, grids, and energy spectra analyzed here.

A phantom study on dose efficiency for orthopedic applications: Comparing slot-scanning radiography using ultra-small-angle tomosynthesis to conventional radiography

PURPOSE

This paper studies the abilities of a twin-robotic x-ray slot-scanning system for orthopedic imaging to reduce dose by scatter rejection compared to conventional digital radiography.

METHODS

We investigate the dose saving capabilities, especially in terms of the signal- and the contrast-to-noise ratio, as well as the scatter-to-primary ratio of the proposed slot-scanning method in comparison to the state-of-the-art method for length-extended imaging. As a baseline, we use x-ray parameters of two clinically established acquisition protocols that provide the same detector entrance dose but are profoundly different in patient dose. To obtain an estimate of the photon-related noise directly from an x-ray image, we implement a Poisson-Gaussian noise model. This model is used to compare the dose efficiency of two settings and combined with the well-known K SNR to determine the transmission parameters. We present a method with an associated measurement protocol, utilizing the robotic capabilities of the used system to automatically obtain quasi-scatter-free ground-truth data with exact geometric correspondence to full-field and slot acquisitions. In total, we investigate two body regions (thoracic spine and lumbar spine) in anterior-posterior view with two patient sizes (BMI = 22 and 30) in two acquisition modes (conventional and slot scan with a flat-panel detector) with and without anti-scatter grid using an anthropomorphic upper-body phantom.

RESULTS

We have shown that it is feasible to combine the proposed approach with the K SNR for the determination of scatter rejection parameters. The use of an anti-scatter grid is indicated for full-field acquisitions allowing for dose savings up to 46% compared to their gridless counterparts. When changing the acquisition mode to the investigated slot scan, the use of an anti-scatter grid has no major impact on the image quality in terms of dose efficiency, in particular for patients with a BMI of 22. However, an increased contrast improvement factor was found. For normal-sized patients, up to 53% of dose can be saved additionally in comparison to full-field acquisitions with grid. Moreover, we could demonstrate that a slot size of 5 cm and air gap of 10 cm is sufficient to achieve scatter-to-primary ratios, which are equal or better compared to those of the full-field acquisitions with a grid.

CONCLUSIONS

We have shown, that the slot-scanning approach is always superior to the conventional full-field acquisition in terms of signal-to-noise and scatter-to-primary ratios. Compared to the state-of-the-art acquisition protocols with a grid, dose savings up to 53% are possible due to the scatter rejection without compromising the SNR. Hence, the use of the slot-scanning method is indicated, especially when it comes to regularly carried-out follow-up acquisitions, for example, in the case of scoliosis monitoring.

A fluence modulation and scatter shielding apparatus for dedicated breast CT: Theory of operation

PURPOSE

To introduce an auxiliary apparatus of fluence modulation and scatter shielding for dedicated breast computed tomography (bCT) and the corresponding patient-specific method of image acquisition.

METHODS

The apparatus is composed of two assemblies, referred to herein as the "Dynamic Fluence Gate" (FG) and "Scatter Shield" (SS), that work in synchrony to form a narrow beam sweeping the entire fan angle coverage of the imaging system during a projection. The apparatus, as a whole, is referred to as FG-SS. FG and SS are pre-patient and post-patient units, respectively. Each is composed of a rotating drum, on top of which are installed two sheets of high x-ray attenuating material, a sensory system, and the constituent robotics. The sheets of each unit are positioned such that an opening - a window Fluence Modulation and Scatter Shielding is formed between them. The rotations of the drums and positioning of the sheets are synchronized and adjusted such that a line of sight is created between the source, FG window, the breast, and the SS window. With line of sight achieved, the narrow beam transitions from the source to the detector. The fluence of the narrow beam during a projection depends on the size, shape, and positioning of the breast. The FG-SS method of imaging is discussed mathematically and demonstrated using computer simulations. A series of Monte Carlo simulations are conducted to evaluate the performance of the system as relates to its impact on the imager's dynamic range, dose distribution to the breast, noise inhomogeneity in reconstructed images, and scatter buildup in projections within small, medium, and large breasts composed of homogeneous medium breast tissue.

RESULTS

Implementation of FG-SS results in near scatter-free projections, reduction in both dose and the required dynamic range of the imager, and equalization of the quantum noise distribution in the reconstructed image. Using the disclosed design, the dynamic range was reduced by factors ranging from 1.6 to 5.5 across the range of breast sizes studied. A reduction in the acquisition of the scattered rays, varying between the factors of 6.1 (in the small breast) and 9.8 (in that large breast) was achieved and consequently, shading artifacts were suppressed. Noise in the CT image was equalized by reducing the overall spatial variation from 29% to <5% in small breast and from 45% to 14% in the large breast. An overall reduction in deposited dose to the breast was achieved - between 26% and 39% depending on the breast size.

CONCLUSIONS

Utilization of the FG-SS apparatus and technique was demonstrated via simulations to result in: significant reductions in dose to the breast, reductions in scatter uptake in projections, reduced required dynamic range of the imager, and homogenizing of quantum noise throughout the reconstructed image.

Physics-based spectral compensation algorithm for x-ray CT with primary modulator

X-ray computed tomography (CT) scatter correction using primary modulator has been continuously developed over the past years, with progress in improving the performance of scatter correction. In this work, we further advance the primary modulator technique towards practical applications where the spectral nonuniformity caused by the modulator continues to be a challenging problem. A physics-based spectral compensation algorithm is proposed to adaptively correct for the spectral nonuniformity, and hence to reduce the resultant ring artifacts on reconstructed CT images. First, an initial spectrum of the CT system without the primary modulator is modeled using an understanding of x-ray CT physics, and optimized by an expectation maximization method; then, the optimized estimation of the initial spectrum is utilized to adaptively calculate the effective modulator thickness from measured transmissions of the primary modulator at each detector element, leading to a set of new spectra that is able to capture the nonuniform spectral distribution of the primary modulator; finally, using the modulator-modeled spectrum, a beam hardening mapping function is generated and beam hardening correction is applied to CT projections. A CatPhan600 phantom and an anthropomorphic thorax phantom were scanned with three different primary modulators to evaluate the approach. For the Catphan phantom, the spectral compensation algorithm efficiently removes the ring (and band) artifacts that otherwise dominate the reconstructed CT image. For the three modulators with nominal copper thickness of 52.5, 105 and 210 [Formula: see text]m, our method reduces the CT number nonuniformity from 147.9, 436.2 and 696.4 Hounsfield units (HU) to 14.6, 26.2 and 13.6 HU, respectively, close to that of the reference image (i.e. 7.5 HU). For the thorax phantom, the ring artifacts are also suppressed significantly on the transaxial image; on the sagittal image, the alternating black-and-white patterns are largely removed, with the CT number nonuniformity being reduced from 282.0 HU to 38.5 HU.

Feasibility of estimating volumetric breast density from mammographic x-ray spectra using a cadmium telluride photon-counting detector

PURPOSE

Mammographic density of glandular breast tissue has a masking effect that can reduce lesion detection accuracy and is also a strong risk factor for breast cancer. Therefore, accurate quantitative estimation of breast density is clinically important. In this study, we investigate experimentally the feasibility of quantifying volumetric breast density with spectral mammography using a CdTe-based photon-counting detector.

METHODS

To demonstrate proof-of-principle, this study was carried out using the single pixel Amptek XR-100T-CdTe detector. The total number of x rays recorded by the detector from a single pencil-beam projection through 50%/50% of adipose/glandular mass fraction-equivalent phantoms was measured. Material decomposition assuming two, four, and eight energy bins was then applied to characterize the inspected phantom into adipose and glandular using log-likelihood estimation, taking into account the polychromatic source, the detector response function, and the energy-dependent attenuation.

RESULTS

Measurement tests were carried out for different doses, kVp settings, and different breast sizes. For dose of 1 mGy and above, the percent relative root mean square (RMS) errors of the estimated breast density was measured below 7% for all three phantom studies. It was also observed that some decrease in RMS errors was achieved using eight energy bins. For 3 and 4 cm thick phantoms, performance at 40 and 45 kVp showed similar performance. However, it was observed that 45 kVp showed better performance for a phantom thickness of 6 cm at low dose levels due to increased statistical variation at lower photon count levels with 40 kVp.

CONCLUSION

The results of the current study suggest that photon-counting spectral mammography systems using CdTe detectors have the potential to be used for accurate quantification of volumetric breast density on a pixel-to-pixel basis, with an RMS error of less than 7%.

Virtual scatter modulation for X-ray CT scatter correction using primary modulator

A new scatter estimation algorithm with a concept of virtual scatter modulation for X-ray scatter correction using primary modulator is proposed to reduce the aliasing errors in the estimated scatter. Virtual scatter modulation can be realized through dividing the measured primary-modulated image by the measured modulation function. After the division, the aggravation of the aliasing of primary due to the non-uniformity of the modulation function is largely transferred to that of scatter. Since scatter in general has less high frequencies than primary does, the aggravation of its aliasing is expected to be weaker, and therefore the overall aliasing can be reduced. A CatPhan©600 phantom and an anthropomorphic thorax phantom are scanned on a tabletop X-ray cone-beam computed tomography system to validate our proposed algorithm. On the Catphan phantom, the oscillations that are clearly observed in the central region of the Catphan scatter profile estimated using the original primary-modulation algorithm, are mostly eliminated with the proposed scatter modulation algorithm, leading to less residual artifacts and better CT number uniformity in the reconstructed image. Compared with 38.9 HU of CT nonuniformity in a selected uniform region when the primary-modulation algorithm is used, the new algorithm significantly reduces it to 4.5 HU, reaching the same level of uniformity as the ground truth reference. On the thorax phantom, overall better CT number uniformity is also achieved.

Design and evaluation of a grid reciprocation scheme for use in digital breast tomosynthesis

This work describes a methodology for efficient removal of scatter radiation during digital breast tomosynthesis (DBT). The goal of this approach is to enable grid image obscuration without a large increase in radiation dose by minimizing misalignment of the grid focal point (GFP) and x-ray focal spot (XFS) during grid reciprocation. Hardware for the motion scheme was built and tested on the dual modality breast tomosynthesis (DMT) scanner, which combines DBT and molecular breast tomosynthesis (MBT) on a single gantry. The DMT scanner uses fully isocentric rotation of tube and x-ray detector for maintaining a fixed tube-detector alignment during DBT imaging. A cellular focused copper prototype grid with 80 cm focal length, 3.85 mm height, 0.1 mm thick lamellae, and 1.1 mm hole pitch was tested. Primary transmission of the grid at 28 kV tube voltage was on average 74% with the grid stationary and aligned for maximum transmission. It fell to 72% during grid reciprocation by the proposed method. Residual grid line artifacts (GLAs) in projection views and reconstructed DBT images are characterized and methods for reducing the visibility of GLAs in the reconstructed volume through projection image flat-field correction and spatial frequency-based filtering of the DBT slices are described and evaluated. The software correction methods reduce the visibility of these artifacts in the reconstructed volume, making them imperceptible both in the reconstructed DBT images and their Fourier transforms.

On image quality metrics and the usefulness of grids in digital mammography

Antiscatter grids are used in digital mammography to reduce the scattered radiation from the breast and improve image contrast. They are, however, imperfect and lead to partial absorption of primary radiation, as well as failing to absorb all scattered radiation. Nevertheless, the general consensus has been that antiscatter grids improve image quality for the majority of breast types and sizes. There is, however, inconsistency in the literature, and recent results show that a substantial image quality improvement can be achieved even for thick breasts if the grid is disposed of. The purpose of this study was to investigate if differences in the considered imaging task and experimental setup could explain the different outcomes. We estimated the dose reduction that can be achieved if the grid were to be removed as a function of breast thickness with varying geometries and experimental conditions. Image quality was quantified by the signal-difference-to-noise ratio (SDNR) measured using an aluminum (Al) filter on blocks of poly(methyl methacrylate) (PMMA), and images were acquired with and without grid at a constant exposure. We also used a theoretical model validated with Monte Carlo simulations. Both theoretically and experimentally, the main finding was that when a large [Formula: see text] Al filter was used, the SDNR values for the gridless images were overestimated up to 25% compared to the values for the small [Formula: see text] filter, and gridless imaging was superior for any PMMA thickness. For the small Al filter, gridless imaging was only superior for PMMAs thinner than 4 cm. This discrepancy can be explained by a different sensitivity to and sampling of the angular scatter spread function, depending on the size of the contrast object. The experimental differences were eliminated either by using a smaller region of interest close to the edge of the large filter or by applying a technique of scatter correction by subtracting the estimated scatter image. These results explain the different conclusions reported in the literature and show the importance of the selection of measurement methods. Since the interesting structures in mammography are below the 1-cm scale, we advocate the use of smaller contrast objects for assessment of antiscatter grid performance.

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.

Characterization of scatter in digital mammography from use of Monte Carlo simulations and comparison to physical measurements

PURPOSE

Monte Carlo simulations were performed with the goal of verifying previously published physical measurements characterizing scatter as a function of apparent thickness. A secondary goal was to provide a way of determining what effect tissue glandularity might have on the scatter characteristics of breast tissue. The overall reason for characterizing mammography scatter in this research is the application of these data to an image processing-based scatter-correction program.

METHODS

mcnpx was used to simulate scatter from an infinitesimal pencil beam using typical mammography geometries and techniques. The spreading of the pencil beam was characterized by two parameters: mean radial extent (MRE) and scatter fraction (SF). The SF and MRE were found as functions of target, filter, tube potential, phantom thickness, and the presence or absence of a grid. The SF was determined by separating scatter and primary by the angle of incidence on the detector, then finding the ratio of the measured scatter to the total number of detected events. The accuracy of the MRE was determined by placing ring-shaped tallies around the impulse and fitting those data to the point-spread function (PSF) equation using the value for MRE derived from the physical measurements. The goodness-of-fit was determined for each data set as a means of assessing the accuracy of the physical MRE data. The effect of breast glandularity on the SF, MRE, and apparent tissue thickness was also considered for a limited number of techniques.

RESULTS

The agreement between the physical measurements and the results of the Monte Carlo simulations was assessed. With a grid, the SFs ranged from 0.065 to 0.089, with absolute differences between the measured and simulated SFs averaging 0.02. Without a grid, the range was 0.28-0.51, with absolute differences averaging -0.01. The goodness-of-fit values comparing the Monte Carlo data to the PSF from the physical measurements ranged from 0.96 to 1.00 with a grid and 0.65 to 0.86 without a grid. Analysis of the data suggested that the nongrid data could be better described by a biexponential function than the single exponential used here. The simulations assessing the effect of breast composition on SF and MRE showed only a slight impact on these quantities. When compared to a mix of 50% glandular/50% adipose tissue, the impact of substituting adipose or glandular breast compositions on the apparent thickness of the tissue was about 5%.

CONCLUSIONS

The findings show agreement between the physical measurements published previously and the Monte Carlo simulations presented here; the resulting data can therefore be used more confidently for an application such as image processing-based scatter correction. The findings also suggest that breast composition does not have a major impact on the scatter characteristics of breast tissue. Application of the scatter data to the development of a scatter-correction software program can be simplified by ignoring the variations in density among breast tissues.

Estimation of scattered radiation in digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is a promising technique to overcome the tissue superposition limitations found in planar 2D x-ray mammography. However, as most DBT systems do not employ an anti-scatter grid, the levels of scattered radiation recorded within the image receptor are significantly higher than that observed in planar 2D x-ray mammography. Knowledge of this field is necessary as part of any correction scheme and for computer modelling and optimisation of this examination. Monte Carlo (MC) simulations are often used for this purpose, however they are computationally expensive and a more rapid method of calculation is desirable. This issue is addressed in this work by the development of a fast kernel-based methodology for scatter field estimation using a detailed realistic DBT geometry. Thickness-dependent scatter kernels, which were validated against the literature with a maximum discrepancy of 4% for an idealised geometry, have been calculated and a new physical parameter (air gap distance) was used to estimate more accurately the distribution of scattered radiation for a series of anthropomorphic breast phantom models. The proposed methodology considers, for the first time, the effects of scattered radiation from the compression paddle and breast support plate, which can represent more than 30% of the total scattered radiation recorded within the image receptor. The results show that the scatter field estimator can calculate scattered radiation images in an average of 80 min for projection angles up to 25° with equal to or less than a 10% error across most of the breast area when compared with direct MC simulations.

Characterization of scatter in digital mammography from physical measurements

PURPOSE

That scattered radiation negatively impacts the quality of medical radiographic imaging is well known. In mammography, even slight amounts of scatter reduce the high contrast required for subtle soft-tissue imaging. In current clinical mammography, image contrast is partially improved by use of an antiscatter grid. This form of scatter rejection comes with a sizeable dose penalty related to the concomitant elimination of valuable primary radiation. Digital mammography allows the use of image processing as a method of scatter correction that might avoid effects that negatively impact primary radiation, while potentially providing more contrast improvement than is currently possible with a grid. For this approach to be feasible, a detailed characterization of the scatter is needed. Previous research has modeled scatter as a constant background that serves as a DC bias across the imaging surface. The goal of this study was to provide a more substantive data set for characterizing the spatially-variant features of scatter radiation at the image detector of modern mammography units.

METHODS

This data set was acquired from a model of the radiation beam as a matrix of very narrow rays or pencil beams. As each pencil beam penetrates tissue, the pencil widens in a predictable manner due to the production of scatter. The resultant spreading of the pencil beam at the detector surface can be characterized by two parameters: mean radial extent (MRE) and scatter fraction (SF). The SF and MRE were calculated from measurements obtained using the beam stop method. Two digital mammography units were utilized, and the SF and MRE were found as functions of target, filter, tube potential, phantom thickness, and presence or absence of a grid. These values were then used to generate general equations allowing the SF and MRE to be calculated for any combination of the above parameters.

RESULTS

With a grid, the SF ranged from a minimum of about 0.05 to a maximum of about 0.16, and the MRE ranged from about 3 to 13 mm. Without a grid, the SF ranged from a minimum of 0.25 to a maximum of 0.52, and the MRE ranged from about 20 to 45 mm. The SF with a grid demonstrated a mild dependence on target/filter combination and kV, whereas the SF without a grid was independent of these factors. The MRE demonstrated a complex relationship as a function of kV, with notable difference among target/filter combinations. The primary source of change in both the SF and MRE was phantom thickness.

CONCLUSIONS

Because breast tissue varies spatially in physical density and elemental content, the effective thickness of breast tissue varies spatially across the imaging field, resulting in a spatially-variant scatter distribution in the imaging field. The data generated in this study can be used to characterize the scatter contribution on a point-by-point basis, for a variety of different techniques.

Evaluation of scatter-to-primary ratio, grid performance and normalized average glandular dose in mammography by Monte Carlo simulation including interference and energy broadening effects

In this work, a computational code for the study of imaging systems and dosimetry in conventional and digital mammography through Monte Carlo simulations is described. The developed code includes interference and Doppler energy broadening for simulation of elastic and inelastic photon scattering, respectively. The code estimates the contribution of scattered radiation to image quality through the spatial distribution of the scatter-to-primary ratio (S/P). It allows the inclusion of different designs of anti-scatter grids (linear or cellular), for evaluation of contrast improvement factor (CIF), Bucky factor (BF) and signal difference-to-noise ratio improvement factor (SIF). It also allows the computation of the normalized average glandular dose, D(g).(N). These quantities were studied for different breast thicknesses and compositions, anode/filter combinations and tube potentials. Results showed that the S/P increases linearly with breast thickness, varying slightly with breast composition or the spectrum used. Evaluation of grid performance showed that the cellular grid provides the highest CIF with smaller BF. The SIF was also greater for the cellular grid, although both grids showed SIF < 1 for thin breasts. Results for D(g).(N) showed that it increases with the half-value layer (HVL) of the spectrum, decreases considerably with breast thickness and has a small dependence on the anode/filter combination. Inclusion of interference effects of breast tissues affected the values of S/P obtained with the grid up to 25%, while the energy broadening effect produced smaller variations on the evaluated quantities.

Scatter correction method for x-ray CT using primary modulation: phantom studies

PURPOSE

Scatter correction is a major challenge in x-ray imaging using large area detectors. Recently, the authors proposed a promising scatter correction method for x-ray computed tomography (CT) using primary modulation. Proof of concept was previously illustrated by Monte Carlo simulations and physical experiments on a small phantom with a simple geometry. In this work, the authors provide a quantitative evaluation of the primary modulation technique and demonstrate its performance in applications where scatter correction is more challenging.

METHODS

The authors first analyze the potential errors of the estimated scatter in the primary modulation method. On two tabletop CT systems, the method is investigated using three phantoms: A Catphan 600 phantom, an anthropomorphic chest phantom, and the Catphan 600 phantom with two annuli. Two different primary modulators are also designed to show the impact of the modulator parameters on the scatter correction efficiency. The first is an aluminum modulator with a weak modulation and a low modulation frequency, and the second is a copper modulator with a strong modulation and a high modulation frequency.

RESULTS

On the Catphan 600 phantom in the first study, the method reduces the error of the CT number in the selected regions of interest (ROIs) from 371.4 to 21.9 Hounsfield units (HU); the contrast to noise ratio also increases from 10.9 to 19.2. On the anthropomorphic chest phantom in the second study, which represents a more difficult case due to the high scatter signals and object heterogeneity, the method reduces the error of the CT number from 327 to 19 HU in the selected ROIs and from 31.4% to 5.7% on the overall average. The third study is to investigate the impact of object size on the efficiency of our method. The scatter-to-primary ratio estimation error on the Catphan 600 phantom without any annulus (20 cm in diameter) is at the level of 0.04, it rises to 0.07 and 0.1 on the phantom with an elliptical annulus (30 cm in the minor axis and 38 cm in the major axis) and with a circular annulus (38 cm in diameter).

CONCLUSIONS

On the three phantom studies, good scatter correction performance of the proposed method has been demonstrated using both image comparisons and quantitative analysis. The theory and experiments demonstrate that a strong primary modulation that possesses a low transmission factor and a high modulation frequency is preferred for high scatter correction accuracy.

Development of a multi-spectral, multi-geometry computational model for X-ray breast imaging

The introduction of novel applications in X-ray breast imaging warrants new research for image acquisition optimisation. A simulation model was developed to investigate the influence of different imaging techniques and acquisition parameters. It was modelled in Monte Carlo N-Particle Extended and contains an X-ray tube with photon production, a breast model and antiscatter grid model. This paper describes the simulation model, compares the results with experimental and literature data and presents the influence of breast and antiscatter grid parameters on scatter radiation.

Evaluation of scatter effects on image quality for breast tomosynthesis

Digital breast tomosynthesis uses a limited number (typically 10-20) of low-dose x-ray projections to produce a pseudo-three-dimensional volume tomographic reconstruction of the breast. The purpose of this investigation was to characterize and evaluate the effect of scattered radiation on the image quality for breast tomosynthesis. In a simulation, scatter point spread functions generated by a Monte Carlo simulation method were convolved over the breast projection to estimate the distribution of scatter for each angle of tomosynthesis projection. The results demonstrate that in the absence of scatter reduction techniques, images will be affected by cupping artifacts, and there will be reduced accuracy of attenuation values inferred from the reconstructed images. The effect of x-ray scatter on the contrast, noise, and lesion signal-difference-to-noise ratio (SDNR) in tomosynthesis reconstruction was measured as a function of the tumor size. When a with-scatter reconstruction was compared to one without scatter for a 5 cm compressed breast, the following results were observed. The contrast in the reconstructed central slice image of a tumorlike mass (14 mm in diameter) was reduced by 30%, the voxel value (inferred attenuation coefficient) was reduced by 28%, and the SDNR fell by 60%. The authors have quantified the degree to which scatter degrades the image quality over a wide range of parameters relevant to breast tomosynthesis, including x-ray beam energy, breast thickness, breast diameter, and breast composition. They also demonstrate, though, that even without a scatter rejection device, the contrast and SDNR in the reconstructed tomosynthesis slice are higher than those of conventional mammographic projection images acquired with a grid at an equivalent total exposure.

The effect of scatter and glare on image quality in contrast-enhanced breast imaging using an a-Si/CsI(TI) full-field flat panel detector

The purpose of this study is to evaluate the performance of an antiscatter grid and its potential benefit on image quality for a full-field digital mammography (FFDM) detector geometry at energies typical for temporal subtraction contrast-enhanced (CE) breast imaging. The signal intensities from primary, scatter, and glare were quantified in images acquired with an a-Si/CsI(T1) FFDM detector using a Rh target and a 0.27 mm Cu filter at tube voltages ranging from 35 to 49 kV. Measurements were obtained at the center of the irradiation region of 20-80 mm thick breast-equivalent phantoms. The phantoms were imaged with and without an antiscatter grid. Based on these data, the performance of the antiscatter grid was determined by calculating the primary and scatter transmission factors (T(P) and T(S)) and Bucky factors (Bf). In addition, glare-to-primary ratios (GPRs) and scatter-to-primary ratios (SPRs) were quantified. The effect of the antiscatter grid on the signal-difference-to-noise ratio (SDNR) was also assessed. It was found that T(P) increases with kV but does not depend on the phantom thickness; T(P) values between 0.81 and 0.84 were measured. T(S) increases with kV and phantom thickness; T(S) values between 0.13 and 0.21 were measured. Bf decreases with kV and increases with phantom thickness; Bf ranges from 1.4 to 2.1. GPR is nearly constant, varying from 0.10 to 0.11. SPR without an antiscatter grid (SPR-) ranges from 0.35 to 1.34. SPR- decreases by approximately 9% from 35 to 49 kV for a given phantom thickness and is 3.5 times larger for an 80 mm thick breast-equivalent phantom than for a 20 mm thick breast-equivalent phantom. SPR with an antiscatter grid (SPR+) ranges from 0.06 to 0.31. SPR+ increases by approximately 23% from 35 to 49 kV for a given phantom thickness; SPR+ is four times larger for an 80 mm breast-equivalent phantom than for a 20 mm breast-equivalent phantom. When imaging a 25 mm PMMA plate at the same mean glandular dose with and without an antiscatter grid, the SDNR is 4% greater with a grid than without. For an 75 mm PMMA plate, the SDNR is 20% greater with a grid. In conclusion, at the higher x-ray energy range used for CE-DM and CE-DBT, an antiscatter grid significantly reduces SPR and improves SDNR. These effects are most pronounced for thick breasts.

Quantification of breast density with dual energy mammography: a simulation study

Breast density, the percentage of glandular breast tissue, has been identified as an important yet underutilized risk factor in the development of breast cancer. A quantitative method to measure breast density with dual energy imaging was investigated using a computer simulation model. Two configurations to measure breast density were evaluated: the usage of monoenergetic beams and an ideal detector, and the usage of polyenergetic beams with spectra from a tungsten anode x-ray tube with a detector modeled after a digital mammography system. The simulation model calculated the mean glandular dose necessary to quantify the variability of breast density to within 1/3%. The breast was modeled as a semicircle 10 cm in radius with equal homogenous thicknesses of adipose and glandular tissues. Breast thicknesses were considered in the range of 2-10 cm and energies in the range of 10-150 keV for the two monoenergetic beams, and 20-150 kVp for spectra with a tungsten anode x-ray tube. For a 4.2 cm breast thickness, the required mean glandular doses were 0.183 microGy for two monoenergetic beams at 19 and 71 keV, and 9.85 microGy for two polyenergetic spectra from a tungsten anode at 32 and 96 kVp with beam filtrations of 50 microm Rh and 300 microm Cu for the low and high energy beams, respectively. The results suggest that for either configuration, breast density can be precisely measured with dual energy imaging requiring only a small amount of additional dose to the breast. The possibility of using a standard screening mammogram as the low energy image is also discussed.

Breast cancer imaging: a perspective for the next decade

Breast imaging is largely indicated for detection, diagnosis, and clinical management of breast cancer and for evaluation of the integrity of breast implants. In this work, a prospective view of techniques for breast cancer detection and diagnosis is provided based on an assessment of current trends. The potential role of emerging techniques that are under various stages of research and development is also addressed. It appears that the primary imaging tool for breast cancer screening in the next decade will be high-resolution, high-contrast, anatomical x-ray imaging with or without depth information. MRI and ultrasonography will have an increasingly important adjunctive role for imaging high-risk patients and women with dense breasts. Pilot studies with dedicated breast CT have demonstrated high-resolution three-dimensional imaging capabilities, but several technological barriers must be overcome before clinical adoption. Radionuclide based imaging techniques and x-ray imaging with intravenously injected contrast offer substantial potential as a diagnostic tools and for evaluation of suspicious lesions. Developing optical and electromagnetic imaging techniques hold significant potential for physiologic information and they are likely to be of most value when integrated with or adjunctively used with techniques that provide anatomic information. Experimental studies with breast specimens suggest that phase-sensitive x-ray imaging techniques can provide edge enhancement and contrast improvement but more research is needed to evaluate their potential role in clinical breast imaging. From the technological perspective, in addition to improvements within each modality, there is likely to be a trend towards multi-modality systems that combine anatomic with physiologic information. We are also likely to transition from a standardized screening, where all women undergo the same imaging exam (mammography), to selection of a screening modality or modalities based an individual-risk or other classification.

The effect of breast compression on mass conspicuity in digital mammography

This study analyzed how the inherent quality of diagnostic information in digital mammography could be affected by breast compression. A digital mammography system was modeled using a Monte Carlo algorithm based on the Penelope program, which has been successfully used to model several medical imaging systems. First, the Monte Carlo program was validated against previous measurements and simulations. Once validated, the Monte Carlo software modeled a digital mammography system by tracking photons through a voxelized software breast phantom, containing anatomical structures and breast masses, and following photons until they were absorbed by a selenium-based flat-panel detector. Simulations were performed for two compression conditions (standard compression and 12.5% reduced compression) and three photon flux conditions (constant flux, constant detector signal, and constant glandular dose). The results showed that reduced compression led to higher scatter fractions, as expected. For the constant photon flux condition, decreased compression also reduced glandular dose. For constant glandular dose, the SdNR for a 4 cm breast was 0.60 +/- 0.11 and 0.62 +/- 0.11 under standard and reduced compressions, respectively. For the 6 cm case with constant glandular dose, the SdNR was 0.50 +/- 0.11 and 0.49 +/- 0.10 under standard and reduced compressions, respectively. The results suggest that if a particular imaging system can handle an approximately 10% increase in total tube output and 10% decrease in detector signal, breast compression can be reduced by about 12% in terms of breast thickness with little impact on image quality or dose.

Comparison of slot scanning digital mammography system with full-field digital mammography system

The purpose of this study was to evaluate and compare microcalcification detectability of two commercial full-field digital mammography (DM) systems. The first unit was a flat panel based DM system (FFDM) which employed an anti-scatter grid method to reject scatter, and the second unit was a charge-coupled device-based DM system (SSDM) which used scanning slot imaging geometry to reduce scatter radiation. Both systems have comparable scatter-to-primary ratios. In this study, 125-160 and 200-250 microm calcium carbonate grains were used to simulate microcalcifications and imaged by both DM systems. The calcium carbonate grains were overlapped with a 5-cm-thick 50% adipose/50% glandular simulated breast tissue slab and an anthropomorphic breast phantom (RMI 165, Gammex) for imaging at two different mean glandular dose levels: 0.87 and 1.74 mGy. A reading study was conducted with seven board certified mammographers with images displayed on review workstations. A five-point confidence level rating was used to score each detection task. Receiver operating characteristic (ROC) analysis was performed and the area under the ROC curve (A(z)) was used to quantify and compare the performances of these two systems. The results showed that with the simulated breast tissue slab (uniform background), the SSDM system resulted in higher A(z)'s than the FFDM system at both MGD levels with the difference statistically significant at 0.87 mGy only. With the anthropomorphic breast phantom (tissue structure background), the SSDM system performed better than the FFDM system at 0.87 mGy but worse at 1.74 mGy. However, the differences were not found to be statistically significant.

Computed tomography with energy-resolved detection: a feasibility study

The feasibility of computed tomography (CT) with energy-resolved x-ray detection has been investigated. A breast CT design with multi slit multi slice (MSMS) data acquisition was used for this study. The MSMS CT includes linear arrays of photon counting detectors separated by gaps. This CT configuration allows for efficient scatter rejection and 3D data acquisition. The energy-resolved CT images were simulated using a digital breast phantom and the design parameters of the proposed MSMS CT. The phantom had 14 cm diameter and 50/50 adipose/glandular composition, and included carcinoma, adipose, blood, iodine and CaCO3 as contrast elements. The x-ray technique was 90 kVp tube voltage with 660 mR skin exposure. Photon counting, charge (energy) integrating and photon energy weighting CT images were generated. The contrast-to-noise (CNR) improvement with photon energy weighting was quantified. The dual energy subtracted images of CaCO3 and iodine were generated using a single CT scan at a fixed x-ray tube voltage. The x-ray spectrum was electronically split into low- and high-energy parts by a photon counting detector. The CNR of the energy weighting CT images of carcinoma, blood, adipose, iodine, and CaCO3 was higher by a factor of 1.16, 1.20, 1.21, 1.36 and 1.35, respectively, as compared to CT with a conventional charge (energy) integrating detector. Photon energy weighting was applied to CT projections prior to dual energy subtraction and reconstruction. Photon energy weighting improved the CNR in dual energy subtracted CT images of CaCO3 and iodine by a factor of 1.35 and 1.33, respectively. The combination of CNR improvements due to scatter rejection and energy weighting was in the range of 1.71-2 depending on the type of the contrast element. The tilted angle CZT detector was considered as the detector of choice. Experiments were performed to test the effect of the tilting angle on the energy spectrum. Using the CZT detector with 20 degrees tilting angle decreased the tailing of the measured x-ray spectrum as compared to a conventional CZT detector. It was concluded that the energy-resolved MSMS CT with tilted angle CZT detector is potentially feasible and could provide a unique combination of photon counting, energy weighting, scatter rejection and single kVp dual energy subtraction CT imaging.

Tomosynthesis and contrast-enhanced digital mammography: recent advances in digital mammography

Digital mammography is more and more replacing conventional mammography. Initial concerns about an inferior image quality of digital mammography have been largely overcome and recent studies even show digital mammography to be superior in women with dense breasts, while at the same time reducing radiation exposure. Nevertheless, an important limitation of digital mammography remains: namely, the fact that summation may obscure lesions in dense breast tissue. However, digital mammography offers the option of so-called advanced applications, and two of these, contrast-enhanced mammography and tomosynthesis, are promising candidates for improving the detection of breast lesions otherwise obscured by the summation of dense tissue. Two techniques of contrast-enhanced mammography are available: temporal subtraction of images acquired before and after contrast administration and the so-called dual-energy technique, which means that pairs of low/high-energy images acquired after contrast administration are subtracted. Tomosynthesis on the other hand provides three-dimensional information on the breast. The images are acquired with different angulations of the X-ray tube while the object or detector is static. Various reconstruction algorithms can then be applied to the set of typically nine to 28 source images to reconstruct 1-mm slices with a reduced risk of obscuring pathology. Combinations of both advanced applications have only been investigated in individual experimental studies; more advanced software algorithms and CAD systems are still in their infancy and have only undergone preliminary clinical evaluation.

Scatter correction method for X-ray CT using primary modulation: theory and preliminary results

An X-ray system with a large area detector has high scatter-to-primary ratios (SPRs), which result in severe artifacts in reconstructed computed tomography (CT) images. A scatter correction algorithm is introduced that provides effective scatter correction but does not require additional patient exposure. The key hypothesis of the algorithm is that the high-frequency components of the X-ray spatial distribution do not result in strong high-frequency signals in the scatter. A calibration sheet with a checkerboard pattern of semitransparent blockers (a "primary modulator") is inserted between the X-ray source and the object. The primary distribution is partially modulated by a high-frequency function, while the scatter distribution still has dominant low-frequency components, based on the hypothesis. Filtering and demodulation techniques suffice to extract the low-frequency components of the primary and hence obtain the scatter estimation. The hypothesis was validated using Monte Carlo (MC) simulation, and the algorithm was evaluated by both MC simulations and physical experiments. Reconstructions of a software humanoid phantom suggested system parameters in the physical implementation and showed that the proposed method reduced the relative mean square error of the reconstructed image in the central region of interest from 74.2% to below 1%. In preliminary physical experiments on the standard evaluation phantom, this error was reduced from 31.8% to 2.3%, and it was also demonstrated that the algorithm has no noticeable impact on the resolution of the reconstructed image in spite of the filter-based approach. Although the proposed scatter correction technique was implemented for X-ray CT, it can also be used in other X-ray imaging applications, as long as a primary modulator can be inserted between the X-ray source and the imaged object.

Tilted angle CZT detector for photon counting/energy weighting x-ray and CT imaging

X-ray imaging with a photon counting/energy weighting detector can provide the highest signal to noise ratio (SNR). Scanning slit/multi-slit x-ray image acquisition can provide a dose-efficient scatter rejection, which increases SNR. Use of a photon counting/energy weighting detector in a scanning slit/multi-slit acquisition geometry could provide highest possible dose efficiency in x-ray and CT imaging. Currently, the most advanced photon counting detector is the cadmium zinc telluride (CZT) detector, which, however, is suboptimal for energy resolved x-ray imaging. A tilted angle CZT detector is proposed in this work for applications in photon counting/energy weighting x-ray and CT imaging. In tilted angle configuration, the x-ray beam hits the surface of the linear array of CZT crystals at a small angle. This allows the use of CZT crystals of a small thickness while maintaining the high photon absorption. Small thickness CZT detectors allow for a significant decrease in the polarization effect in the CZT volume and an increase in count rate. The tilted angle CZT with a small thickness also provides higher spatial and energy resolution, and shorter charge collection time, which potentially enables fast energy resolving x-ray image acquisition. In this work, the major performance parameters of the tilted angle CZT detector, including its count rate, spatial resolution and energy resolution, were evaluated. It was shown that for a CZT detector with a 0.7 mm thickness and 13 degrees tilting angle, the maximum count rate can be increased by 10.7 times, while photon absorption remains >90% at photon energies up to 120 keV. Photon counting/energy weighting x-ray imaging using a tilted angle CZT detector was simulated. SNR improvement due to optimal photon energy weighting was 23% and 14% when adipose contrast element, inserted in soft tissue with 10 cm and 20 cm thickness, respectively, was imaged using 5 energy bins and weighting factors optimized for the adipose. SNR improvement was 42% and 31% when CaCO(3) contrast element, inserted in soft tissue with 10 cm and 20 cm thickness, respectively, was imaged using 5 energy bins and weighting factors optimized for CaCO(3). The SNRs of the photon counting single-kVp dual-energy subtracted images of CaCO(3) and adipose were higher by 2.04 and 2.74 times, respectively, as compared to currently used dual-kVp dual-energy subtracted images. Experiments with a CZT crystal with 2 mm thickness have shown significant decrease in the tailing effect of the CZT pulse spectrum at 59 keV and 122 keV photon energies, when the tilting angle configuration was used. Finally, feasibility of the tilted angle CZT detector for photon counting cone beam breast CT imaging was demonstrated.

Scatter rejection in multislit digital mammography

The scatter to primary ratio (SPR) was measured on a scanning multislit full-field digital mammography system for different thickness of breast equivalent material and different tube voltages. Scatter within the detector was measured separately and was found to be the major source of scatter in the assembly. Measured total SPRs below 6% are reported for breast range 3-7 cm. The performance of the multislit assembly is compared to other imaging geometries with different scatter rejection schemes by using the scatter detective quantum efficiency.

Development and Monte Carlo analysis of antiscatter grids for mammography

Mammography arguably demands the highest fidelity of all x-ray imaging applications, with simultaneous requirements of exceedingly high spatial and contrast resolution. Continuing technical improvements of screen-film and digital mammography systems have led to substantial improvements in image quality, and therefore improvements in the performance of anti-scatter grids are required to keep pace with the improvements in other components of the imaging chain. The development of an air-core honeycomb (cellular) grid using x-ray lithography and electroforming techniques is described, and the production of a 60 mm x 60 mm section of grid is reported. A crossed grid was constructed with 25 microm copper septa, and a period of 550 microm. Monte Carlo and numerical simulation methods were used to analyze the theoretical performance of the fabricated grid, and comparisons with other grid systems (Lorad HTC and carbon fiber interspaced grids) were made over a range of grid ratios. The results demonstrate essentially equivalent performance in terms of contrast improvement factor (CIF) and Bucky factor (BF) between Cu and Au honeycomb grids and the Lorad HTC (itself a copper honeycomb grid). Gold septa improved both CIF and BF performance in higher kVp, higher scatter geometries. The selectivity of honeycomb grids was far better than for linear grids, with a factor of approximately 3.9 improvement at a grid ratio of 5.0. It is concluded that using the fabrication methods described, that practical honeycomb grid structures could be produced for use in mammographic imaging, and that a substantial improvement in scatter rejection would be achieved using these devices.