发明名称 彩色图像超复数空间的自适应水印算法
摘要 本发明属于图像信号处理技术领域,具体为一种彩色图像超复数空间的自适应水印算法。它首先对彩色载体图像进行快速超复数傅氏变换,在超复数频谱实部选择合适频段的基础上,再利用人类视觉系统对彩色载体图像的纹理、边缘和亮度的掩蔽特性,对选择的频段赋予不同的掩蔽强度而嵌入水印,从而在超复数频域内实现了一种彩色载体图像自适应的水印算法。实验结果表明,通过彩色图像的自适应掩蔽,大大提高了超复数频域水印算法的不易感知性和鲁棒性,且抗攻击性能也优于无自适应掩蔽的超复数频域水印算法;和现有文献的比较结果验证了本发明的这些优越性。
申请公布号 CN101246588B 申请公布日期 2011.04.13
申请号 CN200810034890.2 申请日期 2008.03.20
申请人 复旦大学 发明人 江淑红;张建秋;胡波
分类号 G06T1/00(2006.01)I 主分类号 G06T1/00(2006.01)I
代理机构 上海正旦专利代理有限公司 31200 代理人 陆飞;盛志范
主权项 1.一种彩色图像超复数空间的自适应水印算法,其特征在于嵌入数字水印的具体步骤如下:(1)首先对彩色载体图像进行超复数建模:把彩色载体图像f<sub>0</sub>的每一像素的RGB模型分别表示为纯超复数形式:f<sub>0</sub>(x,y)=R(x,y)i+G(x,y)j+B(x,y)k,其中,(x,y)为彩色载体图像中像素的坐标,即像素所在的矩阵行和列的位置,R、G、B分别表示红、绿、蓝三种颜色;(2)采用密钥Key控制的伪随机数发生器,伪随机地选择M<sub>m</sub>×N<sub>m</sub>个整数,取值范围是1-14,作为数字水印的嵌入位置;根据彩色载体图像分块规则,把彩色载体图像f<sub>0</sub>分成8×8的单位小块;(3)计算彩色载体图像的纹理掩蔽因子:计算红、绿、蓝每一色的直方图方差,再取平均得到平均直方图方差,对平均直方图方差进行归一化得到σ<sup>2</sup><sub>u</sub>,用σ<sub>u</sub>的平方根作为每个单位小块的纹理掩蔽的描述子,并把它规整为0~5级,即纹理掩蔽因子为:<img file="FSB00000368571600011.GIF" wi="435" he="69" />其中,round(x)表示对x进行四舍五入;(4)计算彩色载体图像的边缘掩蔽因子:根据超复数彩色边缘检测算法或经典的canny边缘算子,计算彩色载体图像的边缘信息,设P<sub>E</sub>表示8×8单位小块中边缘点的和,则归一化后的5级边缘掩蔽因子为:M<sub>E</sub>=round(5P<sub>E</sub>/max(P<sub>E</sub>));(5)计算彩色载体图像的亮度掩蔽因子:根据彩色载体图像的亮度f<sub>A</sub>=<img file="FSB00000368571600012.GIF" wi="667" he="75" />得到每个8×8单位小块的平均亮度值为:<img file="FSB00000368571600013.GIF" wi="344" he="117" />设<img file="FSB00000368571600014.GIF" wi="81" he="63" />表示彩色载体图像的中等亮度,则亮度掩蔽的描述子由<img file="FSB00000368571600015.GIF" wi="321" he="63" />表示,归一化后的5级亮度掩蔽因子为:M<sub>I</sub>=round(5P<sub>I</sub>/max(P<sub>I</sub>));(6)综合考虑上述彩色载体图像的纹理、边缘和亮度掩蔽特性,彩色载体图像的自适应数字水印掩蔽因子为:J<sub>I</sub>=M<sub>T</sub>-M<sub>E</sub>+M<sub>I</sub>,然后去除M<sub>T</sub>-M<sub>E</sub>+M<sub>I</sub>的最大值和最小值,再把得到的结果规整为0~5级,数值加1,获得最终的1~6级自适应掩蔽因子J<sub>I</sub>;(7)把彩色载体图像f<sub>0</sub>分成8×8的单位小块,对每一个8×8单位小块按照下式进行超复数傅氏变换:F<sup>R</sup>(u,v)=i(real(R<sub>RFT</sub>)+μ·imag(R<sub>RFT</sub>))      +j(real(G<sub>RFT</sub>)+μ·imag(G<sub>RFT</sub>))      +k(real(B<sub>RFT</sub>)+μ·imag(B<sub>RFT</sub>)) 式中real(p)表示取复数p的实部,imag(p)表示取复数p的虚部;p<sub>RFT</sub>表示p的实数傅氏变换;μ为单位虚向量,取强度图像矢量<img file="FSB00000368571600021.GIF" wi="380" he="65" />(8)按照“交叉冗余嵌入”规则,在每个单位小块的超复数傅氏变换F<sup>R</sup>(u,v)的实数部分A(u,v)中的对应嵌入位置,用量化索引调制方法嵌入水印数据,同时修改水印数据对称系数的值;每个单位小块的量化单位Δ为该单位小块的1~6级的自适应掩蔽因子J<sub>I</sub>和最小量化单位Δ<sub>0</sub>的乘积,其中,(u,v)是嵌入水印在频域的坐标;(9)设嵌入水印的频域载体图像为F<sup>R</sup>(u,v)=A(u,v)+iC(u,v)+jD(u,v)+kE(u,v),把所有的嵌入水印的单位小块,按照下式进行超复数傅氏逆变换,得到最终的时空域内的含水印载体图像f<sub>wm</sub>:f<sub>wm</sub>(x,y)=(real(A<sub>IRFT</sub>)+μ·imag(A<sub>IRFT</sub>))         +i(real(C<sub>IRFT</sub>)+μ·imag(C<sub>IRFT</sub>))         +j(real(D<sub>IRFT</sub>)+μ·imag(D<sub>IRFT</sub>))         +k(real(E<sub>IRFT</sub>)+μ·imag(E<sub>IRFT</sub>)) 。
地址 200433 上海市邯郸路220号