Masked smoothing using separable kernels for CT perfusion images

CT perfusion images have a high contrast ratio between voxels representing different anatomy, such as tissue or vessels, which makes image segmentation of tissue and vascular regions relatively easy. However, grey and white matter tissue regions have relatively low values and can suffer from poor signal to noise ratios. While smoothing can improve the image quality of the tissue regions, the inclusion of much higher valued vascular voxels can skew the tissue values. It is thus desirable to smooth tissue voxels separately from other voxel types, as has been previously implemented using mean filter kernels. We created a novel Masked Smoothing method that performs Gaussian smoothing restricted to tissue voxels. Unlike previous methods, it is implemented as a combination of separable kernels and is therefore fast enough to consider for clinical work, even for large kernel sizes.

We compare our Masked Smoothing method to alternatives using Gaussian smoothing on an unaltered image volume and Gaussian smoothing on an image volume with vascular voxels set to zero. Each method was tested on simulation data, collected phantom data, and CT perfusion data sets. We then examined tissue voxels for bias and noise reduction.

Simulation and phantom experiments demonstrate that Masked Smoothing does not bias the underlying tissue value, whereas the other smoothing methods create significant bias. Furthermore, using actual CT perfusion data, we demonstrate significant differences in the calculated CBF and CBV values dependent on the smoothing method used.

The Masked Smoothing is fast enough to allow eventual clinical usage and can remove the bias of tissue voxel values that neighbor blood vessels. Conversely, the other Gaussian smoothing methods introduced significant bias to the tissue voxels.