一种GIS矢量数据图像水印的相似度计算方法技术

本发明专利技术公开了一种GIS矢量数据图像水印的相似度计算方法，属于地理信息版权保护领域。该方法先利用原始水印W对提取的水印W′进行位置校正，使得发生错乱的版权图像的像素回到正确的位置，然后将校正好的水印W″与原始水印W采用字符串近似匹配的动态规划算法进行相似度计算。本发明专利技术的方法能够正确地将提取的图像水印的像素较正到正确的位置，直观反映出数据被篡改的位置，并客观地度量原始水印与提取水印的相似度，在一定程度上提高了水印认证的质量，降低了水印认证的漏检率，完善了地理信息版权保护的理论与方法体系，可应用于GIS矢量数据的版权保护技术与安全传输方面。

【技术实现步骤摘要】

本专利技术属于地理信息版权保护领域，具体涉及一种针对GIS矢量数据图像水印的位置校正与相似度计算方法。
技术介绍
版权保护是版权标记基础上的版权认证、内容保密及访问控制。近年来，以数字水印为代表的版权标记技术日趋成为地理信息版权保护研究与应用的重点。国内外关于如何在矢量数据中嵌入水印已经有了比较成熟的研究，基础算法理论方面已经有了相当数量的文献可供查阅。与“属性明显，位置隐含”的栅格数据相比，矢量数据具有其自身的个性特征“位置明显，属性隐含”、无固定存储顺序的数据组织、多样的数据格式、复杂的投影变换、繁多的空间分析应用，以及精度高、冗余少等。正是由于矢量数据的诸多特性，使得从其中提取的水印信息容易发生移位和错乱，而在水印认证过程中，若采用传统的水印相似度计算方法会造成较高的漏检率，因此有必要提出一种更加合理的水印相似度计算方法。
技术实现思路
本专利技术的目的在于根据GIS矢量数据的特点，提出一种针对GIS矢量数据图像水印的相似度计算方法，从而有效提高版权认证的质量，降低漏检率。为了实现上述目的，本专利技术采取的技术方案主要包括以下两个过程(1)水印位置校正的过程步骤一原始水印W长为Lw，通过嵌入算法将水印W嵌入载体数据D中，对D进行攻击获得含水印的数据D'，从D'中提取出水印W'长为Lw,，记Loc = O;步骤二 通过公式(1)计算两水印W和W'的长度差d ；初始化要校正的水印序列 W"，长度为Lw，W〃 (i) =' ，其中'表示的是被修改的位置，i =0，1，2，-,Lff-I ；d = Lff-Lr I (1)步骤三(a)当原始水印长度大于或等于提取的水印长度时，即Lw > Lw,时i)从提取的水印W'中读取字符W' ⑴，i = 0，1，2，…，Lw-I ；ii)从原始水印W中，从第Loc位开始取d+Ι个字符，若Loc+d+1 > Lw时，取Lw-Loc 个字符，记为W(dis)；iii)将W' (i)与W(dis)进行匹配，若匹配成功，则记录匹配起始位置Pos，进入步骤iv)；若匹配不成功，则i = i+1，Loc = Loc+1，进入步骤ν)；iv)校正水印序列 W"第 Loc+Pos-1 位为W〃(Loc+Pos-1) = W' (i)，i = i+1, Loc = Loc+Pos ；ν)循环i)、iv)的过程，直至原始水印W或提取的水印Ψ读取完成，保存校正后的水印序列W"，此时序列W"是将提取的水印序列W'中的图像像素校正到正确位置后的水印，并且其长度与水印W的长度相等；(b)当原始水印长度小于提取水印的长度时，即Lw < Lr时i)从原始的水印W中读取字符W(i)，i = 0，1，2，…，Lw-I;ii)从提取的水印W'中从第Loc位开始取d+Ι个字符，若Loc+d+1 >LW,时，则取 Lw, -Loc 个字符，记为 W' (dis)；iii)将W(i)与W' (dis)进行匹配，若匹配成功，则记录配起始位置Pos，进入步骤iv)；若匹配不成功则i = i+1，Loc = Loc+1，进入步骤ν)；iv)校正的水印序歹Ij W〃第 i 位为=W" (i) = ff(i), i = i+1, Loc = Loc+Pos ；ν)循环步骤i)至iv)的过程，直至原始水印W或提取的水印W'读取完成，保存校正后的水印序列W"，此时序列W"是将提取水印序列W'中的图像像素校正到正确位置后的水印，并且其长度与水印W的长度相等；(2)相似度计算过程步骤一令校正后的水印序列W"为模式串，原始水印W为匹配串，根据公式(2) 生成距离矩阵 M，其中，i = 0,1,2, · · ·，Lw;j = 0,1,2, ...，Lw;Mi, — i，Mj, —O「00241 M ΛiI(WV) = W(J)) (2)步骤二 设f =Mlw0 ； K表示从校正后的水印序列W"转换为原始水印W的最少编辑次数，当K = O时，表示序列W"与原始水印W完全相同；步骤三读取Mw，判断K与从^的大小，当K>MW时，K=Mw，j = j+l ；步骤四当j < Lw时，循环步骤三，否则进入步骤五；步骤五计算原始水印W与还原后的水印序列W"的相似度Sim= (1-α)*100%, 其中，比值α =1(/1 表示错误水平，通过它可以对模式串W"容忍错误的能力进行度量。本专利技术根据GIS矢量数据的特点，所采用的方法能成功地应用于GIS矢量数据图像水印的校正及相似度计算，一定程度上提高了水印认证的质量，降低了水印认证的漏检率，完善了地理信息版权保护的理论与方法体系，可应用于GIS矢量数据的版权保护技术与安全传输方面。本专利技术的方法同样适用于文本水印。附图说明图1是本专利技术实施例选用的实验数据。图2是本专利技术实施例选用的原始版权标记图像。图3是本专利技术方法的水印位置校正流程图(a)为原始水印长度大于或等于提取水印长度情况下的校正；(b)为原始水印长度小于提取水印长度情况下的校正。图4是本专利技术方法的水印相似度计算流程图。图5是本专利技术实施例中水印位置的校正和相似度计算效果图。具体实施例方式本专利技术采用如下技术思路先利用原始水印W对提取的水印W'进行位置校正，使得发生错乱的版权图像的像素回到正确的位置；将校正好的水印W"与原始水印W采用字符串近似匹配的动态规划算法进行相似度计算。本实例选择某一图像水印，针对图像水印的生成、水印位置的校正、水印的相似性5计算的整个过程，进一步详细说明本专利技术。本实例选择江苏省1 25的shp面图层数据 (如图1)作为载体数据，选择二值图像标志m(如图2)作为版权信息。(1)水印位置校正的过程步骤一原始水印W= “10001110000011100101. · · ”，长为Lw = 6400，通过“最低有效位嵌入算法”嵌入载体数据D中，对D进行删除10%顶点的攻击，获得含水印的数据D'， 从 D'中提取出水印 W' = “10011000001110010110. .· ” 长为 Lw, = 5760，记 Loc = 0 ；步骤二 通过公式(1)计算两水印的长度差d = 640 ；初始化要校正的序列W"，长度为Lw，W〃 (i) =' (1=0，1，2，...，1^-1)，其中'表示的是被修改的位置；步骤三因Lw = 6400，Lr = 5760，满足Lff > Lr条件，执行以下过程(如图3 (a) 所示)a)从提取的水印 W'中读取字符 W' (i)(i =0，1，2，...，LW_1)；b)从原始水印W中，从第Loc位开始取d+Ι个字符，若Loc+641 > Lw时，取Lw-Loc 个字符，记为W (di S)。c)将W' (0)与W(dis)进行匹配，若匹配成功，则记录匹配起始位置Pos = 1，进入步骤d)；d)校正的水印 W"第 Loc+Pos-1 位为(0) = W' (0)，i = 1，Loc = Loc+Pos =1 ；e)循环a)、d)的过程，直至原始水印W读取完或提取的水印W'读取完成，保存校正后的序列W" =“100*11*000001110010110...，，。此时序列W"表示的是将W'中的图像像素校正到正确位置后的水印，它与W的长度相等。(2)本文档来自技高网...

【技术保护点】
１．一种ＧＩＳ矢量数据图像水印的相似度计算方法，包括以下过程：（１）水印位置校正的过程：步骤一：原始水印Ｗ长为ＬＷ，通过嵌入算法将水印Ｗ嵌入载体数据Ｄ中，对Ｄ进行攻击获得含水印的数据Ｄ′，从Ｄ′中提取出水印Ｗ′长为ＬＷ′，记Ｌｏｃ＝０；步骤二：通过公式（１）计算两水印Ｗ和Ｗ′的长度差ｄ；初始化要校正的水印序列Ｗ″，长度为ＬＷ，Ｗ″（ｉ）＝′＊′，其中′＊′表示的是被修改的位置，ｉ＝０，１，２，…，ＬＷ－１；ｄ＝｜ＬＷ－ＬＷ′｜　　　　（１）步骤三：（ａ）当原始水印长度大于或等于提取的水印长度时，即ＬＷ≥ＬＷ′时：ｉ）从提取的水印Ｗ′中读取字符Ｗ′（ｉ），ｉ＝０，１，２，…，ＬＷ－１；ｉｉ）从原始水印Ｗ中，从第Ｌｏｃ位开始取ｄ＋１个字符，若Ｌｏｃ＋ｄ＋１＞ＬＷ时，取ＬＷ－Ｌｏｃ个字符，记为Ｗ（ｄｉｓ）；ｉｉｉ）将Ｗ′（ｉ）与Ｗ（ｄｉｓ）进行匹配，若匹配成功，则记录匹配起始位置Ｐｏｓ，进入步骤ｉｖ）；若匹配不成功，则ｉ＝ｉ＋１，Ｌｏｃ＝Ｌｏｃ＋１，进入步骤ｖ）；ｉｖ）校正水印序列Ｗ″第Ｌｏｃ＋Ｐｏｓ－１位为：Ｗ″（Ｌｏｃ＋Ｐｏｓ－１）＝Ｗ′（ｉ），ｉ＝ｉ＋１，Ｌｏｃ＝Ｌｏｃ＋Ｐｏｓ；ｖ）循环ｉ）、ｉｖ）的过程，直至原始水印Ｗ或提取的水印Ｗ′读取完成，保存校正后的水印序列Ｗ″，此时序列Ｗ″是将提取的水印序列Ｗ′中的图像像素校正到正确位置后的水印，并且其长度与水印Ｗ的长度相等；（ｂ）当原始水印长度小于提取水印的长度时，即ＬＷ＜ＬＷ′时：ｉ）从原始的水印Ｗ中读取字符Ｗ（ｉ），ｉ＝０，１，２，…，ＬＷ－１；ｉｉ）从提取的水印Ｗ′中从第Ｌｏｃ位开始取ｄ＋１个字符，若Ｌｏｃ＋ｄ＋１＞ＬＷ′时，则取ＬＷ′－Ｌｏｃ个字符，记为Ｗ′（ｄｉｓ）；ｉｉｉ）将Ｗ（ｉ）与Ｗ′（ｄｉｓ）进行匹配，若匹配成功，则记录配起始位置Ｐｏｓ，进入步骤ｉｖ）；若匹配不成功则ｉ＝ｉ＋１，Ｌｏｃ＝Ｌｏｃ＋１，进入步骤ｖ）；ｉｖ）校正的水印序列Ｗ″第ｉ位为：Ｗ″（ｉ）＝Ｗ（ｉ），ｉ＝ｉ＋１，Ｌｏｃ＝Ｌｏｃ＋Ｐｏｓ；ｖ）循环步骤ｉ）至ｉｖ）的过程，直至原始水印Ｗ或提取的水印Ｗ′读取完成，保存校正后的水印序列Ｗ″，此时序列Ｗ″是将提取的水印序列Ｗ′中的图像像素校正到正确位置后的水印，并且其长度与水印Ｗ的长度相等；（２）相似度计算过程：步骤一：令校正后的水印序列Ｗ″为模式串，原始水印Ｗ为匹配串，根据公式（２）生成距离矩阵Ｍ，其中，ｉ＝０，１，２，．．．，ＬＷ；ｊ＝０，１，２，．．．，ＬＷ；Ｍｉ，０←ｉ，Ｍｊ，０←０（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｓｕｂ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｓｕｂ）?（ｍｏ）＆ＬｅｆｔＡｒｒｏｗ；（／ｍｏ）?（ｍｆｅｎｃｅｄｏｐｅｎ＝＇｛＇ｃｌｏｓｅ＝＇＇）?（ｍｔａｂｌｅ）?（ｍｔｒ）?（ｍｔｄ）?（ｍｓｕｂ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（／ｍｓｕｂ）?（／ｍｔｄ）?（ｍｔｄ）?（ｍｉ）ｉｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｐ）?（ｍｉ）Ｗ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）＆ｐｒｉｍｅ；（／ｍｏ）?（ｍｏ）＆ｐｒｉｍｅ；（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｓｕｐ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＝（／ｍｏ）?（ｍｉ）Ｗ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｔｄ）?（／ｍｔｒ）?（ｍｔｒ）?（ｍｔｄ）?（ｍｎ）１（／ｍｎ）?（ｍｏ）＋（／ｍｏ）?（ｍｉ）ｍｉｎ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（／ｍｓｕｂ）?（ｍｏ），（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｓｕｂ）?（ｍｏ），（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｊ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?...

【技术特征摘要】

【专利技术属性】
技术研发人员：李莎莎，李安波，周卫，闾国年，
申请(专利权)人：南京师范大学，
类型：发明
国别省市：84