基于svm的嵌入与提取

% designed by superdont
% blog.csdn.net/superdont
clc;
clear all;
o=imread('E:\工作\科研\奇异值\image\lena.bmp');
od=double(o);
w=imread('E:\工作\科研\奇异值\image\w.bmp');
%对分好的块进行svd变换
[u1 d1 v1]=svd(od);
d12=d1+10*w;
[u2 d2 v2]=svd(d12);
ow=u1*d2*v1';
ow=uint8(ow);
figure,imshow(ow);
imwrite(ow,'E:\工作\科研\奇异值\image\lenawatermarked.bmp','bmp');
owr=imread('E:\工作\科研\奇异值\image\lenawatermarked.bmp');
owr=double(owr);
[tu1 td1 tv1]=svd(owr);
td12=u2*td1*v2';
tw=td12-d1
for i=1:512
    for j=1:512
        if tw(i,j)>5
            tw(i,j)=1;
        else
            tw(i,j)=0;
        end
    end
end
% tw=uint8(tw)
   figure,imshow(tw,[]);