lier2_final.m

来自「恒虚警监测的代码」· M 代码 · 共 41 行

%GROUP: Meteor Team
%name:  Zhang Zhi    &&    Sun Xiaoyu
[e-mail]
%deadline 18 nov 2005
clear all,close all;
load lier2.mat
load noise.mat
noise=img(1:75,1:66);
sum=lier2+noise;
angle_diff=0.29;
att=0.4;
filter=min(lier2,noise);
k=zeros(size(filter));
s=angle(filter);
[m,n]=size(filter);
  for i=2:m-1
      for j=2:n-1
          if(abs(s(i,j-1)-s(i,j))<angle_diff)&(abs(s(i-1,j)-s(i,j))<angle_diff)&(abs(s(i,j+1)-s(i,j))<angle_diff)&(abs(s(i+1,j)-s(i,j))<angle_diff)
              k(i,j)=abs(filter(i,j));
          else
              k(i,j)=abs(filter(i,j))*att;
          end;
      end
  end
  for i=2:m-1
      for j=2:n-1
          if(abs(k(i,j-1)-k(i,j))>10)&(abs(k(i-1,j)-k(i,j))>10)&(abs(k(i,j+1)-k(i,j))>10)&(abs(k(i+1,j)-k(i,j))>10)
             k(i,j)=(k(i,j-1)+k(i-1,j)+k(i,j+1)+k(i+1,j))/2;
          end;
      end;
  end;
 result=k;
 endresult=CFARzz(result,16);
 figure(1),
 mesh(endresult);
 figure(2),
 mesh(result);
 figure(3),
 surf(k);
 figure(4),
 surf(endresult);