Conditioning data for calculation of the modulation transfer function

Seven general radiography x-ray detectors with pixel sizes ranging from 175 to 76μm: technical evaluation with the focus on orthopaedic imaging

Aim. Flat panel detectors with small pixel sizes general can potentially improve imaging performance in radiography applications requiring fine detail resolution. This study evaluated the imaging performance of seven detectors, covering a wide range of pixel sizes, in the frame of orthopaedic applications.Material and methods. Pixel sizes ranged from 175 (detector A175) to 76μm (detector G76). Modulation transfer function (MTF) and detective quantum efficiency (DQE) were measured using International Electrotechnical Commission (IEC) RQA3 beam quality. Threshold contrast (CT) and a detectability index (d') were measured at three air kerma/image levels. Rabbit shoulder images acquired at 60 kV, over five air kerma levels, were evaluated in a visual grading study for anatomical sharpness, image noise and overall diagnostic image quality by four radiologists. The detectors were compared to detector E124.Results. The 10% point of the MTF ranged from 3.21 to 4.80 mm-1, in going from detector A175to detector G76. DQE(0.5 mm-1) measured at 2.38μGy/image was 0.50 ± 0.05 for six detectors, but was higher for F100at 0.62. High frequency DQE was superior for the smaller pixel detectors, howeverCTfor 0.25 mm discs correlated best with DQE(0.5 mm-1). Correlation betweenCTand the detectability model was good (R2= 0.964).CTfor 0.25 mm diameter discs was significantly higher for D150and F100compared to E124. The visual grading data revealed higher image quality ratings for detectors D125and F100compared to E124. An increase in air kerma was associated with improved perceived sharpness and overall quality score, independent of detector. Detectors B150, D125, F100and G76, performed well in specific tests, however only F100consistently outperformed the reference detector.Conclusion. Pixel size alone was not a reliable predictor of small detail detectability or even perceived sharpness in a visual grading analysis study.

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.

Impact of an artificial intelligence deep-learning reconstruction algorithm for CT on image quality and potential dose reduction: A phantom study

Background

Recently, computed tomography (CT) manufacturers have developed deep‐learning‐based reconstruction algorithms to compensate for the limitations of iterative reconstruction (IR) algorithms, such as image smoothing and the spatial resolution's dependence on contrast and dose levels.

Purpose

To assess the impact of an artificial intelligence deep‐learning reconstruction (AI‐DLR) algorithm on image quality and dose reduction compared with a hybrid IR algorithm in chest CT for different clinical indications.

Methods

Acquisitions on the CT American College of Radiology (ACR) 464 and CT Torso CTU‐41 phantoms were performed at five dose levels (CTDI_vol: 9.5/7.5/6/2.5/0.4 mGy) used for chest CT conditions. Raw data were reconstructed using filtered backprojection, two levels of IR (iDose⁴ levels 4 (i4) and 7 (i7)), and five levels of AI‐DLR (Precise Image; Smoother, Smooth, Standard, Sharp, Sharper). Noise power spectrum (NPS), task‐based transfer function, and detectability index (d′) were computed: d′‐modeled detection of a soft tissue mediastinal nodule (low‐contrast soft tissue chest nodule within the mediastinum [LCN]), ground‐glass opacity (GGO), or high‐contrast pulmonary (HCP) lesion. The subjective image quality of chest anthropomorphic phantom images was independently evaluated by two radiologists. They assessed image noise, image smoothing, contrast between vessels and fat in the mediastinum for mediastinal images, visual border detection between bronchus and lung parenchyma for parenchymal images, and overall image quality using a commonly used four‐ or five‐point scale.

Results

From Standard to Smoother levels, on average, the noise magnitude decreased (for all dose levels: −66.3% ± 0.5% for mediastinal images and −63.1% ± 0.1% for parenchymal images), the average NPS spatial frequency decreased (for all dose levels: −35.3% ± 2.2% for mediastinal images and −13.3% ± 2.2% for parenchymal images), and the detectability (d′) of the three lesions increased. The opposite pattern was found from Standard to Sharper levels. From Smoother to Sharper levels, the spatial resolution increased for the low‐contrast polyethylene insert and the opposite for the high‐contrast air insert. Compared to the i4 used in clinical practice, d′ values were higher using Smoother (mean for all dose levels: 338.7% ± 29.4%), Smooth (103.4% ± 11.2%), and Standard (34.1% ± 6.6%) levels for the LCN on mediastinal images and Smoother (169.5% ± 53.2% for GGO and 136.9% ± 1.6% for HCP) and Smooth (36.4% ± 22.1% and 24.1% ± 0.9%, respectively) levels for parenchymal images. Radiologists considered the images satisfactory for clinical use at these levels, but adaptation to the dose level of the protocol is required.

Conclusion

With AI‐DLR, the smoothest levels reduced the noise and improved the detectability of chest lesions but increased the image smoothing. The opposite was found with the sharpest levels. The choice of level depends on the dose level and type of image: mediastinal or parenchymal.

Comparison of virtual monoenergetic imaging between a rapid kilovoltage switching dual-energy computed tomography with deep-learning and four dual-energy CTs with iterative reconstruction

Background

To assess the spectral performance of rapid kV switching dual-energy CT (KVSCT-Canon) equipped with a Deep-Learning spectral reconstruction algorithm on virtual-monoenergetic images at low-energy levels and to compare its performances with four other dual-energy CT (DECT) platforms equipped with iterative reconstruction algorithms.

Methods

Two CT phantoms were scanned on five DECT platforms: KVSCT-Canon, fast kV-switching CT (KVSCT-GE), split filter CT, dual-source CT (DSCT), and dual-layer CT (DLCT). The classical parameters of abdomen-pelvic examinations were used for all phantom acquisitions, and a CTDI_vol close to 10 mGy. For KVSCT-Canon, virtual-monoenergetic images were reconstructed with a clinical slice thickness of 0.5 and 1.5 mm to be close to other platforms. Noise power spectrum (NPS) and task-based transfer function (TTF) were evaluated from 40 to 80 keV of virtual-monoenergetic images. A detectability index (d') was computed to model the detection task of two contrast-enhanced lesions as function of keV.

Results

For KVSCT-Canon, the noise magnitude and average NPS spatial frequency (f_av) decreased from 40 to 70 keV and increased thereafter. Similar noise magnitude outcomes were found for KVSCT-GE but the opposite for f_av. For the other DECT platforms, the noise magnitude decreased as the keV increased. For split filter CT, DSCT and DLCT, the f_av values increased from 40 to 80 keV. For all DECT platforms, TTF at 50% (f₅₀) decreased as the keV increased, decreasing spatial resolution. For KVSCT-Canon, d' values peaked at 60 and 70 keV for both simulated lesions and from 50 to 70 keV for KVSCT-GE. d' decreased between 40 and 70 keV for DSCT, DLCT and split filter CT. For KVSCT-Canon, the increase in slice thickness decreases noise magnitude, f_av and f₅₀ and increases d' values. The highest d' values were found for DLCT at 40 and 50 keV and for KVSCT-Canon at 1.5 mm for other keV.

Conclusions

For KVSCT-Canon, the detectability of contrast-enhanced lesions was highest at 60 keV. The highest d' values were found for DLCT at 40 and 50 keV and for KVSCT-Canon at 1.5 mm for other keV.

Validation of a mammographic image quality modification algorithm using 3D-printed breast phantoms

Purpose: To validate a previously proposed algorithm that modifies a mammogram to appear as if it was acquired with different technique factors using realistic phantom-based mammograms. Approach: Two digital mammography systems (an indirect- and a direct-detector-based system) were used to acquire realistic mammographic images of five 3D-printed breast phantoms with the technique factors selected by the automatic exposure control and at various other conditions (denoted by the original images). Additional images under other simulated conditions were also acquired: higher or lower tube voltages, different anode/filter combinations, or lower tube current-time products (target images). The signal and noise in the original images were modified to simulate the target images (simulated images). The accuracy of the image modification algorithm was validated by comparing the target and simulated images using the local mean, local standard deviation (SD), local variance, and power spectra (PS) of the image signals. The absolute relative percent error between the target and simulated images for each parameter was calculated at each sub-region of interest (local parameters) and frequency (PS), and then averaged. Results: The local mean signal, local SD, local variance, and PS of the target and simulated images were very similar, with a relative percent error of 5.5%, 3.8%, 7.8%, and 4.4% (indirect system), respectively, and of 3.7%, 3.8%, 7.7%, and 7.5% (direct system), respectively. Conclusions: The algorithm is appropriate for simulating different technique factors. Therefore, it can be used in various studies, for instance to evaluate the impact of technique factors in cancer detection using clinical images.

Technical Note: A simple method for measuring the slice sensitivity profile of iteratively reconstructed CT images using a non-slanted edge plane

PURPOSE

A method for measuring the slice sensitivity profile (SSP) of computed tomography (CT) images reconstructed with iterative reconstruction (IR) algorithms was reported by the AAPM Task Group 233 (TG233). In this method, the phantom plane edge is slightly slanted with respect to the scan plane to obtain a composite oversampled edge-spread function (ESF). However, it is expected that a fine-sampled ESF can be obtained directly from images reconstructed with a small slice increment without slanting the edge plane. This study aimed to investigate the validity of using a non-slanted edge plane.

METHODS

In the proposed non-slanted edge method, the phantom was positioned so that the plane edge was perpendicular to the longitudinal z-axis, and images were reconstructed with a 1-mm slice thickness and 0.1-mm increment. The mean CT value was obtained in each slice and plotted as a function of slice position along the z-axis, thereby generating the ESF. The SSP was calculated from the ESF by differentiation. In the TG 233-recommended slanted edge method, the SSP was obtained by following the procedure described in the TG233 report. To validate the methodology, we first used filtered back projection (FBP) images to compare SSPs obtained using the non-slanted edge method, slanted edge method, and a standard method using a high-contrast thin object (coin). Next, for two types of IR algorithms, we compared the SSPs obtained using the non-slanted and slanted edge methods.

RESULTS

For the FBP images, the SSP measured using the non-slanted edge method agreed well with SSPs measured using the coin and slanted edge methods. For the IR images, the SSPs measured using the non-slanted and slanted edge methods showed good agreement.

CONCLUSIONS

The non-slanted edge method was demonstrated to be valid. The simplicity and practicality of the method allows routine and accurate determination of the SSP.

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.

Position paper: recommendations for a digital mammography quality assurance program V4.0

In 2001 the ACPSEM published a position paper on quality assurance in screen film mammography which was subsequently adopted as a basis for the quality assurance programs of both the Royal Australian and New Zealand College of Radiologists (RANZCR) and of BreastScreen Australia. Since then the clinical implementation of digital mammography has been realised and it has become evident that existing screen-film protocols were not appropriate to assure the required image quality needed for reliable diagnosis or to address the new dose implications resulting from digital technology. In addition, the advantages and responsibilities inherent in teleradiology are most critical in mammography and also need to be addressed. The current document is the result of a review of current overseas practice and local experience in these areas. At this time the technology of digital imaging is undergoing significant development and there is still a lack of full international consensus about some of the detailed quality control (QC) tests that should be included in quality assurance (QA) programs. This document describes the current status in digital mammography QA and recommends test procedures that may be suitable in the Australasian environment. For completeness, this document also includes a review of the QA programs required for the various types of digital biopsy units used in mammography. In the future, international harmonisation of digital quality assurance in mammography and changes in the technology may require a review of this document. Version 2.0 represented the first of these updates and key changes related to image quality evaluation, ghost image evaluation and interpretation of signal to noise ratio measurements. In Version 3.0 some significant changes, made in light of further experience gained in testing digital mammography equipment were introduced. In Version 4.0, further changes have been made, most notably digital breast tomosynthesis (DBT) testing and QC have been addressed. Some additional testing for conventional projection imaging has been added in order that sites may have the capability to undertake dose surveys to confirm compliance with diagnostic reference levels (DRLs) that may be established at the National or State level. A key recommendation is that dosimetry calculations are now to be undertaken using the methodology of Dance et al. Some minor changes to existing facility QC tests have been made to ensure the suggested procedures align with those most recently adopted by the Royal Australian and New Zealand College of Radiologists and BreastScreen Australia. Future updates of this document may be provided as deemed necessary in electronic format on the ACPSEM's website ( https://www.acpsem.org.au/whatacpsemdoes/standards-position-papers and see also http://www.ranzcr.edu.au/quality-a-safety/radiology/practice-quality-activities/mqap ).

Simulation of images of CDMAM phantom and the estimation of measurement uncertainties of threshold gold thickness

PURPOSE

To demonstrate a method of simulating mammography images of the CDMAM phantom and to investigate the coefficient of variation (CoV) in the threshold gold thickness (tT) measurements associated with use of the phantom.

METHODS

The noise and sharpness of Hologic Dimensions and GE Essential mammography systems were characterized to provide data for the simulation. The simulation method was validated by comparing the tT results of real and simulated images of the CDMAM phantom for three different doses and the two systems. The detection matrices produced from each of 64 images using CDCOM software were randomly resampled to create 512 sets of 8, 16 and 32 images to estimate the CoV of tT. Sets of simulated images for a range of doses were used to estimate the CoVs for a range of diameters and threshold thicknesses.

RESULTS

No significant differences were found for tT or the CoV between real and simulated CDMAM images. It was shown that resampling from 256 images was required for estimating the CoV. The CoV was around 4% using 16 images for most of the phantom but is over double that for details near the edge of the phantom.

CONCLUSIONS

We have demonstrated a method to simulate images of the CDMAM phantom for different systems at a range of doses. We provide data for calculating uncertainties in tT. Any future review of the European guidelines should take into consideration the calculated uncertainties for the 0.1mm detail.

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.

Characteristic image quality of a third generation dual-source MDCT scanner: Noise, resolution, and detectability

PURPOSE

The purpose of this work was to assess the inherent image quality characteristics of a new multidetector computed tomography system in terms of noise, resolution, and detectability index as a function of image acquisition and reconstruction for a range of clinically relevant settings.

METHODS

A multisized image quality phantom (37, 30, 23, 18.5, and 12 cm physical diameter) was imaged on a SOMATOM Force scanner (Siemens Medical Solutions) under variable dose, kVp, and tube current modulation settings. Images were reconstructed with filtered back projection (FBP) and with advanced modeled iterative reconstruction (ADMIRE) with iterative strengths of 3, 4, and 5. Image quality was assessed in terms of the noise power spectrum (NPS), task transfer function (TTF), and detectability index for a range of detection tasks (contrasts of approximately 45, 90, 300, -900, and 1000 HU, and 2-20 mm diameter) based on a non-prewhitening matched filter model observer with eye filter.

RESULTS

Image noise magnitude decreased with decreasing phantom size, increasing dose, and increasing ADMIRE strength, offering up to 64% noise reduction relative to FBP. Noise texture in terms of the NPS was similar between FBP and ADMIRE (<5% shift in peak frequency). The resolution, based on the TTF, improved with increased ADMIRE strength by an average of 15% in the TTF 50% frequency for ADMIRE-5. The detectability index increased with increasing dose and ADMIRE strength by an average of 55%, 90%, and 163% for ADMIRE 3, 4, and 5, respectively. Assessing the impact of mA modulation for a fixed average dose over the length of the phantom, detectability was up to 49% lower in smaller phantom sections and up to 26% higher in larger phantom sections for the modulated scan compared to a fixed tube current scan. Overall, the detectability exhibited less variability with phantom size for modulated scans compared to fixed tube current scans.

CONCLUSIONS

Image quality increased with increasing dose and decreasing phantom size. The CT system exhibited nonlinear noise and resolution properties, especially at very low-doses, large phantom sizes, and for low-contrast objects. Objective image quality metrics generally increased with increasing dose and ADMIRE strength, and with decreasing phantom size. The ADMIRE algorithm could offer comparable image quality at reduced doses or improved image quality at the same dose. The use of tube current modulation resulted in more consistent image quality with changing phantom size.

Estimation of the two-dimensional presampled modulation transfer function of digital radiography devices using one-dimensional test objects

PURPOSE

The modulation transfer function (MTF) of medical imaging devices is commonly reported in the form of orthogonal one-dimensional (1D) measurements made near the vertical and horizontal axes with a slit or edge test device. A more complete description is found by measuring the two-dimensional (2D) MTF. Some 2D test devices have been proposed, but there are some issues associated with their use: (1) they are not generally available; (2) they may require many images; (3) the results may have diminished accuracy; and (4) their implementation may be particularly cumbersome. This current work proposes the application of commonly available 1D test devices for practical and accurate estimation of the 2D presampled MTF of digital imaging systems.

METHODS

Theory was developed and applied to ensure adequate fine sampling of the system line spread function for 1D test devices at orientations other than approximately vertical and horizontal. Methods were also derived and tested for slit nonuniformity correction at arbitrary angle. Techniques were validated with experimental measurements at ten angles using an edge test object and three angles using a slit test device on an indirect-detection flat-panel system [GE Revolution XQ∕i (GE Healthcare, Waukesha, WI)]. The 2D MTF was estimated through a simple surface fit with interpolation based on Delaunay triangulation of the 1D edge-based MTF measurements. Validation by synthesis was also performed with simulated images from a hypothetical direct-detection flat-panel device.

RESULTS

The 2D MTF derived from physical measurements yielded an average relative precision error of 0.26% for frequencies below the cutoff (2.5 mm(-1)) and approximate circular symmetry at frequencies below 4 mm(-1). While slit analysis generally agreed with the results of edge analysis, the two showed subtle differences at frequencies above 4 mm(-1). Slit measurement near 45° revealed radial asymmetry in the MTF resulting from the square pixel aperture (0.2 mm × 0.2 mm), a characteristic which was not necessarily appreciated with the orthogonal 1D MTF measurements. In simulation experiments, both slit- and edge-based measurements resolved the radial asymmetries in the 2D MTF. The average absolute relative accuracy error in the 2D MTF between the DC and cutoff (2.5 mm(-1)) frequencies was 0.13% with average relative precision error of 0.11%. Other simulation results were similar to those derived from physical data.

CONCLUSIONS

Overall, the general availability, acceptance, accuracy, and ease of implementation of 1D test devices for MTF assessment make this a valuable technique for 2D MTF estimation.

Measured and calculated K-fluorescence effects on the MTF of an amorphous-selenium based CCD x-ray detector

PURPOSE

Theoretical reasoning suggests that direct conversion digital x-ray detectors based upon photoconductive amorphous-selenium (a-Se) could attain very high values of the MTF (modulation transfer function) at spatial frequencies well beyond 20 cycles mm(-1). One of the fundamental factors affecting resolution loss, particularly at x-ray energies just above the K-edge of selenium (12.66 keV), is the K-fluorescence reabsorption mechanism, wherein energy can be deposited in the detector at locations laterally displaced from the initial x-ray interaction site. This paper compares measured MTF changes above and below the Se K-edge of a CCD based a-Se x-ray detector with theoretical expectations.

METHODS

A prototype 25 μm sampling pitch (Nyquist frequency = 20 cycles mm(-1), 200 μm thick a-Se layer based x-ray detector, utilizing a specialized CCD readout device (200 × 400 area array), was used to make edge images with monochromatic x-rays above and below the K-edge of Se. A vacuum double crystal monochromator, exposed to polychromatic x-rays from a synchrotron, formed the monochromatic x-ray source. The monochromaticity of the x-rays was 99% or better. The presampling MTF was determined using the slanted edge method. The theory modeling the MTF performance of the detector includes the basic x-ray interaction physics in the a-Se layer as well as effects related to the operation of the CCD and charge trapping at a blocking layer present at the CCD/a-Se interface.

RESULTS

The MTF performance of the prototype a-Se CCD was reduced from the theoretical value prescribed by the basic Se x-ray interaction physics, principally by the presence of a blocking layer. Nevertheless, the K-fluorescence reduction in the MTF was observed, approximately as predicted by theory. For the CCD prototype detector, at five cycles mm(-1), there was a 14% reduction of the MTF, from a value of 0.7 below the K-edge of Se, to 0.6 just above the K-edge.

CONCLUSIONS

The MTF of an a-Se x-ray detector has been measured using monochromatic x-rays above and below the K-edge of selenium. The MTF is poorer above the K-edge by an amount consistent with theoretical expectations.

The effect of angular dose distribution on the detection of microcalcifications in digital breast tomosynthesis

PURPOSE

Substantial effort has been devoted to the clinical development of digital breast tomosynthesis (DBT). DBT is a three-dimensional (3D) x-ray imaging modality that reconstructs a number of thin image slices parallel to a stationary detector plane. Preliminary clinical studies have shown that the removal of overlapping breast tissue reduces image clutter and increases detectability of large, low contrast lesions. However, some studies, as well as anecdotal evidence, suggested decreased conspicuity of small, high contrast objects such as microcalcifications. Several investigators have proposed alternative imaging methods for improving microcalcification detection by delivering half of the total dose to the central view in addition to a separate DBT scan. Preliminary observer studies found possible improvement by either viewing the central projection alone or combining all views with a reconstruction algorithm.

METHODS

In this paper, we developed a generalized imaging theory based on a cascaded linear-system model for DBT to calculate the effect of variable angular dose distribution on the 3D modulation transfer function (MTF) and noise power spectrum (NPS). Using the ideal observer signal-to-noise ratio (SNR), d', as a figure-of-merit (FOM) for a signal embedded in a uniform background, we compared the detectability of objects with different sizes under different imaging conditions (e.g., angular dose distribution and reconstruction filters). Experimental investigation was conducted for three different angular dose schemes (ADS) using a Siemens Novation(TOMO) prototype unit.

RESULTS

Our results show excellent agreement between modeled and experimental measurements of 3D NPS with different angular dose distribution. The ideal observer detectability index for the detection of Gaussian objects with different angular dose distributions depends strongly on the applied reconstruction filter as well as the imaging task. For detection tasks of small calcifications with reconstruction filters used typically in a clinical setting, variable angular dose distribution with more dose delivered to the central views may lead to higher d' than a uniform angular dose distribution.

CONCLUSIONS

The conspicuity of the detection of small calcifications may be improved, under certain imaging conditions, by delivering higher dose toward the central views of a tomosynthesis scan, while also reducing the dose at peripheral angles to keep total administered radiation dose equivalent. The degree of improvement depends on the choice of reconstruction filters as well as the imaging task. The improvement is more substantial for high-frequency imaging tasks and when an aggressive slice-thickness (ST) filter is applied to reduced the high-frequency noise at peripheral angles.

Accurate MTF measurement in digital radiography using noise response

PURPOSE

The authors describe a new technique to determine the system presampled modulation transfer function (MTF) in digital radiography using only the detector noise response.

METHODS

A cascaded-linear systems analysis was used to develop an exact relationship between the two-dimensional noise power spectrum (NPS) and the presampled MTF for a generalized detector system. This relationship was then utilized to determine the two-dimensional presampled MTF. For simplicity, aliasing of the correlated noise component of the NPS was assumed to be negligible. Accuracy of this method was investigated using simulated images from a simple detector model in which the "true" MTF was known exactly. Measurements were also performed on three detector technologies (an x-ray image intensifier, an indirect flat panel detector, and a solid state x-ray image intensifier), and the results were compared using the standard edge-response method. Flat-field and edge images were acquired and analyzed according to guidelines set forth by the International Electrotechnical Commission, using the RQA 5 spectrum.

RESULTS

The presampled MTF determined using the noise-response method for the simulated detector system was in close agreement with the true MTF with an averaged percent difference of 0.3% and a maximum difference of 1.1% observed at the Nyquist frequency (fN). The edge-response method of the simulated detector system also showed very good agreement at lower spatial frequencies (less than 0.5 fN) with an averaged percent difference of 1.6% but showed significant discrepancies at higher spatial frequencies (greater than 0.5 fN) with an averaged percent difference of 17%. Discrepancies were in part a result of noise in the edge image and phasing errors. For all three detector systems, the MTFs obtained using the two methods were found to be in good agreement at spatial frequencies less than 0.5 fN with an averaged percent difference of 3.4%. Above 0.5 fN, differences increased to an average of 20%. Deviations of the experimental results largely followed the trend seen in the simulation results, suggesting that differences between the two methods could be explained as resulting from the inherent inaccuracies of the edge-response measurement technique used in this study. Aliasing of the correlated noise component was shown to have a minimal effect on the measured MTF for the three detectors studied. Systems with significant aliasing of the correlated noise component (e.g., a-Se based detectors) would likely require a more sophisticated fitting scheme to provide accurate results.

CONCLUSIONS

Results indicate that the noise-response method, a simple technique, can be used to accurately measure the MTF of digital x-ray detectors, while alleviating the problems and inaccuracies associated with use of precision test objects, such as a slit or an edge.

Experimental validation of a three-dimensional linear system model for breast tomosynthesis

A three-dimensional (3D) linear model for digital breast tomosynthesis (DBT) was developed to investigate the effects of different imaging system parameters on the reconstructed image quality. In the present work, experimental validation of the model was performed on a prototype DBT system equipped with an amorphous selenium (a-Se) digital mammography detector and filtered back-projection (FBP) reconstruction methods. The detector can be operated in either full resolution with 85 microm pixel size or 2 x 1 pixel binning mode to reduce acquisition time. Twenty-five projection images were acquired with a nominal angular range of +/- 20 degrees. The images were reconstructed using a slice thickness of 1 mm with 0.085 x 0.085 mm in-plane pixel dimension. The imaging performance was characterized by spatial frequency-dependent parameters including a 3D noise power spectrum (NPS) and in-plane modulation transfer function (MTF). Scatter-free uniform x-ray images were acquired at four different exposure levels for noise analysis. An aluminum (Al) edge phantom with 0.2 mm thickness was imaged to measure the in-plane presampling MTF. The measured in-plane MTF and 3D NPS were both in good agreement with the model. The dependence of DBT image quality on reconstruction filters was investigated. It was found that the slice thickness (ST) filter, a Hanning window to limit the high-frequency components in the slice thickness direction, reduces noise aliasing and improves 3D DQE. An ACR phantom was imaged to investigate the effects of angular range and detector operational modes on reconstructed image quality. It was found that increasing the angular range improves the MTF at low frequencies, resulting in better detection of large-area, low-contrast mass lesions in the phantom. There is a trade-off between noise and resolution for pixel binning and full resolution modes, and the choice of detector mode will depend on radiation dose and the targeted lesion.

The optimal optical readout for the x-ray light valve--document scanners

The x-ray light valve (XLV) is a novel, potentially low-cost, x-ray detector that converts an x-ray exposure into an optical image stored in a liquid crystal cell. This optical image is then transferred from the liquid crystal cell to a computer through an optical-to-digital imaging readout system. Previously, CCD-based cameras were used for the optical readout, but recently it was proposed that an inexpensive optical scanner, such as an office document scanner, is a better match to the optical properties of the XLV. A methodology for characterizing a document scanner's ability to produce medical quality images from the XLV is outlined and tested on a particular scanner (Canon LiDE 30). This scanner was shown to have key characteristics of a medical device-a linear response, dynamic range sufficient for chest radiography (although not mammography) in a single pass, and an MTF and NPS that exceed the requirements for all medical applications of the scanner. This combination of criteria shows that a document scanner can be used as a digitization method for the XLV.

Singular value description of a digital radiographic detector: theory and measurements

The H operator represents the deterministic performance of any imaging system. For a linear, digital imaging system, this system operator can be written in terms of a matrix, H, that describes the deterministic response of the system to a set of point objects. A singular value decomposition of this matrix results in a set of orthogonal functions (singular vectors) that form the system basis. A linear combination of these vectors completely describes the transfer of objects through the linear system, where the respective singular values associated with each singular vector describe the magnitude with which that contribution to the object is transferred through the system. This paper is focused on the measurement, analysis, and interpretation of the H matrix for digital x-ray detectors. A key ingredient in the measurement of the H matrix is the detector response to a single x ray (or infinitestimal x-ray beam). The authors have developed a method to estimate the 2D detector shift-variant, asymmetric ray response function (RRF) from multiple measured line response functions (LRFs) using a modified edge technique. The RRF measurements cover a range of x-ray incident angles from 0 degree (equivalent location at the detector center) to 30 degrees (equivalent location at the detector edge) for a standard radiographic or cone-beam CT geometric setup. To demonstrate the method, three beam qualities were tested using the inherent, Lu/Er, and Yb beam filtration. The authors show that measures using the LRF, derived from an edge measurement, underestimate the system's performance when compared with the H matrix derived using the RRF. Furthermore, the authors show that edge measurements must be performed at multiple directions in order to capture rotational asymmetries of the RRF. The authors interpret the results of the H matrix SVD and provide correlations with the familiar MTF methodology. Discussion is made about the benefits of the H matrix technique with regards to signal detection theory, and the characterization of shift-variant imaging systems.

The effect of breast compression on mass conspicuity in digital mammography

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

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.

Interim recommendations for a digital mammography quality assurance program

In 2001 the ACPSEM published a position paper on quality assurance in screen film mammography which was subsequently adopted as a basis for the quality assurance programs of both the Royal Australian and New Zealand College of Radiologists (RANZCR) and of BreastScreen Australia. Since then the clinical implementation of digital mammography has been realised and it has become evident that existing screen-film protocols were not appropriate to assure the required image quality needed for reliable diagnosis or to address the new dose implications resulting from digital technology. In addition, the advantages and responsibilities inherent in teleradiology are most critical in mammography and also need to be addressed. The current document is the result of a review of current overseas practice and local experience in these areas. At this time the technology of digital imaging is undergoing significant development and there is still a lack of full international consensus about some of the detailed Quality Control tests that should be included in quality assurance (QA) programs. This document describes the current status in digital mammography QA and recommends test procedures that may be suitable in the Australasian environment. For completeness, this document also includes a review of the QA programs required for the various types of digital biopsy units used in mammography. In the future, international harmonisation of digital quality assurance in mammography and changes in the technology may require a review of this document. Accordingly, updates of this document will be provided as deemed necessary in electronic format on the ACPSEM's website (see http://www.acpsem.org.au/au/subgroup/radiology/RadiologySG_index.html).

The dependence of the modulation transfer function on the blocking layer thickness in amorphous selenium x-ray detectors

Blocking layers are used to reduce leakage current in amorphous selenium detectors. The effect of the thickness of the blocking layer on the presampling modulation transfer function (MTF) and on dark current was experimentally determined in prototype single-line CCD-based amorphous selenium (a-Se) x-ray detectors. The sampling pitch of the detectors evaluated was 25 microm and the blocking layer thicknesses varied from 1 to 51 microm. The blocking layers resided on the signal collection electrodes which, in this configuration, were used to collect electrons. The combined thickness of the blocking layer and a-Se bulk in each detector was approximately 200 microm. As expected, the dark current increased monotonically as the thickness of the blocking layer was decreased. It was found that if the blocking layer thickness was small compared to the sampling pitch, it caused a negligible reduction in MTF. However, the MTF was observed to decrease dramatically at spatial frequencies near the Nyquist frequency as the blocking layer thickness approached or exceeded the electrode sampling pitch. This observed reduction in MTF is shown to be consistent with predictions of an electrostatic model wherein the image charge from the a-Se is trapped at a characteristic depth within the blocking layer, generally near the interface between the blocking layer and the a-Se bulk.

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

Oblique incidence of x rays on an imaging detector causes blurring that reduces spatial resolution. For simple projection imaging this effect is small and often ignored. However, for breast tomosynthesis, the incidence angle can be larger (>20 degrees), leading to increased blur for some of the projections. The modulation transfer function (MTF) is measured for a typical phosphor-coupled flat-panel detector versus angular incidence of the x-ray beam for two x-ray spectra: 26 kV Mo/Mo and 40 kV Rh/Al. At an incidence angle of 40 degrees the MTF at 5 mm(-1) falls by 35% and 40% for each spectrum, respectively (and 65%/80% at 8 mm(-1)). Increasing the detector absorber thickness to achieve improved quantum efficiency will cause the blurring effect due to beam obliquity to become greater. The impact of this blur is likely to cause misregistration and increased relative noise in tomosynthesis reconstructed images.

Retrospective analysis of a detector fault for a full field digital mammography system

This paper describes objective and subjective image quality measurements acquired as part of a routine quality assurance (QA) programme for an amorphous selenium (a-Se) full field digital mammography (FFDM) system between August-04 and February-05. During this period, the FFDM detector developed a fault and was replaced. A retrospective analysis of objective image quality parameters (modulation transfer function (MTF), normalized noise power spectrum (NNPS) and detective quantum efficiency (DQE)) is presented to try and gain a deeper understanding of the detector problem that occurred. These measurements are discussed in conjunction with routine contrast-detail (c-d) results acquired with the CDMAM (Artinis, The Netherlands) test object. There was significant reduction in MTF over this period of time indicating an increase in blurring occurring within the a-Se converter layer. This blurring was not isotropic, being greater in the data line direction (left to right across the detector) than in the gate line direction (chest wall to nipple). The initial value of the 50% MTF point was 6 mm(-1); for the faulty detector the 50% MTF points occurred at 3.4 mm(-1) and 1.0 mm(-1) in the gate line and data line directions, respectively. Prior to NNPS estimation, variance images were formed of the detector flat field images. Spatial distribution of variance was not uniform, suggesting that the physical blurring process was not constant across the detector. This change in variance with image position implied that the stationarity of the noise statistics within the image was limited and that care would be needed when performing objective measurements. The NNPS measurements confirmed the results found for the MTF, with a strong reduction in NNPS as a function of spatial frequency. This reduction was far more severe in the data line direction. A somewhat tentative DQE estimate was made; in the gate line direction there was little change in DQE up to 2.5 mm(-1) but at the Nyquist frequency the DQE had fallen to approximately 35% of the original value. There was severe attenuation of DQE in the data line direction, the DQE falling to less than 0.01 above approximately 3.0 mm(-1). C-d results showed an increase in threshold contrast of approximately 25% for details less than 0.2 mm in diameter, while no reduction in c-d performance was found at the largest detail diameters (1.0 mm and above). Despite the detector fault, the c-d curve was found to pass the European protocol acceptable c-d curve.

A comparison between objective and subjective image quality measurements for a full field digital mammography system

This paper presents pre-sampling modulation transfer function (MTF), normalized noise power spectrum (NNPS) and detective quantum efficiency (DQE) results for an amorphous selenium (a-Se) full field digital mammography system. MTF was calculated from the image of an angled 0.5 mm thick Cu edge, acquired without additional beam filtration. NNPS data were acquired at detector air-kerma levels ranging from 9.1 microGy to 331 microGy, using a standard mammography x-ray spectrum of 28 kV, Mo/Mo target/filter combination and 4 cm of PMMA additional filtration. Prior to NNPS estimation, the image statistics were assessed using a variance image. This method was able to easily identify a detector artefact and should prove useful in routine quality assurance (QA) measurements. Detector DQE, calculated from the NNPS and MTF data, dropped to 0.3 for low detector air-kerma settings but reached an approximately constant value of 0.6 above 50 microGy at the detector. Subjective image quality data were also obtained at these detector air-kerma settings using the CDMAM contrast-detail (c-d) test object. The c-d data reflected the trend seen in DQE, with threshold contrast increasing at low detector air-kerma values. The c-d data were then compared against predictions made using two established models, the Rose model and a standard signal detection theory model. Using DQE(0), the Rose model gave results within approximately 15% on average for all the detector air-kerma values studied and for detail diameters down to 0.2 mm. Similar agreement was also found between the measured c-d data and the signal detection theory results, which were calculated using an ideal human visual response function and a system magnification of unity. The use of full spatial frequency DQE improved the agreement between the calculated and observer results for detail sizes below 0.13 mm.

Comparison of full-field digital mammography to screen-film mammography with respect to contrast and spatial resolution in tissue equivalent breast phantoms

To determine if the improved contrast resolution of full-field digital mammography (FFDM) with reduced spatial resolution allows for superior or equal phantom object detection compared with screen-film mammography (SFM). Tissue equivalent breast phantoms simulating an adipose to glandular ratio of 50/50,30/70, and 20/80 were imaged according to each manufacturers' recommendation with four full-field digital mammography units (Fuji, Sectra, Fischer, and General Electric) and a screen-film mammography unit (MammoMatII 2000, Siemens, Munich, Germany). A total of 20 images were obtained in both hard- and soft-copy formats. For the purpose of soft-copy display, the screen-film hard-copy images were digitized with a 50 microm micron scanner. Six radiologists, experts in breast imaging, and three physicists, experts in scoring mammography phantoms, participated in a reader study where each reader scored each phantom for visibility of line-pairs and for 24 objects (fibers, clusters of specks, and masses). The data were recorded, entered into a database, and analyzed by a mixed-effect model. The limiting spatial resolution in line-pairs per millimeter visible with the digital units was less, regardless of display modality used, than that provided by the screen-film unit. The difference was statistically significant for the General Electric (p < 0.01) and Fuji digital mammography units (p = 0.03). With respect to the number of visible objects, a statistically significant higher number could be detected with the screen-film unit as compared to the Fischer (p < 0.01) and Sectra (p < 0.01) digital mammography units, but there was no significant difference between the other digital units and screen film. Overall, there was significantly better performance on the 50/50 phantom than with the 30/70 and 20/80 phantoms (p = 0.01, p < 0.01) for object visibility. For the digital mammography units, soft-copy display performed better than hard-copy display for the Fischer and Sectra images, but worse for Fuji and General Electric. In addition, soft-copy display of digitized screen-film images was significantly better than hard-copy display (p =0.02) of the original screen films for object visibility, but worse for spatial resolution. The higher contrast resolution of the FFDM units tested did not result in improved detection of line-pair resolution or objects in the phantoms tested versus screen-film mammography. The phantom performance of a digital mammography unit seems to be influenced by the type of detection task (line-pair resolution versus object visibility), the display modality (soft-copy versus hard-copy) chosen to score the phantoms, and the parenchymal pattern composition of the phantom.

Validation of MTF measurement for digital mammography quality control

The modulation transfer function (MTF) describes the spatial resolution properties of imaging systems. In this work, the accuracy of our implementation of the edge method for calculating the presampled MTF was examined. Synthetic edge images with known MTF were used as gold standards for determining the robustness of the edge method. These images simulated realistic data from clinical digital mammography systems, and contained intrinsic system factors that could affect the MTF accuracy, such as noise, scatter, and flat-field nonuniformities. Our algorithm is not influenced by detector dose variations for MTF accuracy up to 1/2 the sampling frequency. We investigated several methods for noise reduction, including truncating the supersampled line spread function (LSF), windowing the LSF, applying a local exponential fit to the LSF, and applying a monotonic constraint to the supersampled edge spread function. Only the monotonic constraint did not introduce a systematic error; the other methods could result in MTF underestimation. Overall, our edge method consistently computed MTFs which were in good agreement with the true MTF. The edge method was then applied to images from a commercial storage-phosphor based digital mammography system. The calculated MTF was affected by the size (sides of 2.5, 5, or 10 cm) and the composition (lead or tungsten) of the edge device. However, the effects on the MTF were observed only with regard to the low frequency drop (LFD). Scatter nonuniformity was dependent on edge size, and could lead to slight underestimation of LFD. Nevertheless, this negative effect could be minimized by using an edge of 5 cm or larger. An edge composed of lead is susceptible to L-fluorescence, which causes overestimation of the LFD. The results of this work are intended to underline the need for clear guidelines if the MTF is to be given a more crucial role in acceptance tests and routine assessment of digital mammography systems: the MTF algorithm and edge object test tool need to be publicly validated.

Image quality measurements and metrics in full field digital mammography: an overview

This paper gives an overview of test procedures developed to assess the performance of full field digital mammography systems. We make a distinction between tests of the individual components of the imaging chain and global system tests. Most tests are not yet fully standardised. Where possible, we illustrate the test methodologies on a selenium flat-panel system.

Evaluation of the imaging properties of an amorphous selenium-based flat panel detector for digital fluoroscopy

The imaging performance of an amorphous selenium (a-Se) flat-panel detector for digital fluoroscopy was experimentally evaluated using the spatial frequency dependent modulation transfer function (MTF), noise power spectrum (NPS), and detective quantum efficiency (DQE). These parameters were investigated at beam qualities and exposures within the range typical of gastrointestinal fluoroscopic imaging (approximately 0.1 - 10 microR, 75 kV). The investigation does not take into consideration the detector cover, which in clinical use will lower the DQE measured here by its percent attenuation. The MTF was found to be less than the expected aperture response and the NPS was not white which together indicate presampling blurring. The cause of this blurring was attributed to charge trapping at the interface between two different layers of the a-Se. The effect on the DQE was also consistent with presampling blur, which reduces the aliasing in the NPS and thereby reduces the spatial frequency dependence of the DQE. (The DQE was independent of spatial frequency from 0.12 to 0.73 mm(-1) due to antialiasing of the NPS.) Moreover, the first zero of the measured MTF and the aperture response appeared at the same spatial frequency (6.66 mm(-1) for a pixel of 150 microm). Hence, the geometric fill factor (77%) was increased to an effective fill factor of 99 +/- 1%. A large scale ( approximately 32 pixels) correlation in the noise due to the configuration of the readout electronics caused increased noise power in the gate line NPS at low spatial frequency (< 0.1 mm(-1)). The DQE (f = 0) was exposure independent over a large range of exposures but became exposure dependent at low exposures due to the electronic noise.