【正文】 gne_source,watermark_source,watermarked_source)%N=512。%先處理水印圖像(128x128) ，如下：arnold_image=Arnold(watermark_source,10,0)。%讀入水印圖像并且進(jìn)行 arnold 變換，10 表示進(jìn)行的次數(shù)（可作密鑰）%arnold_image 是 uint8 的格式%[ca,ch,cv,cd]=dwt2(arnold_image,39。haar39。)。%進(jìn)行一級(jí) dwt 變換%然后處理載體圖像，如下：origne_image=imread(origne_source)。%讀入載體圖像dims=ndims(origne_image)。if(dims==2)[image_rol,image_col]=size(origne_image)。if(image_rol512) if(image_col512) disp(39。載體圖像的高度和寬度都必須大于 51239。)。 return。 endendorigne_image2=double(origne_image)。part_origne_image2=origne_image2。[LL1,LH1,HL1,HH1]=dwt2(part_origne_image2,39。haar39。)。%對(duì)載體圖像先進(jìn)行一級(jí) dwt 變換（320x256）[LL2,LH2,HL2,HH2]=dwt2(LL1,39。haar39。)。%對(duì) LL1 進(jìn)行 dwt 變換（160x128）[LL3,LH3,HL3,HH3]=dwt2(LL2,39。haar39。)。%變成 80x64，水印圖像為 64x64%水印嵌入過(guò)程，如下：%嵌入分兩個(gè)部分%a1=。a2=。arnold_image=double(arnold_image)。[UW,SW,VW]=svd(arnold_image)。%對(duì)水印作 svd 分解[UL,SL,VL]=svd(LL3)。%對(duì) 3 級(jí)低頻分量作 svd 分解 80x64 [UH,SH,VH]=svd(HH3)。%對(duì) 3 級(jí)高頻分量作 svd 分解SHH=SH。U=UW。V=VW。%128x128g1=。g2=。g3=。%由于第一個(gè)奇異值很大，改變一點(diǎn)，圖像失真很多，所以強(qiáng)度不應(yīng)很大[a,b]=size(LL3)。col=min(a,b)。SL(1,1)=SL(1,1)+g3*SW(1,1)。for i=2:col SL(i,i)=SL(i,i)+g1*SW(i,i)。endSH(1,1)=SH(1,1)+g3*SW(65,65)。for i=2:col SH(i,i)=SH(i,i)+g2*SW(i+64,i+64)。endLL3=UL*SL*VL39。HH3=UH*SH*VH39。LL2=idwt2(LL3,LH3,HL3,HH3,39。haar39。)。%小波逆變換29LL1=idwt2(LL2,LH2,HL2,HH2,39。haar39。)。watermarked_image1=idwt2(LL1,LH1,HL1,HH1,39。haar39。)。watermarked_image=uint8(watermarked_image1)。imwrite(watermarked_image,watermarked_source)。figuresubplot(1,2,1)。imshow(origne_image)。title(39。原始圖像39。)。subplot(1,2,2)。imshow(watermarked_image)。title(39。嵌入水印后的圖像 39。)。end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%文件名：%水印提取函數(shù)%參數(shù)說(shuō)明：%watermarked_image嵌入水印的圖像%origne_source原始圖像function water_get=out(watermarked_source,origne_source,UW,VW)origne_image=imread(origne_source)。dims=ndims(origne_image)。if(dims~=2) disp(39。輸入的不是灰度圖像39。)。end[o_rol,o_col]=size(origne_image)。watermarked_image2=imread(watermarked_source)。[w_rol,w_col]=size(watermarked_image2)。origne_image1=double(origne_image)。watermarked_image1=double(watermarked_image2)。 watermarked_image(i+temp_rol,j+temp_col)=watermarked_image2(i,j)。part_origne_image=origne_image1。part_watermarked_image=watermarked_image1。[o_LL1,o_LH1,o_HL1,o_HH1]=dwt2(part_origne_image,39。haar39。)。%對(duì)載體圖像先進(jìn)行 dwt 變換320x256[o_LL2,o_LH2,o_HL2,o_HH2]=dwt2(o_LL1,39。haar39。)。%160*126[o_LL3,o_LH3,o_HL3,o_HH3]=dwt2(o_LL2,39。haar39。)。%80*64[w_LL1,w_LH1,w_HL1,w_HH1]=dwt2(part_watermarked_image,39。haar39。)。[w_LL2,w_LH2,w_HL2,w_HH2]=dwt2(w_LL1,39。haar39。)。[w_LL3,w_LH3,w_HL3,w_HH3]=dwt2(w_LL2,39。haar39。)。SW1=zeros(128)。%%ch=(w_LH3o_LH3)/a1。cv=(w_HL3o_HL3)/a2。cd=(w_HH3o_HH3)/a3。[o_UL,o_SL,o_VL]=svd(o_LL3)。%對(duì)原始圖像分量作 svd 分解[o_UH,o_SH,o_VH]=svd(o_HH3)。%對(duì)原始圖像分量作 svd 分解[w_UL,w_SL,w_VL]=svd(w_LL3)。%對(duì)嵌入水印分量作 svd 分解[w_UH,w_SH,w_VH]=svd(w_HH3)。%對(duì)嵌入水印分量作 svd 分解%[w_ULH,w_SLH,w_VLH]=svd(w_LH3)。%[w_UHL,w_SHL,w_VHL]=svd(w_LL3)。[a,b]=size(w_LL3)。col=min(a,b)。if(col64) col=64。endtemp=fix((64col)/2)。if(temp0) temp=0。endg1=。g2=。g3=。SW1(1,1)=(w_SL(1,1)o_SL(1,1))/g3。for i=2:col SW1(i,i)=(w_SL(i,i)o_SL(i,i))/g1。附錄30endSW1(65,65)=(w_SH(1,1)o_SH(1,1))/g3。for j=2:col SW1(j+64,j+64)=(w_SH(j,j)o_SH(j,j))/g2。end SW=SW1。%SW(1,1)=。SW(1,1)=SW(1,1)。water_image=UW*SW*VW39。water_image=uint8(water_image)。imwrite(water_image,39。39。)。water_temp=Arnold(39。39。,10,1)。water_get=uint8(water_temp)。%figure。data=imread(39。39。)。subplot(1,2,1)。imshow(data)。title(39。原水印圖像39。)。subplot(1,2,2)。imshow(water_temp)。title(39。提取的水印圖像39。)。imwrite(water_temp,39。39。)。%%%%%%%%%%%性能分析%%%%%%%%%%%%%%%%%%% 攻擊實(shí)驗(yàn) 測(cè)試魯棒性 %%%%%%%%%%%disp(39。對(duì)嵌入水印的圖像的攻擊實(shí)驗(yàn)，請(qǐng)輸入選擇項(xiàng)：39。)。disp(39。1添加白噪聲 39。)。disp(39。2高斯低通濾波 39。)。disp(39。3JPEG 壓縮39。)。disp(39。4圖像剪切 39。)。disp(39。5旋轉(zhuǎn)變換 39。)。disp(39。6直接檢測(cè)水印 39。)。disp(39。7中值濾波 39。)。disp(39。8馬賽克攻擊 39。)。disp(39。其他不攻擊 39。)。d=input(39。請(qǐng)輸入選擇（ 18）:39。)。origne_image_name=39。39。 a1=imread(39。39。)。 im_prime=imread(origne_image_name)。 water_prime=imread(39。39。)。 a1=double(a1)。 switch d case 6 [M,N]=size(a1)。 water_get=out(39。39。,origne_image_name,U,V)。 psnr_arr=psnr(a1,im_prime,M,N)。 nc_arr=nc(water_get,water_prime)。 disp(39。嵌入水印后的圖像：39。)。 disp(39。峰值信噪比：psnr=39。)。disp(psnr_arr)。 disp(39。相關(guān)性：nc=39。)。disp(nc_arr)。 subplot(1,2,1)。imshow(a1,[])。title(39。未受攻擊的含水印圖像39。)。 subplot(1,2,2)。imshow(water_get)。title(39。水印圖像39。)。 M1=a1。 M_1=uint8(M1)。 imwrite(M_1,39。39。,39。bmp39。)。 case 1 WImage2=a1。 count=0。 for i=0:20:100 count=count+1。 arr1(count)=i。31 noise0=i*randn(size(WImage2))。 WImage2=WImage2+noise0。%subplot(2,3,4)。 %imshow(WImage2,[])。 %title(39。加入白噪聲后圖像 39。)。 M1=WImage2。 [M2,N2]=size(WImage2)。 M_1=uint8(M1)。 name=[39。whitenoise_39。 int2str(i) 39。.bmp39。]。 psnr_arr1(count)=psnr(WImage2,im_prime,M2,N2)。 imwrite(M_1,name,39。bmp39。)。 water_get=out(name,origne_image_name,U,V)。 nc_arr1(count)=nc(water_get,water_prime)。 ber_arr1(count)=ber(M1,a1)。 end arr1=0:20:100。 name=[39。whitenoise_39。 int2str(40) 39。.bmp39。]。 water_get=out(name,origne_image_name,U,V)。 one_water_image=imread(name)。 figure。 subplot(1,2,1)。imshow(one_water_image)。title(39。加入白噪聲－４０后的圖像39。)。 subplot(1,2,2)。imshow(water_get)。title(39。加噪后提取的水印39。)。 imwrite(water_get,39。39。)。