enbed.m

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

%--------将源图像逐行存储到矩阵S_I里
A=OpenBitmap('6.tif');
k=25;
l=32;
S_CI=DvdBptSubBp(A,k,l);
%--------将源图像逐行存储到矩阵S_CI里
%------------调用fastICA算法
[icasig, A_matrix, W_matrix] =FASTICA(S_CI);
%--------------------------------------------------------------------------

%---Step1------------------------------------------------------------------
%icasig         %估算的独立分量(行向量) 
%W_matrix       %估算的分离矩阵  
%A_matrix       %估算的混合矩阵
%--------------------将水印逐列存储到向量W里。-------------------------------
W=WMToV('wartermark.bmp');
%--------------------------------------------------------------------------

%----Step2-------将icasig最大的主成份icasig(1,:)和水印组和起来---------------
%---计算方差最大的主成份----------------------------------------------------
index=max_cov(icasig);
S_W=[icasig(index,:)
     W];
SS = [0.8 0.1
    0.9 0.2];
%--SS是混合矩阵-------------------------------------------------------------
%--------------------------------------------------------------------------


%-----Step3:Y=A*S_W--------------------------------------------------------
Y=SS*S_W;
y1=Y(1,:);
y2=Y(2,:);
%--------------------------------------------------------------------------



%----Step4 用y1代替icasig最大的主成份icasig(index,:)--------------------
icasig(index,:)=y1;
S_CI=inv(W_matrix)*icasig;%---或者S_CI=A_matrix*icasig;
SubBPtBP(A,k,l,S_CI);
%--------------------------------------------------------------------------

imshow(A);