glwdnew.m

来自「在图像处理中对于比较大地图像」· M 代码 · 共 44 行

%map=ww(6300,7300,2300,3300);
function [linewid,linedis,se]=glwdnew(mapbw)
%cof=histcof(map);
%mapen=mapenhance(map);
%mapbw=im2bw(mapen,cof);
sumc=sum(mapbw);
sumr=sum(mapbw');
idxc=find(sumc==0);
idxr=find(sumr==0);
if length(idxc)/length(sumc)>length(idxr)/length(sumr)
    se=ones(1,3); sumo=sumc;
else
    se=ones(3,1); sumo=sumr;
end
counts=sum(sum(double(mapbw)));
meanv=counts/length(sumo);
ll=sumo>meanv*.5;
rlong=length(ll);
rr1=ll(1:rlong-1); rr2=ll(2:rlong);
rout=xor(rr1,rr2);
ridx=find(rout);

if ll(1)==1
    ridx=ridx(2:length(ridx));
end
if ll(rlong)==1
    ridx=ridx(1:length(ridx)-1);
end
ridl=length(ridx);
a1=ridx(1:2:ridl-1);
b1=ridx(2:2:ridl);
wids=b1-a1;
linewid=round(median(wids));

ridout=round((a1+1+b1)/2);
long=length(ridout);
rdis=ridout(2:long)-ridout(1:long-1);
rdiss=round(median(rdis));
linedis=rdiss-linewid;