基于四元数PHFM的三图像零水印构造方法和检测方法技术

本发明专利技术公开了基于四元数PHFM的三图像零水印构造方法和检测方法，属于数字多媒体防伪和信息安全保护领域，解决的问题是如何实现同时对三幅同类的数字医学图像进行版权保护；其方法包括构造上述三幅医学图像的四元数极谐Fourier矩；基于选取的矩的幅值计算二维特征矩阵；构造感知哈希序列；构造Logistic混沌映射序列并得到混沌序列；生成零水印图像。其检测方法包括构造三幅待检测的医学图像的四元数极谐Fourier矩；基于选取的矩的幅值计算二维特征矩阵；构造感知哈希序列；构造Logistic混沌映射序列并得到混沌序列；提取水印图像。该方法可对三幅医学图像进行版权保护，且不会对原始医学图像进行修改。

【技术实现步骤摘要】
基于四元数PHFM的三图像零水印构造方法和检测方法
本专利技术涉及数字多媒体防伪和信息安全保护领域，具体地说是基于四元数PHFM的三图像零水印构造方法和检测方法。
技术介绍
当前医院医疗数字信息化建设进程加快，诸如计算机断层(CT)成像、超声成像(US)、核磁共振成像(NMRI)等数字医学影像作为医生判断患者健康状况的重要参考依据，这些图像在网络传输时很可能会遭受到被信息截获或被故意篡改等安全问题。增强数字医疗信息系统的健壮性，尤其有效保证医学图像的可靠性、可用性和保密性已成为目前亟待解决的问题。数字水印技术是有效解决该问题的技术手段之一。利用数字水印技术将具有标志性或有意义的数字信息嵌入到医学图像中，其特有的鲁棒性和安全性保证图像在经历数据信息处理后，仍能够可靠完整地提取出水印标志，实现其版权保护、完整性认证等作用。然而嵌入的水印信息会破坏医学图像的完整性，由于医学图像包含了一些重要的病理信息，任何较小的改动都会影响医生对病理的判断，故而设计一种无损的医学图像版权保护方式至关重要。在实际应用中可能需要对多幅同样的数字医学影像分别进行版权保护，例如，在CT图像中，身体同一部位通常会被采集成多幅图像，上述多幅图像均需进行版权保护，但现有的零水印方案主要是面向一幅图像进行版权保护，虽然可以重复用于多幅图像，但是在重复操作中会耗费大量时间、占用大量的存储空间。如何实现对多幅同样的数字医学影像进行版权保护，是需要解决的技术问题。
技术实现思路
本专利技术的技术任务是提供一种基于四元数PHFM的三图像零水印构造方法和检测方法，来解决如何同时实现对三幅同类医学图像进行版权保护的问题。本专利技术的技术任务是按以下方式实现的：基于四元数PHFM(英文全称为PolarHarmonicFourierMoments，中文翻译为极谐傅里叶矩)的三图像零水印构造方法，用于同时使用三幅医学图像构造零水印图像，三幅医学图像的大小均为N×N，水印图像的大小为P×Q，包括如下步骤：S100、计算上述三幅医学图像的四元数极谐Fourier矩；S200、基于密钥选取P×Q个四元数极谐Fourier矩，并基于选取的四元数极谐Fourier矩的幅值计算二维特征矩阵；S300、对二维特征矩阵逐列比较列向量中每个幅值和列向量幅值平均值的大小，构造感知哈希序列；S400、构造Logistic混沌映射序列并将其转换为二值序列，对二值序列和感知哈希序列进行异或操作得到混沌序列；S500、将混沌序列与水印图像进行异或操作，生成零水印图像。进一步的，步骤S100中，以三幅医学图像作为四元数的三个虚部，计算上述三幅医学图像的四元数极谐Fourier矩。进一步的，步骤S200包括如下分步骤：S210、根据密钥从上述四元数极谐Fourier矩中随机选取P×Q个矩；S220、计算上述选取的四元数极谐Fourier矩的幅值作为图像特征；S230、对选取的四元数极谐Fourier矩的幅值进行升维计算得到二维特征矩阵V，V＝(vi,j,0≤i≤P,0≤j≤Q)其中，vi,j表示二维特征矩阵V在坐标(i,j)处的值。进一步的，步骤S300包括如下分步骤：S310、对于二维特征矩阵逐列计算每列向量对应的列向量幅值平均值；S320、对于二维矩阵逐列进行如下计算：将当前列向量的每个幅值逐一与当前列向量的列向量幅值平均值进行大小比较，并生成序列元素；其中，当幅值大于或等于列向量幅值平均值时，序列元素记为1，当幅值小于列向量幅值平均值时，序列元素记为0；S330、将上述序列元素扫描成序列，即得到感知哈希序列H，H＝(hi,j,0＜i＜P,0＜j＜Q)其中，hi,j表示感知哈希序列H在坐标(i,j)处的值。进一步的，步骤S400包括如下分步骤：S410、通过初始值为x0的Logistic混沌映射生成长度为P×Q的混沌序列L；S420、将混沌序列L转换为二值序列L'，其中，Threshold为混沌序列L的均值；S430、对二值序列L'与哈希序列H进行异或操作生成混沌序列H'，基于四元数PHFM的三图像零水印检测方法，同时对三幅医学图像进行水印检测，三幅医学图像的大小均为N×N，包括如下步骤：L100、计算上述三幅待检测医学图像的四元数极谐Fourier矩；L200、基于密钥选取P×Q个四元数极谐Fourier矩，并基于选取的四元数极谐Fourier矩的幅值计算二维特征矩阵；L300、对二维特征矩阵逐列比较列向量中每个幅值和列向量幅值平均值的大小，构造感知哈希序列；L400、构造Logistic混沌映射序列并将其转换为二值序列，对二值序列和感知哈希序列进行异或操作得到混沌序列；L500、将混沌序列与通过上述任一项基于四元数PHFM的三图像零水印构造方法构造的零水印图像进行异或操作，提取水印图像。进一步的，步骤L100中，以三幅待检测医学图像作为四元数的三个虚部，构造上述三幅医学图像的四元数极谐Fourier矩。进一步的，步骤L200包括如下分步骤：L210、根据密钥从上述四元数极谐Fourier矩中随机选取P×Q个四元数极谐Fourier矩；L220、计算上述选取的四元数极谐Fourier矩的幅值作为图像特征；L230、对四元数极谐Fourier矩的幅值进行升维计算得到二维特征矩阵其中，表示二维特征矩阵在坐标(i,j)处的值。进一步的，步骤L300包括如下分步骤：L310、对于二维特征矩阵逐列计算每列向量对应的列向量幅值平均值；L320、对于二维矩阵逐列进行如下计算：将当前列向量的每个幅值逐一与当前列向量的列向量幅值平均值进行大小比较，并生成序列元素；其中，当幅值大于或等于列向量幅值平均值时，序列元素记为1，当幅值小于列向量幅值平均值时，序列元素记为0；L330、将上述序列元素扫描成序列，即得到感知哈希序列其中，表示感知哈希序列在坐标(i,j)处的值。进一步的，步骤L400包括如下分步骤：L410、通过初始值为x0的Logistic混沌映射生成长度为P×Q的混沌序列L420、将混沌序列转换为二值序列其中，Threshold为混沌序列的均值；L430、对二值序列与哈希序列进行异或操作生成混沌序列本专利技术的一种基于四元数PHFM的三图像零水印构造方法和检测方法具有以下优点：1、该方法中实现了同时对三幅同类的医学图像进行版权保护，有效提升了算法的效率，节省了存储空间；2、首次构造了三幅不同图像的四元数极谐Fourier矩；3、该方法中水印图像并没有嵌入到医学图像载体中，而是与载体相关联，确切的说是与图像特征向量相关联，没有对原始的医学图像进行任何修改，保证了医学图像的完整性。附图说明下面结合附图对本专利技术进一步说明。附图1为实施例1基于四元数PHFM的三图像零水印构造方法的流程框图；附图2为实施例2基于四元数PHFM的三图像零水印检测方法的流程框图。具体实施方式参照说明书附图和具体实施例对本专利技术的基于四元数PHFM的三图像零水印构造方法和检测方法作以下详细地说明。实施例1：如附图1所示，本专利技术的基于四元数PHFM的三图像零水印构造方法，用于同时使用三幅医学图像构造零水印图像，三幅医学图像大小均为N×N，水印图像W的大小为P×Q，包括如下步骤：S100、计算上述三幅医学图像的四元数极谐Fourier矩本文档来自技高网...

【技术保护点】
1.基于四元数PHFM的三图像零水印构造方法，其特征在于用于同时使用三幅医学图像构造零水印图像，三幅医学图像的大小均为N×N，水印图像的大小为P×Q，包括如下步骤：S100、计算上述三幅医学图像的四元数极谐Fourier矩；S200、基于密钥选取P×Q个四元数极谐Fourier矩，并基于选取的四元数极谐Fourier矩的幅值计算二维特征矩阵；S300、对二维特征矩阵逐列比较列向量中每个幅值和列向量幅值平均值的大小，构造感知哈希序列；S400、构造Logistic混沌映射序列并将其转换为二值序列，对二值序列和感知哈希序列进行异或操作得到混沌序列；S500、将混沌序列与水印图像进行异或操作，生成零水印图像。

【技术特征摘要】
1.基于四元数PHFM的三图像零水印构造方法，其特征在于用于同时使用三幅医学图像构造零水印图像，三幅医学图像的大小均为N×N，水印图像的大小为P×Q，包括如下步骤：S100、计算上述三幅医学图像的四元数极谐Fourier矩；S200、基于密钥选取P×Q个四元数极谐Fourier矩，并基于选取的四元数极谐Fourier矩的幅值计算二维特征矩阵；S300、对二维特征矩阵逐列比较列向量中每个幅值和列向量幅值平均值的大小，构造感知哈希序列；S400、构造Logistic混沌映射序列并将其转换为二值序列，对二值序列和感知哈希序列进行异或操作得到混沌序列；S500、将混沌序列与水印图像进行异或操作，生成零水印图像。2.根据权利要求1所述的基于四元数PHFM的三图像零水印构造方法，其特征在于步骤S100中，以三幅医学图像作为四元数的三个虚部，计算上述三幅医学图像的四元数极谐Fourier矩。3.根据权利要求2所述的基于四元数PHFM的三图像零水印构造方法，其特征在于步骤S200包括如下分步骤：S210、根据密钥从上述四元数极谐Fourier矩中随机选取P×Q个矩；S220、计算上述选取的四元数极谐Fourier矩的幅值作为图像特征；S230、对四元数极谐Fourier矩的幅值进行升维计算得到二维特征矩阵V，V＝(vi,j,0≤i≤P,0≤j≤Q)其中，vi,j表示二维特征矩阵V在坐标(i,j)处的值。4.根据权要求3所述的基于四元数PHFM的三图像零水印构造方法，其特征在于步骤S300包括如下分步骤：S310、对于二维特征矩阵逐列计算每列向量对应的列向量幅值平均值；S320、对于二维矩阵逐列进行如下计算：将当前列向量的每个幅值逐一与当前列向量的列向量幅值平均值进行大小比较，并生成序列元素；其中，当幅值大于或等于列向量幅值平均值时，序列元素记为1，当幅值小于列向量幅值平均值时，序列元素记为0；S330、将上述序列元素扫描成序列，即得到感知哈希序列H，H＝(hi,j,0＜i＜P,0＜j＜Q)其中，hi,j表示感知哈希序列H在坐标(i,j)处的值。5.根据权要求4所述的基于四元数PHFM的三图像零水印构造方法，其特征在于步骤S400包括如下分步骤：S410、通过初始值为x0的Logistic混沌映射生成长度为P×Q的混沌序列L；S420、将混沌序列L转换为二值序列L'，其中，Threshold为L的均值；S430、对二值序列L'与哈希序列H进行异或操作生成混沌序列H'，6....

【专利技术属性】
技术研发人员：王春鹏，夏之秋，马宾，李健，
申请(专利权)人：齐鲁工业大学，
类型：发明
国别省市：山东,37