Quantitative Comparison of Sinc-Approximating Kernels for Medical Image Interpolation

Abstract. Interpolation is required in many medical image processing operations. From sampling theory, it follows that the ideal interpolation kernel is the sinc function, which is of infinite extent. In the attempt to obtain practical and computationally efficient image processing al-gorithms, many sinc-approximating interpolation kernels have been de-vised. In this paper we present the results of a quantitative comparison of 84 different sinc-approximating kernels, with spatial extents ranging from 2 to 10 grid points in each dimension. The evaluation involves the application of geometrical transformations to medical images from dif-ferent modalities (CT, MR, and PET), using the different kernels. The results show very clearly that, of all kernels with a spatial extent of 2 grid points, the linear interpolation kernel performs best. Of all kernels with an extent of 4 grid points, the cubic convolution kernel is the best (28 %- 75 % reduction of the errors as compared to linear interpolation). Even better results (44 %- 95 % reduction) are obtained with kernels of larger extent, notably the Welch, Cosine, Lanczos, and Kaiser windowed sinc kernels. In general, the truncated sinc kernel is one of the worst performing kernels.

CAVASS: A Computer-Assisted Visualization and Analysis Software System

The Medical Image Processing Group at the University of Pennsylvania has been developing (and distributing with source code) medical image analysis and visualization software systems for a long period of time. Our most recent system, 3DVIEWNIX, was first released in 1993. Since that time, a number of significant advancements have taken place with regard to computer platforms and operating systems, networking capability, the rise of parallel processing standards, and the development of open-source toolkits. The development of CAVASS by our group is the next generation of 3DVIEWNIX. CAVASS will be freely available and open source, and it is integrated with toolkits such as Insight Toolkit and Visualization Toolkit. CAVASS runs on Windows, Unix, Linux, and Mac but shares a single code base. Rather than requiring expensive multiprocessor systems, it seamlessly provides for parallel processing via inexpensive clusters of work stations for more time-consuming algorithms. Most importantly, CAVASS is directed at the visualization, processing, and analysis of 3-dimensional and higher-dimensional medical imagery, so support for digital imaging and communication in medicine data and the efficient implementation of algorithms is given paramount importance

Edge Detection by Adaptive Splitting II. The Three-Dimensional Case

In Llanas and Lantarón, J. Sci. Comput. 46, 485–518 (2011) we proposed an algorithm (EDAS-d) to approximate the jump discontinuity set of functions defined on subsets of ℝ d . This procedure is based on adaptive splitting of the domain of the function guided by the value of an average integral. The above study was limited to the 1D and 2D versions of the algorithm. In this paper we address the three-dimensional problem. We prove an integral inequality (in the case d=3) which constitutes the basis of EDAS-3. We have performed detailed computational experiments demonstrating effective edge detection in 3D function models with different interface topologies. EDAS-1 and EDAS-2 appealing properties are extensible to the 3D cas

Urinary sodium excretion and cardiovascular disease mortality.

Comment on: Stolarz-Skrzypek K, Kuznetsova T, Thijs L, Tikhonoff V, Seidlerová J, Richart T, Jin Y, Olszanecka A, Malyutina S, Casiglia E, Filipovsk J, Kawecka-Jaszcz K, Nikitin Y, Staessen JA; European Project on Genes in Hypertension (EPOGH) Investigators. Fatal and nonfatal outcomes, incidence of hypertension, and blood pressure changes in relation to urinary sodium excretion. JAMA. 2011 May 4;305(17):1777-85. PMID: 21540421

Quantitative comparison of sinc-approximating kernels for medical image interpolation

Interpolation is required in many medical image processing operations. From sampling theory, it follows that the ideal interpolation kernel is the sinc function, which is of infinite extent. In the attempt to obtain practical and computationally efficient image processing algorithms, many sinc-approximating interpolation kernels have been devised. In this paper we present the results of a quantitative comparison of 84 different sinc-approximating kernels, with spatial extents ranging from 2 to 10 grid points in each dimension. The evaluation involves the application of geometrical transformations to medical images from different modalities (CT, MR, and PET), using the different kernels. The results show very clearly that, of all kernels with a spatial extent of 2 grid points, the linear interpolation kernel performs best. Of all kernels with an extent of 4 grid points, the cubic convolution kernel is the best (28% – 75% reduction of the errors as compared to linear interpolation). Even better results (44% – 95% reduction) are obtained with kernels of larger extent, notably the Welch, Cosine, Lanczos, and Kaiser windowed sinc kernels. In general, the truncated sinc kernel is one of the worst performing kernels