• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 29
    • 下载费用:30 金币  

    重庆时时彩推算方法: 基于三维DWTDCT感知哈希和混沌的体数据水印.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201410145427.0

    申请日:

    2014.04.10

    公开号:

    CN103942747A

    公开日:

    2014.07.23

    当前法律状态:

    撤回

    有效性:

    无权

    法律详情: 发明专利申请公布后的视为撤回IPC(主分类):G06T 1/00申请公布日:20140723|||实质审查的生效IPC(主分类):G06T 1/00申请日:20140410|||公开
    IPC分类号: G06T1/00 主分类号: G06T1/00
    申请人: 海南大学
    发明人: 李京兵; 黄梦醒; 周又玲; 陈晓枫; 王英杰
    地址: 570228 海南省??谑忻览记嗣翊蟮?8号
    优先权:
    专利代理机构: 代理人:
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410145427.0

    授权公告号:

    ||||||

    法律状态公告日:

    2017.06.30|||2014.08.27|||2014.07.23

    法律状态类型:

    发明专利申请公布后的视为撤回|||实质审查的生效|||公开

    摘要

    本发明公开了一种基于三维DWT-DCT感知哈希和混沌的体数据水印技术,属于多媒体信号处理领域。本发明的步骤是:首先利用Logistic?Map生成混沌序列,利用混沌序列对水印进行置乱;然后进行水印的嵌入,对医用体数据进行全局3DDWT-DCT变换,选取前4×4×2个系数,再进行3D-IDCT变换,然后在反变换系数中通过感知哈希提取一个具有鲁棒特性的感知哈希值,并将水印序列与该感知哈希值相关联得到一个二值密钥序列,并将该二值密钥序列存于第三方;再通过对待测体数据进行三维DWT-DCT体数据的感知哈希值的提取,并与存于第三方的二值序列相关联来进行水印的提取,最后利用混沌序列对水印进行还原。

    权利要求书

    权利要求书
    1.  基于三维DWT-DCT感知哈希和混沌的体数据水印,其特征在于:对体数据进行三维DWT-DCT变换,选取前4×4×2个系数,再进行3D-IDCT变换,然后在反变换系数中提取一个鲁棒的感知哈希值,并将水印序列与感知哈希值相关联,实现了医学体数据数字水印的抗几何攻击和常规攻击,该体数据数字水印实现方法共分四大部分,共计八个步骤:
    第一部分是水印的混沌置乱:利用LogisticMap性质产生混沌序列对水印进行置乱,得到混沌置乱的水印BW(i,j);
    1)由逻辑初始值x0通过LogisticMap生成混沌序列X(j);
    2)将混沌序列X(j)中的值按照从小到大的顺序排列,再根据X(j)中各个值排序前后的位置变化对水印像素的空间位置进行置乱,得到混沌置乱的水印BW(i,j);
    第二部分:水印的嵌入
    3)通过三维DWT-DCT,求得原始体数据的一个感知哈希值H(j);即:对原始体数据进行三维小波变换,然后对近似系数进行三维DCT变换,选取低频的前4×4×2变换系数进行三维反DCT变换,求取反变换后系数的平均值,然后将每个反变换后的系数与平均值进行比较,进行二值量化处理,大于或等于平均值,记为1;小于平均值,记为0,得到体数据的感知哈希值H(j);
    4)利用密码学Hash函数性质和水印序列W(i,j),生成含水印信息的二值密钥序列Key(i,j),
    保存Key(i,j),下面提取水印时要用到,通过把Key(i,j)作为密钥向第三方申请,以获得对原始医学体数据的所有权;
    第三部分:水印的提取
    5)对于待测体数据,按照步骤3)的方法,求出待测体数据的感知哈希值H’(j);
    6)利用存在于第三方的二值逻辑密钥序列Key(i,j)和待测体数据的感知哈希值H’(j),提取出水印
    第四部分水印的还原;
    7)由逻辑初始值x0通过LogisticMap生成混沌序列X(j);
    8)对混沌序列X(j)中的值由小到大进行排序,根据X(j)中各个值排序前后的位置变化,对提取的水印像素的空间位置进行还原,得到还原的水印W’(i,j);
    将W(i,j)和W’(i,j)进行归一化相关系数计算,来确定是否有水印嵌入。

    说明书

    说明书基于三维DWT-DCT感知哈希和混沌的体数据水印
    技术领域
    本发明涉及一种基于三维DWT-DCT感知哈希和混沌的体数据数字水印技术,是一种多媒体数据?;し椒?,属于多媒体信号处理领域。
    技术背景
    利用互联网可以实现远程医疗诊断,医学体数据在网络上进行远程传输时,记录在医学图片上的病人的个人信息,很容易被泄露。如何?;せ颊叩母鋈艘?、使得病人的CT、MRI等医用图像上的个人信息、病人的电子病历等数据不被泄露,这一问题随着互联网的普及变得日益严重。现今的加密方法和访问控制已很难满足医学体数据信息安全的要求,若把个人信息作为数字水印嵌入在医学图片中,即医学图像数字水印技术(Medical Image Watermarking,简称MIW)可以有效地解决这一问题。
    数字水印技术最初是用于互联网上的数字媒体的版权?;?,现在利用数字水印的不可见性、鲁棒性等特点,可以把病人的个人信息隐藏在其医学图像中,以保证它在互联网上的安全传输。医学图像数字水印的出现,使得远程医疗诊断、远程手术所需的相关病人资料在互 联网上传输时,可以有效的?;げ∪说囊?,避免病人的资料被篡改。
    目前对于抗几何攻击的体数据多数字水印算法的研究较少。而体数据在医学图像中大量存在,如:CT、MRI图像都是由切片组成的体数据,因此研究如何在体数据中嵌入多数字水印有较大意义,并且对于医学体数据,一般是不允许修改其内容的。这又为在体数据中嵌入多水印提高了难度。
    另外,下一代图像压缩标准JPEG2000是基于小波变换的。因此,对利用三维DWT-DCT感知哈希,在体数据中嵌入水印的工作有较大意义。
    发明内容
    本发明的目的是提出一种基于三维DWT-DCT感知哈希和混沌置乱,实现既能抗击几何攻击又能抗击常规攻击的体数据水印嵌入与提取方法,它有较高的鲁棒性,并且水印的嵌入不影响原始体数据的体素数据值,是一种零水印方案,并且水印在嵌入前进行了混沌置乱从而提高了水印嵌入的安全性。
    为了实现上述目的,本发明是这样进行的:应用Logistic Map产生混沌序列对水印进行混沌置乱和还原,提高水印的安全性;首先对医用体数据进行全局3D-DWT变换,得到“近似系数”和“细节系数”,这类似二维图像的小波变换,“近似系数”代表体数据的低频特性,反映的是体数据的主要外部轮廓;“细节系数”代表体数据的高频特性,反映的是体数据的高频细节信息。小波变换本身的抗几何攻击能力较差,因此,我们先对体数据进行三维小波变换(DWT),然后 再对反映低频特性的“近似系数”进行全局余弦变换(DCT),选取前4×4×2个系数,再进行反变换,3D-IDCT,求取反变换后系数的平均值,然后将每个反变换后的系数与平均值进行比较,进行二值量化处理,大于或等于平均值,记为1;小于平均值,记为0,得到体数据的感知哈希值。感知哈希具有鲁棒性和不可感知性,并将水印与提取的感知哈希值相关联,实现体数据的鲁棒数字水印技术。本发明所采用的方法包括水印的混沌置乱、水印的嵌入、提取和还原四大部分,第一部分为水印的混沌置乱,包括:(1)通过Logistic Map产生混沌序列X(j);(2)根据X(j)对水印进行置乱,得到混杂置乱的水印BW(i,j);第二部分为水印的嵌入,包括:(3)通过三维DWT-DCT感知哈希算法,得到原始体数据的鲁棒感知哈希H(j),(4)利用密码学HASH函数性质,嵌入混沌置乱后的水印BW(i,j),生成含水印信息的二值密钥序列Key(i,j),然后将二值逻辑序列Key(i,j)存在第三方;第三部分为水印的提取,包括:(5)求出待测体数据的感知哈希值H’(j),(6)利用存在于第三方的二值逻辑密钥序列Key(i,j)和待测体数据的感知哈希值H’(j),提取出水印BW’(i,j);第四部分为水印的还原,包括:(7)应用LogisticMap,得到相同的混沌序列X(j),(8)通过X(j)对水印进行还原。
    现对本发明的方法进行详细说明如下:
    首先选择一幅有意义的二值图像作为要嵌入医学体数据的水印,记为W={w(i,j)|w(i,j)=0,1;1≤i≤M1,1≤j≤N1};同时,选取Matlab中自带的一个MRI体数据作为原始医学体数据,表示为: F={f(i,j,k)|f(i,j,k)∈R;1≤i≤M,1≤j≤N,1≤k≤P}。其中,f(i,j,k)表示原始医学体数据的体素(Voxel)数据值,这类似二维图像中的像素灰度值,方便起见,设M1=M2,M=N。
    第一部分:对水印的混沌置乱
    1)通过Logistic Map生成混沌序列;
    由初始值x0通过Logistic Map混沌系统生成混沌序列X(j)。
    2)得到混沌置乱的水??;
    首先,将原始水印转化为二值水印W(i,j),然后,将混沌序列X(j)中的值按照从小到大的顺序进行排序,最后,根据X(j)中各个值排序前后的位置变化对水印像素的位置空间进行置乱,得到混沌置乱的水印BW(i,j)。
    第二部分:水印的嵌入
    3)通过三维DWT-DCT感知哈希算法,得到原始体数据的一个鲁棒感知哈希值H(j);
    先对原始体数据F(i,j,k)进行三维小波变换,得到逼近子图系数FAL,再对逼近子图FAL进行全局三维DCT变换,得到DWT-DCT系数矩阵FD(i,j,k),在系数矩阵FD(i,j,k)中选取前4×4×2个系数FD4(i,j,k),再对选取出的系数矩阵FD4(i,j,k)进行三维反DCT变换,得到反变换后的系数FID(i,j,k),求取反变换后系数的平均值,然后将每个反变换后的系数与平均值进行比较,进行二值量化处理,大于或等于平均值,记为1;小于平均值,记为0,得到体数据的感知哈希值H(j)。主要过程描述如下:
    FAL(i,j,k)=DWT3(F(i,j,k))
    FD4(i,j,k)=DCT3(FAL(i,j,k))
    FID(i,j,k)=IDCT3(FD4(i,j,k))
    H(j)=BINARY(FID(i,j,k))
    4)利用HASH函数,嵌入水??;生成含水印信息的二值密钥序Key(i,j);
    Key(i,j)=H(j)⊕BW(i,j)]]>
    Key(i,j)是由体数据的感知哈希值H(j)和水印图像BW(i,j),通过密码学常用的Hash函数生成。保存Key(i,j),在下面提取水印时要用到。通过将Key(i,j)作为密钥向第三方申请,以获得医学体数据的所有权和使用权,达到版权?;さ哪康?。并且水印的嵌入不影响原始医学体数据的质量,是一种零水印方案。
    第三部分:水印的提取
    5)求出待测数据的感知哈希值H’(j);
    设待测体数据为F’(i,j,k),经过三维小波变换得到逼近子图系数FA’L,再对逼近子图FA’L进行全局三维DCT变换,得到系数矩阵FD’(i,j,k),选取低频部分的前4x4x2系数矩阵,然后进行反DCT变换,再按上述步骤3)类似的方法,求得待测体数据的感知哈希值H’(j);
    FA’L(i,j,k)=DWT3(F’(i,j,k))
    FD’4(i,j,k)=DCT3(FA’L(i,j,k))
    FID’(i,j,k)=IDCT3(FD’4(i,j,k))
    H’(j)=BINARY(FID’(i,j,k))
    6)在待测体数据中提取出水印BW’(i,j);
    BW,(i,j)=Key(i,j)⊕H,(j)]]>
    根据在嵌入水印时生成的逻辑密钥序列Key(i,j)和待测体数据的感知哈希值H’(j),利用Hash函数性质可以提取出待测体数据中含有的水印BW’(i,j)。
    第四部分:水印的还原
    7)通过Logistic Map生成混沌序列;
    由与上文step1相同的初始值x0通过Logistic Map混沌系统生成相同的混沌序列X(j);
    8)还原提取的水??;
    首先将混沌序列X(j)中的值按照从小到大进行排序,然后根据X(j)中各个值排序前后的位置变化对水印像素的位置空间进行还原得到还原的水印W’(i,j)。
    再根据W(i,j)和W’(i,j)的相关程度来判别是否有水印嵌入,从而确认待测体数据的所有权和病患信息的安全性问题。
    本发明与现有的医学水印技术比较有以下优点:
    首先,由于本发明是基于三维DWT-DCT感知哈希算法的数字水印技术,通过后面的实验数据证实,该水印不仅有较强的抗常规攻击能力,而且有较强的抗几何攻击能力;其次,嵌入的水印是经过Logistic Map混沌置乱的,使得水印信息变得杂乱无章,提高了水印信息的安全性;最后,水印的嵌入不影响原始体数据的体素数据值, 是一种零水印技术,更好得?;ち艘窖迨?。这个特性,尤其是在医疗图像处理等方面具有很高的实用价值,使用范围广,并且可实现大水印的嵌入与提取。
    以下我们从理论基础和试验数据说明:
    1)三维离散小波变换(DWT)
    三维小波变换的一层分解过程如图1所示,图1中的L、H分别表示医学体数据经过低频和高频滤波之后得到的低频成分和高频成分,与二维图像的小波变换类似,医学体数据经过三维小波变换后,被分解成一个代表体数据低频特性的“近似系数”LLL1(低频三维子带),和代表该体数据的高频信息的“细节系数”(高频三维子带),下标“1”表示三维DWT的第一层分解;一个体数据的三维小波变换(两层)的三维示意图见图4,图2为体数据的一个切片,图3为体数据的三维成像,图4为体数据的三维小波变换(两层)。
    2)三维离散余弦变换(3D-DCT)
    三维DCT变换公式如下:
    对应大小为M×N×P体数据,三维离散余弦正变换(DCT)公式如下:
    F(u,v,w)=c(u)c(v)c(w)[Σx=0M-1Σy=0N-1Σp=0P-1f(x,y,z)·cos(2x+1)2Mcos(2y+1)2Ncos(2z+1)2P]u=0,1,...,M-1;v=0,1,...,N-1;w=0,1,...,P-1;]]>
    式中
    c(u)=1/Mu=02/Mu=1,2,...,M-1]]>
    c(v)=1/Nv=02/Nv=1,2,...,N-1]]>
    c(w)=1/Pw=02/Pw=1,2,...,P-1]]>
    这里,f(x,y,z)是体数据V在(x,y,z)处的体素(voxel)数据值,F(u,v,w)是该体素数据对应的3D-DCT变换系数。
    三维离散余弦反变换(IDCT)公式如下:
    f(x,y,z)=Σu=0M-1Σv=0N-1Σw=0P-1[c(u)c(v)c(w)F(u,v,w)cos(2x+1)2Mcos(2y+1)2Ncos(2z+1)2P]x=0,1,...,M-1;y=0,1,...,N-1;z=0,1,...,P-1]]>
    其中,(x,y,z)为空间域采样值;(u,v,w)为频率域采样值。医学体数据的获得可通过CT和MRI(Magnetic Resnane Iamge,磁共振成像),体数据(Volume data)是由许多层的切片(slice)组成,每个切片是一个二维图像,大小为M×N,切片的层数为P。
    3)Logistic Map
    混沌是一种貌似无规则的运动,指在确定性系统中出现的类似随机的过程。因此,有了它的初始值和参数,我们就能够生成这个混沌系统。Logistic Map是最著名的一种混沌系统之一,它是由以下公式给出的非线性映射:
    xk+1=μxk(1-xk)
    其中,0≤μ≤4为增长参数,xk∈(0,1)为系统变量,k是迭代次数?;煦缍ο低车难芯抗ぷ髦赋?,当增长参数3.569945≤μ≤4时,Logistic Map工作于混沌状态??梢钥吹匠跏贾涤幸桓鲂⌒〉牟煌岬贾禄煦缧蛄械南灾钜?。因此,以上序列是一个理想的密钥序列。 本文中设定μ=4,混沌序列由不同的初始值x0产生。
    4)体数据的一个感知哈希函数的选取方法
    目前大部分水印算法抗几何攻击能力差的主要原因是:人们将数字水印嵌入在体素或变换系数中,体数据的轻微几何变换,常?;岬贾绿逅厥葜祷虮浠幌凳档耐蝗唤洗蟊浠?。这样嵌入在体数据中的多水印便被轻易攻击。如果能够找到一个反映体数据几何特点的感知哈希值,当体数据发生小的几何变换时,该感知哈希值不会发生明显的突变,然后我们把要嵌入的数字水印和该体数据的感知哈希值相关联,就可以较好地解决水印的鲁棒性问题。我们选取一些常规攻击和几何攻击的实验数据见表1所示。表1中用作测试的原图是图5,是matlab中自带的一个MRI体数据的一个切片(取第十个),表1中“第1列”显示的是体数据受到攻击的类型,受到常规攻击后的该切片图像见图6至图8,常规攻击对应的三维成像见图9至图12;受到几何攻击后的切片图像见图13至图16,其对应的三维成像见图17至图20。表1的“第2列”表示的是体数据受到攻击后的峰值信噪比(PSNR);表1的“第3列”到“第10列”表示的是从反DCT变换后变换系数的实部中任意选取“F(1,1,4)、F(1,2,1)”等八个系数值。表1的“第11列”是DWT-DCT感知哈希算法二值量化处理后求出来的平均值。对于常规攻击或是几何攻击,这些系数值F(1,1,4)、F(1,2,1)等可能发生一些变换,但是它与平均值的大小关系仍然不变,我们将大于或等于平均值,记为1;小于平均值,记为0,那么对于原始体数据来说,系数值F(1,1,4)、F(1,2,1)等对应的哈希值序列为:“00010101”,具体 见表1的第12列,观察该列可以发现,无论常规攻击还是几何攻击该符号序列和原始体数据的保持相似,与原始体数据归一化相关系数都较大,为1.0,见表1“第13列”(方便起见这里取了8个三维反DCT变换系数符号)。
    表1 基于DWT-DCT体数据感知哈希值受不同攻击后的变化值

                                 *DWT-DCT与感知哈希处理后系数单位1.0e+002
    为了进一步证明按上述方法提取的感知哈希值,具有不可碰撞性,即对于不同的体数据,感知哈希值不同;我们又把不同的测试对象(见图21至图27),通过三维DWT-DCT感知哈希算法对它们进行处理。从统计学角度,这里取了前8×8×4个DWT-DCT系数。并且求出每个体数据的感知哈希值序列相互之间的相关系数,计算结果如表2所示,从图21至图27求取的感知哈希值分别为Ha至Hg。
    从表2可以看出,首先感知哈希值自身之间的相关系数最大,为1.00;其次,图26和图27之间的相关系数也较大为0.69,而这两个图是形状相似的两个肝的体数据;图21和图22,相关系数为0.62, 也较大,在表中为第三大相关系数,而这两个图都是人体的头部,也比较相似。其它感知哈希值之间的相关系数值较小,这与我们人眼观察到的相符合,这说明按该发明的方法提取的感知哈希值,反映了体数据的主要外形特征,有较好的鲁棒性和不可感知性。
    表2 不同体数据感知哈希值之间的相关系数(向量长度256bit)
     HaHbHcHdHeHfHgHa1.000.620.46-0.10-0.510.110.14Hb0.621.000.260.08-0.510.180.18Hc0.460.261.00-0.05-0.37-0.12-0.14Hd-0.100.08-0.051.00-0.100.090.12He-0.51-0.51-0.37-0.101.00-0.24-0.21Hf0.110.18-0.120.09-0.241.000.69Hg0.140.18-0.140.12-0.210.691.00
    综上所述,我们利用三维DWT-DCT感知哈希值,得到体数据多水印的嵌入与提取方法。
    附图说明
    图1是三维小波变换示意图(一层)。
    图2是原始体数据的一个切片。
    图3是原始体数据对应的三维成像。
    图4是对原始体数据进行三维小波变换(两层)的结果显示。
    图5是原始体数据的一个切片(默认是体数据的第10个切片)。
    图6是经过10%的高斯噪声干扰后的切片图像。
    图7是经过JPEG压缩(压缩质量为2%)后的切片图像。
    图8是经过中值滤波后的切片图像(滤波参数为[5x5]))。
    图9是原始体数据对应的三维成像。
    图10是体数据受强度为10%的高斯干扰后对应的三维成像。
    图11是JPEG压缩(压缩质量为2%)后对应的三维成像。
    图12是经过中值滤波后对应的三维成像(滤波参数为[5x5])。
    图13是经过顺时旋转20度的切片图像。
    图14是经过缩放0.5倍的切片图像。
    图15是垂直方向下移10%的切片图像。
    图16是Z轴方向剪切10%后的第一个切片图像。
    图17是顺时旋转20度的三维成像。
    图18是缩放系数为0.5的三维成像。
    图19是垂直方向下移10%的三维成像。
    图20是Z轴方向剪切10%的三维成像。
    图21是体数据MRI_1的三维成像。
    图22是体数据MRI_2的三维成像。
    图23是体数据MRI_3的三维成像。
    图24是体数据Teddybear的三维成像。
    图25是体数据Tooth的三维成像。
    图26是体数据Liver_1的三维成像。
    图27是体数据Liver_2的三维成像。
    图28是原始水印图像。
    图29是混沌置乱后的水印图像。
    图30是不加干扰时的水印切片。
    图31是不加干扰时的三维重建图。
    图32是不加干扰时提取的水印。
    图33是高斯噪声干扰后的切片图像(高斯噪声强度10%)。
    图34是高斯噪声干扰后的三维重建图(高斯噪声强度10%)。
    图35是高斯噪声干扰后提取的水印图像(高斯噪声强度10%)。
    图36是JPEG压缩后的切片图像(压缩质量参数为5%)。
    图37是JPEG压缩后的体数据三维成像(压缩质量参数为5%)。图38是JPEG压缩后提取的水印图像(压缩质量参数为5%)。
    图39是中值滤波后的切片图像(滤波参数为[5x5],滤波次数为10次)。
    图40是中值滤波后的体数据的三维成像(滤波参数为[5x5],滤波次数为10次)。
    图41是中值滤波后提取的水印图像(滤波参数为[5x5],滤波次数为10次)。
    图42是顺时旋转20度后的切片图像。
    图43是顺时旋转20度后体数据的三维成像。
    图44是顺时旋转20度后提取的水印图像。
    图45是缩放系数为0.5的切片图像。
    图46是缩放系数为0.5的三维成像。
    图47是缩放系数为0.5时提取的水印图像。
    图48是垂直下移10%的切片图像。
    图49是垂直下移10%的体数据对应的三维成像。
    图50是垂直下移10%后提取的水印图像。
    图51是X轴方向剪切20%后,体数据的第一个切片图像。
    图52是X轴方向剪切20%后,体数据的三维成像。
    图53是X轴方向剪切20%后,提取的水印图像。
    图54是扭曲频率因子为13时的切片图像。
    图55是扭曲频率因子为13时的体数据的三维成像。
    图56是扭曲频率因子为13时提取的水印图像。
    具体实施方式
    下面结合附图对本发明作进一步说明,选择一幅有意义的二值图像作为原始多水印,记为W={w(i,j)|w(i,j)=0,1;1≤i≤M1,1≤j≤N1},见图28,这里水印的大小为32×32。通过Logistic Map混沌置乱后的水印,见图29,可以明显看到水印发生了很大的变化,安全性提高。仿真平台为Matlab2010a,原始医学体数据的一个切片见图5,是取自matlab中自带的磁共振三维图像体数据(MRI.mat),体数据的大小为128x128x27,见图9。原始体数据表示为F(i,j,k),其中1≤i,j≤128;1≤k≤27,对应的三维DWT-DCT系数矩阵为FD(i,j,k),其中1≤i,j≤128;1≤k≤27??悸堑铰嘲粜院鸵淮涡郧度攵嗨〉娜萘课颐侨∏?×4×2个系数。再进行3D-IDCT变换,并求取反变换后系数的平均值,然后通过平均值对反变换后的系数进行二值量化处理,得到体数据的感知哈希值H(j)。通过水印算法检测出W’(i,j)后,我们通过计算归一化相关系数NC(Normalized Cross Correlation)来判断是 否有水印嵌入。
    不加任何干扰时,图30是不加干扰时的切片图像(这里默认选择第十个切片,测试用体数据共由27个切片组成);
    图31是不加干扰时的体数据三维成像;
    图32是不加干扰时提取的水印,可以看到NC=1.00,可以准确得提取水印。
    下面我们通过具体实验来判断该数字水印方法的抗常规攻击能力和抗几何攻击能力。
    先测试该水印算法抗常规攻击的能力。
    (1)加入高斯噪声
    使用imnoise()函数在水印图像中加入高斯噪声。
    表3是水印抗高斯噪声干扰的实验数据。从中可以看到,当高斯噪声强度高达25%时,水印体数据的PSNR降至0.08dB,这时提取的水印,相关系数NC=1.00,仍能准确得提取水印。这说明采用该发明有好的抗高斯噪声能力。
    图33是高斯噪声强度10%时的切片图像,在视觉上已很模糊;
    图34是对应的体数据三维成像,在视觉上已很模糊,体数据的PSNR=3.30dB,较低;
    图35是提取的水印,能准确得提取水印,NC=1.00。
    表3 水印抗高斯噪声干扰数据
    噪声强度(%)13510152025PSNR(dB)12.538.056.023.301.820.810.08NC1.001.001.001.001.001.001.00
    (2)JPEG压缩处理
    采用图像压缩质量百分数作为参数对水印体数据进行JPEG压缩;表4为水印体数据抗JPEG压缩实验数据。当压缩质量仅为2%,这时压缩质量较低,仍然可以提取出水印,NC=1.00。
    图36是压缩质量为5%的切片图像,该图已经出现方块效应;
    图37是对应的体数据三维成像,该图已经出现立体方块效应,此
    时体数据的PSNR=17.61dB;
    图38是提取的水印,NC=1.00,可以准确提取水印。
    表4 水印抗JPEG压缩实验数据
    压缩质量(%)2481020406080PSNR(dB)16.5717.8220.2121.2023.1025.0626.6129.31NC1.001.001.001.001.001.001.001.00
    (3)中值滤波处理
    表5为水印体数据抗中值滤波能力,从表中看出,当中值滤波参数为[7x7],滤波重复次数为20时,仍然可以测得水印的存在,NC=1.00。
    图39是中值滤波参数为[5x5],滤波重复次数为10的切片图像,图像已出现模糊;
    图40是对应的体数据三维成像,这时耳朵等轮廓已不太分明;
    图41是提取的水印,NC=1.00,可以准确提取水印。
    表5 水印抗中值滤波实验数据

    水印抗几何攻击能力
    (1)旋转变换
    表6为水印抗旋转攻击实验数据。从表中可以看到当水印体数据顺时旋转35°时,NC=0.71,仍然可以提取水印。
    图42是顺时旋转20度的水印切片图像;
    图43是相应的体数据三维成像,这时,水印体数据的信噪比较低,PSNR=12.44dB;
    图44是提取的水印,NC=0.95,可以准确地提取水印。
    表6 水印抗旋转攻击实验数据
    顺时旋转5度10度15度20度25度30度35度PSNR(dB)16.5413.9712.9812.4412.0411.6811.33NC1.001.000.950.950.710.710.71
    (2)缩放变换
    表7为水印体数据抗缩放攻击实验数据,从表7可以看到当水印体数据缩放因子小至0.2时,相关系数NC=1.00,可准确提取出水印。
    图45是缩放后的水印切片图像(缩放因子为0.5);
    图46是缩放攻击后,体数据对应的三维成像(缩放因子为0.5);
    图47是缩放攻击后,提取的水印,NC=1.00,可以准确得提取出水印。
    表7 水印抗缩放攻击实验数据
    缩放因子0.20.50.81.01.22.04.0NC0.721.001.001.001.001.001.00
    (3)平移变换
    表8是水印抗平移变换实验数据。从表中得知当水平或垂直移动10%时,NC值都高于0.5,可以准确提取水印,故该水印方法有较强 的抗平移变换能力。
    图48是切片垂直下移10%的图像;
    图49是体数据的每个切片垂直下移10%后,对应的三维成像,这时PSNR=10.85dB,信噪比较低;
    图50是提取的水印,可以准确提取水印,NC=1.00。
    表8 水印抗平移变换实验数据

    (4)剪切攻击
    表9为水印抗剪切攻击实验数据,从表中可以看到,当从X轴方向剪切,剪切量为60%时,仍然可以提取水印,NC=0.86,说明该水印算法有较强的抗剪切攻击能力。
    图51是按X轴方向剪切20%后,第一个切片图像;
    图52是按X轴方向剪切20%后对应的三维成像,可以发现,剪切攻击的效果明显;顶部相对原图的三维成像,切去了一块。
    图53是提取的水印,可以准确得提取水印,NC=1.00。
    表9 水印抗剪切攻击实验数据
    X轴剪切(%)2468102030405060NC1.001.001.001.001.001.001.001.000.930.86
    (5)扭曲攻击
    表10为水印抗扭曲攻击实验数据,当扭曲因子为24时,这时体数据的信噪比较低PSNR=9.68dB,但这时NC=0.87,仍然可以提取 水印。当扭曲因子为3时,这时体数据的信噪比较低PSNR=10.13dB,但这时NC1=0.73,此时NC值相对较低,但是仍然可以提取多水印。并且从表10中发现,当扭曲因子较低时,对体数据的低频特性影响较大,所以NC值较??;而当扭曲因子较大时,对体数据的高频特性影响较大,即对体数据的外部轮廓影响较小,所以NC值较大;表中的数据与我们在前面对体数据的中低频系数的分析一致。
    图54是扭曲攻击后的切片图像(扭曲因子为13);
    图55是扭曲攻击后对应的体数据三维成像,PSNR=9.83dB,信噪比较低;
    图56是提取的水印,NC=0.87,可以较为准确地提取水印。
    表10 水印抗扭曲攻击实验数据
    扭曲频率因子357913172024PSNR(dB)10.1310.169.899.589.839.869.689.68NC0.730.870.870.870.870.870.870.87

    关 键 词:
    基于 三维 DWTDCT 感知 混沌 数据 水印
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:基于三维DWTDCT感知哈希和混沌的体数据水印.pdf
    链接地址://www.4mum.com.cn/p-6157228.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03