comatch.m

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

function objout=comatch(map,mod)
[maph,mapw]=size(map);
[modh,modw]=size(mod);

marea=modh*modw;
objout=[];
for ii=1:maph-modh+1
    for kk=1:mapw-modw+1
        mappart=map(ii:ii+modh-1,kk:kk+modw-1);
        logout=~xor(mappart,mod);
        matchout=sum(sum(logout))/marea;
        if matchout>.9
            if ~isempty(objout)
                rr=abs(ii-objout(:,1)); cc=abs(kk-objout(:,2));
                andout=(rr<modh-1)&(cc<modw-1);
                idx=find(andout==1);
                if~isempty(idx)
                    if objout(idx,3)<matchout
                        objout(idx,:)=[ii kk matchout];
                    end
                else
                    objout=[objout;ii kk matchout];
                end
            else
                objout=[objout;ii kk matchout];
            end
        end
    end
end