togridnew.m

在图像处理中对于比较大地图像

%map=ww(2500,3500,1500,2500);
function [rindx,cindx]=togrid(mapenbw,linewid,linedis,se)
%cof=histcof(map);
%mapbw=im2bw(map,cof);
%mapen=mapenhance(map);
%mapenbw=im2bw(mapen,cof);

timeo=round(linedis/2)-2;
timec=round((linewid+linedis)/2)-2;
if timeo<1,timeo=1;end
if timec<1,timec=1;end
    
mapbwor=dilate(mapenbw,se,timeo);
mapbwor=~mapbwor;
mapbwor=dilate(mapbwor,se,timeo);
mapbwor=~mapbwor;
map1=erode(mapbwor,se,timec);
mapo1=dilate(map1,se,timec);

mapbwoc=dilate(mapenbw,se',timeo);
mapbwoc=~mapbwoc;
mapbwoc=dilate(mapbwoc,se',timeo);
mapbwoc=~mapbwoc;
map2=erode(mapbwoc,se',timec);
mapo2=dilate(map2,se',timec);

[m,n]=size(se);
if m>n
    cs=sum(mapo1);rs=sum(mapo2');
else
    rs=sum(mapo1');cs=sum(mapo2);
end
lldis=linewid+linedis;
rindx=fgrid(rs,lldis);
cindx=fgrid(cs,lldis);