dctdwt)图像数字水印的方案毕业论文基于dct和dwt的算法有详细程序

dctdwt)图像数字水印的方案毕业论文基于dct和dwt的算法有详细程序内容摘要：

加各种噪声： imnoise()可以对图像加入各种噪声，如白噪声、椒盐噪声等，加入噪声是对水印鲁棒性考验的一种常见的攻击； 放大 /缩小： imresize()可以以指定的插值方法来对图像进行放大和缩小。 第 四 章： 图像数字水印技术的实现 基于离散余弦变法（ DCT）实现数字水印技术 ① 打开原始及水印图像 ： subplot(2,2,1) I=uigetfile(39。 *.bmp39。 ,39。 打开原始彩 色图像文件 39。 )。 RGB=imread(I)。 image(RGB)。 title(39。 原始彩色图像 39。 )。 subplot(2,2,2) I=uigetfile(39。 *.bmp39。 ,39。 打开水印灰度图像文件 39。 )。 imshow(I)。 title(39。 灰度水印图像 39。 )。 subplot(2,2,3) H=imread(I)。 J=dct2(H)。 imshow(log(abs(J)),[]),colorbar。 title(39。 水印图像经 DCT 变换后能量分布情况 39。 ) 运行结果： ② 水印全过程： 0%水印加入程序 Q=input(39。 请输入放缩因子的值 (建议小于 1)： Q＝ 39。 ) subplot(2,3,1) RGB=imread(39。 浙江台州学院 39。 ,39。 jpg39。 )。 imshow(RGB)。 title(39。 原始图像 39。 )。 subplot(2,3,2) N=dct2(RGB(:,:,3))。 imshow(log(abs(N)),[]),colorbar。 title(39。 Y 分量能量分布 39。 )。 subplot(2,3,4) I=imread(39。 lena139。 ,39。 bmp39。 )。 imshow(I)。 title(39。 灰度水印图像 39。 )。 subplot(2,3,5) M=dct2(I)。 imshow(log(abs(M)),[]),colorbar。 title(39。 水印能量分布 39。 )。 subplot(2,3,6) J=M(1:128,1:128)。 J(128:364,128:400)=0。 J=rot90(J)。 J=rot90(J)。 J(365:600,401:750)=0。 J=rot90(J)。 J=rot90(J)。 N=N+Q*J。 K=idct2(N)。 RGB(:,:,3)=K。 imshow(RGB)。 title(39。 加入水印后图像 39。 )。 %水印提取程序 subplot(2,3,3) RGB1=imread(39。 浙江台州学院 39。 ,39。 jpg39。 )。 N=dct2(RGB(:,:,3))。 M=dct2(RGB1(:,:,3))。 M=(NM)/Q。 B=idct2(M(236:365,350:401))。 Y=mat2gray(B)。 imshow(Y)。 title(39。 提取的水印图像 39。 ) 运行结果： ③ 水印全过程（经剪切检测水印） %水印加入程序 Q=input(39。 请输入放缩因子的值 (建议小于 1)： Q＝ 39。 ) subplot(3,3,1) RGB=imread(39。 MM39。 ,39。 jpg39。 )。 imshow(RGB)。 title(39。 原始图像 39。 )。 subplot(3,3,2) imshow(RGB(:,:,3))。 title(39。 B 分量 39。 )。 subplot(3,3,3) N=dct2(RGB(:,:,3))。 imshow(log(abs(N)),[]),colorbar。 title(39。 B 分量能量分布 39。 )。 subplot(3,3,4) I=imread(39。 lena139。 ,39。 bmp39。 )。 imshow(I)。 title(39。 灰度水印图像 39。 )。 subplot(3,3,5) M=dct2(I)。 imshow(log(abs(M)),[]),colorbar。 title(39。 水印能量分布 39。 )。 subplot(3,3,7) J=M(1:128,1:128)。 J(128:464,128:364)=0。 J=rot90(J)。 J=rot90(J)。 J(465:800,365:600)=0。 J=rot90(J)。 J=rot90(J)。 N=N+Q*J。 K=idct2(N)。 RGB(:,:,3)=K。 imshow(RGB)。 title(39。 加入水印后图像 39。 )。 subplot(3,3,8) I=imcrop(RGB,[1 1 598 798])。 imshow(I)。 subplot(3,3,9) %水印提取程序 subplot(3,3,6) RGB1=imread(39。 MM39。 ,39。 jpg39。 )。 J=RGB1(:,:,3)。 X=J(1:799,1:599)。 N=dct2(I(:,:,3))。 M=dct2(X)。 M=(NM)/Q。 B=idct2(M(337:464,237:364))。 Y=mat2gray(B)。 imshow(Y)。 title(39。 经放缩后提取的水印图像 39。 ) 运行结果： ④ 水印全过程（经空域压缩检测水印） 程序源代码 %水印加入程序 Q=input(39。 请输入放缩因子的值 (建议小于 1)： Q＝ 39。 ) P=input(39。 请输入您所希望的图像放缩系数值 (建议取值不要小于 ):P=39。 ) subplot(3,3,1) RGB=imread(39。 浙江台州学院 39。 ,39。 jpg39。 )。 imshow(RGB)。 title(39。 原始图像 39。 )。 subplot(3,3,2) imshow(RGB(:,:,3))。 title(39。 B 分量 39。 )。 subplot(3,3,3) N=dct2(RGB(:,:,3))。 imshow(log(abs(N)),[]),colorbar。 titl。