抗投影变换的矢量地理数据水印算法
作者简介:佟德宇(1989-),男,江苏高邮人,博士生,研究方向为地理数据安全。E-mail:tdyforweb@163.com
收稿日期: 2015-07-27
要求修回日期: 2015-09-11
网络出版日期: 2016-08-10
基金资助
国家社科基金重大项目(11&ZD162)
国家自然科学基金项目(41301413)
江苏省自然科学基金项目(BK20130903)
公益性行业科研专项(201512021、201512019)
A Watermarking Algorithm Resisting to Projection Transformation for Vector Geographic Data
Received date: 2015-07-27
Request revised date: 2015-09-11
Online published: 2016-08-10
Copyright
在矢量地理数据水印算法研究中,以往研究较多考虑增删点、裁剪、平移、旋转等攻击方式,而对抗投影变换的攻击方式研究较少。由于投影变换在GIS中具有重要意义,故本文提出了一种抗投影变换的矢量地理数据水印算法。水印嵌入前,对待嵌矢量地理数据利用道格拉斯算法进行压缩,使用四叉树分块选取特征点,提取并保存特征点及其属性信息;水印嵌入采用坐标映射和量化机制以增强水印算法的鲁棒性;检测水印时,将待检测数据与原始特征点进行属性信息匹配,匹配成功的同名点采用二元三次多项式进行最小二乘法拟合,根据拟合的多项式系数对待检测数据进行投影变换,最终实现水印信息的提取。实验结果表明,本方法能抵抗投影变换攻击、增删点、几何变换以及它们的复合攻击,具有较好的可行性和实用性。
佟德宇 , 任娜 , 朱长青 , 林威 . 抗投影变换的矢量地理数据水印算法[J]. 地球信息科学学报, 2016 , 18(8) : 1037 -1042 . DOI: 10.3724/SP.J.1047.2016.01037
Among the researches on watermarking algorithms for vector geographic data, more attentions have been paid on the attacking methods such as increasing points, deleting points, translation, rotation and so on. However, there is scarce research on watermarking algorithms resisting the attack of projection transformation. As a common data processing method in geography, projection transformation would affect and ruin the watermark embedded by normal ways in vector geographic data. Therefore this paper proposes an anti-projection watermarking algorithm for this kind of data based on the correspondence points matching. First, feature points of the data are extracted using Douglas-Peucker compression and quad-tree partition. And the points are stored with their attribute information together. Then the watermark is embedded based on coordinate mapping and quantization mechanism, which would enhance the robustness of the watermarking algorithm. Data containing watermark can now be distributed and shared. When it is needed to detect the watermark from the suspect data, the detected one will be compared and matched with the points stored according to their attribute information. If the corresponding points are matched, the project transformation of them can be realized by the bivariate cubic polynomial. When calculating the coefficients of the polynomial, the least square method and QR decomposition are used to improve its accuracy. After the calculation, the polynomial will help to transform the projection coordinate system from the detected data to the original data. And the watermark will be able to be detected after the transformation if enough accuracy is ensured. Experiments have been conducted to prove the proposed algorithm is robust against the attacks of projection transformation, adding data, deleting data, ordinary geometric transform and their compound attacks. Results of the experiments show that the algorithm has a good feasibility and can be taken into practice easily.
Fig.1 General view of the experimental data图1 实验数据示意图 |
Tab.1 Result of the algorithm robustness after the project transformation表1 投影变换攻击鲁棒性实验结果 |
实验编号 | 实验1 | 实验2 | 实验3 | 实验4 | 实验5 | 实验6 |
---|---|---|---|---|---|---|
实验数据 | 第1组 | 第1组 | 第2组 | 第2组 | 第3组 | 第3组 |
原始坐标系 | WGS84 | 墨卡托 | 兰伯特 | WGS84 | 高斯克吕格 | 兰伯特 |
投影攻击后坐标系 | 墨卡托 | WGS84 | WGS84 | 兰伯特 | 兰伯特 | 高斯克吕格 |
水印精度/m | 0.1 | 0.1 | 0.01 | 0.01 | 0.001 | 0.001 |
X坐标误差最大值 | 3.6 | 4.6 | 5.6×10-3 | 4.9×10-3 | 4.3×10-3 | 3.3×10-3 |
X坐标误差最小值 | 4.3×10-8 | 9.3×10-8 | 6.0×10-8 | 2.5×10-8 | 5.7×10-9 | 6.1×10-9 |
X坐标误差平均值 | 0.12 | 0.14 | 1.5×10-3 | 2.4×10-3 | 7.9×10-4 | 4.7×10-4 |
Y坐标误差最大值 | 2.5 | 3.2 | 0.32 | 0.27 | 5.6×10-3 | 4.3×10-3 |
Y坐标误差最小值 | 6.6×10-8 | 3.8×10-8 | 1.3×10-8 | 8.0×10-8 | 3.5×10-8 | 1.9×10-8 |
Y坐标误差平均值 | 0.14 | 0.15 | 2.6×10-3 | 1.9×10-3 | 5.1×10-4 | 3.7×10-4 |
自相关系数 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 |
检测结果 | 成功 | 成功 | 成功 | 成功 | 成功 | 成功 |
运行时间/s | 6.78 | 6.96 | 3.28 | 3.04 | 11.90 | 10.29 |
Tab.2 Result of the algorithm robustness after the compound attack表2 复合攻击鲁棒性实验结果 |
实验7 | 实验8 | 实验9 | 实验10 | 实验11 | 实验12 | 实验13 | |
---|---|---|---|---|---|---|---|
实验数据 | 第2组 | 第1组 | 第1组 | 第2组 | 第3组 | 第2组 | 第3组 |
复合攻击方式 | 压缩20%+实验3 | 增加30%+实验1 | 删除30%+实验2 | 平移500 m+实验4 | 旋转120°+实验5 | 压缩20%+平移500 m+实验4 | 增加30%+旋转60°+实验6 |
自相关系数 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 |
检测结果 | 成功 | 成功 | 成功 | 成功 | 成功 | 成功 | 成功 |
运行时间/s | 2.98 | 7.23 | 6.22 | 3.03 | 11.56 | 2.86 | 13.45 |
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
[
|
[4] |
|
[5] |
[
|
[6] |
|
[7] |
[
|
[8] |
[
|
[9] |
[
|
[10] |
|
[11] |
|
[12] |
[
|
[13] |
|
[14] |
|
[15] |
[
|
[16] |
[
|
[17] |
[
|
[18] |
[
|
/
〈 |
|
〉 |