 |
|
||||||||||||||||||||||||||||||||||
|
|||||||||||||||||||||||||||||||||||
 |
|||||||||||||||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
 |
|
Dieses Dokument ist leider nur in Englisch verfügbar.
Abstract | Pictures | 3D Convolution | Wavelets | Morphology | Results | Papers Hardware Accelerated FiltersAbstractMany volume analysis operations used for image enhancement, data processing, or feature extraction use three dimensional convolutions, morphological operators, or wavelet analysis. These three categories of filters are usually performed by the main processor, which is usually bound by main memory bandwidth due to the huge data sets nowadays. Graphics memory, on the other hand, has much higher bandwidth and can be accessed in a more parallel manner. Additionally, the cache strategies used in these architectures usually suite much more the way memory has to be addressed for linear filters. By exploiting these hardware systems we are able to accelerate time consuming filtering steps. Pictures
3D ConvolutionFor most interactive 3D visualization techniques the data set has to be loaded to some dedicated texture memory located in the graphics subsystem. In order to be able to interactively filter these data sets, algorithms working directly on the graphics hardware are of additional benefit, because the data set does not have to be reloaded after filtering. Most linear filters can be decomposed into a seperable convolution. Current graphics hardware is able to perform 2D convolutions. In order to extend this filter type to 3D, we have to decompose the 3D convolution into a 2D one perpendicular to the z axis, and into a 1D convolution along the z axis. The convolution itself is performed in the frame buffer while slices of the volume are read back into texture memory. Wavelet AnalysisWavelets are a widely used multiresolution analysis. Figure 2 depicts the simplest wavelet of all, the so-called Haar wavelet. Using this mother wavelet, a data set can be decomposed lossless into a set of wavelet coefficients by performing a wavelet decomposition. The reverse step - recreating the data set from the coefficients - is called reconstruction.
Both operators, decomposition and reconstruction, can benefit greatly from the high memory bandwidth of modern graphics cards. However, the OpenGL pipeline does not completely fit the necessary operator sequence. Thus the invocation sequence of OpenGL commands with their according parameters is - while not time consuming - rather complex. The presented papers discuss in detail these aspects. During a typical decomposition step the data set has to be filtered several times. Due to the lack of floating point data types in the graphics hardware, accuracy is quite an issue. However, wavelets seem to be relatively stable in that aspect. Morphological OperatorsThe disadvantage of wavelet transformation and other filters, that are based on linear combinations of the input data, is the fact that the analyzed structures become distorted. This can be seen especially on lower frequency scales, accomplished by large filter kernels. One class of special non-linear filters that do not flatten the contours of the original data set are the morphological operators. This filter type computes the minimum respectively maximum of pixel values within a given scope. Before the actual maxima are calculated, the values of the so-called structuring element are added to the pixel values. As one can see in Figure~3 for binary operators, the erosion operator cuts away parts of the boundary of the analyzed image X. The amount that is removed is defined by the structuring element Y, which is typically much smaller than the input image. The dilation operator on the other hand enlarges the input data. By combining these basic operators to more complex ones we get the so-called opening and closing operators. The opening operator breaks up small bridges between connected regions while the closing operator tends to fill small gaps in solid components. Morphological operators on gray scale images are defined accordingly.
As usual we assume that the given data set is already stored in texture memory. When filtering a volume data set with a so-called diamond shaped morphological operator, we sweep over the volume in three passes. In the first two passes we perform the filtering along the x and y axes. These steps only access data from within one single texture, which contains data perpendicular to the z axis, as defined in the previous section. In the third pass we filter along the z axis. As the data along this axis is spread in multiple textures, triangles textured with several different images have to be rendered. ResultsBy using standard OpenGL hardware the time needed for typical filtering applications can be significantly reduced, as it is shown in Table 1 to 3.
Papers and Technical Reports
|
|
Letzte Änderung Sep 12 2000 (wwwcgi) | © Universität Stuttgart | Impressum |