CN108629190B - 地理信息数据脱密方法 - Google Patents

地理信息数据脱密方法 Download PDF

Info

Publication number
CN108629190B
CN108629190B CN201810249182.4A CN201810249182A CN108629190B CN 108629190 B CN108629190 B CN 108629190B CN 201810249182 A CN201810249182 A CN 201810249182A CN 108629190 B CN108629190 B CN 108629190B
Authority
CN
China
Prior art keywords
remote sensing
decryption
sensing image
transformation
vector data
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.)
Expired - Fee Related
Application number
CN201810249182.4A
Other languages
English (en)
Other versions
CN108629190A (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.)
Xi'an Transmission And Transformation Project Environmental Influence Control Technology Center Co ltd
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shaanxi Electric Power Co Ltd
Original Assignee
Xi'an Transmission And Transformation Project Environmental Influence Control Technology Center Co ltd
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shaanxi Electric Power Co Ltd
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 Xi'an Transmission And Transformation Project Environmental Influence Control Technology Center Co ltd, State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Shaanxi Electric Power Co Ltd filed Critical Xi'an Transmission And Transformation Project Environmental Influence Control Technology Center Co ltd
Priority to CN201810249182.4A priority Critical patent/CN108629190B/zh
Publication of CN108629190A publication Critical patent/CN108629190A/zh
Application granted granted Critical
Publication of CN108629190B publication Critical patent/CN108629190B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F21/00—Security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
    • G06F21/60—Protecting data
    • G06F21/602—Providing cryptographic facilities or services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F21/00—Security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
    • G06F21/60—Protecting data
    • G06F21/62—Protecting access to data via a platform, e.g. using keys or access control rules
    • G06F21/6209—Protecting access to data via a platform, e.g. using keys or access control rules to a single file or object, e.g. in a secure envelope, encrypted and accessed using a key, or with access control rules appended to the object itself
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • G06T3/147—Transformations for image registration, e.g. adjusting or mapping for alignment of images using affine transformations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G06T7/337—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10032—Satellite or aerial image; Remote sensing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Computer Security & Cryptography (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Bioethics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及一种地理信息数据脱密方法，包括如下步骤：获取目标区域的基准遥感影像；获取目标区域的多个子遥感影像；获取目标区域的矢量数据；通过坐标平移变换，对基准遥感影像进行数据脱密处理，实现基准遥感影像的位置脱密；利用多个子遥感影像上的控制点以及基准遥感影像上的同名控制点，将多个子遥感影像进行配准，完成每个子遥感影像的数据脱密；利用脱密后的子遥感影像和矢量数据上的同名控制点将矢量数据进行空间变换，进而完成矢量数据的空间位置脱密。本发明以基准遥感影像脱密、子遥感影像脱密以及矢量数据脱密展开，提高了脱密精度，且同时实现了脱密后的影像拼接。

Description

地理信息数据脱密方法
技术领域
本发明涉及遥感影像测绘技术领域，具体涉及一种地理信息数据脱密方 法。
背景技术
目前，国家电网公司电网环境保护工作存在部分问题，一定程度上制约 了输变电工程的可研、环评、核准、竣工验收的快速推进，并对工程建成后 的环保、水保验收等方面产生了较大影响，由于缺乏环境敏感区域分布有关 数据，在输变电工程选址选线和建设阶段，部分工程在环境敏感区域避让方 面采取的措施和管理支撑手段不足；输变电工程环境影响预测技术在自动化、 信息化程度方面不高，在选址选线阶段和路径调整时无法快速、及时和准确 评估输变电工程对环境敏感点的影响。
在输变电工程规划、科研以及环评等前期工作中，进行环境敏感区识别， 判断项目可能存在的环境风险时，主要通过现场搜索、踏勘等手段进行，个 别重大项目，会采用遥感、航拍手段进行校核，但由于缺乏敏感区边界信息， 识别中仍存在易遗漏、识别不全的问题。
随着信息技术的快速发展，信息安全形势更不容忽视，以智能终端、 Web2.0、物联网和云计算为代表的新技术，一方面带来了生活和工作的便利 性，另一方面，这些技术的普遍应用同时也带来了更多的安全风险。
测绘地理数据是当今社会化地理数据的基础来源，随着GIS等技术手段的应 用和互联网的发展，地图和地理信息服务等测绘成果被大众广泛熟知和应 用。如何既有效的满足社会对于测绘数据和地理数据的需要，又能保证国家 保密数据的安全性，是现在测绘面临的问题之一。公众版地图是指在国家基 本比例尺地图基础上经过面向国家安全和公众使用双重需求的整合处理而 形成的，它主要供国民经济建设各部门进行规划、设计、科研等使用，亦可 作为社会公众从事与地理信息有关活动的基础底图。公众版地图具有公开性、可用性、权威性等特性。
由于地图公开和电网地理信息保密的双重需求，需要对电网地理信息进 行脱密处理。输变电工程点多面广，线路涉及重要及敏感区域广阔，有水系 的天然水体和水库等人工附属结构；陆地交通、水上交通、航空运输等交通 设施；各种起着基础和主导作用的地形要素及地貌；涉及到国防工程和军事 设施的三角点、独立天文点、发电厂、油、气井、变电所、油库、煤气库等， 以及军事基地等保密单位。电网线路以及重要敏感区域必须进行脱密技术处 理，才能保障电网地理信息及线路涉及敏感区域的信息安全。
发明内容
为保障电网地理信息及线路涉及敏感区域的信息安全，本发明提供了一 种地理信息数据脱密方法，包括如下步骤：
步骤S1：获取目标区域的基准遥感影像；
步骤S2：获取目标区域的多个子遥感影像；
步骤S3：获取目标区域的矢量数据；
步骤S4：通过坐标平移变换，对基准遥感影像进行数据脱密处理，实现 基准遥感影像的位置脱密；
步骤S5：利用多个子遥感影像上的控制点以及基准遥感影像上的同名控 制点，将多个子遥感影像进行配准，完成每个子遥感影像的数据脱密；
步骤S6：利用脱密后的子遥感影像的控制点以及步骤S3中矢量数据上 的同名控制点将矢量数据进行空间变换，进而完成矢量数据的空间位置脱密。
其中，还包括：
步骤S7：对矢量数据进行几何形状脱密；
步骤S8：删除矢量数据中的涉密属性信息。
其中，所述步骤S7中，对矢量数据进行的几何形状脱密包括面平滑处 理及面简化处理；
所述步骤S8中，所删除的涉密属性信息包括军事设施、天然水体、人 工附属结构、交通枢纽、变电站以及换流站。
其中，所述步骤S7中，进行面平滑处理及面简化处理时，依据地理信 息数据脱密原则，对矢量属性表中的重点地物及敏感地物进行处理。
其中，所述步骤S5中，所选择的控制点及同名控制点选自房屋拐角、 道路交叉口或河流交叉口，每个子遥感影像中选取至少6个控制点，并均匀 分布；
所述步骤S6中，所选择的控制点及同名控制点同样选自房屋拐角、道 路交叉口或河流交叉口，并均匀分布。
其中，所述步骤S5中，通过三阶多项式函数或样条曲线函数实现子遥 感影像的数据脱密。
其中，所述三阶多项式函数的变换方法为：
Figure BDA0001605785850000041
Figure BDA0001605785850000042
其中，aij，bij为多项式系数，N为多项式次数，x、y以及X、Y分别 为脱密后及脱密前的影像坐标。
其中，所述步骤S5还包括：对每个子遥感图像进行重采样，并且，重 采样的方法采用三次卷积内插法。
其中，所述步骤S6中，通过射影变换函数完成多个子遥感图像与矢量 数据的融合。
其中，所述步骤S5中，子遥感影像脱密后的位置偏移中误差大于50；同 样的，所述步骤S6中，矢量数据脱密后的位置偏移中误差也大于50。
本发明提供的地理信息数据脱密方法，以基准遥感影像脱密、子遥感影 像脱密以及矢量数据脱密展开，提高了脱密精度，且同时实现了脱密后的影 像拼接、影像平移以及敏感数据的删除。
附图说明
图1：基准遥感影像坐标平移原理图；
图2：基准遥感影像平移实施例图；
图3：一阶多项式变换RMS表；
图4：二阶多项式变换RMS表；
图5：三阶多项式变换RMS表；
图6：样条函数变换RMS表；
图7：待脱密的矢量数据与脱密后的子遥感影像的位置关系图；
图8：矢量数据控制点的全局图；
图9：矢量数据控制点的局部图；
图10：仿射变换RMS表；
图11：射影变换RMS表；
图12：相似变换RMS表；
图13：校正后的矢量数据图；
图14：矢量数据待平滑处理的区域要素图；
图15：矢量数据面平滑处理结果对比图；
图16：矢量数据面简化前的区域图；
图17：矢量数据面简化后的局部放大图。
具体实施方式
为了对本发明的技术方案及有益效果有更进一步的了解，下面结合附图 详细说明本发明的技术方案及其产生的有益效果。
一、基准遥感影像的脱密
对基准遥感影像的脱密通过坐标平移变换进行，如图1所示，坐标系XOY 与坐标系X’O’Y’相应的坐标轴彼此平行，并且具有相同的正向。坐标系X’O’Y’ 是由坐标系XOY平行移动而得到的。设P点在坐标系XOY中的坐标为(x,y)， 在X’O’Y’中坐标为(x’,y’)，而(a，b)是O’在坐标系XOY中的坐标，于是： x＝x’+a,y＝y’+a，即为坐标系平移后坐标关系式。
图2为本发明的基准遥感影像平移实施例图，本发明中，选用北京市主 城区为研究区，通过对基准遥感影像进行平移，实现对基准遥感影像脱密处 理。
二、子遥感影像的脱密
本发明的一较佳实施例中，生成了四个子遥感图像，通过将每个子遥感 图像以及平移后的基准遥感影像加载到GIS软件中，选择二者的同名控制点 点进行地理配准，地理配准提供了多种影像变换方法，本发明通过对多种方 法的误差进行对比，选择了脱密程度最好的三阶多项式函数或样条曲线函 数。
1、子遥感影像空间脱密(配准)方法
(1)多项式变换
在数学方法上,对于不同二维笛卡儿坐标系统间的空间转换,通常采用 的是二元n次多项式,多项式数学表达式为：
Figure BDA0001605785850000061
Figure BDA0001605785850000062
式中：aij，bij为多项式系数，N为多项式次数。
二元n次多项式将不同坐标系统下的对应点坐标联系起来，(x,y)和 (X,Y)分别对应不同坐标系统中的像元坐标。这是一种多项式数字模拟坐 标变换的方法，一旦有了该多项式，就可以从一个坐标系统推算出另一个坐 标系统中的对应点坐标。
不同的二元n次多项式，反映了几何畸变的遥感图像与无几何畸变的遥 感图像间的像元坐标的对应关系，其中哪种多项式是最佳的空间变换模拟式， 能达到图像间坐标的完全配准，是需要考虑和分析的。在二元n次多项式数 字模拟中，从提高几何校正精度的角度考虑，需要兼顾的因素主要有引起几 何畸变的原因和产生数学运算误差因素。
多项式中n值的选择与几何畸变的复杂程度密切相关。
当n＝1，上述的坐标空间变换成为二元一次多项式，可以进行线性的坐 标变换，解决比例尺、中心移动、歪斜等方面的几何畸变，适用于第2级别 以上的遥感数据。n值的不同选择，可以得到不同的空间变换式。
当n≥2，上述的坐标空间变换成为二元非线性多项式，解决遥感器偏航、 俯仰、滚动等因素引起的几何畸变。从理论上讲，n值越大，越能校正复杂 的几何畸变，但计算量也相对要大。实际应用中n值通常取小于等于3。
最后是控制点GCP数目的确定，从数学运算上来说，一次多项式变换， 存在6个系数要计算，需要GCP的最少数目是3。二次多项式变换，有12个系 数需要计算，GCP最少数目是6。n次多项式，GCP的最小数目为(n+1)(n+2) /2。但在实际应用中，采用最小GCP数目几何校正效果往往不好。
(2)样条曲线函数变换
样条曲线函数变换实际上是一种橡皮页变换方法，并对局部精度(而非 全局精度)进行优化。它基于样条函数-一种可维护相邻多项式之间的连续 性和平滑度的分段多项式
Figure BDA0001605785850000081
其中：
j＝1，2，…，N；
N为点数。
δj是通过求解线性方程组而获得的系数；
rj是点(x,y)到第j点之间的距离；
T(x,y)＝a1+a2x+a3y；
其中：
ai是通过求解线性方程组而获得的系数。
Figure BDA0001605785850000082
其中：
R是点与样本之间的距离；
τ2是权重参数；
ko是修正贝塞尔函数；
C是大小等于0.577215的常数。
出于计算目的，输出栅格的整个空间被划分为大小相等的块或区域。x 方向和y方向上的区域数相等，并且这些区域的形状均为矩形。将输入点数 据集中的总点数除以指定的点数值可以确定区域数。如果数据的分布不太均 匀，则这些区域包含的点数可能会明显不同，而点数值只是粗略的平均值。 如果任何一个区域中的点数小于八，则该区域将会扩大到至少包含八个点。
样条曲线函数可将源控制点准确地变换至目标控制点；但不能保证距控 制点的像素距离是准确的。当控制点很重要并且需要进行精确配准的时候就 会用到此变换。添加更多的控制点会提高样条函数变换的总体精度。样条函 数至少需要10个控制点。
(3)变换误差
影像分辨率控制点不可能是绝对同一位置，有的控制点选择误差较大不 但不会提高校正精度，反而会降低影像校正质量，因此要求出所有地面控制 点的均方根误差(RMS),RMS值越大，控制点误差越大，则该控制点需要剔 除。通过调节控制点的位置和数目，多项式系数值和控制点的RMS都在变化， 通常当RMS值小于1时，精度控制在一个像元大小，控制点精度较高。控制 点均方根误差数学表达式为：
Figure BDA0001605785850000091
式中：x，y是地面控制点在原影像上的坐标，x’，y’是变换后影像控制 点坐标。
2、选择控制点
本发明中，利用GIS软件，加载待变换的覆盖实验区4景融合后子遥感影 像(由多光谱影像和全色影像融合生成)。
对4个子遥感影像进行地理配准，首先在待配准影像上选择明显的房屋 拐角、道路交叉口或河流交叉口作为控制点，然后选择基准遥感影像同名点 为控制点进行连接。一般每个待配准的遥感影像选择6个以上控制点，并均 匀分布。
3、分析不同变换函数影像配准残差和RMS误差对比
遥感影像配准最重要的过程是变换函数算法选择，分析不同变换函数对 配准效果影响，以便于选择最合适的变换函数算法。变换函数可利用一阶多 项式、二阶多项式、三阶多项式和样条函数。图3-图6分别是一阶多项式变 换RMS表、二阶多项式变换RMS表、三阶多项式变换RMS表和样条函数变换RMS 表。通过对比一阶多项式、二阶多项式、三阶多项式和样条函数RMS误差表， 发现一阶多项式残差最大16.96，最小1.15，总误差9.55，二阶多项式残差 最大14.53，最小0.88，总误差8.28，三阶多项式残差最大5.18，最小0.07， 总误差2.90，样条函数统计为残差和总误差都为0。残差和RMS总误差三阶多 项式变换和样条函数变换较小，说明利用三阶多项式变换和样条曲线变换配 准更加准确。
4、选择变换函数
变换函数越复杂，脱密程度就越高。考虑到脱密效果和配准精度，样条 函数变换和多项式变换中三阶多项式变换函数已经比较复杂，且配准精度较 高，能够满足空间位置脱密后无缝衔接，变换函数选择样条函数和三阶多项 式。
5、重采样
重采样有三种算法，分别为最近邻、双线性和双三次卷积。最近邻法计 算简单、运算量小，且不破坏原始影像的灰度信息，但是视觉效应差，几何 精度只能达到±0.5个像元，采样后图像的不连续性较突出；双线性内插法 克服了最近邻法使图像不连续的缺点，计算较三次卷积内插法简单，不足的 是图像边缘受到平滑作用，会出现轮廓模糊现象；三次卷积内插法会使图像 边缘增强，但是计算量太大，适合高精度的影像处理要求，因此，本发明采 用三次卷积内插法作为重采样方法。
重复步骤2-5，对每个子遥感影像执行配准，完成子遥感影像的位置脱 密。
6、子遥感影像脱密精度计算
子遥感影像数据进行脱密处理后，利用GIS中的测量功能和excel计算 功能计算数据脱密是否满足国家规定的地图公开脱密位置精度要求。
本发明通过在脱密前后的子遥感影像中选择若干在整个项目区域特征 明显同名点位置，测出脱密前后的位置偏移距离，计算出上述4个子遥感影 像的偏移距离中误差分别为2217.6m、2211.0m、2238.5m以及2207.7m，满 足国家对公共版地图的上线脱密要求。
三、矢量数据的脱密
本发明中，同样将已脱密的子遥感图像及矢量数据加载到GIS软件中， 选择二者的同名控制点进行空间校正(脱密)，空间校正提供了多种方法， 本发明通过对多种方法的误差进行对比，选择了脱密效果最好的射影变换函 数。
1、矢量数据的脱密方法
对矢量数据采用空间变换法进行位置脱密。矢量数据的空间变换法常用 的变换方法有仿射变换、相似变换和射影变换。其中，仿射变换可对图层坐 标进行缩放、倾斜、旋转和平移；相似变换可以缩放、旋转和平移数据，但 不会单独对轴进行缩放，也不会产生任何倾斜，可使变换后的要素保持原有 的横纵比，保持要素的相对形状；射影变换基于更加复杂的公式，对图层坐 标的变换形式更全面。
(1)仿射变换
仿射变换由它对仿射坐标系的作用所完全决定。因为仿射坐标系由原点 和各坐标轴上的单位点所决定，所以一个平面的仿射变换由不共线的三点和 它们的对应点所完全决定。借助于仿射变换对于仿射坐标系的作用，得到仿 射变换的坐标表示。在定义仿射变换时，采用的是纯几何的方式，直到研究 了仿射变换的简单性质以后，才导出了仿射变换的坐标表示。正交变换和相 似变换也是这样处理的。其实，也可以从仿射变换的坐标表示出发，定义仿 射变换，即仿射变换是由下列公式定义的：
Figure BDA0001605785850000121
其中，
Figure BDA0001605785850000122
相似变换可以分解为放缩、平移、旋转和翻转变换的复合。相似变换是 仿射变换的一种特殊情况，也就是在仿射变换中去除错位变换这个因子后的 结果。
(2)射影变换
由有限次中心射影的积定义的两条直线间的一一对应变换称为一维射影 变换。由有限次中心射影的积定义的两个平面之间的一一对应变换称为二维 射影变换。
在射影变换方法中，设在点场π与π，上各建立了齐次射影坐标系， x＝(x1,x2,x3)，x’＝(x’1,x’2,x’3)分别是其上点的射影坐标，如果映射
φ:π→π’可以表示成如下公式：
Figure BDA0001605785850000131
|A|＝|aij|≠0,ρ≠0
则称φ是π到π’的二维射影对应；特别地，若π＝π’，则称φ是π上 的二维射影变换。
均方根误差(root-mean-square error)亦称标准误差，其定义为
Figure BDA0001605785850000132
在有限测量次数中，均方根误差常用下式表示:
Figure BDA0001605785850000133
公式中，n为测量次数；d为一组测量值与平均值的偏差。
2、选择控制点
本发明中，同于子遥感影像的脱密，对矢量数据与子遥感影像进行配准 时，同样于矢量数据及子遥感影像上选择同样位置的明显的房屋拐角、道路 交叉口或河流交叉口作为控制点，一般每个待配准的子遥感影像选择6个以 上控制点，并均匀分布。图7所示为待脱密的矢量数据与脱密后的子遥感影 像的位置关系图，图8及图9分别为矢量数据控制点的全局图和局部图。
3、分析不同变换函数矢量配准残差和RMS误差对比
变换函数可利用仿射变换、射影变换和相似变换。图10-图12分别是仿 射变换RMS表、射影变换RMS表和相似变换RMS表。通过对比仿射变换、射 影变换和相似变换函数RMS误差表，从三种变换方法RMS误差表中可知，仿 射变换RMS误差为2.058，射影变换RMS误差为1.985，相似变换RMS误差 为2.316,射影变换RMS误差最小。
4、选择变换函数
仿射变换、射影变换和相似变换三种变换方法中，射影变换函数相对于 仿射变换和相似变换更复杂，脱密程度最高。在配准精度方面，射影变换也 高于仿射变换和相似变换，能够满足空间位置脱密后无缝衔接，故本发明中， 变换函数选择射影变换函数。
如图13所示，为本发明校正后的矢量数据图，采用射影变换函数对矢 量数据进行校正，既实现了位置脱密，又能较好的与遥感影像契合。
5、矢量数据脱密精度计算
矢量数据脱密处理之后，利用GIS中的测量功能和excel计算功能计算 数据脱密是否满足国家规定地图公开脱密位置精度要求。
(1)分别通过在脱密前后的矢量数据中选择若干在整个项目区域特征 明显同名点位置(如道路交叉口，河流交叉口等)，用GIS中的测量功能测 量出脱密前后同名点位置偏移距离；
(2)将统计的偏移距离利用excel按照中误差计算公式：
Figure BDA0001605785850000141
计算出的中误差m，当m大于50m时满足脱密要求，中误差m通过对n个同 名点偏移距离求平方和的平均值并开根号求得，计算得到2225.596m，满足 国家对公众版地图上线脱密要求。
6、矢量数据局部要素脱密实施方案
对矢量数据的地理信息数据进行了空间位置脱密后，矢量数据中的重要 地物和敏感地物还需要进行几何形状脱密，本发明中，通过对局部地物数据 进行平滑、简化等处理，以对图形特征进行钝化，降低重点和敏感地理数据 的几何精度。
(1)面平滑实施方法
依据地理信息数据脱密原则，根据矢量属性表，找出重点地物和敏感地 物，如水库、变电所、油库、煤气库、涉及国防工程和军事设施等，选择这 些需要平滑处理的区域要素，如图14所示。
可利用ARCGIS软件中面平滑工具进行处理，得到处理结果，如图15所 示。
(2)面简化实施方法
与面平滑处理类似，依据地理信息数据脱密原则，根据矢量属性表，找 出重点地物和敏感地物，选择这些需要简化处理的区域要素。
可利用ARCGIS软件中面平滑工具进行处理。简化容差中，最大允许偏 移量可以根据需要进行调整，本发明中设置为10米，最小面积根据需要进 行设置，这里设置为100平方米。如图16及图17所示，分别为面简化前的 区域图以及面简化后的局部放大图。
7、矢量数据属性脱密实施方案
根据公开地图要求，矢量数据牵涉到保密信息的需要删除，防止信息泄 露，本发明中，通过在矢量属性表中查询有军事设施、天然水体、人工附属 结构、交通枢纽、变电站、换流站等需要删除的属性信息，可利用GIS软件、 编写程序等方法将涉密属性信息删除，
8、地理信息数据脱密原则
完全删除地理数据中含有涉密的地理要素或像素及其他属性信息。
(1)水系的天然水体和人工附属结构都是涉密的重要信息，因此应删 除水库库容和大坝的高度、河宽、水深、流速、底质和岸质属性等。
(2)居民地是重要的地物要素，从涉密角度看，保密单位的地理位置、 分布特征、编制及部署均应进行保密处理，按不依比例尺的居民地符号式样 表示或不表示。
(3)交通运输条件是涉密重要组成部分。陆地交通、水上交通、航空 运输详细信息对其影响非常重要，应删除交通枢纽的地理位置、分布和特征， 桥梁、隧道的位置、长度、宽度、高度、性质、载重量、运输能力、周围地 形状况等不易表达的信息。
(4)在任何形态的地形中，起着基础和主导作用的地形要素是地貌。 地形景观与特征、图幅的最高点、主要山峰等制高点的地理位置及高程、冲 沟比高、陡崖比高等影响较大的地貌信息应删除。
(5)对潜力目标进行定位和破坏，已成为现代战争的一个重要特点。 总体而言，国防工程中的地理位置、分布特征等属性，均不得在公开地图产 品上表示。如三角点、独立天文点不表示，发电厂按不依比例尺的居民地符 号表示，油、气井、变电所、油库、煤气库、涉及国防工程和军事设施的山 洞不表示，高压输电线、油气管道不表示等。
本发明的地理信息数据脱密方法，以基准遥感影像脱密、子遥感影像脱 密以及矢量数据脱密展开，其有益效果如下：
(1)通过平移方法将基准遥感影像进行平移，使其平移距离远远大于 国家规定对公众版上线地图数据精度要求，为待脱密遥感影像位置脱密奠定 了精度要求。
(2)通过采用地理配准的方法对待脱密子遥感影像进行配准，既保障 了待脱密子遥感影像无缝拼接问题，又实现了脱密。通过对比不同影像配准 算法，根据算法越复杂，脱密越安全的原则，选择了相对其它变换方法更复 杂的三次多项式变换方法和样条函数变换方法，提高了配准影像精度。
(3)通过以脱密后的子遥感影像为基准进行矢量数据脱密，采用空间 校正方法对矢量进行校正，既实现了矢量脱密，又能保障与脱密后子遥感影 像较好的契合。通过选择校正误差更小的射影变换函数，提高了矢量数据的 脱密程度。
(4)通过利用平滑和简化方法对矢量要素中的重点区域或地物进行几 何精度处理，使其几何形状发生畸变，从而达到矢量数据局部要素脱密的目 的。
(5)通过删除敏感地物要素属性，从而实现矢量数据属性脱密，结合 使用降低敏感地物几何形状方法，能够达到更好的脱密效果。
虽然本发明已利用上述较佳实施例进行说明，然其并非用以限定本发明 的保护范围，任何本领域技术人员在不脱离本发明的精神和范围之内，相对 上述实施例进行各种变动与修改仍属本发明所保护的范围，因此本发明的保 护范围以权利要求书所界定的为准。

Claims (10)

1.一种地理信息数据脱密方法，其特征在于包括如下步骤：
步骤S1：获取目标区域的基准遥感影像；
步骤S2：获取目标区域的多个子遥感影像；
步骤S3：获取目标区域的矢量数据；
步骤S4：通过坐标平移变换，对基准遥感影像进行数据脱密处理，实现基准遥感影像的位置脱密；
步骤S5：利用多个子遥感影像上的控制点以及平移后的基准遥感影像上的同名控制点，将多个子遥感影像进行配准，完成每个子遥感影像的数据脱密；
步骤S6：利用脱密后的子遥感影像的控制点以及步骤S3中矢量数据上的同名控制点将矢量数据进行空间变换，进而完成矢量数据的空间位置脱密。
2.如权利要求1所述的地理信息数据脱密方法，其特征在于还包括：
步骤S7：对矢量数据进行几何形状脱密；
步骤S8：删除矢量数据中的涉密属性信息。
3.如权利要求2所述的地理信息数据脱密方法，其特征在于：
所述步骤S7中，对矢量数据进行的几何形状脱密包括面平滑处理及面简化处理；
所述步骤S8中，所删除的涉密属性信息包括军事设施、天然水体、人工附属结构、交通枢纽、变电站以及换流站。
4.如权利要求3所述的地理信息数据脱密方法，其特征在于：所述步骤S7中，进行面平滑处理及面简化处理时，依据地理信息数据脱密原则，对矢量属性表中的重点地物及敏感地物进行处理。
5.如权利要求1所述的地理信息数据脱密方法，其特征在于：
所述步骤S5中，所选择的控制点及同名控制点选自房屋拐角、道路交叉口或河流交叉口，每个子遥感影像中选取至少6个控制点，并均匀分布；
所述步骤S6中，所选择的控制点及同名控制点同样选自房屋拐角、道路交叉口或河流交叉口，并均匀分布。
6.如权利要求1所述的地理信息数据脱密方法，其特征在于：所述步骤S5中，通过二元N次多项式或样条曲线函数实现子遥感影像的数据脱密。
7.如权利要求6所述的地理信息数据脱密方法，其特征在于：如果采用二元N次多项式，则二元N次多项式的变换方法为：
Figure FDA0002436490220000021
Figure FDA0002436490220000022
其中，aij，bij为多项式系数，N为多项式次数，x、y以及X、Y分别为脱密后及脱密前的影像坐标。
8.如权利要求6所述的地理信息数据脱密方法，其特征在于：所述步骤S5还包括：对每个子遥感图像进行重采样，并且，重采样的方法采用三次卷积内插法。
9.如权利要求1所述的地理信息数据脱密方法，其特征在于：所述步骤S6中，通过射影变换函数完成多个子遥感图像与矢量数据的融合。
10.如权利要求1-9中任一项所述的地理信息数据脱密方法，其特征在于：所述步骤S5中，子遥感影像脱密后的位置偏移中误差大于50；同样的，所述步骤S6中，矢量数据脱密后的位置偏移中误差也大于50。
CN201810249182.4A 2018-03-23 2018-03-23 地理信息数据脱密方法 Expired - Fee Related CN108629190B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810249182.4A CN108629190B (zh) 2018-03-23 2018-03-23 地理信息数据脱密方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810249182.4A CN108629190B (zh) 2018-03-23 2018-03-23 地理信息数据脱密方法

Publications (2)

Publication Number Publication Date
CN108629190A CN108629190A (zh) 2018-10-09
CN108629190B true CN108629190B (zh) 2020-08-25

Family

ID=63696319

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810249182.4A Expired - Fee Related CN108629190B (zh) 2018-03-23 2018-03-23 地理信息数据脱密方法

Country Status (1)

Country Link
CN (1) CN108629190B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109597870A (zh) * 2018-12-07 2019-04-09 国网安徽省电力有限公司黄山供电公司 一种地理数据可用性的评定方法
CN109885775B (zh) * 2019-03-05 2020-03-17 重庆工商大学融智学院 一种实现生态环境空间大数据集成共享的方法
CN110427450A (zh) * 2019-08-09 2019-11-08 深圳四维集思技术服务有限公司 基于互联网地图的油气管线数据处理方法和系统
CN111028178B (zh) * 2019-12-20 2022-04-29 武汉大学 一种基于深度学习的遥感影像数据自动化几何纠正方法
CN112749241A (zh) * 2020-07-28 2021-05-04 广西华蓝岩土工程有限公司 一种涉密宗地图界址点大地坐标脱密处理方法
CN117195249B (zh) * 2023-08-04 2024-04-09 武汉大学 基于差分隐私的矢量地理数据脱密处理方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101685539A (zh) * 2009-06-29 2010-03-31 中国测绘科学研究院 一种遥感影像在线正射纠正方法和系统
CN103067159A (zh) * 2012-12-28 2013-04-24 南京师范大学 一种gis矢量数据可逆脱密方法
CN104361590A (zh) * 2014-11-12 2015-02-18 河海大学 一种控制点自适应分布的高分辨率遥感影像配准方法
CN205318382U (zh) * 2016-01-18 2016-06-15 耿翀 基于arcgis的脱密地理信息互联网发布引擎系统
CN106778347A (zh) * 2016-12-22 2017-05-31 南京师范大学 一种基于三角函数的矢栅地理数据可逆脱密方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6768942B1 (en) * 2000-09-18 2004-07-27 Navigation Technologies Corp. Navigation system with decryption functions and secure geographic database

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101685539A (zh) * 2009-06-29 2010-03-31 中国测绘科学研究院 一种遥感影像在线正射纠正方法和系统
CN103067159A (zh) * 2012-12-28 2013-04-24 南京师范大学 一种gis矢量数据可逆脱密方法
CN104361590A (zh) * 2014-11-12 2015-02-18 河海大学 一种控制点自适应分布的高分辨率遥感影像配准方法
CN205318382U (zh) * 2016-01-18 2016-06-15 耿翀 基于arcgis的脱密地理信息互联网发布引擎系统
CN106778347A (zh) * 2016-12-22 2017-05-31 南京师范大学 一种基于三角函数的矢栅地理数据可逆脱密方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
拓扑保形的矢量地理数据几何脱密方法;吕海洋，等;《中国矿业大学学报》;20170531(第3期);第648-654页 *

Also Published As

Publication number Publication date
CN108629190A (zh) 2018-10-09

Similar Documents

Publication Publication Date Title
CN108629190B (zh) 地理信息数据脱密方法
Yue et al. High-accuracy surface modelling and its application to DEM generation
Uggla et al. Geographic capabilities and limitations of Industry Foundation Classes
KR100686287B1 (ko) 공간정보/위치정보의 정밀 위치정보 변환을 위한 변동량모델링 방법
Doytsher A rubber sheeting algorithm for non-rectangular maps
CN111696156B (zh) 一种免控制点的遥感影像坐标转换方法
Tong et al. Area‐preservation simplification of polygonal boundaries by the use of the structured total least squares method with constraints
Javadnejad et al. An assessment of UAS-based photogrammetry for civil integrated management (CIM) modeling of pipes
Di Benedetto et al. Integration of LiDAR data into a regional topographic database for the generation of a 3D city model
Tereshchuk et al. Efficiency of application of satellite technology when performing land and cadastral works in settlements
Tong et al. A least squares-based method for adjusting the boundaries of area objects
Arias et al. ORTHOIMAGE‐BASED DOCUMENTATION OF ARCHAEOLOGICAL STRUCTURES: THE CASE OF A MEDIAEVAL WALL IN PONTEVEDRA, SPAIN
CN115409962B (zh) 虚幻引擎内构建坐标系统的方法、电子设备和存储介质
Ma et al. Low‐Altitude Photogrammetry and Remote Sensing in UAV for Improving Mapping Accuracy
CN115409960A (zh) 基于虚幻引擎的模型构建方法、电子设备和存储介质
CN109239745A (zh) 一种高分辨率卫星影像及有理多项式参数变换方法
Tong et al. Optimal selection of virtual control points with planar constraints for large‐scale block adjustment of satellite imagery
Crosilla et al. Fifty years of the Friuli Venezia Giulia regional technical map: a best practice in the Italian cartographic context
Moscovici et al. Mathematical cartography based on georeferencing maps
US9958273B2 (en) Computer-implemented map of assets conversion method and system
Mohammed et al. Spatial Correction of Cadastral Maps Based on Satellite Imagery
Oh et al. Large Scale Topographic Mapping Using Stereo Pair Imagery from CAS500-1 Satellite
Siriba et al. Mutual linear feature matching and alignment designed for geometric accuracy enhancement of graphical cadastral datasets
Chen et al. Optimized patch backprojection in orthorectification for high resolution satellite images
Elshehaby et al. Assessment of cartographic potential of EgyptSat-1 satellite image (case study in flat areas)

Legal Events

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

Granted publication date: 20200825

Termination date: 20210323

CF01 Termination of patent right due to non-payment of annual fee