im1_path = which('lamp.jpg'); im1 = imread(im1_path); figure(1) colormap(gray(256)) clf image(im1) title('Original Image') y = nldifc(im1,7,1,20,linspace(20,100,10),10,1,2,'aos', 'grad', 'dfstep', 5, 'alt1','imscale', 'norm', 'inf'); figure(1) subplot(1,2,1) image(im1) title('Original Image') subplot(1,2,2) image(im1) title('Nonlinear diffusion') figure(2) colormap(hot) y = nldifc(im1,7,1,20,linspace(20,100,10),10,1,2,'aos', 'grad', 'dfstep', 5, 'alt1','imscale', 'norm', 'inf');
http://www.mathworks.com/matlabcentral/fileexchange/3710-nonlinear-diffusion-toolbox