Usually the goal of volume visualization is to get meaningful, twodimensional pictures from threedimensional data sets (for instance from CT or MRI imaging). A useful concept for that purpose is the gradient, which can be interpreted as a normal to a surface of equal values (an isosurface) passing through the point of interest. Consequently, the gradient can be used to shade surfaces, since most lighting models (for instance Gouraud shading [5] or Phong shading [16]) make use of it or it can be used for classification in ray casting techniques [7]. Möller et al. even concede gradient reconstruction having a greater impact on image quality than function reconstruction itself [12].
From signal processing theory we know that the ideal reconstruction
filter is the sinc function [15] given by
Sampling and reconstruction are often investigated in frequency domain. A function can be transformed by means of the Fourier transformation from the spatial domain to the frequency domain, where additional information about its behavior can be obtained (reconstruction filters, for instance, can be compared to the ideal one, i.e., the sinc filter). An excellent introduction to this topic is given by Blinn [2,3]. In order to reconstruct the derivative of the function it is possible to reconstruct the function and derive it discretely. However, we chose another approach which is to derive the ideal reconstruction filter and use it as a derivative reconstruction filter directly.
Deriving the sinc function (Eq. 1) yields the cosc
function [1] (note that the definition is not
analogous):
In this paper we will compare various filters obtained by bounding the ideal derivative reconstruction filter given by Eq. 3 to some finite extend. The simplest way to do this is to truncate the cosc function, which is, anyway, not the best way because of discontinuities at the edges. So special functions (adversely called windows) were designed to get a smoother transition to zero. Their definitions and a discussion of their frequency responses is given in Section 3 after a brief overview of some work already done in this area in Section 2. In Section 4 the test scenario used to compare the windowed cosc filters is described and the result are given and discussed in Section 5.
