CN103065278A - 一种针对shp线面图层完整性的多级认证方法 - Google Patents

一种针对shp线面图层完整性的多级认证方法 Download PDF

Info

Publication number
CN103065278A
CN103065278A CN2012105665002A CN201210566500A CN103065278A CN 103065278 A CN103065278 A CN 103065278A CN 2012105665002 A CN2012105665002 A CN 2012105665002A CN 201210566500 A CN201210566500 A CN 201210566500A CN 103065278 A CN103065278 A CN 103065278A
Authority
CN
China
Prior art keywords
data
watermark
watermark information
sequence
classification
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012105665002A
Other languages
English (en)
Other versions
CN103065278B (zh
Inventor
李安波
吴赛松
闾国年
吴雪荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Normal University
Original Assignee
Nanjing Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN201210566500.2A priority Critical patent/CN103065278B/zh
Publication of CN103065278A publication Critical patent/CN103065278A/zh
Application granted granted Critical
Publication of CN103065278B publication Critical patent/CN103065278B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Editing Of Facsimile Originals (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种针对shp线面图层完整性的多级认证方法，属于地理信息完整性认证领域。该方法包括两个过程：（1）水印生成与嵌入过程：参数设置，生成混沌序列，矢量数据分级，矢量数据分组，然后进行各分组数据的水印生成，分组数据水印嵌入；（2）水印的提取与检测过程：参数设置，生成混沌序列，矢量数据分级，矢量数据分组，然后提取各分组数据原始水印信息，获取各分组数据校验水印信息，将各分组数据H i ， j 的原始水印信息与校验水印信息比对，基于比对结果，给出待认证数据的完整性等级。本发明通过数据分级认证和认证链技术，较好地克服了现有技术中的缺点，不但能对shp线面图层的完整性进行篡改定位，而且可以进一步确定篡改的等级和程度。

Description

一种针对shp线面图层完整性的多级认证方法
技术领域
       本发明属于地理信息完整性认证领域，具体涉及一种支持shp线面图层篡改定位和篡改程度识别的完整性多级认证技术。
背景技术
       矢量数字地图具有精度高、无损缩放、易于存储、易于发布的优点。与此同时矢量数字地图也极容易被编辑篡改。因此很多情况下需要对矢量数字地图的完整性进行认证，即明确矢量数字地图是否遭到篡改。当前对矢量数字地图完整性认证主要采用的是数字签名技术和易碎水印技术。其中，数字签名技术虽可以很好的满足对矢量数字地图完整性进行认证的要求，但是存在不能进行篡改定位的缺点；专利《一种基于混沌映射的shp线面图层易碎水印技术》，申请号：201210142896.8；公开号：CN102663673A中虽克服了数字签名技术不能进行篡改定位的缺点，但是存在着只支持单级认证的不足。在当前地理数据安全分级管理的情况下，难以满足数据多级完整性认证的需求。
发明内容
       本发明目的：针对GIS矢量数据的特点，基于数据分级与混沌映射技术，提出一种适用于shp线面图层完整性多级认证方法。
       为了实现上述目的，本发明采取的技术方案：
       一种针对shp线面图层完整性的多级认证方法，包括如下过程：
       （一）水印生成与嵌入过程
       步骤11：参数设置
       设原始线面图层数据为R，用户设置数据分级阈值                                               、数据分组点个数M和混沌密钥
Figure 2012105665002100002DEST_PATH_IMAGE004
，将设置的参数保存密钥文件KeyFile中；
       步骤12：生成混沌序列
       a）生成混沌序列：基于（1）式所示的Logistic方程、设置的参数，生成长度为M的混沌序列，
                                         （1）
       b)二值化混沌序列：根据公式（2）将序列
Figure 2012105665002100002DEST_PATH_IMAGE008
量化得到一维二值混沌位置序列S，其中，
Figure 2012105665002100002DEST_PATH_IMAGE010
，
Figure 2012105665002100002DEST_PATH_IMAGE012
；
                         （2）
       步骤13：矢量数据分级
       对线面要素进行分级，根据用户设置的分级阈值，使用道格拉斯—普克算法，对R中每个线面要素基于（3）式进行分级，其中，d表示道格拉斯—普克算法中数据压缩的极差，经分级后
Figure 2012105665002100002DEST_PATH_IMAGE016
，
Figure 2012105665002100002DEST_PATH_IMAGE018
代表第i级数据集，简称数据集
Figure 132742DEST_PATH_IMAGE018
；
Figure 2012105665002100002DEST_PATH_IMAGE020
                （3）
       步骤14：矢量数据分组
       对经过步骤13分级后的数据进一步进行分组处理，即对每级数据集
Figure 151514DEST_PATH_IMAGE018
按先后顺序进一步进行数据分组，每个分组数据中点个数为M，数据集
Figure 2012105665002100002DEST_PATH_IMAGE022
分组结果为
Figure 2012105665002100002DEST_PATH_IMAGE024
，其中，?count(
Figure 750991DEST_PATH_IMAGE022
) /M?，count()表示
Figure 163835DEST_PATH_IMAGE022
中点个数，??表示向下取整；
      步骤15：各分组数据
Figure 2012105665002100002DEST_PATH_IMAGE028
水印生成
       数据经过分级、分组处理后，进行各分组数据的水印生成，各分组数据
Figure 87929DEST_PATH_IMAGE028
的水印生成方法如下：分组数据
Figure 971440DEST_PATH_IMAGE028
与下一个分组数据
Figure 2012105665002100002DEST_PATH_IMAGE030
组成水印生成链，其中，最后一组数据与第一组数据组成水印生成链，
       a）特征值的计算
       特征值计算选取的是长度因子，根据公式（4）计算
Figure 271971DEST_PATH_IMAGE028
和中每条线面要素的折线长度
Figure 2012105665002100002DEST_PATH_IMAGE032
，
Figure 2012105665002100002DEST_PATH_IMAGE034
表示分组数据
Figure 871897DEST_PATH_IMAGE028
的点坐标，
Figure 2012105665002100002DEST_PATH_IMAGE036
（4）
       b）基于特征值的水印信息生成
       ①将Lj与Lj+1组成认证连Lˊ，
Figure 2012105665002100002DEST_PATH_IMAGE038
，并采用MD5哈希函数对Lˊ进行加密，
       ②对加密结果进行进制转换，MD5哈希函数加密结果为32位十六进制序列，无论是内容上还是长度上都需进一步转化，根据转换后序列长度最接近M的原则，见（5）式，确定转换后的进制，其中，N表示转换后进制；
    N=
Figure 2012105665002100002DEST_PATH_IMAGE040
                              （5）
       步骤16：分组数据
Figure 445967DEST_PATH_IMAGE028
水印嵌入
       对步骤15生成水印信息，先按级别再按分组依次嵌入原始线面图层数据R中，各分组数据
Figure 550189DEST_PATH_IMAGE028
水印嵌入过程如下：
       根据步骤12生成的位置混沌序列S，和公式（6）将水印
Figure 2012105665002100002DEST_PATH_IMAGE042
嵌入
Figure 644047DEST_PATH_IMAGE028
点坐标中，嵌入后记为
Figure 2012105665002100002DEST_PATH_IMAGE044
，
Figure 2012105665002100002DEST_PATH_IMAGE046
表示嵌入前
Figure 909943DEST_PATH_IMAGE028
点坐标；
Figure 2012105665002100002DEST_PATH_IMAGE048
表示嵌入后点坐标，
Figure 2012105665002100002DEST_PATH_IMAGE050
表示水印
Figure 890941DEST_PATH_IMAGE042
序列第i位；
Figure 2012105665002100002DEST_PATH_IMAGE052
             （6）
       （二）水印的提取与检测过程
       步骤21：参数设置
       设待认证线面图层数据为H，读密钥文件KeyFile，获取相关参数设置；
       步骤22：生成混沌序列
       a）生成混沌序列：利用公式（1）的Logistic方程、获取的相关参数，生成长度为M的沌序列，
       b）二值化混沌序列：根据公式（2）将序列
Figure 104884DEST_PATH_IMAGE008
量化得到一维二值混沌位置序列S，
Figure 807261DEST_PATH_IMAGE010
，
Figure 903393DEST_PATH_IMAGE012
；
       步骤23：矢量数据分级
       对线面要素进行分级，根据从KeyFile文件中读取的分级阈值，使用道格拉斯—普克算法，对H中每个线面要素进行分级，基于（3）式，经分级后
Figure 2012105665002100002DEST_PATH_IMAGE054
，
Figure 2012105665002100002DEST_PATH_IMAGE056
代表第i级数据集，简称数据集
Figure 801948DEST_PATH_IMAGE056
；
       步骤24：矢量数据分组
       对经过步骤23分级后的数据进一步进行分组处理，即对每级数据集
Figure 667136DEST_PATH_IMAGE056
按数据的先后顺序进行数据分组，每个分组数据中点个数为M，数据集Hi分组结果为
Figure 2012105665002100002DEST_PATH_IMAGE058
，其中，
Figure 9255DEST_PATH_IMAGE026
?count(Hi) /M?，count(Hi)表示Hi中点个数，??表示向下取整；
      步骤25：提取各分组数据原始水印信息
       原始水印信息
Figure 2012105665002100002DEST_PATH_IMAGE062
是先按级别再按分组依次提取的，各分组数据
Figure 310793DEST_PATH_IMAGE060
的原始水印信息具体方法如下：
       a）根据步骤22生成混沌序列S，获取原始水印信息嵌入位置，
       b）根据公式（7）从j组点坐标中提取原始水印信息
Figure 29350DEST_PATH_IMAGE042
，（
Figure 2012105665002100002DEST_PATH_IMAGE064
）表示j组中第r个点坐标，表示混沌序列S第r位；
Figure 2012105665002100002DEST_PATH_IMAGE068
                 （7）
       步骤26：获取各分组数据校验水印信息
       根据步骤15，先按级别再按分组依次生成各分组数据
Figure 996486DEST_PATH_IMAGE060
校验水印信息
Figure 2012105665002100002DEST_PATH_IMAGE070
 ；
       步骤27：各分组数据
Figure 50898DEST_PATH_IMAGE060
的原始水印信息与校验水印信息比对
       数据的完整性认证方法是先按级别再按分组依次循环比对原始水印信息与校验水印信息是否一致，若原始水印信息与校验水印信息一致，则待认证数据完整；否则，待认证数据遭到篡改，各分组数据
Figure 838726DEST_PATH_IMAGE060
原始水印信息与校验水印信息比对方法如下：
       依次比较分组数据
Figure 412927DEST_PATH_IMAGE060
提取的原始水印信息
Figure 2012105665002100002DEST_PATH_IMAGE072
与获取的校验水印信息
Figure 2012105665002100002DEST_PATH_IMAGE074
是否一致。若
Figure 96849DEST_PATH_IMAGE072
=
Figure 123711DEST_PATH_IMAGE074
，则分级数据集Hi中j组数据完整；若
Figure 698917DEST_PATH_IMAGE072
≠
Figure 127625DEST_PATH_IMAGE074
，且
Figure 2012105665002100002DEST_PATH_IMAGE076
≠
Figure 2012105665002100002DEST_PATH_IMAGE078
则分级数据集Hi遭到篡改，并对篡改进行定位于j+1组数据；
       步骤28：数据完整性分级认证
       基于步骤27的比对结果，依据（8）式给出待认证数据的完整性等级，A级表示数据未遭到篡改，完整性好，可信等级最高，B级表示数据遭到轻微篡改，完整性比A级差；C级表示数据遭到一般篡改，数据完整性比B级差；D级表示数据遭到严重篡改，数据完整性比C级差；E级表示数据遭到非常严重的篡改，数据完整性最差，可信等级最低。
Figure DEST_PATH_IMAGE080
        （8）
本发明通过数据分级认证和认证链技术，较好地克服了现有技术中的缺点，不但能对shp线面图层的完整性进行篡改定位，而且可以进一步确定篡改的等级和程度。
附图说明
       图1是本发明水印生成与嵌入流程图；
       图2是本发明水印提取与检测流程图；
       图3是本发明实施实例采用的实验数据；
       图4是实施实例中原始数据与篡改数据的局部对比图，（a）图为原始数据局部截图，（b）图、（c）图、（d）图为篡改数据的检测结果局部截图。其中，（b）图中点表示遭到篡改的一级点、（c）图中点为遭到篡改的二级点、（d）图中点为遭到篡改的三级点。
具体实施方式
       下面结合附图与实施实例，做进一步的详细说明。
       本实例选择一shp线图层数据，针对数据分级、数据分组、水印生成与嵌入、水印提取与检测的整个过程（面图层与线图层的方法一致），进一步详细说明本发明。本实例选择南京市管线图层数据（如图3所示）作为实验数据。       
（一）水印生成与嵌入流程（如图1）
       步骤11：参数设置
       数据被分为三级，分级阈值分别为1600、600。矢量数据每个数据组包含点个数M设置为100个。依据（5）式确定转换后的进制为三进制。Logistic混沌序列初值Q0设为0.31，混沌系数k设为3.92。并将设置的参数保存至密钥文件KeyFile中。
       步骤12：生成混沌序列
       a）基于（1）式Logistic方程、密钥[k,Q0]、M，生成长度为100的混沌序列。
       b）基于（2）式将序列
Figure 716869DEST_PATH_IMAGE008
量化得到一个二值序列，
Figure DEST_PATH_IMAGE084
。
       步骤13：数据分级
       针对原始线图层R进行分级，根据设定的分级阈值将R分为三级，得到
Figure DEST_PATH_IMAGE086
。
Figure 949136DEST_PATH_IMAGE018
表示第i级数据集，简称数据集
Figure 344345DEST_PATH_IMAGE018
。本实施实例中，数据集
Figure DEST_PATH_IMAGE088
共有1246个点，数据集
Figure DEST_PATH_IMAGE090
中共有779个点，数据集
Figure DEST_PATH_IMAGE092
中共有1379个点。
       步骤14：数据分组
       对经过步骤13分级后的数据集进一步进行分组处理，即对每个数据集Ri进行数据分组。令每个分组数据中点个数为M，得到
Figure DEST_PATH_IMAGE094
，
Figure 565242DEST_PATH_IMAGE028
表示数据集
Figure 387705DEST_PATH_IMAGE022
第j组数据。本实施例中，分为12组，
Figure 306168DEST_PATH_IMAGE090
分为7组，分为13组。
       步骤15：各分组数据
Figure 702831DEST_PATH_IMAGE028
水印生成
       根据数据分级、数据分组情况先按级别再按分组生成水印信息。以数据集
Figure 191582DEST_PATH_IMAGE088
中各分组数据为例说明。
a）  数据集
Figure 663014DEST_PATH_IMAGE088
共有12组数据。根据公式（4）求出每组点列的折线长度，其中，
b） 
Figure DEST_PATH_IMAGE098
。
       b）将相邻两组折线长度
Figure DEST_PATH_IMAGE100
与
Figure DEST_PATH_IMAGE102
组成水印生成链G，最后一组
Figure DEST_PATH_IMAGE104
与第一组
Figure DEST_PATH_IMAGE106
组成水印生成链。
Figure DEST_PATH_IMAGE108
。
       c）采用MD5哈希函数对进行加密，并将加密结果转为三进制整数序列作为i组水印值
Figure DEST_PATH_IMAGE112
。其中，第1组
Figure DEST_PATH_IMAGE114
，第12组
Figure DEST_PATH_IMAGE116
。
       步骤16：各分组数据
Figure 232405DEST_PATH_IMAGE028
水印嵌入
       先按级别再按分组完成水印的嵌入，各分组数据
Figure 396670DEST_PATH_IMAGE028
嵌入方法一致。以分组数据
Figure DEST_PATH_IMAGE118
的水印嵌入为例。根据步骤12生成的位置序列S和公式（6）将该组对应的水印信息嵌入坐标中。
       
Figure DEST_PATH_IMAGE120
，
Figure DEST_PATH_IMAGE122
表示水印嵌入前坐标；
Figure DEST_PATH_IMAGE124
，
Figure DEST_PATH_IMAGE126
表示水印嵌入后坐标。其中，根据位置序列S0=0取第1组第一个点的横坐标
Figure DEST_PATH_IMAGE128
、第1组水印
Figure DEST_PATH_IMAGE130
的第1位
Figure DEST_PATH_IMAGE132
，嵌入后横坐标
Figure DEST_PATH_IMAGE134
。
       （二）水印提取与检测过程（如图2）
        针对上述过程嵌入水印后的线图层数据，随机选取三条线要素进行了适当的点坐标移动编辑，得到受攻击后数据H，以下的水印提取和检测都针对该攻击后数据H进行。
       步骤21：读取密钥文件KeyFile，获取相关参数。根据（1）式Logistic混沌方程、获取的密钥参数[k，Q0]、分组数据包含点个数M，生成位置序列S。其具体方法同所述过程（一）步骤12。
       步骤22：根据分级阈值与所述过程（一）中步骤13的方法，将受攻击后地图数据H分为三级。得
Figure DEST_PATH_IMAGE136
。表示第i级数据集，记为数据集
Figure 435033DEST_PATH_IMAGE138
。
       步骤23：对分级后数据集进一步进行分组，各分组数据中包含点个数M。具体方法同所述过程（一）中步骤14。
Figure DEST_PATH_IMAGE140
分为12组，分为7组，
Figure DEST_PATH_IMAGE144
分为13组。
       步骤24：原始水印信息的提取
       先按级别再按分组依次提取各分组数据原始水印信息，其中，，
Figure DEST_PATH_IMAGE150
，
Figure DEST_PATH_IMAGE152
。各分组数据原始水印信息
Figure DEST_PATH_IMAGE154
提取方法一致。以
Figure DEST_PATH_IMAGE156
中各分组数据原始水印信息提取为例。
       根据位置序列S，循环从数据集
Figure 196370DEST_PATH_IMAGE140
坐标中提取水印信息
Figure 378DEST_PATH_IMAGE156
。此过程为嵌入的逆过程。其中第1组提取水印
Figure DEST_PATH_IMAGE158
，第13组
Figure DEST_PATH_IMAGE160
。具体方法为：
       根据位置序列S，从各分组数据顶点坐标中提取水印，若则从横坐标提取水印否则从纵坐标提取水印。以第1组为例，提取
Figure DEST_PATH_IMAGE164
。
Figure DEST_PATH_IMAGE166
，读取第1组第一个点坐标横坐标
Figure 650671DEST_PATH_IMAGE134
。则第1组原始水印序列第1位为2。
       步骤25：根据所述过程（一）步骤15，先按级别再按组循环生成各分组数据校验水印信息
Figure DEST_PATH_IMAGE168
，其中，，
Figure DEST_PATH_IMAGE172
，
Figure DEST_PATH_IMAGE174
。其中，
Figure 401590DEST_PATH_IMAGE158
，
Figure DEST_PATH_IMAGE176
。
       步骤26：依次比较各级别下各分组数据提取原始水印信息与生成的校验水印信息
Figure DEST_PATH_IMAGE180
。若
Figure 86518DEST_PATH_IMAGE154
=
Figure DEST_PATH_IMAGE182
，则判定数据组j与数据组j+1未遭到篡改。若
Figure 327006DEST_PATH_IMAGE154
≠
Figure 277645DEST_PATH_IMAGE180
且 
Figure DEST_PATH_IMAGE184
≠
Figure DEST_PATH_IMAGE186
，则j+1组数据遭到篡改。
       例如：
Figure DEST_PATH_IMAGE188
≠
Figure DEST_PATH_IMAGE190
，
Figure DEST_PATH_IMAGE192
≠说明第2级数据集中第6组数据遭到篡改。
       步骤27：数据完整性认证。依据（8）式，数据遭到轻微篡改、数据完整等级为C。 
              （三）实验分析
       由上述实例（图4）可知：实例表明该方案不但能对shp线面图层的完整性进行篡改定位，而且可以进一步确定篡改的等级和程度，较好满足了进行shp线面图层的完整性进行多级认证的需求。
       本发明实施实例中仅以shp格式的线图层进行易碎水印的生成、嵌入、提取与检测过程。该方法也可以适用于shp格式的面图层数据。

Claims (1)

1.一种针对shp线面图层完整性的多级认证方法，包括如下过程：
    （一）水印生成与嵌入过程
    步骤11：参数设置
    设原始线面图层数据为R，用户设置数据分级阈值                                               
Figure 2012105665002100001DEST_PATH_IMAGE002
、数据分组点个数M和混沌密钥，将设置的参数保存密钥文件KeyFile中；
    步骤12：生成混沌序列
    a）生成混沌序列：基于（1）式所示的Logistic方程、设置的参数，生成长度为M的混沌序列，
Figure 2012105665002100001DEST_PATH_IMAGE006
                           （1）
    b)二值化混沌序列：根据公式（2）将序列
Figure 2012105665002100001DEST_PATH_IMAGE008
量化得到一维二值混沌位置序列S，其中，
Figure 2012105665002100001DEST_PATH_IMAGE010
，
Figure 2012105665002100001DEST_PATH_IMAGE012
；
Figure DEST_PATH_IMAGE014
                         （2）
    步骤13：矢量数据分级
    对线面要素进行分级，根据用户设置的分级阈值，使用道格拉斯—普克算法，对R中每个线面要素基于（3）式进行分级，其中，d表示道格拉斯—普克算法中数据压缩的极差，经分级后
Figure DEST_PATH_IMAGE016
，
Figure DEST_PATH_IMAGE018
代表第i级数据集，简称数据集
Figure 100506DEST_PATH_IMAGE018
；
Figure DEST_PATH_IMAGE020
                （3）
    步骤14：矢量数据分组
    对经过步骤13分级后的数据进一步进行分组处理，即对每级数据集
Figure 220909DEST_PATH_IMAGE018
按先后顺序进一步进行数据分组，每个分组数据中点个数为M，数据集
Figure DEST_PATH_IMAGE022
分组结果为
Figure DEST_PATH_IMAGE024
，其中，
Figure DEST_PATH_IMAGE026
?count(
Figure 922017DEST_PATH_IMAGE022
) /M?，count(
Figure 51647DEST_PATH_IMAGE022
)表示
Figure 334861DEST_PATH_IMAGE022
中点个数，??表示向下取整；
      步骤15：各分组数据
Figure DEST_PATH_IMAGE028
水印生成
    数据经过分级、分组处理后，进行各分组数据的水印生成，各分组数据
Figure 875433DEST_PATH_IMAGE028
的水印生成方法如下：分组数据与下一个分组数据
Figure DEST_PATH_IMAGE030
组成水印生成链，其中，最后一组数据与第一组数据组成水印生成链，
       a）特征值的计算
    特征值计算选取的是长度因子，根据公式（4）计算
Figure 279049DEST_PATH_IMAGE028
和
Figure 682349DEST_PATH_IMAGE030
中每条线要素的折线长度
Figure DEST_PATH_IMAGE032
，
Figure DEST_PATH_IMAGE034
表示分组数据
Figure 128243DEST_PATH_IMAGE028
的点坐标，
（4）
    b）基于特征值的水印信息生成
    ①将Lj与Lj+1组成认证连Lˊ，
Figure DEST_PATH_IMAGE038
，并采用MD5哈希函数对Lˊ进行加密，
    ②对加密结果进行进制转换，根据转换后序列长度最接近M的原则，见（5）式，确定转换后的进制，其中，N表示转换后进制；
    N=
Figure DEST_PATH_IMAGE040
                              （5）
    步骤16：分组数据
Figure 554676DEST_PATH_IMAGE028
水印嵌入
    对步骤15生成水印信息，先按级别再按分组依次嵌入原始线面图层数据R中，各分组数据
Figure 26109DEST_PATH_IMAGE028
水印嵌入过程如下：
       根据步骤12生成的位置混沌序列S，和公式（6）将水印嵌入
Figure 267603DEST_PATH_IMAGE028
点坐标中，嵌入后记为，
Figure DEST_PATH_IMAGE046
表示嵌入前
Figure 900710DEST_PATH_IMAGE028
点坐标；
Figure DEST_PATH_IMAGE048
表示嵌入后
Figure 611177DEST_PATH_IMAGE028
点坐标，
Figure DEST_PATH_IMAGE050
表示水印
Figure 355142DEST_PATH_IMAGE042
序列第i位；
Figure DEST_PATH_IMAGE052
             （6）
    （二）水印的提取与检测过程
    步骤21：参数设置
    设待认证线面图层数据为H，读密钥文件KeyFile，获取相关参数设置；
    步骤22：生成混沌序列
    a）生成混沌序列：利用公式（1）的Logistic方程、获取的相关参数，生成长度为M的沌序列，
    b）二值化混沌序列：根据公式（2）将序列
Figure 716722DEST_PATH_IMAGE008
量化得到一维二值混沌位置序列S，，
Figure 249651DEST_PATH_IMAGE012
；
    步骤23：矢量数据分级
    对线面要素进行分级，根据从KeyFile文件中读取的分级阈值，使用道格拉斯—普克算法，对H中每个线面要素进行分级，基于（3）式，经分级后
Figure DEST_PATH_IMAGE054
，
Figure DEST_PATH_IMAGE056
代表第i级数据集，简称数据集；
    步骤24：矢量数据分组
    对经过步骤23分级后的数据进一步进行分组处理，即对每级数据集
Figure 747815DEST_PATH_IMAGE056
按数据的先后顺序进行数据分组，每个分组数据中点个数为M，数据集Hi分组结果为
Figure DEST_PATH_IMAGE058
，其中，?count(Hi) /M?，count(Hi)表示Hi中点个数，??表示向下取整；
      步骤25：提取各分组数据
Figure DEST_PATH_IMAGE060
原始水印信息
    原始水印信息是先按级别再按分组依次提取的，各分组数据的原始水印信息具体方法如下：
    a）根据步骤22生成混沌序列S，获取原始水印信息嵌入位置，
    b）根据公式（7）从j组点坐标中提取原始水印信息
Figure 488958DEST_PATH_IMAGE042
，（）表示j组中第r个点坐标，
Figure DEST_PATH_IMAGE066
表示混沌序列S第r位；
                 （7）
    步骤26：获取各分组数据
Figure 779125DEST_PATH_IMAGE060
校验水印信息
    根据步骤15，先按级别再按分组依次生成各分组数据
Figure 456094DEST_PATH_IMAGE060
校验水印信息
Figure DEST_PATH_IMAGE070
 ；
    步骤27：各分组数据
Figure 346558DEST_PATH_IMAGE060
的原始水印信息与校验水印信息比对
比对方法如下：依次比较分组数据
Figure 501596DEST_PATH_IMAGE060
提取的原始水印信息
Figure DEST_PATH_IMAGE072
与获取的校验水印信息
Figure DEST_PATH_IMAGE074
是否一致，若=
Figure 759719DEST_PATH_IMAGE074
，则分级数据集Hi中j组数据完整；若
Figure 419371DEST_PATH_IMAGE072
≠
Figure 643679DEST_PATH_IMAGE074
，且
Figure DEST_PATH_IMAGE076
≠
Figure DEST_PATH_IMAGE078
则分级数据集Hi遭到篡改，并对篡改进行定位于j+1组数据；
    步骤28：数据完整性分级认证
    基于步骤27的比对结果，给出待认证数据的完整性等级。
CN201210566500.2A 2012-12-24 2012-12-24 一种针对shp线面图层完整性的多级认证方法 Expired - Fee Related CN103065278B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210566500.2A CN103065278B (zh) 2012-12-24 2012-12-24 一种针对shp线面图层完整性的多级认证方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210566500.2A CN103065278B (zh) 2012-12-24 2012-12-24 一种针对shp线面图层完整性的多级认证方法

Publications (2)

Publication Number Publication Date
CN103065278A true CN103065278A (zh) 2013-04-24
CN103065278B CN103065278B (zh) 2015-08-05

Family

ID=48107895

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210566500.2A Expired - Fee Related CN103065278B (zh) 2012-12-24 2012-12-24 一种针对shp线面图层完整性的多级认证方法

Country Status (1)

Country Link
CN (1) CN103065278B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103559678A (zh) * 2013-10-30 2014-02-05 南京师范大学 一种shp线面图层数据的置乱与还原方法
CN104050622A (zh) * 2014-06-02 2014-09-17 鲁东大学 基于三进制编码的彩色图像盲水印方法
CN104462886A (zh) * 2014-11-28 2015-03-25 重庆市地理信息中心 一种基于矢量空间数据对象存储顺序的数字水印方法
CN104751065A (zh) * 2015-03-09 2015-07-01 哈尔滨工程大学 一种基于地理信息产品文件的加密方法
CN105550971A (zh) * 2015-12-15 2016-05-04 南京师范大学 一种基于距离比值调制的矢量地理数据密级标识方法
CN108053358A (zh) * 2017-11-22 2018-05-18 明鉴方寸（北京）科技有限公司 一种水印图生成方法、装置以及存储设备
CN112308755A (zh) * 2020-10-30 2021-02-02 平安普惠企业管理有限公司 基于Hybrid混合型App页面水印的图片检测方法、系统及装置
CN114647824A (zh) * 2022-05-23 2022-06-21 南京信息工程大学 一种神经网络主动保护方法、系统、存储介质及计算设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7171561B2 (en) * 2002-10-17 2007-01-30 The United States Of America As Represented By The Secretary Of The Air Force Method and apparatus for detecting and extracting fileprints
CN102184244A (zh) * 2011-05-18 2011-09-14 武汉大学 一种适合于polygon类型、高鲁棒的地理数据库水印方法
CN102663673A (zh) * 2012-05-10 2012-09-12 南京师范大学 一种基于混沌映射的shp线面图层易碎水印技术

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7171561B2 (en) * 2002-10-17 2007-01-30 The United States Of America As Represented By The Secretary Of The Air Force Method and apparatus for detecting and extracting fileprints
CN102184244A (zh) * 2011-05-18 2011-09-14 武汉大学 一种适合于polygon类型、高鲁棒的地理数据库水印方法
CN102663673A (zh) * 2012-05-10 2012-09-12 南京师范大学 一种基于混沌映射的shp线面图层易碎水印技术

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李安波等: ""空间数据文件的实时版权保护"", 《地球信息科学学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103559678A (zh) * 2013-10-30 2014-02-05 南京师范大学 一种shp线面图层数据的置乱与还原方法
CN103559678B (zh) * 2013-10-30 2016-05-04 南京师范大学 一种shp线面图层数据的置乱与还原方法
CN104050622A (zh) * 2014-06-02 2014-09-17 鲁东大学 基于三进制编码的彩色图像盲水印方法
CN104462886A (zh) * 2014-11-28 2015-03-25 重庆市地理信息中心 一种基于矢量空间数据对象存储顺序的数字水印方法
CN104462886B (zh) * 2014-11-28 2017-10-31 重庆市地理信息中心 一种基于矢量空间数据对象存储顺序的数字水印方法
CN104751065A (zh) * 2015-03-09 2015-07-01 哈尔滨工程大学 一种基于地理信息产品文件的加密方法
CN105550971A (zh) * 2015-12-15 2016-05-04 南京师范大学 一种基于距离比值调制的矢量地理数据密级标识方法
CN105550971B (zh) * 2015-12-15 2018-07-20 南京师范大学 一种基于距离比值调制的矢量地理数据密级标识方法
CN108053358A (zh) * 2017-11-22 2018-05-18 明鉴方寸（北京）科技有限公司 一种水印图生成方法、装置以及存储设备
CN112308755A (zh) * 2020-10-30 2021-02-02 平安普惠企业管理有限公司 基于Hybrid混合型App页面水印的图片检测方法、系统及装置
CN114647824A (zh) * 2022-05-23 2022-06-21 南京信息工程大学 一种神经网络主动保护方法、系统、存储介质及计算设备

Also Published As

Publication number Publication date
CN103065278B (zh) 2015-08-05

Similar Documents

Publication Publication Date Title
CN103065278A (zh) 一种针对shp线面图层完整性的多级认证方法
CN101458810B (zh) 一种基于对象属性特征的矢量地图水印方法
CN102110281B (zh) 一种针对shape点图层数据的差值扩展水印嵌入与提取方法
CN108648133B (zh) 一种结合块旋转和马赛克的无嵌入伪装方法
Wang et al. Reversible fragile watermarking for locating tampered blocks in 2D vector maps
CN104636764B (zh) 一种图像隐写分析方法以及其装置
CN104036531B (zh) 一种基于矢量量化以及二叉查找树的信息隐藏方法
CN113129198B (zh) 一种零水印生成方法、系统及版权的侵权比对方法、系统
CN103747271A (zh) 一种基于混合感知哈希的视频篡改检测方法和装置
CN104063731A (zh) 一种采用数字水印技术的二维码防伪印刷及验证方法
CN103325082A (zh) 一种基于lsd平面的矢量地图可逆信息隐藏方法
CN105488434A (zh) 一种基于标记的矢量地图完整性认证方法
CN102184244B (zh) 一种适合于polygon类型、高鲁棒的地理数据库水印方法
CN103986721A (zh) 一种面向ip语音的最低有效位隐写方法
CN110430334B (zh) 一种基于扇形构图单元叠加的马赛克拼图伪装及恢复方法
Swain A data hiding technique by mixing MFPVD and LSB substitution in a pixel
Neyman et al. Reversible fragile watermarking based on difference expansion using manhattan distances for 2d vector map
CN103903217A (zh) 一种基于顶点插入的矢量地图完整性认证方法
CN113902135A (zh) 样本数据量联合扩充方法、装置、设备、系统及存储介质
CN103236265A (zh) 一种针对MP3Stegz的隐写检测方法
CN108876693A (zh) 一种基于角度差分的矢量地理空间数据数字水印方法
CN102663673B (zh) 一种基于混沌映射的shp线面图层易碎水印方法
CN102509058B (zh) 基于冗余位替换的点类型gis矢量数据伪装与还原方法
Cao et al. Iterative embedding-based reversible watermarking for 2D-vector maps
CN106599727A (zh) 一种矢量地图完整性认证方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150805

Termination date: 20171224