按 Ctrl+C 复制代码clear all;close all;clc;k=15; %导热系数,控制平滑lambda=0.15; %控制平滑N=20; %迭代次数img=double(imread('lena.jpg'));imshow(img,[]);[m n]=size(img);imgn=zeros(m,n);for i=1:N for p=2:m-1 for q=2:n-1 %当前像素的散度,对四个方向分别求偏导,局部不同方向上的变化量, %如果变化较多,就证明是边界,想方法保留边界 NI=img(p-1,q)-img(p,q); SI=img(p+1,q)-img(p,q); EI=img(p,q-1)-img(p,q); WI=img(p,q+1)-img(p,q); %四个方向上的导热系数,该方向变化越大,求得的值越小,从而达到保留边界的目的 cN=exp(-NI^2/(k*k)); cS=exp(-SI^2/(k*k)); cE=exp(-EI^2/(k*k)); cW=exp(-WI^2/(k*k)); imgn(p,q)=img(p,q)+lambda*(cN*NI+cS*SI+cE*EI+cW*WI); %扩散后的新值 end end img=imgn; %整个图像扩散完毕,用已扩散图像的重新扩散。endfigure;imshow(imgn,[]);clear all;close all;clc;k=15; 你的当前访问异常,请进行认证后继续阅读剩余内容。 提交