Observation of super-resolution in digital breast tomosynthesis

Exploring advanced 2D acquisitions in breast tomosynthesis: T-shaped and Pentagon geometries

In this study, we investigate the performance of advanced 2D acquisition geometries - Pentagon and T-shaped - in digital breast tomosynthesis (DBT) and compare them against the conventional 1D geometry. Unlike the conventional approach, our proposed 2D geometries also incorporate anterior projections away from the chest wall. Implemented on the Next-Generation Tomosynthesis (NGT) prototype developed by X-ray Physics Lab (XPL), UPenn, we utilized various phantoms to compare three geometries: a Defrise slab phantom with alternating plastic slabs to study low-frequency modulation; a Checkerboard breast phantom (a 2D adaptation of the Defrise phantom design) to study the ability to reconstruct the fine features of the checkerboard squares; and the 360° Star-pattern phantom to assess aliasing and compute the Fourier-spectral distortion (FSD) metric that assesses spectral leakage and the contrast transfer function. We find that both Pentagon and T-shaped scans provide greater modulation amplitude of the Defrise phantom slabs and better resolve the squares of the Checkerboard phantom against the conventional scan. Notably, the Pentagon geometry exhibited a significant reduction in aliasing of spatial frequencies oriented in the right-left (RL) medio-lateral direction, which was corroborated by a near complete elimination of spectral leakage in the FSD plot. Conversely T-shaped scan redistributes the aliasing between both posteroanterior (PA) and RL directions thus maintaining non-inferiority against the conventional scan which is predominantly affected by PA aliasing. The results of this study underscore the potential of incorporating advanced 2D geometries in DBT systems, offering marked improvements in imaging performance over the conventional 1D approach.

Impact of super-resolution and image acquisition on the detection of calcifications in digital breast tomosynthesis

OBJECTIVES

A virtual clinical trial (VCT) method is proposed to determine the limit of calcification detection in tomosynthesis.

METHODS

Breast anatomy, focal findings, image acquisition, and interpretation (n = 14 readers) were simulated using screening data (n = 660 patients). Calcifications (0.2-0.4 mm3) were inserted into virtual breast phantoms. Digital breast tomosynthesis (DBT) acquisitions were simulated assuming various acquisition geometries: source motion (continuous and step-and-shoot), detector element size (140 and 70 µm), and reconstructed voxel size (35-140 µm). VCT results were estimated using multiple-reader multiple-case analyses and d' statistics. Signal-to-noise (SNR) analyses were also performed using BR3D phantoms.

RESULTS

Source motion and reconstructed voxel size demonstrated significant changes in the performance of imaging systems. Acquisition geometries that use 70 µm reconstruction voxel size and step-and-shoot motion significantly improved calcification detection. Comparing 70 with 100 µm reconstructed voxel size for step-and-shoot, the ΔAUC was 0.0558 (0.0647) and d' ratio was 1.27 (1.29) for 140 µm (70 µm) detector element size. Comparing step-and-shoot with a continuous motion for a 70 µm reconstructed voxel size, the ΔAUC was 0.0863 (0.0434) and the d' ratio was 1.40 (1.19) for 140 µm (70 µm) detector element. Small detector element sizes (e.g., 70 µm) did not significantly improve detection. The SNR results with the BR3D phantom show that calcification detection is dependent upon reconstructed voxel size and detector element size, supporting VCT results with comparable agreement (ratios: d' = 1.16 ± 0.11, SNR = 1.34 ± 0.13).

CONCLUSION

DBT acquisition geometries that use super-resolution (smaller reconstructed voxels than the detector element size) combined with step-and-shoot motion have the potential to improve the detection of calcifications.

CLINICAL RELEVANCE

Calcifications may not always be discernable in tomosynthesis because of differences in acquisition and reconstruction methods. VCTs can identify strategies to optimize acquisition and reconstruction parameters for calcification detection in tomosynthesis, most notably through super-resolution in the reconstruction.

KEY POINTS

• Super-resolution improves calcification detection and SNR in tomosynthesis; specifically, with the use of smaller reconstruction voxels. • Calcification detection using step-and-shoot motion is superior to that using continuous tube motion. • A detector element size of 70 µm does not provide better detection than 140 µm for small calcifications at the threshold of detectability.

Non-Isocentric Geometry for Next-Generation Tomosynthesis With Super-Resolution

Our lab at the University of Pennsylvania (UPenn) is investigating novel designs for digital breast tomosynthesis. We built a next-generation tomosynthesis system with a non-isocentric geometry (superior-to-inferior detector motion). This paper examines four metrics of image quality affected by this design. First, aliasing was analyzed in reconstructions prepared with smaller pixelation than the detector. Aliasing was assessed with a theoretical model of r -factor, a metric calculating amplitudes of alias signal relative to input signal in the Fourier transform of the reconstruction of a sinusoidal object. Aliasing was also assessed experimentally with a bar pattern (illustrating spatial variations in aliasing) and 360°-star pattern (illustrating directional anisotropies in aliasing). Second, the point spread function (PSF) was modeled in the direction perpendicular to the detector to assess out-of-plane blurring. Third, power spectra were analyzed in an anthropomorphic phantom developed by UPenn and manufactured by Computerized Imaging Reference Systems (CIRS), Inc. (Norfolk, VA). Finally, calcifications were analyzed in the CIRS Model 020 BR3D Breast Imaging Phantom in terms of signal-to-noise ratio (SNR); i.e., mean calcification signal relative to background-tissue noise. Image quality was generally superior in the non-isocentric geometry: Aliasing artifacts were suppressed in both theoretical and experimental reconstructions prepared with smaller pixelation than the detector. PSF width was also reduced at most positions. Anatomic noise was reduced. Finally, SNR in calcification detection was improved. (A potential trade-off of smaller-pixel reconstructions was reduced SNR; however, SNR was still improved by the detector-motion acquisition.) In conclusion, the non-isocentric geometry improved image quality in several ways.

Impact of Tomosynthesis Acquisition on 3D Segmentations of Breast Outline and Adipose/Dense Tissue with AI: A Simulation-Based Study

Digital breast tomosynthesis (DBT) reconstructions introduce out-of-plane artifacts and false-tissue boundaries impacting the dense/adipose and breast outline (convex hull) segmentations. A virtual clinical trial method was proposed to segment both the breast tissues and the breast outline in DBT reconstructions. The DBT images of a representative population were simulated using three acquisition geometries: a left-right scan (conventional, I), a two-directional scan in the shape of a "T" (II), and an extra-wide range (XWR, III) left-right scan at a six-times higher dose than I. The nnU-Net was modified including two losses for segmentation: (1) tissues and (2) breast outline. The impact of loss (1) and the combination of loss (1) and (2) was evaluated using models trained with data simulating geometry I. The impact of the geometry was evaluated using the combined loss (1&2). The loss (1&2) improved the convex hull estimates, resolving 22.2% of the false classification of air voxels. Geometry II was superior to I and III, resolving 99.1% and 96.8% of the false classification of air voxels. Geometry III (Dice = (0.98, 0.94)) was superior to I (0.92, 0.78) and II (0.93, 0.74) for the tissue segmentation (adipose, dense, respectively). Thus, the loss (1&2) provided better segmentation, and geometries T and XWR improved the dense/adipose and breast outline segmentations relative to the conventional scan.

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.

Achieving Isotropic Super-Resolution with a Non-Isocentric Acquisition Geometry in a Next-Generation Tomosynthesis System

We have constructed a prototype next-generation tomosynthesis (NGT) system that supports a non-isocentric acquisition geometry for digital breast tomosynthesis (DBT). In this geometry, the detector gradually descends in the superior-to-inferior direction. The aim of this work is to demonstrate that this geometry offers isotropic super-resolution (SR), unlike clinical DBT systems which are characterized by anisotropies in SR. To this end, a theoretical model of a sinusoidal test object was developed with frequency exceeding the alias frequency of the detector. We simulated two geometries: (1) a conventional geometry with a stationary detector, and (2) a non-isocentric geometry. The input frequency was varied over the full 360° range of angles in the plane of the object. To investigate whether SR was achieved, we calculated the Fourier transform of the reconstruction. The amplitude of the tallest peak below the alias frequency was measured relative to the peak at the input frequency. This ratio (termed the r-factor) should approach zero to achieve high-quality SR. In the conventional geometry, the r-factor was minimized (approaching zero) if the orientation of the frequency was parallel with the source motion, yet exceeded unity (prohibiting SR) in the orientation perpendicular to the source motion. However, in the non-isocentric geometry, the r-factor was minimized (approaching zero) for all orientations of the frequency, meaning SR was achieved isotropically. In summary, isotropic SR in DBT can be achieved using the non-isocentric acquisition geometry supported by the NGT system.

Effect of Dose Level on Radiologists' Detection of Microcalcifications in Digital Breast Tomosynthesis: An Observer Study with Breast Phantoms

OBJECTIVES

To compare radiologists' sensitivity, confidence level, and reading efficiency of detecting microcalcifications in digital breast tomosynthesis (DBT) at two clinically relevant dose levels.

MATERIALS AND METHODS

Six 5-cm-thick heterogeneous breast phantoms embedded with a total of 144 simulated microcalcification clusters of four speck sizes were imaged at two dose modes by a clinical DBT system. The DBT volumes at the two dose levels were read independently by six MQSA radiologists and one fellow with 1-33 years (median 12 years) of experience in a fully-crossed counter-balanced manner. The radiologist located each potential cluster and rated its conspicuity and his/her confidence that the marked location contained a cluster. The differences in the results between the two dose modes were analyzed by two-tailed paired t-test.

RESULTS

Compared to the lower-dose mode, the average glandular dose in the higher-dose mode for the 5-cm phantoms increased from 1.34 to 2.07 mGy. The detection sensitivity increased for all speck sizes and significantly for the two smaller sizes (p <0.05). An average of 13.8% fewer false positive clusters was marked. The average conspicuity rating and the radiologists' confidence level were higher for all speck sizes and reached significance (p <0.05) for the three larger sizes. The average reading time per detected cluster reduced significantly (p <0.05) by an average of 13.2%.

CONCLUSION

For a 5-cm-thick breast, an increase in average glandular dose from 1.34 to 2.07 mGy for DBT imaging increased the conspicuity of microcalcifications, improved the detection sensitivity by radiologists, increased their confidence levels, reduced false positive detections, and increased the reading efficiency.

Estimated risk of radiation-induced cancer following breast screening employing tomosynthesis and digital mammography

The objective of this study was to estimate the risk of radiation-induced injury to the organs due to ionizing radiation following breast screening recommendations employing digital breast tomosynthesis (DBT) and digital mammography (DM). Using the Monte Carlo method, absorbed doses in the tissues and organs were calculated on an adult female phantom, considering two-view craniocaudal (CC) and mediolateral oblique (MLO) projections for each breast. The results showed differences in the total effective risk due to DM (CC + MLO) and DBT (CC + MLO) examinations in Brazil, ranging from 20.73 cases 10^-5(DM) to 27.19 cases 10^-5(DBT). Significant differences were also observed in the total effective risk of cancer incidence in the lungs due to DM (CC + MLO) and DBT (CC + MLO), ranging from 1.75×10^-01cases 10^-5(DM) to 1.76×10^-01cases 10^-5(DBT). The results indicate that the total effective risk of incidence should be considered as an additional parameter for the evaluation of DBT or DBT + DM program performance.

Development and Evaluation of the Fourier Spectral Distortion Metric

A spatial resolution metric is presented for tomosynthesis. The Fourier spectral distortion metric (FSD) was developed to evaluate specific resolution properties of different imaging techniques for digital tomosynthesis using a star pattern image to plot modulation in the frequency domain. The FSD samples the spatial resolution of a star-pattern image tangentially over an acute angle and for a range of spatial frequencies in a 2D image or 3D image reconstruction slice. The FSD graph portrays all frequencies present in a star pattern quadrant. In addition to the fundamental input frequency of the star pattern, the FSD graph shows spectral leakage, square wave harmonics, and residual noise. The contrast transfer function (CTF) is obtained using the FSD graph. The CTF is analogous to the modulation transfer function (MTF), but it is not normalized to unity at zero spatial frequency. Unlike the MTF, this metric separates the fundamental input-frequency from the other signals in the Fourier domain. This metric helps determine optimal image reconstruction parameters, the in-plane limit of spatial resolution with respect to aliased signals, and a threshold criterion for an image to support super resolution and reduce aliasing artifacts. Various sampling parameters were evaluated to optimize this metric and ascertain measurement accuracy. The FSD adequately compares resolution properties of 2D images and 3D image reconstruction slices for various x ray imaging modes without suppressing aliased signals.

Development of Magnification Tomosynthesis for Superior Resolution in Diagnostic Mammography

Our previous work showed that digital breast tomosynthesis (DBT) supports super-resolution (SR). Clinical systems are not yet designed to optimize SR; this can be demonstrated with a high-frequency line-resolution pattern. SR is achieved if frequencies are oriented laterally, but not if frequencies are oriented in the perpendicular direction; i.e., the posteroanterior (PA) direction. We are developing a next-generation tomosynthesis (NGT) prototype with new trajectories for the x-ray source. This system is being designed to optimize SR not just for screening, but also for diagnostic mammography; specifically, for magnification DBT (M-DBT). SR is not achieved clinically in magnification mammography, since the acquisition is 2D. The aim of this study is to investigate SR in M-DBT, and analyze how anisotropies differ from screening DBT (S-DBT). We have a theoretical model of a high-frequency sinusoidal test object. First, a conventional scanning motion (directed laterally) was simulated. In the PA direction, SR was not achieved in either S-DBT or M-DBT. Next, the scanning motion was angled relative to the lateral direction. This motion introduces submillimeter offsets in source positions in the PA direction. Theoretical modeling demonstrated that SR was achieved in M-DBT, but not in S-DBT, in the PA direction. This work shows that, with the use of magnification, anisotropies in SR are more sensitive to small offsets in the source motion, leading to insights into how to design M-DBT systems.

Super-Resolution in Digital Breast Tomosynthesis: Limitations of the Conventional System Design and Strategies for Optimization

Our previous work explored the use of super-resolution as a way to improve the visibility of calcifications in digital breast tomosynthesis. This paper demonstrates that there are anisotropies in super-resolution throughout the reconstruction, and investigates new motion paths for the x-ray tube to suppress these anisotropies. We used a theoretical model of a sinusoidal test object to demonstrate the existence of the anisotropies. In addition, high-frequency test objects were simulated with virtual clinical trial (VCT) software developed for breast imaging. The simulated objects include a lead bar pattern phantom as well as punctate calcifications in a breast-like background. In a conventional acquisition geometry in which the source motion is directed laterally, we found that super-resolution is not achievable if the frequency is oriented in the perpendicular direction (posteroanteriorly). Also, there are positions, corresponding to various slices above the breast support, at which super-resolution is inherently not achievable. The existence of these anisotropies was validated with VCT simulations. At locations predicted by theoretical modeling, the bar pattern phantom showed aliasing, and the spacing between individual calcifications was not properly resolved. To show that super-resolution can be optimized by re-designing the acquisition geometry, we applied our theoretical model to the analysis of new motion paths for the x-ray tube; specifically, motions with more degrees of freedom and with more rapid pulsing (submillimeter spacing) between source positions. These two strategies can be used in combination to suppress the anisotropies in super-resolution.

Proposing Rapid Source Pulsing for Improved Super-Resolution in Digital Breast Tomosynthesis

Our previous work showed that digital breast tomosynthesis (DBT) systems are capable of super-resolution, or subpixel resolution relative to the detector. Using a bar pattern phantom, it is possible to demonstrate that there are anisotropies in super-resolution throughout the reconstruction. These anisotropies are lessened in acquisition geometries with narrow spacing between source positions. This paper demonstrates that by re-arranging the source positions in the scan, the anisotropies can be minimized even further. To this end, a theoretical model of the reconstruction of a high-frequency sinusoidal test object was developed from first principles. We modeled the effect of clustering additional source positions around each conventional source position in fine increments (submillimeter). This design can be implemented by rapidly pulsing the source during a continuous sweep of the x-ray tube. It is shown that it is not possible to eliminate the anisotropies in a conventional DBT system with uniformly-spaced source positions, even if the increments of spacing are narrower than those used clinically. However, super-resolution can be achieved everywhere if the source positions are re-arranged in clusters with submillimeter spacing. Our previous work investigated a different approach for optimizing super-resolution through the use of detector motion perpendicular to the breast support. The advantage of introducing rapid source pulsing is that detector motion is no longer required; this mitigates the need for a thick detector housing, which may be cumbersome for patient positioning.

Effect of source blur on digital breast tomosynthesis reconstruction

PURPOSE

Most digital breast tomosynthesis (DBT) reconstruction methods neglect the blurring of the projection views caused by the finite size or motion of the x-ray focal spot. This paper studies the effect of source blur on the spatial resolution of reconstructed DBT using analytical calculation and simulation, and compares the influence of source blur over a range of blurred source sizes.

METHODS

Mathematically derived formulas describe the point spread function (PSF) of source blur on the detector plane as a function of the spatial locations of the finite-sized source and the object. By using the available technical parameters of some clinical DBT systems, we estimated the effective source sizes over a range of exposure time and DBT scan geometries. We used the CatSim simulation tool (GE Global Research, NY) to generate digital phantoms containing line pairs and beads at different locations and imaged with sources of four different sizes covering the range of potential source blur. By analyzing the relative contrasts of the test objects in the reconstructed images, we studied the effect of the source blur on the spatial resolution of DBT. Furthermore, we simulated a detector that rotated in synchrony with the source about the rotation center and calculated the spatial distribution of the blurring distance in the imaged volume to estimate its influence on source blur.

RESULTS

Calculations demonstrate that the PSF is highly shift-variant, making it challenging to accurately implement during reconstruction. The results of the simulated phantoms demonstrated that a typical finite-sized focal spot (~0.3 mm) will not affect the reconstructed image resolution if the x-ray tube is stationary during data acquisition. If the x-ray tube moves during exposure, the extra blur due to the source motion may degrade image resolution, depending on the effective size of the source along the direction of the motion. A detector that rotates in synchrony with the source does not reduce the influence of source blur substantially.

CONCLUSIONS

This study demonstrates that the extra source blur due to the motion of the x-ray tube during image acquisition substantially degrades the reconstructed image resolution. This effect cannot be alleviated by rotating the detector in synchrony with the source. The simulation results suggest that there are potential benefits of modeling the source blur in image reconstruction for DBT systems using continuous-motion acquisition mode.

Digital Breast Tomosynthesis: Physics, Artifacts, and Quality Control Considerations

As digital breast tomosynthesis (DBT) becomes widely used, radiologists must understand the basic principles of (a) image acquisition, (b) artifacts, and (c) quality control (QC) that are specific to DBT. Standard acquisition parameters common to both full-field digital mammography (FFDM) and DBT are combinations of x-ray tube voltage, current exposure time, and anode target and filter combinations. Image acquisition parameters specific to DBT include tube motion, sweep angle, and number of projections. Continuous tube motion or x-ray emission decreases imaging time but leads to focal spot blurring when compared with step-and-shoot techniques. The sweep angle and number of projections determines resolution. Wider sweep angles allow greater out-of-plane (z-axis) resolution, improving visualization of masses and architecture distortion. A greater number of projections increases in-plane or x-y axis resolution, improving visualization of microcalcifications. Artifacts related to DBT include blurring-ripple, truncation, and loss of skin and superficial tissue resolution. Motion artifacts are difficult to recognize because of inherent out-of-plane blurring. To maintain optimal image quality and an "as low as reasonably achievable" (ALARA) radiation dose, regular QC must be performed. DBT is considered a new imaging modality; therefore, breast imaging facilities are required to obtain a separate certification in addition to that in FFDM, and all personnel (radiologists, technologists, and medical physicists) are mandated to complete initial DBT training and maintain appropriate continuing medical education credits. ^©RSNA, 2019.

Non-stationary model of oblique x-ray incidence in amorphous selenium detectors: I. Point spread function

PURPOSE

In previous work, a theoretical model of the point spread function (PSF) for oblique x-ray incidence in amorphous selenium (a-Se) detectors was proposed. The purpose of this paper is to develop a complementary model that includes two additional features. First, the incidence angle and the directionality of ray incidence are calculated at each position, assuming a divergent x-ray beam geometry. This approach allows the non-stationarity of the PSF to be modeled. Second, this paper develops a framework that is applicable to a digital system, unlike previous work which did not model the presence of a thin-film transistor (TFT) array.

METHODS

At each point on the detector, the incidence angle and the ray incidence direction are determined using ray tracing. Based on these calculations, an existing model for the PSF of the x-ray converter (Med Phys. 1995;22:365-374) is generalized to a non-stationary model. The PSF is convolved with the product of two rectangle functions, which model the sampling of the TFT array. The rectangle functions match the detector element (del) size in two dimensions.

RESULTS

It is shown that the PSF can be calculated in closed form. This solution is used to simulate a digital mammography (DM) system at two x-ray energies (20 and 40 keV). Based on the divergence of the x-ray beam, the direction of ray incidence varies with position. Along this direction, the PSF is broader than the reference rect function matching the del size. The broadening is more pronounced with increasing obliquity. At high energy, the PSF deviates more strongly from the reference rect function, indicating that there is more blurring. In addition, the PSF is calculated along the polar angle perpendicular to the ray incidence direction. For this polar angle, the shape of the PSF is dependent upon whether the ray incidence direction is parallel with the sides of the detector. If the ray incidence direction is parallel with either dimension, the PSF is a perfect rectangle function, matching the del size. However, if the ray incidence direction is at an oblique angle relative to the sides of the detector, the PSF is not rectangular. These results illustrate the non-stationarity of the PSF.

CONCLUSIONS

This paper demonstrates that an existing model of the PSF of a-Se detectors can be generalized to include the effects of non-stationarity and digitization. The PSF is determined in closed form. This solution offers the advantage of shorter computation time relative to approaches that use numerical methods. This model is a tool for simulating a-Se detectors in future work, such as in virtual clinical trials with computational phantoms.

Can Breast Compression Be Reduced in Digital Mammography and Breast Tomosynthesis?

OBJECTIVE

The objective of this study was to investigate the impact of decreasing breast compression during digital mammography and breast tomosynthesis (DBT) on perceived pain and image quality.

MATERIALS AND METHODS

In this two-part study, two groups of women with prior mammograms were recruited. In part 1, subjects were positioned for craniocaudal (CC) and mediolateral oblique (MLO) views, and four levels of compression force were applied to evaluate changes in breast thickness, perceived pain, and relative tissue coverage. No imaging was performed. In part 2, two MLO DBT images of one breast of each patient were acquired at standard and reduced compression. Blurring artifacts and tissue coverage were judged by three breast imaging radiologists, and compression force, breast thickness, relative tissue coverage, and perceived pain were recorded.

RESULTS

Only the first reduction in force was feasible because further reduction resulted in inadequate breast immobilization. Mean force reductions of 48% and 47% for the CC and MLO views, respectively, resulted in a significantly reduced perceived pain level, whereas the thickness of the compressed breast increased by 0.02 cm (CC view) and 0.09 (MLO view, part 1 of the study) and 0.38 cm (MLO view, part 2 of the study), respectively, with no change in tissue coverage or increase in motion blurring.

CONCLUSION

Mammography and DBT acquisitions may be possible using half of the compression force used currently, with a significant and substantial reduction in perceived pain with no clinically significant change in breast thickness and tissue coverage.

Segmented separable footprint projector for digital breast tomosynthesis and its application for subpixel reconstruction

PURPOSE

Digital forward and back projectors play a significant role in iterative image reconstruction. The accuracy of the projector affects the quality of the reconstructed images. Digital breast tomosynthesis (DBT) often uses the ray-tracing (RT) projector that ignores finite detector element size. This paper proposes a modified version of the separable footprint (SF) projector, called the segmented separable footprint (SG) projector, that calculates efficiently the Radon transform mean value over each detector element. The SG projector is specifically designed for DBT reconstruction because of the large height-to-width ratio of the voxels generally used in DBT. This study evaluates the effectiveness of the SG projector in reducing projection error and improving DBT reconstruction quality.

METHODS

We quantitatively compared the projection error of the RT and the SG projector at different locations and their performance in regular and subpixel DBT reconstruction. Subpixel reconstructions used finer voxels in the imaged volume than the detector pixel size. Subpixel reconstruction with RT projector uses interpolated projection views as input to provide adequate coverage of the finer voxel grid with the traced rays. Subpixel reconstruction with the SG projector, however, uses the measured projection views without interpolation. We simulated DBT projections of a test phantom using CatSim (GE Global Research, Niskayuna, NY) under idealized imaging conditions without noise and blur, to analyze the effects of the projectors and subpixel reconstruction without other image degrading factors. The phantom contained an array of horizontal and vertical line pair patterns (1 to 9.5 line pairs/mm) and pairs of closely spaced spheres (diameters 0.053 to 0.5 mm) embedded at the mid-plane of a 5-cm-thick breast tissue-equivalent uniform volume. The images were reconstructed with regular simultaneous algebraic reconstruction technique (SART) and subpixel SART using different projectors. The resolution and contrast of the test objects in the reconstructed images and the computation times were compared under different reconstruction conditions.

RESULTS

The SG projector reduced the projector error by 1 to 2 orders of magnitude at most locations. In the worst case, the SG projector still reduced the projection error by about 50%. In the DBT reconstructed slices parallel to the detector plane, the SG projector not only increased the contrast of the line pairs and spheres but also produced more smooth and continuous reconstructed images, whereas the discrete and sparse nature of the RT projector caused artifacts appearing as patterned noise. For subpixel reconstruction, the SG projector significantly increased object contrast and computation speed, especially for high subpixel ratios, compared with the RT projector implemented with accelerated Siddon's algorithm. The difference in the depth resolution among the projectors is negligible under the conditions studied. Our results also demonstrated that subpixel reconstruction can improve the spatial resolution of the reconstructed images, and can exceed the Nyquist limit of the detector under some conditions.

CONCLUSIONS

The SG projector was more accurate and faster than the RT projector. The SG projector also substantially reduced computation time and improved the image quality for the tomosynthesized images with and without subpixel reconstruction.

Digital breast tomosynthesis: lessons learned from early clinical implementation

The limitations of mammography are well known and are partly related to the fact that with conventional imaging, the three-dimensional volume of the breast is imaged and presented in a two-dimensional format. Because normal breast tissue is similar in x-ray attenuation to some breast cancers, clinically relevant malignancies may be obscured by normal overlapping tissue. In addition, complex areas of normal tissue may be perceived as suspicious. The limitations of two-dimensional breast imaging lead to low sensitivity in detecting some cancers and high false-positive recall rates. Although mammographic screening has been shown to reduce breast cancer deaths by approximately 30%, controversy exists over when and how often screening mammography should occur. Digital breast tomosynthesis (DBT) is rapidly being implemented in breast imaging clinics around the world as early clinical data demonstrate that it may address some of the limitations of conventional mammography. With DBT, multiple low-dose x-ray images are acquired in an arc and reconstructed to create a three-dimensional image, thus minimizing the impact of overlapping breast tissue and improving lesion conspicuity. Early studies of screening DBT have shown decreased false-positive callback rates and increased rates of cancer detection (particularly for invasive cancers), resulting in increased sensitivity and specificity. In our clinical practice, we have completed more than 2 years of using two-view digital mammography combined with two-view DBT for all screening and select diagnostic imaging examinations (over 25,000 patients). Our experience, combined with previously published data, demonstrates that the combined use of DBT and digital mammography is associated with improved outcomes for screening and diagnostic imaging. Online supplemental material is available for this article.

Digital Breast Tomosynthesis in the Diagnostic Setting: Indications and Clinical Applications

Digital breast tomosynthesis (DBT) is an emerging technology used in diagnostic breast imaging to evaluate potential abnormalities. In DBT, the compressed breast tissue is imaged in a quasi-three-dimensional manner by performing a series of low-dose radiographic exposures and using the resultant projection image dataset to reconstruct cross-sectional in-plane images in standard mammographic views. Improved visualization of breast detail at diagnostic DBT allows improved characterization of findings, including normal structures and breast cancer. This technology reduces the summation of overlapping breast tissue, which can mimic breast cancer, and provides improved detail of noncalcified mammographic findings seen in breast cancer. It also assists in lesion localization and determining mammographic extent of disease in women with known or suspected breast cancer. The authors review the potential uses, benefits, and limitations of DBT in the diagnostic setting and discuss how radiologists can best use DBT to characterize lesions, localize potential abnormalities, and evaluate the extent of known or suspected breast cancer. The authors' experience shows that DBT can be implemented effectively in the diagnostic workflow to evaluate and localize potential lesions more efficiently. DBT may potentially replace conventional supplemental mammography at diagnostic workup and obviate ultrasonography in select cases.

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.

Oblique reconstructions in tomosynthesis. I. Linear systems theory

PURPOSE

By convention, slices in a tomosynthesis reconstruction are created on planes parallel to the detector. It has not yet been demonstrated that slices can be generated along oblique directions through the same volume, analogous to multiplanar reconstructions in computed tomography (CT). The purpose of this work is to give a proof-of-principle justification for oblique reconstructions in tomosynthesis, which acquires projection images over a smaller angular range than CT.

METHODS

To investigate the visibility of individual frequencies in an oblique reconstruction, a theoretical framework is developed in which the reconstruction of a sinusoidal input is calculated. The test frequency is pitched at an angle in a 2D parallel-beam acquisition geometry. Reconstructions are evaluated along the pitch of the object. The modulation transfer function (MTF) is calculated from the relative signal at various test frequencies. The MTF determines whether modulation is within detectable limits in oblique reconstructions. In the previous linear systems (LS) model [B. Zhao and W. Zhao, "Three-dimensional linear system analysis for breast tomosynthesis," Med. Phys. 35(12), 5219-5232 (2008)], the MTF was calculated only in reconstructed slices parallel to the detector. This paper generalizes the MTF calculation to reconstructed slices at all possible pitches. Unlike the previous LS model, this paper also analyzes the effect of object thickness on the MTF. A second test object that is considered is a rod whose long axis is pitched similar to the sinusoidal input. The rod is used to assess whether the length of an object can be correctly estimated in oblique reconstructions.

RESULTS

To simulate the conventional display of the reconstruction, slices are first created along a 0° pitch. This direction is perpendicular to the rays of the central projection. The authors show that the input frequency of a pitched sinusoidal object cannot be determined using these slices. By changing the pitch of the slice to match the object, it is shown that the input frequency is properly resolved. To prove that modulation is preserved in pitched slices, the MTF is also calculated. Modulation is within detectable limits over a broad range of pitches if the object is thin, but is detectable over a narrower range of pitches if the object is thick. Turning next to the second test object, it is shown that the length of a pitched rod can be correctly determined in oblique reconstructions. Concordant with the behavior of the MTF, the length estimate is accurate over a broad range of pitches if the object is thin, but is correct over a narrower range of pitches if the object is thick.

CONCLUSIONS

This work justifies the feasibility of oblique reconstructions in tomosynthesis. It is demonstrated that pitched test objects are most easily visualized with oblique reconstructions instead of conventional reconstructions. In order to achieve high image quality over a broad range of pitches, the object must be thin. By analyzing the effect of reconstruction pitch and object thickness on image quality, this paper generalizes the previous LS model for tomosynthesis.

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.

Proposing a New Velocity Profile for Continuous X-Ray Tube Motion in Digital Breast Tomosynthesis

In digital breast tomosynthesis (DBT), a 3D image of the breast is generated from x-ray projections at various angles. There are two mechanisms for acquiring projection images in DBT, step-and-shoot motion and continuous tube motion. The benefit of continuous tube motion is shorter scan time and hence less patient motion; the trade-off is focal spot blurring. To minimize focal spot blurring in a system with continuous tube motion, this study proposes a new velocity profile for the x-ray tube during the scan. Unlike existing systems for which the x-ray tube has constant angular velocity, we investigate a smoothly-varying tube velocity that approaches zero during each projection and is larger between projections. With this unique design, the filtered backprojection reconstruction of a sinusoidal test object was calculated, and modulation was determined at various frequencies. It is shown that the newly proposed tube velocity yields increased modulation in the reconstruction relative to a conventional system with continuous tube motion. The modulation in the re-designed system differs minimally from an analogous step-and-shoot system operated with the same scan time. This improvement in image quality was validated with reconstructions of microcalcifications in computer breast phantoms. It is known that continuous tube motion reduces the contrast of microcalcifications relative to step- and-shoot systems; we show that the newly proposed tube motion increases the contrast of microcalcifications compared to conventional continuous tube motion. In conclusion, this work proposes a strategy for optimizing the velocity of tube motion in DBT.