In-plane Modulation Transfer Function Research Articles

Task-based detectability in anatomical background in digital mammography, digital breast tomosynthesis and synthetic mammography

Objective. Determining the detectability of targets for the different imaging modalities in mammography in the presence of anatomical background noise is challenging. This work proposes a method to compare the image quality and detectability of targets in digital mammography (DM), digital breast tomosynthesis (DBT) and synthetic mammography. Approach. The low-frequency structured noise produced by a water phantom with acrylic spheres was used to simulate anatomical background noise for the different types of images. A method was developed to apply the non-prewhitening observer model with eye filter (NPWE) in these conditions. A homogeneous poly(methyl) methacrylate phantom with a 0.2 mm thick aluminium disc was used to calculate 2D in-plane modulation transfer function (MTF), noise power spectrum (NPS), noise equivalent quanta, and system detective quantum efficiency for 30, 50 and 70 mm thicknesses. The in-depth MTFs of DBT volumes were determined using a thin tungsten wire. The MTF, system NPS and anatomical NPS were used in the NPWE model to calculate the threshold gold thickness of the gold discs contained in the CDMAM phantom, which was taken as reference. Main results. The correspondence between the NPWE model and the CDMAM phantom (linear Pearson correlation 0.980) yielded a threshold detectability index that was used to determine the threshold diameter of spherical microcalcifications and masses. DBT imaging improved the detection of masses, which depended mostly on the reduction of anatomical background noise. Conversely, DM images yielded the best detection of microcalcification s. Significance. The method presented in this study was able to quantify image quality and object detectability for the different imaging modalities and levels of anatomical background noise.

Effect of optical blurring of X-ray source on breast tomosynthesis image quality: Modulation transfer function, anatomical noise power spectrum, and signal detectability perspectives

We investigated the effect of the optical blurring of X-ray source on digital breast tomosynthesis (DBT) image quality using well-designed DBT simulator and table-top experimental systems. To measure the in-plane modulation transfer function (MTF), we used simulated sphere phantom and Teflon sphere phantom and generated their projection data using two acquisition modes (i.e., step-and-shoot mode and continuous mode). After reconstruction, we measured the in-plane MTF using reconstructed sphere phantom images. In addition, we measured the anatomical noise power spectrum (aNPS) and signal detectability. We constructed simulated breast phantoms with a 50% volume glandular fraction (VGF) of breast anatomy using the power law spectrum and inserted spherical objects with 1 mm, 2 mm, and 5 mm diameters as breast masses. Projection data were acquired using two acquisition modes, and in-plane breast images were reconstructed using the Feldkamp-Davis-Kress (FDK) algorithm. For the experimental study, we used BR3D breast phantom with 50% VGF and obtained projection data using a table-top experimental system. To compare the detection performance of the two acquisition modes, we calculated the signal detectability using the channelized Hotelling observer (CHO) with Laguerre-Gauss (LG) channels. Our results show that spatial resolution of in-plane image in continuous mode was degraded due to the optical blurring of X-ray source. This blurring effect was reflected in aNPS, resulting in large β values. From a signal detectability perspective, the signal detectability in step-and-shoot mode is higher than that in continuous mode for small spherical signals but not large spherical signals. Although the step-and-shoot mode has disadvantage in terms of scan time compared to the continuous mode, scanning in step-and-shoot mode is better for detecting small signals, indicating that there is a tradeoff between scan time and image quality.

In-plane image quality and NPWE detectability index in digital breast tomosynthesis

A rigorous 2D analysis of signal and noise transfer applied to reconstructed planes in digital breast tomosynthesis (DBT) is necessary for system characterization and optimization. This work proposes a method for assessing technical image quality and system detective quantum efficiency (DQEsys) for reconstructed planes in DBT. Measurements of 2D in-plane modulation transfer function (MTF) and noise power spectrum (NPS) were made on five DBT systems using different acquisition parameters, reconstruction algorithms and plane spacing. This work develops the noise equivalent quanta (NEQ), DQEsys and detectability index (d’) calculated using a non-prewhitening model observer with eye filter (NPWE) for reconstructed DBT planes. The images required for this implementation were acquired using a homogeneous test object of thickness 40 mm poly(methyl) methacrylate plus 0.5 mm Al; 2D MTF was calculated from an Al disc of thickness 0.2 mm and diameter 50 mm positioned within the phantom. The radiant contrast of the MTF disc and the air kerma at the system input were used as normalization factors. The NPWE detectability index was then compared to the in-plane contrast-detail (c-d) threshold measured using the CDMAM phantom.The MTF and NPS measured on the different systems showed a strong anisotropy, consistent with the cascaded models developed in the literature for DBT. Detectability indices calculated from the measured MTF and NPS successfully predicted changes in c-d detectability for details between 0.1 mm and 2.0 mm, for DBT plane spacings between 0.5 mm and 10 mm, and for air kerma values at the system input between 157 µGy and 1170 μGy. The linear Pearson correlation between the detectability index and threshold gold thickness of the CDMAM phantom was −0.996. The method implements a parametric means of assessing the technical image quality of reconstructed DBT planes, providing valuable information for optimization of DBT systems.

Development of a scanner-specific simulation framework for photon-counting computed tomography

The aim of this study was to develop and validate a simulation platform that generates photon-counting CT images of voxelized phantoms with detailed modeling of manufacturer-specific components including the geometry and physics of the x-ray source, source filtrations, anti-scatter grids, and photon-counting detectors. The simulator generates projection images accounting for both primary and scattered photons using a computational phantom, scanner configuration, and imaging settings. Beam hardening artifacts are corrected using a spectrum and threshold dependent water correction algorithm. Physical and computational versions of a clinical phantom (ACR) were used for validation purposes. The physical phantom was imaged using a research prototype photon-counting CT (Siemens Healthcare) with standard (macro) mode, at four dose levels and with two energy thresholds. The computational phantom was imaged with the developed simulator with the same parameters and settings used in the actual acquisition. Images from both the real and simulated acquisitions were reconstructed using a reconstruction software (FreeCT). Primary image quality metrics such as noise magnitude, noise ratio, noise correlation coefficients, noise power spectrum, CT number, in-plane modulation transfer function, and slice sensitivity profiles were extracted from both real and simulated data and compared. The simulator was further evaluated for imaging contrast materials (bismuth, iodine, and gadolinium) at three concentration levels and six energy thresholds. Qualitatively, the simulated images showed similar appearance to the real ones. Quantitatively, the average relative error in image quality measurements were all less than 4% across all the measurements. The developed simulator will enable systematic optimization and evaluation of the emerging photon-counting computed tomography technology.

Inverse filtering approach to measure directional in-plane modulation transfer function using a sphere phantom for a digital tomosynthesis system.

We propose a method to measure the directional in-plane modulation transfer function (MTF) of a digital tomosynthesis system using a sphere phantom. To assess the spatial resolution of an in-plane image of the tomosynthesis system, projection data of a sphere phantom were generated within a limited data acquisition range of 40°, and reconstructed by the FDK algorithm. To measure the in-plane MTF, we divided the Fourier transform of the reconstructed sphere phantom by that of the ideal sphere phantom, and then performed plane integral along the fz-direction. When dividing, small values in the denominator can introduce estimation errors, and these errors were reduced by the proposed method. To evaluate the performance of the proposed method, the in-plane MTF estimated by simulation and experimental data was compared to the ideal in-plane MTF generated by computer simulations using a point object. For quantitative evaluation, we measured frequency values at half-maximum and full-maximum of the directional in-plane MTF along the three different directions (i.e., f0° -, f30° -, and f60° -directions) and compared them with those of the ideal in-plane MTF. Although the sphere phantom has been regarded as an inappropriate object due to the anisotropic characteristics of tomosynthesis image, our results show that the proposed method has a reliable estimation performance, demonstrating the sphere phantom is still suitable for measuring the directional in-plane MTF for a digital tomosynthesis system.

Three-dimensional cascaded system analysis of a 50 µm pixel pitch wafer-scale CMOS active pixel sensor x-ray detector for digital breast tomosynthesis

High-resolution, low-noise x-ray detectors based on the complementary metal-oxide-semiconductor (CMOS) active pixel sensor (APS) technology have been developed and proposed for digital breast tomosynthesis (DBT). In this study, we evaluated the three-dimensional (3D) imaging performance of a 50 µm pixel pitch CMOS APS x-ray detector named DynAMITe (Dynamic Range Adjustable for Medical Imaging Technology). The two-dimensional (2D) angle-dependent modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE) were experimentally characterized and modeled using the cascaded system analysis at oblique incident angles up to 30°. The cascaded system model was extended to the 3D spatial frequency space in combination with the filtered back-projection (FBP) reconstruction method to calculate the 3D and in-plane MTF, NNPS and DQE parameters. The results demonstrate that the beam obliquity blurs the 2D MTF and DQE in the high spatial frequency range. However, this effect can be eliminated after FBP image reconstruction. In addition, impacts of the image acquisition geometry and detector parameters were evaluated using the 3D cascaded system analysis for DBT. The result shows that a wider projection angle range (e.g. ±30°) improves the low spatial frequency (below 5 mm−1) performance of the CMOS APS detector. In addition, to maintain a high spatial resolution for DBT, a focal spot size of smaller than 0.3 mm should be used. Theoretical analysis suggests that a pixelated scintillator in combination with the 50 µm pixel pitch CMOS APS detector could further improve the 3D image resolution. Finally, the 3D imaging performance of the CMOS APS and an indirect amorphous silicon (a-Si:H) thin-film transistor (TFT) passive pixel sensor (PPS) detector was simulated and compared.

Further results on the evaluation of the performance of a digital breast tomosynthesis system in the clinical environment

Introduction Digital breast tomosynthesis (DBT) is spreading throughout medical centres worldwide. However, there is a wide variability of technical characteristics among DBT systems that need to be assessed for quality control (QC) implementation. Purpose After evaluating the technical performance of three commercial DBT systems in a previous work, we pursue to research a fourth system with notorious differences with the previous ones. Materials and methods We evaluate a General Electric SenoClaire DBT system with angular range (AR) of 25°, 9 projections, step-and-shoot X-ray tube motion and iterative reconstruction. Tests for assessing the image quality of projections and reconstructed planes were: contrast-to-noise ratio (CNR); spatial resolution, via the in-plane modulation transfer function (MTF) and the in-depth point spread function (Z-PSF), and artefact spread function (ASF). All was implemented in custom-designed phantoms according to literature. Results In the reconstructed slices and slabs, CNR has a very weak dependency on PMMA thickness, as opposed to projections, possibly due to the iterative algorithm. The tube travel MTF is 30% larger respect to the measured average from continuous tube movement systems, due to the step-and-shoot setup, while the chest-wall nipple resolution remains similar. The Z-PSF is narrower than systems with larger AR due to the iterative reconstruction but the ASF is similar to systems with a smaller AR, because of a greater angular spacing between projections. Conclusion We observe in the clinic some theoretical results when comparing with previously studied DBT systems, and the results could be useful as a guide for future QC procedures. Disclosure None to disclose.

Investigation of spatial resolution characteristics of an in vivo microcomputed tomography system

The spatial resolution characteristics of an in vivo microcomputed tomography (CT) system was investigated in the in-plane (x–y), cross plane (z) and projection imaging modes. The microCT system utilized in this study employs a flat panel detector with a 127µm pixel pitch, a microfocus x-ray tube with a focal spot size ranging from 5–30µm, and accommodates three geometric magnifications (M) of 1.72, 2.54 and 5.10. The in-plane modulation transfer function (MTF) curves were measured as a function of the number of projections, geometric magnification (M), detector binning and reconstruction magnification (MRecon). The in plane cutoff frequency (10% MTF) ranged from 2.31lp/mm (M=1.72, 2×2 binning) to 12.56lp/mm (M=5.10, 1×1 binning) and a bar pattern phantom validated those measurements. A slight degradation in the spatial resolution was observed when comparing the image reconstruction with 511 and 918 projections, whose effect was visible at the lower frequencies. Small value of MRecon has little or no impact on the in-plane spatial resolution owning to a stable system. Large value of MRecon has implications on the spatial resolution and it was evident when comparing the bar pattern images reconstructed with MRecon=1.25 and 2.5. The cross plane MTF curves showed that the spatial resolution increased as the slice thickness decreased. The cutoff frequencies in the projection imaging mode yielded slightly higher values as compared to the in-plane and cross plane modes at all the geometric magnifications (M). At M=5.10, the cutoff resolution of the projection and cross plane on an ultra-high contrast resolution bar chip phantom were 14.9lp/mm and 13–13.5lp/mm. Due to the finite focal spot size of the x-ray tube, the detector blur and the reconstruction kernel functions, the system׳s spatial resolution does not reach the limiting spatial resolution as defined by the Nyquist׳s detector criteria with an ideal point source. The geometric magnification employed in the microCTs provide a tradeoff between field of view and spatial resolution for a wide range of applications.

The effect of averaging adjacent planes for artifact reduction in matrix inversion tomosynthesis

Matrix inversion tomosynthesis (MITS) uses linear systems theory and knowledge of the imaging geometry to remove tomographic blur that is present in conventional backprojection tomosynthesis reconstructions, leaving in-plane detail rendered clearly. The use of partial-pixel interpolation during the backprojection process introduces imprecision in the MITS modeling of tomographic blur, and creates low-contrast artifacts in some MITS planes. This paper examines the use of MITS slabs, created by averaging several adjacent MITS planes, as a method for suppressing partial-pixel artifacts. Human chest tomosynthesis projection data, acquired as part of an IRB-approved pilot study, were used to generate MITS planes, three-plane MITS slabs (MITSa3), five-plane MITS slabs (MITSa5), and seven-plane MITS slabs (MITSa7). These were qualitatively examined for partial-pixel artifacts and the visibility of normal and abnormal anatomy. Additionally, small (5 mm) subtle pulmonary nodules were simulated and digitally superimposed upon human chest tomosynthesis projection images, and their visibility was qualitatively assessed in the different reconstruction techniques. Simulated images of a thin wire were used to generate modulation transfer function (MTF) and slice-sensitivity profile curves for the different MITS and MITS slab techniques, and these were examined for indications of partial-pixel artifacts and frequency response uniformity. Finally, mean-subtracted, exposure-normalized noise power spectra (ENNPS) estimates were computed and compared for MITS and MITS slab reconstructions, generated from 10 sets of tomosynthesis projection data of an acrylic slab. The simulated in-plane MTF response of each technique was also combined with the square root of the ENNPS estimate to yield stochastic signal-to-noise ratio (SNR) information about the different reconstruction techniques. For scan angles of 20° and 5 mm plane separation, seven MITS planes must be averaged to sufficiently remove partial-pixel artifacts. MITSa7 does appear to subtly reduce the contrast of high-frequency "edge" information, but the removal of partial-pixel artifacts makes the appearance of low-contrast, fine-detail anatomy even more conspicuous in MITSa7 slices. MITSa7 also appears to render simulated subtle 5 mm pulmonary nodules with greater visibility than MITS alone, in both the open lung and regions overlying the mediastinum. Finally, the MITSa7 technique reduces stochastic image variance, though the in-plane stochastic SNR (for very thin objects which do not span multiple MITS planes) is only improved at spatial frequencies between 0.05 and 0.20 cycles∕mm. The MITSa7 method is an improvement over traditional single-plane MITS for thoracic imaging and the pulmonary nodule detection task, and thus the authors plan to use the MITSa7 approach for all future MITS research at the authors' institution.

Measurements of system sharpness for two digital breast tomosynthesis systems

The aim of this work was to propose system sharpness parameters for digital breast tomosynthesis (DBT) systems that include the influence of focus size and focus motion for use in quality assurance protocols. X-ray focus size was measured using a multiple pinhole test object, while detector presampling modulation transfer function (MTF) was measured from projection images of a 10 cm × 10 cm, 1 mm thick steel edge, for the Siemens Inspiration and Hologic Selenia Dimensions DBT systems. The height of the edge above the table was then varied from 1 to 78 mm. The MTF expected from theory for the projection images was calculated from the measured detector MTF, focus size MTF and focus motion MTF and was compared against measured curves. Two methods were used to measure the in-plane MTF in the DBT volume: a tungsten wire of diameter 25 µm and an Al edge 0.2 mm thick, both imaged with a 15 mm thick poly(methyl methacrylate) (PMMA) plate. The in-depth point spread function (PSF) was measured using an angled tungsten wire. The full 3D MTF was estimated with a 0.5 mm diameter aluminium bead held in a 45 mm thick PMMA phantom, with the bead 15 and 65 mm above the table. Inspiration DBT projection images are saved at native detector resolution (85 µm), while the Dimensions re-bins projections to 140 µm pixels (2 × 2 binning); both systems used 2 × 2 binning of projection data before reconstruction. The 50% point for the MTF (MTF0.50) measured in the DBT projection images for the tube-travel direction fell as a function of height above the table from 3.60 to 0.90 mm−1 for the Inspiration system and from 2.50 to 1.20 mm−1 for the Dimensions unit. The maximum deviation of measured MTF0.50 from the calculated value was 13%. MTF0.50 measured in-plane (tube-travel direction) fell as a function of height above the table from 1.66 to 0.97 mm−1 for the Inspiration system and from 2.21 to 1.31 mm−1 for the Dimensions system. The full-width half-maximum for the in-depth PSF was 3.0 and 5.9 mm for the Inspiration and Dimensions systems, respectively. There was no difference in the 3D MTF curves, sectioned in the tube-travel direction, for bead heights of 15 and 65 mm above the table. A 25 µm tungsten wire held within a 15 mm thick PMMA plate was found to be a suitable test object for measurement of in-plane MTF. Evaluation of MTF as a function of height above the table, both in the projection images and in the reconstructed planes, provides important information on the impact of focus size and focus motion on the DBT system's imaging performance.

Evaluation of Image Quality of Multiplanar Reconstruction Images: Effect of Z-increment of Original Axial Images

The aim of our study was to evaluate the image quality of multiplanar reconstruction images (MPRs) focusing on the effect of z-increment of original axial images using signal to noise ratio (SNR) measurement in in- plane and longitudinal directions. SNRs of MPRs were calculated using modulation transfer function (MTF) and noise power spectrum (NPS). We scanned a bead phantom with a diameter of 0.1 mm and a water phantom with a diameter of 250 mm for calculating MTF and NPS using a MDCT with 0.5 s per rotation, 1.0 pitch and 64 x 0.6 mm collimation, and 50 mm field of view. Axial images for generating MPRs were reconstructed with standard kernel (B40), and 1.00 mm slice width. Coronal images were generated from two datasets with 0.1 mm and 0.5 mm z-increments of axial images respectively. For measuring the SNRs, the MTFs and NPSs in in-plane and longitudinal directions of each dataset were calculated from coronal bead images and coronal uniform noise images, respectively. Differences of MTF, NPS, and SNR were compared in in-plane and longitudinal directions. The MTF of longitudinal direction of the dataset with 0.1 mm z-increment was higher than the dataset with 0.5 mm z-increment. 10% MTFs of longitudinal direction with 0.1 mm and 0.5 mm z-increments were 0.75 cycles/mm and 0.68 cycles/mm, respectively. Conversely, the NPS of longitudinal direction of the dataset with a 0.1 mm z-increment was lower than the dataset with a 0.5 mm z-increment. As a consequence, the SNRs of longitudinal direction had relatively no difference between the datasets. In in-plane direction, MTFs, NPSs and SNRs had no differences between the datasets. A tradeoff relationship was indicated between spatial resolution and noise characteristic in the longitudinal direction due to the effects of different z-increment of original axial images used in generating MPRs. MPR using 0.5 mm z-increment of axial images had comparable SNR to MPR using 0.1 mm z-increment of axial images in our experimental condition. CLINICAL RELEVANCE/APPLICATION: Using 0.5 mm z-increment of original axial images for generating MPRs is effective for reducing the data volume, reconstruction time and transfer time without reducing image quality.

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.

TH-E-L100J-05: A Three-Dimensional Linear Model for Investigating the Image Quality of Breast Tomosynthesis

Purpose: To build a three dimensional (3D) linear model for breast tomosynthesis in order to investigate the effects of different imaging system design parameters on the reconstructed tomosynthesis image quality. Method and Materials: The 3D model incorporates the detector performance, imaging geometry and reconstruction filters for the filtered backprojection (FBP) method. To validate the model, experiments were performed on a prototype breast tomosynthesis system equipped with an amorphous selenium (a-Se) digital mammography detector. The detector can be operated in full resolution with 85 micron pixel size or binning mode to reduce acquisition time. Twenty—five projection images were acquired with an angular range of ±20°. The images were reconstructed using a slice thickness of 1 mm with 0.085 mm × 0.085 mm in-plane pixel dimension. The 3D noise power spectrum (NPS) was computed using reconstructed scatter free uniform images and an edge phantom was imaged to measure the in-plane modulation transfer function (MTF). An ACR mammography accreditation phantom was imaged to demonstrate the effects of detector operation modes and image acquisition geometry, i.e., angular range, on reconstructed image quality. Results: The measured in-plane MTF and 3D NPS are both in good agreement with the model. Our results showed that both NPS and in-plane MTF depend on the reconstruction filters and angular range. Increasing angular range helps improve the MTF at low frequencies, resulting in better detection of the large-area, low-contrast mass lesions in the phantom. Pixel binning cuts the acquisition time in half at the cost of decreased high frequency response. This results in reduced visibility of small calcification specks. Conclusion: A 3D linear system model for breast tomosynthesis was developed and validated with experimental measurements. The model can be used to predict the effects of system design parameters on reconstructed image quality. *Research sponsored by Siemens AG.

Novel ultrahigh resolution data acquisition and image reconstruction for multi‐detector row CT

We present and evaluate a special ultrahigh resolution mode providing considerably enhanced spatial resolution both in the scan plane and in the z-axis direction for a routine medical multi-detector row computed tomography (CT) system. Data acquisition is performed by using a flying focal spot both in the scan plane and in the z-axis direction in combination with tantalum grids that are inserted in front of the multi-row detector to reduce the aperture of the detector elements both in-plane and in the z-axis direction. The dose utilization of the system for standard applications is not affected, since the grids are moved into place only when needed and are removed for standard scanning. By means of this technique, image slices with a nominal section width of 0.4 mm (measured full width at half maximum=0.45 mm) can be reconstructed in spiral mode on a CT system with a detector configuration of 32 x 0.6 mm. The measured 2% value of the in-plane modulation transfer function (MTF) is 20.4 lp/cm, the measured 2% value of the longitudinal (z axis) MTF is 21.5 lp/cm. In a resolution phantom with metal line pair test patterns, spatial resolution of 20 lp/cm can be demonstrated both in the scan plane and along the z axis. This corresponds to an object size of 0.25 mm that can be resolved. The new mode is intended for ultrahigh resolution bone imaging, in particular for wrists, joints, and inner ear studies, where a higher level of image noise due to the reduced aperture is an acceptable trade-off for the clinical benefit brought about by the improved spatial resolution.

A computer simulation platform for the optimization of a breast tomosynthesis system

In breast tomosynthesis there is a compromise between resolution, noise, and acquisition speed for a given glandular dose. The purpose of the present work is to develop a simulation platform to investigate the potential imaging performance for the many possible tomosynthesis system configurations. The simulation platform was used to investigate the dependence of image blur and signal difference to noise ratio (SDNR) for several different tomosynthesis acquisition configurations. Simulated projections of a slanted thin tungsten wire placed in different object planes were modified according to the detector's modulation transfer function (MTF), with or without pixel binning. In addition, the focal spot blur (FSB), which depends on the location of the wire, the system geometry, the source-detector movement speed, and the exposure time, was also incorporated into the projections. Both expectation maximization (EM) and filtered back projection (FBP) based algorithms were used for 3D image reconstruction. The in-plane MTF was calculated from the reconstructed image of the wire. To evaluate the noise performance, simulated noiseless projections of calcification and tumor in uniform breast tissue were modified with the noise power spectrum (NPS) calculated from a cascaded linear system model for the detector for a given x-ray dose. The SDNR of the reconstructed images was calculated with different tomosynthesis configurations, e.g., pixel binning, view number, and angular range. Our results showed that for a source-to-imager distance (SID) of 66 cm, pixel binning (2 x 2) caused more degradation to the in-plane MTF than the blur caused by the moving focal spot and reconstruction. The in-depth resolution can be improved by increasing the angular range.

Performance evaluation of local ROI algorithms for exact ROI reconstruction in spiral cone-beam computed tomography

Investigates the performance of exact reconstruction algorithms for spiral cone-beam computed tomography. The authors compare two different approaches of exact algorithms: the first one is the Radon algorithm with explicit calculation of Radon data and subsequent direct inversion, the second approach is a reformulation of the Radon algorithm as a three-dimensional (3-D) filtered backprojection (FBP) method. Both algorithms are able to reconstruct a region of interest (ROI) of a long object by using the local ROI technique. The goal of this study is to identify which algorithm has the better performance with respect to image quality and practicability. The properties of the algorithms are analyzed. Image quality is evaluated regarding artifacts, spatial as well as contrast resolution, and noise properties. The 3-D point spread function (PSF) is examined for violations of isotropy and spatial invariance. In-plane and axial resolutions are derived from the in-plane modulation transfer function (MTF) and the slice sensitivity profile (SSP), respectively. The results show that for the two implementations investigated in this study, the image quality is similar, with a slight advantage of the FBP case. The algorithms yield images almost free of artifacts but with a poor spatial resolution. Regarding practicability, however, the FBP approach seems to be a more convenient and more flexible method because of the sequential processing of the projections and the easy control of image characteristics.