Local filtration based scatter correction for cone-beam CT using primary modulation

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 practical and robust method for beam blocker-based cone beam CT scatter correction

Objective. In the traditional beam-blocker based cone beam CT (CBCT) scatter correction, the scatter measured in the region shaded by lead strips was multiplied by a correction factor to directly represent the scatter in the unblocked region. The correction factor optimization is a tedious process and lacks an objective stop criterion. To skip the optimization process, an indirect scatter estimation method was developed and validated in phantom imaging.Approach.A beam-blocker made of lead strips was mounted between the x-ray source and object for scatter estimation. The primary signal between lead strips in the blocked region was first calculated by subtracting the measured scatter, and then used to calculate the scatter signal in the unblocked region corresponding to the same attenuation path. The calculated scatter signal was smoothed via local filtration and used to correct the measured projection in the unblocked region. Finally, the CBCT was reconstructed via Feldkamp-Davis-Kress algorithm. A Catphan and a head phantom were used to verify the performance of the proposed method in both full- and half-blocker scenarios, and with and without a bow-tie filter.Main Results. For scans without the bow-tie filter, the CT number error was reduced to 3.97±2.27 and 5.51±3.90 HU in the full- and half-blocker scenarios, respectively, for the Catphan, and to 4.01±2.18 and 7.97 ± 4.05 HU for the head phantom. When the bow-tie filter was applied, the CT number error was reduced to 2.29±1.42 and 6.72±0.77 HU in the full- and half-blocker scenarios, respectively, for the Catphan, and 2.35±1.25 and 4.96 ± 1.89 HU for the head phantom.Significance. The proposed method effectively avoids the influence of the inserted beam blocker itself on the scatter intensity estimation, and proves a more practical and robust way for the beam-blocker based scatter correction in CBCT scanning.

Development and validation of a scatter-corrected CBCT image-guided method for cervical cancer brachytherapy

Background and purpose

Multiple patient transfers have a nonnegligible impact on the accuracy of dose delivery for cervical cancer brachytherapy. We consider using on-site cone-beam CT (CBCT) to resolve this problem. However, CBCT clinical applications are limited due to inadequate image quality. This paper implements a scatter correction method using planning CT (pCT) prior to obtaining high-quality CBCT images and evaluates the dose calculation accuracy of CBCT-guided brachytherapy for cervical cancer.

Materials and methods

The CBCT of a self-developed female pelvis phantom and five patients was first corrected using empirical uniform scatter correction in the projection domain and further corrected in the image domain. In both phantom and patient studies, the CBCT image quality before and after scatter correction was evaluated with registered pCT (rCT). Model-based dose calculation was performed using the commercial package Acuros®BV. The dose distributions of rCT-based plans and corrected CBCT-based plans in the phantom and patients were compared using 3D local gamma analysis. A statistical analysis of the differences in dosimetric parameters of five patients was also performed.

Results

In both phantom and patient studies, the HU error of selected ROIs was reduced to less than 15 HU. Using the dose distribution of the rCT-based plan as the baseline, the γ pass rate (2%, 2 mm) of the corrected CBCT-based plan in phantom and patients all exceeded 98% and 93%, respectively, with the threshold dose set to 3, 6, 9, and 12 Gy. The average percentage deviation (APD) of D90 of HRCTV and D2cc of OARs was less than 1% between rCT-based and corrected CBCT-based plans.

Conclusion

Scatter correction using a pCT prior can effectively improve the CBCT image quality and CBCT-based cervical brachytherapy dose calculation accuracy, indicating promising prospects in both simplified brachytherapy processes and accurate brachytherapy dose delivery.

Planning CT-guided robust and fast cone-beam CT scatter correction using a local filtration technique

PURPOSE

Cone-beam CT (CBCT) has been widely utilized in image-guided radiotherapy. Planning CT (pCT)-aided CBCT scatter correction could further enhance image quality and extend CBCT application to dose calculation and adaptive planning. Nevertheless, existing pCT-based approaches demand accurate registration between pCT and CBCT, leading to limited imaging performance and increased computational cost when large anatomical discrepancies exist. In this work, we proposed a robust and fast CBCT scatter correction method using local filtration technique and rigid registration between pCT and CBCT (LF-RR).

METHODS

First of all, the pCT was rigidly registered with CBCT, then forward projection was performed on registered pCT to create scatter-free projections. The raw scatter signals were obtained via subtracting the scatter-free projections from the measured CBCT projections. Based on frequency and intensity threshold criteria, reliable scatter signals were selected from the raw scatter signals, and further filtered for global scatter estimation via local filtration technique. Finally, corrected CBCT was reconstructed with the projections generated by subtracting the scatter estimation from the raw CBCT projections using FDK algorithm. The LF-RR method was evaluated via comparison with another pCT-based scatter correction method based on Median and Gaussian filters (MG method).

RESULTS

Proposed method was first validated on an anthropomorphic pelvis phantom, and showed satisfied performance on scatter removal even when anatomical mismatches were intentionally created on pCT. The quantitative analysis was further performed on four clinical CBCT images. Compared with the uncorrected CBCT, CBCT corrected by MG with rigid registration (MG-RR), MG with deformable registration (MG-DR), and LF-RR reduced the CT number error from 79 ± 35 to 25 ± 18 , 17 ± 13 and 7 ± 3 HU for adipose and from 115 ± 61 to 36 ± 22 , 30 ± 24 , 7 ± 3 HU for muscle, respectively. After correction, the spatial non-uniformity (SNU) of CBCT corrected with MG-RR, MG-DR and LF-RR was 51 ± 13 , 60 ± 21 , and 21 ± 9 HU for adipose, and 50 ± 22 , 57 ± 41 , and 25 ± 6 HU for muscle, respectively. Meanwhile, the contrast-to-noise ratio (CNR) between muscle and adipose was increased by a factor of 2.70, 2.89 and 2.56, respectively. Using the LF-RR method, the scatter correction of 656 projections can be finished within 10 s and the corrected volumetric images (200 slices) can be obtained within 2 min.

CONCLUSION

We developed a fast and robust pCT-based CBCT scatter correction method which exploits the local-filtration technique to promote the accuracy of scatter estimation and is resistant to pCT-to-CBCT registration uncertainties. Both phantom and patient studies showed the superiority of the proposed correction in imaging accuracy and computational efficiency, indicating promisingfuture clinical application.

Impact of bowtie filter and detector collimation on multislice CT scatter profiles: A simulation study

PURPOSE

To investigate via Monte Carlo simulations, the impact of scan subject size, antiscatter grid (ASG), collimator size, and bowtie filter on the distribution of scatter radiation in a typical realistically modeled third generation 16 slice diagnostic computed tomography (CT) scanner.

METHODS

Full radiation transport was simulated with Geant4 in a realistic CT scanner geometric model, including the imaging phantom, bowtie filter (BTF), collimators and detector assembly, except for the ASGs. An analytical method was employed to quantify the probable transmission through the ASG of each photon intersecting the detector array. Normalized scatter profiles (NSP) and scatter-to-primary-ratio (SPR) profiles were simulated for 90 and 140 kVp beams for different size phantoms and slice thicknesses. The impact of CT scatter on the reconstructed attenuation coefficient factor was also studied as were the modulating effects of phantom- and patient-tissue heterogeneities on scatter profiles. A method to characterize the relative spatial frequency content of sinogram signals was developed to assess the latter.

RESULTS

For the 21.4-cm diameter phantom, NSP and SPR increase linearly with collimator opening for both tube potentials, with the 90 kVp scan exhibiting slightly larger NSP and SPR. The BTF modestly modulates scatter under the phantom center, reducing the prominent off-axis lobes by factors of 1.1-1.3. The ASG reduces scatter on the central axis NSP threefold, and reduces scatter at the detectors outside the phantom shadow by factors of 25 to 500. For the phantoms with diameters of 27 and 32 cm, the scatter increases roughly three- and fourfold, respectively, demonstrating that scatter monotonically increases with phantom size, despite deployment of the ASG and BTF. In the absence of a scan subject, the ASG reduces the signal profile arising photons scattered by the BTF. Without ASG, the in-air scatter profile is relatively flat compared to the scatter profile when the ASG is present. For both 90 and 140 kVp photon spectra, the calculated attenuation coefficient decreases linearly with increasing collimation size. For both homogeneous and heterogeneous objects, NSPs are dominated by low spatial frequency content compared to the primary signal. However, the SPR, which quantifies the local magnitude of nonlinear detector response and is dominated by the high frequency content of the primary profile, can contribute strongly to high-spatial frequency streaking artifacts near high-density structures in reconstructed image artifacts.

CONCLUSION

Public-domain Monte Carlo codes, Geant-4 in particular, is a feasible method for characterizing CT detector response to scattered- and off-focal radiation. Our study demonstrates that the ASG substantially reduces the scatter radiation and reshapes scatter-radiation profiles and affects the accuracy with which the detector array can measure narrow-beam attenuation due its inability to distinguish between true uncollided primary and narrow-angle coherently scattered photons. Hence, incorporating the impact of detector array collimation into the forward-projection signal formation models used by iterative reconstruction algorithms is necessary to use CT for accurately characterizing material properties. While tissue heterogeneities exercise a modest influence on local NPS shape and magnitude, they do not add significant high spatial frequency content.

Physics-based spectral compensation algorithm for x-ray CT with primary modulator

X-ray computed tomography (CT) scatter correction using primary modulator has been continuously developed over the past years, with progress in improving the performance of scatter correction. In this work, we further advance the primary modulator technique towards practical applications where the spectral nonuniformity caused by the modulator continues to be a challenging problem. A physics-based spectral compensation algorithm is proposed to adaptively correct for the spectral nonuniformity, and hence to reduce the resultant ring artifacts on reconstructed CT images. First, an initial spectrum of the CT system without the primary modulator is modeled using an understanding of x-ray CT physics, and optimized by an expectation maximization method; then, the optimized estimation of the initial spectrum is utilized to adaptively calculate the effective modulator thickness from measured transmissions of the primary modulator at each detector element, leading to a set of new spectra that is able to capture the nonuniform spectral distribution of the primary modulator; finally, using the modulator-modeled spectrum, a beam hardening mapping function is generated and beam hardening correction is applied to CT projections. A CatPhan600 phantom and an anthropomorphic thorax phantom were scanned with three different primary modulators to evaluate the approach. For the Catphan phantom, the spectral compensation algorithm efficiently removes the ring (and band) artifacts that otherwise dominate the reconstructed CT image. For the three modulators with nominal copper thickness of 52.5, 105 and 210 [Formula: see text]m, our method reduces the CT number nonuniformity from 147.9, 436.2 and 696.4 Hounsfield units (HU) to 14.6, 26.2 and 13.6 HU, respectively, close to that of the reference image (i.e. 7.5 HU). For the thorax phantom, the ring artifacts are also suppressed significantly on the transaxial image; on the sagittal image, the alternating black-and-white patterns are largely removed, with the CT number nonuniformity being reduced from 282.0 HU to 38.5 HU.

Fast shading correction for cone-beam CT via partitioned tissue classification

The quantitative use of cone beam computed tomography (CBCT) in radiation therapy is limited by severe shading artifacts, even with system embedded correction. We recently proposed effective shading correction methods, using planning CT (pCT) as prior information to estimate low-frequency errors in either the projection domain or image domain. In this work, we further improve the clinical practicality of our previous methods by removing the requirement of prior pCT images. Clinical CBCT images are typically composed of a limited number of tissues. By utilizing the low frequency characteristic of shading distribution, we first generate a 'shading-free' template image by enforcing uniformity on CBCT voxels of the same tissue type via a technique named partitioned tissue classification. Only a small subset of voxels in the template image are used in the correction process to generate sparse samples of shading artifacts. Local filtration, a Fourier transform based algorithm, is employed to efficiently process the sparse errors to compute a full-field distribution of shading artifacts for CBCT correction. We evaluate the method's performance using an anthropomorphic pelvis phantom and 6 pelvis patients. The proposed method improves the image quality of CBCT for both phantom and patients to a level matching that of pCT. On the pelvis phantom, the signal non-uniformity (SNU) is reduced from 12.11% to 3.11% and 8.40% to 2.21% on fat and muscle, respectively. The maximum CT number error is reduced from 70 to 10 HU and 73 to 11 HU on fat and muscle, respectively. On patients, the average SNU is reduced from 9.22% to 1.06% and 11.41% to 1.67% on fat and muscle, respectively. The maximum CT number error is reduced from 95 to 9 HU and 88 to 8 HU on fat and muscle, respectively. The typical processing time for one CBCT dataset is about 45 s on a standard PC.

Evaluation and Clinical Application of a Commercially Available Iterative Reconstruction Algorithm for CBCT-Based IGRT

PURPOSE

We have quantitatively evaluated the image quality of a new commercially available iterative cone-beam computed tomography reconstruction algorithm over standard cone-beam computed tomography image reconstruction results.

METHODS

This iterative cone-beam computed tomography reconstruction pipeline uses a finite element solver (AcurosCTS)-based scatter correction and a statistical (iterative) reconstruction in addition to a standard kernel-based correction followed by filtered back-projection-based Feldkamp-Davis-Kress cone-beam computed tomography reconstruction. Standard full-fan half-rotation Head, half-fan full-rotation Head, and standard Pelvis cone-beam computed tomography protocols have been investigated to scan a quality assurance phantom via the following image quality metrics: uniformity, HU constancy, spatial resolution, low contrast detection, noise level, and contrast-to-noise ratio. An anthropomorphic head phantom was scanned for verification of noise reduction. Clinical patient image data sets for 5 head/neck patients and 5 prostate patients were qualitatively evaluated.

RESULTS

Quality assurance phantom study results showed that relative to filtered back-projection-based cone-beam computed tomography, noise was reduced from 28.8 ± 0.3 HU to a range between 18.3 ± 0.2 and 5.9 ± 0.2 HU for Full-Fan Head scans, from 14.4 ± 0.2 HU to a range between 12.8 ± 0.3 and 5.2 ± 0.3 HU for Half-Fan Head scans, and from 6.2 ± 0.1 HU to a range between 3.8 ± 0.1 and 2.0 ± 0.2 HU for Pelvis scans, with the iterative cone-beam computed tomography algorithm. Spatial resolution was marginally improved while results for uniformity and HU constancy were similar. For the head phantom study, noise was reduced from 43.6 HU to a range between 24.8 and 13.0 HU for a Full-Fan Head and from 35.1 HU to a range between 22.9 and 14.0 HU for a Half-Fan Head scan. The patient data study showed that artifacts due to photon starvation and streak artifacts were all reduced, and image noise in specified target regions were reduced to 62% ± 15% for 10 patients.

CONCLUSION

Noise and contrast-to-noise ratio image quality characteristics were significantly improved using the iterative cone-beam computed tomography reconstruction algorithm relative to the filtered back-projection-based cone-beam computed tomography method. These improvements will enhance the accuracy of cone-beam computed tomography-based image-guided applications.

Single-Scan Dual-Energy CT Using Primary Modulation

Compared with conventional computed tomography (CT), dual-energy CT (DECT) provides better material differentiation but requires projection data acquired with two different effective x-ray spectra, limiting DECT applications to specialized scanners. We propose a hardware-based method, known as PM-DECT, which utilizes primary beam modulation to enable single-scan DECT on a conventional CT scanner. PM-DECT inserts an attenuation sheet with a spatially varying pattern-primary beam modulator-between the x-ray source and imaged object. During a CT scan, the modulator selectively hardens the x-ray beam, thereby increasing the average photon energy at specific detector pixel locations. Thus, PM-DECT simultaneously acquires high and low energy data at each projection angle. From the sparse projection data, high and low energy CT images are jointly reconstructed and simultaneously decomposed into basis materials via an iterative CT reconstruction algorithm with gradient weighting and an improved version of similarity based regularization. Studies on Catphan 600 and anthropomorphic head phantoms demonstrate that PM-DECT retains a high level of spatial resolution compared with conventional CT scans. Electron density values calculated from decomposed images indicate a limited error of 1.12% for PM-DECT. Comparison against a two-scan DECT technique shows that PM-DECT's image reconstruction from sparse data sets contributes only 0.66% error. By granting the opportunity for high-quality single-scan DECT on conventional CT scanners via limited hardware modification, PM-DECT has the potential to liberate DECT from specialized scanners, extending clinical availability, and implementation.

Optimization of the geometry and speed of a moving blocker system for cone-beam computed tomography scatter correction

PURPOSE

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.

METHODS

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).

RESULTS

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.

CONCLUSIONS

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.

4D cone-beam computed tomography (CBCT) using a moving blocker for simultaneous radiation dose reduction and scatter correction

Four-dimensional (4D) x-ray cone-beam computed tomography (CBCT) is important for a precise radiation therapy for lung cancer. Due to the repeated use and 4D acquisition over a course of radiotherapy, the radiation dose becomes a concern. Meanwhile, the scatter contamination in CBCT deteriorates image quality for treatment tasks. In this work, we propose the use of a moving blocker (MB) during the 4D CBCT acquisition ('4D MB') and to combine motion-compensated reconstruction to address these two issues simultaneously. In 4D MB CBCT, the moving blocker reduces the x-ray flux passing through the patient and collects the scatter information in the blocked region at the same time. The scatter signal is estimated from the blocked region for correction. Even though the number of projection views and projection data in each view are not complete for conventional reconstruction, 4D reconstruction with a total-variation (TV) constraint and a motion-compensated temporal constraint can utilize both spatial gradient sparsity and temporal correlations among different phases to overcome the missing data problem. The feasibility simulation studies using the 4D NCAT phantom showed that 4D MB with motion-compensated reconstruction with 1/3 imaging dose reduction could produce satisfactory images and achieve 37% improvement on structural similarity (SSIM) index and 55% improvement on root mean square error (RMSE), compared to 4D reconstruction at the regular imaging dose without scatter correction. For the same 4D MB data, 4D reconstruction outperformed 3D TV reconstruction by 28% on SSIM and 34% on RMSE. A study of synthetic patient data also demonstrated the potential of 4D MB to reduce the radiation dose by 1/3 without compromising the image quality. This work paves the way for more comprehensive studies to investigate the dose reduction limit offered by this novel 4D MB method using physical phantom experiments and real patient data based on clinical relevant metrics.

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.

A Feasibility Study of Low-Dose Single-Scan Dual-Energy Cone-Beam CT in Many-View Under-Sampling Framework

A single-scan dual-energy low-dose cone-beam CT (CBCT) imaging technique that exploits a multi-slit filter is proposed in this paper. The multi-slit filter installed between the x-ray source and the scanned object is reciprocated during a scan. The x-ray beams through the slits would generate relatively low-energy x-ray projection data, while the filtered beams would make high-energy projection data. An iterative image reconstruction algorithm that uses an adaptive-steepest-descent method to minimize image total-variation under the constraint of data fidelity was applied to reconstructing the image from the low-energy projection data. Since the high-energy projection data suffer from a substantially high noise level due to the beam filtration, we have developed a new algorithm that exploits the joint sparsity between the low- and high-energy CT images for image reconstruction of the high-energy CT image. The proposed image reconstruction algorithm uses a gradient magnitude image (GMI) of the low-energy CT image by regularizing the difference of GMIs of the low- and high-energy CT images to be minimized. The feasibility of the proposed technique has been demonstrated by the use of various phantoms in the experimental CBCT setup. Furthermore, based on the proposed dual-energy imaging, a material differentiation was performed and its potential utility has been shown. The proposed imaging technique produced promising results for its potential application to a low-dose single-scan dual-energy CBCT.

Fast shading correction for cone beam CT in radiation therapy via sparse sampling on planning CT

PURPOSE

The image quality of cone beam computed tomography (CBCT) is limited by severe shading artifacts, hindering its quantitative applications in radiation therapy. In this work, we propose an image-domain shading correction method using planning CT (pCT) as prior information which is highly adaptive to clinical environment.

METHOD

We propose to perform shading correction via sparse sampling on pCT. The method starts with a coarse mapping between the first-pass CBCT images obtained from the Varian TrueBeam system and the pCT. The scatter correction method embedded in the Varian commercial software removes some image errors but the CBCT images still contain severe shading artifacts. The difference images between the mapped pCT and the CBCT are considered as shading errors, but only sparse shading samples are selected for correction using empirical constraints to avoid carrying over false information from pCT. A Fourier-Transform-based technique, referred to as local filtration, is proposed to efficiently process the sparse data for effective shading correction. The performance of the proposed method is evaluated on one anthropomorphic pelvis phantom and 17 patients, who were scheduled for radiation therapy. (The codes of the proposed method and sample data can be downloaded from https://sites.google.com/view/linxicbct) RESULTS: The proposed shading correction substantially improves the CBCT image quality on both the phantom and the patients to a level close to that of the pCT images. On the phantom, the spatial nonuniformity (SNU) difference between CBCT and pCT is reduced from 74 to 1 HU. The root of mean square difference of SNU between CBCT and pCT is reduced from 83 to 10 HU on the pelvis patients, and from 101 to 12 HU on the thorax patients. The robustness of the proposed shading correction is fully investigated with simulated registration errors between CBCT and pCT on the phantom and mis-registration on patients. The sparse sampling scheme of our method successfully avoids false structures in the corrected CBCT even when the maximum registration error is as high as 8 mm.

CONCLUSION

We develop an effective shading correction algorithm for CBCT readily implementable on clinical data as a software plug-in without modifications of current imaging hardware and protocol. The algorithm is directly applied on the output images from a commercial CBCT scanner with high computational efficiency and negligible memory burden.

Virtual scatter modulation for X-ray CT scatter correction using primary modulator

A new scatter estimation algorithm with a concept of virtual scatter modulation for X-ray scatter correction using primary modulator is proposed to reduce the aliasing errors in the estimated scatter. Virtual scatter modulation can be realized through dividing the measured primary-modulated image by the measured modulation function. After the division, the aggravation of the aliasing of primary due to the non-uniformity of the modulation function is largely transferred to that of scatter. Since scatter in general has less high frequencies than primary does, the aggravation of its aliasing is expected to be weaker, and therefore the overall aliasing can be reduced. A CatPhan©600 phantom and an anthropomorphic thorax phantom are scanned on a tabletop X-ray cone-beam computed tomography system to validate our proposed algorithm. On the Catphan phantom, the oscillations that are clearly observed in the central region of the Catphan scatter profile estimated using the original primary-modulation algorithm, are mostly eliminated with the proposed scatter modulation algorithm, leading to less residual artifacts and better CT number uniformity in the reconstructed image. Compared with 38.9 HU of CT nonuniformity in a selected uniform region when the primary-modulation algorithm is used, the new algorithm significantly reduces it to 4.5 HU, reaching the same level of uniformity as the ground truth reference. On the thorax phantom, overall better CT number uniformity is also achieved.