一种基于栅格地图特征的小波域零水印方法技术

本发明专利技术公开了一种基于栅格地图特征的小波域零水印方法。现有的水印的嵌入过程复杂化，消耗计算时间太长，不利于实际应用。本发明专利技术包括水印信息嵌入和水印信息检测。水印信息嵌入是指提取原始栅格地图的零水印信息，保存入库。水印信息检测是指将待检测栅格地图中提取的零水印信息与水印信息库中的水印信息操作得到最后水印信息。本发明专利技术具有较好的隐蔽性，可以克服可逆数字水印中存在的安全漏洞。

【技术实现步骤摘要】

本专利技术属于多媒体信息安全
，特别是针对栅格数字地图数据，提出一种 基于栅格数字地图特征的零水印方法。
技术介绍
目前已经提出了许多数字水印方法，这些方法一般可以分为两类频域水印法和 空域水印法。这些嵌入数字水印的方法都对数据的空域信息或其频域信息做了一定的修改 来嵌入水印信息，为了不让人眼发觉人为修改的痕迹，很多方法采用了基于HVS (人类视觉 系统)的视觉掩膜。用加视觉掩膜的方法在一定程度上解决了水印可感知性和鲁棒性之间 的矛盾。但是加视觉掩膜使得水印的嵌入过程复杂化，消耗计算时间太长，不利于实际应 用，而且如果别有用心者了解加视觉掩膜的方法，不能排除有能力的人用方法篡改含有水 印信息的地图，同样使人眼难以发觉的细小的改动，从而导致水印信息检测失败。因此，这 就使得数字水印的安全性受到了限制。所以最好的办法是运用零水印来解决。所谓的零水 印就是不对原始图像或数据进行修改的一种新的数字水印技术。它是利用图像或者数据的 重要特征来构造水印信息，而不是来修改这些特征。所以零水印技术很好地解决了不可见 数字水印的可感知性和鲁棒性之间的矛盾。特别是针对栅格地图的数据要求，由于栅格的 地理空间数据越来越被社会运用，它是描述关于人类赖以生存的地球的重要信息，是国家 基础性、公益性设施，是国家基础设施建设和地球科学研究的支撑性成果，是国家经济、国 防建设中不可缺少的资源，对社会持续发展起着重要作用，因此零水印的应用会对于栅格 地图数据的版权保护起到至关重要的作用。因此零水印非常适用于栅格地图的安全保护， 即不修改原始地图的内容要能满足版权的保护。
技术实现思路
本专利技术针对现有技术的不足，提供了一种基于二维栅格数字地图特征的零水印方法。本专利技术方法包括水印信息嵌入和水印信息检测。水印信息嵌入所述的水印信息嵌入是指提取原始栅格地图的零水印信息，保存入库，具体包括以下 步骤步骤(1)调制水印信息，二值水印图结合密钥和运用Arnold变换把水印信息置乱，得 到二进制比特流w，具体是设定水印图像为NXN的二值图，则水印信息长度为/=NXN，采 用离散化的Arnold变换并结合密钥得到最终水印信息比特流W=(Wl，w2, . . . , Wi, . . . , W1), 其中Wi e {0，1}，0彡Kl ；同时生成两个长度为1的二进制比特流X=(Xl，X2,... ,Xi,..., X1)和y=(yi，y2，...⑶，……)，用来保存特征数据，其中Xi=Yi=O, 0 ^ i<7 ；步骤( 对栅格地图数据进行分块处理根据水印信息的长度Λ把栅格地图数据分成 1 土夬，且每块大小为MXM方形小块；然后进行以下步骤步骤a、对每个MXM的小块数据进行二维sym2小波分解(sym2为Symlets小波系中的 一种)，分解级数为2;步骤b、对第i块小波系数中的LIA和Hffii系数(LIA代表第2级小波分解中的逼近 系数，Hffii代表第2级小波分解中的对角线细节系数)分别进行统计取LI^中的最小值与 最大值，均分成十个等间隔的区间，把LI^中的元素分到此十个区间之内，并返回每个范围 内元素的个数作为一个一维向量，记为V1，再取Hffii中的最小值与最大值，均分成十个等间 隔的区间，把Hffii中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个 一维向量，记为V2;步骤c、分别对向量V1和V2的元素进行从大到小的排序，记为V1 ’和V2'，获得V1'排在向 量前面的两个元素ApB1',以及他们在向量V1中的位置HA1、HB1，获得V2’排在向量前面的两 个元素A2’、B2’，以及他们在向量V2中的位置HA2、HB2，然后做如下处理 如果 abs (Hai-Hbi) >T，则令 Xi=I ； 如果 abs (Hai-Hbi) < T，则令 Xi=O ； 如果 abs (Ha2-Hb2) >T，则令 yi=l ； 如果 abs (Ha2-Hb2) < T，则令 Yi=O ； 其中，abs ()代表取绝对值，T代表阈值。步骤( 生成零水印信息把置乱后的水印信息二进制比特流w和X进行异或操 作得到零水印信息I1, 二进制比特流w和y进行异或操作得到零水印信息I2 ；步骤(4)保存零水印信息把生成的零水印信息I1和I2保存到水印信息库中，以备检 测时需要。水印信息检测所述的水印信息检测是指将待检测栅格地图中提取的零水印信息与水印信息库中的 水印信息操作得到最后水印信息，具体步骤如下步骤⑴根据待检测栅格地图数据，取出水印信息库中的水印信息IJPI2，并生成两个 长度为 7 的二进制比特流 α=(αι，α2，···，Qi,..., CI1)禾口 3=(3i，β2，···，^i,..., β D，用来保存特征数据，其中α严1=0，0彡iO;；步骤( 对待检测栅格地图数据进行分块处理根据水印信息的长度，把栅格地图数 据分成每块大小为MXM的小块；然后进行一下步骤步骤a、对每个MXM的小块数据进行二维sym2小波分解(sym2为Symlets小波系中的 一种)，分解级数为2;步骤b、对第i块小波系数中的LI^和Hffii系数(LIA代表第2级小波分解中的逼近 系数，Hffii代表第2级小波分解中的对角线细节系数)分别进行统计取LI^中的最小值与 最大值，均分成十个等间隔的区间，把LI^中的元素分到此十个区间之内，并返回每个范围 内元素的个数作为一个一维向量，记为V3，再取Hffii中的最小值与最大值，均分成十个等间 隔的区间，把Hffii中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个 一维向量，记为V4;步骤c、分别对向量V3和V4的元素进行从大到小的排序，记为V3’和V4',获得V3'排在向 量前面的两个元素C3'、D3'，以及他们在向量V3中的位置H。3、HD3，获得V4'排在向量前面的两 个元素C4’、D4',以及他们在向量V4中的位置Hc4、Hd4，然后做如下处理5如果 abs (Hg3-Hd3) >T，则令 Qi=I ； 如果 abs (Hc3-Hd3) < T，则令 α i=0 ； 如果 abs (Hc4-HJ >T,则令 ^i=I ； 如果 abs (Hc4-H 4) < T 则令 β 尸0 ； 其中，abs()代农収绝对值，T代表阈值。步骤( 得到检测信息把α和I1进行异或操作，得到W1 ’，把β和I2进行异或 操作，得到W2'；步骤⑷反置乱水印信息把二进制比特流Wp W；结合密钥进行Arnold反变换提取 出最终的水印图像I1 ’和I2'。本专利技术的有益效果为该方法不需要对原始栅格地图数据做任何修改，根据栅格地图数据图像块颜色差别比 较明显的特征，对栅格地图数据做小波变换，并在小波域系数上提取零水印信息。该水印方 法，一方面避免把水印作为冗余信息嵌入在地图数据中引起数据改变的缺点，具有较好的 隐蔽性，可以很好地解决数字水印的不可感知性和鲁棒性之间的矛盾，也可以克服可逆数 字水印中存在的安全漏洞，是一种天然的盲水印系统，有很大的实用价值；另一方面能够抵 抗多种常见的攻击，比如缩放、平移、旋转、裁剪，具有较强的抵抗力和稳健性。具体实施例方式本专利技术方本文档来自技高网...

【技术保护点】
１．一种基于栅格地图特征的小波域零水印方法，其特征在于该方法包括水印信息嵌入和水印信息检测；所述的水印信息嵌入是指提取原始栅格地图的零水印信息，保存入库，具体包括以下步骤：步骤（１）调制水印信息，二值水印图结合密钥和运用Ａｒｎｏｌｄ变换把水印信息置乱，得到二进制比特流ｗ，具体是：设定水印图像为Ｎ×Ｎ的二值图，则水印信息长度为ｌ＝Ｎ×Ｎ，采用离散化的Ａｒｎｏｌｄ　变换并结合密钥得到最终水印信息比特流ｗ＝（ｗ１，ｗ２，．．．，ｗｉ，．．．，ｗｌ），其中ｗｉ∈｛０，１｝，０≤ｉ（ｌ；同时生成两个长度为ｌ的二进制比特流ｘ＝（ｘ１，ｘ２，．．．，ｘｉ，．．．，ｘｌ）和ｙ＝（ｙ１，ｙ２，．．．，ｙｉ，．．．，ｙｌ），用来保存特征数据，其中ｘｉ＝ｙｉ＝０，０≤ｉ（ｌ；步骤（２）对栅格地图数据进行分块处理：根据水印信息的长度ｌ，把栅格地图数据分成ｌ块，且每块大小为Ｍ×Ｍ方形小块；然后进行以下步骤：步骤ａ、对每个Ｍ×Ｍ的小块数据进行二维ｓｙｍ２小波分解，分解级数为２；步骤ｂ、对第ｉ块小波系数中的ＬＬ２ｉ和ＨＨ２ｉ系数分别进行统计：取ＬＬ２ｉ中的最小值与最大值，均分成十个等间隔的区间，把ＬＬ２ｉ中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个一维向量，记为Ｖ１，再取ＨＨ２ｉ中的最小值与最大值，均分成十个等间隔的区间，把ＨＨ２ｉ中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个一维向量，记为Ｖ２；步骤ｃ、分别对向量Ｖ１和Ｖ２的元素进行从大到小的排序，记为Ｖ１’和Ｖ２’，获得Ｖ１’排在向量前面的两个元素Ａ１’、Ｂ１’，以及他们在向量Ｖ１中的位置ＨＡ１、ＨＢ１，获得Ｖ２’排在向量前面的两个元素Ａ２’、Ｂ２’，以及他们在向量Ｖ２中的位置ＨＡ２、ＨＢ２，然后做如下处理：如果ａｂｓ（ＨＡ１－ＨＢ１））Ｔ，则令ｘｉ＝１；如果ａｂｓ（ＨＡ１－ＨＢ１）　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　Ｔ，则令ｘｉ＝０；如果ａｂｓ（ＨＡ２－ＨＢ２））Ｔ，则令ｙｉ＝１；如果ａｂｓ（ＨＡ２－ＨＢ２）Ｔ，则令ｙｉ＝０；其中，ａｂｓ（）代表取绝对值，Ｔ代表阈值；步骤（３）生成零水印信息：把置乱后的水印信息二进制比特流ｗ和ｘ进行异或操作得到零水印信息Ｉ１，二进制比特流ｗ和ｙ进行异或操作得到零水印信息Ｉ２；步骤（４）保存零水印信息：把生成的零水印信息Ｉ１和Ｉ２保存到水印信息库中，以备检测时需要；所述的水印信息检测是指将待检测栅格地图中提取的零水印信息与水印信息库中的水印信息操作得到最后水印信息，具体步骤如下：步骤（５）根据待检测栅格地图数据，取出水印信息库中的水印信息Ｉ１和Ｉ２，并生成两个长度为ｌ的二进制比特流α＝（α１，α２，．．．，αｉ，．．．，αｌ）和β＝（β１，β２，．．．，βｉ，．．．，βｌ），用来保存特征数据，其中αｉ＝βｉ＝０，０≤ｉ（ｌ；；步骤（６）对待检测栅格地图数据进行分块处理：根据水印信息的长度，把栅格地图数据分成每块大小为Ｍ×Ｍ的小块；然后进行以下步骤：步骤ｄ、对每个Ｍ×Ｍ的小块数据进行二维ｓｙｍ２小波分解，分解级数为２；步骤ｅ、对第ｉ块小波系数中的ＬＬ２ｉ和ＨＨ２ｉ系数分别进行统计：取ＬＬ２ｉ中的最小值与最大值，均分成十个等间隔的区间，把ＬＬ２ｉ中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个一维向量，记为Ｖ３，再取ＨＨ２ｉ中的最小值与最大值，均分成十个等间隔的区间，把ＨＨ２ｉ中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个一维向量，记为Ｖ４；步骤ｆ、分别对向量Ｖ３和Ｖ４的元素进行从大到小的排序，记为Ｖ３’和Ｖ４’，获得Ｖ３’排在向量前面的两个元素Ｃ３’、Ｄ３’，以及他们在向量Ｖ３中的位置ＨＣ３、ＨＤ３，获得Ｖ４’排在向量前面的两个元素Ｃ４’、Ｄ４’，以及他们在向量Ｖ４中的位置ＨＣ４、ＨＤ４，然后做如下处理：如果ａｂｓ（ＨＣ３－ＨＤ３））Ｔ，则令αｉ＝１；如果ａｂｓ（ＨＣ３－ＨＤ３）Ｔ，则令αｉ＝０；如果ａｂｓ（ＨＣ４－ＨＤ４））Ｔ，则令βｉ＝１；如果ａｂｓ（ＨＣ４－ＨＤ４）Ｔ，则令βｉ＝０；步骤（７）得到检测信息：把α和Ｉ１进行异或操作，得到ｗ１’，把β和Ｉ２进行异或操作，得到ｗ２’；步骤（８）反置乱水印信息：把二进制比特流ｗ１’、ｗ２’结合密钥进行Ａｒｎｏｌｄ反变换提取出最终的水印图像Ｉ１’和Ｉ２’。...

【技术特征摘要】

【专利技术属性】
技术研发人员：王勋，凌云，章志勇，林传峰，黄定军，
申请(专利权)人：浙江工商大学，
类型：发明
国别省市：86