Resolution at oblique incidence angles of a flat panel imager for breast tomosynthesis

Optimizing angular range in digital breast tomosynthesis: A phantom study investigating lesion detection across varied breast density and thickness

PURPOSE

To determine the optimal angular range (AR) for digital breast tomosynthesis (DBT) systems that provides highest lesion visibility across various breast densities and thicknesses.

METHOD

A modular DBT phantom, consisting of tissue-equivalent adipose and glandular modules, along with a module embedded with test objects (speckles, masses, fibers), was used to create combinations simulating different breast thicknesses, densities, and lesion locations. A prototype DBT system operated at four ARs (AR±7.5°, AR±12.5°, AR±19°, and AR±25°) to acquire 11 projection images for each combination, with separate fixed doses for thin and thick combinations. Three blinded radiologists independently assessed lesion visibility in reconstructed images; assessments were averaged and compared using linear mixed models.

RESULTS

Speckle visibility was highest with AR±7.5° or AR±12.5°, decreasing with wider ARs in all density and thickness combinations. The difference between AR±7.5° and AR±12.5° was not statistically significant, except for the tube-side speckles in thin-fatty combinations (5.83 [AR±7.5°] vs. 5.39 [AR±12.5°], P = 0.019). Mass visibility was not affected by AR in thick combinations, while AR±12.5° exhibited the highest mass visibility for both thin-fatty and thin-dense combinations (P = 0.032 and 0.007, respectively). Different ARs provided highest fiber visibility for different combinations; however, AR±12.5° consistently provided highest or comparable visibility. AR±12.5° showed highest overall lesion visibility for all density and thickness combinations.

CONCLUSIONS

AR±12.5° exhibited the highest overall lesion visibility across various phantom thicknesses and densities using a projection number of 11.

Performance evaluation of digital breast tomosynthesis systems: physical methods and experimental data

Digital breast tomosynthesis (DBT) has become a well-established breast imaging technique, whose performance has been investigated in many clinical studies, including a number of prospective clinical trials. Results from these studies generally point to non-inferiority in terms of microcalcification detection and superior mass-lesion detection for DBT imaging compared to digital mammography (DM). This modality has become an essential tool in the clinic for assessment and ad-hoc screening but is not yet implemented in most breast screening programmes at a state or national level. While evidence on the clinical utility of DBT has been accumulating, there has also been progress in the development of methods for technical performance assessment and quality control of these imaging systems. DBT is a relatively complicated ‘pseudo-3D’ modality whose technical assessment poses a number of difficulties. This paper reviews methods for the technical performance assessment of DBT devices, starting at the component level in part one and leading up to discussion of system evaluation with physical test objects in part two. We provide some historical and basic theoretical perspective, often starting from methods developed for DM imaging. Data from a multi-vendor comparison are also included, acquired under the medical physics quality control protocol developed by EUREF and currently being consolidated by a European Federation of Organisations for Medical Physics working group. These data and associated methods can serve as a reference for the development of reference data and provide some context for clinical studies.

Endorectal digital prostate tomosynthesis (endoDPT): a proof-of-concept study

In this study we present endorectal digital prostate tomosynthesis (endoDPT), a proposed method of high resolution prostate imaging using an endorectal x-ray sensor and an external x-ray source. endoDPT may be useful for visualizing the fine detail of small structures such as low dose rate brachytherapy (LDRBT) seeds that are difficult to visualize with current methods. The resolution of endoDPT was characterized through measurement of the modulation transfer function (MTF) and artifact spread function (ASF) in computational and physical phantoms. The qualitative resolution of endoDPT was assessed relative to computed tomography (CT) through imaging of LDRBT seeds implanted inex vivocanine prostates. The x-ray sensor MTF reached 10% at 11.50 mm^-1, the reconstruction algorithm MTF reached a maximum at 7.08 mm^-1, and the ASF was 2.5 mm (full-width at half-maximum). Fine structures in LDRBT seeds like the 0.05 mm thick shell were visible with endoDPT but not CT. All endoDPT images exhibited an overshoot artifact. The measured MTFs were consistent with other studies using similar x-ray sensors and demonstrated improved resolution compared to digital breast tomosynthesis; this result was due to the smaller endoDPT x-ray sensor detection element size and quantitatively demonstrates the high resolution of endoDPT. The ASF demonstrated worse depth resolution compared to in-plane resolution, due to partial angular sampling; partial angular sampling also caused the observed overshoot artifact in the endoDPT images. However, endoDPT still was able to visualize fine structures such as the LDRBT seed shell to a much higher degree than CT. This high-resolution visualization may be useful for improvements in patient specific LDRBT dosimetry. Overall, these results indicate endoDPT is capable of high in-plane spatial resolution and is thus well poised for optimization and studies assessing clinical utility.

MTF of columnar phosphors with a homogenous part: an analytical approach

A method for the theoretical estimation of the MTF of columnar phosphors with a homogeneous part at the end used in X-ray imaging has been developed. This method considers the light transport inside the scintillator through an analytical modelling, the optical photon beams distribution on the scintillator-optical sensor interface, and uses the definition of the PSF and a Gauss fitted LSF to estimate the MTF of an indirect detector. This method was applied to a columnar CsI:Tl scintillator and validated against experimental results found in literature, and a good agreement was observed. It was found that, by increasing the pixel size of the optical detector and the thickness of the scintillator, the MTF decreased as expected. This method may be used in evaluating the performance of the columnar phosphors used in medical imaging, given their physical and geometrical characteristics.Graphical abstract (a) Side view of a part of the scintillator where five crystal columns with homogeneous ends attached to an optical sensor is shown. (b) Propagation of two random optical photon beams emitted from point K with different angles of emission is shown. All the symbols are explained analytically in the text.

The role of digital breast tomosynthesis in breast cancer screening: a manufacturer- and metrics-specific analysis

Aim

Digital Breast Tomosynthesis (DBT), with or without Digital Mammography (DM) or Synthetic Mammography (SM), has been introduced or is under consideration for its introduction in breast cancer screening in several countries, as it has been shown that it has advantages over DM. Despite this there is no agreement on how to implement DBT in screening, and in many cases there is a lack of official guidance on the optimum usage of each commercially available system. The aim of this review is to carry out a manufacturer-specific summary of studies on the implementation of DBT in breast cancer screening.

Methods

An exhaustive literature review was undertaken to identify clinical observer studies that evaluated at least one of five common metrics: sensitivity, specificity, area under the curve (AUC) of the receiver-operating characteristics (ROC) analysis, recall rate and cancer detection rate. Four common DBT implementation methods were discussed in this review: (1) DBT, (2) DM with DBT, (3) 1-view DBT with or without 1-view DM or 2-view DM and (4) DBT with SM.

Results

A summary of 89 studies, selected from a database of 677 studies, on the assessment of the implementation of DBT in breast cancer screening is presented in tables and discussed in a manufacturer- and metric-specific approach. Much more studies were carried out using some DBT systems than others. For one implementation method of DBT by one manufacturer there is a shortage of studies, for another implementation there are conflicting results. In some cases, there is a strong agreement between studies, making the advantages and disadvantages of each system clear.

Conclusion

The optimum implementation method of DBT in breast screening, in terms of diagnostic benefit and patient radiation dose, for one manufacturer does not necessarily apply to other manufacturers.

The threshold detectable mass diameter for 2D-mammography and digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is currently under consideration for replacement of, or combined use with 2D-mammography in national breast screening programmes. To investigate the potential benefits that DBT can bring to screening, the threshold detectable lesion diameters were measured for different forms of DBT in comparison to 2D-mammography. The aim of this study was to compare the threshold detectable mass diameters obtained with narrow angle (15°/15 projections) and wide angle (50°/25 projections) DBT in comparison to 2D-mammography. Simulated images of 60 mm thick compressed breasts were produced with and without masses using a set of validated image modelling tools for 2D-mammography and DBT. Image processing and reconstruction were performed using commercial software. A series of 4-alternative forced choice (4AFC) experiments was conducted for signal detection with the masses as targets. The threshold detectable mass diameter was found for each imaging modality with a mean glandular dose of 2.5 mGy. The resulting values of the threshold diameter for 2D-mammography (10.2 ± 1.4 mm) were found to be larger (p < 0.001) than those for narrow angle DBT (6.0 ± 1.1 mm) and wide angle DBT (5.6 ± 1.2 mm). There was no significant difference between the threshold diameters for wide and narrow angle DBT. Implications for the introduction of DBT alone or in combination with 2D-mammography in breast cancer screening are discussed.

Nonstationary model of oblique x-ray incidence in amorphous selenium detectors: II. Transfer functions

PURPOSE

One limitation of experimental techniques for quantifying resolution and noise in detectors is that the measurement is made in a region-of-interest (ROI). With theoretical modeling, these properties can be measured at a point, allowing for quantification of spatial anisotropy. This paper calculates nonstationary transfer functions for amorphous selenium (a-Se) detectors in breast imaging. We use this model to demonstrate the performance advantage of a "next-generation" tomosynthesis (NGT) system, which is capable of x-ray source motion with more degrees of freedom than a clinical tomosynthesis system.

METHODS

Using Swank's formulation, the optical transfer function (OTF) and presampled noise power spectra (NPS) are determined based on the point spread function derived in Part 1. The modulation transfer function (MTF) is found from the normalized modulus of the OTF. To take into account the presence of digitization, the presampled NPS is convolved with a two-dimensional comb function, for which the period along each direction is the reciprocal of the detector element size. The detective quantum efficiency (DQE) is then determined from combined knowledge of the OTF and NPS.

RESULTS

First, the model is used to demonstrate the loss of image quality due to oblique x-ray incidence. The MTF is calculated along various polar angles, corresponding to different orientations of the input frequency. The MTF is independent of the incidence angle if the polar angle is perpendicular to the ray incidence direction. However, along other polar angles, oblique incidence results in MTF degradation at high frequencies. The MTF degradation is most substantial along the ray incidence direction. Unlike the MTF, the normalized NPS (NNPS) is independent of the incidence angle. To measure the relative signal-to-noise, the DQE is also calculated. Oblique incidence yields high-frequency DQE degradation, which is more pronounced than the MTF degradation. This arises because the DQE is proportionate with the square of the MTF. Ultimately, this model is used to evaluate how the image quality varies over the detector area. For various projection images, we calculate the variation in the incidence angle over this area. With the NGT system, the source can be positioned in such a way that this variation is minimized, and hence the DQE exhibits less anisotropy. To achieve this improvement in the image quality, the source needs to have a component of motion in the posteroanterior (PA) direction, which is perpendicular to the conventional direction of source motion in tomosynthesis.

CONCLUSIONS

In a-Se detectors, the DQE at high frequencies is degraded due to oblique incidence. The DQE degradation is more pronounced than the MTF degradation. This model is used to quantify the spatial variation in DQE over the detector area. The use of PA source motion is a strategy for minimizing this variation and thus improving the image quality.

Cascaded systems analysis of shift-variant image quality in slit-scanning breast tomosynthesis

PURPOSE

Digital breast tomosynthesis (DBT) is becoming an important part of breast cancer screening and diagnosis. Compared to two-dimensional mammography, tomosynthesis introduces limited three-dimensional (3D) resolution, but maintains high in-plane resolution, low dose, and allows for similar clinical protocols. The scanning motion and oblique projections of tomosynthesis acquisitions introduce shift-variance to the image quality, in addition to effects such as source blurring and geometric magnification. Shift-variant detector response caused by oblique incidence has been extensively studied previously and is most easily mitigated by letting the source and detector move in sync. In addition, conical reconstruction grids, that is, a grid aligned with the central tomosynthesis projection, have been proposed to compensate for magnification effects. This paper introduces a shift-variant cascaded systems model for tomosynthesis and validates it against measurements. As an example, the model was used to investigate the shift-variance of a tomosynthesis system.

METHODS

The shift-variant cascaded systems model was validated on a slit-scanning photon-counting DBT system, with synchronous source-detector movement, using simple back-projection in a conical reconstruction volume. The modulation transfer function (MTF), normalized noise-power spectrum (NNPS), and detective quantum efficiency (DQE) were used as figures of merit. Simulations were performed for single points while measurements were done over a finite volume, assuming local shift invariance. To investigate the full extent of shift-variance, 80 locations across the volume were simulated, and the MTF and DQE at 2.5 lp/mm were calculated as a function of position.

RESULTS

The simulated metrics generally agreed well with their corresponding measurements. The frequency at 50% MTF along the scan direction showed a relatively small variation, ranging from 2.1 to 2.4 lp/mm for the different locations. The frequency at 50% MTF along the chest-mammilla direction showed a larger variation, ranging from 2.9 to 3.8 lp/mm. All points exhibited a similarly shaped NNPS but the noise magnitude varied with slice height. The zero-frequency DQE in reconstructed slices was lower than that of the projections, an effect likely caused by noise-aliasing increasing the zero-frequency noise.

CONCLUSIONS

A shift-variant cascaded systems model has been developed for slit-scanning tomosynthesis using simple back-projection. The model was successfully validated against measurements. Even though the study was performed on a slit-scanning system, several parts of the framework can be applied and extended to other tomosynthesis geometries. The conical reconstruction grid has low variation in image quality in the scan direction where the 3D information is acquired, but source and geometric magnification still dominate in the slit direction, causing a larger variation in image quality. We conclude that image quality is close to shift-invariant in the scan direction, but not in the height and chest-mammilla directions, and we recommend that small measurement volumes are used when measuring image quality in these directions to minimize the effects of shift variance.

Analysis of Volume Overestimation Artifacts in the Breast Outline Segmentation in Tomosynthesis

In digital breast tomosynthesis (DBT), the reconstruction is calculated from x-ray projection images acquired over a small range of angles. One step in the reconstruction process is to identify the pixels that fall outside the shadow of the breast, to segment the breast from the background (air). In each projection, rays are back-projected from these pixels to the focal spot. All voxels along these rays are identified as air. By combining these results over all projections, a breast outline can be determined for the reconstruction. This paper quantifies the accuracy of this breast segmentation strategy in DBT. In this study, a physical phantom modeling a breast under compression was analyzed with a prototype next-generation tomosynthesis (NGT) system described in previous work. Multiple wires were wrapped around the phantom. Since the wires are thin and high contrast, their exact location can be determined from the reconstruction. Breast parenchyma was portrayed outside the outline defined by the wires. Specifically, the size of the phantom was overestimated along the posteroanterior (PA) direction; i.e., perpendicular to the plane of conventional source motion. To analyze how the acquisition geometry affects the accuracy of the breast outline segmentation, a computational phantom was also simulated. The simulation identified two ways to improve the segmentation accuracy; either by increasing the angular range of source motion laterally or by increasing the range in the PA direction. The latter approach is a unique feature of the NGT design; the advantage of this approach was validated with our prototype system.

Modelling the transport of optical photons in scintillation detectors for diagnostic and radiotherapy imaging

Computational modelling of radiation transport can enhance the understanding of the relative importance of individual processes involved in imaging systems. Modelling is a powerful tool for improving detector designs in ways that are impractical or impossible to achieve through experimental measurements. Modelling of light transport in scintillation detectors used in radiology and radiotherapy imaging that rely on the detection of visible light plays an increasingly important role in detector design. Historically, researchers have invested heavily in modelling the transport of ionizing radiation while light transport is often ignored or coarsely modelled. Due to the complexity of existing light transport simulation tools and the breadth of custom codes developed by users, light transport studies are seldom fully exploited and have not reached their full potential. This topical review aims at providing an overview of the methods employed in freely available and other described optical Monte Carlo packages and analytical models and discussing their respective advantages and limitations. In particular, applications of optical transport modelling in nuclear medicine, diagnostic and radiotherapy imaging are described. A discussion on the evolution of these modelling tools into future developments and applications is presented.

Task-Based Modeling of a 5k Ultra-High-Resolution Medical Imaging System for Digital Breast Tomosynthesis

High-resolution, low-noise X-ray detectors based on CMOS active pixel sensor (APS) technology have demonstrated superior imaging performance for digital breast tomosynthesis (DBT). This paper presents a task-based model for a high-resolution medical imaging system to evaluate its ability to detect simulated microcalcifications and masses as lesions for breast cancer. A 3-D cascaded system analysis for a 50- [Formula: see text] pixel pitch CMOS APS X-ray detector was integrated with an object task function, a medical imaging display model, and the human eye contrast sensitivity function to calculate the detectability index and area under the ROC curve (AUC). It was demonstrated that the display pixel pitch and zoom factor should be optimized to improve the AUC for detecting small microcalcifications. In addition, detector electronic noise of smaller than 300 e^- and a high display maximum luminance (>1000 cd/cm ²) are desirable to distinguish microcalcifications of [Formula: see text] in size. For low contrast mass detection, a medical imaging display with a minimum of 12-bit gray levels is recommended to realize accurate luminance levels. A wide projection angle range of greater than ±30° in combination with the image gray level magnification could improve the mass detectability especially when the anatomical background noise is high. On the other hand, a narrower projection angle range below ±20° can improve the small, high contrast object detection. Due to the low mass contrast and luminance, the ambient luminance should be controlled below 5 cd/ [Formula: see text]. Task-based modeling provides important firsthand imaging performance of the high-resolution CMOS-based medical imaging system that is still at early stage development for DBT. The modeling results could guide the prototype design and clinical studies in the future.

Characterisation of noise and sharpness of images from four digital breast tomosynthesis systems for simulation of images for virtual clinical trials

In-depth evaluation of the noise and sharpness characteristics of FujiFilm Innovality, GE SenoClaire, Hologic Selenia Dimensions and Siemens Inspiration digital breast tomosynthesis (DBT) systems was performed with the intention of improving image simulation for virtual clinical trials. Noise power spectra (NPS) and modulation transfer function curves (MTF) were measured for planar modes and for the first and central projections for DBT modes. In DBT mode, the x-ray beam was blocked for the projections before the central projection in order to remove the influence of lag and ghosting from the previous images. A quadratic fit between the NPS and linearised pixel value gave the noise coefficients for planar and DBT imaging modes. The spatial frequencies corresponding to an MTF of 0.5 (MTF0.5) were calculated from the MTF measurements made on the breast support and at 40 mm above the breast support. This was done for the first and the central projections. The percentage of signal carried over from the first projection to subsequent images (lag) was measured using a slit. The noise associated with lag was also evaluated. The DBT modes typically had lower electronic noise coefficients but higher structural noise coefficients compared to the respective planar mode MTF0.5 measured 40 mm above the table was between 6% and 47% lower for continuous scanning systems compared to 1% lower for step and shoot systems. For wide angle DBT, the MTF0.5 of the first projection was 18% (FujiFilm) and 28% (Siemens) lower than for the central projection. Lag in the second projection was 2.2%, 0.3%, 0.8% for the FujiFilm, GE and Hologic systems respectively. In all cases, the noise associated with lag was negligible. Current modelling frameworks for virtual clinical trials of breast DBT systems need to be adapted to account for signals from lag and variations in the MTF at wide angles.

The effect of system geometry and dose on the threshold detectable calcification diameter in 2D-mammography and digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is under consideration to replace or to be used in combination with 2D-mammography in breast screening. The aim of this study was the comparison of the detection of microcalcification clusters by human observers in simulated breast images using 2D-mammography, narrow angle (15°/15 projections) and wide angle (50°/25 projections) DBT. The effects of the cluster height in the breast and the dose to the breast on calcification detection were also tested. Simulated images of 6 cm thick compressed breasts were produced with and without microcalcification clusters inserted, using a set of image modelling tools for 2D-mammography and DBT. Image processing and reconstruction were performed using commercial software. A series of 4-alternative forced choice (4AFC) experiments was conducted for signal detection with the microcalcification clusters as targets. Threshold detectable calcification diameter was found for each imaging modality with standard dose: 2D-mammography: 2D-mammography (165  ±  9 µm), narrow angle DBT (211  ±  11 µm) and wide angle DBT (257  ±  14 µm). Statistically significant differences were found when using different doses, but different geometries had a greater effect. No differences were found between the threshold detectable calcification diameters at different heights in the breast. Calcification clusters may have a lower detectability using DBT than 2D imaging.

Digital Breast Tomosynthesis: State of the Art

This topical review on digital breast tomosynthesis (DBT) is provided with the intent of describing the state of the art in terms of technology, results from recent clinical studies, advanced applications, and ongoing efforts to develop multimodality imaging systems that include DBT. Particular emphasis is placed on clinical studies. The observations of increase in cancer detection rates, particularly for invasive cancers, and the reduction in false-positive rates with DBT in prospective trials indicate its benefit for breast cancer screening. Retrospective multireader multicase studies show either noninferiority or superiority of DBT compared with mammography. Methods to curtail radiation dose are of importance. (©) RSNA, 2015.

The effect of amorphous selenium detector thickness on dual-energy digital breast imaging

PURPOSE

Contrast enhanced (CE) imaging techniques for both planar digital mammography (DM) and three-dimensional (3D) digital breast tomosynthesis (DBT) applications requires x-ray photon energies higher than the k-edge of iodine (33.2 keV). As a result, x-ray tube potentials much higher (>40 kVp) than those typical for screening mammography must be utilized. Amorphous selenium (a-Se) based direct conversion flat-panel imagers (FPI) have been widely used in DM and DBT imaging systems. The a-Se layer is typically 200 μm thick with quantum detective efficiency (QDE) >87% for x-ray energies below 26 keV. However, QDE decreases substantially above this energy. To improve the object detectability of either CE-DM or CE-DBT, it may be advantageous to increase the thickness (dSe) of the a-Se layer. Increasing the dSe will improve the detective quantum efficiency (DQE) at the higher energies used in CE imaging. However, because most DBT systems are designed with partially isocentric geometries, where the gantry moves about a stationary detector, the oblique entry of x-rays will introduce additional blur to the system. The present investigation quantifies the effect of a-Se thickness on imaging performance for both CE-DM and CE-DBT, discussing the effects of improving photon absorption and blurring from oblique entry of x-rays.

METHODS

In this paper, a cascaded linear system model (CLSM) was used to investigate the effect of dSe on the imaging performance (i.e., MTF, NPS, and DQE) of FPI in CE-DM and CE-DBT. The results from the model are used to calculate the ideal observer signal-to-noise ratio, d', which is used as a figure-of-merit to determine the total effect of increasing dSe for CE-DM and CE-DBT.

RESULTS

The results of the CLSM show that increasing dSe causes a substantial increase in QDE at the high energies used in CE-DM. However, at the oblique projection angles used in DBT, the increased length of penetration through a-Se introduces additional image blur. The reduced MTF and DQE at high spatial frequencies lead to reduced two-dimensional d'. These losses in projection image resolution may subsequently result in a decrease in the 3D d', but the degree of which is largely dependent on the DBT reconstruction algorithm. For a filtered backprojection (FBP) algorithm with spectral apodization and slice-thickness filters, which dominate the blur for reconstructed images at oblique angles, the effect of oblique entry of x-rays on 3D d' is minimal. Thus, increasing dSe results in an improvement in d' for both CE-DM and CE-DBT with typical FBP reconstruction parameters.

CONCLUSIONS

Increased dSe improves CE breast imaging performance by increasing QDE of detectors at higher energies, e.g., 49 kVp. Although there is additional blur in the oblique angled projections of a DBT scan, the overall 3D d' for DBT is not degraded because the dominant source blur at these angles results from the reconstruction filters of the employed FBP algorithm.

Characterization of a high-energy in-line phase contrast tomosynthesis prototype

PURPOSE

In this research, a high-energy in-line phase contrast tomosynthesis prototype was developed and characterized through quantitative investigations and phantom studies.

METHODS

The prototype system consists of an x-ray source, a motorized rotation stage, and a CMOS detector with a pixel pitch of 0.05 mm. The x-ray source was operated at 120 kVp for this study, and the objects were mounted on the rotation stage 76.2 cm (R1) from the source and 114.3 cm (R2) from the detector. The large air gap between the object and detector guarantees sufficient phase-shift effects. The quantitative evaluation of this prototype included modulation transfer function and noise power spectrum measurements conducted under both projection mode and tomosynthesis mode. Phantom studies were performed including three custom designed phantoms with complex structures: a five-layer bubble wrap phantom, a fishbone phantom, and a chicken breast phantom with embedded fibrils and mass structures extracted from an ACR phantom. In-plane images of the phantoms were acquired to investigate their image qualities through observation, intensity profile plots, edge enhancement evaluations, and/or contrast-to-noise ratio calculations. In addition, the robust phase-attenuation duality (PAD)-based phase retrieval method was applied to tomosynthesis for the first time in this research. It was utilized as a preprocessing method to fully exhibit phase contrast on the angular projection before reconstruction.

RESULTS

The resolution and noise characteristics of this high-energy in-line phase contrast tomosynthesis prototype were successfully investigated and demonstrated. The phantom studies demonstrated that this imaging prototype can successfully remove the structure overlapping in phantom projections, obtain delineate interfaces, and achieve better contrast-to-noise ratio after applying phase retrieval to the angular projections.

CONCLUSIONS

This research successfully demonstrated a high-energy in-line phase contrast tomosynthesis prototype. In addition, the PAD-based method of phase retrieval was combined with tomosynthesis imaging for the first time, which demonstrated its capability in significantly improving the contrast-to-noise ratios in the images.

Digital breast tomosynthesis: studies of the effects of acquisition geometry on contrast-to-noise ratio and observer preference of low-contrast objects in breast phantom images

The effect of acquisition geometry in digital breast tomosynthesis was evaluated with studies of contrast-to-noise ratios (CNRs) and observer preference. Contrast-detail (CD) test objects in 5 cm thick phantoms with breast-like backgrounds were imaged. Twelve different angular acquisitions (average glandular dose for each ~1.1 mGy) were performed ranging from narrow angle 16° with 17 projection views (16d17p) to wide angle 64d17p. Focal slices of SART-reconstructed images of the CD arrays were selected for CNR computations and the reader preference study. For the latter, pairs of images obtained with different acquisition geometries were randomized and scored by 7 trained readers. The total scores for all images and readings for each acquisition geometry were compared as were the CNRs. In general, readers preferred images acquired with wide angle as opposed to narrow angle geometries. The mean percent preferred was highly correlated with tomosynthesis angle (R = 0.91). The highest scoring geometries were 60d21p (95%), 64d17p (80%), and 48d17p (72%); the lowest scoring were 16d17p (4%), 24d9p (17%) and 24d13p (33%). The measured CNRs for the various acquisitions showed much overlap but were overall highest for wide-angle acquisitions. Finally, the mean reader scores were well correlated with the mean CNRs (R = 0.83).

Development and validation of a modelling framework for simulating 2D-mammography and breast tomosynthesis images

Planar 2D x-ray mammography is generally accepted as the preferred screening technique used for breast cancer detection. Recently, digital breast tomosynthesis (DBT) has been introduced to overcome some of the inherent limitations of conventional planar imaging, and future technological enhancements are expected to result in the introduction of further innovative modalities. However, it is crucial to understand the impact of any new imaging technology or methodology on cancer detection rates and patient recall. Any such assessment conventionally requires large scale clinical trials demanding significant investment in time and resources. The concept of virtual clinical trials and virtual performance assessment may offer a viable alternative to this approach. However, virtual approaches require a collection of specialized modelling tools which can be used to emulate the image acquisition process and simulate images of a quality indistinguishable from their real clinical counterparts. In this paper, we present two image simulation chains constructed using modelling tools that can be used for the evaluation of 2D-mammography and DBT systems. We validate both approaches by comparing simulated images with real images acquired using the system being simulated. A comparison of the contrast-to-noise ratios and image blurring for real and simulated images of test objects shows good agreement ( < 9% error). This suggests that our simulation approach is a promising alternative to conventional physical performance assessment followed by large scale clinical trials.

Digital breast tomosynthesis: observer performance of clustered microcalcification detection on breast phantom images acquired with an experimental system using variable scan angles, angular increments, and number of projection views

PURPOSE

To investigate the dependence of microcalcification cluster detectability on tomographic scan angle, angular increment, and number of projection views acquired at digital breast tomosynthesis ( DBT digital breast tomosynthesis ).

MATERIALS AND METHODS

A prototype DBT digital breast tomosynthesis system operated in step-and-shoot mode was used to image breast phantoms. Four 5-cm-thick phantoms embedded with 81 simulated microcalcification clusters of three speck sizes (subtle, medium, and obvious) were imaged by using a rhodium target and rhodium filter with 29 kV, 50 mAs, and seven acquisition protocols. Fixed angular increments were used in four protocols (denoted as scan angle, angular increment, and number of projection views, respectively: 16°, 1°, and 17; 24°, 3°, and nine; 30°, 3°, and 11; and 60°, 3°, and 21), and variable increments were used in three (40°, variable, and 13; 40°, variable, and 15; and 60°, variable, and 21). The reconstructed DBT digital breast tomosynthesis images were interpreted by six radiologists who located the microcalcification clusters and rated their conspicuity.

RESULTS

The mean sensitivity for detection of subtle clusters ranged from 80% (22.5 of 28) to 96% (26.8 of 28) for the seven DBT digital breast tomosynthesis protocols; the highest sensitivity was achieved with the 16°, 1°, and 17 protocol (96%), but the difference was significant only for the 60°, 3°, and 21 protocol (80%, P < .002) and did not reach significance for the other five protocols (P = .01-.15). The mean sensitivity for detection of medium and obvious clusters ranged from 97% (28.2 of 29) to 100% (24 of 24), but the differences fell short of significance (P = .08 to >.99). The conspicuity of subtle and medium clusters with the 16°, 1°, and 17 protocol was rated higher than those with other protocols; the differences were significant for subtle clusters with the 24°, 3°, and nine protocol and for medium clusters with 24°, 3°, and nine; 30°, 3°, and 11; 60°, 3° and 21; and 60°, variable, and 21 protocols (P < .002).

CONCLUSION

With imaging that did not include x-ray source motion or patient motion during acquisition of the projection views, narrow-angle DBT digital breast tomosynthesis provided higher sensitivity and conspicuity than wide-angle DBT digital breast tomosynthesis for subtle microcalcification clusters.

Carbon nanotube electron field emitters for x-ray imaging of human breast cancer

For imaging human breast cancer, digital breast tomosynthesis (DBT) has been shown to improve image quality and breast cancer detection in comparison to two-dimensional (2D) mammography. Current DBT systems have limited spatial resolution and lengthy scan times. Stationary DBT (s-DBT), utilizing an array of carbon nanotube (CNT) field emission x-ray sources, provides increased spatial resolution and potentially faster imaging than current DBT systems. This study presents the results of detailed evaluations of CNT cathodes for x-ray breast imaging tasks. The following were investigated: high current, long-term stability of CNT cathodes for DBT; feasibility of using CNT cathodes to perform a 2D radiograph function; and cathode performance through several years of imaging. Results show that a breast tomosynthesis system using CNT cathodes could run far beyond the experimentally tested lifetime of one to two years. CNT cathodes were found capable of producing higher currents than typical DBT would require, indicating that the s-DBT imaging time can be further reduced. The feasibility of using a single cathode of the s-DBT tube to perform 2D mammography in 4 s was demonstrated. Over the lifetime of the prototype s-DBT system, it was found that both cathode performance and transmission rate were stable and consistent.

Optimization of configuration parameters in a newly developed digital breast tomosynthesis system

The purpose of the present work was to investigate the effects of variable projection-view (PV) and angular dose (AD) distributions on the reconstructed image quality for improving microcalcification detection. The PV densities at central and peripheral sites were varied through the distribution of 21 PVs acquired over ± 25° angular range. To vary the AD distribution, 7 PVs in the central region were targeted with two, four and six times the peripheral dose, and the number of central PVs receiving four times the peripheral dose was increased from 3 to 11. The contrast-to-noise ratio (CNR) for in-focus plane quality and the full width at half maximum (FWHM) of artifact spread function (ASF) for resolution in the z-direction were used. Although the ASF improved with increasing PV densities at two peripheral sites, the CNRs were inferior to those obtained with other subsets. With increasing PV density in the central area, the vertical resolution decreased but the CNR increased. Although increasing the central PV or AD concentrations improved image quality, excessive central densities reduced image quality by increasing noise in peripheral views.

Oblique reconstructions in tomosynthesis. II. Super-resolution

PURPOSE

In tomosynthesis, super-resolution has been demonstrated using reconstruction planes parallel to the detector. Super-resolution allows for subpixel resolution relative to the detector. The purpose of this work is to develop an analytical model that generalizes super-resolution to oblique reconstruction planes.

METHODS

In a digital tomosynthesis system, a sinusoidal test object is modeled along oblique angles (i.e., "pitches") relative to the plane of the detector in a 3D divergent-beam acquisition geometry. To investigate the potential for super-resolution, the input frequency is specified to be greater than the alias frequency of the detector. Reconstructions are evaluated in an oblique plane along the extent of the object using simple backprojection (SBP) and filtered backprojection (FBP). By comparing the amplitude of the reconstruction against the attenuation coefficient of the object at various frequencies, the modulation transfer function (MTF) is calculated to determine whether modulation is within detectable limits for super-resolution. For experimental validation of super-resolution, a goniometry stand was used to orient a bar pattern phantom along various pitches relative to the breast support in a commercial digital breast tomosynthesis system.

RESULTS

Using theoretical modeling, it is shown that a single projection image cannot resolve a sine input whose frequency exceeds the detector alias frequency. The high frequency input is correctly visualized in SBP or FBP reconstruction using a slice along the pitch of the object. The Fourier transform of this reconstructed slice is maximized at the input frequency as proof that the object is resolved. Consistent with the theoretical results, experimental images of a bar pattern phantom showed super-resolution in oblique reconstructions. At various pitches, the highest frequency with detectable modulation was determined by visual inspection of the bar patterns. The dependency of the highest detectable frequency on pitch followed the same trend as the analytical model. It was demonstrated that super-resolution is not achievable if the pitch of the object approaches 90°, corresponding to the case in which the test frequency is perpendicular to the breast support. Only low frequency objects are detectable at pitches close to 90°.

CONCLUSIONS

This work provides a platform for investigating super-resolution in oblique reconstructions for tomosynthesis. In breast imaging, this study should have applications in visualizing microcalcifications and other subtle signs of cancer.

A virtual trial framework for quantifying the detectability of masses in breast tomosynthesis projection data

PURPOSE

Digital breast tomosynthesis (DBT) is a promising breast cancer screening tool that has already begun making inroads into clinical practice. However, there is ongoing debate over how to quantitatively evaluate and optimize these systems, because different definitions of image quality can lead to different optimal design strategies. Powerful and accurate tools are desired to extend our understanding of DBT system optimization and validate published design principles.

METHODS

The authors developed a virtual trial framework for task-specific DBT assessment that uses digital phantoms, open-source x-ray transport codes, and a projection-space, spatial-domain observer model for quantitative system evaluation. The authors considered evaluation of reconstruction algorithms as a separate problem and focused on the information content in the raw, unfiltered projection images. Specifically, the authors investigated the effects of scan angle and number of angular projections on detectability of a small (3 mm diameter) signal embedded in randomly-varying anatomical backgrounds. Detectability was measured by the area under the receiver-operating characteristic curve (AUC). Experiments were repeated for three test cases where the detectability-limiting factor was anatomical variability, quantum noise, or electronic noise. The authors also juxtaposed the virtual trial framework with other published studies to illustrate its advantages and disadvantages.

RESULTS

The large number of variables in a virtual DBT study make it difficult to directly compare different authors' results, so each result must be interpreted within the context of the specific virtual trial framework. The following results apply to 25% density phantoms with 5.15 cm compressed thickness and 500 μm(3) voxels (larger 500 μm(2) detector pixels were used to avoid voxel-edge artifacts): 1. For raw, unfiltered projection images in the anatomical-variability-limited regime, AUC appeared to remain constant or increase slightly with scan angle. 2. In the same regime, when the authors fixed the scan angle, AUC increased asymptotically with the number of projections. The threshold number of projections for asymptotic AUC performance depended on the scan angle. In the quantum- and electronic-noise dominant regimes, AUC behaviors as a function of scan angle and number of projections sometimes differed from the anatomy-limited regime. For example, with a fixed scan angle, AUC generally decreased with the number of projections in the electronic-noise dominant regime. These results are intended to demonstrate the capabilities of the virtual trial framework, not to be used as optimization rules for DBT.

CONCLUSIONS

The authors have demonstrated a novel simulation framework and tools for evaluating DBT systems in an objective, task-specific manner. This framework facilitates further investigation of image quality tradeoffs in DBT.

Effective DQE (eDQE) for monoscopic and stereoscopic chest radiography imaging systems with the incorporation of anatomical noise

PURPOSE

Stereoscopic chest biplane correlation imaging (stereo∕BCI) has been proposed as an alternative modality to single view chest x-ray (CXR). The metrics effective modulation transfer function (eMTF), effective normalized noise power spectrum (eNNPS), and effective detective quantum efficiency (eDQE) have been proposed as clinically relevant metrics for assessing clinical system performance taking into consideration the magnification and scatter effects. This study compared the metrics eMTF, eNNPS, eDQE, and detectability index for stereo∕BCI and single view CXR under isodose conditions at two magnifications for two anthropomorphic phantoms of differing sizes.

METHODS

Measurements for the eMTF were taken for two phantom sizes with an opaque edge test device using established techniques. The eNNPS was measured at two isodose conditions for two phantoms using established techniques. The scatter was measured for two phantoms using an established beam stop method. All measurements were also taken at two different magnifications with two phantoms. A geometrical phantom was used for comparison with prior results for CXR although the results for an anatomy free phantom are not expected to vary for BCI.

RESULTS

Stereo∕BCI resulted in improved metrics compared to single view CXR. Results indicated that magnification can potentially improve the detection performance primarily due to the air gap which reduced scatter by ∼20%. For both phantoms, at isodose, eDQE(0) for stereo∕BCI was ∼100 times higher than that for CXR. Magnification at isodose improved eDQE(0) by ∼10 times for stereo∕BCI. Increasing the dose did not improve eDQE. The detectability index for stereo∕BCI was ∼100 times better than single view CXR for all conditions. The detectability index was also not improved with increased dose.

CONCLUSIONS

The findings indicate that stereo∕BCI with magnification may improve detectability of subtle lung nodules compared to single view CXR. Results were improved with magnification for the smaller phantom but not for the larger phantom. The effective DQE and the detectability index did not improve with increasing dose.

Observation of super-resolution in digital breast tomosynthesis

PURPOSE

Digital breast tomosynthesis (DBT) is a 3D x-ray imaging modality in which tomographic sections of the breast are generated from a limited range of tube angles. Because oblique x-ray incidence shifts the image of an object in subpixel detector element increments with each increasing projection angle, it is demonstrated that DBT is capable of super-resolution (i.e., subpixel resolution).

METHODS

By convention, DBT reconstructions are performed on planes parallel to the breast support at various depths of the breast volume. In order for resolution in each reconstructed slice to be comparable to the detector, the pixel size should match that of the detector elements; hence, the highest frequency that can be resolved in the plane of reconstruction is the alias frequency of the detector. This study considers reconstruction grids with much smaller pixelation to visualize higher frequencies. For analytical proof of super-resolution, a theoretical framework is developed in which the reconstruction of a high frequency sinusoidal input is calculated using both simple backprojection (SBP) and filtered backprojection. To study the frequency spectrum of the reconstruction, its Fourier transform is also determined. The experimental feasibility of super-resolution was investigated by acquiring images of a bar pattern phantom with frequencies higher than the detector alias frequency.

RESULTS

Using analytical modeling, it is shown that the central projection cannot resolve frequencies exceeding the detector alias frequency. The Fourier transform of the central projection is maximized at a lower frequency than the input as evidence of aliasing. By contrast, SBP reconstruction can resolve the input, and its Fourier transform is correctly maximized at the input frequency. Incorporating filters into the reconstruction smoothens pixelation artifacts in the spatial domain and reduces spectral leakage in the Fourier domain. It is also demonstrated that the existence of super-resolution is dependent on position in the reconstruction and on the directionality of the input frequency. Consistent with the analytical results, experimental reconstructions of bar patterns showed visibility of frequencies greater than the detector alias frequency. Super-resolution was present at positions predicted from analytical modeling.

CONCLUSIONS

This work demonstrates the existence of super-resolution in DBT. Super-resolution has the potential to impact the visualization of fine structural details in the breast, such as microcalcifications and other subtle signs of cancer.

A review of breast tomosynthesis. Part I. The image acquisition process

Mammography is a very well-established imaging modality for the early detection and diagnosis of breast cancer. However, since the introduction of digital imaging to the realm of radiology, more advanced, and especially tomographic imaging methods have been made possible. One of these methods, breast tomosynthesis, has finally been introduced to the clinic for routine everyday use, with potential to in the future replace mammography for screening for breast cancer. In this two part paper, the extensive research performed during the development of breast tomosynthesis is reviewed, with a focus on the research addressing the medical physics aspects of this imaging modality. This first paper will review the research performed on the issues relevant to the image acquisition process, including system design, optimization of geometry and technique, x-ray scatter, and radiation dose. The companion to this paper will review all other aspects of breast tomosynthesis imaging, including the reconstruction process.

Implementation and evaluation of an expectation maximization reconstruction algorithm for gamma emission breast tomosynthesis

PURPOSE

We are developing a dual modality tomosynthesis breast scanner in which x-ray transmission tomosynthesis and gamma emission tomosynthesis are performed sequentially with the breast in a common configuration. In both modalities projection data are obtained over an angular range of less than 180° from one side of the mildly compressed breast resulting in incomplete and asymmetrical sampling. The objective of this work is to implement and evaluate a maximum likelihood expectation maximization (MLEM) reconstruction algorithm for gamma emission breast tomosynthesis (GEBT).

METHODS

A combination of Monte Carlo simulations and phantom experiments was used to test the MLEM algorithm for GEBT. The algorithm utilizes prior information obtained from the x-ray breast tomosynthesis scan to partially compensate for the incomplete angular sampling and to perform attenuation correction (AC) and resolution recovery (RR). System spatial resolution, image artifacts, lesion contrast, and signal to noise ratio (SNR) were measured as image quality figures of merit. To test the robustness of the reconstruction algorithm and to assess the relative impacts of correction techniques with changing angular range, simulations and experiments were both performed using acquisition angular ranges of 45°, 90° and 135°. For comparison, a single projection containing the same total number of counts as the full GEBT scan was also obtained to simulate planar breast scintigraphy.

RESULTS

The in-plane spatial resolution of the reconstructed GEBT images is independent of source position within the reconstructed volume and independent of acquisition angular range. For 45° acquisitions, spatial resolution in the depth dimension (the direction of breast compression) is degraded with increasing source depth (increasing distance from the collimator surface). Increasing the acquisition angular range from 45° to 135° both greatly reduces this depth dependence and improves the average depth dimension resolution from 10.8 to 4.8 mm. The 135° acquisition results in a near-isotropic, spatially uniform 3D resolution of approximately 4.3 mm full width at half maximum. Background nonuniformity (cupping) artifacts arise primarily from angular incompleteness for small angular range acquisition but primarily from gamma ray attenuation at larger angular range. However, background artifacts can be largely eliminated if both prior information regularization and AC are applied. An artificial decrease in lesion voxel value with increasing lesion depth can also be substantially reduced through a combination of AC and RR. In experiments using compressible gelatin breast phantoms, lesion contrast and SNR are about 2.6-8.8 times and 2.3-5.6 times higher, respectively, in GEBT than in planar breast scintigraphy depending on the acquisition angle, the gamma camera trajectory, and the lesion location. In addition, the strong reduction in lesion contrast and SNR with increasing lesion depth that is observed in planar breast scintigraphy can be largely overcome in GEBT.

CONCLUSIONS

The authors have demonstrated a promising EM-based reconstruction scheme for use in GEBT. Compared to planar breast scintigraphy GEBT provides superior and less position-dependent lesion contrast, lesion SNR, and spatial resolution as well as more accurate quantification of lesion-to-background activity concentration ratio.

Quantification of scattered radiation in projection mammography: four practical methods compared

PURPOSE

Four different practical methodologies of quantifying scattered radiation for two different digital mammographic systems are compared. The study considered both grid in and grid out geometries for two different antiscatter grid types, a typical linear grid and a cellular grid design. The aim was to find quick and reproducible methods that could be used in place of the beam stop technique.

METHODS

The scatter to primary ratio (SPR) and the scatter fraction (SF) were used to quantify scattered radiation as a function of poly(methyl methacrylate) (PMMA) thickness, grid position, and beam quality. The four scatter estimation methods applied were (1) the beam stop method, (2) a hybrid method that combined measured detector (scatter-free) modulation transfer function (MTF) data and a Monte Carlo simulation of the scatter point spread function, (3) from the low frequency drop data taken from the system MTF, and (4) from the edge spread function (ESF) measured in the presence of PMMA. Repeatability error was assessed for all methods.

RESULTS

SPR results acquired with the beam stop method ranged from 0.052 to 0.187 for the system with linear grid and from 0.012 to 0.064 for the cellular grid system, as PMMA thickness was increased from 20 to 80 mm. With the grid removed, beam stop SPR was similar for both systems, ranging between 0.268 and 1.124, for corresponding MTF thicknesses. The direct MTF method had a maximum difference of 24% from the beam stop SPR and SF data for all conditions except the cellular grid in geometry, where maximum difference in SPR was 0.044 (164%). The ESF technique gave large differences from the beam stops for both grid geometries but agreement was within 21% for the grid out geometry. Repeatability error with beam stops was between 1% and 5% for the grid out geometries, while for the grid in cases it was 13% and 87% for the linear and cellular grids, respectively. Repeatability error for the direct MTF method applied to both systems and grid geometries ranged between 3% and 12%.

CONCLUSIONS

All three alternative methods to the beam stop technique gave reasonable estimates of SPR without grid, with a maximum difference of 24% (mean difference 8%). For the grid in geometry, the direct MTF method gave a maximum difference of 24% for the linear grid system, while maximum percentage difference was 119% (absolute difference of 0.042) for the system with the cellular grid, where SPR values were low. Except for cases where the SPR is very low, the direct MTF method offers a quick and reproducible alternative to the beam stop technique.

Optimization of phosphor-based detector design for oblique x-ray incidence in digital breast tomosynthesis

PURPOSE

In digital breast tomosynthesis (DBT), a volumetric reconstruction of the breast is generated from a limited range of x-ray projections. One trade-off of DBT is resolution loss in the projections due to non-normal (i.e., oblique) x-ray incidence. Although degradation in image quality due to oblique incidence has been studied using empirical data and Monte Carlo simulations, a theoretical treatment has been lacking. The purpose of this work is to extend Swank's calculations of the transfer functions of turbid granular phosphors to oblique incidence. The model is ultimately used as a tool for optimizing the design of DBT detectors.

METHODS

A quantum-limited system and 20 keV x-rays are considered. Under these assumptions, the modulation transfer function (MTF) and noise power spectra (NPS) are derived using the diffusion approximation to the Boltzmann equation to model optical scatter within the phosphor. This approach is applicable to a nonstructured scintillator such as gadolinium oxysulfide doped with terbium (Gd(2)O(2)S:Tb), which is commonly used in breast imaging and which can reasonably approximate other detector materials. The detective quantum efficiency (DQE) is then determined from the Nishikawa formulation, where it is written as the product of the x-ray quantum detection efficiency, the Swank factor, and the Lubberts fraction. Transfer functions are calculated for both front- and back-screen configurations, which differ by positioning the photocathode at the exit or entrance point of the x-ray beam, respectively.

RESULTS

In the front-screen configuration, MTF and DQE are found to have considerable angular dependence, while NPS is shown to vary minimally with projection angle. As expected, the high frequency MTF and DQE are degraded substantially at large angles. By contrast, all transfer functions for the back-screen configuration have the advantage of significantly less angular dependence. Using these models, we investigated the possibility for optimizing the design of DBT detectors. As an example optimization strategy, the phosphor thickness which maximizes the DQE at a fixed frequency is analyzed. This work demonstrates that the optimal phosphor thickness for the front-screen is angularly dependent, shifting to lower thickness at higher angles. Conversely, the back-screen is not optimized by a single thickness but instead attains reasonably high DQE values over a large range of thicknesses. Although the back-screen configuration is not suited for current detectors using a glass substrate, it may prove to be preferred in future detectors using newly proposed plastic thin-film transistor (TFT) substrates.

CONCLUSIONS

Using the diffusion approximation to the Boltzmann equation to model the spread of light in a scintillator, this paper develops an analytical model of MTF, NPS, and DQE for a phosphor irradiated obliquely. The model is set apart from other studies on oblique incidence in being derived from first principles. This work has applications in the optimization of DBT detector design.

Analysis of parenchymal texture with digital breast tomosynthesis: comparison with digital mammography and implications for cancer risk assessment

PURPOSE

To correlate the parenchymal texture features at digital breast tomosynthesis (DBT) and digital mammography with breast percent density (PD), an established breast cancer risk factor, in a screening population of women.

MATERIALS AND METHODS

This HIPAA-compliant study was approved by the institutional review board. Bilateral DBT images and digital mammograms from 71 women (mean age, 54 years; age range, 34-75 years) with negative or benign findings at screening mammography were retrospectively collected from a separate institutional review board-approved DBT screening trial (performed from July 2007 to March 2008) in which all women had given written informed consent. Parenchymal texture features of skewness, coarseness, contrast, energy, homogeneity, and fractal dimension were computed from the retroareolar region. Principal component analysis (PCA) was applied to obtain orthogonal texture components. Mammographic PD was estimated with software. Correlation analysis and multiple linear regression with generalized estimating equations were performed to determine the association between texture features and breast PD. Regression was adjusted for age to determine the independent association of texture to breast PD when age was also considered as a predictor variable.

RESULTS

Texture feature correlations to breast PD were stronger with DBT than with digital mammography. Statistically significant correlations (P < .001) were observed for contrast (r = 0.48), energy (r = -0.47), and homogeneity (r = -0.56) at DBT and for contrast (r = 0.26), energy (r = -0.26), and homogeneity (r = -0.33) at digital mammography. Multiple linear regression analysis of PCA texture components as predictors of PD also demonstrated significantly stronger associations with DBT. The association was strongest when age was also considered as a predictor of PD (R² = 0.41 for DBT and 0.28 for digital mammography; P < .001).

CONCLUSION

Parenchymal texture features are more strongly correlated to breast PD in DBT than in digital mammography. The authors' long-term hypothesis is that parenchymal texture analysis with DBT will result in quantitative imaging biomarkers that can improve the estimation of breast cancer risk.

Oblique incidence effects in direct x-ray detectors: a first-order approximation using a physics-based analytical model

PURPOSE

The authors describe the modifications to a previously developed analytical model of indirect CsI:Tl-based detector response required for studying oblique x-ray incidence effects in direct semiconductor-based detectors. This first-order approximation analysis allows the authors to describe the associated degradation in resolution in direct detectors and compare the predictions to the published data for indirect detectors.

METHODS

The proposed model is based on a physics-based analytical description developed by Freed et al. ["A fast, angle-dependent, analytical model of CsI detector response for optimization of 3D x-ray breast imaging systems," Med. Phys. 37(6), 2593-2605 (2010)] that describes detector response functions for indirect detectors and oblique incident x rays. The model, modified in this work to address direct detector response, describes the dependence of the response with x-ray energy, thickness of the transducer layer, and the depth-dependent blur and collection efficiency.

RESULTS

The authors report the detector response functions for indirect and direct detector models for typical thicknesses utilized in clinical systems for full-field digital mammography (150 microm for indirect CsI:Tl and 200 microm for a-Se direct detectors). The results suggest that the oblique incidence effect in a semiconductor detector differs from that in indirect detectors in two ways: The direct detector model produces a sharper overall PRF compared to the response corresponding to the indirect detector model for normal x-ray incidence and a larger relative increase in blur along the x-ray incidence direction compared to that found in indirect detectors with respect to the response at normal incidence angles.

CONCLUSIONS

Compared to the effect seen in indirect detectors, the direct detector model exhibits a sharper response at normal x-ray incidence and a larger relative increase in blur along the x-ray incidence direction with respect to the blur in the orthogonal direction. The results suggest that the oblique incidence effect in direct detectors can be considered to be caused mostly by the geometry of the path where the x-ray beam and its secondary particles deposit energy in the semiconductor layer.

Penalized maximum likelihood reconstruction for improved microcalcification detection in breast tomosynthesis

We examined the application of an iterative penalized maximum likelihood (PML) reconstruction method for improved detectability of microcalcifications (MCs) in digital breast tomosynthesis (DBT). Localized receiver operating characteristic (LROC) psychophysical studies with human observers and 2-D image slices were conducted to evaluate the performance of this reconstruction method and to compare its performance against the commonly used Feldkamp FBP algorithm. DBT projections were generated using rigorous computer simulations that included accurate modeling of the noise and detector blur. Acquisition dose levels of 0.7, 1.0, and 1.5 mGy in a 5-cm-thick compressed breast were tested. The defined task was to localize and detect MC clusters consisting of seven MCs. The individual MC diameter was 150 μm. Compressed-breast phantoms derived from CT images of actual mastectomy specimens provided realistic background structures for the detection task. Four observers each read 98 test images for each combination of reconstruction method and acquisition dose. All observers performed better with the PML images than with the FBP images. With the acquisition dose of 0.7 mGy, the average areas under the LROC curve (A(L)) for the PML and FBP algorithms were 0.69 and 0.43, respectively. For the 1.0-mGy dose, the values of A(L) were 0.93 (PML) and 0.7 (FBP), while the 1.5-mGy dose resulted in areas of 1.0 and 0.9, respectively, for the PML and FBP algorithms. A 2-D analysis of variance applied to the individual observer areas showed statistically significant differences (at a significance level of 0.05) between the reconstruction strategies at all three dose levels. There were no significant differences in observer performance for any of the dose levels.

A comparative analysis of OTF, NPS, and DQE in energy integrating and photon counting digital x-ray detectors

PURPOSE

One of the benefits of photon counting (PC) detectors over energy integrating (EI) detectors is the absence of many additive noise sources, such as electronic noise and secondary quantum noise. The purpose of this work is to demonstrate that thresholding voltage gains to detect individual x rays actually generates an unexpected source of white noise in photon counters.

METHODS

To distinguish the two detector types, their point spread function (PSF) is interpreted differently. The PSF of the energy integrating detector is treated as a weighting function for counting x rays, while the PSF of the photon counting detector is interpreted as a probability. Although this model ignores some subtleties of real imaging systems, such as scatter and the energy-dependent amplification of secondary quanta in indirect-converting detectors, it is useful for demonstrating fundamental differences between the two detector types. From first principles, the optical transfer function (OTF) is calculated as the continuous Fourier transform of the PSF, the noise power spectra (NPS) is determined by the discrete space Fourier transform (DSFT) of the autocovariance of signal intensity, and the detective quantum efficiency (DQE) is found from combined knowledge of the OTF and NPS. To illustrate the calculation of the transfer functions, the PSF is modeled as the convolution of a Gaussian with the product of rect functions. The Gaussian reflects the blurring of the x-ray converter, while the rect functions model the sampling of the detector.

RESULTS

The transfer functions are first calculated assuming outside noise sources such as electronic noise and secondary quantum noise are negligible. It is demonstrated that while OTF is the same for two detector types possessing an equivalent PSF, a frequency-independent (i.e., "white") difference in their NPS exists such that NPS(PC) > or = NPS(EI) and hence DQE(PC) < or = DQE(EI). The necessary and sufficient condition for equality is that the PSF is a binary function given as zero or unity everywhere. In analyzing the model detector with Gaussian blurring, the difference in NPS and DQE between the two detector types is found to increase with the blurring of the x-ray converter. Ultimately, the expression for the additive white noise of the photon counter is compared against the expression for electronic noise and secondary quantum noise in an energy integrator. Thus, a method is provided to determine the average secondary quanta that the energy integrator must produce for each x ray to have superior DQE to a photon counter with the same PSF.

CONCLUSIONS

This article develops analytical models of OTF, NPS, and DQE for energy integrating and photon counting digital x-ray detectors. While many subtleties of real imaging systems have not been modeled, this work is illustrative in demonstrating an additive source of white noise in photon counting detectors which has not yet been described in the literature. One benefit of this analysis is a framework for determining the average secondary quanta that an energy integrating detector must produce for each x ray to have superior DQE to competing photon counting technology.

A statistical, task-based evaluation method for three-dimensional x-ray breast imaging systems using variable-background phantoms

PURPOSE

For the last few years, development and optimization of three-dimensional (3D) x-ray breast imaging systems, such as digital breast tomosynthesis (DBT) and computed tomography, have drawn much attention from the medical imaging community, either academia or industry. However, there is still much room for understanding how to best optimize and evaluate the devices over a large space of many different system parameters and geometries. Current evaluation methods, which work well for 2D systems, do not incorporate the depth information from the 3D imaging systems. Therefore, it is critical to develop a statistically sound evaluation method to investigate the usefulness of inclusion of depth and background-variability information into the assessment and optimization of the 3D systems.

METHODS

In this paper, we present a mathematical framework for a statistical assessment of planar and 3D x-ray breast imaging systems. Our method is based on statistical decision theory, in particular, making use of the ideal linear observer called the Hotelling observer. We also present a physical phantom that consists of spheres of different sizes and materials for producing an ensemble of randomly varying backgrounds to be imaged for a given patient class. Lastly, we demonstrate our evaluation method in comparing laboratory mammography and three-angle DBT systems for signal detection tasks using the phantom's projection data. We compare the variable phantom case to that of a phantom of the same dimensions filled with water, which we call the uniform phantom, based on the performance of the Hotelling observer as a function of signal size and intensity.

RESULTS

Detectability trends calculated using the variable and uniform phantom methods are different from each other for both mammography and DBT systems.

CONCLUSIONS

Our results indicate that measuring the system's detection performance with consideration of background variability may lead to differences in system performance estimates and comparisons. For the assessment of 3D systems, to accurately determine trade offs between image quality and radiation dose, it is critical to incorporate randomness arising from the imaging chain including background variability into system performance calculations.

A fast, angle-dependent, analytical model of CsI detector response for optimization of 3D x-ray breast imaging systems

PURPOSE

Accurate models of detector blur are crucial for performing meaningful optimizations of three-dimensional (3D) x-ray breast imaging systems as well as for developing reconstruction algorithms that faithfully reproduce the imaged object anatomy. So far, x-ray detector blur has either been ignored or modeled as a shift-invariant symmetric function for these applications. The recent development of a Monte Carlo simulation package called MANTIS has allowed detailed modeling of these detector blur functions and demonstrated the magnitude of the anisotropy for both tomosynthesis and breast CT imaging systems. Despite the detailed results that MANTIS produces, the long simulation times required make inclusion of these results impractical in rigorous optimization and reconstruction algorithms. As a result, there is a need for detector blur models that can be rapidly generated.

METHODS

In this study, the authors have derived an analytical model for deterministic detector blur functions, referred to here as point response functions (PRFs), of columnar CsI phosphor screens. The analytical model is x-ray energy and incidence angle dependent and draws on results from MANTIS to indirectly include complicated interactions that are not explicitly included in the mathematical model. Once the mathematical expression is derived, values of the coefficients are determined by a two-dimensional (2D) fit to MANTIS-generated results based on a figure-of-merit (FOM) that measures the normalized differences between the MANTIS and analytical model results averaged over a region of interest. A smaller FOM indicates a better fit. This analysis was performed for a monochromatic x-ray energy of 25 keV, a CsI scintillator thickness of 150 microm, and four incidence angles (0 degrees, 15 degrees, 30 degrees, and 45 degrees).

RESULTS

The FOMs comparing the analytical model to MANTIS for these parameters were 0.1951 +/- 0.0011, 0.1915 +/- 0.0014, 0.2266 +/- 0.0021, and 0.2416 +/- 0.0074 for 0 degrees, 15 degrees, 30 degrees, and 45 degrees, respectively. As a comparison, the same FOMs comparing MANTIS to 2D symmetric Gaussian fits to the zero-angle PRF were 0.6234 +/- 0.0020, 0.9058 +/- 0.0029, 1.491 +/- 0.012, and 2.757 +/- 0.039 for the same set of incidence angles. Therefore, the analytical model matches MANTIS results much better than a 2D symmetric Gaussian function. A comparison was also made against experimental data for a 170 microm thick CsI screen and an x-ray energy of 25.6 keV. The corresponding FOMs were 0.3457 +/- 0.0036, 0.3281 +/- 0.0057, 0.3422 +/- 0.0023, and 0.3677 +/- 0.0041 for 0 degrees, 15 degrees, 30 degrees, and 45 degrees, respectively. In a previous study, FOMs comparing the same experimental data to MANTIS PRFs were found to be 0.2944 +/- 0.0027, 0.2387 +/- 0.0039, 0.2816 +/- 0.0025, and 0.2665 +/- 0.0032 for the same set of incidence angles.

CONCLUSIONS

The two sets of derived FOMs, comparing MANTIS-generated PRFs and experimental data to the analytical model, demonstrate that the analytical model is able to reproduce experimental data with a FOM of less than two times that comparing MANTIs and experimental data. This performance is achieved in less than one millionth the computation time required to generate a comparable PRF with MANTIS. Such small computation times will allow for the inclusion of detailed detector physics in rigorous optimization and reconstruction algorithms for 3D x-ray breast imaging systems.

Task-based assessment of breast tomosynthesis: effect of acquisition parameters and quantum noise

PURPOSE

Tomosynthesis is a promising modality for breast imaging. The appearance of the tomosynthesis reconstructed image is greatly affected by the choice of acquisition and reconstruction parameters. The purpose of this study was to investigate the limitations of tomosynthesis breast imaging due to scan parameters and quantum noise. Tomosynthesis image quality was assessed based on performance of a mathematical observer model in a signal-known exactly (SKE) detection task.

METHODS

SKE detectability (d') was estimated using a prewhitening observer model. Structured breast background was simulated using filtered noise. Detectability was estimated for designer nodules ranging from 0.05 to 0.8 cm in diameter. Tomosynthesis slices were reconstructed using iterative maximum-likelihood expectation-maximization. The tomosynthesis scan angle was varied between 15 degrees and 60 degrees, the number of views between 11 and 41 and the total number of x-ray quanta was infinity, 6 X 10(5), and 6 x 10(4). Detectability in tomosynthesis was compared to that in a single projection.

RESULTS

For constant angular sampling distance, increasing the angular scan range increased detectability for all signal sizes. Large-scale signals were little affected by quantum noise or angular sampling. For small-scale signals, quantum noise and insufficient angular sampling degraded detectability. At high quantum noise levels, angular step size of 3 degrees or below was sufficient to avoid image degradation. At lower quantum noise levels, increased angular sampling always resulted in increased detectability. The ratio of detectability in the tomosynthesis slice to that in a single projection exhibited a peak that shifted to larger signal sizes when the angular range increased. For a given angular range, the peak shifted toward smaller signals when the number of views was increased. The ratio was greater than unity for all conditions evaluated.

CONCLUSION

The effect of acquisition parameters on lesion detectability depends on signal size. Tomosynthesis scan angle had an effect on detectability for all signals sizes, while quantum noise and angular sampling only affected the detectability small-scale signals.

Calculation of OTF, NPS, and DQE for Oblique X-Ray Incidence on Turbid Granular Phosphors

Digital breast tomosynthesis (DBT) is an imaging modality in which tomographic sections of the breast are generated from a limited range of x-ray tube angles. One drawback of DBT is resolution loss in the oblique projection images. The purpose of this work is to extend Swank's formulation of the transfer functions of turbid granular phosphors to oblique x-ray incidence, using the diffusion approximation to the Boltzmann equation to model the spread of light in the phosphor. As expected, the modulation transfer function (MTF) and noise power spectra (NPS) are found to decrease with projection angle regardless of frequency. By contrast, the dependence of detective quantum efficiency (DQE) on projection angle is frequency dependent. DQE increases with projection angle at low frequencies, and only decreases with projection angle at high frequencies. Importantly, the x-ray quantum detection efficiency (AQ ) and the Swank information factor (AS ) are also found to be angularly dependent.

Experimental validation of Monte Carlo (MANTIS) simulated x-ray response of columnar CsI scintillator screens

PURPOSE

MANTIS is a Monte Carlo code developed for the detailed simulation of columnar CsI scintillator screens in x-ray imaging systems. Validation of this code is needed to provide a reliable and valuable tool for system optimization and accurate reconstructions for a variety of x-ray applications. Whereas previous validation efforts have focused on matching of summary statistics, in this work the authors examine the complete point response function (PRF) of the detector system in addition to relative light output values.

METHODS

Relative light output values and high-resolution PRFs have been experimentally measured with a custom setup. A corresponding set of simulated light output values and PRFs have also been produced, where detailed knowledge of the experimental setup and CsI:Tl screen structures are accounted for in the simulations. Four different screens were investigated with different thicknesses, column tilt angles, and substrate types. A quantitative comparison between the experimental and simulated PRFs was performed for four different incidence angles (0 degrees, 15 degrees, 30 degrees, and 45 degrees) and two different x-ray spectra (40 and 70 kVp). The figure of merit (FOM) used measures the normalized differences between the simulated and experimental data averaged over a region of interest.

RESULTS

Experimental relative light output values ranged from 1.456 to 1.650 and were in approximate agreement for aluminum substrates, but poor agreement for graphite substrates. The FOMs for all screen types, incidence angles, and energies ranged from 0.1929 to 0.4775. To put these FOMs in context, the same FOM was computed for 2D symmetric Gaussians fit to the same experimental data. These FOMs ranged from 0.2068 to 0.8029. Our analysis demonstrates that MANTIS reproduces experimental PRFs with higher accuracy than a symmetric 2D Gaussian fit to the experimental data in the majority of cases. Examination of the spatial distribution of differences between the PRFs shows that the main reason for errors between MANTIS and the experimental data is that MANTIS-generated PRFs are sharper than the experimental PRFs.

CONCLUSIONS

The experimental validation of MANTIS performed in this study demonstrates that MANTIS is able to reliably predict experimental PRFs, especially for thinner screens, and can reproduce the highly asymmetric shape seen in the experimental data. As a result, optimizations and reconstructions carried out using MANTIS should yield results indicative of actual detector performance. Better characterization of screen properties is necessary to reconcile the simulated light output values with experimental data.

Optimization of the acquisition geometry in digital tomosynthesis of the breast

Digital tomosynthesis of the breast continues to be intensively studied as an adjunct or replacement of conventional mammography. Although many of the acquisition parameters found in tomosynthesis imaging are also found in conventional mammography and therefore most of the traditional values from mammography have been used in the former, two acquisition geometry parameters, the angular range covered during acquisition and the number of projections per projection set, are applicable only to tomosynthesis. Therefore, in the preclinical and clinical studies reported on tomosynthesis of the breast, a wide variety of values have been used for these two parameters. In this study, 63 different combinations of angular range and number of projections were evaluated using computer simulation methods to characterize how these two parameters affect reconstruction quality and to find which of these combinations maximize it. For this, a computer simulation of a digital tomosynthesis system that included empirically determined system response characteristics was developed and used to generate 9450 different breast tissue volume reconstructions. These reconstructions were analyzed objectively using metrics for in-plane lesion visibility and vertical resolution in the form of the contrast-to-noise ratio and artifact spread function (ASF). It was found that although maximizing the angular range covered does always increase the vertical resolution in tomosynthesis, increasing the number of projections in the acquisition set beyond a relatively low threshold does not further improve the vertical resolution. This threshold value for the minimal number of projections needed to minimize the ASF was found to vary proportionally with angular range. For example, for a 60 degrees angular range, the threshold number of projections was found to be 13. Given the clear inverse relationship between the number of projections and in-plane reconstruction quality under a limited total glandular dose condition, the optimum acquisition geometry in tomosynthesis imaging of the breast is that which maximizes the angular range while maintaining the number of projections close to the threshold values found. Finally, of the 63 acquisition geometries studied, the one that resulted in the highest quality reconstruction, considering both in-plane quality and vertical resolution, consisted of the acquisition of 13 projections over a 60 degrees angular range.

Effect of oblique X-ray incidence in flat-panel computed tomography of the breast

We quantify the variation in resolution due to anisotropy caused by oblique X-ray incidence in indirect flat-panel detectors for computed tomography breast imaging systems. We consider a geometry and detector type utilized in breast computed tomography (CT) systems currently being developed. Our methods rely on mantis, a combined X-ray, electron, and optical Monte Carlo transport open source code. The physics models are the most accurate available in general-purpose Monte Carlo packages in the diagnostic energy range. We consider maximum-obliquity angles of 10 ( degrees ) and 13 ( degrees ) at the centers of the 30 and 40 cm detector edges, respectively, and 16 ( degrees ) at the corner of the detector. Our results indicate that blur is asymmetric and that the resolution properties vary significantly with the angle (or location) of incidence. Our results suggest that the asymmetry can be as high as a factor of 2.6 between orthogonal directions. Anisotropy maps predicted by mantis provide an understanding of the effect that such variations have on the imaging system and allow more accurate modeling and optimization of breast CT systems. These maps of anisotropy across the detector could lead to improved reconstruction and help motivate physics-based strategies for computer detection of breast lesions.

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.

Investigation of physical image characteristics and phenomenon of edge enhancement by phase contrast using equipment typical for mammography

A technique called phase contrast mammography (PCM) has only recently been applied in clinical examination. In this application, PCM images are acquired at a 1.75 x magnification using an x-ray tube for clinical use, and then reduced to the real size of the object by image processing. The images showed enhanced object edges; reportedly, this enhancement occurred because of the refraction of x rays through a cylindrical object. The authors measured the physical image characteristics of PCM to compare the image characteristics of PCM with those of conventional mammography. More specifically, they measured the object-edge-response characteristics and the noise characteristics in the spatial frequency domain. The results revealed that the edge-response characteristics of PCM outperformed those of conventional mammography. In addition, the characteristics changed with the object-placement conditions and the object shapes. The noise characteristics of PCM were better than those of conventional mammography. Subsequently, to verify why object edges were enhanced in PCM images, the authors simulated image profiles that would be obtained if the x rays were refracted and totally reflected by using not only a cylindrical substance but also a planar substance as the object. So, they confirmed that the object edges in PCM images were enhanced because x rays were refracted irrespective of the object shapes. Further, they found that the edge enhancements depended on the object shapes and positions. It was also proposed that the larger magnification than 1.75 in the commercialized system might be more suitable for PCM. Finally, the authors investigated phase-contrast effects to breast tissues by the simulation and demonstrated that PCM would be helpful in the diagnoses of mammography.

Fundamental x-ray interaction limits in diagnostic imaging detectors: spatial resolution

The practice of diagnostic x-ray imaging has been transformed with the emergence of digital detector technology. Although digital systems offer many practical advantages over conventional film-based systems, their spatial resolution performance can be a limitation. The authors present a Monte Carlo study to determine fundamental resolution limits caused by x-ray interactions in four converter materials: Amorphous silicon (a-Si), amorphous selenium, cesium iodide, and lead iodide. The "x-ray interaction" modulation transfer function (MTF) was determined for each material and compared in terms of the 50% MTF spatial frequency and Wagner's effective aperture for incident photon energies between 10 and 150 keV and various converter thicknesses. Several conclusions can be drawn from their Monte Carlo study. (i) In low-Z (a-Si) converters, reabsorption of Compton scatter x rays limits spatial resolution with a sharp MTF drop at very low spatial frequencies (< 0.3 cycles/mm), especially above 60 keV; while in high-Z materials, reabsorption of characteristic x rays plays a dominant role, resulting in a mid-frequency (1-5 cycles/mm) MTF drop. (ii) Coherent scatter plays a minor role in the x-ray interaction MTF. (iii) The spread of energy due to secondary electron (e.g., photoelectrons) transport is significant only at very high spatial frequencies. (iv) Unlike the spread of optical light in phosphors, the spread of absorbed energy from x-ray interactions does not significantly degrade spatial resolution as converter thickness is increased. (v) The effective aperture results reported here represent fundamental spatial resolution limits of the materials tested and serve as target benchmarks for the design and development of future digital x-ray detectors.

The x-ray light valve: a low-cost, digital radiographic imaging system--spatial resolution

An x-ray light valve (XLV) coupled with an optical scanner has the potential to meet the need for a low-cost, high quality digital imaging system for general radiography. The XLV/scanner concept combines three well-established, and hence, low-cost technologies: An amorphous selenium (a-Se) layer as an x-ray-to-charge transducer, a liquid crystal (LC) cell as an analog display, and an optical scanner for image digitization. The XLV consists of an a-Se layer and LC cell in a sandwich structure which produces an optical image in the LC layer upon x-ray exposure. The XLV/scanner system consists of an XLV in combination with an optical scanner for image readout. Here, the effect of each component on the spatial resolution of an XLV/scanner system is investigated. A theoretical model of spatial resolution of an XLV is presented based on calculations of the modulation transfer function (MTF) for a-Se and a LC cell. From these component MTFs, the theoretical MTF of the XLV is derived. The model was validated by experiments on a prototype XLV/scanner system. The MTF of the scanner alone was obtained by scanning an optical test target and the MTF of the XLV/scanner system was measured using x rays. From the measured MTF of the scanner, the theoretical MTF of the XLV/scanner system was established and compared with the experimental results. Good general agreement exists between experimental and theoretical results in the frequency range of interest for general radiography, although the theoretical curves slightly overstate the measured MTFs. The experimental MTF of the XLV was compared with the MTF of two clinical systems and was shown to have the capability to exceed the resolution of flat-panel detectors. From this, the authors can conclude that the XLV has an adequate resolution for general radiography. The XLV/scanner also has the potential to eliminate aliasing while maintaining a MTF that exceeds that of a flat-panel imager.

Anniversary paper. Development of x-ray computed tomography: the role of medical physics and AAPM from the 1970s to present

The AAPM, through its members, meetings, and its flagship journal Medical Physics, has played an important role in the development and growth of x-ray tomography in the last 50 years. From a spate of early articles in the 1970s characterizing the first commercial computed tomography (CT) scanners through the "slice wars" of the 1990s and 2000s, the history of CT and related techniques such as tomosynthesis can readily be traced through the pages of Medical Physics and the annals of the AAPM and RSNA/AAPM Annual Meetings. In this article, the authors intend to give a brief review of the role of Medical Physics and the AAPM in CT and tomosynthesis imaging over the last few decades.

Imaging performance of an amorphous selenium digital mammography detector in a breast tomosynthesis system

In breast tomosynthesis a rapid sequence of N images is acquired when the x-ray tube sweeps through different angular views with respect to the breast. Since the total dose to the breast is kept the same as that in regular mammography, the exposure used for each image of tomosynthesis is 1/N. The low dose and high frame rate pose a tremendous challenge to the imaging performance of digital mammography detectors. The purpose of the present work is to investigate the detector performance in different operational modes designed for tomosynthesis acquisition, e.g., binning or full resolution readout, the range of view angles, and the number of views N. A prototype breast tomosynthesis system with a nominal angular range of +/-25 degrees was used in our investigation. The system was equipped with an amorphous selenium (a-Se) full field digital mammography detector with pixel size of 85 microm. The detector can be read out in full resolution or 2 x 1 binning (binning in the tube travel direction). The focal spot blur due to continuous tube travel was measured for different acquisition geometries, and it was found that pixel binning, instead of focal spot blur, dominates the detector modulation transfer function (MTF). The noise power spectrum (NPS) and detective quantum efficiency (DQE) of the detector were measured with the exposure range of 0.4-6 mR, which is relevant to the low dose used in tomosynthesis. It was found that DQE at 0.4 mR is only 20% less than that at highest exposure for both detector readout modes. The detector temporal performance was categorized as lag and ghosting, both of which were measured as a function of x-ray exposure. The first frame lags were 8% and 4%, respectively, for binning and full resolution mode. Ghosting is negligible and independent of the frame rate. The results showed that the detector performance is x-ray quantum noise limited at the low exposures used in each view of tomosynthesis, and the temporal performance at high frame rate (up to 2 frames per second) is adequate for tomosynthesis.

Automated detection of microcalcification clusters for digital breast tomosynthesis using projection data only: a preliminary study

Digital breast tomosynthesis (DBT) is a promising modality for breast imaging in which an anisotropic volume image of the breast is obtained. We present an algorithm for computerized detection of microcalcification clusters (MCCs) for DBT. This algorithm operates on the projection views only. Therefore it does not depend on reconstruction, and is computationally efficient. The algorithm was developed using a database of 30 image sets with microcalcifications, and a control group of 30 image sets without visible findings. The patient data were acquired on the first DBT prototype at Massachusetts General Hospital. Algorithm sensitivity was estimated to be 0.86 at 1.3 false positive clusters, which is below that of current MCC detection algorithms for full-field digital mammography. Because of the small number of patient cases, algorithm parameters were not optimized and one linear classifier was used. An actual limitation of our approach may be that the signal-to-noise ratio in the projection images is too low for microcalcification detection. Furthermore, the database consisted of predominantly small MCC. This may be related to the image quality obtained with this first prototype.