function?[xedgenoise]=adp_spa_denoise_new3（trM1）
%基于個(gè)數(shù)判斷的脈沖噪聲的中值濾波
[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