matlab 基于噪声点判断的中值滤波法

function [xedgenoise]=adp_spa_denoise_new3（trM1）
%基于个数判断的脉冲噪声的中值滤波
[depwide]=size（t）;
x=ones（size（t））;
edge=zeros（size（t））;
noise=zeros（size（t））;

for i=2:dep-1
    for j=2:wide-1
       f1=t（i-1j-1）-t（ij）;
       f2=t（i-1j）-t（ij）;
       f3=t（i-1j+1）-t（ij）;
      
       f4=t（ij-1）-t（ij）;
       f5=t（ij）-t（ij）;
       f6=t（ij+1）-t（ij）;
       
       f7=t（i+1j-1）-t（ij）;
       f8=t（i+1j）-t（ij）;
       f9=t（i+1j+1）-t（ij）;
       
       f=[f1 f2 f3 f4 f5 f6 f7 f8 f9];
       f_mean=0;
       f_var=0;
       
       for k=1:9
           f_mean=f_mean+f（k）;
       end
       
       f_mean=f_mean/9;
       
       for k=1:9
           f_var=f_var+（f