www.pudn.com > mulscal_edge.rar > local_max_dujiaoxian.m

function [out]=local_max_dujiaoxian(cd,n) 
 
[row,col]=size(cd); 
m1=(n+1)/2; 
m2=(n-1)/2; 
out(row,col)=0; 
 
for irow=m1:(row-m1) 
    for jcol=m1:(col-m1) 
        lefttop_h=irow-m2;%×óÉϽǺá×ø±ê 
        lefttop_v=jcol-m2;%×óÉϽÇ×Ý×ø±ê 
        sub_c=sub_m(cd,lefttop_h,lefttop_v,n); 
        if max_c(sub_c,cd(irow,jcol)) 
           out(irow,jcol)=cd(irow,jcol); 
        end 
    end 
end