Dataset Information

Scatter correction for full-fan volumetric CT using a stationary beam blocker in a single full scan.

ABSTRACT: Applications of volumetric CT (VCT) are hampered by shading and streaking artifacts in the reconstructed images. These artifacts are mainly due to strong x-ray scatter signals accompanied with the large illumination area within one projection, which lead to CT number inaccuracy, image contrast loss and spatial nonuniformity. Although different scatter correction algorithms have been proposed in literature, a standard solution still remains unclear. Measurement-based methods use a beam blocker to acquire scatter samples. These techniques have unrivaled advantages over other existing algorithms in that they are simple and efficient, and achieve high scatter estimation accuracy without prior knowledge of the imaged object. Nevertheless, primary signal loss is inevitable in the scatter measurement, and multiple scans or moving the beam blocker during data acquisition are typically employed to compensate for the missing primary data. In this paper, we propose a new measurement-based scatter correction algorithm without primary compensation for full-fan VCT. An accurate reconstruction is obtained with one single-scan and a stationary x-ray beam blocker, two seemingly incompatible features which enable simple and efficient scatter correction without increase of scan time or patient dose.Based on the CT reconstruction theory, we distribute the blocked data over the projection area where primary signals are considered approximately redundant in a full scan, such that the CT image quality is not degraded even with primary loss. Scatter is then accurately estimated by interpolation and scatter-corrected CT images are obtained using an FDK-based reconstruction algorithm.The proposed method is evaluated using two phantom studies on a tabletop CBCT system. On the Catphan©600 phantom, our approach reduces the reconstruction error from 207 Hounsfield unit (HU) to 9 HU in the selected region of interest, and improves the image contrast by a factor of 2.0 in the high-contrast regions. On an anthropomorphic head phantom, the reconstruction error is reduced from 97 HU to 6 HU in the soft tissue region and image spatial nonuniformity decreases from 27% to 5% after correction.Our method inherits the main advantages of measurement-based methods while avoiding their shortcomings. It has the potential to become a practical scatter correction solution widely implementable on different VCT systems.

SUBMITTER: Niu T

PROVIDER: S-EPMC3215690 | biostudies-other | 2011 Nov

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:X-ray scatter is a significant barrier to image quality improvements in cone-beam computed tomography (CBCT). A moving blocker-based strategy was previously proposed to simultaneously estimate scatter and reconstruct the complete volume within the field of view (FOV) from a single CBCT scan. A blocker consisting of lead stripes is inserted between the X-ray source and the imaging object, and moves back and forth along the rotation axis during gantry rotation. While promising results were obtained in our previous studies, the geometric design and moving speed of the blocker were set empirically. The goal of this work is to optimize the geometry and speed of the moving block system.Performance of the blocker was examined through Monte Carlo (MC) simulation and experimental studies with various geometry designs and moving speeds. All hypothetical designs employed an anthropomorphic pelvic phantom. The scatter estimation accuracy was quantified by using lead stripes ranging from 5 to 100 pixels on the detector plane. An iterative reconstruction based on total variation minimization was used to reconstruct CBCT images from unblocked projection data after scatter correction. The reconstructed image was evaluated under various combinations of lead strip width and interspace (ranging from 10 to 60 pixels) and different moving speed (ranging from 1 to 30 pixels per projection).MC simulation showed that the scatter estimation error varied from 0.8% to 5.8%. Phantom experiment showed that CT number error in the reconstructed CBCT images varied from 13 to 35. Highest reconstruction accuracy was achieved when the strip width was 20 pixels and interspace was 60 pixels and the moving speed was 15 pixels per projection.Scatter estimation can be achieved in a large range of lead strip width and interspace combinations. The moving speed does not have a very strong effect on reconstruction result if it is above 5 pixels per projection. Geometry design of the blocker affected image reconstruction accuracy more. The optimal geometry of the blocker has a strip width of 20 pixels and an interspace three times the strip width, which means 25% detector is covered by the blocker, while the optimal moving speed is 15 pixels per projection.

Project description:X-ray scatter results in a significant degradation of image quality in computed tomography (CT), representing a major limitation in cone-beam CT (CBCT) and large field-of-view diagnostic scanners. In this work, a novel scatter estimation and correction technique is proposed that utilizes peripheral detection of scatter during the patient scan to simultaneously acquire image and patient-specific scatter information in a single scan, and in conjunction with a proposed compressed sensing scatter recovery technique to reconstruct and correct for the patient-specific scatter in the projection space.The method consists of the detection of patient scatter at the edges of the field of view (FOV) followed by measurement based compressed sensing recovery of the scatter through-out the projection space. In the prototype implementation, the kV x-ray source of the Varian TrueBeam OBI system was blocked at the edges of the projection FOV, and the image detector in the corresponding blocked region was used for scatter detection. The design enables image data acquisition of the projection data on the unblocked central region of and scatter data at the blocked boundary regions. For the initial scatter estimation on the central FOV, a prior consisting of a hybrid scatter model that combines the scatter interpolation method and scatter convolution model is estimated using the acquired scatter distribution on boundary region. With the hybrid scatter estimation model, compressed sensing optimization is performed to generate the scatter map by penalizing the L1 norm of the discrete cosine transform of scatter signal. The estimated scatter is subtracted from the projection data by soft-tuning, and the scatter-corrected CBCT volume is obtained by the conventional Feldkamp-Davis-Kress algorithm. Experimental studies using image quality and anthropomorphic phantoms on a Varian TrueBeam system were carried out to evaluate the performance of the proposed scheme.The scatter shading artifacts were markedly suppressed in the reconstructed images using the proposed method. On the Catphan©504 phantom, the proposed method reduced the error of CT number to 13 Hounsfield units, 10% of that without scatter correction, and increased the image contrast by a factor of 2 in high-contrast regions. On the anthropomorphic phantom, the spatial nonuniformity decreased from 10.8% to 6.8% after correction.A novel scatter correction method, enabling unobstructed acquisition of the high frequency image data and concurrent detection of the patient-specific low frequency scatter data at the edges of the FOV, is proposed and validated in this work. Relative to blocker based techniques, rather than obstructing the central portion of the FOV which degrades and limits the image reconstruction, compressed sensing is used to solve for the scatter from detection of scatter at the periphery of the FOV, enabling for the highest quality reconstruction in the central region and robust patient-specific scatter correction.

Project description:PURPOSE:Scatter is a major factor degrading the image quality of cone beam computed tomography (CBCT). Conventional scatter correction strategies require handcrafted analytical models with ad hoc assumptions, which often leads to less accurate scatter removal. This study aims to develop an effective scatter correction method using a residual convolutional neural network (CNN). METHODS:A U-net based 25-layer CNN was constructed for CBCT scatter correction. The establishment of the model consists of three steps: model training, validation, and testing. For model training, a total of 1800 pairs of x-ray projection and the corresponding scatter-only distribution in nonanthropomorphic phantoms taken in full-fan scan were generated using Monte Carlo simulation of a CBCT scanner installed with a proton therapy system. An end-to-end CNN training was implemented with two major loss functions for 100 epochs with a mini-batch size of 10. Image rotations and flips were randomly applied to augment the training datasets during training. For validation, 200 projections of a digital head phantom were collected. The proposed CNN-based method was compared to a conventional projection-domain scatter correction method named fast adaptive scatter kernel superposition (fASKS) method using 360 projections of an anthropomorphic head phantom. Two different loss functions were applied for the same CNN to evaluate the impact of loss functions on the final results. Furthermore, the CNN model trained with full-fan projections was fine-tuned for scatter correction in half-fan scan by using transfer learning with additional 360 half-fan projection pairs of nonanthropomorphic phantoms. The tuned-CNN model for half-fan scan was compared with the fASKS method as well as the CNN-based method without the fine-tuning using additional lung phantom projections. RESULTS:The CNN-based method provides projections with significantly reduced scatter and CBCT images with more accurate Hounsfield Units (HUs) than that of the fASKS-based method. Root mean squared error of the CNN-corrected projections was improved to 0.0862 compared to 0.278 for uncorrected projections or 0.117 for the fASKS-corrected projections. The CNN-corrected reconstruction provided better HU quantification, especially in regions near the air or bone interfaces. All four image quality measures, which include mean absolute error (MAE), mean squared error (MSE), peak signal-to-noise ratio (PSNR), and structural similarity (SSIM), indicated that the CNN-corrected images were significantly better than that of the fASKS-corrected images. Moreover, the proposed transfer learning technique made it possible for the CNN model trained with full-fan projections to be applicable to remove scatters in half-fan projections after fine-tuning with only a small number of additional half-fan training datasets. SSIM value of the tuned-CNN-corrected images was 0.9993 compared to 0.9984 for the non-tuned-CNN-corrected images or 0.9990 for the fASKS-corrected images. Finally, the CNN-based method is computationally efficient - the correction time for the 360 projections only took less than 5 s in the reported experiments on a PC (4.20 GHz Intel Core-i7 CPU) with a single NVIDIA GTX 1070 GPU. CONCLUSIONS:The proposed deep learning-based method provides an effective tool for CBCT scatter correction and holds significant value for quantitative imaging and image-guided radiation therapy.

Project description:A human observer study was performed for a signal detection task for the case of fan-beam x-ray computed tomography. Hotelling observer (HO) performance was calculated for the same detection task without the use of efficient channels. By considering the full image covariance produced by the filtered backprojection (FBP) algorithm and avoiding the use of channels in the computation of HO performance, the authors establish an absolute upper bound on signal detectability. Therefore, this study serves as a baseline for relating human and ideal observer performance in the case of fan-beam CT.Eight human observers participated in a two-alternative forced choice experiment where the signal of interest was a small simulated ellipsoid in the presence of independent, identically distributed Gaussian detector noise. Theoretical performance of the HO, which is equivalent to the ideal observer in this case (see Sec. 13.2.12 in Barrett and Myers [Foundations of Image Science (Wiley, Hoboken, NJ, 2004)], was also computed and compared to the performance of the human observers. In addition to a reference FBP implementation, two FBP implementations with inherent loss of HO signal detectability (e.g., by apodizing the ramp filter) were also investigated. Each of these latter two implementations takes the form of a discrete-to-discrete linear operator (i.e., a matrix), which has a nontrivial null-space resulting in the loss of detectability.Estimated observer detectability index (d(A)) values for the human observers and SNR values for the HO were obtained. While Hanning filtering in the FBP implementation with a cutoff frequency of 1/4 of the Nyquist frequency reduces HO SNR (due to the reconstruction matrix's nontrivial null-space), this filtering was shown to consistently improve human observer performance. By contrast, increasing the image pixel size was seen to have a comparable effect on both the HO and the human observers, degrading performance.These results, which characterize HO and human observer performance for a signal detection task in fan-beam FBP noise, form a basis for applying model observer metrics to fan-beam CT when knowledge of the full image-domain noise statistics is important. Further, by calculating HO performance without relying on channels, these results are particularly relevant when an information theoretic approach is considered, e.g., in optimization of the image reconstruction algorithm with respect to preservation of signal detectability. Finally, the HO (which is here equivalent to the ideal observer) provides an absolute upper bound on detection performance, and our results therefore provide insight into the performance of human observers relative to the optimum for two different cases wherein ideal observer performance is compromised through degradation of the data quality. In one case (regularization), human performance is improved to practically ideal performance, and in the other (larger pixel size), ideal and human observer performance are approximately degraded equivalently.

Dataset Information

Scatter correction for full-fan volumetric CT using a stationary beam blocker in a single full scan.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets