会员体验
专利管家(专利管理)
工作空间(专利管理)
风险监控(情报监控)
数据分析(专利分析)
侵权分析(诉讼无效)
联系我们
交流群
官方交流:
QQ群: 891211   
微信请扫码    >>>
现在联系顾问~
首页 / 专利库 / 电脑安全 / 密码学 / 一种基于DFT可抗几何攻击的医学图像多重水印方法

一种基于DFT可抗几何攻击的医学图像多重水印方法

阅读:1033发布:2021-03-02

IPRDB可以提供一种基于DFT可抗几何攻击的医学图像多重水印方法专利检索,专利查询,专利分析的服务。并且本发明涉及一种基于DFT可抗几何攻击的医学图像多重水印方法,是先进行多水印的嵌入,包括:(1)对原始医学图像进行全图DFT变换,在变换系数中提取该图的一个重要视觉特征的向量;(2)利用该特征向量和要嵌入的多个水印通过Hash函数得到对应的二值逻辑序列,并将其存于第三方;然后进行多水印提取,包括:(3)对待测医学图像进行全图DFT变换,找到待测图像的一个视觉特征向量;(4)利用Hash函数性质和存于第三方的二值逻辑序列来提取水印。本发明将医学图像的视觉特征向量、加密技术和第三方的概念有机结合起来,解决了多重水印嵌入的问题,有较强的抗几何攻击能力和抗常规攻击能力,以保护医学图像的版权和病患信息的隐秘性。,下面是一种基于DFT可抗几何攻击的医学图像多重水印方法专利的具体信息内容。

1.一种基于DFT可抗几何攻击的医学图像多重水印方法,其特征在于:基于全局DFT及抗几何攻击的特征向量的提取,并将水印技术、密码学中的Hash函数特性和“第三方”概念有机结合起来,实现了在医学图像中多重数字水印的嵌入,该方法共分两个部分,共计四个步骤:第一部分是多重水印嵌入:通过对多重水印的嵌入操作,得到相应的二值逻辑序列kKey(j);

1)对原始医学图像进行全局DFT,在变换系数中,利用低中频系数的符号序列来得到该图的抗几何攻击的向量V(j);

k

2)利用Hash函数和要嵌入的多重水印W(j),k=0,1,2,...,n;得到二值逻辑序列kKey(j),

k k

保存Key(j),下面提取水印时要用到,通过把Key(j)作为密钥向第三方申请,以获得对原始医学图像的所有权和使用权;

k

第二部分是多重水印提取:通过二值逻辑序列Key(j)和待测医学图像的抗几何攻击k的特征向量V’(j),提取出多重水印W’(j);

3)对待测医学图像进行全局DFT;在变换系数中,根据低中频系数的符号提取出待测图像的视觉特征向量V’(j);

k

4)利 用 Hash 函 数 性 质,和 存 在 第 三 方 的 Key(j),提 取 出 水 印,k k将W(j)和W’(j)进行归一化相关系数计算,来确定医学图像的所有权。

说明书全文

一种基于DFT可抗几何攻击的医学图像多重水印方法

技术领域

[0001] 本发明属于多媒体信号处理领域,涉及一种基于DFT变换和图像视觉特征的医学图像多重数字水印技术,具体是一种基于DFT可抗几何攻击的医学图像多重水印方法。

背景技术

[0002] 近年来,随着计算机科学技术和多媒体通信技术的迅速发展,远程医疗日益普及,但在互联网上进行医学图像传输时,病人的个人信息容易泄漏,利用数字水印的不可见性和鲁棒性可以较好的解决这个问题,即把病人的个人信息作为数字水印嵌入在医学图像中。
[0003] 目前对医学图像数字水印领域的研究主要集中在空间域和变换域(DFT、DFT和DWT)两个方面,它们分别通过改变空间域的某些象素的灰度或变换域的一些系数的值来嵌入水印。
[0004] 鉴于对医学图像病灶区保护的特殊性要求,一般的医学数字水印方法常选择将水印信息嵌入到图像的非感兴趣区域(Region of Non-interest,RONI)。医学图像中的感兴趣区域ROI(Region of Interest)指的是那些包含重要病理特征或诊疗信息的病灶区,若在该区域嵌入水印,则有可能造成错误的诊断。但往往人们在寻找ROI时,要花费很长的时间与精力,并且一旦选择有误,则有可能干扰医生的诊断。
[0005] 在医学数字水印研究领域,至今为止抗击几何攻击和多重水印的嵌入仍是一个比较难以解决的课题,至于同时能有效抵抗常规攻击和几何攻击的多重水印嵌入方法研究,目前尚未见报道,尚属空白。而实际应用中,医学图像数字水印常常同时受到这两种攻击。

发明内容

[0006] 本发明的目的是提供一种基于DFT可抗几何攻击的医学图像多重水印方法,将医学图像的视觉特征向量、加密技术和第三方的概念有机结合起来,不需要进行感兴趣区域的选取,具有很理想的鲁棒性和不可见性,有效地解决了多重水印嵌入的问题,同时解决了医学图像应用中出现的抗击几何攻击和抗击常规攻击问题,以保护医学图像的版权和病患信息的隐秘性。
[0007] 为了实现上述目的,本发明是这样进行的:基于全图DFT变换,在DFT变换系数中,提取一个抗几何攻击的医学图像视觉特征向量,并将水印技术与密码学有机结合起来,实现了多重数字水印的抗几何和常规攻击。本发明所采用的方法包括水印嵌入和水印提取两大部分,第一部分为多重水印嵌入方法,包括:(1)通过进行全图DFT变换,得到图像的k一个视觉特征向量V(j);(2)根据要嵌入的多重水印W(j),k=1,2,...,n;n表示嵌入的水印个数和在医学图像中提取的特征向量V(j),通过Hash函数运算,生成二值逻辑序列k k
Key(j),然后将二值逻辑序列Key(j)存在第三方。第二部分为多重水印提取方法,包括:
k(
(3)求出待测图像的视觉特征向量V’(j);(4)利用已存在第三方的二值逻辑序列Key j)k
和待测医学图像的特征向量V’(j),提取出多重水印W’(j)。
[0008] 现对本发明的方法进行详细说明如下:
[0009] 首先用Wk(j)表示要嵌入的多重水印,Wk(j)={wk(j)|w(j)=0,1;1≤j≤L,1≤k≤n},L代表要嵌入的一个水印长度,n为嵌入水印的个数。原始图像记为F={f(i,j)|f(i,j)∈R;1≤i≤N1,1≤j≤N2)},其中,f(i,j)表示原始医学图像的像素灰度值,设N1=N2=N。
[0010] 第一部分:水印嵌入方法
[0011] 1)通过进行全图DFT变换,得到图像的视觉特征向量V(j)。
[0012] 先对原图F(i,j)进行全图DFT变换,得到DFT系数矩阵FF(i,j),再对DFT系数矩阵FF(i,j),在低中频系数中,取前L个值,并通过DFT系数符号运算得到该图像的视觉特征向量V(j),方便起见,这里一个复数看成实部、虚部(虚部只看系数)两个系数,当系数值为“正”时我们用“1”表示(含系数值为“0”的情况),系数为负时用“0”表示,主要过程描述如下:
[0013] FF(i,j)=DFT2(F(i,j))
[0014] V(j)=-Sign(FF(i,j))
[0015] 2)根据水印Wk(j)和图像的视觉特征向量V(j)生成一个二值逻辑序列Keyk(j)。
[0016]
[0017] Keyk(j)是由图像的视觉特征向量V(j)和水印Wk(j),通过密码学常用的Hash函k k数生成。保存Key(j),在以后提取水印时需用。通过将Key(j)作为密钥向第三方申请,以获得原图像的查看和所有权,从而达到保护医学图像的目的。
[0018] 第二部分:水印提取方法
[0019] 3)求出待测图像的视觉特征向量V’(j)。
[0020] 设待测图像为F’(i,j),经过全图DFT变换后得到DFT系数矩阵为FF’(i,j),按上述Step1方法,求得待测图像的视觉特征向量V’(j);
[0021] FF’(i,j)=DFT2(F’(i,j))
[0022] V’(j)=-Sign(FF’(i,j))
[0023] 4)在待测图像中提取出水印Wk’(j)。
[0024] 根据在嵌入水印时生成的Keyk(j)和待测图像的视觉特征向量V’(j),利用Hashk性质可以提取出待测图像的水印W’(j)。
[0025]
[0026] 再根据Wk(j)和Wk’(j)的相关程度来判别待测图像的所有权和隐患信息的安全性问题。
[0027] 本发明与现有的医学水印技术比较有以下优点:
[0028] 由于本发明是基于DFT变换的多数字水印技术,有较强的抗几何攻击能力和抗常规攻击能力;不需要人为的进行感兴趣区域的选取,提高另外水印的嵌入速度,并且该发明解决了多重水印嵌入的问题;嵌入的多重水印是一种零水印,不影响原始医学图像质量,在医疗方面具有很高的实用价值。
[0029] 以下从理论基础和试验数据说明:
[0030] 1)离散傅里叶变换
[0031] 二维离散傅里叶正变换(DFT)公式如下:
[0032]
[0033] u=0,1,Λ,M-1;v=0,1,Λ,N-1;
[0034] 二维离散傅里叶反变换(IDFT)公式如下:
[0035]
[0036] x=0,1,Λ,M-1;y=0,1,Λ,N-1
[0037] 其中x,y为空间域采样值;u,v为频率域采样值,通常数字图像用像素方阵表示,即M=N
[0038] 从上面的公式可知,DFT的系数符号是和分量的相位有关的。
[0039] 2)医学图像视觉主要特征向量的选取方法
[0040] 目前大部分医学图像水印算法抗几何攻击能力差的主要原因是:人们将数字水印嵌入在像素或变换系数中,医学图像的轻微几何变换,常常导致像素值或变换系数值的有较大变化。这样便会使嵌入的水印很轻易的就受到攻击。如果能够找到反映图像几何特点的视觉特征向量,那么当图像发生小的几何变换时,该图像的视觉特征值不会发生明显的突变。Hayes研究表明对图像特征而言,相位比幅度更重要。我们对大量的全图DFT数据(低中频)经过观察发现,当对一个医学图像进行常见的几何变换时,低中频系数大小可能发生一些变化,但其系数符号基本保持不变,我们选取一些实验数据见表1所示。表1中用作测试一个切片图像(128x128)。表中第1列显示的是医学图像受到攻击的类型,受到常规攻击后的医学图像见图1(b)-(d),受到几何攻击后的医学图像见图2(a)-(d)。第3列到第6列,这是在DFT系数矩阵中取的FF(1,1)-FF(1,5),共5x2=10个低中频系数(这里把一个复数,看成实部和虚部两个系数)。其中系数F(1,1)表示医学图像的直流分量值。对于常规攻击,这些低中频系数值FF(1,1)-FF(1,5)基本保持不变,和原始医学图像值近似相等;对于几何攻击,部分系数有较大变化,但是我们可以发现,医学图像在受到几何攻击时,部分DFT低中频系数的大小发生了变化但其符号基本没有发生变化。我们将傅里叶系数(这里复数看成实部和虚部两个系数值),正值和零用“1”表示,负值用“0”表示,那么对于原始医学图像来说,DFT系数矩阵中的FF(1,1)-FF(1,5)系数,对应的系数符号序列为:“1100001111”,见表1的第7列,观察该列可以发现,无论常规攻击还是几何攻击该符号序列和原始医学图像能保持相似,与原始医学图像归一化相关系数都较大(见第8列),(方便起见这里取了5个DFT系数符号)。
[0041] 为了进一步证明全图的DFT变换系数符号序列是属于该图的一个视觉重要特征,又把不同的测试图像,见图3(a)-(g),进行全图DFT变换,得到对应的DFT系数FF(1,1)-FF(4,4),并且求出与原图的符号序列的相关系数,计算结果如表2所示。从表2可以看出,不同医学图像之间,符号序列相差较大,相关度较小,小于0.5。
[0042] 这更加说明DFT系数的符号序列可以反映该医学图像的主要视觉特征。当水印图像受到一定程度的常规攻击和几何攻击后,该向量基本不变,这也符合DFT“有很强的提取图像特征”能力。
[0043] 表1图像全图DFT变换低中频部分系数及受不同攻击后的变化值
[0044]
[0045] *DFT变换系数单位1.00e+003
[0046] 表2不同医学图像特征向量的相关系数(向量长度32bit)
[0047]Pa Pb Pc Pd Pe Pf Pg
Pa 1.00 0.38 0.25 -0.18 0.12 -0.26 0.00
Pb 0.38 1.00 0.38 -0.11 -0.12 0.14 -0.13
Pc 0.25 0.38 1.00 -0.01 -0.25 0.24 0.13
Pd -0.18 -0.11 -0.01 1.00 0.25 0.09 0.27
Pe 0.12 -0.12 -0.25 0.25 1.00 -0.01 0.38
Pf -0.26 0.14 0.24 0.09 -0.01 1.00 0.26
Pg 0.00 -0.13 0.13 0.27 0.38 0.26 1.00
[0048] 3)水印嵌入的位置和一次性嵌入的长度
[0049] 根据人类视觉特性(HVS),低中频信号对人的视觉影响较大,代表着医学图像的主要特征。因此所选取的医学图像的视觉特征向量是低中频系数的符号,低中频系数的个数选择与进行全图DFT变换的原始医学图像的大小、以及一次性嵌入的信息量和要求的鲁棒性有关,L值越小,一次性嵌入的信息量越少,但鲁棒性越高。在后面的试验中,选取L的长度为32。
[0050] 综上所述,通过对医学图像的全局DFT系数的分析,利用DFT低中频系数的符号序列得到一种取得医学图像的一个抗几何攻击的特征向量的方法,利用该特征向量和Hash函数、“第三方”概念实现了在医学图像中嵌入多水印的方法。经过实验证明,该方法实现了多水印的嵌入,并且水印的嵌入不影响医学图像的内容,并且有较好的鲁棒性。

附图说明

[0051] 图1(a)是原始医学图像。
[0052] 图1(b)是经过高斯干扰的图像。
[0053] 图1(c)是经过JPEG攻击的图像。
[0054] 图1(d)是经过中值滤波的图像。
[0055] 图2(a)是经过旋转变换的图像。
[0056] 图2(b)是经过缩放2.0的图像。
[0057] 图2(c)是经过缩放0.5的图像。
[0058] 图2(d)是经过垂直移动的图像。
[0059] 图3(a)是标准测试图MRI_1。
[0060] 图3(b)是标准测试图MRI_2。
[0061] 图3(c)是标准测试图MRI_3。
[0062] 图3(d)是标准测试图Engine。
[0063] 图3(e)是标准测试图Head。
[0064] 图3(f)是标准测试图Teddy bear。
[0065] 图3(g)是标准测试图Mri_1back。
[0066] 图4(a)没有干扰时的水印图像。
[0067] 图4(b)没有干扰时的水印检测。
[0068] 图5(a)有高斯干扰时的水印图像(高斯干扰强度为3%)。
[0069] 图5(b)有高斯干扰时的水印检测。
[0070] 图6(a)JPEG压缩后的水印图像(压缩质量为4%)。
[0071] 图6(b)JPEG压缩后的水印检测。
[0072] 图7(a)中值滤波后的水印图像(经过[3,3]的20次滤波)。
[0073] 图7(b)中值滤波后的水印检测。
[0074] 图8(a)旋转20度后的水印图像。
[0075] 图8(b)旋转20度后的水印检测。
[0076] 图9(a)缩放因子为4.0的水印图像。
[0077] 图9(b)缩放因子为4.0的水印检测。
[0078] 图10(a)缩放因子为0.5的水印图像。
[0079] 图10(b)缩放因子为0.5的水印检测。
[0080] 图11(a)移动3%后的图像。(水平左移)
[0081] 图11(b)移动3%后的水印检测。
[0082] 图12(a)剪切6%的水印图像。
[0083] 图12(b)剪切6%的水印检测。

具体实施方式

[0084] 下面结合附图对本发明作进一步说明使用1000组独立的二值伪随机序列(取值为+1或0),每组序列长度为32bit,在这1000组数据中,任抽取三组(这里选择第300组、第500组、第700组),作为嵌入的三个水印序列,即我们嵌入了总长度为32x3=96bit的水印序列。实验所用的原始医学图像,是一大脑的个切片的图像(128x128)见图4(a)。设原图表示为F(i,j),其中1≤i≤128,1≤j≤128,对应的全图DFT系数矩阵为FF(i,j),选择低中频系数Y(j),1≤j≤L,第一个值Y(1)代表图像的直流分量,然后由低到高的频率顺序排列。考虑到鲁棒性和一次性嵌入水印的容量,我们选择中低频的4x4=16个复数系数做特征向量(这里把一个复数,看成实部和虚部两个系数),则共有16x2=32个低中频k系数,即L=32。嵌入的多水印W由k个子水印W(j)组成,本例中子水印的个数k取3,这k
里子水印记为W(j),1≤j≤32,1≤k≤3;选取的DFT系数矩阵为FF(i,j),1≤i≤4,k k k
1≤j≤4。通过水印算法检测出W′(j)后,再通过算W(j)和W′(j)的归一化相关系k
数NC(Normalized Cross Correlation),方便起见,用NC1、NC2和NC3表示与提取出的三个水印对应的三个相关系数,用于判断是否有水印嵌入。
[0085] 图4(a)是不加干扰时的水印图像;
[0086] 图4(b)不加干扰时,水印检测器的输出,可以看到NC1=1.00,NC2=1.00,NC3=1.00,明显检测到水印的存在。
[0087] 下面我们通过具体试验来判断该数字水印方法的抗常规攻击能力和抗几何攻击能力鲁棒性。
[0088] 先测试该水印算法抗常规攻击的能力。
[0089] (1)加入高斯噪声
[0090] 使用imnoise()函数在水印图像中加入高斯噪音。
[0091] 图5(a)为当高斯噪声强度为3%时的水印图像,在视觉上已很模糊;
[0092] 图5(b)水印检测器的输出,能很明显的检测到水印的存在,NC1=1.00,NC2=1.00,NC3=1.00。
[0093] 表3是水印抗高斯干扰时的检测数据。从实验数据可以看到,当高斯噪声强度高达为25%时,水印图像PSNR降至0.13dB,这时检测水印,相关系数NC1=0.90,NC2=0.93,NC3=0.93,仍能检测出水印的存在.这说明采用该发明有好的抗高斯噪声能力。
[0094] 表3水印抗高斯噪声干扰数据
[0095]噪声强度(%) 1 3 5 10 15 20 25
PSNR(dB) 12.43 7.94 6.01 3.25 1.85 0.84 0.13
NC1 1.00 1.00 1.00 1.00 0.96 0.90 0.90
NC2 1.00 1.00 1.00 1.00 0.95 0.93 0.93
NC3 1.00 1.00 1.00 1.00 0.93 0.93 0.93
[0096] (2)JPEG压缩处理
[0097] 采用图像压缩质量百分数作为参数对水印图像进行JPEG压缩;
[0098] 图6(a)是压缩质量为4%的图像,该图已经出现方块效应;
[0099] 图6(b)是水印检测器的响应,NC1=1.00,NC2=1.00,NC3=1.00,检测效果明显。
[0100] 表4为水印图像抗JPEG的试验数据。当压缩质量为很差,压缩质量为2%时,仍然可以测得水印的存在,NC1=0.75,NC2=0.75,NC3=0.75。
[0101] 表4水印抗JPEG压缩的实验数据
[0102]压缩质量(%) 2 4 8 10 20 40 60 80
PSNR(dB) 16.32 17.61 19.99 20.98 23.04 25.06 26.52 29.27
NC1 0.75 1.00 1.00 0.88 1.00 1.00 1.00 1.00
NC2 0.75 1.00 1.00 0.90 1.00 1.00 1.00 1.00
NC3 0.75 1.00 1.00 0.89 1.00 1.00 1.00 1.00
[0103] (3)中值滤波处理
[0104] 图7(a)是中值滤波参数为[3x3],滤波重复次数为20的医学图像,图像已出现模糊;
[0105] 图7(b)是水印检测器的响应,NC1=0.83,NC2=0.82,NC3=0.85,检测效果明显。
[0106] 表5为水印图像抗中值滤波能力,从表中看出,当中值滤波参数为[7x7],滤波重复次数为20时,仍然可以测得水印的存在,NC1=0.75,NC2=0.75,NC3=0.75。
[0107] 表5水印抗中值滤波实验数据
[0108]
[0109] 水印抗几何攻击能力:
[0110] (1)旋转变换
[0111] 图8(a)是水印图像旋转20°,这时水印图像的PSNR=12.38dB,信噪比很低;
[0112] 图8(b)为检测的水印图像,可以明显检测到水印的存在NC1=0.89,NC2=0.87,NC3=0.88。
[0113] 表6为水印抗旋转攻击试验数据。表中可以看到当水印图像旋转30°时,NC1=0.69,NC2=0.68,NC3=0.71,仍然可以检测到水印存在;Pitas等人提出的抗几何攻击算法,把水印嵌入DFT幅度谱的圆环中,只能抵抗不大于3度的旋转。
[0114] 表6水印抗旋转攻击实验数据
[0115]
[0116] (2)缩放变换
[0117] 图9(a)是当缩放因子4.0时的水印图像,这时中心图像比原图的要大;
[0118] 图9(b)为水印检测结果,可以检测到水印的存在,NC1=1.00,NC2=1.00,NC3=1.00。
[0119] 图10(a)为缩放因子为0.5的水印图像,这时中心图像比原图的要小好多;
[0120] 图10(b)是水印检测结果,可以明显检测到水印的存在NC1=1.00,NC2=1.00,NC3=1.00。
[0121] 表7为水印缩放攻击试验数据,从表8可以看到当水印图像缩放因子小至0.2时,相关系数NC1=0.75,NC2=0.75,NC3=0.75,仍可测得水印。Pereira等采用的在DFT中置入模板的方法,只能抵御缩放因子不小于0.65的缩放,说明该发明有较强的抗缩放能力。
[0122] 表7水印缩放攻击实验数据
[0123]缩放因子 0.2 0.5 0.8 1.00 1.2 2.0 4.0
NC1 0.75 1.00 0.94 1.00 1.00 1.00 1.00
NC2 0.75 1.00 0.95 1.00 1.00 1.00 1.00
NC3 0.75 1.00 0.94 1.00 1.00 1.00 1.00
[0124] (3)平移变换
[0125] 图11(a)为图像水平左移3%的情况,这时PSNR=12.28dB,信噪比很低;
[0126] 图11(b)为水印检测器输出,可以明显检测到水印的存在NC1=0.94,NC2=0.93,NC3=0.94。
[0127] 表8是水印抗平移攻击试验数据。从表中得知当水平或垂直移动10%,仍然可以检测到水印的存在,故该数字水印有较强的抗平移能力。
[0128] 表8水印抗平移攻击实验数据
[0129]
[0130] (4)剪切试验
[0131] 图12(a)为对水印图像按Y轴方向进行剪切6%的情况,这时顶部相对于原始医学图像,已被剪切掉一部分了;
[0132] 图12(b)为其水印检测情况,可以明显检测到水印的存在,NC1=0.87,NC2=0.90,NC3=0.87。
[0133] 表9为水印抗切割试验数据,从表中试验数据可以得知,该算法有一定的抗剪切能力。
[0134] 表9水印抗剪切攻击实验数据(按Y轴方向剪切)
[0135]
[0136] 通过以上的实验说明,该水印的嵌入方法,有较强的抗常规攻击能力和几何攻击能力,并且水印的嵌入不影响医学图像的值,是一种零水印。
高效检索全球专利

IPRDB是专利检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,专利查询、专利分析

电话:13651749426

侵权分析

IPRDB的侵权分析产品是IPRDB结合多位一线专利维权律师和专利侵权分析师的智慧,开发出来的一款特色产品,也是市面上唯一一款帮助企业研发人员、科研工作者、专利律师、专利分析师快速定位侵权分析的产品,极大的减少了用户重复工作量,提升工作效率,降低无效或侵权分析的准入门槛。

立即试用