3 years ago

Tissue classification based on 3D local intensity ... - IEEE Xplore

Tissue classification based on 3D local intensity ... - IEEE Xplore


SATO ET AL.: TISSUE CLASSIFICATION BASED ON 3D LOCAL INTENSITY STRUCTURES FOR VOLUME RENDERING 165 Fig. 3. Plots of normalized p responses of local structure filters for corresponding local models, S fh …0; r †; f g, where r is continuously varied and f ˆ i s i1 ( 1 ˆ 1, s ˆ 2 , and i ˆ 1; 2; 3; 4). See Appendix A for the theoretical derivations of the response curves shown in (b)-(d). (a) Response of the edge filter for the edge model ( ˆ edge). (b) Response of the line filter for the line model ( ˆ line). (c) Response of the blob filter for the blob model ( ˆ blob). (d) Response of the sheet filter for the sheet model ( ˆ sheet). where 2 2fsheet; line; blobg. S int can also be extended so as to be combined with Gaussian blur, in which case we represent it as S int ff; f g. 3.3 Implementation Our 3D local structure filtering methods described above assume that volume data with isotropic voxels are used as input data. However, voxels in medical volume data are usually anisotropic since they generally have lower resolution along the third directionÐi.e., the direction orthogonal to the slice planeÐthan within slices. Rotational invariant feature extraction becomes more intuitive in a space where the sample distances are uniform. That is, structures of a particular size can be detected on the same scale independent of the direction when the signal sampling is isotropic. We therefore introduce a preprocessing procedure for 3D local structure filtering in which we perform interpolation to make each voxel isotropic. Linear and spline-ong>basedong> interpolation methods are often used, but blurring is inherently involved in these approaches. Because, as noted above, the original volume data is inherently blurrier in the third direction, further degradation of the data in that direction should be avoided. For this reason, we opted to employ sinc interpolation so as not to introduce any additional blurring. After Gaussian-shaped slopes are added at the beginning and end of each profile in the third direction to avoid unwanted Gibbs ringing (see Appendix B), sinc interpolation is performed by zero-filled expansion in the frequency domain [22], [23]. The sinc interpolation and 3D local structure filtering were implemented on a Sun Enterprise server with multi- CPUs using multithreaded programming. With eight CPUs (168 MHz), interpolation and single-scale 3D filtering for a 256 256 128 volume were performed in about five minutes. Separable implementation was used for efficient 3D Gaussian derivative convolution. For example, the computation of the second derivatives of Gaussian in the Hessian matrix was implemented using three separate convolutions with one-dimensional kernels as represented by f x i y j z k…x; @ 2 f†ˆ @x i @y j @z G…x; f† f…x† k ˆ di dx i G…x; d j f† dy j G…y; d k f† dz k G…z; f†f…x† ; …20† where i, j, and k are nonnegative integers satisfying i ‡ j ‡ k ˆ 2 and 4 f was used as the radius of the kernel [10]. Using this decomposition, the amount of computation needed can be reduced from O…n 3 † to O…3n†, where n is the kernel diameter. Details of the implementation of eigenvalue computation for the Hessian matrix are described in [10].

166 IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS, VOL. 6, NO. 2, APRIL-JUNE 2000 Fig. 4. Multichannel tissue ong>classificationong>. (a) Example using the max (OR) operations. Discrete ong>classificationong> …p† (upper left) is obtained by combining …p i † in different channels using the max (OR) operation. The resultant multidimensional opacity function …p† (right) is obtained by the min of the discrete ong>classificationong> and a continuous one-dimensional (1D) opacity function 0 …S int † (lower left). (b) Example using the min (AND) 4 MULTICHANNEL TISSUE CLASSIFICATION BASED ON LOCAL STRUCTURES 4.1 Representation of Classification Strategy A multidimensional feature space can be defined whose axes correspond to the measures of similarity to different local structures selected from S 0;1 , S 2 , and M 2 , where 0;1 2fint; edgeg and 2 2fsheet; line; blobg. Let p ˆ …p 1 ;p 2 ; ...;p m † (m is the number of features) be a feature vector in the multidimensional feature space, which consists of multichannel values of the measures of similarity to different local structures at each voxel. The opacity function …p† and color function c…p† ˆ…R…p†;G…p†;B…p†† for volume rendering are given as functions of the multidimensional variable p. We consider the following problem: Given tissue classes t ˆ…t 1 ;t 2 ; ...;t n † (n is the number of tissue classes), how can we determine the tissue opacity function j …p† and tissue color c j …p† for tissue t j ? We use a 1D weight function …p i †, whose range is ‰0; 1Š, of feature p i as a primitive function to specify the opacity function j …p† for tissue t j . We consider a generalized combination rule of the 1D weight functions to obtain the opacity function …p† formally given by …p† ˆmax min 1;i…p i †; min 2;i…p i †; ...; min ri…p i † ; 1im 1im 1im …21† where r is the number of terms min 1im ki …p i † and ki …p i † denotes the weight function for feature p i in the kth term. The max and min operations in (21) respectively correspond to union and intersection operations in fuzzy theory [24]. Suppose ki …p i † has a box-shaped function given by where ki …p i †ˆ…p i ; L ki ;H ki †; …x; L; H† ˆ 1; L x

An Ensemble Classification-Based Approach Applied ... - IEEE Xplore
Neural network based classification of road ... - IEEE Xplore
The Occlusion Spectrum for Volume Classification and ... - IEEE Xplore
Social relationship classification based on interaction ... - IEEE Xplore
Feature-based Modulation Classification Using ... - IEEE Xplore
Tissue Characterization of Equine Tendons With ... - IEEE Xplore
Gabor Volume based Local Binary Pattern for Face ... - IEEE Xplore
Local Intensity Order Pattern for Feature Description - IEEE Xplore
Locally adaptive metric nearest-neighbor classification ... - IEEE Xplore
High-Resolution Determination of Soft Tissue ... - IEEE Xplore
Graphics Cluster Based Visualization of 3D Medical ... - IEEE Xplore
Tissue and Species Differences in the Application of ... - IEEE Xplore
Decision Fusion for the Classification of Urban ... - IEEE Xplore
An Adaptive Approach for Image Classification - IEEE Xplore
Automatic classification of shoeprints for use in ... - IEEE Xplore
A Statistical Image-Based Approach for the 3D ... - IEEE Xplore
A classification and comparison framework for ... - IEEE Xplore
Vision-Based Localization for Leader–Follower ... - IEEE Xplore
Vision Based Vehicle Localization for Autonomous ... - IEEE Xplore
Graph-Based Local Kernel Regression for Image Editing - IEEE Xplore
Visualization by Proxy: A Novel Framework for ... - IEEE Xplore
A Novel, Fast, and Complete 3D Segmentation of ... - IEEE Xplore
Multiscale Poisson Intensity and Density Estimation - IEEE Xplore
Object-Based Shadow Extraction and Correction of ... - IEEE Xplore
A Robust Descriptor based on Weber's Law - IEEE Xplore
A Probabilistic Shape-based Segmentation Using ... - IEEE Xplore
Volume Conduction in an Anatomically Based Surface ... - IEEE Xplore
Modeling of Skin Tissue Ablation by Nanosecond ... - IEEE Xplore
Estimation of surface normals for use in 3D rendering ... - IEEE Xplore