X-ray scatter correction method for dedicated breast computed tomography

Single-shot quantitative x-ray imaging using a primary modulator and dual-layer detector

BACKGROUND

Conventional x-ray imaging and fluoroscopy have limitations in quantitation due to several challenges, including scatter, beam hardening, and overlapping tissues. Dual-energy (DE) imaging, with its capability to quantify area density of specific materials, is well-suited to address such limitations, but only if the dual-energy projections are acquired with perfect spatial and temporal alignment and corrected for scatter.

PURPOSE

In this work, we propose single-shot quantitative imaging (SSQI) by combining the use of a primary modulator (PM) and dual-layer (DL) detector, which enables motion-free DE imaging with scatter correction in a single exposure.

METHODS

The key components of our SSQI setup include a PM and DL detector, where the former enables scatter correction for the latter while the latter enables beam hardening correction for the former. The SSQI algorithm allows simultaneous recovery of two material-specific images and two scatter images using four sub-measurements from the PM encoding. The concept was first demonstrated using simulation of chest x-ray imaging for a COVID patient. For validation, we set up SSQI geometry on our tabletop system and imaged acrylic and copper slabs with known thicknesses (acrylic: 0-22.5 cm; copper: 0-0.9 mm), estimated scatter with our SSQI algorithm, and compared the material decomposition (MD) for different combinations of the two materials with ground truth. Second, we imaged an anthropomorphic chest phantom containing contrast in the coronary arteries and compared the MD with and without SSQI. Lastly, to evaluate SSQI in dynamic applications, we constructed a flow phantom that enabled dynamic imaging of iodine contrast.

RESULTS

Our simulation study demonstrated that SSQI led to accurate scatter correction and MD, particularly for smaller focal blur and finer PM pitch. In the validation study, we found that the root mean squared error (RMSE) of SSQI estimation was 0.13 cm for acrylic and 0.04 mm for copper. For the anthropomorphic phantom, direct MD resulted in incorrect interpretation of contrast and soft tissue, while SSQI successfully distinguished them quantitatively, reducing RMSE in material-specific images by 38%-92%. For the flow phantom, SSQI was able to perform accurate dynamic quantitative imaging, separating contrast from the background.

CONCLUSIONS

We demonstrated the potential of SSQI for robust quantitative x-ray imaging. The integration of SSQI is straightforward with the addition of a PM and upgrade to a DL detector, which may enable its widespread adoption, including in techniques such as radiography and dynamic imaging (i.e., real-time image guidance and cone-beam CT).

Evaluation of the feasibility of a multisource CBCT for maxillofacial imaging

Objective. The aim of this study was to investigate the feasibility of improving the image quality and accuracy of cone beam computed tomography (CBCT) by replacing the conventional wide cone angle x-ray tube with a distributed x-ray source array positioned in the axial direction.Approach. The multisource CBCT (ms-CBCT) design was experimentally simulated using a benchtop scanner with a carbon nanotube x-ray tube and a flat-panel detector. The source was collimated and translated in the axial direction to simulate a source array with a reduced cone angle for each beam. An adjacent scatter ratio subtraction (ASRS) method was implemented for residual scatter reduction. Several phantoms were imaged using the ms-CBCT and conventional CBCT configurations under otherwise similar conditions. The Requirements of the ms-CBCT design on the x-ray source and detector were evaluated.Main results. Compared to the conventional CBCT, the ms-CBCT design with 8 sources and ASRS significantly improved the image quality and accuracy, including: (1) reducing the cupping artifact from 15% to 3.5%; (2) reducing the spatial nonuniformity of the CT Hounsfield unit values from 38.0 to 9.2; (3) improving the contrast-to-noise ratio of the low contrast objects (acrylic and low density polyethylene inserts) against the water-equivalent background by ∼20% and (4) reducing the root-mean-square error of the HU values by 70%, from 420.1 to 124.4. The imaging dose and scanning time used by the current clinical CBCT for maxillofacial imaging can be achieved by current source and detector technologies.Significance. The ms-CBCT design significantly reduces the scatter and improves the image quality and accuracy compared to the conventional CBCT.

A novel beam stopper-based approach for scatter correction in digital planar radiography

X-ray scatter in planar radiography degrades the contrast resolution of the image, thus reducing its diagnostic utility. Antiscatter grids partially block scattered photons at the cost of increasing the dose delivered by two- to four-fold and posing geometrical restrictions that hinder their use for other acquisition settings, such as portable radiography. The few software-based approaches investigated for planar radiography mainly estimate the scatter map from a low-frequency version of the image. We present a novel method for scatter correction in planar imaging based on direct patient measurements. Samples from the shadowed regions of an additional partially obstructed projection acquired with a beam stopper placed between the X-ray source and the patient are used to estimate the scatter map. Evaluation with simulated and real data showed an increase in contrast resolution for both lung and spine and recovery of ground truth values superior to those of three recently proposed methods. Our method avoids the biases of post-processing methods and yields results similar to those for an antiscatter grid while removing geometrical restrictions at around half the radiation dose. It can be used in unconventional imaging techniques, such as portable radiography, where training datasets needed for deep-learning approaches would be very difficult to obtain.

Deep learning for x-ray scatter correction in dedicated breast CT

BACKGROUND

Accurate correction of x-ray scatter in dedicated breast computed tomography (bCT) imaging may result in improved visual interpretation and is crucial to achieve quantitative accuracy during image reconstruction and analysis.

PURPOSE

To develop a deep learning (DL) model to correct for x-ray scatter in bCT projection images.

METHODS

A total of 115 patient scans acquired with a bCT clinical system were segmented into the major breast tissue types (skin, adipose, and fibroglandular tissue). The resulting breast phantoms were divided into training (n = 110) and internal validation cohort (n = 5). Training phantoms were augmented by a factor of four by random translation of the breast in the image field of view. Using a previously validated Monte Carlo (MC) simulation algorithm, 12 primary and scatter bCT projection images with a 30-degree step were generated from each phantom. For each projection, the thickness map and breast location in the field of view were also calculated. A U-Net based DL model was developed to estimate the scatter signal based on the total input simulated image and trained single-projection-wise, with the thickness map and breast location provided as additional inputs. The model was internally validated using MC-simulated projections and tested using an external data set of 10 phantoms derived from images acquired with a different bCT system. For this purpose, the mean relative difference (MRD) and mean absolute error (MAE) were calculated. To test for accuracy in reconstructed images, a full bCT acquisition was mimicked with MC-simulations and then assessed by calculating the MAE and the structural similarity (SSIM). Subsequently, scatter was estimated and subtracted from the bCT scans of three patients to obtain the scatter-corrected image. The scatter-corrected projections were reconstructed and compared with the uncorrected reconstructions by evaluating the correction of the cupping artifact, increase in image contrast, and contrast-to-noise ratio (CNR).

RESULTS

The mean MRD and MAE across all cases (min, max) for the internal validation set were 0.04% (-1.1%, 1.3%) and 2.94% (2.7%, 3.2%), while for the external test set they were -0.64% (-1.6%, 0.2%) and 2.84% (2.3%, 3.5%), respectively. For MC-simulated reconstruction slices, the computed SSIM was 0.99 and the MAE was 0.11% (range: 0%, 0.35%) with a single outlier slice of 2.06%. For the three patient bCT reconstructed images, the correction increased the contrast by a mean of 25% (range: 20%, 30%), and reduced the cupping artifact. The mean CNR increased by 0.32 after scatter correction, which was not found to be significant (95% confidence interval: [-0.01, 0.65], p = 0.059). The time required to correct the scatter in a single bCT projection was 0.2 s on an NVIDIA GeForce GTX 1080 GPU.

CONCLUSION

The developed DL model could accurately estimate scatter in bCT projection images and could enhance contrast and correct for cupping artifact in reconstructed patient images without significantly affecting the CNR. The time required for correction would allow its use in daily clinical practice, and the reported accuracy will potentially allow quantitative reconstructions.

Cone-beam breast CT using an offset detector: effect of detector offset and image reconstruction algorithm

Objective.A dedicated cone-beam breast computed tomography (BCT) using a high-resolution, low-noise detector operating in offset-detector geometry has been developed. This study investigates the effects of varying detector offsets and image reconstruction algorithms to determine the appropriate combination of detector offset and reconstruction algorithm.Approach.Projection datasets (300 projections in 360°) of 30 breasts containing calcified lesions that were acquired using a prototype cone-beam BCT system comprising a 40 × 30 cm flat-panel detector with 1024 × 768 detector pixels were reconstructed using Feldkamp-Davis-Kress (FDK) algorithm and served as the reference. The projection datasets were retrospectively truncated to emulate cone-beam datasets with sinograms of 768×768 and 640×768 detector pixels, corresponding to 5 cm and 7.5 cm lateral offsets, respectively. These datasets were reconstructed using the FDK algorithm with appropriate weights and an ASD-POCS-based Fast, total variation-Regularized, Iterative, Statistical reconstruction Technique (FRIST), resulting in a total of 4 offset-detector reconstructions (2 detector offsets × 2 reconstruction methods). Signal difference-to-noise ratio (SDNR), variance, and full-width at half-maximum (FWHM) of calcifications in two orthogonal directions were determined from all reconstructions. All quantitative measurements were performed on images in units of linear attenuation coefficient (1/cm).Results.The FWHM of calcifications did not differ (P > 0.262) among reconstruction algorithms and detector formats, implying comparable spatial resolution. For a chosen detector offset, the FRIST algorithm outperformed FDK in terms of variance and SDNR (P < 0.0001). For a given reconstruction method, the 5 cm offset provided better results.Significance.This study indicates the feasibility of using the compressed sensing-based, FRIST algorithm to reconstruct sinograms from offset-detectors. Among the reconstruction methods and detector offsets studied, FRIST reconstructions corresponding to a 30 cm × 30 cm with 5 cm lateral offset, achieved the best performance. A clinical prototype using such an offset geometry has been developed and installed for clinical trials.

A novel hardware duo of beam modulation and shielding to reduce scatter acquisition and dose in cone-beam breast CT

PURPOSE

In cone-beam breast CT, scattered photons form a large portion of the acquired signal, adversely impacting image quality throughout the frequency response of the imaging system. Prior simulation studies provided proof of concept for utilization of a hardware solution to prevent scatter acquisition. Here, we report the design, implementation, and characterization of an auxiliary apparatus of fluence modulation and scatter shielding that does indeed lead to projections with a reduced level of scatter.

METHODS

An apparatus was designed for permanent installation within an existing cone-beam CT system. The apparatus is composed of two primary assemblies: a "Fluence Modulator" (FM) and a "Scatter Shield" (SS). The design of the assemblies enables them to operate in synchrony during image acquisition, converting the sourced x-rays into a moving narrow beam. During a projection, this narrow beam sweeps the entire fan angle coverage of the imaging system. As the two assemblies are contingent on one another, their joint implementation is described in the singular as apparatus FM-SS. The FM and the SS assemblies are each comprised a metal housing, a sensory system, and a robotic system. A controller unit handles their relative movements. A series of comparative studies were conducted to evaluate the performance of a cone-beam CT system in two "modes" of operation: with and without FM-SS installed, and to compare the results of physical implementation with those previously simulated. The dynamic range requirements of the utilized detector in the cone-beam CT imaging system were first characterized, independent of the mode of operation. We then characterized and compared the spatial resolution of the imaging system with, and without, FM-SS. A physical breast phantom, representative of an average size breast, was developed and imaged. Actual differences in signal level obtained with, versus without, FM-SS were then compared to the expected level gains based on previously reported simulations. Following these initial assessments, the scatter acquisition in each projection in both modes of operation was investigated. Finally, as an initial study of the impact of FM-SS on radiation dose in an average size breast, a series of Monte Carlo simulations were coupled with physical measurements of air kerma, with and without FM-SS.

RESULTS

With implementation of FM-SS, the detector's required dynamic range was reduced by a factor of 5.5. Substantial reduction in the acquisition of the scattered rays, by a factor of 5.1 was achieved. With the implementation of FM-SS, deposited dose was reduced by 27% in the studied breast.

CONCLUSIONS

The disclosed implementation of FM-SS, within a cone-beam breast CT system, results in reduction of scatter-components in acquired projections, reduction of dose deposit to the breast, and relaxation of requirements for the detector's dynamic range. Controlling or correcting for patient motion occurring during image acquisition remains an open problem to be solved prior to practical clinical usage of FM-SS cone-beam breast CT.

Flat-panel conebeam CT in the clinic: history and current state

Research into conebeam CT concepts began as soon as the first clinical single-slice CT scanner was conceived. Early implementations of conebeam CT in the 1980s focused on high-contrast applications where concurrent high resolution ( < 200    μ m ), for visualization of small contrast-filled vessels, bones, or teeth, was an imaging requirement that could not be met by the contemporaneous CT scanners. However, the use of nonlinear imagers, e.g., x-ray image intensifiers, limited the clinical utility of the earliest diagnostic conebeam CT systems. The development of consumer-electronics large-area displays provided a technical foundation that was leveraged in the 1990s to first produce large-area digital x-ray detectors for use in radiography and then compact flat panels suitable for high-resolution and high-frame-rate conebeam CT. In this review, we show the concurrent evolution of digital flat panel (DFP) technology and clinical conebeam CT. We give a brief summary of conebeam CT reconstruction, followed by a brief review of the correction approaches for DFP-specific artifacts. The historical development and current status of flat-panel conebeam CT in four clinical areas-breast, fixed C-arm, image-guided radiation therapy, and extremity/head-is presented. Advances in DFP technology over the past two decades have led to improved visualization of high-contrast, high-resolution clinical tasks, and image quality now approaches the soft-tissue contrast resolution that is the standard in clinical CT. Future technical developments in DFPs will enable an even broader range of clinical applications; research in the arena of flat-panel CT shows no signs of slowing down.

Advancements towards the implementation of clinical phase-contrast breast computed tomography at Elettra

Breast computed tomography (BCT) is an emerging application of X-ray tomography in radiological practice. A few clinical prototypes are under evaluation in hospitals and new systems are under development aiming at improving spatial and contrast resolution and reducing delivered dose. At the same time, synchrotron-radiation phase-contrast mammography has been demonstrated to offer substantial advantages when compared with conventional mammography. At Elettra, the Italian synchrotron radiation facility, a clinical program of phase-contrast BCT based on the free-space propagation approach is under development. In this paper, full-volume breast samples imaged with a beam energy of 32 keV delivering a mean glandular dose of 5 mGy are presented. The whole acquisition setup mimics a clinical study in order to evaluate its feasibility in terms of acquisition time and image quality. Acquisitions are performed using a high-resolution CdTe photon-counting detector and the projection data are processed via a phase-retrieval algorithm. Tomographic reconstructions are compared with conventional mammographic images acquired prior to surgery and with histologic examinations. Results indicate that BCT with monochromatic beam and free-space propagation phase-contrast imaging provide relevant three-dimensional insights of breast morphology at clinically acceptable doses and scan times.

Phantom-based study exploring the effects of different scatter correction approaches on the reconstructed images generated by contrast-enhanced stationary digital breast tomosynthesis

Stationary digital breast tomosynthesis (sDBT) is an emerging technology in which the single rotating x-ray tube is replaced by a fixed array of multiple carbon nanotube-enabled sources, providing a higher spatial and temporal resolution. As such, sDBT offers a promising platform for contrast-enhanced (CE) imaging. However, given the minimal enhancement above background with standard operational tube settings and iodine dosing, CE breast imaging requires additional acquisition steps to isolate the iodine signal, using either temporal or dual energy subtraction (TS or DES) protocols. Also, correcting for factors that limit contrast is critical, and scatter and noise pose unique challenges during tomosynthesis. This phantom-based study of CE sDBT compared different postacquisition scatter correction approaches on the quality of the reconstructed image slices. Beam-pass collimation was used to sample scatter indirectly, from which an interpolated scatter map was obtained for each projection image. Scatter-corrected projections provided the information for reconstruction. Comparison between the application of different scatter maps demonstrated the significant effect that processing has on the contrast-to-noise ratio and feature detectability ([Formula: see text]) in the final displayed images and emphasized the critical importance of scatter correction during DES.

Characterization of CT Hounsfield Units for 3D acquisition trajectories on a dedicated breast CT system

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is important, especially when using multiple sources, novel detector and non-traditional trajectories. Dedicated clinical breast CT (BCT) systems therefore should be similarly evaluated. In this study, uniform cylindrical phantoms filled with various uniform density fluids were used to characterize differences in HU values between simple circular and complex 3D (saddle) orbits. Based on ACR recommendations, the HU accuracy, center-to-edge variability within a slice, and overall variability within the reconstructed volume were characterized for simple and complex acquisitions possible on a single versatile BCT system. Results illustrate the statistically significantly better performance of the saddle orbit, especially close to the chest and nipple regions of what would clinically be a pendant breast volume. The incomplete cone beam acquisition of a simple circular orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values throughout the reconstructed volumes. This study further establishes the advantages of using 3D acquisition trajectories for breast CT as well as other applications by demonstrating the robustness of HU values throughout large reconstructed volumes.

The role of off-focus radiation in scatter correction for dedicated cone beam breast CT

PURPOSE

Dedicated cone beam breast CT (CBBCT) suffers from x-ray scatter contamination. We aim to identify the source of the significant difference between the scatter distributions estimated by two recent methods proposed by our group and to investigate its effect on CBBCT image quality.

METHOD

We recently proposed two novel methods of scatter correction for CBBCT, using a library based (LB) technique and a forward projection (FP) model. Despite similar enhancement on CBBCT image qualities, these two methods obtain very different scatter distributions. We hypothesize that the off-focus radiation (OFR) is the contributor and results in nontrivial signals in x-ray projections, which is ignored in the scatter estimation via the LB method. Experiments using a thin wire test tool are designed to study the effect of OFR on CBBCT spatial resolution by measuring the point spread function (PSF) and the modulation transfer function (MTF). A narrow collimator setting is used to suppress the OFR-induced signals. In addition, "PSFs" and "MTFs" are measured on clinical CBBCT images obtained by the LB and FP methods using small calcifications as point sources. The improvement of spatial resolution achieved by suppressing OFR in the wire experiment as well as in the clinical study is quantified by the improvement ratios of PSFs and spatial frequencies at different MTF values. Our hypothesis that OFR causes the imaging difference between the FP and LB methods is verified if these ratios obtained from experimental and clinical data are consistent.

RESULTS

In the wire experiment, the results show that suppression of OFR increases the maximum signal of the PSF by about 14% and reduces the full-width-at-half-maximum (FWHM) by about 12.0%. Similar improvement on spatial resolution is achieved by the FP method compared with the LB method in the patient study. The improvement ratios of spatial frequencies at different MTF values without OFR match very well in both studies at a level of around 16%, with an average root-mean-square difference of 0.47%.

CONCLUSION

The results of the wire experiment and the clinical study indicate that the main difference between the LB and FP methods is whether the OFR-induced signals are included after scatter correction. Our study further shows that OFR significantly affects the image spatial resolution of CBBCT, indicating that the visualization of micro-calcifications is susceptible to OFR contamination. Our finding is therefore important in further improvement of diagnostic performance of CBBCT.

Implementation and CT sampling characterization of a third-generation SPECT-CT system for dedicated breast imaging

Stand-alone cone beam computed tomography (CT) and single-photon emission computed tomography (SPECT) systems capable of complex acquisition trajectories have previously been developed for breast imaging. Fully three-dimensional (3-D) motions of SPECT systems provide views into the chest wall and throughout the entire volume. The polar tilting capability of the CBCT system has shown improvement in sampling close to the chest wall, while eliminating cone beam artifacts. Here, a single hybrid SPECT-CT system, with each individual modality capable of independently traversing complex trajectories around a common pendant breast volume, was developed. We present the practical implementation of this design and preliminary results of the CT system. The fully 3-D SPECT was nested inside the suspended CT gantry and oriented perpendicular to the CT source-detector pair. Both subsystems were positioned on a rotation stage, with the combined polar and azimuthal motions enabling spherical trajectories. Six trajectories were used for initial evaluation of the tilt capable CT system. The developed system can achieve polar tilt angles with a [Formula: see text] positioning error and no hysteresis. Initial imaging results demonstrate that additional off-axis projection views of various geometric resolution phantoms facilitate more complete sampling, more consistent attenuation value recovery, and markedly improved reconstructions. This system could have various applications in diagnostic or therapeutic breast imaging.

Transmission characteristics of a two dimensional antiscatter grid prototype for CBCT

AIM

High fraction of scattered radiation in cone-beam CT (CBCT) imaging degrades CT number accuracy and visualization of low contrast objects. To suppress scatter in CBCT projections, we developed a focused, two-dimensional antiscatter grid (2DASG) prototype. In this work, we report on the primary and scatter transmission characteristics of the 2DASG prototype aimed for linac mounted, offset detector geometry CBCT systems in radiation therapy, and compared its performance to a conventional one-dimensional ASG (1DASG).

METHODS

The 2DASG is an array of through-holes separated by 0.1 mm septa that was fabricated from tungsten using additive manufacturing techniques. Through-holes' focusing geometry was designed for offset detector CBCT in Varian TrueBeam system. Two types of ASGs were evaluated: (a) a conventional 1DASG with a grid ratio of 10, (b) the 2DASG prototype with a grid ratio of 8.2. To assess the scatter suppression performance of both ASGs, Scatter-to-primary ratio (SPR) and scatter transmission fraction (Ts) were measured using the beam stop method. Scatter and primary intensities were modulated by varying the phantom thickness between 10 and 40 cm. Additionally, the effect of air gap and bow tie (BT) filter on SPR and Ts were evaluated. Average primary transmission fraction (T_P ) and pixel specific primary transmission were also measured for both ASGs. To assess the effect of transmission characteristics on projection image signal-to-noise ratio (SNR), SNR improvement factor was calculated. Improvement in contrast to noise ratio (CNR) was demonstrated using a low contrast object.

RESULTS

In comparison to 1DASG, 2DASG reduced SPRs by a factor of 3 to 6 across the range of phantom setups investigated. Ts values for 1D and 2DASGs were in the range of 21 to 29%, and 5 to 14% respectively. 2DASG continued to provide lower SPR and Ts at increased air gap and with BT filter. Tp of 1D and 2DASGs were 70.6% and 84.7% respectively. Due to the septal shadow of the 2DASG, its pixel specific primary transmission values varied between 32.5% and 99.1%. With respect to 1DASG, 2DASG provided up to factor of 1.7 more improvement in SNR across the SPR range investigated. Moreover, 2DASG provided improved visualization of low contrast objects with respect to 1DASG and NOASG setups.

CONCLUSIONS

When compared to a conventional 1DASG, 2DASG prototype provided noticeably lower SPR and Ts values, indicating its superior scatter suppression performance. 2DASG also provided 19% higher average primary transmission that was attributed to the absence of interseptal spacers and optimized grid geometry. Our results indicate that the combined effect of lower scatter and higher primary transmission provided by 2DASG may potentially translate into more accurate CT numbers and improved contrast resolution in CBCT images.

Characterization of X-ray scattering for various phantoms and clinical breast geometries using breast CT on a dedicated hybrid system

OBJECTIVE

The purpose of this study was to utilize a dedicated breast CT system using a 2D beam stop array to physically evaluate the scatter to primary ratios (SPRs) of different geometric phantoms and prospectively acquired clinical patient data.

METHODS

Including clinically unrealizable compositions of 100% glandular and 100% fat, projection images were acquired using three geometrically different phantoms filled with fluids simulating breast tissue. The beam stop array method was used for measuring scatter in projection space, and creating the scatter corrected primary images. 2D SPRs were calculated. Additionally, a new figure of merit, the 3D normalized scatter contribution (NSC) volumes were calculated.

RESULTS

The 2D SPR values (0.52-1.10) were primarily dependent on phantom geometry; a secondary dependence was due to their uniform density; 2D SPRs were low frequency and smoothly varying in the uniformly filled phantoms. SPRs of clinical patient data followed similar trends as phantoms, but with noticeable deviations and high frequency components due to the heterogeneous distribution of glandular tissue. The maximum measured patient 2D SPRs were all <0.6, even for the largest diameter breast. These results demonstrate modest scatter components with changing object geometries and densities; the 3D NSC volumes with higher frequency components help visualize scatter distribution throughout the reconstructed image volumes. Furthermore, the SPRs in the heterogeneous clinical breast cases were underestimated by the equivalent density, uniformly filled phantoms.

CONCLUSIONS

These results provide guidance on the use of uniformly distributed density and differently shaped phantoms when considering simulations. They also clearly demonstrate that results from patients can vary considerably from 2D SPRs of uniformly simulated phantoms.

Estimating scatter from sparsely measured primary signal

Scatter radiation severely degrades the image quality. Measurement-based scatter correction methods sample the scatter signal at sparsely distributed points, from which the scatter profile is estimated and deterministically removed from the projection image. The estimation of the scatter profile is generally done through a spline interpolation and the resulting scatter profile is quite smooth. Consequently, the noise is intact and the signal-to-noise ratio is reduced in the projection image after scatter correction, leading to image artifacts and increased noise in the reconstruction images. We propose a simple and effective method, referred to as filtered scatter-to-primary ratio ([Formula: see text]-SPR) estimation, to estimate the scatter profile using the sparsely sampled scatter signal. Using the primary sampling device and the stationary digital tomosynthesis systems previously developed in our lab, we evaluated and compared the [Formula: see text]-SPR method in comparison with existing methods in terms of contrast ratio, signal difference-to-noise ratio, and modulation transfer function. A significant improvement in image quality is observed in both the projection and the reconstruction images using the proposed method.

Library based x-ray scatter correction for dedicated cone beam breast CT

PURPOSE

The image quality of dedicated cone beam breast CT (CBBCT) is limited by substantial scatter contamination, resulting in cupping artifacts and contrast-loss in reconstructed images. Such effects obscure the visibility of soft-tissue lesions and calcifications, which hinders breast cancer detection and diagnosis. In this work, we propose a library-based software approach to suppress scatter on CBBCT images with high efficiency, accuracy, and reliability.

METHODS

The authors precompute a scatter library on simplified breast models with different sizes using the geant4-based Monte Carlo (MC) toolkit. The breast is approximated as a semiellipsoid with homogeneous glandular/adipose tissue mixture. For scatter correction on real clinical data, the authors estimate the breast size from a first-pass breast CT reconstruction and then select the corresponding scatter distribution from the library. The selected scatter distribution from simplified breast models is spatially translated to match the projection data from the clinical scan and is subtracted from the measured projection for effective scatter correction. The method performance was evaluated using 15 sets of patient data, with a wide range of breast sizes representing about 95% of general population. Spatial nonuniformity (SNU) and contrast to signal deviation ratio (CDR) were used as metrics for evaluation.

RESULTS

Since the time-consuming MC simulation for library generation is precomputed, the authors' method efficiently corrects for scatter with minimal processing time. Furthermore, the authors find that a scatter library on a simple breast model with only one input parameter, i.e., the breast diameter, sufficiently guarantees improvements in SNU and CDR. For the 15 clinical datasets, the authors' method reduces the average SNU from 7.14% to 2.47% in coronal views and from 10.14% to 3.02% in sagittal views. On average, the CDR is improved by a factor of 1.49 in coronal views and 2.12 in sagittal views.

CONCLUSIONS

The library-based scatter correction does not require increase in radiation dose or hardware modifications, and it improves over the existing methods on implementation simplicity and computational efficiency. As demonstrated through patient studies, the authors' approach is effective and stable, and is therefore clinically attractive for CBBCT imaging.

Bowtie filters for dedicated breast CT: theory and computational implementation

PURPOSE

To design bowtie filters with improved properties for dedicated breast CT to improve image quality and reduce dose to the patient.

METHODS

The authors present three different bowtie filters designed for a cylindrical 14-cm diameter phantom with a uniform composition of 40/60 breast tissue, which vary in their design objectives and performance improvements. Bowtie design #1 is based on single material spectral matching and produces nearly uniform spectral shape for radiation incident upon the detector. Bowtie design #2 uses the idea of basis material decomposition to produce the same spectral shape and intensity at the detector, using two different materials. Bowtie design #3 eliminates the beam hardening effect in the reconstructed image by adjusting the bowtie filter thickness so that the effective attenuation coefficient for every ray is the same. All three designs are obtained using analytical computational methods and linear attenuation coefficients. Thus, the designs do not take into account the effects of scatter. The authors considered this to be a reasonable approach to the filter design problem since the use of Monte Carlo methods would have been computationally intensive. The filter profiles for a cone-angle of 0° were used for the entire length of each filter because the differences between those profiles and the correct cone-beam profiles for the cone angles in our system are very small, and the constant profiles allowed construction of the filters with the facilities available to us. For evaluation of the filters, we used Monte Carlo simulation techniques and the full cone-beam geometry. Images were generated with and without each bowtie filter to analyze the effect on dose distribution, noise uniformity, and contrast-to-noise ratio (CNR) homogeneity. Line profiles through the reconstructed images generated from the simulated projection images were also used as validation for the filter designs.

RESULTS

Examples of the three designs are presented. Initial verification of performance of the designs was done using analytical computations of HVL, intensity, and effective attenuation coefficient behind the phantom as a function of fan-angle with a cone-angle of 0°. The performance of the designs depends only weakly on incident spectrum and tissue composition. For all designs, the dynamic range requirement on the detector was reduced compared to the no-bowtie-filter case. Further verification of the filter designs was achieved through analysis of reconstructed images from simulations. Simulation data also showed that the use of our bowtie filters can reduce peripheral dose to the breast by 61% and provide uniform noise and CNR distributions. The bowtie filter design concepts validated in this work were then used to create a computational realization of a 3D anthropomorphic bowtie filter capable of achieving a constant effective attenuation coefficient behind the entire field-of-view of an anthropomorphic breast phantom.

CONCLUSIONS

Three different bowtie filter designs that vary in performance improvements were described and evaluated using computational and simulation techniques. Results indicate that the designs are robust against variations in breast diameter, breast composition, and tube voltage, and that the use of these filters can reduce patient dose and improve image quality compared to the no-bowtie-filter case.

X-ray scatter correction method for dedicated breast computed tomography: improvements and initial patient testing

A previously proposed x-ray scatter correction method for dedicated breast computed tomography was further developed and implemented so as to allow for initial patient testing. The method involves the acquisition of a complete second set of breast CT projections covering 360° with a perforated tungsten plate in the path of the x-ray beam. To make patient testing feasible, a wirelessly controlled electronic positioner for the tungsten plate was designed and added to a breast CT system. Other improvements to the algorithm were implemented, including automated exclusion of non-valid primary estimate points and the use of a different approximation method to estimate the full scatter signal. To evaluate the effectiveness of the algorithm, evaluation of the resulting image quality was performed with a breast phantom and with nine patient images. The improvements in the algorithm resulted in the avoidance of introduction of artifacts, especially at the object borders, which was an issue in the previous implementation in some cases. Both contrast, in terms of signal difference and signal difference-to-noise ratio were improved with the proposed method, as opposed to with the correction algorithm incorporated in the system, which does not recover contrast. Patient image evaluation also showed enhanced contrast, better cupping correction, and more consistent voxel values for the different tissues. The algorithm also reduces artifacts present in reconstructions of non-regularly shaped breasts. With the implemented hardware and software improvements, the proposed method can be reliably used during patient breast CT imaging, resulting in improvement of image quality, no introduction of artifacts, and in some cases reduction of artifacts already present. The impact of the algorithm on actual clinical performance for detection, diagnosis and other clinical tasks in breast imaging remains to be evaluated.

A breast-specific, negligible-dose scatter correction technique for dedicated cone-beam breast CT: a physics-based approach to improve Hounsfield Unit accuracy

The purpose of this research was to develop a method to correct the cupping artifact caused from x-ray scattering and to achieve consistent Hounsfield Unit (HU) values of breast tissues for a dedicated breast CT (bCT) system. The use of a beam passing array (BPA) composed of parallel-holes has been previously proposed for scatter correction in various imaging applications. In this study, we first verified the efficacy and accuracy using BPA to measure the scatter signal on a cone-beam bCT system. A systematic scatter correction approach was then developed by modeling the scatter-to-primary ratio (SPR) in projection images acquired with and without BPA. To quantitatively evaluate the improved accuracy of HU values, different breast tissue-equivalent phantoms were scanned and radially averaged HU profiles through reconstructed planes were evaluated. The dependency of the correction method on object size and number of projections was studied. A simplified application of the proposed method on five clinical patient scans was performed to demonstrate efficacy. For the typical 10-18 cm breast diameters seen in the bCT application, the proposed method can effectively correct for the cupping artifact and reduce the variation of HU values of breast equivalent material from 150 to 40 HU. The measured HU values of 100% glandular tissue, 50/50 glandular/adipose tissue, and 100% adipose tissue were approximately 46, -35, and -94, respectively. It was found that only six BPA projections were necessary to accurately implement this method, and the additional dose requirement is less than 1% of the exam dose. The proposed method can effectively correct for the cupping artifact caused from x-ray scattering and retain consistent HU values of breast tissues.

Personalized estimates of radiation dose from dedicated breast CT in a diagnostic population and comparison with diagnostic mammography

This study retrospectively analyzed the mean glandular dose (MGD) to 133 breasts from 132 subjects, all women, who participated in a clinical trial evaluating dedicated breast CT in a diagnostic population. The clinical trial was conducted in adherence to a protocol approved by institutional review boards and the study participants provided written informed consent. Individual estimates of MGD to each breast from dedicated breast CT was obtained by combining x-ray beam characteristics with estimates of breast dimensions and fibroglandular fraction from volumetric breast CT images, and using normalized glandular dose coefficients. For each study participant and for the breast corresponding to that imaged with breast CT, an estimate of the MGD from diagnostic mammography (including supplemental views) was obtained from the DICOM image headers for comparison. This estimate uses normalized glandular dose coefficients corresponding to a breast with 50% fibroglandular weight fraction. The median fibroglandular weight fraction for the study cohort determined from volumetric breast CT images was 15%. Hence, the MGD from diagnostic mammography was corrected to be representative of the study cohort. Individualized estimates of MGD from breast CT ranged from 5.7 to 27.8 mGy. Corresponding to the breasts imaged with breast CT, the MGD from diagnostic mammography ranged from 2.6 to 31.6 mGy. The mean (± inter-breast SD) and the median MGD (mGy) from dedicated breast CT exam were 13.9 ± 4.6 and 12.6, respectively. For the corresponding breasts, the mean (± inter-breast SD) and the median MGD (mGy) from diagnostic mammography were 12.4 ± 6.3 and 11.1, respectively. Statistical analysis indicated that at the 0.05 level, the distributions of MGD from dedicated breast CT and diagnostic mammography were significantly different (Wilcoxon signed ranks test, p = 0.007). While the interquartile range and the range (maximum-minimum) of MGD from dedicated breast CT was lower than diagnostic mammography, the median MGD from dedicated breast CT was approximately 13.5% higher than that from diagnostic mammography. The MGD for breast CT is based on a 1.45 mm skin layer and that for diagnostic mammography is based on a 4 mm skin layer; thus, favoring a lower estimate for MGD from diagnostic mammography. The median MGD from dedicated breast CT corresponds to the median MGD from four to five diagnostic mammography views. In comparison, for the same 133 breasts, the mean and the median number of views per breast during diagnostic mammography were 4.53 and 4, respectively. Paired analysis showed that there was approximately equal likelihood of receiving lower MGD from either breast CT or diagnostic mammography. Future work will investigate methods to reduce and optimize radiation dose from dedicated breast CT.

Modelling the physics in the iterative reconstruction for transmission computed tomography

There is an increasing interest in iterative reconstruction (IR) as a key tool to improve quality and increase applicability of x-ray CT imaging. IR has the ability to significantly reduce patient dose; it provides the flexibility to reconstruct images from arbitrary x-ray system geometries and allows one to include detailed models of photon transport and detection physics to accurately correct for a wide variety of image degrading effects. This paper reviews discretization issues and modelling of finite spatial resolution, Compton scatter in the scanned object, data noise and the energy spectrum. The widespread implementation of IR with a highly accurate model-based correction, however, still requires significant effort. In addition, new hardware will provide new opportunities and challenges to improve CT with new modelling.

Modelling the physics in the iterative reconstruction for transmission computed tomography

There is an increasing interest in iterative reconstruction (IR) as a key tool to improve quality and increase applicability of x-ray CT imaging. IR has the ability to significantly reduce patient dose; it provides the flexibility to reconstruct images from arbitrary x-ray system geometries and allows one to include detailed models of photon transport and detection physics to accurately correct for a wide variety of image degrading effects. This paper reviews discretization issues and modelling of finite spatial resolution, Compton scatter in the scanned object, data noise and the energy spectrum. The widespread implementation of IR with a highly accurate model-based correction, however, still requires significant effort. In addition, new hardware will provide new opportunities and challenges to improve CT with new modelling.

Fast parallel algorithms for the x-ray transform and its adjoint

PURPOSE

Iterative reconstruction methods often offer better imaging quality and allow for reconstructions with lower imaging dose than classical methods in computed tomography. However, the computational speed is a major concern for these iterative methods, for which the x-ray transform and its adjoint are two most time-consuming components. The speed issue becomes even notable for the 3D imaging such as cone beam scans or helical scans, since the x-ray transform and its adjoint are frequently computed as there is usually not enough computer memory to save the corresponding system matrix. The purpose of this paper is to optimize the algorithm for computing the x-ray transform and its adjoint, and their parallel computation.

METHODS

The fast and highly parallelizable algorithms for the x-ray transform and its adjoint are proposed for the infinitely narrow beam in both 2D and 3D. The extension of these fast algorithms to the finite-size beam is proposed in 2D and discussed in 3D.

RESULTS

The CPU and GPU codes are available at https://sites.google.com/site/fastxraytransform. The proposed algorithm is faster than Siddon's algorithm for computing the x-ray transform. In particular, the improvement for the parallel computation can be an order of magnitude.

CONCLUSIONS

The authors have proposed fast and highly parallelizable algorithms for the x-ray transform and its adjoint, which are extendable for the finite-size beam. The proposed algorithms are suitable for parallel computing in the sense that the computational cost per parallel thread is O(1).