www.pudn.com > ridge_extract.rar > coordinate.m, change:2007-05-08,size:399b

function a=coordinate(centerx,centery,I) 
M=double(I); 
[r,c]=size(M); 
[n,nc]=size(centerx); 
% for i=1:r 
%     for j=1:c 
%        for k=1:n 
%         radiu =sqrt((i-centerx(k))^2+(j-centery(k))^2); 
%         if  (radiu=10)&(radiu>=6) 
%             M(i,j)=0; 
%         end 
%        end 
% end 
% end  
for i=1:3 
    M(centerx+i-1,:)=0; 
end 
for j=1:3 
    M(:,centery+j-1)=0; 
end 
a=M;