【導(dǎo)讀】等作品，可以進(jìn)行學(xué)術(shù)交流。而一些具有特殊意義的數(shù)字信息，如涉及司法訴訟、政府機(jī)要等信息，更是。圖像、視頻制品中迅速得到廣泛的研究和發(fā)展。建立產(chǎn)品所有權(quán)、辨認(rèn)購買者或提供數(shù)字產(chǎn)品的一些額外信息。像與原始圖像沒有絲毫差別。對特定的攻擊進(jìn)行設(shè)計(jì)。水印技術(shù)的安全性是其最重要的特性，由于它的商業(yè)性，其算法必須公開，算法的安全性完全取決于密鑰，而不對算法進(jìn)行保密。到了保護(hù)產(chǎn)權(quán)的作用。據(jù)是否被修改過或是否經(jīng)過特殊的處理。所發(fā)送的數(shù)據(jù)中來給傳統(tǒng)的信號發(fā)射系統(tǒng)升級。這一運(yùn)用的目的是提取對原始信號進(jìn)行處理過的信息。數(shù)字水印技術(shù)便是二十世紀(jì)九十年代初出現(xiàn)的為解決這。2020年，王向陽等提出了一種DCT域自適應(yīng)彩色圖像二維數(shù)字水。其他DCT域數(shù)字水印算法，采用的多是基于DCT的8×8圖像塊。數(shù)字水印技術(shù)是近十年來橫跨信。保護(hù)及其沖突問題的解決。盲水印的魯棒性較強(qiáng)，而其應(yīng)用受到存儲(chǔ)成本的限制。不過，目前學(xué)術(shù)界研究的

　　

【正文】 tmark_1(i,j)=0。 end end end %%%%%%%%%% 計(jì)算 NC（歸一化相關(guān)系數(shù)） %%%%%%%% g_mark=double(tmark_1)。 o_mark=double(omark)。 [m,n]=size(g_mark)。 nc_0=0。 淮北煤炭師范學(xué)院學(xué)士畢業(yè)論文 基于變換域的數(shù)字水印技術(shù)的研究 26 nc_1=0。 nc_2=0。 for i=1:m for j=1:n nc_0=nc_0+g_mark(i,j)*o_mark(i,j)。 nc_1=nc_1+o_mark(i,j)*o_mark(i,j)。 nc_2=nc_2+g_mark(i,j)*g_mark(i,j)。 end end NC=nc_0/sqrt(nc_1*nc_2)。 %%%%%% 顯示提取水印 %%%%%%%% subplot(2,3,5),imshow(39。39。),title(‘原始水印圖像’ )。 subplot(2,3,6),imshow(tmark_1,[])。 name=‘提取得水印圖像’ 。 title(strcat(num2str(name),39。NC=39。,num2str(NC)))。 附錄 B（基于 FFT 域的數(shù)字水印算法程序） %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%% 嵌入水印算法 %%%%%%%%%%%%%% clear all。 clc。 %%%%%% 讀取水印 圖象 %%%%%%% omark=double(imread(39。39。))/255。 mo=size(omark,1)。 no=size(omark,2)。 %%%%% 以下生成水印信息 %%%%%% mark=omark。 alpha=200。 R1=randn(1,8)。 rand1=round(R1)。 R2=randn(1,8)。 rand2=round(R2)。 %%%%%%% 讀取宿主圖像 %%%%%%% cimage=imread(39。39。)。 淮北煤炭師范學(xué)院學(xué)士畢業(yè)論文 基于變換域的數(shù)字水印技術(shù)的研究 27 figure(1)。 subplot(2,3,1),imshow(cimage,[]),title(39。原始的宿主圖像 39。)。 [mc,nc]=size(cimage)。 cimage=double(cimage)。 fft0=blkproc(cimage,[8,8],39。fft239。)。 %%%%%%% 嵌入水印 %%%%%%%%% fft_abs=abs(fft0)。 fft_angle=angle(fft0)。 fft1=fft_abs。 for ii=1:mo for jj=1:no x=(ii1)*8。y=(jj1)*8。 if mark(ii,jj)==1 k=rand1。 else k=rand2。 end fft1(x+1,y+8)=fft_abs(x+1,y+8)+alpha*k(1)。 fft1(x+2,y+7)=fft_abs(x+2,y+7)+alpha*k(2)。 fft1(x+3,y+6)=fft_abs(x+3,y+6)+alpha*k(3)。 fft1(x+4,y+5)=fft_abs(x+4,y+5)+alpha*k(4)。 fft1(x+5,y+4)=fft_abs(x+5,y+4)+alpha*k(5)。 fft1(x+6,y+3)=fft_abs(x+6,y+3)+alpha*k(6)。 fft1(x+7,y+2)=fft_abs(x+7,y+2)+alpha*k(7)。 fft1(x+8,y+1)=fft_abs(x+8,y+1)+alpha*k(8)。 end end fft_abs_0=fft1。 fft_abs_1=(fft_abs_0.*exp(i*fft_angle))。 wimage_0=blkproc(fft_abs_1,[8,8],39。ifft239。)。 wimage=abs(wimage_0)。 %%%%%% 顯示嵌入水印圖象 %%%%% wimage_1=uint8(wimage)。 imwrite(wimage_1,39。39。,39。bmp39。)。 subplot(2,3,2),imshow(wimage,[]),title(39。嵌入水印后的圖象 39。)。 %%%%%%%%%%%%% 進(jìn)行各種攻擊 %%%%%%%%%%%% 淮北煤炭師范學(xué)院學(xué)士畢業(yè)論文 基于變換域的數(shù)字水印技術(shù)的研究 28 disp(39。1加入白噪聲 39。)。 disp(39。2高斯低通濾波 39。)。 disp(39。3剪切圖像 39。)。 disp(39。4旋轉(zhuǎn)圖象 39。)。 disp(39。5直接檢測 39。)。 begin=input(39。請選擇攻擊 （ 15） 39。)。 switch begin %%%%%%%%加入白噪聲 %%%%%%%% case 1 Aimage1=wimage。 Wnoise=12*randn(size(Aimage1))。 Aimage1=Aimage1+Wnoise。 subplot(2,3,4),imshow(Aimage1,[]),title(39。加入白噪聲后的圖象 39。)。 att=Aimage1。 imwrite(att,39。39。)。 %%%%%% 高斯低通濾波 %%%%%%%% case 2 Aimage2=wimage。 H=fspecial(39。gaussian39。,[4,4],)。 Aimage2=imfilter(Aimage2,H)。 subplot(2,3,4),imshow(Aimage2,[]),title(39。高斯低通濾波后的圖象 39。)。 att=Aimage2。 imwrite(att,39。39。) %%%%%%% 剪切攻擊 %%%%%%%%% case 3 Aimage3=wimage。 Aimage3(1:128,1:128)=256。 subplot(2,3,4),imshow(Aimage3,[]),title(39。剪切后的圖象 39。)。 att=Aimage3。 imwrite(att,39。39。)。 %%%%%%% 旋轉(zhuǎn)攻擊 %%%%%%%%% case 4 Aimage4=wimage。 Aimage4=imrotate(Aimage4,10,39。bilinear39。,39。crop39。)。 Aimage_4=mat2gray(Aimage4)。 淮北煤炭師范學(xué)院學(xué)士畢業(yè)論文 基于變換域的數(shù)字水印技術(shù)的研究 29 subplot(2,3,4),imshow(Aimage_4,[]),title(39。旋轉(zhuǎn) 10度后的圖象 39。)。 att=Aimage_4。 imwrite(att,39。39。)。 %%%%%%% 沒有受到攻擊 %%%%%%% case 5 subplot(2,3,4),imshow(wimage,[]),title(39。直接提取的圖象 39。)。 att=wimage。 imwrite(att,39。39。)。 end %%%%%%%%%%%%%%% 提取水印 %%%%%%%%%%%%% tmark_0=att。 tmark_0=blkproc(tmark_0,[8,8],39。fft239。)。 tmark_0=abs(tmark_0)。 pass=zeros(1,8)。 for ii=1:mo for jj=1:no x=(ii1)*8。y=(jj1)*8。 pass(1)=tmark_0(x+1,y+8)。 pass(2)=tmark_0(x+2,y+7)。 pass(3)=tmark_0(x+3,y+6)。 pass(4)=tmark_0(x+4,y+5)。 pass(5)=tmark_0(x+5,y+4)。 pass(6)=tmark_0(x+6,y+3)。 pass(7)=tmark_0(x+7,y+2)。 pass(8)=tmark_0(x+8,y+1)。 if (corr2(pass,rand1)corr2(pass,rand2)) tmark_1(ii,jj)=1。 else tmark_1(ii,jj)=0。 end end end %%%%%%%%%% 計(jì)算 NC（歸一化相關(guān)系數(shù)） %%%%%%%% g_mark=double(tmark_1)。 o_mark=double(omark)。 淮北煤炭師范學(xué)院學(xué)士畢業(yè)論文 基于變換域的數(shù)字水印技術(shù)的研究 30 [m,n]=size(g_mark)。 nc_0=0。 nc_1=0。 nc_2=0。 for ii=1:m for jj=1:n nc_0=nc_0+g_mark(ii,jj)*o_mark(ii,jj)。 nc_1=nc_1+o_mark(ii,jj)*o_mark(ii,jj)。 nc_2=nc_2+g_mark(ii,jj)*g_mark(ii,jj)。 end end NC=nc_0/sqrt(nc_1*nc_2)。 %%%%%% 顯示提取水印圖象 %%%%% subplot(2,3,5),imshow(39。39。),title(39。原始水印圖象 39。)。 subplot(2,3,6),imshow(tmark_1,[])。 name=39。提取的水印圖像 39。 title(strcat(num2str(name),39。NC=39。,num2str(NC)))。 致 謝 經(jīng)過了四年的努力，我終于完成了我的學(xué)業(yè)，并順利完成了學(xué)士畢業(yè)論文，在這期間得到了諸多老師的指導(dǎo)與幫助，還有許 多 親人朋友 的關(guān)心與支持，在這里一并表示衷心的感謝！ 首先，感謝竇德召老師，他在我 寫畢業(yè)論文期間給我提供了許多資料，并指導(dǎo)我寫完了畢業(yè)論文，老師嚴(yán)謹(jǐn)樸實(shí)的治學(xué)態(tài)度，誨人不倦的教學(xué)精神，和藹可親的處世之道給我留下了很深刻的映像，讓我受益無窮。另外，還要感謝圖書館電子資源管理室，它在我查閱資料方面也給我提供了許多幫助！ 感謝所有關(guān)心過我的人 ，謹(jǐn) 以我的努力作為回報(bào)！