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.