【導(dǎo)讀】顯出來，這使得數(shù)字水印技術(shù)成為熱點問題。本論文研究的是基于離散余弦變換。試，包括壓縮攻擊、加噪聲攻擊、旋轉(zhuǎn)攻擊等。提取結(jié)果與測試結(jié)果達到了預(yù)期。最后還分別從嵌入方法、輸出結(jié)果、提取方法等方面比較了可見水印與。不可見水印的同異。

　　

【正文】 t Overview IEEE Signal Processing Magazine[J], 2021,17(5):2046. 2 向德生 ,楊格蘭 ,熊岳山 . 數(shù)字水印技術(shù)研究 [J]. 計算機工程與設(shè)計 ,2021,26(2):326328,334 3 李永全 . 基于 MATLAB 的 DCT 域數(shù)字水印技術(shù)實現(xiàn) [J]. 信息技術(shù) ,2021,(4):6671 4 胡彥，陳昭 . MATLAB 在數(shù)字水印中的應(yīng)用 [J]. 計算機工程，2021,(7):184186 5 周熠 . 數(shù)字水印及其方法研究 [J]. 通信技術(shù) ,2021,(10):124125 6 許文麗 . 基于版權(quán)保護的圖像數(shù)字水印研究 [D]. 西安電子科技大學(xué)，2021:2931 7 杜江 . 信息隱藏與數(shù)字水印技術(shù)研究 [D]. 西安電子科技大學(xué)，2021:6068 8 Ejima, M. and A. Myazaki. On the evaluation of performance of digital watermarking in the frequency domain[C]. Proc. of the IEEE Int. Conf. on Image Processing,2021,2: 546549. 9 任曉揚 ,韓勇 . 基于 DCT數(shù)字水印算法的 Matlab 實現(xiàn) [J]. 儀器儀表用戶 ,2021,(1)： 116117 10 黃曉 ,李純厚 ,李毅 . 基于 DCT 數(shù)字水印的 Matlab 仿真實現(xiàn) [J]. 電子技術(shù)， 2021,(11)： 5859 11 歐陽琦 ,李炳法 ,林 春薔 ,龐堅 ,高東妮 . 基于 DCT 域的圖像數(shù)字水印技術(shù) [J]. 信息技術(shù) , 2021,(9):2023 12 徐奔 . 基于 MATLAB 的 DCT 域數(shù)字圖像水印技術(shù) [J]. 計算機安全 ,2021,(8):3537 13 楊帆等 . 數(shù)字圖像處理與分析 [M]. 北京航空航天大學(xué)出版社，2021:4950 14 曹華 . 圖像和視頻水印嵌入新方法研究 [D]. 華中科技大學(xué)， 2021:3035 15 Nikolaidis,A. and . Asymptotically optimal detection for additive watermarking in the DCT and DWT domains [C]. IEEE Trans. Image Processing, 2021,2(10): 563571. 河北工業(yè)大學(xué) 城市學(xué)院 2021 屆本科畢業(yè)設(shè)計說明書 24 16 徐先傳 ,張琦 . 一種基于 DCT 域的數(shù)字水印方案研究 [J]. 微計算機信息 , 2021,(12):5961 17 徐金東 ,黎洪松 ,倪夢瑩 . 一種基于關(guān)系的 DCT 域數(shù)字水印改進算法 [J]. 北京師范大學(xué)學(xué)報 (自然科學(xué)版 ), 2021,(1):5759 18 曾高榮 . 數(shù)字水印 的魯棒性評測建模和算法研究 [D]. 北京交通大學(xué)，2021:2535 19 楊樹國 . 魯棒圖像數(shù)字水印技術(shù)的研究 [D]. 哈爾濱工程大學(xué)， 2021:26 20 Voloshynovskiy, S., S. Pereira and T. Pun. Attacks on Digital Watermarks: Classification, EstimationBased Attacks, and Benchmarks[J]. Comm. Magazine, 2021, 39(8): 118126. 河北工業(yè)大學(xué) 城市學(xué)院 2021 屆本科畢業(yè)設(shè)計說明書 25 致謝 首先對我的指導(dǎo)老師王楊老師表示深深地敬意和由衷的感謝，感謝王老師在畢業(yè)設(shè)計期間給予我的深深教誨、無私的關(guān)懷和悉心的指導(dǎo)。王老師淵博的學(xué)識、敏銳的思維、嚴謹?shù)闹螌W(xué)態(tài)度和孜孜不倦的工作作風(fēng)，都令我受益匪淺，為我以后踏入社會打下了堅實的基礎(chǔ)。通過這次畢業(yè)設(shè)計，我不僅在圖像處理方面學(xué)到了更深的知識，并且在數(shù)字水印及 MATLAB 軟件應(yīng)用方面有了比較大的提高，培養(yǎng)了自己的自學(xué)能力，這些進步對在以后的工作學(xué)習(xí)中起到了很好的幫助。 同時，感謝我周圍的同學(xué)們，是他們營造 出積極向上的學(xué)習(xí)環(huán)境和溫暖如家的生活氛圍，給我鼓勵和關(guān)懷，并且還要感謝所有在論文完成期間支持我，幫助我、鼓勵我的朋友。 向參加論文評審、答辯的專家和老師表示衷心的感謝！ 河北工業(yè)大學(xué) 城市學(xué)院 2021 屆本科畢業(yè)設(shè)計說明書 26 附錄 A 以下是不可見水印的嵌入程序，由于基本原理相同，只是在嵌入的位置有所不同，所以省去了可見水印的嵌入程序。 %嵌入前的準備 Y1=imread(39。39。)。 %原圖像 Y2=imread(39。39。)。 %待嵌入的水印圖像 figure(1)。 subplot(2,1,1)。imshow(Y1)。 subplot(2,1,2)。imshow(Y2)。 G1=rgb2gray(Y1)。 %將彩色圖像變成灰度圖像 G2=rgb2gray(Y2)。 figure(2)。 subplot(2,1,1)。imshow(G1)。 subplot(2,1,2)。imshow(G2)。 D1=dct2(G1)。 %進行 DCT變換 D2=dct2(G2)。 %嵌入水印 [x,y]=size(D1)。 %求出原圖像 DCT域的數(shù)組的行列數(shù) [m,n]=size(D2)。 %求出帶嵌入水印 DCT域的數(shù)組的行列數(shù) x1=round(x/300)。 %假設(shè)前 1/300后為中低頻 D0=zeros(x,y)。 %創(chuàng)建全零數(shù)組，為嵌入水印做準備 i=1。 %i為待嵌入水印的變換域的行坐標，初始值為 1 j=1。 %j為待嵌入水印的變換域的列坐標，初始值為 1 for u=x1:1:x for v=1:1:y D0(u*x+v)=D2(i,j)。 %將水印 DCT域的值寫到零數(shù)組的中頻段，為把水印 %嵌入原圖的中頻段做準備 if j==n %j=n時表示到達了一行的行尾 河北工業(yè)大學(xué) 城市學(xué)院 2021 屆本科畢業(yè)設(shè)計說明書 27 if i==m %i=m表示到達了最后一列 break。 % j=n并且 i=m時代表到達了最后一個元素，退出內(nèi)層循環(huán) else %若不是最后一個元素，僅是行尾，但不是列尾 j=1。 %列指向列首 i=i+1。 %行指向下一行 end else j=j+1。 %不是一行的末尾元素，則指向下一個元素 end end if j==n amp。amp。 i==m %既是列尾，又是行尾 break。 %則為最后一個元素，退出外層循環(huán) end end OD=D0*+D1。 %包含水印頻域信息的零數(shù)組與原圖像的變換域相疊加， 。 %寫入文件 O=idct2(OD)。 UO=uint8(round(O))。 %轉(zhuǎn)換為無符號整數(shù)，為保存嵌入水印后的圖像做準備 imwrite(UO,39。39。,39。jpg39。) %將嵌入水印后的圖像寫入文件保存 河北工業(yè)大學(xué) 城市學(xué)院 2021 屆本科畢業(yè)設(shè)計說明書 28 附錄 B 以下是不可見水印的提取程序，由于基本原理相同，只是提取的位置有所不同，所以省去了可見水印的提取程序。 %水印的提取 前的準備 clc。 clear。 Y1=imread(39。39。)。 %原圖像 Y2=imread(39。39。)。 %已嵌入水印后的圖像 G1=rgb2gray(Y1)。 %將彩色圖像變成灰度圖像 G2=Y2。 D1=dct2(G1)。 %進行 DCT變換 D2=dct2(G2)。 %提取水印 D0=D2D1。 %頻域相減，提取水印 [x,y]=size(D0)。 %求出相減后的行列數(shù) [m,n]=size(D3)。 %求出水印 DCT域的行列數(shù) x1=round(x/300)。 %只有嵌入者知道水印是從何處開始嵌入的 u=x1。 %u是水印提取后在頻域的行起點位置 v=1。 %v是提取后在頻域的列號，初始值為 1 OD=zeros(m,n)。 %構(gòu)建與水印的變換域相同大小的零數(shù)組，用來裝載水印 for i=1:1:m for j=1:1:n OD(i,j)=D0(u*x+v).*。 %將處于中頻的水印調(diào)整到從零頻開始，并增強 if v==y v=1。 %如果為最后一列，換行 u=u+1。 else v=v+1。 %如果不是最后一列，不換行，但要指向下一行 河北工業(yè)大學(xué) 城市學(xué)院 2021 屆本科畢業(yè)設(shè)計說明書 29 end %因為原圖的變換域一定大于水印的變換域 ，所以 不會出現(xiàn)最后一 %行最后一列的情況 end end O=idct2(OD)。 h=fspecial(39。average39。,[4 4])。 %均值濾波 O=imfilter(O,h)。 h=fspecial(39。unsharp39。,)。 %銳化 O=imfilter(O,h,39。replicate39。)。 imshow(O,[0,255])。