CN105302874A - 基于地理云数据的空间匹配方法 - Google Patents
基于地理云数据的空间匹配方法 Download PDFInfo
- Publication number
- CN105302874A CN105302874A CN201510648550.9A CN201510648550A CN105302874A CN 105302874 A CN105302874 A CN 105302874A CN 201510648550 A CN201510648550 A CN 201510648550A CN 105302874 A CN105302874 A CN 105302874A
- Authority
- CN
- China
- Prior art keywords
- spatial
- cloud data
- point
- matching
- scan model
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明揭示了一种基于地理云数据的空间匹配方法,包括:步骤1,建立分布式空间数据库,并采用ICP算法进行点云数据配准,在点云配准过程中,采用多点的最小二乘拟合法进行平差,求取转换矩阵;步骤2,进行同名实体的目标匹配。本发明能够快速查询和调用地理信息空间位置信息和属性信息,并通过进一步运用同名实体匹配方法来实现空间数据的精准匹配;同时本发明还采用最小二乘法拟合,通过三维空间中的多个点求取转换矩阵,计算方法简单,可以有效提高配准精度。
Description
技术领域
本发明涉及一种空间匹配技术,尤其是涉及一种基于地理云数据的空间匹配方法。
背景技术
随着信息社会的到来,人类获取空间信息的能力较过去有了明显提高,面对海量信息,人们对于信息的需求要求发生巨大变化,对信息的广泛性、精准性、快速性及综合性要求越来越高。
传统的GIS(GeographicInformationSystem,地理信息系统)数据库往往都是单独存在,信息孤岛现象严重,各部门根据自己的具体情况和应用需求不同,从而导致GIS数据在结构数据,数据模型、数据格式上存在差异,数据难以共享。虽然现阶段WebGIS(网络地理信息系统)也能够实现异构地理空间数据的共享需求,但现阶段基于语义匹配的研究还有待进一步完善。
众所周知,数据的采集和更新是GIS工程中最耗时耗力的一项工作,许多GIS应用系统存在数据更新周期过长等问题,这是由于空间数据标准化发展的滞后,许多空间数据在数据模型和格式以及比例尺,地图投影等方面存在差异,为了充分利用已有信息,加快数据采集的高效性,这就需要基于地理数据库的空间匹配技术,其定义就是将同一实体不同数据库的空间位置信息和属性信息相互匹配或者同一地区两个以上不同来源的数据库合并生成更具完善的地理信息数据库。
基于目标空间匹配主要是通过对目标几何、拓扑和语义进行相似性度量,识别出同一地区的空间数据集中的同义地物,从而建立两个空间数据集之间的相互关联。匹配技术在计算机识别、模式识别、图像分析、图像理解等领域已经运用,现有的匹配方法主要基于多种类型的图像特征来进行图像匹配技术,缺乏完善的协同智能化的方法。
目前,解决不同点云数据集配准问题的直接方法有在被测物体上贴固定球和使用高精度的旋转工作台,固定球会使部分物体不可见,一般应贴在相对较平坦的区域流动式光学扫描仪Atos是在物体上或固定物体的工具上贴特制的圆形标签,其作用相当于固定球。根据前后2个视角观察的3个或3个以上不共线的公共标签来对数据进行配准,利用旋转工作台可直接对测量数据进行配准,但是物体底部和下部的数据点是无法采集的。
针对复杂的匹配问题,一方面设计更加高效的搜索算法引入更多的图结构和语义信息降低匹配问题复杂度,另一方面采用近似算法允许匹配结果中存在一定概率的错误保证算法达到多项式级的时间复杂度。在实际应用中,由于图数据的形式、结构、规模以及包含信息不同又产生了适合不同应用场景的扩展图匹配问题。
点云配准是将2个或2个以上坐标系中的大容量3维空间数据点集转换到统一坐标系统中的数学计算过程。ICP算法是一种迭代算法,具有很高的匹配精度。但由算法描述可知,它存在计算量较大且迭代过程可能无法收敛到全局最优解的缺陷。ICP算法最耗费时间的步骤是求解邻近点对的过程,因为它采用的是全局搜索。为适应不同的环境并克服ICP算法自身的部分缺陷,近十几年,研究学者为改善ICP(IterativeClosestPoint,迭代最近点算法)算法效率和精度,对搜索最近点过程和加快收敛速度方面进行了改进。
发明内容
本发明的目的在于克服现有技术的缺陷,提供一种基于地理云数据的空间匹配方法,采用分布式网络一体化集成管理,及对ICP算法的改进,以实现地理信息数据的实时更新管理,区域地物实时快速精准匹配。
为实现上述目的,本发明提出如下技术方案:一种基于地理云数据的空间匹配方法,包括:
建立分布式空间数据库,并采用ICP算法进行点云数据配准;
进行同名实体的目标匹配;
其中,所述点云数据配准过程包括:
在两个测试单元内分别得到各自的公共点三维空间坐标OAXAYAZA、OBXBYBZB,且两者具有的坐标变换关系为:
采用多点的最小二乘拟合法进行平差,求取转换矩阵MBA;
采用比较点间欧式距离的方法进行点对匹配;
所述同名实体的目标匹配过程包括:
集中搜集与待匹配对象存在匹配关系的候选匹配对象;
对待匹配对象与所述候选匹配对象进行相似性评价,以确定它们之间的匹配关系。
优选地,所述分布式空间数据库包括分布式多空间数据管理系统,所述分布式多空间数据管理系统存在于多个空间数据库管理系统和非空间数据库管理系统之上。
优选地,所述分布式多空间数据管理系统结构采用Client/Server结构。
优选地,采用多点的最小二乘拟合法进行平差,求取转换矩阵MBA的过程包括:
采用4个或4个以上的点,将两个坐标变换关系写成为:
式中,XA和XB包括两个坐标系下4个或4个以上点的齐次坐标,则转换矩阵MBA的解为:
优选地,采用比较点间欧式距离的方法进行点对匹配的过程包括:
初始化迭代,迭代次数为k,设置欧式距离均方差阀值r以及最大迭代次数Kmax;
对待匹配扫描模型中的每个点搜索其在参考扫描模型中欧式距离最近的点,生成整个扫描模型的邻近点对集合Yk;
由所述参考扫描模型中的扫描点和匹配点对集合,计算扫描模型间的变换矩阵Pk=(PRk|PTk)T和匹配点对间的欧式距离误差dk;
根据变换矩阵Pk=(PRk|PTk)T变换当前扫描模型Sr中的所有扫描点位置;
计算变换后,参考扫描模型与当前扫描模型间的对应点对间的欧式距离误差dk+1,并计算第k和k+1次迭代中误差变化量,当该变化量小于欧式距离均方差阀值r或迭代次数k大于kmax时,停止迭代。
优选地,转换矩阵MBA中共有12个未知量,且旋转矩阵R:
满足的正交条件为:
本发明的有益效果是:
1、本发明主要采用建立分布式空间数据库系统,实现分布式计算,从而快速查询和调用地理信息空间位置信息和属性信息,通过进一步运用同名实体匹配方法来实现空间数据的精准匹配。
2、采用最小二乘法拟合,通过三维空间中的多个点求取转换矩阵,计算方法简单,可以有效提高配准精度。
附图说明
图1是本发明基于地理云数据的空间匹配方法的流程示意图;
图2是本发明建立的分布式空间数据库的原理示意图;
图3a和图3b是本发明模式图P在数据图G中存在的两个与其匹配的子图的原理示意图。
具体实施方式
下面将结合本发明的附图,对本发明实施例的技术方案进行清楚、完整的描述。
如图1所示,本发明所揭示的一种基于地理云数据的空间匹配方法,包括:
步骤一,建立分布式空间数据库,并采用ICP算法进行点云数据配准。
结合图1和图2所示,本发明建立的分布式多空间数据库的体系结构采用了Client/Serve结构,它是完全独立的参与空间数据库的管理系统,各自管理各自的数据库。分布式多空间数据库管理系统是在这些独立的空间数据库管理系统和非空间数据库管理系统之上运行的一层软件,它负责管理全局的控制信息,包括全局模式、全局元数据、提交和控制涉及不止一个参与空间数据库的全局查询和全局事务。多空间数据库系统没有对参与空间数据库系统做出任何改动,全局用户可以透明地访问分布的异构的空间数据源。
多空间数据库系统是已经存在的若干个空间数据库之上,为全局用户提供一个统一的存储空间数据环境,并且规定了本数据由本地拥有和管理,采用的分割式组织方式,按照某种逻辑划分在分布式数据库中。
步骤一中,点云数据配准过程具体包括以下步骤:
在两个测试单元内分别得到各自的公共点三维空间坐标0AxAYAzA、0BxBYBzB,设测量单元l的坐标系为0AxAYAzA,测量单元2的坐标系为0BxBYBzB,且二者具有的坐标变换关系为:
即
式中,转换矩阵MBA中共有12个未知量,而旋转矩阵
满足的正交条件为:
所以独立变量只有6个。因此,如果已知3个点在2个坐标系下的三维坐标即可计算出二者的转换矩阵,进而实现三维数据的配准。现在考虑在已知 和i之3时,如何求解转换矩阵MBA。若采用3点在2测量单元坐标系的坐标值,要充分利用旋转矩阵的正交条件,求解非线性方程,计算量较大,且引入计算误差。所以本发明采用多点的最小二乘拟合法进行平差,如采用4个或4个以上的点,写成为:
式中XA和XB包括2坐标系下4个或4个以上点的齐次坐标,则矩阵MBA的解为:
所以本阀门采用最小二乘法拟合,通过三维空间中的多个点求取转换矩阵,计算方法简单,可以有效提高配准精度。
ICP算法在进行点对匹配时,设扫描匹配过程中的参考扫描模型为Sr,待匹配的当前扫描模型为Sc,模型中的数据点个数为N,参考扫描模型Sr和当前扫描模型间Sc的变换矩阵为P=(p0,p1,p2,p3,p4,p5,p6)T,其中,PR=(p0,p1,p2,p3)T四元组表示旋转偏移量,可转化为3×3的旋转矩阵R(PR),PT=(p4,p5,p6)T表示平移偏移量,这样,模型间的变换矩阵可以表示为P=(PR|PT)T。
本发明中,ICP算法在进行点对匹配时采用的是比较点间欧式距离的方法,具体过程包括:
第一步:初始化迭代,P0=(1,0,0,0,0,0,0)T,迭代次数k,设置欧式距离均方差阀值r以及最大迭代次数Kmax。
第二步:迭代步骤:
(1)对待匹配扫描模型中的每个点搜索其在参考扫描模型中欧式距离最近的点,生成整个扫描模型的邻近点对集合Yk。
(2)由参考模型中的扫描点和匹配点对集合计算扫描模型间的变换矩阵Pk=(PRk|PTk)T和匹配点对间的欧式距离误差dk。
(3)根据变换矩阵Pk=(PRk|PTk)T变换当前扫描模型Sr中的所有扫描点位置。
(4)计算变换后,参考扫描模型与当前扫描模型间的对应点对间的欧式距离误差dk+1,并计算第k和k+1次迭代中误差变化量。当该变化量小于欧式距离均方差阀值r或迭代次数k大于kmax时,停止迭代。
精确子图匹配,已知一个数据图G=(V,E,L)和一个模式图P=(Vp,Ep,Lp),精确图匹配问题要求在模式图和数据子图Gsub=(Vsub,Esub,Lsub)(Gsub属于G)之间建立一个双射函数f:
精确图匹配问题需要在数据图中找到所有满足上述条件的子图作为模式图的匹配结果。图3a和图3b分别展示了模式图P在数据图G中存在的两个与其匹配的子图M1=[(v1,u3),(v2,u5),(v3,u9),(v4,u6)]和M2=[(v1,u6),(v2,u5),(v3,u9),(v4,u3)]。
步骤二,同名实体的目标匹配方法的建立
同名实体的匹配过程主要包含两个关键步骤:1、确定匹配候选集2、基于相似性评价的匹配判断。第一步负责在其他数据集中搜集与待匹配对象可能存在匹配关系的候选匹配对象,第二步负责对待匹配对象与候选匹配对象进行相似性评价,以确定他们之间的匹配关系。
本发明基于地理信息系统技术,以云计算技术为核心,采用分布式网络一体化集成管理,实现地理信息数据的实时更新管理。实现网络资源整合,智能终端实时传送,区域地物实时快速精准匹配。
该平台构建的系统可以整合政府各部门数据资源,实现网格化资源管理,形成较为完善的智慧城市管理信息化体系。系统依靠对业务流程信息的搜集,并整合各部门资源信息,应用大数据融合,挖掘数据价值,以实现服务方式的转变,资源配置的优化,业务流程的改进,从而驱动服务模式的创新,使决策科学有据。
本发明的技术内容及技术特征已揭示如上,然而熟悉本领域的技术人员仍可能基于本发明的教示及揭示而作种种不背离本发明精神的替换及修饰,因此,本发明保护范围应不限于实施例所揭示的内容,而应包括各种不背离本发明的替换及修饰,并为本专利申请权利要求所涵盖。
Claims (6)
1.一种基于地理云数据的空间匹配方法,其特征在于,包括:
建立分布式空间数据库,并采用ICP算法进行点云数据配准;
进行同名实体的目标匹配;
其中,所述点云数据配准过程包括:
在两个测试单元内分别得到各自的公共点三维空间坐标OAXAYAZA、OBXBYBZB,且两者具有的坐标变换关系为:
其中,MBA为转换矩阵,R为旋转矩阵;
采用多点的最小二乘拟合法进行平差,求取转换矩阵MBA;
采用比较点间欧式距离的方法进行点对匹配;
所述同名实体的目标匹配过程包括:
集中搜集与待匹配对象存在匹配关系的候选匹配对象;
对待匹配对象与所述候选匹配对象进行相似性评价,以确定它们之间的匹配关系。
2.根据权利要求1所述的基于地理云数据的空间匹配方法,其特征在于,所述分布式空间数据库包括分布式多空间数据管理系统,所述分布式多空间数据管理系统存在于多个空间数据库管理系统和非空间数据库管理系统之上。
3.根据权利要求2所述的基于地理云数据的空间匹配方法,其特征在于,所述分布式多空间数据管理系统结构采用Client/Server结构。
4.根据权利要求1所述的基于地理云数据的空间匹配方法,其特征在于,采用多点的最小二乘拟合法进行平差,求取转换矩阵MBA的过程包括:
采用4个或4个以上的点,将两个坐标变换关系写成为:
式中,XA和XB包括两个坐标系下4个或4个以上点的齐次坐标,则转换矩阵MBA的解为:
5.根据权利要求1所述的基于地理云数据的空间匹配方法,其特征在于,采用比较点间欧式距离的方法进行点对匹配的过程包括:
初始化迭代,迭代次数为k,设置欧式距离均方差阀值r以及最大迭代次数Kmax;
对待匹配扫描模型中的每个点搜索其在参考扫描模型中欧式距离最近的点,生成整个扫描模型的邻近点对集合Yk;
由所述参考扫描模型中的扫描点和匹配点对集合,计算扫描模型间的变换矩阵Pk=(PRk|PTk)T和匹配点对间的欧式距离误差dk;
根据变换矩阵Pk=(PRk|PTk)T变换当前扫描模型Sr中的所有扫描点位置;
计算变换后,参考扫描模型与当前扫描模型间的对应点对间的欧式距离误差dk+1,并计算第k和k+1次迭代中误差变化量,当该变化量小于欧式距离均方差阀值r或迭代次数k大于kmax时,停止迭代。
6.根据权利要求1所述的基于地理云数据的空间匹配方法,其特征在于,所述转换矩阵MBA中共有12个未知量,且旋转矩阵R:
满足的正交条件为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510648550.9A CN105302874A (zh) | 2015-10-09 | 2015-10-09 | 基于地理云数据的空间匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510648550.9A CN105302874A (zh) | 2015-10-09 | 2015-10-09 | 基于地理云数据的空间匹配方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105302874A true CN105302874A (zh) | 2016-02-03 |
Family
ID=55200144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510648550.9A Pending CN105302874A (zh) | 2015-10-09 | 2015-10-09 | 基于地理云数据的空间匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105302874A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679126A (zh) * | 2017-09-21 | 2018-02-09 | 北京星闪世图科技有限公司 | 激光三维点云数据存储和管理方法及其系统 |
CN108389221A (zh) * | 2018-01-30 | 2018-08-10 | 深圳市菲森科技有限公司 | 三维图像的扫描方法及系统 |
CN109100730A (zh) * | 2018-05-18 | 2018-12-28 | 北京师范大学-香港浸会大学联合国际学院 | 一种多车协同快速建图方法 |
CN112484729A (zh) * | 2020-11-11 | 2021-03-12 | 深圳市优必选科技股份有限公司 | 导航地图切换方法、装置、终端设备及存储介质 |
CN112529945A (zh) * | 2020-11-17 | 2021-03-19 | 西安电子科技大学 | 一种多视角三维isar散射点集配准方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005068936A1 (en) * | 2004-01-16 | 2005-07-28 | Geotango International, Corp. | System,computer program and method for 3d object measurement, modeling and mapping from single imagery |
CN104330074A (zh) * | 2014-11-03 | 2015-02-04 | 熊友谊 | 一种智能测绘平台及其实现方法 |
-
2015
- 2015-10-09 CN CN201510648550.9A patent/CN105302874A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005068936A1 (en) * | 2004-01-16 | 2005-07-28 | Geotango International, Corp. | System,computer program and method for 3d object measurement, modeling and mapping from single imagery |
CN104330074A (zh) * | 2014-11-03 | 2015-02-04 | 熊友谊 | 一种智能测绘平台及其实现方法 |
Non-Patent Citations (5)
Title |
---|
于静等: ""大规模图数据匹配技术综述"", 《计算机研究与发展》 * |
吴禄慎等: ""基于特征点的改进ICP三维点云配准技术"", 《南昌大学学报·工科版》 * |
翟仁建: ""基于全局一致性评价的多尺度矢量空间数据匹配方法研究"", 《中国博士学位论文全文数据库基础科学辑(月刊)》 * |
邬伦等: ""分布式多空间数据库系统的集成技术"", 《地理学与国土研究》 * |
陈三清: ""点云数据的配准方法研究"", 《光盘技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679126A (zh) * | 2017-09-21 | 2018-02-09 | 北京星闪世图科技有限公司 | 激光三维点云数据存储和管理方法及其系统 |
CN108389221A (zh) * | 2018-01-30 | 2018-08-10 | 深圳市菲森科技有限公司 | 三维图像的扫描方法及系统 |
CN108389221B (zh) * | 2018-01-30 | 2021-10-01 | 深圳市菲森科技有限公司 | 三维图像的扫描方法及系统 |
CN109100730A (zh) * | 2018-05-18 | 2018-12-28 | 北京师范大学-香港浸会大学联合国际学院 | 一种多车协同快速建图方法 |
CN109100730B (zh) * | 2018-05-18 | 2022-05-24 | 北京师范大学-香港浸会大学联合国际学院 | 一种多车协同快速建图方法 |
CN112484729A (zh) * | 2020-11-11 | 2021-03-12 | 深圳市优必选科技股份有限公司 | 导航地图切换方法、装置、终端设备及存储介质 |
CN112529945A (zh) * | 2020-11-17 | 2021-03-19 | 西安电子科技大学 | 一种多视角三维isar散射点集配准方法 |
CN112529945B (zh) * | 2020-11-17 | 2023-02-21 | 西安电子科技大学 | 一种多视角三维isar散射点集配准方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105302874A (zh) | 基于地理云数据的空间匹配方法 | |
CN102867066B (zh) | 数据汇总装置和数据汇总方法 | |
CN105718628B (zh) | 面向公差技术的装配体几何要素误差传递关系图表示和构建方法 | |
CN103778191B (zh) | 一种顾及空间邻近关系的矢量等高线数据划分方法 | |
CN108536923A (zh) | 一种基于建筑cad图的室内拓扑地图生成方法及系统 | |
Huang et al. | Feature matching in cadastral map integration with a case study of Beijing | |
CN108627798A (zh) | 基于线性判别分析和梯度提升树的wlan室内定位算法 | |
CN113656477A (zh) | 一种国土空间多源异构数据的校验与融合方法 | |
Yu et al. | High-precision camera pose estimation and optimization in a large-scene 3D reconstruction system | |
CN112580428A (zh) | 一种配电网设计方法及装置 | |
Guo et al. | Line-based 3d building abstraction and polygonal surface reconstruction from images | |
Xu et al. | Fast and accurate registration of large scene vehicle-borne laser point clouds based on road marking information | |
CN107562693B (zh) | 土地利用/覆被要素多特征变化转移矩阵矢量化提取方法 | |
Wang et al. | Integration of 3DGIS and BIM and its application in visual detection of concealed facilities | |
CN106767405B (zh) | 相位映射辅助三维成像系统快速对应点匹配的方法及装置 | |
Cao et al. | MHSAN: Multi-view hierarchical self-attention network for 3D shape recognition | |
KR101063827B1 (ko) | 한국토지정보시스템 연속지적도와 수치지형도의 기하학적 지도 변환을 위한 반자동화된 공액점 쌍 추출방법 | |
CN102254093A (zh) | 基于泰森多边形的连通域统计相关算法 | |
CN110019632A (zh) | 一种基于形状树的地理要素形状分析方法 | |
Ma et al. | Complex buildings orientation recognition and description based on vector reconstruction | |
CN106652032B (zh) | 一种基于Linux集群平台的DEM并行等高线生成方法 | |
Li et al. | Efficient and precise visual location estimation by effective priority matching-based pose verification in edge-cloud collaborative IoT | |
Li et al. | A cluster-based 3D reconstruction system for large-scale scenes | |
CN105447875A (zh) | 一种电子地形图自动几何校正方法 | |
Dalyot et al. | Integrating network structures of different geometric representations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160203 |
|
RJ01 | Rejection of invention patent application after publication |