基于2x2块的Type-1水印提取算法(Matlab实现)

function rw = wm_extr_type1(loca)

% 基于2x2块的Type-1水印提取算法
% 输入参数：
%   loca: 包含嵌入水印的图像
% 输出参数：
%   rw: 提取的水印向量

[mm,nn] = size(loca);
index = 1;
rw = zeros(1,mm*nn/4); %预分配水印向量

% 取出2*2的块
for k = 1:2:mm-1
    for kk = 1:2:nn-1
        tem = loca(k:k+1,kk:kk+1);
        a11 = tem(1,1);
        a12 = tem(1,2);
        a21 = tem(2,1);
        a22 = tem(2,2);
        % -------a11+a21 >= a12+a22则水印为0，否则为1-------
        if a11+a21 >= a12+a22
            rw(index) = 0;
        else
            rw(index) = 1;
        end
        index = index + 1;
    end    
end
end