wm_embedtt.m

这是<<Matlab在图像处理中的应用>>这本书的源码

function [icasig,mixedsig,img2,s12]=wm_embedtt(features,img)
[Mc,Nc]=size(img);
randn('state',10);
Iw=randn(30);
Ic=transform(Iw);
[mc,nc]=size(Ic);%
len=size(features,1);
features_e=cell(1,len);
rt=(mc-1)/2;
fea_e=cell(1,len);
img2=img;
for i=1:len
for j=-rt:rt
  for k=-rt:rt
   rf1=round(features(i,1));
   rf2=round(features(i,2));
   fea_e{1,i}(j+rt+1,k+rt+1)=img(rf2+j,rf1+k); 
  end
end
s1{1,i}=double(fea_e{1,i});
end

for jj=1:len
    for l=1:mc
        for k=1:nc
            rf1=round(features(jj,1));
            rf2=round(features(jj,2));
         end
    end
end
alp=17;
s2=alp*Ic;
alp1=0.75;bet1=0.25;al1=0.34;
for i=1:len
   s11{1,i}=s1{1,i}+alp1*s2;%nvf?
end

for jj=1:len
    for l=1:mc
        for k=1:nc
            rf1=round(features(jj,1));
            rf2=round(features(jj,2));
        img2(rf2-rt-1+k,rf1-rt-1+l)=s11{1,jj}(k,l);
        end
    end
end
figure,imshow(img2);
imwrite(img2,'watermarked.bmp');