没有合适的资源?快使用搜索试试~ 我知道了~
基于matlab图像数字水印的方案文档.doc
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 140 浏览量
2023-10-22
23:13:48
上传
评论
收藏 1.28MB DOC 举报
温馨提示
试读
16页
在数据库中存储在国际互联网上传输的水印图像一般会被压缩,有时达到很高的压缩比。因此,数字水印算法所面临的第一个考验就是压缩。JPEG和EZW(Embedded Zero-Tree Wavelet)压缩是最常见的两种压缩方法。JPEG是基于离散余弦变换域的压缩方法,而EZW是基于小波变换域的压缩方法。前人的研究证明采用与压缩算法相同的变换域水印方法,对于压缩的稳健性较强。因此,我研究图像文件水印算法主要集中在变换域算法及利用人眼视觉特性上。
资源推荐
资源详情
资源评论
图像数字水印的方案
3.1 图像数字水印的技术方案
在数据库中存储在国际互联网上传输的水印图像一般会被压缩,有时达到很高的压
缩比。因此,数字水印算法所面临的第一个考验就是压缩。JPEG 和 EZW(Embedded
Zero-Tree Wavelet)压缩是最常见的两种压缩方法。JPEG 是基于离散余弦变换域的压
缩方法,而 EZW 是基于小波变换域的压缩方法。前人的研究证明采用与压缩算法相同
的变换域水印方法,对于压缩的稳健性较强。因此,我研究图像文件水印算法主要集中
在变换域算法及利用人眼视觉特性上。
数字水印的嵌入要求即要考虑视觉透明性,又要保证嵌入水印后图像的稳健性,这
两个方面存在着矛盾。保证视觉透明性,就要将水印嵌入到人眼不敏感区,也就是嵌入
到图像的高频分量中。而多数图像处理方法对于图像高频部分的损坏程度较高,如有损
压缩、高频滤波等。水印很容易在经历图像处理的过程中丢失。这样,则无法保证图像
数字水印的稳健性。如果要获得很好的稳健性,数字水印应加在人眼敏感的低频部分,
图像的大部分能量集中在低频部分,如果对于低频部分进行处理,水印固然会失去,而
图像也没有了利用价值,然而,水印的嵌入会对图像的质量有非常大的影响,这又无法
保证视觉透明性。
数字水印算法的实现基本分为三个部分:宿主图像的变换,水印的嵌入和水印的检
测,分别描述如下。
3.2 基于 DCT 域的图像数字水印技术
离散余弦变换(Discrete Cosine Transform)属于正交变换图像编码方法中的一种。
正交变换图像编码始于 1968 年。当时安德鲁斯(Andrews)等人发现大多数自然图像
的高频分量相对幅度较低,可完全舍弃或者只用少数码字编码,提出不对图像本身编码,
只对其二维傅立叶(DFT)系数进行编码和传输。但 DFT 是一种正交变换,运算量很
大,常常使实时处理发生困难,第二年他们就用 Walsh-Hadamard 变换(WHT)取代 DFT
可以使运算量明显减少,这是因为 WHT 变换只有加减法而无需乘法。但是更有意义的
是离散余弦变换和离散正旋变换的出现,它们具有快速算法,精确度高。其中最重要的
是 1974 年提出的 DCT,因为其变换矩阵的基向量很近似于托伯利兹矩阵的特征向量,
而托伯利兹矩阵又体现了人类语言及图像信号的相关性。因此,DCT 常常被认为是语
音与图像信号变换的准最佳变换。
图像是二维的,所以在研究时主要用到二维 DCT,以及二维 IDCT 来对图像进行
处理。
3.2.1 离散余弦变换(DCT)的定义
数字图像 X(m,n)是具有 M 行 N 列的一个矩阵。为了同时减弱或去除图像数据
相关性,可以运用二维 DCT,将图像从空间域转换到 DCT 变换域。
根据定义,二维离散余弦变换(DCT)定义如下:
N
ln
M
km
nmXlckc
MN
lkY
M
m
N
n
2
)12(
cos
2
)12(
cos),()()(
2
),(
1
0
1
0
��
��
�
� �
�
�
�
�
式中:m,k=0,1,…,M-1; n,
l
=0,1,…,N-1。
�
�
�
��
�
�
1,,2,11
021
)(
Mk
k
kc
�
其中函数
�
�
�
��
�
�
1,,2,11
021
)(
Nk
k
lc
�
二维逆离散余弦变换(IDCT)的定义如下:
� �
�
�
�
�
��
�
1
0
1
0
2
)12(
cos
2
)12(
cos),()()(
2
),(
M
K
N
L
N
ln
M
km
lkYlckc
MN
nmX
��
式中:m,k=0,1,…,M-1; n,
l
=0,1,…N-1。
3.2.2 离散余弦变换的特点
在基于 DCT 的变换编码中,图像是先经分块(8×8 或 16×16)后再经 DCT,这
种变换是局部的,只反映了图像某一部分的信息。当然也可以对整幅图像的特点,但是
运算速度比分块 DCT 要慢。图像经 DCT 后,得到的 DCT 图像有三个特点:
一是系数值全部集中到 0 值附近(从直方图统计的意义上),动态范围很小,这说
明用较小的量化比特数即可表示 DCT 系数;
二是 DCT 变换后图像能量集中在图像的低频部分,即 DCT 图像中不为零的系数大
部分集中在一起(左上角),因此编码效率很高。
三是没有保留原图像块的精细结构,从中反映不了原图像块的边缘、轮廓等信息,
这一特点是由 DCT 缺乏时局域性造成的。
如下左图 3—1 是原始图像经过 DCT 变换后的系数图像为图 3—2。两条线划分出
图像的低频、中频和高频分别所在的矩形区域。可以看出,图像 DCT 变换后大部分参
数接近于零,只有左上角的低频部分有较大的数值,中频部分参数值相对较小,而大部
分高频参数值非常小,接近于零。
MM
MMM
MM
MM
MM
MM
MMM
MMM
MMM
MMM
…
MMM
MMM
….
MMM
MMM
…
MMM
K
…
图 3—1:原图像 图 3—2:变换后的系数图像
3.2.3 离散余弦变换的数字水印算法
根据离散余弦变换后的参数性质,本文采用了以 ZigZag 方式重排变换域系数的方
法,选出中频分量,用数字水印序列对其进行非线性调制。水印检测时,待检测图像仍
按比方式选择变换域系数,与待水印进行相关运算,与阈值比较来判断是否所含水印。
离散余弦域的数字水印算法的具体实现分为三步:宿主图像的变换,数字水印的嵌
入,数字水印的检测。
3.2.3.1 宿主图像的 DCT 变换
对于 N×N 大小的 256 灰度级的宿主图像 I 进行 N×N 二维离散余弦变换(DCT)。
以 ZigZag 方 式 对 于 DCT 变 换 后 的 图 像 频 率 系 数 重 新 排 列 成 一 维 向 量 Y={y
1
,
y
2
,…y
N×N
}.
并取出序列中第 L+1 到 L+M 的中频系数部分,得到 Y
L
={ Y
L+1
, Y
L+2
,…, Y
L+M
}
3.2.3.2 数字水印的嵌入
假设数字水印 W 为一服从标准正态分布的随机实数序列,用数字序列表示为
W={W
1
,W
2
,…W
M
}。用 W 对 Y 序列中第 L+1 到 L+M 的中频系数部分的值进行
修改,按以下公式进行:
经过修改的系数序列 Y′ ={ Y
1
′, Y
2
′,... Y′
N×N
} 以 ZigZag 逆变换形式重组,
再进行 N×N DCT 逆变换,得到嵌有数字水印的图像 I′。
3.2.3.3 数字水印的检测
待检测的可能含有水印的图像 I"
。假设 I"未损失大量信息,可以近似认为 I"= I
′。在此假设下可以运用统计的方法来检测水印。
(1)待检水印域待检图像中频系数相关性的测定
同样对 I′进行 DCT 变换,以 ZigZag 方式将 DCT 系数排成一维向量 Y "= { Y
1
", Y
2
",... Y
N×N
"}。由于假设 I"=I′,则 Y"= Y ′。
取出 Y"(等于 Y′)中第 L+1 到 L+M 的中频系数部分 Y
L
"={ Y
L+1
", Y
L+2
"
’
,...
Y
L+M
"}。假设待检测的数字水印 X={X
1
,X
2
,... X
M
}为一符合标准正态分布的实数伪随
机序列。则可以通过待检水印与图像中频系数作相关运算来判断是否所加入了水印。只
有在待检水印为所加入的水印时,才能得到较大的相关值。否则相关值很小,接近于零。
��
��
�
�
���� ))()((
1
)(
1
),'cov(
2
1
XiWYXY
M
XY
M
XYZ
iiLiiL
M
i
iiLL
�
用符号 E 表示数学期望,得到:
'
' 2
, 0
i i
i i i i
y y i L i L m
y y y w L i L M
a a
ì
= £ ³ +
í
= + < £ + >
î
或
(2)阈值的确定
根据中心极限定理,参照水印匹配与不匹配两种情况得到阈值为
2/)(
2
Lz
YET
�
�
。
由于原始图像难以得到,因此从实用性出发,阈值定义为:
�
�
�
��
M
i
iLLZ
Y
M
YET
1
2
2/)(
�
�
综上所述,满足
5.0
2
2
�
�
�
Z
Z
T
TZ 或
时,则表明检测到匹配水印。否则,未检测到
匹配的水印。
3.3 MATLAB 工具简介
3.3.1. 简介
Matlab 是当前在国内外十分流行的工程设计和系统仿真软件包。它是 MathWorks
公司于 1982 年推出的一套高性能的数值计算和可视化软件,它集数值分析、矩阵运算、
信号处理和图形显示于一体,构成了一人方便的、界面友好的用户环境。
Matlab 的推出得到了各个领域专家、学者的广泛关注,其强大的扩展功能为各个领域
的应用提供了基础。由各个专家学者相继推出了 MATLAB 工具箱,其中的信号处理
(signal processing)、控制系统(control system)、神经网络(neural network)、图像处理(image
processing)、鲁棒控制(robust control)、非线性系统控制设计(nonlinear system control
design)、系统辨识(system identification)、最优化(optimization)、模糊逻辑(fuzzy logic)、
小波(wavelet)、通信(communication)、统计(statistics)等工具箱,这些工具箱给各个领域
的研究和工程应用提供了有力的工具,借助于这些“巨人肩上的工具”,各个层次的研究
人员可直观、方便地进行分析、计算及设计工作,从而大大地节省了时间。
3.3.2. MATLAB 研究数字水印的优点
① 集成了 DCT、DWT 等函数有丰富的小波函数和处理函数,这不仅方便了研究人员,
而且使源程序简洁明了、易实现。
② 强大的数学运算功能。能够方便、高效地实现音频、视频中的大量矩阵运算。
③ 提供了图像处理工具箱、小波分析工具箱、数字信号处理工具箱。用来编制跨数字
图像处理技术、数字信号处理等多学科的数字水印技术是非常好的选择。
④ MATLAB 与目前最强大的编程工具——Visual C++具有良好的接口。
2
( ) ( )
( ) 0
( ) 0
E z E Y X W
E z X W
E z
a
ì
= =
ï
= ¹
í
ï
=
î
没有水印存在
剩余15页未读,继续阅读
资源评论
等天晴i
- 粉丝: 3434
- 资源: 10万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功