[发明专利]用于RAW格式图像的CWT-QIM数字盲水印算法

权利要求书

1.用于RAW格式图像的CWT-QIM数字盲水印算法，其特征在于，包括水印的嵌入过程和水印的提取过程，具体步骤如下：

(1)水印的嵌入过程

步骤a.水印图像的生成

假设待嵌入RAW图像的大小为M×N，原始水印信息为一维向量W，根据光栅扫描顺序，原始水印信息W描述为一维向量的形式：

W＝{w₀,w₁,…,w_l-1}            式1

其中它的每个元素w_i满足w_i∈{-1，1}，l大小为m、n分别是DCT分块的块尺寸与水印嵌入的DT-CWT的变换层数；根据水印信息的长度利用随机数产生器对W生成一个长为l的随机数构成的一维向量，记为K；

K＝{k₀,k₁,…,k_l-1}             式2

将得到的l个水印信号w_i与密钥k_i实行线性调制，得到一个调制后的含水印的信号即：

其中，是由l个组成的一维向量；

步骤b.选择嵌入区域

将原始RAW图像插值提取大小为的B通道的图像矩阵，对B通道的图像矩阵使用n次DT-CWT变换，而每次DT-CWT变换会产生一个低频子带和6个不同方向的高频子带，6个不同方向分别为±15°、±45°和±75°，选择两个互相垂直、方向为±45°、大小为的高频子带作为嵌入的载体，划分两个子带的取模矩阵，对大小为的子带等比例的划分为S块大小为m×m且不重复的方块矩阵，其中S大小如式4所示：

步骤c.DCT变换

通过对取模后获得矩阵按照m×m大小的块连续划分整个对整个的每个块按式5、式6去分别计算它的DCT系数矩阵D_k；

采样提取D_k的直流系数d_k，具体如式7描述：

d_k＝D_k(1，1)         式7

其中k代表块的编号，将每块的直流系数d_k按式8合成矩阵此时，是大小为的直流系数矩阵，代表子带每个块最稳定不易修改的能量分布；

步骤d.基于矢量量化调制的嵌入

对两个子带分别进行分块DCT，得到两个直流矩阵与相比得到比例系数矩阵L，其数学描述如式9所示：

将坐标轴等间隔划分成两个区间集，区间间隔为Δ；根据水印值是0或1，动态调整待嵌入水印的系数值L_x，y，数学上的表示为式10、式11；

其中L′描述了修改后的L值，也就是水印嵌入后±45°方向的新比例矩阵；

将经过n次DT-CWT变换的分割为S块大小m×m的矩阵对块矩阵整块的放大倍得到矩阵；数学上的表示为式12：

从上到下、从左向右依次的修改S个矩阵，按照式13的排列修改原子带方向系数矩阵：

利用得到的系数矩阵再结合未修改的其他层系数矩阵B^L、通过逆DT-CWT变换得到B通道矩阵，至此得到含水印图像；如式14所示：

(2)水印的提取过程

信息接受方从接收到的含水印图像中插值提取出B通道系数矩阵，根据嵌入方共享的子带方向d、嵌入层数n、分块大小m、量化步长Δ，可以得到嵌入的水印图案，进而由接收方利用密钥矩阵K来准确恢复水印图像，K＝{k₁,k₂,…,k_q}；假设X为一幅待提取水印的RGB图像，尺寸为M×N，接受方根据已知的嵌入参数信息通过相应的逆操作进而获得水印图像其向量表示为具体实现步骤如下：

步骤a.将插值提取的B通道系数矩阵经过n层DT-CWT变换得到B^L、提取嵌入方向d的高频子带与作为水印检测区域，这里d取1、2；

步骤b.分别对与进行大小为m×m的DCT分块，由式5得到相应DCT系数矩阵D₁、D₂，再经式7、式8从系数矩阵D₁、D₂提取稳定的直流系数矩阵

步骤c.根据式9计算比例系数矩阵L，记为l_i，i＝1，2...，q，L代表着q组变换域比例系数，然后，根据QIM逆变换式15得出位序列其中符号函数mod(x,y)定义为x对y的取模函数，round(x)定义为找寻对向量x中每个元素朝最近的方向取整数部分，并返回与x同维的向量；

步骤d.最后用与嵌入时相同密钥K解析水印得到原始水印W₀。

2.根据权利要求1所述的用于RAW格式图像的CWT-QIM数字盲水印算法，其特征在于，在光栅扫描原始水印信息的过程中采用纠错编码技术。