X-ray fluoroscopy noise modeling for filter design

Abstract

Fluoroscopy is an invaluable tool in various medical practices such as catheterization or image-guided surgery. Patient's screen for prolonged time requires substantial reduction in X-ray exposure: The limited number of photons generates relevant quantum noise. Denoising is essential to enhance fluoroscopic image quality and can be considerably improved by considering the peculiar noise characteristics. This study presents analytical models of fluoroscopic noise to express the variance of noise as a function of gray level, a practical method to estimate the parameters of the models and a possible application to improve the performance of noise filtering. Quantum noise is modeled as a Poisson distribution and results strongly signal-dependent. However, fluoroscopic devices generally apply gray-level transformations (i.e., logarithmic-mapping, gamma-correction) for image enhancement. The resulting statistical transformations of the noise were analytically derived. In addition, a characterization of the statistics of noise for fluoroscopic image differences was offered by resorting to Skellam distribution. Real fluoroscopic sequences of a simple step-phantom were acquired by a conventional fluoroscopic device and were utilized as actual noise measurements to compare with. An adaptive spatio-temporal filter based on the local conditional average of similar pixels has been proposed. The gray-level differences between the local pixel and the neighboring pixels have been assumed as measure of similarity. Filter performance was evaluated by using real fluoroscopic images of a step phantom and acquired during a pacemaker implantation. The comparison between experimental data and the analytical derivation of the relationship between noise variance and mean pixel intensity (noise-parameter models) were presented relatively to raw-images, after applying logarithmic-mapping or gamma-correction and for difference images. Results have confirmed a great agreement (adjusted R-squared values > 0.8). Clipping effects of real sensors were also addressed. A fine image restoration has been obtained by using a conditioned spatio-temporal average filter based on the noise statistics previously estimated. Fluoroscopic noise modeling is useful to design effective procedures for noise estimation and image filtering. In particular, filter performance analysis has showed that the knowledge of the noise model and the accurate estimate of noise characteristics can significantly improve the image restoration, especially for edge preserving. Fluoroscopic image enhancement can support further X-ray exposure reduction, medical image analysis and automated object identification (i.e., surgery tools, anatomical structures).