CN104751399A - 一种基于信息分存的涉密线状地理要素伪装与还原方法 - Google Patents

Abstract

本发明公开了一种基于信息分存的涉密线状地理要素伪装与还原方法，包括如下过程：(1)针对线类型GIS矢量数据伪装过程，包括线要素分级、分段、平移旋转处理、生成新图层等步骤；(2)针对线类型GIS矢量数据还原过程，包括线要素分段、旋转平移处理、线要素重组还原等过程。本发明的方法具有较高的隐蔽性和不受容量限制的优势，为矢量地理数据的信息伪装提供了新思路，在隐藏秘密信息存在性的基础上有效保障秘密数据在数据传输、隐藏通信中的安全性。

Description

一种基于信息分存的涉密线状地理要素伪装与还原方法

技术领域

本发明属于地理信息安全领域，具体涉及一种基于信息分存的涉密线状地理要素伪装与还原的方法。

背景技术

地理信息的安全防护涉及国家安全和社会利益，是当前急需解决的国家和社会重大需求问题。在对秘密GIS矢量数据信息内容隐藏的基础上，进一步隐藏秘密信息的存在，结合分存技术从形态特征进行伪装。

目前，信息分存技术主要应用于图像领域，将图像信息分为具有一定可视效果的n幅图像，如果知道图像信息中的m(m<n)幅子图像，则该图像可以得到恢复，如果图像少于m幅，则图像无法得到恢复。借鉴分存技术在图像领域的应用，将分存技术应用到矢量地理数据信息伪装中来，伪装后数据同样具有可视效果，并且对伪装后的数据进行分图层存储，以增加算法的安全性。

发明内容

本发明的目的在于：基于信息分存伪装技术，提出一种针对GIS线图层数据中线状地理要素的伪装与还原方法，能有效保障涉密线数据在数据传输、隐藏通信中的安全性。

为了实现上述目的，本发明采取的技术方案主要包括：

一种基于信息分存的涉密线状地理要素伪装与还原方法，具体过程如下：

(一)针对线类型GIS矢量数据伪装过程

步骤11：加载一幅矢量地图数据R，输入初始平移量c作为用户密钥参数，c∈[0,100]，存储到密钥文件key中；

步骤12：读取数据R中所有线要素L＝{l_i|i＝1,2,…,m}，并为每个线要素中的点按顺序设置序号，m为线要素的个数；

步骤13：利用道格拉斯-普克算法，将线要素l_i分为两级l_i1和l_i2，可根据安全性需要对一级线要素l_i1进行适当的伪装处理；

步骤14：对线要素l_i2进行分段处理，

(a)线要素l_i1中的点集合记为P_i1＝{p_j|j＝1,2,…,n},n为线要素l_i1中的点个数；

(b)依次读取点集合P_i1中相邻两个点p_j、p_j+1的序号I_j、I_j+1，在线要素l_i2的点集合中遍历序号位于区间(I_j,I_j+1)之间的所有点，将这些点分为一段；

(c)循环步骤(b)，直至将线要素l_i2分为(n-1)个分段，记为S＝{s_k|k＝1,2,…,n-1}；

步骤15：确定当前线要素l_i的基线l的直线方程，

(a)假设点集合P_i1的首尾点坐标分别记为p_b(x₁,y₁),p_e(x_n,y_n)，根据公式(1)计算基线斜率k：

$k=\left\{\begin{array}{c}\&infin;,(x_n=x_1)\\ \frac{y_n-y_1}{x_n-x_1},(x_n\&NotEqual;x_1)\end{array}\right.---\left(1\right)$

(b)假设p_b点经平移后的点p_b'在基线l上，p_b'坐标记为(x₁',y₁')，平移公式如(2)所示，其中c为用户指定参数，

$\left\{\begin{array}{c}{x_1}^{\&prime;}=x_1+c\\ {y_1}^{\&prime;}=y_1+c\end{array}\right.---\left(2\right)$

(c)基线l的方程可表示为公式(3)：

$\left\{\begin{array}{c}x={x_1}^{\&prime;},(x_n=x_1)\\ y=k(x-{x_1}^{\&prime;})+{y_1}^{\&prime;},(x_n\&NotEqual;x_1)\end{array}\right.---\left(3\right)$

步骤16：对点集合P_i1做平移旋转变换使其落在基线l上,处理后的点作为线要素l_i2中每段点集合的旋转基点；平移后的点集合记为P_i1'＝{p_j'|j＝1,2,…,n}，旋转后的点集合记为P_i1″＝{p_j″|j＝1,2,…,n}，

(a)根据公式(4)、(5)计算当前分段s_k对应的横坐标平移偏移量h_j和纵坐标平移偏移量v_j，

$h_j=\left\{\begin{array}{c}c,(j=1)\\ {x_j}^{\&prime;\&prime;}-x_j(j>1)\end{array}\right.---\left(4\right)$

$v_j=\left\{\begin{array}{c}c,(j=1)\\ {y_j}^{\&prime;\&prime;}-y_j(j>1)\end{array}\right.---\left(5\right)$

其中H＝{h_j|j＝1,2,…,n}为横坐标平移偏量集合，V＝{v_j|j＝1,2,…,n}为纵坐标平移偏量集合，(x_j″,y_j″)为p_j″中的点坐标；

(b)根据公式(6)计算平移后的点集合P_i1'中的点坐标p_j'(x_j',y_j')，

$\left\{\begin{array}{c}{x_j}^{\&prime;}=x_j+h_j\\ {y_j}^{\&prime;}=y_j+v_j\end{array}\right.---\left(6\right)$

(c)根据公式(1)计算当前段直线p_jp_j+1斜率k_j，根据公式(7)确定基线l和当前段直线p_jp_j+1之间的夹角α_j；当前段直线p_jp_j+1按逆时针旋转与基线重合的旋转角度α_j'如公式(8)所示，旋转角度集合记为A＝{α_j'|j＝1,2,…,n}，

${\mathrm{\&alpha;}}_j=\left\{\begin{array}{c}0,(k_j\&RightArrow;\&infin;,k\&RightArrow;\&infin;)\\ \mathrm{\&pi;}/2,({\mathrm{kk}}_j=-1)\\ \arctan k-\mathrm{\&pi;}/2,(k_j\&RightArrow;\&infin;)\\ \mathrm{\&pi;}/2-{\arctan k}_j,(k\&RightArrow;\&infin;)\\ \arctan \left(\frac{k-k_j}{1+k*k_j}\right),({\mathrm{kk}}_j\&NotEqual;-1)\end{array}\right.---\left(7\right)$

(d)根据公式(9)、(10)计算旋转后的点集合p_j″中的点坐标(x_j″,y_j″)，

${x_j}^{\&prime;\&prime;}=\left\{\begin{array}{c}{x_j}^{\&prime;},(j=1)\\ ({x_j}^{\&prime;}-{x_{j-1}}^{\&prime;\&prime;})*\cos {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}-({y_j}^{\&prime;}-{y_{j-1}}^{\&prime;\&prime;})*\sin {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}+{x_{j-1}}^{\&prime;\&prime;},(j>1)\end{array}\right.---\left(9\right)$

${y_j}^{\&prime;\&prime;}=\left\{\begin{array}{c}{y_j}^{\&prime;},(j=1)\\ ({y_j}^{\&prime;}-{y_{j-1}}^{\&prime;\&prime;})*\cos {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}-({x_j}^{\&prime;}-{x_{i-1}}^{\&prime;\&prime;})*\sin {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}+{y_{j-1}}^{\&prime;\&prime;},(j>1)\end{array}\right.---\left(10\right)$

(e)循环处理(a)-(d)，直至段S中的各段处理完毕；

步骤17：对二级线要素l_i2中所有段S＝{s_k|k＝1,2,…,n-1}进行伪装处理；假设段s_k的点集合记为SP＝{sp_w|w＝1,2,…,t},其中t为该段中点的个数，

(a)对段s_k进行平移处理，段s_k中的点集合SP平移后的集合记为SP'＝{sp_w'|w＝1,2,…,t},横纵坐标对应的平移偏量分别为h_w和v_w，根据公式(6)做相应的平移变换；

(b)对平移后的每段进行旋转处理，SP'旋转后的集合记为SP″＝{sp_w″|w＝1,2,…,t}，每段对应的旋转角度为α_W',旋转基点为点集合SP″中对应的点sp_w″,根据公式(9)、(10)做相应的旋转变换；

(c)循环执行(a)和(b)，直到所有段完成伪装处理；

步骤18：生成新图层数据R1和R2，将线要素l_i1中的点集合按原始顺序组成的线要素添加到数据R1中，将线要素l_i2中经伪装变换后的点集合按照原始顺序组成的线要素存放到数据R2中；

步骤19：重复步骤12-18，直到所有线要素完成伪装处理为止；

(二)针对线类型GIS矢量数据还原过程

步骤21：打开伪装后的矢量数据R1和R2，打开密钥文件key，读取用户密钥参数c；

步骤22：分别读取数据R1和数据R2中的线要素L_i1＝{l_i1|i＝1,2,…,m}，L_i2＝{l_i2|i＝1,2,…,m}，其中，m为数据R1中线要素的个数；

步骤23：对线要素l_i2进行分段处理，

(a)线要素L_i1中还原后的点集合记为P_i1＝{p_j|j＝1,2,…,n},n为线要素l_i1中的点个数；

(b)依次读取点集合P_i1中相邻两个点p_j、p_j+1的序号I_j、I_j+1，在线要素l_i2点集合中遍历序号位于区间(I_j,I_j+1)之间的所有点，将这些点分为一段；

(c)循环步骤(b)，直至将线要素l_i2分为(n-1)个分段，记为S＝{s_k|k＝1,2,…,n-1}；

步骤24：确定当前线要素l_i的基线l的直线方程，

(a)假设点集合P_i1的首尾点坐标分别记为p_b(x₁,y₁),p_e(x_n,y_n)，根据公式(1)计算基线斜率k；

(b)假设p_b点经平移后的点p_b'在基线l上，p_b'坐标记为(x₁',y₁')，平移公式如(2)所示，其中c为用户指定参数；

(c)基线l的方程可表示为公式(3)；

步骤25：对点集合P_i1做平移旋转变换使其落在基线l上,处理后的点作为每段点集合的旋转基点，

(a)根据公式(4)、(5)计算当前分段s_k对应的横坐标平移偏移量h_j和纵坐标平移偏移量v_j；

(b)根据公式(6)计算平移后的点集合P_i1'中的点坐标p_j'(x_j',y_j')；

(c)根据公式(1)计算当前段直线p_jp_j+1斜率k_j，根据公式(7)确定基线l和当前段直线P_jP_j+1之间的夹角α_j，当前段直线p_jp_j+1按逆时针旋转与基线重合的旋转角度α_j'由公式(8)计算得出，旋转角度集合记为A＝{α_j'|j＝1,2,…,n}；

(d)根据公式(9)、(10)计算旋转后的点集合_pj″中的点坐标(x_j″,y_j″)；

(e)循环处理(a)-(d)，直至段S中各段处理完毕；

步骤26：对二级线要素l_i2中所有段S＝{s_k|k＝1,2,…,n-1}进行还原处理；假设段s_k的点集合记为SP＝{sp_w|w＝1,2,…,t},sp_w坐标记为(x_w,y_w),其中t为该段中点的个数，

(a)根据公式(9)、(10)对每段进行旋转处理：点集合SP中各点旋转后的点集合记为SP'＝{sp_w'|w＝1,2,…,t},其中每段对应的旋转角度为A'＝{α_j″|j＝1,2,…,n}，由公式(11)计算得出，

α_j″＝2π-α_j'   (11)

(b)对旋转后的每段进行平移操作：点集合SP'平移后的集合记为SP″＝{sp_w″|w＝1,2,…,t},横纵坐标对应的平移偏量分别为-h_j和-v_j，根据公式(9)做相应的平移变换；

(c)循环执行(a)和(b)，直到完成所有段的还原处理；

步骤27：将线要素l_i2中还原后的点集合与线要素l_i1中的点集合按照坐标序号从小到大的顺序重组成为线要素l_i，完成当前线要素的还原处理；

步骤28：重复步骤22-27，直到完成每个线要素的还原处理。

本发明结合分存技术从形态特征进行伪装，使伪装后的矢量数据仍然是有意义的数据，从而降低攻击者的注意力，提高伪装的隐蔽性，并且避免了传统伪装中容量的问题，具有较高的隐蔽性和不受容量限制的优势，为矢量地理数据信息伪装提供了新思路。本发明在隐藏秘密信息存在性的基础上，能有效保障秘密数据在数据传输、隐藏通信中的安全性，在一定程度上满足了GIS矢量数据的隐藏通信和安全传输需求。

附图说明

图1是本发明方法的流程图，(a)为伪装流程图，(b)为还原流程图；

图2是本发明原始数据，(a)为原始数据，(b)一级原始数据,(c)为二级原始数据,(d)为二级伪装数据；

图3是本发明原始二级数据和伪装数据对比图，(a)为二级伪装与原始数据对比图，(b)为局部细节对比图，其中虚线表示原始二级数据，实线表示伪装二级数据。

具体实施方式

本实施例选择shp格式矢量数据，针对数据的读取、伪装、还原的整个过程，进一步详细说明本发明。本实例选择中国河流图。

(一)伪装过程

步骤11：加载中国河流图，输入初始平移量c＝0,存储到密钥文件key中。

步骤12：读取当前线要素l₁，线要素的个数m＝31。

步骤13：利用道格拉斯-普克算法，将线要素l₁分为两级l₁₁和l₁₂。

步骤14：对l₁₂进行分段处理。

(a)l₁₁中的点集合记为P₁₁＝{(435481.944441,4164230.4092540005)，(435481.21390171145,4164252.6066547194)，(435466.944441,4164275.4092540005)},n为3。

(b)依次读取P₁₁中相邻两个点p₁、p₂、p₃的序号0、7、30，在l₁₂点集合中遍历序号位于(0,7)之间的所有点，将序号为1至6的点分为一段。

(c)循环步骤(b)，直至将l₁₂分为2个分段，记为S＝{s_k|k＝1,2}。

步骤15：确定当前线要素l_i的基线l的直线方程。

(a)根据公式(1)计算基线斜率k＝-3.0。

(b)首点经平移后的点p_b'在基线l上，根据公式(2)计算得到p_b'坐标(435481.944441,4164230.4092540005)。

(c)基线l的方程可表示为y＝-3(x-435481.944441)+4164230.4092540005

步骤16：对点集合P₁₁做平移旋转变换使其落在基线l上,处理后的点作为每段点集合的旋转基点。

(a)根据公式(4)、(5)计算当前分段s₁对应的横纵坐标平移偏移量分别为0,0。

(b)根据公式(6)计算平移后的点坐标p_j'(435481.944441,4164230.4092540005)。

(c)根据公式(1)计算当前段直线斜率k₁＝-30.384951318763715，根据公式(7)确定基线l和当前段直线的夹角α_j＝0.28885140051280989，根据当前段直线按顺时针旋转与基线重合的旋转角度α_j'＝6.572036707692396。

(d)根据公式(9)、(10)计算旋转后的P₁″点坐标(435481.944441,4164230.4092540005)

(e)循环处理(a)-(d)，其中H₁＝{0，-5.2926956237060949},V₁＝{0,-0.12769598234444857}，旋转角度集合A＝{6.572036707692396，-0.23741104309531336}旋转基点集合P₁₁″＝{(435481.944441,4164230.4092540005)，(435475.92120608775,4164252.4789587371)}。

步骤17：对二级点l₁₂中所有段进行伪装处理。s₁的点集合记为SP＝{(435482.05815706245,4164243.7518112813)，(435482.03656739253,416426.8110503652)，…,(435481.50921598467,4164251.7062506238)},其中t＝6。

(a)对每段进行平移处理。根据公式(6)，h₁＝0，v₁＝0，s₁中的点集合SP平移后的集合记为SP'＝{(435482.05815706245,4164243.7518112813)，(435482.03656739253,416426.8110503652)，…,(435481.50921598467,4164251.7062506238)}

(b)对平移后的每段进行旋转处理。根据公式(9)、(10)，SP'旋转后的集合记为SP″＝{(435479.2527998961,4164244.2314438336)，(435478.36067612254,4164247.1577939587)，…,(435476.46076735476,4164251.6999775474)}。

(c)循环执行(a)(b)，直到所有段完成伪装处理。

步骤18：生成新图层R1、R2，将l₁中的点集合按原始顺序组成的线要素添加到R1中，将L2中经伪装变换后的点集合按照原始顺序组成的线要素存放到R2中。

步骤19：重复步骤12-18，直到所有线要素完成伪装处理。

(二)还原过程

步骤21：打开伪装后的矢量数据R1和R2，打开密钥文件key，读取用户密钥参数c＝0。

步骤22：分别读取R1和R2中的当前线要素l₁₁，l₁₂＝{(435482.05815706245,4164243.7518112813)，(435482.03656739253,416426.8110503652)，…,(435481.50921598467,4164251.7062506238)，…},其中点个数为28。

步骤23：对线要素l₁₂进行分段处理。

(a)L_i1中的点集合记为P₁₁＝{(435481.944441,4164230.4092540005)，(435481.21390171145,4164252.6066547194)，(435466.944441,4164275.4092540005)},n＝3。

(b)依次读取P₁₁中相邻点的序号0、7、30，在l_i2点集合中遍历序号位于区间(0,7)、(7,30)之间的所有点，将这些点分为一段。

(c)循环步骤(b)，直至将l₁₂分为2个分段，记为S＝{s_k|k＝1,2}。

步骤24：确定当前线要素l_i的基线l的直线方程。

(a)l₁₁中的点集合记为P₁₁＝{(435481.944441,4164230.4092540005)，(435481.21390171145,4164252.6066547194)，(435466.944441,4164275.4092540005)},n为3。

(b)依次读取P₁₁中相邻两个点p₁、p₂、p₃的序号0、7、30，在l₁₂点集合中遍历序号位于(0,7)之间的所有点，将序号为1至6的点分为一段。

(c)循环步骤(b)，直至将l₁₂分为2个分段，记为S＝{s_k|k＝1,2}。

步骤25：对点集合P_i1做平移旋转变换使其落在基线l上,处理后的点作为每段点集合的旋转基点。

(a)根据公式(4)、(5)计算当前分段s_k对应的横纵坐标平移偏移量H₁＝{0，-5.2926956237060949},V₁＝{0,-0.12769598234444857}。

(b)根据公式(6)计算平移后的点集合P_i1'＝{(435481.944441,4164230.4092540005)，(435475.92120608775,4164252.4789587371)。}

(c)由公式(8)计算得出旋转角度集合为A＝{6.572036707692396，-0.23741104309531336}}。

(d)根据公式(9)、(10)计算旋转后的点集合p_j″＝{(435481.944441,4164230.4092540005)，(435475.92120608775,4164252.4789587371)。}

(e)循环处理(a)-(d)，直至S中各段处理完毕。

步骤26：对二级点l_i2中所有段S＝{s_k|k＝1,2,…,n-1}进行还原处理。

(a)根据公式(9)、(10)对每段进行旋转处理：由公式(11)计算得出每段对应的旋转角度，点集合SP中各点旋转后的点集合记为SP'＝{(435482.05815706245,4164243.7518112813)，(435482.03656739253,416426.8110503652)，…,(435481.50921598467,4164251.7062506238)}。

(b)对旋转后的每段进行平移操作：由于对应的横纵坐标偏移量分别为0、0，根据公式(9)，点集合SP'平移后的集合记为SP″＝{(435482.05815706245,4164243.7518112813)，(435482.03656739253,416426.8110503652)，…,(435481.50921598467,4164251.7062506238)}。

(c)循环执行(a)(b)，直到完成所有段的还原处理。

步骤27：将l_i2中还原后的点集合与l_i1中的点集合按照坐标序号从小到大的顺序重组成为线要素l₁，完成当前线要素的还原处理。

步骤28：重复步骤22-27，直到完成每个线要素的还原处理。

Claims (2)

1.一种基于信息分存的涉密线状地理要素伪装与还原方法，其特征在于，具体过程如下：

(一)针对线类型GIS矢量数据伪装过程

步骤11：加载一幅矢量地图数据R，输入初始平移量c作为用户密钥参数，c∈[0,100]，存储到密钥文件key中；

步骤12：读取数据R中所有线要素L＝{l_i|i＝1,2,…,m}，并为每个线要素中的点按顺序设置序号，m为线要素的个数；

步骤13：利用道格拉斯-普克算法，将线要素l_i分为两级l_i1和l_i2；

步骤14：对线要素l_i2进行分段处理，

(a)线要素l_i1中的点集合记为P_i1＝{p_j|j＝1,2,…,n},n为线要素l_i1中的点个数；

(b)依次读取点集合P_i1中相邻两个点p_j、p_j+1的序号I_j、I_j+1，在线要素l_i2的点集合中遍历序号位于区间(I_j,I_j+1)之间的所有点，将这些点分为一段；

(c)循环步骤(b)，直至将线要素l_i2分为(n-1)个分段，记为S＝{s_k|k＝1,2,…,n-1}；

步骤15：确定当前线要素l_i的基线l的直线方程，

(a)假设点集合P_i1的首尾点坐标分别记为p_b(x₁,y₁),p_e(x_n,y_n)，根据公式(1)计算基线斜率k：

$k=\left\{\begin{array}{c}\&infin;,(x_n=x_1)\\ \frac{y_n-y_1}{x_n-x_1},(x_n\&NotEqual;x_1)\end{array}\right.---\left(1\right)$

(b)假设p_b点经平移后的点p_b'在基线l上，p_b'坐标记为(x₁',y₁')，平移公式如(2)所示，其中c为用户指定参数，

$\left\{\begin{array}{c}{x_1}^{\&prime;}=x_1+c\\ {y_1}^{\&prime;}=y_1+c\end{array}\right.---\left(2\right)$

(c)基线l的方程可表示为公式(3)：

$\left\{\begin{array}{c}x={x_1}^{\&prime;},(x_n=x_1)\\ y=k(x-{x_1}^{\&prime;})+{y_1}^{\&prime;},(x_n\&NotEqual;x_1)\end{array}\right.---\left(3\right)$

步骤16：对点集合P_i1做平移旋转变换使其落在基线l上,处理后的点作为线要素l_i2中每段点集合的旋转基点；平移后的点集合记为P_i1'＝{p_j'|j＝1,2,…,n}，旋转后的点集合记为P_i1″＝{p_j″|j＝1,2,…,n}，

(a)根据公式(4)、(5)计算当前分段s_k对应的横坐标平移偏移量h_j和纵坐标平移偏移量v_j，

$h_j=\left\{\begin{array}{c}c,(j=1)\\ {x_j}^{\&prime;\&prime;}-x_j(j>1)\end{array}\right.---\left(4\right)$

$v_j=\left\{\begin{array}{c}c,(j=1)\\ {y_j}^{\&prime;\&prime;}-y_j(j>1)\end{array}\right.---\left(5\right)$

其中H＝{h_j|j＝1,2,…,n}为横坐标平移偏量集合，V＝{v_j|j＝1,2,…,n}为纵坐标平移偏量集合，(x_j″,y_j″)为p_j″中的点坐标；

(b)根据公式(6)计算平移后的点集合P_i1'中的点坐标p_j'(x_j',y_j')，

$\left\{\begin{array}{c}{x_j}^{\&prime;}=x_j+h_j\\ {y_j}^{\&prime;}=y_j+v_j\end{array}\right.---\left(6\right)$

(c)根据公式(1)计算当前段直线p_jp_j+1斜率k_j，根据公式(7)确定基线l和当前段直线p_jp_j+1之间的夹角α_j；当前段直线p_jp_j+1按逆时针旋转与基线重合的旋转角度α_j'如公式(8)所示，旋转角度集合记为A＝{α_j'|j＝1,2,…,n}，

${\mathrm{\&alpha;}}_j=\left\{\begin{array}{c}0,(k_j\&RightArrow;\&infin;,k\&RightArrow;\&infin;)\\ \mathrm{\&pi;}/2,(k{k}_j=-1)\\ \arctan k-\mathrm{\&pi;}/2,(k_j\&RightArrow;\&infin;)\\ \mathrm{\&pi;}/2-\arctan k_j,(k\&RightArrow;\&infin;)\\ \mathrm{arc}\tan \left(\frac{k-k_j}{1+k*k_j}\right),(k{k}_j\&NotEqual;-1)\end{array}\right.---\left(7\right)$

(d)根据公式(9)、(10)计算旋转后的点集合p_j″中的点坐标(x_j″,y_j″)，

${x_j}^{\&prime;\&prime;}=\left\{\begin{array}{c}{x_j}^{\&prime;},(j=1)\\ ({x_j}^{\&prime;}-{x_{j-1}}^{\&prime;\&prime;})*\cos {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}-({y_j}^{\&prime;}-{y_{j-1}}^{\&prime;\&prime;})*\sin {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}+{x_{j-1}}^{\&prime;\&prime;},(j>1)\end{array}\right.---\left(9\right)$

${y_j}^{\&prime;\&prime;}=\left\{\begin{array}{c}{y_j}^{\&prime;},(j=1)\\ ({y_j}^{\&prime;}-{y_{j-1}}^{\&prime;\&prime;})*\cos {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}+({x_j}^{\&prime;}-{x_{j-1}}^{\&prime;\&prime;})*\sin {{\mathrm{\&alpha;}}_{j-1}}^{\&prime;}+{y_{j-1}}^{\&prime;\&prime;},(j>1)\end{array}\right.---\left(10\right)$

(e)循环处理(a)-(d)，直至段S中的各段处理完毕；

步骤17：对二级线要素l_i2中所有段S＝{s_k|k＝1,2,…,n-1}进行伪装处理；假设段s_k的点集合记为SP＝{sp_w|w＝1,2,…,t},其中t为该段中点的个数，

(a)对段s_k进行平移处理，段s_k中的点集合SP平移后的集合记为SP'＝{sp_w'|w＝1,2,…,t},横纵坐标对应的平移偏量分别为h_w和v_w，根据公式(6)做相应的平移变换；

(b)对平移后的每段进行旋转处理，SP'旋转后的集合记为SP″＝{sp_w″|w＝1,2,…,t}，每段对应的旋转角度为α_W',旋转基点为点集合SP″中对应的点sp_w″,根据公式(9)、(10)做相应的旋转变换；

(c)循环执行(a)和(b)，直到所有段完成伪装处理；

步骤18：生成新图层数据R1和R2，将线要素l_i1中的点集合按原始顺序组成的线要素添加到数据R1中，将线要素l_i2中经伪装变换后的点集合按照原始顺序组成的线要素存放到数据R2中；

步骤19：重复步骤12-18，直到所有线要素完成伪装处理为止；

(二)针对线类型GIS矢量数据还原过程

步骤21：打开伪装后的矢量数据R1和R2，打开密钥文件key，读取用户密钥参数c；

步骤22：分别读取数据R1和数据R2中的线要素L_i1＝{l_i1|i＝1,2,…,m}，L_i2＝{l_i2|i＝1,2,…,m}，其中，m为数据R1中线要素的个数；

步骤23：对线要素l_i2进行分段处理，

(a)线要素L_i1中还原后的点集合记为P_i1＝{p_j|j＝1,2,…,n},n为线要素l_i1中的点个数；

(b)依次读取点集合P_i1中相邻两个点p_j、p_j+1的序号I_j、I_j+1，在线要素l_i2点集合中遍历序号位于区间(I_j,I_j+1)之间的所有点，将这些点分为一段；

(c)循环步骤(b)，直至将线要素l_i2分为(n-1)个分段，记为S＝{s_k|k＝1,2,…,n-1}；

步骤24：确定当前线要素l_i的基线l的直线方程，

(a)假设点集合P_i1的首尾点坐标分别记为p_b(x₁,y₁),p_e(x_n,y_n)，根据公式(1)计算基线斜率k；

(b)假设p_b点经平移后的点p_b'在基线l上，p_b'坐标记为(x₁',y₁')，平移公式如(2)所示，其中c为用户指定参数；

(c)基线l的方程可表示为公式(3)；

步骤25：对点集合P_i1做平移旋转变换使其落在基线l上,处理后的点作为每段点集合的旋转基点，

(a)根据公式(4)、(5)计算当前分段s_k对应的横坐标平移偏移量h_j和纵坐标平移偏移量v_j；

(b)根据公式(6)计算平移后的点集合P_i1'中的点坐标p_j'(x_j',y_j')；

(c)根据公式(1)计算当前段直线p_jp_j+1斜率k_j，根据公式(7)确定基线l和当前段直线P_jP_j+1之间的夹角α_j，当前段直线p_jp_j+1按逆时针旋转与基线重合的旋转角度α_j'由公式(8)计算得出，旋转角度集合记为A＝{α_j'|j＝1,2,…,n}；

(d)根据公式(9)、(10)计算旋转后的点集合p_j″中的点坐标(x_j″,y_j″)；

(e)循环处理(a)-(d)，直至段S中各段处理完毕；

步骤26：对二级线要素l_i2中所有段S＝{s_k|k＝1,2,…,n-1}进行还原处理；假设段s_k的点集合记为SP＝{sp_w|w＝1,2,…,t},sp_w坐标记为(x_w,y_w),其中t为该段中点的个数，

(a)根据公式(9)、(10)对每段进行旋转处理：点集合SP中各点旋转后的点集合记为SP'＝{sp_w'|w＝1,2,…,t},其中每段对应的旋转角度为A'＝{α_j″|j＝1,2,…,n}，由公式(11)计算得出，

α_j″＝2π-α_j'                        (11)

(b)对旋转后的每段进行平移操作：点集合SP'平移后的集合记为SP″＝{sp_w″|w＝1,2,…,t},横纵坐标对应的平移偏量分别为-h_j和-v_j，根据公式(9)做相应的平移变换；

(c)循环执行(a)和(b)，直到完成所有段的还原处理；

步骤27：将线要素l_i2中还原后的点集合与线要素l_i1中的点集合按照坐标序号从小到大的顺序重组成为线要素l_i，完成当前线要素的还原处理；

步骤28：重复步骤22-27，直到完成每个线要素的还原处理。

2.根据权利要求1所述的一种基于信息分存的涉密线状地理要素伪装与还原方法，其特征在于，所述步骤13中可根据安全性需要对一级线要素l_i1进行伪装处理。