getlhw.m

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

map=wwx([1000,2000,2000,3000]);
cof =histcofnew(map);
mapbw=im2bw(map,cof);
% ss=strel('line',1,90);
for ll=1:13
    ss=strel('line',2*ll+1,0);
    ermap=imdilate(mapbw,ss);
    immap=imdilate(~ermap,ss);
    clmap=~immap;
    cou(ll)=sum(sum(clmap));
 
end
discou=cou(1:length(cou)-1)-cou(2:length(cou));
discou1=cou(2:length(cou))-cou(1:length(cou)-1);








% map=wwx([1000,2000,2000,3000]);
% cof =histcofnew(map);
% mapbw=im2bw(map,cof);
% % ss=strel('line',1,90);
% for ll=1:13
%     ss=strel('line',2*ll+1,0);
%     ermap=imerode(mapbw,ss);
%     immap=imdilate(ermap,ss);
%     cou(ll)=sum(sum(immap));
%  
% end
% discou=cou(1:length(cou)-1)-cou(2:length(cou));