SATD算法Matlab代码

%//SATD模板匹配算法-哈达姆变换（hadamard）
clear all;
close all;
%%
src=double（imread（‘lena.jpg‘））;%//长宽相等的
mask=double（rgb2gray（imread（‘lena_mask.jpg‘）））;%//长宽相等的
M=size（src1）;%//搜索图大小
N=size（mask1）;%//模板大小
%%
hdm_matrix=hadamard（N）;%//hadamard变换矩阵
hdm=zeros（M-NM-N）;%//保存SATD值
for i=1:M-N
    for j=1:M-N
        temp=（src（i:i+N-1j:j+N-1）-mask）/256;
        sw=（hdm_matrix*temp*hdm_matrix）/256;