[发明专利]一种基于阈值分割和直方图均衡的可逆水印方法

权利要求书

1.一种基于阈值分割和直方图均衡的可逆水印方法，其特征在于，包括：

A、水印嵌入阶段：使用高斯滤波器对载体医疗图像进行平滑处理，抑制图像噪声；使用阈值分割算法对处理后的图像的梯度幅值矩阵进行处理，获取图像前景区的轮廓；对横纵双向扫描获得的两个二值矩阵进行与运算，进而分割图像的前景区和背景区；将前景区作为水印嵌入区域，使用直方图移位技术，对前景区的直方图中像素值介于峰值点和零值点之间的直方柱进行移位，实现直方图均衡，在移位的同时，进行水印循环嵌入，为了在之后实现信息提取和载体无失真恢复，少量的附加信息随水印一同嵌入，随着嵌入循环次数的增加，每次循环中嵌入的附加信息逐渐增长，当图像像素为峰值点的数量小于附加信息的长度时，嵌入循环结束；

B、水印提取、载体无失真恢复阶段：利用峰值点和零值点从载密图像的前景区提取信息，利用提取到的信息进行图像恢复，提取到的信息除去附加信息便是本轮所提取的水印，进行多次循环提取，当某一次提取信息的附加信息显示为最后一轮，提取终止。

2.根据权利要求1所述的一种基于阈值分割和直方图均衡的可逆水印方法，其特征在于，步骤A进一步包括以下步骤：

A1、构造一个模板大小为m×n的高斯滤波器，所述m，n为大于零的正整数，采用高斯滤波器对图像进行平滑降噪处理，并计算处理后图像的梯度幅值矩阵；

A2、使用阈值分割算法对梯度幅值矩阵进行处理，得到一个与原图像等大小的二值矩阵BW；

A3、构造与原图像等大小的二值矩阵BW1和BW2，横向扫描BW的每一行，获取每一行中值为‘1’的最左侧的坐标left和最右侧的坐标right，将BW1对应行处于left和right之间的值设置为‘1’，两侧的值设置为‘0’，纵向扫描BW的每一列，获取每一列中值为‘1’的最上方的坐标top和最下方的坐标bottom，将BW2对应列处于top和bottom之间的值设置为‘1’，令两侧的值为‘0’；

A4、构造与原图像等大小的二值矩阵bw，同步扫描BW1和BW2，若两个矩阵中同一位置处的值都是‘1’，则令bw中对应坐标下的值为‘1’，否则设置为‘0’，载体图像中对应于bw矩阵中值为‘1’的位置的像素为前景区，其他像素为背景区；

A5、由载体图像的前景区像素构建直方图，获取直方图中的峰值点Mx和零值点Mn，并依据所选峰值点Mx和零值点Mn确定本轮嵌入过程中处于峰值点Mx和零值点Mn之间的直方柱的移动方向，当峰值点小于零值点时，直方柱向右移动，当峰值点大于零值点时，直方柱向左移动；

A6、横向扫描图像，使用标记序列S对像素值p等于零值点Mn或其相邻像素值Mn-1或Mn+1的像素进行记录，若直方柱向右移动，当像素值p为Mn时，向S中添加标记‘0’，当像素值p为Mn-1时，向S中添加标记‘1’，反之，若直方柱向左移动，当像素值p为Mn时，向S中添加标记‘0’，当像素值p为Mn+1时，向S中添加标记‘1’；

A7、将上一轮嵌入过程用到的峰值点Mx和零值点Mn分别转换成8位二进制的表示形式，并将其组合成16位的xn，首轮xn以16个‘0’填充，最后一轮的峰值点和零值点转换成二进制的表示后，使用LSB替换的方式嵌入本轮载体图像最后一行的前16个像素的LSB，并将原像素的LSB连接在附加信息的xn之后；

A8、每轮嵌入先进行附加信息的嵌入，再进行水印的嵌入，按照横向扫描的顺序同步扫描载体图像I和二值矩阵bw的每一行，若直方柱向右移动，对于前景区像素，当其像素值p大于Mx且小于Mn时，将其值修改为p+1，当像素值p小于Mx或大于等于Mn时，其值不变，当像素值p等于Mx时，将其值修改为p+m，m为当前要嵌入的一个信息比特，若直方柱向左移动，对于前景区像素，当其像素值p大于Mn且小于Mx时，将其值修改为p-1，当像素值p小于等于Mn或大于Mx时，其值不变，当像素值p等于Mx时，将其值修改为p-m，m为当前要嵌入的一个信息比特；

A9、重复A5～A8，直到某一轮循环可嵌入的比特数量，即图像像素为峰值点的个数，小于附加信息的比特数时，该轮的上一轮即为嵌入的最后一轮循环，嵌入终止。

3.根据权利要求1所述的一种基于阈值分割和直方图均衡的可逆水印方法，其特征在于，步骤B进一步包括以下步骤：

B1、首轮提取，获取载密图像最后一行前16个像素的LSB，每8位一组，分别构成该轮峰值点Mx和零值点Mn；

B2、按照横向扫描的顺序同步扫描载密图像I和二值矩阵bw的每一行，从载密图像的前景区提取信息，若峰值点Mx小于零值点Mn，当像素值p等于Mx时，提取信息比特‘0’，当像素值p等于Mx+1时，提取信息比特‘1’，若峰值点Mx大于零值点Mn，当像素值p等于Mx时，提取信息比特‘0’，当像素值p等于Mx-1时，提取信息比特‘1’；

B3、每一轮提取的信息的前16位，每8位一组，分别构成下一轮提取的峰值点Mx和零值点Mn，若该轮为首轮提取，则用提取到的信息的第17～32位以LSB替换方式还原载密图像最后一行的前16个像素的LSB；

B4、余下的信息比特R由标记序列S和水印比特构成，即首轮提取S时从所提取的信息比特的第33位开始，其余轮次提取S时从所提取的信息比特的第17位开始；

B5、按照横向扫描的顺序同步扫描当前载密图像I和二值矩阵bw的每一行，还原前景区像素，若峰值点Mx小于零值点Mn，当像素值p'大于Mx且小于Mn时，将其值修改为p′-1，当像素值p'等于Mn且当前S标记为‘1”时，将其值修改为p′-1，当像素值p'等于Mn且当前S标记为‘0’时，其值不变，当像素值p'小于等于Mx或大于Mn时，其值不变，若峰值点Mx大于零值点Mn，当像素值p'大于Mn且小于Mx时，将其值修改为p′+1，当像素值p'等于Mn且当前S标记为‘1’时，将其值修改为p′+1，当像素值p'等于Mn且当前S标记为‘0’时，其值不变，当像素值p'小于Mn或大于等于Mx时，其值不变；

B6、除去用于像素还原的标记序列S的比特，R中剩余的比特为该轮提取到的水印比特；

B7、重复B2～B6，直到某一次提取的信息的前16位全部为‘0’，则该轮为提取的最后一轮，提取终止，将提取的水印信息输出。