stable_unstable_diffusion.m

超声图像实时消噪实验程序

a = imread('cameraman.tif');
a = im2double( a );
niter = 70;

smoothkernel = inline('[ 0 dt 0; dt (1-4*dt) dt; 0 dt 0 ]','dt');
subplot(1,3,1); imshow(a); title('Original Image');
b=a;
subplot(1,3,2); 
for i=1:niter
	b = conv2br( b, smoothkernel(0.24) );
end
imshow(b);title('dt = 0.24');

c=a;
subplot(1,3,3);
for i=1:niter
	c = conv2br( c, smoothkernel(0.26) );
end
imshow(c);title('dt = 0.26');