www.gusucode.com > IPCV_Eval_Kit_R2019a_0ce6858工具箱matlab程序源码 > IPCV_Eval_Kit_R2019a_0ce6858/code/demo_files/I2_07_2_2_3DFilter_fspecial3.m
%% MRI夋憸偺3師尦僈僂僔傾儞僼傿儖僞儕儞僌 clear; clc; close all; % 弶婜壔 %% MRI夋憸偺儘乕僪 load mristack % 壜帇壔 intensity = [0 20 40 120 220 1024]; alpha = [0 0 0 0 0.38 0.5]; color = ([0 0 0; 43 0 0; 103 37 20; 199 155 97; 216 213 201; 255 255 255]) ./ 255; queryPoints = linspace(min(intensity),max(intensity),256); alphamap = interp1(intensity,alpha,queryPoints)'; colormap = interp1(intensity,color,queryPoints); figure; volshow(mristack,'Colormap',colormap,... 'Alphamap',alphamap,'ScaleFactors',[1 1 10]); %% 3D懭墌僼傿儖僞乕偵傛傞僗儉乕僕儞僌 H = fspecial3('ellipsoid',[7 7 3]); volSmooth = imfilter(mristack,H,'replicate'); figure; volshow(volSmooth,'Colormap',colormap,... 'Alphamap',alphamap,'ScaleFactors',[1 1 10]); %% 3D僜乕儀儖僄僢僕僼傿儖僞乕偵傛傞僄僢僕専弌 H = fspecial3('sobel','Y'); edgesHor = imfilter(mristack,H,'replicate'); alpha = [0 0 0 0 0 0.5]; queryPoints = linspace(min(intensity),max(intensity),256); alphamap = interp1(intensity,alpha,queryPoints)'; figure; volshow(edgesHor,'Colormap',colormap,... 'Alphamap',alphamap,'ScaleFactors',[1 1 10]); %% Copyright 2018 The MathWorks, Inc.