matlab-(含教程)基于matlab的非线性滤波和各向异性扩散

function [ I2 ] = anisotropic_diff_lma( I,lambda,k,typeg ) [r,c]=size(I); I2=zeros(size(I)); GNI=0;GSI=0;GEI=0;GWI=0; for i=2:r-1 for j=2:c-1 %% for diffusion in north direction GNI=I(i-1,j)-I(i,j); [cN] = ConductionCoeff( GNI,k,typeg ); %% for diffusion in south direction GSI=I(i+1,j)-I(i,j); cS = ConductionCoeff( GNI,k,typeg ); %% for diffusion in east direction GEI=I(i,j+1)-I(i,j); cE = ConductionCoeff( GEI,k,typeg ); %% for diffusion in west direction GWI=I(i,j-1)-I(i,j); cW = ConductionCoeff( GWI,k,typeg ); %% anisotropic (discrete) implementation I2(i,j) = I(i,j) + (((cN*GNI)+(cS*GSI)+(cE*GEI)+(cW*GWI))*lambda); if I2(i,j)>=1 I2(i,j)=1; end if I2(i,j)<=0 I2(i,j)=0; end end end end