six.rar.rar_six_steganography_steganography matlab_watermark cod

clear all clc; %保存开始时间 start_time=cputime; figure(1); %读出原始图像 subplot(1,2,1); input=imread('image.jpg'); imshow(input); title('原始图像'); %读出水印图像 subplot(1,2,2); watermarked_image=imread('watermarked.bmp'); %watermarked_image=imread('dwtnoise.bmp'); imshow(watermarked_image,[]); title('水印图像'); %三色分离 input=double(input); watermarked_image=double(watermarked_image); inputr=input(:,:,1); watermarked_imager=watermarked_image(:,:,1); inputg=input(:,:,2); watermarked_imageg=watermarked_image(:,:,2); inputb=input(:,:,3); watermarked_imageb=watermarked_image(:,:,3); %水印图像R的分解 [Cwr,Swr]=wavedec2(watermarked_imager,2,'haar'); %图像R的分解 [Cr,Sr]=wavedec2(inputr,2,'haar'); %水印图像G的分解 [Cwg,Swg]=wavedec2(watermarked_imageg,2,'haar'); %图像G的分解 [Cg,Sg]=wavedec2(inputg,2,'haar'); %水印图像B的分解 [Cwb,Swb]=wavedec2(watermarked_imageb,2,'haar'); %图像B的分解 [Cb,Sb]=wavedec2(inputb,2,'haar'); %提取水印的小波系数 r=0.08; for k=0:3 whr(k+1,:)=Cwr(1+size(Cwr,2)/4+k*size(Cwr,2)/16:... size(Cwr,2)/4+(k+1)*size(Cwr,2)/16)-... Cr(1+size(Cr,2)/4+k*size(Cr,2)/16:... size(Cr,2)/4+(k+1)*size(Cr,2)/16); wvr(k+1,:)=Cwr(1+size(Cwr,2)/2+k*size(Cwr,2)/16:... size(Cwr,2)/2+(k+1)*size(Cwr,2)/16)-... Cr(1+size(Cr,2)/2+k*size(Cr,2)/16:... size(Cr,2)/2+(k+1)*size(Cr,2)/16); wdr(k+1,:)=Cwr(1+3*size(Cwr,2)/4+k*size(Cwr,2)/16:... 3*size(Cwr,2)/4+(k+1)*size(Cwr,2)/16)-... Cr(1+3*size(Cr,2)/4+k*size(Cr,2)/16:... 3*size(Cr,2)/4+(k+1)*size(Cr,2)/16); end whr=(whr(1,:)+whr(2,:)+whr(3,:)+whr(4,:))/(4*r); wvr=(wvr(1,:)+wvr(2,:)+wvr(3,:)+wvr(4,:))/(4*r); wdr=(wdr(1,:)+wdr(2,:)+wdr(3,:)+wdr(4,:))/(4*r); war=(Cwr(1:size(Cwr,2)/16)-Cr(1:size(Cr,2)/16))/r; g=0.1; for k=0:3 whg(k+1,:)=Cwg(1+size(Cwg,2)/4+k*size(Cwg,2)/16:... size(Cwg,2)/4+(k+1)*size(Cwg,2)/16)-Cg(1+... size(Cg,2)/4+k*size(Cg,2)/16:... size(Cg,2)/4+(k+1)*size(Cg,2)/16); wvg(k+1,:)=Cwg(1+size(Cwg,2)/2+k*size(Cwg,2)/16:... size(Cwg,2)/2+(k+1)*size(Cwg,2)/16)-Cg(1+... size(Cg,2)/2+k*size(Cg,2)/16:... size(Cg,2)/2+(k+1)*size(Cg,2)/16); wdg(k+1,:)=Cwg(1+3*size(Cwg,2)/4+k*size(Cwg,2)/16:... 3*size(Cwg,2)/4+(k+1)*size(Cwg,2)/16)-... Cg(1+3*size(Cg,2)/4+k*size(Cg,2)/16:... 3*size(Cg,2)/4+(k+1)*size(Cg,2)/16); end whg=(whg(1,:)+whg(2,:)+whg(3,:)+whg(4,:))/(4*g); wvg=(wvg(1,:)+wvg(2,:)+wvg(3,:)+wvg(4,:))/(4*g); wdg=(wdg(1,:)+wdg(2,:)+wdg(3,:)+wdg(4,:))/(4*g); wag=(Cwg(1:size(Cwg,2)/16)-Cg(1:size(Cg,2)/16))/g; b=0.155; for k=0:3 whb(k+1,:)=Cwb(1+size(Cwb,2)/4+k*size(Cwb,2)/16:... size(Cwb,2)/4+(k+1)*size(Cwb,2)/16)-... Cb(1+size(Cb,2)/4+k*size(Cb,2)/16:... size(Cb,2)/4+(k+1)*size(Cb,2)/16); wvb(k+1,:)=Cwb(1+size(Cwb,2)/2+k*size(Cwb,2)/16:... size(Cwb,2)/2+(k+1)*size(Cwb,2)/16)-... Cb(1+size(Cb,2)/2+k*size(Cb,2)/16:... size(Cb,2)/2+(k+1)*size(Cb,2)/16); wdb(k+1,:)=Cwb(1+3*size(Cwb,2)/4+k*size(Cwb,2)/16:... 3*size(Cwb,2)/4+(k+1)*size(Cwb,2)/16)-... Cb(1+3*size(Cb,2)/4+k*size(Cb,2)/16:... 3*size(Cb,2)/4+(k+1)*size(Cb,2)/16); end whb=(whb(1,:)+whb(2,:)+whb(3,:)+whb(4,:))/(4*b); wvb=(wvb(1,:)+wvb(2,:)+wvb(3,:)+wvb(4,:))/(4*b); wdb=(wdb(1,:)+wdb(2,:)+wdb(3,:)+wdb(4,:))/(4*b); wab=(Cwb(1:size(Cwb,2)/16)-Cb(1:size(Cb,2)/16))/b; %重构水印图像 cwr=[war,whr,wvr,wdr]; swr(:,1)=[sqrt(size(war,2)),sqrt(size(war,2)),2*sqrt(size(war,2))]; swr(:,2)=[sqrt(size(war,2)),sqrt(size(war,2)),2*sqrt(size(war,2))]; wr=waverec2(cwr,swr,'haar'); figure(10),imshow(uint8(wr)); cwg=[wag,whg,wvg,wdg]; swg(:,1)=[sqrt(size(wag,2)),sqrt(size(wag,2)),2*sqrt(size(wag,2))]; swg(:,2)=[sqrt(size(wag,2)),sqrt(size(wag,2)),2*sqrt(size(wag,2))]; wg=waverec2(cwg,swg,'haar'); figure(11),imshow(uint8(wg)); cwb=[wab,whb,wvb,wdb]; swb(:,1)=[sqrt(size(wab,2)),sqrt(size(wab,2)),2*sqrt(size(wab,2))]; swb(:,2)=[sqrt(size(wab,2)),sqrt(size(wab,2)),2*sqrt(size(wab,2))]; wb=waverec2(cwb,swb,'haar'); figure(12),imshow(uint8(wb)); %三色叠加 temp=size(wr); pic=zeros(temp(1),temp(2),3); for i=1:temp(1); for j=1:temp(2); pic(i,j,1)=wr(i,j); pic(i,j,2)=wg(i,j); pic(i,j,3)=wb(i,j); end end output=uint8(round(pic)); %转化为uint8 watermark_image_uint8=uint8(output); imwrite(watermark_image_uint8,'extract.bmp','bmp'); figure(2); subplot(1,2,1); %imshow(watermark_image_uint8); imshow(uint8(wb)); title('提取水印'); subplot(1,2,2); imshow('palm101256.bmp'); title('原始水印'); %显示时间 elapsed_time=cputime-start_time;