dvdbptsubbp.m

ICA是比较新的用于图像水印的方法

function Z=DvdBptSubBp(A,k,l)
%将位图A分割成(m/k)*(n/l)个子图，每个子图大小是k*l
%该函数返回一个矩阵,每一行都是一个子图


%------------------------------------------
dim=size(A);  %dim是由矩阵A的维数组成的向量
m_subbp=floor(dim(1)/k);
n_subbp=floor(dim(2)/l);
%------------------------------------------


%------------------------------------------
%构造一个矩阵，每一行代表一个子图
Z=ones(m_subbp*n_subbp,k*l);
%------------------------------------------


%------------------------------------------
%开始给矩阵Z附值
for i=1:m_subbp*n_subbp
row=ceil(i/n_subbp);%第i个子图所在行
column=n_subbp-(row*n_subbp-i);%第i个子图所在列
temp=A(((row-1)*k+1):row*k,((column-1)*l+1):column*l);
Z(i,:)=(temp(:))';
end
%------------------------------------------