www.pudn.com > mulscal_edge.rar > ada_thr.asv

%对矩阵取自适应阈值 
function out=ada_thr(c,n) 
 
[row,col]=size(c); 
 
m1=(n+1)/2;%n为奇数是模板的大小 
m2=(n-1)/2; 
temp_c=c; 
for irow=m1:(row-m1) 
    for jcol=m1:(col-m1) 
        if c(irow,jcol)>0 
           lefttop_h=irow-m2;%左上角横坐标 
           lefttop_v=jcol-m2;%左上角纵坐标 
           sub_c=sub_m(c,lefttop_h,lefttop_v,n); 
                 mean_array=mean(sub_c); 
              adapt_thr=mean(mean_array)*a0; 
              temp_c(irow,jcol)=wthresh(c(irow,jcol),'h',adapt_thr); 
           end 
        end 
    end 
end 
out=temp_c;