测绘学报 | 奚旭：基于区间映射和最大扰动区域的矢量地图可逆水印算法

原创测绘学报智绘科服

本文内容来源于《测绘学报》2022年第11期（审图号GS京（2022）1229号）

基于区间映射和最大扰动区域的矢量地图可逆水印算法

奚旭¹, 张新长²

1. 苏州科技大学地理科学与测绘工程学院，江苏苏州 215009;
2. 广州大学地理科学与遥感学院，广东广州 510006

基金项目：国家自然科学基金(42101420;42071441)

摘要：针对传统矢量地图可逆水印算法往往因片面追求水印容量而容易暴露水印的问题，本文设计了一种扰动程度可控的矢量地图可逆水印算法。首先，基于量化索引调制的思想设计了坐标点的状态区间轴，利用水印信息调制坐标点的状态值实现水印的嵌入；然后，利用最大扰动区域方法限定坐标点状态区间的范围，确保水印信息的嵌入对数据的扰动控制在合理范围内，不会对数据质量造成较大影响；最后，将二进制水印信息十进制化，利用十进制数在最大扰动区域范围内调制坐标点的状态值，实现水印信息的大容量嵌入。试验结果表明，本文提出的矢量地图可逆水印算法不仅容量大，而且扰动程度可控，对常规几何攻击、坐标点攻击等具备较强的稳健性，能够较好地处理水印容量、不可见性及稳健性之间的相互制约关系。

关键词：矢量地图可逆水印区间映射最大扰动区域

引文格式：奚旭, 张新长. 基于区间映射和最大扰动区域的矢量地图可逆水印算法[J]. 测绘学报，2022，51(11)：2379-2389. DOI: 10.11947/j.AGCS.2022.20210552 复制到剪切板

XI Xu, ZHANG Xinchang. Reversible watermarking for vector maps based on interval mapping and maximum perturbation region[J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(11): 2379-2389. DOI: 10.11947/j.AGCS.2022.20210552 复制到剪切板

阅读全文：http://xb.chinasmp.com/article/2022/1001-1595/20221114.htm

引言

矢量地图具有极高的社会经济价值，在数字经济、智慧城市、虚拟现实等新兴发展理念中具有重要基础功能^[1^-3^]。在当前大数据时代，矢量地图的非法传播、倒卖与篡改等行为越来越猖狂，给数据合法拥有方造成了难以估量的损失^[4^-5^]。数字水印^[6^-7^]作为信息安全领域的前沿技术被引用到矢量地图中，为地理空间数据的版权保护提供了技术支持。

目前，矢量地图水印的研究成果已经非常丰富，根据传统的嵌入方式主要分为空间域水印算法和频率域水印算法^[5^, 8^]，这两种算法通过直接或间接修改坐标值的方式实现水印嵌入。由于水印信息的嵌入会对坐标点造成扰动，因此，在设计水印方案的同时设计数据质量控制措施是矢量地图数字水印研究领域的一个重要方向^[9^]。文献[10]用相对位置精度和绝对位置精度作为约束条件，在线要素的周长和面要素的面积中嵌入水印，确保坐标点的变化在允许容差范围内。文献[11]利用多边形相似变换，基于量化索引调制的方法，将水印信息调制于多边形要素的空间关系的特征集合中，有效保持了要素的几何形状。文献[12]将水印嵌入至扩大倍数离散傅里叶变换的(discrete Fourier transform，DFT)相位和幅度上，在一定程度上减少了数据扰动等。以上研究均设计了质量控制方法，但仍然会在不同程度上降低数据保真度，因此，具有无损性质的水印算法越来越受关注。

可逆水印^[13^]算法在正确提取水印的同时可以恢复原始数据，保证数据质量不受影响，具有无损特性，是矢量地图数字水印研究的一个热点方向^[14^-16^]。文献[17]基于离散余弦变换(discrete cosine transform，DCT)设计了矢量地图的可逆水印算法，这是可逆水印首次在矢量地图中的应用，但该方法是基于频率域变换实现的水印嵌入，在水印容量、可逆性及不可见性方面的表现都存在不足。为了提升可逆水印算法在矢量地图中的应用能力，国内外学者在提升算法性能方面做出了很多工作。文献[18—19]分别设计了基于差值扩张的矢量地图可逆水印算法，通过修改相邻顶点间的距离值实现水印嵌入，这类方法在水印容量方面表现出较好的性能，但嵌入水印时未考虑要素的几何特征，导致存在较多的要素形态扭曲情况。文献[20]对差值扩张算法进行了改进，通过修改顶点坐标与线要素的重心坐标的差值实现水印的嵌入，在提升水印容量方面效果显著，但数据的扰动情况并未得到改善。文献[21]以单个要素中特征点集的平均坐标为水印嵌入域，设计了一种可迭代嵌入的矢量地图可逆水印方案，该方法显著提升了水印容量，但迭代嵌入对要素的扰动情况介绍较少。文献[22—24]将地理坐标转换为极坐标，以极坐标为嵌入域设计了多种矢量地图可逆水印算法，并以QR码为水印信息载体，通过提升荷载能力的方式实现了水印信息的大容量嵌入，且采用QR码可减少对数据的扰动，但在不可见性的定量控制方面仍关注较少。综上所述，提升水印容量是矢量地图可逆水印算法研究的一个热点方向，水印容量的提升可以增大可逆水印的恢复提取概率，但也意味着载体数据需要牺牲更多的数据精度来提供水印嵌入空间。虽然可逆水印算法可以在提取水印后实现数据恢复，但过度扰动不仅会影响数据的可用性，也可能会暴露水印信息的存在，从而降低算法的安全性。因此，关注可逆算法下数据的扰动问题，设计科学合理的定量化可控扰动措施，对于提升矢量地图可逆水印算法的实用性具有重要意义。

针对上述问题，本文通过建立坐标点的状态区间轴，以量化索引调制的思想调制坐标点的状态值，将水印信息映射到状态区间中，状态区间根据矢量地图中坐标点的最大扰动区域(maximum perturbation region, MPR)进行范围限制，从而控制坐标点的合理移动范围，实现扰动程度可控的矢量地图可逆水印算法。

1 可控扰动的区间映射可逆水印算法

1.1 数据预处理方法

1.1.1 矢量地图预处理

数据预处理的主要目的是提取在预处理操作下仍然稳健的顶点，将这些顶点用于水印信息的嵌入域，可提升水印方案的整体稳健性。本文的数据预处理环节主要包括特征点的提取与重排序两部分。为了增强水印算法抵抗数据压缩、坐标点攻击的能力，在嵌入水印前，采用Douglas-Peucker算法^[25^]对矢量地图进行压缩简化，设置压缩阈值τ，提取简化后的坐标点p_i(i=1, 2, …, n)作为待嵌入水印的特征点。得到这些特征点后进行重排序，使基于该排序的水印算法获得抵抗数据重排序、旋转等攻击的能力，基本步骤如下：①计算所有特征点到地图几何中心的距离d_i(i=1, 2, …，n)；②根据该距离从小到大的顺序对特征点进行重排序，获取重排序后的特征点集p_i^re(i=1, 2, …, n)，后续将基于该点集进行水印嵌入处理。

1.1.2 水印信息预处理

本文提出的水印算法依赖水印信息调制坐标点在状态区间范围内的状态值，从而实现水印信息的嵌入。为了提升水印容量，将二进制流水印信息转为十进制数调制坐标状态值，利用十进制数较大的信息承载能力提升水印容量。设待嵌入的水印是长度为l的二进制流信息W^binary={w_i^binary|w_i^binary∈{0, 1}, i∈[0, l-1]}，设c为每个坐标值携带的水印位，n_w为待嵌入水印的坐标值个数，则需要保证l=n_w×c。坐标点的x、y坐标同时携带水印，则嵌入长度为l的水印需要n_v(n_v=n_w/2)个坐标点。将二进制流W^binary转换为十进制数组成的一维数组W^decimal，通过c来分割W^binary，确保每一个坐标值正好嵌入一个水印信息，W^decimal={w_j^decimal|w_j^decimal∈{0, 2^c-1}, j∈[0, n_w-1]}，其中

(1)

在嵌入水印时，将W^decimal中的n_w个水印序列按照设定的顺序嵌入n_w个坐标值中。

在水印提取阶段，设(W^decimal)′是从n_w个坐标值中提取出来的水印信息组成的数据，(W^decimal)′={(w_j^decimal)′|(w_j^decimal)′∈(0, 2^c-1), j∈(0, n′_w-1)}，通过式(2)将十进制流水印信息转换为二进制流(W^binary)′={(w_i^binary)′|(w_i^binary)′∈{0, 1}, i∈[0, l-1]}

(2)

1.2 水印嵌入方法

水印嵌入的基本思路如下：①建立一个可包含所有特征点的状态坐标轴，并将特征点按横坐标的值置于轴内；②设定合适区间等分状态坐标轴，使各个点在状态区间内；③继续将状态区间等分，得到状态值子区间；④利用水印信息依次调制重排序后的特征点的状态值，实现水印嵌入。

对区间的状态值(也称为序列号)作如下定义：假设x^left、x和x^right(x^left≤x≤x^right)是基于x轴排序序列中的3个坐标，将x^left和x^right之间的区间等分为m(m≥2)个子区间，那么x坐标所在位置的子区间序列号就是区间状态值s。图 1中所示x的区间状态值为3。

图 1 m个等分区间时x的区间状态值

Fig. 1 The interval state value of x at m equal intervals

图选项

定义Q为包含x^left、x和x^right所在的区间，即Q=(x^left, x, x^right)；x^center为区间Q的中点，即x^center=(x^left+x^right)/2。设区间个数m=2^c+1，则Q被等分为2^c+1个子区间，如图 2所示。

图 2 x在Q中的区间状态值

Fig. 2 The interval state value of x in Q

图选项

定义w为x需要携带的水印，其中0≤w≤2^c。x在嵌入水印后，其数值将映射到水印w所在的区间，因为0≤w≤2^c，所以x的数值依然会在区间[x^left, x^right)中。区间Q内的任意位置的x通过这种映射方法都不会超过Q的范围。

定义r为x位于区间Q中的系数，当x < x^center时，r取值0；当x≥x^center时，r取值1。r值的设定规则如图 3所示。

图 3 r值的设定规则

Fig. 3 The setting rule of r

图选项

在水印嵌入过程中，需要嵌入水印的x在区间Q的系数r可以通过式(3)进行计算

(3)

将x进行水印嵌入后的值记为x′，x′在Q中的区间状态值记为s′，s′也就是x′应该落在Q中的子区间的序号，其计算公式为

(4)

因为w取值范围为0≤w≤2^c，由式(3)可知，r的取值为0或1，所以s′的取值范围为0≤s′ < 2^c+1，x′也不会超过Q的范围。通过式(5)移动x到第s′个区间

(5)

即获得嵌入水印的x′的值。其中，

，是待嵌入水印的坐标点到系数区间最左端的距离(图 4)。可以将k转换为

(6)

图 4 x经过区间映射嵌入水印

Fig. 4 The watermark embedded in x through interval mapping

图选项

定义区间Q的子区间的等间距l_s为

(7)

因为

，根据式(5)可知，x′将位于s′这个子区间内，如图 4所示。

通过对x区间状态值的映射可以实现水印的嵌入，如果区间Q不变，可以提取出x′的区间状态值，然后通过算法可以恢复原始的x。

在水印的提取阶段，x′在Q的区间状态值s′，可以通过x′与Q的子区间间隔l_s求出

(8)

将式(8)转换为

(9)

则水印x′携带的水印值w可以通过式(10)被提取出来

(10)

最后，原始x可以通过式(11)进行恢复

(11)

在上述算法的实现过程中，坐标的区间状态值与水印的信息位进行了一一映射，坐标点所在的区间位置也记录了对应水印信息的索引信息，在数学原理上不仅能够对坐标进行水印嵌入，而且可以精确地提取出水印并恢复数据。

1.3 基于MPR的扰动控制

以上水印嵌入算法中，因为x只在x^left和x^right之间移动，要复现水印算法只需要构建固定的x^left和x^right。如果还能控制x^left和x^right的差值，那么就能将x控制在一定范围内，从而控制水印嵌入对数据的扰动程度。文献[26]提出的MPR方法已证实坐标点在MPR内移动不会使数据发生拓扑变化，数据的几何特征也能得到较好维持，本文采用MPR方法，对x的区间范围进行限定。MPR的计算如图 5所示：获取原始地图G的所有坐标点，基于这些点构建Delaunay三角网T，用t_i表示T中的三角形，δ(t_i^j)为t_i中内切圆的半径；得到与坐标点v_i邻接的所有三角形内切圆后，其中半径最小的内切圆就是v_i的MPR。以v_i为圆心，R_i为半径的圆就是坐标点v_i的MPR

(12)

图 5 基于Delaunay三角网的最大扰动区域计算方法

Fig. 5 The maximum perturbation region was calculated based on Delaunay triangular mesh

图选项

式中，r_i=min_{1≤j≤d_i}δ(t_i)；t_i^j为与坐标点v_i邻接的所有三角形；d_i为t_i^j的数量。

计算得到所有坐标点的MPR后，取最小半径R_min作为x^left和x^right的差值，此时x^left=x-R_min，x^right=x+R_min。为了在水印提取时能重构状态区间轴，在水印嵌入前，首先找到坐标轴中最值x_max和x_min，然后自定义设置坐标轴大区间L_x的最小值和最大值，要求X_min≤x_min、X_max≥x_max。本文采用图层空间域的最值作为X_min和X_max，可实现全盲检测。提取时，只需要输入X_min和X_max作为密钥重构新的坐标轴即可。定义L_X(X_min, X, X_max)为包含X_min、X和X_max的区间，将区间L_X等分为D_x(D_x≥2)个子区间，每个子区间的长度为2R_min，则有

(13)

对X={x_i|i∈[1, n]}中的每一个待嵌入水印的坐标点x_i计算左右相邻的虚拟坐标值x_i^left和x_i^right，并组成待嵌入水印的区间Q(x_i^left, x_i, x_i^right)。定义λ_{x_i}(0≤λ_{x_i}≤D_x-1)为x_i在L_x中的区间序列号，可以通过式(14)求出

(14)

根据λ_{x_i}可以相继求出x_i的相邻虚拟坐标x_i^left和x_i^right

(15)

图 6为整个X坐标区间及x的虚拟坐标。

图 6 虚拟坐标

Fig. 6 Virtual coordinates

图选项

同样地，也可以在y坐标中嵌入水印，在2R_min基础上得到两个相邻虚拟坐标y_i^bottom和y_i^top，并以同样的方式嵌入水印。

1.4 水印嵌入流程

基于区间映射与MPR扰动控制的矢量地图可逆水印算法嵌入流程如图 7所示，具体步骤如下。

图 7 水印嵌入流程

Fig. 7 Flowchart of watermark embedding

图选项

(1) 数据预处理：对原始矢量地图G，采用Douglas-Peucker算法，设置阈值τ，提取特征点p_i；计算各特征点至地图几何中心的距离d_i，按距离值大小重排序；将二值水印图像I进行Arnold置乱，选择置乱加密的图像I^A读取为二进制流水印信息W^binary，并转换为十进制数W^decimal。

(2) 构建区间状态轴：根据图层域的最值X_min和X_max建立一个可包含所有特征点x坐标的坐标轴L_X，并将特征点集按坐标值大小置于L_X中的相应位置，得到区间状态轴。

(3) 获取区间状态值：根据原始矢量地图计算所有坐标点的MPR，取最小MPR的直径2R_min作为区间状态轴的区间范围Q；利用水印位c继续将区间Q划分为2^c+1个子区间，根据特征点所处的子区间位置获取每个特征点的区间状态值s。

(4) 嵌入水印信息：根据1.2节的方法，用水印信息W^decimal依次调制重排序后的特征点的状态值，使所有特征点得到新的区间状态值s′，并获得含水印信息的特征点p′_i。

(5) 将含水印特征点与非特征点一起写入文件，得到含水印矢量地图G′。

1.5 水印提取与检测

水印提取是水印嵌入的逆过程，基本步骤如下。

(1) 重构状态区间轴：采用Douglas-Peucker算法预设的阈值τ压缩含水印矢量地图G′，获取特征点集合p′_i，并根据各个点至几何中心的距离重现排序；根据图层域的最值X_min和X_max，以及MPR的最小半径R_min重构状态区间轴，将各个特征点置于坐标轴内，得到状态值s′。

(2) 水印提取与数据恢复：根据1.2节提出的水印算法进行水印信息提取并同时恢复坐标值，对提取出的十进制水印数组(W^decimal)′进行二进制流转换，并读取为二进制图像(I^A)′；通过Arnold逆置乱，恢复得到水印图像I′。

(3) 水印检测：将提取到的多个水印图像I′与原始水印图像I进行对比验证，用归一化相关系数(normalized correlation，NC)衡量图像之间的相似性^[27^]，具体公式如下

(16)

式中，w_i, j与w′_i, j分别为原始水印图像与提取水印图像在第(i, j)位的水印信息。计算提取的水印信息具体误差时，用比特误差率(bit error rate，BER)来表示

(17)

式中，L为水印信息的长度。水印容量是衡量被保护数据含有水印信息量大小的主要指标，矢量地图的基本计量为顶点，因此水印容量指标定义为

(18)

在水印提取时，要对密钥进行解析，按照顺序对嵌入的水印点进行空间查找，从而提取出水印信息。一旦水印点被删除，查找不到，采用空值占位，也不会影响后续顶点的顺序。

2 试验与分析

2.1 试验数据处理

2.1.1 原始矢量地图

本文选取涵盖点、线、面3种类型的矢量地图作为试验数据，如图 8所示，分别为交通站台点、道路、水体及建筑物轮廓图。表 1统计了预处理之后各个数据的基本信息，其中线、面数据的压缩阈值设置为2 m，点数据不进行压缩简化。通过计算，得到各个数据的MPR半径，确定了各数据的坐标点活动范围。可以看出，点密度越大的数据，其MPR范围越小，嵌入水印时可扰动的空间也越小。

图 8 原始矢量地图

Fig. 8 Original vector maps

图选项

表 1 原始矢量地图的数据信息

Tab. 1 Data information of original vector maps

表选项

2.1.2 水印信息

试验嵌入的水印信息是一张大小为64×64像素的二值图像(图 9(a))，图像内含有版权信息“地测学院”。为了提升安全性，采用Arnold算法^[28^]对其进行加密，在Arnold公式中所有参数取1的情况下(Arnold算法的公式参考文献[28])，得到本文水印图像的置乱周期为48次。图 9展示了水印图像在置乱10次、20次、30次及48次时的结果，其中第10次、20次和30次置乱的水印图像为难以识别的噪声图，具有较好的隐蔽性，在第48次置乱时恢复成原始状态。本文选择第20次置乱的图像作为加密态的水印信息用于水印嵌入与提取。

图 9 水印图像

Fig. 9 Watermark images

图选项

2.2 水印容量分析

通过对二进制流水印信息的转化，得到c的最大可取数值为38，即每个坐标值最多可携带38 bit的水印码。为了防止水印嵌入后的坐标值在保存时被数据格式精度约减，产生较大误差，一般取较小值。综合考虑shape数据的存储特征及可嵌入量，将c值设置为8，即每个坐标值可携带8 bit信息。在可逆水印算法中，水印容量越大，提取出完整水印信息的概率越大。表 2对比了近几年较热的大容量水印算法以及本文算法下的水印容量。在本文算法中，交通站点数据为非压缩情况下的有效水印嵌入量，可达到16 bit/点；水体数据中仅为0.578 8 bit/点，主要是因为水体数据的点密度很大，压缩后得到的特征点数量很少，因此水印容量明显小于其他数据。通过对比分析，本文算法在水印容量方面几乎仅次于文献[32]，在有效水印嵌入量方面更是达到了16 bit/点，水印容量性能突出，本文算法的大容量特点主要得益于水印信息的十进制化操作以及坐标点状态值位置调制的包容性。

表 2 水印容量对比

Tab. 2 Comparison of watermark capacity

表选项

2.3 不可见性分析

从拓扑质量、视觉观测和误差分析这3个方面进行水印不可见性分析。对含水印矢量地图拓扑检查发现不存在要素重叠及容差错误，因此水印嵌入未导致拓扑变更，图形质量得到较好维持。图 10展示了道路图在水印嵌入前、后的对比情况，在不放大的情况几乎难以察觉坐标点发生了偏移，放大后可见坐标点发生了细微改变，但要素的形状不存在明显的形变，要素的几何特征维持良好。表 3统计了水印嵌入后与原始数据的误差值，在所有数据中，最大的误差改变量为0.000 851 32 m。根据国标《1∶500 1∶1000 1∶2000外业数字测图技术规程GB14912—2005》中的地形图精度平面误差要求，1∶500地形图的点位误差为±0.15 m、1∶1000地形图的点位误差为±0.30 m、1∶2000地形图的点位误差为±0.60 m，由此可见，水印嵌入对数据的扰动非常小，不会影响数据的可用性。

图 10 矢量地图嵌入水印前、后的叠置对比

Fig. 10 Superposition comparison of watermarked vector map and original vector map

图选项

表 3 水印嵌入后的精度误差统计

Tab. 3 Accuracy error statistics after watermark embedding

表选项

2.4 稳健性分析

为了检验本文算法的稳健性，分别对含水印矢量地图进行几何攻击、坐标点攻击以及裁剪攻击等，并与相关的水印算法作对比分析。限于篇幅，以下试验均在道路数据上进行。

2.4.1 几何攻击

本次试验的几何攻击主要包括旋转、平移、缩放及简化攻击。表 4的统计结果显示，本文算法以及对比算法在不同强度的旋转、平移及缩放下，都可以提取到完整的水印图像，具有良好的抗几何攻击能力。在这些攻击下，数据未发生扭曲，得到的最小MPR的半径与原始数据的R_min是一致的，特征点与原始数据也一致，因此可获取完整水印信息。本文算法预设的压缩阈值为2 m，在简化攻击时，不超过2 m的压缩阈值时可得到完整水印图像，超过2 m时，随着特征点的减少，可获取的水印信息也不断减少；而对比算法中只要有简化攻击，就无法提取出完整水印图像，且随着压缩阈值的提升，甚至无法完成提取试验。

表 4 几何攻击试验结果

Tab. 4 Experiment results of geometric attacks

表选项

2.4.2 坐标点攻击

坐标点攻击包括随机增点和随机删点，试验强度为原数据的10%~60%逐步递增。图 11所示为本文算法与文献[21]算法和文献[22]算法在不同强度坐标点攻击下得到的水印图像BER的变化曲线。图 11(a)中，随机增点攻击对本文算法和文献[21]算法没有任何影响，始终可以得到BER值为0的完整水印图像；图 11(b)中，随着删点比例的逐步增加，3种算法下可得到的水印图像的BER值都越来越高，本文算法与文献[22]算法结果比较接近，即便在60%的删点比例下，仍可以得到BER值低于10%的水印图像，展现出良好的抗删点攻击能力。综上，本文算法在抗增点和删点攻击均展现出优越的性能，具有较强的抗坐标点攻击的能力。

图 11 坐标点攻击下水印图像的BER变化曲线

Fig. 11 BER variation curves of the extracted watermark images under point attacks

图选项

2.4.3 裁剪攻击

表 5统计了在不同裁剪比例及裁剪位置情况下，本文算法及对比算法可提取到的最优水印图像。本文算法中，含水印特征点被删除时，该点所带的水印信息自动以空值占位，因此提取的水印图像不存在噪点。3种算法都随着裁剪尺度的增加，可提取的水印图像质量越来越差，但都展现出较强的稳健性，在30%、40%及50%裁剪比例下均可以获取可识别的水印图像。本文算法因有空值占位，在同一裁剪比例下，得到的NC值与BER值要明显要好于对比算法的结果；但在60%的裁剪尺度下，得到的水印可识别度较差，难以有效验证版权，而文献[21—22]得到的结果虽也难以完成验证功能，但仍有一定辨识度。总体而言，本文算法对裁剪攻击具有良好的稳健性，在50%的裁剪尺度下仍可提取出可靠的水印图像，足以胜任常规使用。

表 5 裁剪攻击试验结果

Tab. 5 Experiment results of cropping attacks

表选项

2.5 可逆性分析

从理论以及数学的角度，本文算法是可以完美无损地恢复原始矢量地图的，但对恢复后的数据与原始数据对比发现，恢复后的数据与原始数据存在一些极小的误差。主要原因在于受限于部分数据的MPR半径过小及GIS软件的精度位存储长度的限制，可供运算的结果在存储时会在精度位上进行四舍五入的约简，从而导致恢复后的数据与初始数据有一些微小误差。误差结果见表 6，可以看出，这些误差对数据质量的影响很小，且仅限于数量不多的特征点，因此，本文算法具有较好的可逆性。

表 6 数据恢复后与原始数据的误差统计

Tab. 6 Error statistics between recovered data and original data

表选项

3 结论

本文基于量化索引调制的思想设计了一种矢量地图区间映射的可逆水印算法，利用水印信息位调制坐标点的区间状态值的方式嵌入水印信息，通过区间位与水印信息位的一一映射实现水印信息的可逆提取与盲检测。在水印嵌入阶段，将二值水印信息进行十进制转化，并用于坐标区间状态值的调制，利用十进制数相较二进制数更强的信息承载能力实现水印的大容量嵌入；坐标的区间状态值范围采用最大扰动区域方法进行限制，将坐标扰动控制在合理范围内，从而控制了水印的隐蔽性。结合特征点提取，基于空间特征重排序等数据预处理，本文提出的可逆水印算法对常规几何攻击、坐标点攻击以及裁剪攻击等具有较强的稳健性，能较好地平衡水印方案的容量、不可见性及稳健性之间的相互制约关系。

作者简介

第一作者简介：奚旭(1990—)，男，博士，讲师，研究方向为空间数据安全防护。E-mail：xixu@usts.edu.cn

通信作者：张新长, E-mail：zhangxc@gzhu.edu.cn

初审：张艳玲

复审：宋启凡

终审：金君