findmask.m

来自「用小波实现图像去噪的一个程序」· M 代码 · 共 28 行

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%% 寻找掩膜，这里写成一个函数，方便调用%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


function [W,Cor,M,k,pw,NCor]= findmask(W,Cor,M,k,NCor,weight)
   
     %%%%%%%%%%%%%%%%%% 检测掩膜%%%%%%%%%%%%%%%%%%%
     [m,n]=size(W);
      for i=1:1:m
          for j=1:1:n
              if (abs(NCor(i,j))>=weight*abs(W(i,j)))&(M(i,j)==0)
                  Cor(i,j)=0;
                  W(i,j)=0;
                  M(i,j)=1;
                  k=k+1;
              end
          end
      end
      
      %%%%%%%%%%%%%%%%%%顺便计算一下能量，使得判断的时候有用，呵呵！
      pw=sum(sum(W.*W));
      PCor=sum(sum(Cor.*Cor));
      normalization=sqrt(pw/(PCor+eps));
      NCor=Cor*normalization;