资源描述：
《matlab实现的小波变换彩色图像水印嵌入和提取程序代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、functiondwtglclearallclc;%保存开始时间start_time=cputime;figure(1);%读出原始图像subplot(2,2,1);input=imread('image.bmp');imshow(input);title('原始图像');%读出水印subplot(2,2,2);water=imread('watermark.bmp');imshow(water);title('水印');%三色分离input=double(input);water=double(water);inputr=input(:,:,1);
2、waterr=water(:,:,1);inputg=input(:,:,2);waterg=water(:,:,2);inputb=double(input(:,:,3));waterb=double(water(:,:,3));%系数r大，增加鲁棒性，r小增加透明性r=0.04;%水印R的分解[Cwr,Swr]=wavedec2(waterr,1,'haar');%图像R的分解[Cr,Sr]=wavedec2(inputr,2,'haar');%水印的嵌入Cr(1:size(Cwr,2)/16)=...Cr(1:size(Cwr,2)/16)+r
3、*Cwr(1:size(Cwr,2)/16);k=0;whilek<=size(Cr,2)/size(Cwr,2)-1Cr(1+size(Cr,2)/4+k*size(Cwr,2)/4:size(Cr,2)/4+...(k+1)*size(Cwr,2)/4)=Cr(1+size(Cr,2)/4+...k*size(Cwr,2)/4:size(Cr,2)/4+(k+1)*size(Cwr,2)/4)+...r*Cwr(1+size(Cwr,2)/4:size(Cwr,2)/2);Cr(1+size(Cr,2)/2+k*size(Cwr,2)/4:siz
4、e(Cr,2)/2+...(k+1)*size(Cwr,2)/4)=Cr(1+size(Cr,2)/2+...k*size(Cwr,2)/4:size(Cr,2)/2+(k+1)*size(Cwr,2)/4)+...r*Cwr(1+size(Cwr,2)/2:3*size(Cwr,2)/4);Cr(1+3*size(Cwr,2)/4+k*size(Cwr,2)/4:3*size(Cwr,2)/4+...(k+1)*size(Cwr,2)/4)=Cr(1+3*size(Cr,2)/4+...k*size(Cwr,2)/4:3*size(Cr,2)/4+
5、(k+1)*size(Cwr,2)/4)+...r*Cwr(1+3*size(Cwr,2)/4:size(Cwr,2));k=k+1;end;Cr(1:size(Cwr,2)/4)=Cr(1:size(Cwr,2)/4)+r*Cwr(1:size(Cwr,2)/4);g=0.02;%水印G的分解[Cwg,Swg]=WAVEDEC2(waterg,1,'haar');%图像G的分解[Cg,Sg]=WAVEDEC2(inputg,2,'haar');%水印的嵌入Cg(1:size(Cwg,2)/16)=...Cg(1:size(Cwg,2)/16)+g*
6、Cwg(1:size(Cwg,2)/16);k=0;whilek<=size(Cg,2)/size(Cwg,2)-1Cg(1+size(Cg,2)/4+k*size(Cwg,2)/4:size(Cg,2)/4+...(k+1)*size(Cwg,2)/4)=Cg(1+size(Cg,2)/4+...k*size(Cwg,2)/4:size(Cg,2)/4+(k+1)*size(Cwg,2)/4)+...g*Cwg(1+size(Cwg,2)/4:size(Cwg,2)/2);Cg(1+size(Cg,2)/2+k*size(Cwg,2)/4:size
7、(Cg,2)/2+...(k+1)*size(Cwg,2)/4)=Cg(1+size(Cg,2)/2+...k*size(Cwg,2)/4:size(Cg,2)/2+(k+1)*size(Cwg,2)/4)+...g*Cwg(1+size(Cwg,2)/2:3*size(Cwg,2)/4);Cg(1+3*size(Cg,2)/4+k*size(Cwg,2)/4:3*size(Cg,2)/4+...(k+1)*size(Cwg,2)/4)=Cg(1+3*size(Cg,2)/4+...k*size(Cwg,2)/4:3*size(Cg,2)/4+(k+
8、1)*size(Cwg,2)/4)+...g*Cwg(1+3*size(Cwg,2)/4:size(Cwg,
