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.