基于数字全息和傅里叶梅林变换的数字水印方法技术

本发明专利技术公开了一种基于数字全息和傅里叶梅林变换的数字水印嵌入方法，包括：将二值水印图像通过共轭对称延拓和离散傅里叶变换(DFT)，得到二维离散实值函数F(μ,ν)，将该实值函数映射至0～255生成数字水印全息图；在载体图像I离散傅里叶变换的中频区域中嵌入数字水印全息图的矩形块，在幅度谱的中心对称的相应块嵌入相同的矩形块；将新的幅度谱和原始相位谱进行逆离散傅里叶变换，得到含水印图像I

【技术实现步骤摘要】
基于数字全息和傅里叶梅林变换的数字水印方法
本专利技术属于数字水印
，具体地涉及一种基于数字全息和傅里叶梅林变换的抗几何攻击的数字水印嵌入和提取方法，尤其适用于打印扫描过程对水印图像的攻击。
技术介绍
在数字水印
，常用的攻击包括噪声、滤波、剪切、尺度变化、IBM(Internationalbusinessmachines)攻击、有损压缩编码及几何攻击，其中几何攻击包括对数字图像的旋转、缩放和平移变换，这些攻击方法的直接后果是破坏了数字水印检测器与数字水印的同步性，从而使检测器检测不到图像中的数字水印信号。而现有技术中普通数字图像水印技术对抗攻击鲁棒性差以及水印信息容量小，不能很好地抵抗常用的攻击手段。几何变换后水印嵌入的初始位置发生了空间上的变化，而现有大多数算法都会从水印嵌入的原始位置中提取水印，所以水印的提取在几何攻击后往往失效，因此在数字水印的鲁棒性研究中抵抗几何攻击相对抵抗其他攻击手段而言更具有挑战性。傅里叶-梅林变换是基于几何不变域的典型水印算法。O’Ruanaidh等首次利用傅里叶-梅林变换(Fourier-Mellin)，在图像中寻找平移变换、尺寸变换和旋转变换不变域。傅里叶-梅林变换过程是将数极坐标变换(LPT)和离散傅里叶变换相结合，利用Fourier变换幅值的空间平移不变性来抵抗平移攻击；再利用对数极坐标系中能将笛卡尔坐标系中旋转和缩放变化转化为平移的性质。D.Zheng提出将水印嵌入在图像的中频区域,并给出了为避免水印提取时穷尽搜索而采用的计算水印嵌入点的相位关联性算法，但Zheng算法在水印嵌入时需要选择指定的嵌入点，给嵌入过程带来了更多的计算量。
技术实现思路
为了解决上述技术问题，本专利技术目的是：提供了一种基于数字全息和傅里叶梅林变换的数字水印嵌入和提取方法，既简化了嵌入方法，又提高嵌入质量，该方法在鲁棒性和不可见性之间能达到很好的平衡。本专利技术的技术方案是：一种基于数字全息和傅里叶梅林变换的数字水印嵌入方法，包括以下步骤：S01：将二值水印图像通过共轭对称延拓和离散傅里叶变换(DFT)，得到二维离散实值函数F(μ,ν)，将该实值函数映射至0～255生成数字水印全息图；S02：在载体图像I离散傅里叶变换的中频区域中嵌入数字水印全息图的矩形块，在幅度谱的中心对称的相应块嵌入相同的矩形块；S03：将新的幅度谱和原始相位谱进行逆离散傅里叶变换，得到含水印图像IW。优选的，所述步骤S01中的二维离散实值函数F(μ,ν)由下列步骤得到：对二值水印图像w0(x,y)进行调制，调制后的水印图像为式中是由高斯随机数决定，通过调制的水印图像经过傅里叶变换，同参考光相干涉，将干涉产生的强度分布场作为嵌入的傅里叶变换全息图；设物光波f0(m,n)可用幅度分布函数A1(m,n)和相位分布函数表示，其中A1(m,n)表示图像信息，其公式为：其中，m＝1,2,…M-1；n＝1,2,…N/2-1，由物光波f0(m,n)共轭对称延拓后得到的复振幅分布为：式中，上标“*”表示复共轭，对μ为缩放比例，进行二维离散傅里叶变换为：将式(1)、(2)代入式(3)，可得优选的，所述步骤S02具体为：对I进行离散傅里叶变换，得到幅度谱FA(u,v)和相位谱FP；对幅度谱FA(u,v)进行对数极坐标变换(LPT)得到FA(r,θ)；将数字水印全息图的矩形块对称嵌入到FA(r,θ)的中频区域，嵌入公式为：当RL＜r＜RH时，式中：RL为极坐标系中最小半径；RH为极坐标系中最大半径；RH-RL为水印图的宽度；α为嵌入强度，W为水印信号；对进行逆对数极坐标变换(ILPT)得到优选的，所述嵌入强度α为2％。本专利技术还公开了一种基于数字全息和傅里叶梅林变换的数字水印提取方法，包括以下步骤：S01：对载体图像I和含水印图像IW分别作傅里叶变换和对数极坐标变换，得到对数极坐标下的幅度谱分别为FA(r,θ)和S02：提取水印信号W，当RL＜r＜RH时，与现有技术相比，本专利技术的优点是：将水印以全息图的方式隐藏在载体图像离散傅里叶变换(discreteFouriertransform,DFT)中频幅度谱中，在中心对称的相应块做相同的嵌入。该方法利用全息图的不可撕毁性，并巧用log-polar变换将在DFT域幅度谱的中频圆环转化为LPM域的矩形块，只须嵌入半个圆环即可，这样既简化了嵌入方法，又提高嵌入质量。仿真实验表明，对加水印图像进行几何攻击(如：剪切、缩放、旋转)、高斯滤波和JPEG压缩后，该方法在鲁棒性和不可见性之间能达到很好的平衡。附图说明下面结合附图及实施例对本专利技术作进一步描述：图1是本专利技术生成数字水印全息图的流程图；图2(a)是二值水印图像；图2(b)是数字水印全息图；图3是本专利技术水印嵌入方法的流程图；图4(a)为水印图嵌入中频区域的示意图；图4(b)是进行逆对数极坐标变换得到的示意图；图5为本专利技术水印提取方法的流程图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明了，下面结合具体实施方式并参照附图，对本专利技术进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本专利技术的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本专利技术的概念。实施例：本专利技术包括数字水印全息图的生成以及水印的嵌入与提取。水印的生成如图1所示，设二值水印图像为w0(x,y)，为取得较平稳的傅里叶谱，先对水印图像进行调制，调制后的水印图像表示为式中是由高斯随机数决定的。通过调制的水印图像经过傅里叶变换，同参考光相干涉，干涉产生的强度分布场就是用来嵌入的傅里叶变换全息图。设物光波f0(m,n)可用幅度分布函数A1(m,n)和相位分布函数表示，其中A1(m,n)表示图像信息，其公式为：其中，m＝1,2,…M-1；n＝1,2,…N/2-1。由物光波f0(m,n)共轭对称延拓后得到的复振幅分布为：式中，上标表示复共轭。式(2)的二维离散Fourier变换为将式(1)、(2)代入式(3)，可得因此，一个二维离散物光波复振幅f0(m,n)经共轭对称延拓和DFT变换后可得到二维离散实值函数F(μ,ν)，将此实值函数映射至0～255即可可生成数字全息图。由式(4)可见，全息图中的任意一点都不包含原始物光波的振幅和相位信息，具有不可撕毁性。图2(a)为128X128原始水印图像，图2(b)为二值化得到的计算数字全息图。设载体图像为512pixel×512pixelLena图像I，水印图像为128pixel×128pixel的全息图像，则水印的嵌入算法步骤如图3所示：步骤1：对I做离散傅里叶变换，得到幅度谱FA和相位谱FP。步骤2：对幅度谱FA(u,v)进行对数极坐标变换(LPT)得到FA(r,θ)。由LPT的性质可知：FA(r,θ)分辨率u＝ercosθ,v＝ersinθ0≤θ＜2πr∈R2(5)步骤3：考虑到离散傅里叶变换幅度谱的中心对称性，将水印图对称嵌入到FA(r,θ)的中频区域，见图4(a)，嵌入公式为：当RL＜r＜RH时，式中：RL为极坐标系中最小半径；RH为极坐标系中最大半径；RH-RL为水印图的宽度；α为嵌入强度，α越大水印的鲁棒性越好，但不可见性会降低。对于数字图像水印而言，几乎很少技术能做到嵌入水印完全不可见性，本文档来自技高网...

【技术保护点】
一种基于数字全息和傅里叶梅林变换的数字水印嵌入方法，其特征在于，包括以下步骤：S01：将二值水印图像通过共轭对称延拓和离散傅里叶变换(DFT)，得到二维离散实值函数F(μ,ν)，将该实值函数映射至0～255生成数字水印全息图；S02：在载体图像I离散傅里叶变换的中频区域中嵌入数字水印全息图的矩形块，在幅度谱的中心对称的相应块嵌入相同的矩形块；S03：将新的幅度谱和原始相位谱进行逆离散傅里叶变换，得到含水印图像I

【技术特征摘要】
1.一种基于数字全息和傅里叶梅林变换的数字水印嵌入方法，其特征在于，包括以下步骤：S01：将二值水印图像通过共轭对称延拓和离散傅里叶变换(DFT)，得到二维离散实值函数F(μ,ν)，将该实值函数映射至0～255生成数字水印全息图；S02：在载体图像I离散傅里叶变换的中频区域中嵌入数字水印全息图的矩形块，在幅度谱的中心对称的相应块嵌入相同的矩形块；S03：将新的幅度谱和原始相位谱进行逆离散傅里叶变换，得到含水印图像IW。2.根据权利要求1所述的基于数字全息和傅里叶梅林变换的数字水印嵌入方法，其特征在于，所述步骤S01中的二维离散实值函数F(μ,ν)由下列步骤得到：对二值水印图像w0(x,y)进行调制，调制后的水印图像为式中是由高斯随机数决定，通过调制的水印图像经过傅里叶变换，同参考光相干涉，将干涉产生的强度分布场作为嵌入的傅里叶变换全息图；设物光波f0(m,n)可用幅度分布函数A1(m,n)和相位分布函数表示，其中A1(m,n)表示图像信息，其公式为：其中，m＝1,2,…M-1；n＝1,2,…N/2-1，由物光波f0(m,n)共轭对称延拓后得到的复振幅分布为：式中，上标“*”表示复共轭，对μ为缩放比例，进行二维离散傅里叶变换为：

【专利技术属性】
技术研发人员：陈维霞，龚冬梅，顾济华，陈大庆，顾涵，
申请(专利权)人：苏州大学，常熟理工学院，
类型：发明
国别省市：江苏,32