最简单的匹配算法matalb代码-绝对误差和算法（SAD）

本文链接：https://blog.youkuaiyun.com/Miracle_zk/article/details/87925099

博客详细解析了从网上获取的一个SAD算法Matlab代码存在的两个主要问题，并提供了修复方案。一是uint8数据类型在进行加减运算时被限制在0-255范围内，二是rectangle函数中长宽标识的错误。作者通过调整代码解决了这些问题，实现了模板匹配功能。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

从网上搜索了一个SAD的matlab代码，发现不能用。

唉唉，这是注定不让我轻松啊，

经过我认真的思考，发现了错误所在。

错误所在：

1、这个错误隐藏性很强，uint8的数据加减被限制在了 0-255 导致19-50=0的情况出现。

2、rectangle的长宽标识是反的。rectangle('position',[Y,X,N,N],'edgecolor','r');

效果：

模版

搜索图

上代码：


%%

%绝对误差和算法（SAD）

clear all;
close all;

%%

src=imread('Fig_rose.tif');

[a b d]=size(src);%行数（垂直分辨率）、列数（水平分辨率）、层数

if d==3
    src=rgb2gray(src);
end

mask=imread('mask.tif');

[m n d]=size(mask);

if d==3
    mask=rgb2gray(mask);
end

%%

N=n;%模板尺寸，默认模板为正方形

M=a;%代搜索图像尺寸，默认搜索图像为正方形

%%

dst=zeros(M-N+1,M-N+1);

for i=1:M-N+1         %子图选取，每次滑动一个像素

    for j=1:M-N+1
        
        temp=src(i:i+N-1,j:j+N-1);%当前子图

        dst(i,j)=sum(sum(abs(double(temp)-double(mask))));
    end
end


abs_min=min(min(dst));

[X,Y]=find(dst==abs_min);%最小值的位置

figure;

imshow(mask);title('模板');

figure;

imshow(src);

hold on;
rectangle('position',[Y,X,N,N],'edgecolor','r');%从点x，y开始绘制长宽为N-1的正方形

hold off;title('搜索图');

%%
%imwrite(uint8(src(200:260,440:500)),'mask.tif');