CN104809689B - 一种基于轮廓的建筑物点云模型底图配准方法 - Google Patents

一种基于轮廓的建筑物点云模型底图配准方法 Download PDF

Info

Publication number
CN104809689B
CN104809689B CN201510250777.8A CN201510250777A CN104809689B CN 104809689 B CN104809689 B CN 104809689B CN 201510250777 A CN201510250777 A CN 201510250777A CN 104809689 B CN104809689 B CN 104809689B
Authority
CN
China
Prior art keywords
point cloud
cloud model
building
line segment
lines
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
CN201510250777.8A
Other languages
English (en)
Other versions
CN104809689A (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.)
Shenzhen Research Institute Beijing Institute Of Technology
Beijing Institute of Technology BIT
Original Assignee
Shenzhen Research Institute Beijing Institute Of Technology
Beijing Institute of Technology BIT
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 Shenzhen Research Institute Beijing Institute Of Technology, Beijing Institute of Technology BIT filed Critical Shenzhen Research Institute Beijing Institute Of Technology
Priority to CN201510250777.8A priority Critical patent/CN104809689B/zh
Publication of CN104809689A publication Critical patent/CN104809689A/zh
Application granted granted Critical
Publication of CN104809689B publication Critical patent/CN104809689B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Processing Or Creating Images (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及一种基于轮廓的建筑物点云模型底图配准方法，属于计算机视觉技术领域。本发明方法通过对三维重建得到的建筑物点云模型进行结构分析，自动检测当前模型的姿态，最终实现点云模型坐标的归一化，并在此基础上得到当前建筑物的点云俯视轮廓线条图；对重建建筑所在区域的遥感图像进行轮廓分析，自动检测得到底图的轮廓线条；根据重建建筑的轮廓线条估计当前建筑的结构拐点，然后借助轮廓匹配信息实现点云模型与卫星底图间的空间对齐，得到从点云模型到卫星底图间的映射矩阵，自动实现建筑物点云模型在卫星底图上的放置。对比现有技术，本发明方法能够准确对点云模型和底图进行自动化配准，有效降低配准所需人工工作量，降低配准成本。

Description

一种基于轮廓的建筑物点云模型底图配准方法
技术领域：
本发明属于计算机视觉技术领域，涉及一种基于轮廓的建筑物点云模型底图配准方法，具体涉及一种点云模型结构分析、图像轮廓提取及匹配的建筑物点云模型和底图配准方法。
背景技术：
在三维重建中，围绕建筑物拍摄多角度图像序列，能够重建得到较好的点云模型。但是，这些点云模型各自处于独立的坐标系中，不能直接构成建筑群。而独立的建筑物模型场景展示能力较弱，不能有效表达整个建筑群的空间位置关系。因此，需要将建筑群内独立的建筑物点云模型有机融合为一个整体，才能实现整个三维场景的直观展示，体现空间展示的优点。直接从图像序列重建得到的建筑物点云模型仅由一系列具有形状约束的空间离散点的构成，该模型的方向、尺度等几何信息仅在当前局部坐标系下起作用，相对世界坐标系只能保证整体形状一致性。传统的点云模型展示形式为各建筑物模型独立显示，并未将整个大场景联合起来，空间描述能力较弱，不能很好发挥三维展示的特点。一些地理信息系统(GIS)也可将三维模型放置于虚拟地球等地理平台上进行展示，不过放置的过程需要提供大量人工交互以实现空间对准，整个过程自动化程度很低，不利于实际应用。尤其是面对大量模型需要配准的情况，巨大的交互工作量需要消耗大量的时间和精力。
发明内容：
本发明的目的是针对以上问题，提供一种基于轮廓的建筑物点云模型底图配准方法，通过建筑物点云模型结构分析、底图轮廓提取和基于轮廓匹配的建筑物点云模型与底图配准相关技术，实现点云模型与底图之间的空间对齐，最终实现点云模型在底图上的自动放置。
本发明方法的思想是通过对三维重建得到的建筑物点云模型进行结构分析，自动检测当前模型的姿态，最终实现点云模型坐标的归一化，并在此基础上得到当前建筑物的点云俯视轮廓线条图；对重建建筑所在区域的遥感图像进行轮廓分析，自动检测得到底图的轮廓线条；根据重建建筑的轮廓线条估计当前建筑的结构拐点，然后借助轮廓匹配信息实现点云模型与卫星底图间的空间对齐，得到从点云模型到卫星底图间的映射矩阵，自动实现建筑物点云模型在卫星底图上的放置。
本发明的目的是通过以下技术方案实现的：
一种基于轮廓的建筑物点云模型底图配准方法，包括以下步骤：
步骤一、通过以下步骤进行建筑物点云模型俯视轮廓线条提取：
步骤1.1对建筑物点云模型进行姿态估计；
步骤1.2通过以下步骤根据建筑物点云模型z轴方向对其俯视轮廓进行检测：
步骤1.2.1对点云模型在世界坐标系下进行姿态矫正，然后对姿态矫正后的点云模型在XOY平面上投影得到当前建筑物点云模型的俯视图；
步骤1.2.2在俯视图中根据四叉树分割方法检测出模型的轮廓；
步骤1.2.3通过改进的LSD算法对建筑物点云模型俯视轮廓图做边缘检测，得到建筑物的俯视边界线条，即建筑物的俯视轮廓：
步骤二、通过改进的LSD算法对建筑物遥感影像进行基于线特征的轮廓提取；
步骤三、通过以下步骤进行基于轮廓匹配的建筑物点云模型与底图配准：
步骤3.1通过以下步骤筛选经过步骤2.3后得到的所有线段：
步骤3.1.1针对点云模型俯视图lo和卫星底图lb，根据现代化建筑俯视图自身的轮廓特性，将属于当前建筑轮廓所有的线条聚类为两个近似垂直的簇，同时将不属于这两个簇的线段视为噪声线条，予以剔除；
步骤3.1.2求出lo中所有线条的平均长度，如果存在线条长度小于平均长度的倍，则剔除；为预定的点云模型俯视图噪声线条阈值；
步骤3.1.3求出lb中所有线条的平均长度，如果存在线条长度小于平均长度的倍，则剔除；为预定的卫星底图噪声线条阈值；
步骤3.2对经步骤3.1处理后得到的线段通过下述过程找出鲁棒性高的线段交点，即轮廓拐点：
对于检测得到的两个近似垂直的线段簇和取C1中的线段其两个端点分别为和那么所在直线和C2中所有线段所在直线均有交点，对于C2中任意一条线段首先判定和中哪一个点距离所在直线最近；若距离所在直线较近，则将线段和所在直线的交点记为线段在端的交点，反之则记为线段在端的交点；线段在端的交点的能量函数为
和分别在C2中检测得到的最小能量值对应的交点即为当前线段的端点，也是需要保留的备选轮廓拐点；
步骤3.3通过以下步骤找到匹配轮廓拐点，然后计算仿射矩阵：
步骤3.3.1通过对两个图上所有的拐点建立特征点向量并依次进行映射找出所有可能的仿射矩阵，即候选仿射矩阵；
步骤3.3.2通过以下过程找出最佳仿射矩阵：
对于某一候选仿射矩阵H，假设图像lo中一条线段在图像lb中的线段为定义线段和线段的能量函数为
E1、E2、Em和分别为线段和线段的左端点、中点、右端点，d(.)表示两个点之间的距离，对于lo中任一条线段从图像lb中找到与其能量函数最小的即为它的映射线段，线段的能量值定义为
其中，N4为lb中轮廓线的数目；
当前仿射矩阵H的能量值为
其中，N2为lo中轮廓线的数目；
在所有仿射矩阵中，找出能量值最小的仿射矩阵即为最佳仿射矩阵；
步骤3.4通过仿射矩阵将建筑物点云模型映射到卫星底图。
一种基于轮廓的建筑物点云模型底图配准装置，包括轮廓线条提取模块、点云模型与底图配准模块；其中轮廓线条提取模块与点云模型与底图配准模块相连；
轮廓线条提取模块用于提取点云模型和卫星底图中建筑物的轮廓线条；
点云模型与底图配准模块用于根据点云模型与卫星底图中的建筑物轮廓线条将建筑物的点云模型正确放置到卫星底图上；
作为优选，所述轮廓线条提取模块包含了点云模型轮廓线条提取单元和卫星底图轮廓线条提取单元；点云模型轮廓线条提取单元用于提取点云模型中建筑物的轮廓线条，卫星底图轮廓线条提取单元用于提取卫星底图中的建筑物轮廓线条。
有益效果：
对比现有技术，本发明方法能够准确对点云模型和底图进行自动化配准，有效降低配准所需人工工作量，降低配准成本。
附图说明：
图1为本发明实施例一种基于轮廓的建筑物点云模型底图配准方法流程示意图；
图2为本发明实施例线段之间的纵向距离概念示意图；
图3为本发明实施例一种基于轮廓的建筑物点云模型底图配准装置结构示意图。
具体实施方式
下面将结合附图和实施例对本发明加以详细说明，同时也叙述了本发明技术方案解决的技术问题及有益效果，需要指出的是，所描述的实施例仅旨在便于对本发明的理解，而对其不起任何限定作用。
实施例1
一种基于轮廓的建筑物点云模型底图配准方法，包括以下步骤：
步骤一、建筑物点云模型边界线条提取
步骤1.1建筑物点云模型姿态估计
在本实施例中利用开源点云数据处理库PCL(Point Cloud Library)对点云模型中点的法线进行估计，具体步骤如下：
a)使用pcl的StatisticalOutlierRemoval滤波器移除偏远点；
b)使用pcl的NormalEstimation方法通过点附近的点云信息拟合切面的方法估计该点的法向量；
然后对所有点的法线以角度θ进行聚类，因为常规的现代化建筑物的各相邻外表面自身就存在互相垂直的空间结构关系，所以存在三个分类的向量在相互近似垂直的同时它们所在聚类的三维点数目应该是较多的，所以先对所有聚类按照三维点个数从大到小排列，然后从前往后依次遍历所有可能的组合情况，直到找到第一组三个互相近似垂直的向量d1,d2,d3，也即三个向量的夹角都在80度到100度之间变动，这时候将向量d1,d2,d3分别视为建筑物点云模型的垂直正方向和建筑物墙面的2个法线向量，然后从d1,d2,d3中选取与模型所在世界坐标系坐标轴z轴夹角最小的向量作为点云模型的正方向。
作为优选，θ＝10；
步骤1.2建筑物点云模型俯视轮廓检测
步骤1.2.1得到点云模型在世界坐标系下的俯视图
在检测到点云模型的z轴方向后，需要将点云模型归一化到统一的世界坐标系下，以便得到整个建筑物点云模型的俯视图，进而检测得到该建筑的俯视轮廓线条，具体做法如下:
利用上一步骤中得到的三个向量d1,d2,d3与世界坐标系中3个正方向进行比较，求出检测到的建筑物点云所在的坐标系与世界坐标系间的旋转角度分别为α,β,γ，那么点云模型各点(x,y,z)矫正后的坐标为
矫正建筑物点云模型的坐标系后，将新坐标系下的点云模型在XOY平面上投影便得到当前建筑物点云模型的俯视图；
步骤1.2.2从俯视图中检测出模型的轮廓
本方法采用一种自顶向下的分割方法对模型俯视投影区域进行四叉树分割，具体如下：
首先提取模型所在的最外层矩形区域(即当前投影模型的包络)，并针对这个矩形区域进行四叉树分割条件判断；如果满足分割条件，则将该矩形区域分割为4个子区域，然后分别对四个子区域进行下一步分割条件判断和分割处理，如此迭代循环，直至不能再进行细分截止；
所述四叉树分割判断条件主要包含以下几个因素：当前区域内点云密度，当区域内点云密度ρ小于预设阈值ρmin时，该区域便不再进行分割；投影到当前区域内所有点的坐标z值的标准差SD，当标准差小于预设的阈值SDmin时，该区域便不再进行分割；
作为优选，ρmin为整体点云密度的0.35倍，SDmin为坐标z值最大值与最小值之和的0.45倍。
在得到四叉树分割的结果后，针对分割后的区域进行初步分类，将整个点云模型俯视图分割为以下四种区域：无点区域、墙面点区域、屋顶点区域和地面点区域；分类时主要考虑到的因素是各分区区域内点的密度信息和z值分布信息；例如墙面点区域在整个z轴分布密度较高，屋顶点区域仅在z值较大的某个平面分布较密，地面点在z值较小的平面分布比较密集等；
接下来，对于初步分类后的墙面点区域根据相邻关系将墙面进行拼接融合，得到相对准确的最表层三维点区域；然后根据这些三维点的俯视投影，初步勾勒出整栋建筑的俯视轮廓图。
步骤1.3通过改进的LSD算法对建筑物点云模型俯视轮廓图做边缘检测，得到建筑物的俯视边界线条，即建筑物的俯视轮廓：
由于图像数据存在诸如形变、模糊或遮挡等质量问题，导致LSD算法的提取结果往往存在大量的断裂和噪声，无法较为完整地勾勒出建筑物的完整线特征信息；因此，本实施例对传统的LSD算法进行了改进，通过以下过程获取建筑物的俯视边界线条：
a)使用传统LSD算法提取图像中的线段
b)对于上一步骤得到的线段进行以下五种判定：
i.噪声剔除
根据预设的阈值THRl，对步骤2.1得到的所有线段进行比较，删除其
中长度小于THRl的线段；
作为优选，THRl＝10；
ii.斜率一致性判定
当两条线段l1和l2之间的夹角θ小于预设的夹角阈值THRθ时，便认为这
两条线段是平行的；
作为优选，THRθ＝15；
iii.平行线过滤
对于被判定为平行的线段对l1和l2，若两线段间距离大于一定阈值且有足够大的平行部分，则认为这两条线段平行且足够分离，不进行合并，这样可以保证一些平行的建筑线特征不会被合并在一起，具体可通过如下过程：
假设两个线段的长度分别为L1和L2，a为两条线段平行且投影相互重叠的部分长度，设定比例阈值THRp，当小于THRp时，认为两条线段是同一条长线段断裂成的子线段，需要将线段l1和l2进行合并，反之则表示这两条线段为两条独立且平行的线条，直接保留；
作为优选，THRp＝0.4；
iv.横向距离判定
当两条线段l1和l2之间互相投影的距离足够近时，可以判定它们处于同一条直线上，Li为线段l1上各采样点到l2所在直线间的距离，Gj为线段l2上各采样点到l1所在直线之间的距离。可以得到以下公式
式中，N1和N2分别为线段l1和l2上点的数目。设定横向距离阈值THRc，当两条线段投影距离满足max(dL,dG)＜THRc时，则判定两条线段在同一条直线上；
作为优选，THRc＝4.5；
v.纵向距离判定
当检测到两条线段l1和l2处于同一条直线上时，这两条线段可能出现以下两种情况：两条线段原本就是同一条长线段断裂成为的子线段，这时候线段l1和l2的纵向距离较近；线段l1和l2不属于同一条轮廓线，仅仅因为两条线段所处的局部区域线条处于同一直线上，此时两条线段间纵向距离一般较大，极端情况为两条线段分别处于图像的两端，具体可通过下述过程完成：
L1和L2分别为两条线段l1和l2的长度，l1和l2间的纵向距离为dmin，其具体含义请参考图1，定义两条线段同属一条线段的隶属度为设定一个隶属度阈值THRz，当满足时表示这两条线段是同一条线段断裂成的两个子线段，需要进行拼接，反之则表示这个线段对只是两条距离较远但近似共线的独立线段；
作为优选，THRz＝0.1；
c)利用上述判定方式拟合新的线条
通过上述5种判定，可以在利用LSD算法检测得到初始线段后，利用线条的长度和几何分布等特征分别做噪声线段剔除和线段可拼接性检测，然后对可拼接的线段利用最小二乘算法做融合处理：对于检测到的可拼接线段，直接利用两条线段中包含的所有像素直接拟合出一条新的线段，并在拟合得到新的线条后重复步骤b)和步骤c)再次进行判定和拟合，直至不能再检测出可拼接线段时结束；
步骤二、基于线特征的建筑物遥感影像轮廓提取
使用上述改进的LSD算法得到遥感建筑物图像的轮廓；
步骤三、基于轮廓匹配的建筑物点云模型与底图配准
因为点云模型的建筑物轮廓与底图中建筑物轮廓存在旋转、平移、缩放关系，所以需要得到两个图的映射关系，而通过线条是无法确定两张图的仿射矩阵的，所以需要在提取出一系列准确匹配的特征对。针对初步提取到的建筑物轮廓线条，在经过一定处理可以得到更加鲁棒的建筑物轮廓线条，然后通过搜索一系列鲁棒性高的线段交点，即本文所需的建筑物轮廓结构拐点，就可以利用它们进行建筑物匹配；
步骤3.1筛选上一步骤中得到的所有线段
步骤3.1.1针对点云模型俯视图lo和卫星底图lb，根据现代化建筑俯视图自身的轮廓特性，将属于当前建筑轮廓所有的线条聚类为两个近似垂直的簇，也即两个簇各自的平均线条方向夹角在80到100度之间，同时将不属于这两个簇的线段视为噪声线条，予以剔除；
步骤3.1.2求出lo中所有线条的平均长度，如果存在线条长度小于平均长度的倍，则剔除；
作为优选，
步骤3.1.3求出lb中所有线条的平均长度，如果存在线条长度小于平均长度的倍，则剔除；
作为优选，
步骤3.2通过上一步骤得到的线段找出鲁棒性高的线段交点，即轮廓拐点
对于检测得到的两个近似垂直的线段簇和取C1中的线段其两个端点分别为和那么所在直线和C2中所有线段所在直线均有交点，对于C2中任意一条线段首先判定和中哪一个点距离所在直线最近。若距离所在直线较近，则将线段和所在直线的交点记为线段在端的交点，反之则记为线段在端的交点。设线段在端的交点的能量函数为
那么和分别在C2中检测得到的最小能量值对应的交点即为当前线段的端点，也是需要保留的备选轮廓拐点；
步骤3.3找到匹配轮廓拐点计算仿射矩阵
步骤3.3.1找出所有可能的仿射矩阵
假设点云模型俯视图lo最终检测得到的轮廓拐点为生成的轮廓线为卫星底图lb最终检测到的轮廓拐点为生成的轮廓线为对于模型俯视图中的任意一个拐点得到与它相接的两个线段的另外两个端点这三个点组成一个特征点向量在卫星底图上任一拐点同样找到与它相接的两个线段的另外两个端点，组成两个特征点向量 然后通过特征向量分别与的对应关系就可以得到两个仿射矩阵，遍历两个图上所有的拐点可以得到所有可能的仿射矩阵；
步骤3.3.2找出最佳仿射矩阵
对于某一候选仿射矩阵H，假设图像lo中一条线段在图像lb中的线段为定义线段和线段的能量函数为
E1、E2、Em和分别为线段和线段的左端点、中点、右端点，d(.)表示两个点之间的距离，本例中使用欧氏距离，对于lo中任一条线段从图像lb中找到与其能量函数最小的即为它的映射线段，线段的能量值定义为
当前仿射矩阵H的能量值为
在所有仿射矩阵中，找出能量值最小的仿射矩阵即为所求，然后通过仿射矩阵就可以把点云模型映射到卫星底图。
实施例2
如图3所示为一种基于轮廓的建筑物点云模型底图配准装置结构示意图，由图中可以看出该装置包括轮廓线条提取模块、点云模型与底图配准模块；
轮廓线条提取模块用于提取点云模型和卫星底图中建筑物的轮廓线条；
点云模型与底图配准模块用于根据点云模型与卫星底图中的建筑物轮廓线条将建筑物的点云模型正确放置到卫星底图上；
作为优选，所述轮廓线条提取模块包含了点云模型轮廓线条提取单元和卫星底图轮廓线条提取单元；点云模型轮廓线条提取单元用于提取点云模型中建筑物的轮廓线条，卫星底图轮廓线条提取单元用于提取卫星底图中的建筑物轮廓线条。
以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于轮廓的建筑物点云模型底图配准方法，其特征在于，包括以下步骤：
步骤一、通过以下步骤进行建筑物点云模型俯视轮廓线条提取：
步骤1.1对建筑物点云模型进行姿态估计；
步骤1.2通过以下步骤根据建筑物点云模型z轴方向对其俯视轮廓进行提取：
步骤1.2.1对点云模型在世界坐标系下进行姿态矫正，然后对姿态矫正后的点云模型在XOY平面上投影得到当前建筑物点云模型的俯视图；
步骤1.2.2在俯视图中根据四叉树分割方法检测出模型的轮廓；
步骤1.2.3通过如下改进的LSD算法对建筑物点云模型俯视轮廓图做边缘检测，得到建筑物的俯视边界线条，即建筑物的俯视轮廓：
步骤1.2.3.1使用传统LSD算法提取图像中的线段；
步骤1.2.3.2对于上一步骤得到的线段进行以下五种判定：
(1)噪声剔除：根据预设的阈值THRl，对所有线段进行比较，删除其中长度小于THRl的线段；
(2)斜率一致性判定：当两条线段l1和l2之间的夹角θ小于预设的夹角阈值THRθ时，便认为这两条线段是平行的；
(3)平行线过滤：对于被判定为平行的线段对l1和l2，当小于THRp时，认为两条线段是同一条长线段断裂成的子线段，需要将线段l1和l2进行合并，反之则表示这两条线段为两条独立且平行的线条，直接保留；其中，l1和l2的长度分别为L1和L2，a为两条线段平行且投影相互重叠的部分长度，THRp为预设的比例阈值；
(4)横向距离判定：当两条线段l1和l2之间互相投影的距离足够近时，可以判定它们处于同一条直线上，Li为线段l1上各采样点到l2所在直线间的距离，Gj为线段l2上各采样点到l1所在直线之间的距离，根据下述公式：
其中，N1和N2分别为线段l1和l2上采样点的数目；以及预设的横向距离阈值THRc，当两条线段投影距离满足max(dL,dG)＜THRc时，则判定两条线段在同一条直线上；
(5)纵向距离判定：当检测到两条线段l1和l2处于同一条直线上时，根据下述两条线段同属一条线段的隶属度计算公式以及预设的隶属度阈值THRz判断，当时表示这两条线段是同一条线段断裂成的两个子线段，需要进行拼接；反之则表示这个线段对只是两条距离较远但近似共线的独立线段：
其中，dmin为l1和l2间的纵向距离；
步骤1.2.3.3利用上述判定拟合新的线条；
步骤1.2.3.4对经过步骤1.2.3.3得到的线条重复步骤1.2.3.2和步骤1.2.3.3过程，直到不能再检测出可拼接的线段时为止；
步骤二、通过改进的LSD算法对建筑物遥感影像进行基于线特征的轮廓提取；
步骤三、通过以下步骤进行基于轮廓匹配的建筑物点云模型与底图配准：
步骤3.1通过以下步骤筛选经过前面步骤后得到的所有线段：
步骤3.1.1针对点云模型俯视图lo和卫星底图lb，根据现代化建筑俯视图自身的轮廓特性，将属于当前建筑轮廓所有的线条聚类为两个近似垂直的簇，同时将不属于这两个簇的线段视为噪声线条，予以剔除；
步骤3.1.2求出lo中所有线条的平均长度，如果存在线条长度小于平均长度的倍，则剔除；为预定的点云模型俯视图噪声线条阈值；
步骤3.1.3求出lb中所有线条的平均长度，如果存在线条长度小于平均长度的倍，则剔除；为预定的卫星底图噪声线条阈值；
步骤3.2对经步骤3.1处理后得到的线段通过下述过程找出鲁棒性高的线段交点：
对于检测得到的两个近似垂直的线段簇和取C1中的线段其两个端点分别为和那么所在直线和C2中所有线段所在直线均有交点，对于C2中任意一条线段首先判定和中哪一个点距离所在直线最近；若距离所在直线较近，则将线段和所在直线的交点记为线段在端的交点，反之则记为线段在端的交点；线段在端的交点的能量函数为
和分别在C2中检测得到的最小能量值对应的交点即为当前线段的端点，也是需要保留的备选轮廓拐点；
步骤3.3通过以下步骤找到匹配轮廓拐点以及最佳仿射矩阵：
步骤3.3.1通过对两个图上所有的备选轮廓拐点建立特征点向量并依次进行映射找出所有可能的仿射矩阵，即候选仿射矩阵；
步骤3.3.2通过以下过程找出最佳仿射矩阵：
对于某一候选仿射矩阵H，假设图像lo中一条线段在图像lb中的线段为定义线段和线段的能量函数为
E1、E2、Em和分别为线段和线段的左端点、中点、右端点，d(.)表示两个点之间的距离，对于lo中任一条线段从图像lb中找到与其能量函数最小的即为它的映射线段，线段的能量值定义为
其中，N4为lb中轮廓线的数目；
当前仿射矩阵H的能量值为
其中，N2为lo中轮廓线的数目；
在所有仿射矩阵中，能量值最小的仿射矩阵即为最佳仿射矩阵；
步骤3.4通过最佳仿射矩阵将建筑物点云模型映射到卫星底图。
2.根据权利要求1所述的一种基于轮廓的建筑物点云模型底图配准方法，其特征在于：所述步骤1.1进一步通过以下步骤完成姿态估计：
步骤1.1.1采用开源点云数据处理库PCL对点云模型中点的法线进行估计，具体步骤如下：
a)使用pcl的StatisticalOutlierRemoval滤波器移除偏远点；
b)使用pcl的NormalEstimation方法通过点附近的点云信息拟合切面的方法估计该点的法向量；
步骤1.1.2对所有点的法线以角度θ进行聚类，然后对所有聚类按照三维点个数从大到小排列，接下来从前往后依次遍历所有可能的组合情况，直到找到第一组三个互相近似垂直的向量d1,d2,d3，也即三个向量的夹角都在80度到100度之间变动，这时候将向量d1,d2,d3分别视为建筑物点云模型的垂直正方向和建筑物墙面的2个法线向量；
步骤1.1.3从d1,d2,d3中选取与模型所在世界坐标系坐标轴z轴夹角最小的向量作为点云模型的正方向。
3.根据权利要求2所述的一种基于轮廓的建筑物点云模型底图配准方法，其特征在于：θ＝10。
4.根据权利要求1所述的一种基于轮廓的建筑物点云模型底图配准方法，其特征在于：所述步骤1.2.2进一步通过以下过程完成：
首先，提取模型所在的最外层矩形区域，并针对这个矩形区域进行四叉树分割条件判断；如果满足分割条件，则将该矩形区域分割为4个子区域，然后分别对四个子区域进行下一步分割条件判断和分割处理，如此迭代循环，直至不能再进行细分截止；
其次，对分割后的区域根据各分区区域内点的密度信息和z值分布信息进行初步分类，将整个点云模型俯视图分割为以下四种区域：无点区域、墙面点区域、屋顶点区域和地面点区域；
再次，对于初步分类后的墙面点区域根据相邻关系将墙面进行拼接融合，得到相对准确的最表层三维点区域；然后根据这些三维点的俯视投影，初步勾勒出整栋建筑的俯视轮廓图；
最后，通过改进的LSD算法对建筑物点云模型俯视轮廓图做边缘检测，得到建筑物的俯视边界线条。
5.根据权利要求4所述的一种基于轮廓的建筑物点云模型底图配准方法，其特征在于：所述四叉树分割判断条件包含以下几个因素：当前区域内点云密度，当区域内点云密度ρ小于预设阈值ρmin时，该区域便不再进行分割；投影到当前区域内所有点的坐标z值的标准差SD，当标准差小于预设的阈值SDmin时，该区域便不再进行分割。
6.根据权利要求5所述的一种基于轮廓的建筑物点云模型底图配准方法，其特征在于：ρmin为整体点云密度的0.35倍，SDmin为坐标z值最大值与最小值之和的0.45倍。
7.根据权利要求1所述的一种基于轮廓的建筑物点云模型底图配准方法，其特征在于：所述THRl＝10，THRθ＝15，THRp＝0.4，THRc＝4.5，THRz＝0.1。
8.根据权利要求1所述的一种基于轮廓的建筑物点云模型底图配准方法，其特征在于：所述
CN201510250777.8A 2015-05-15 2015-05-15 一种基于轮廓的建筑物点云模型底图配准方法 Expired - Fee Related CN104809689B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510250777.8A CN104809689B (zh) 2015-05-15 2015-05-15 一种基于轮廓的建筑物点云模型底图配准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510250777.8A CN104809689B (zh) 2015-05-15 2015-05-15 一种基于轮廓的建筑物点云模型底图配准方法

Publications (2)

Publication Number Publication Date
CN104809689A CN104809689A (zh) 2015-07-29
CN104809689B true CN104809689B (zh) 2018-03-30

Family

ID=53694494

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510250777.8A Expired - Fee Related CN104809689B (zh) 2015-05-15 2015-05-15 一种基于轮廓的建筑物点云模型底图配准方法

Country Status (1)

Country Link
CN (1) CN104809689B (zh)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107154016B (zh) * 2016-03-01 2019-02-26 腾讯科技（深圳）有限公司 立体图像中目标对象的拼接方法和装置
CN106384192A (zh) * 2016-08-31 2017-02-08 南京品淳通信科技有限公司 一种基于gis的3d指挥调度系统以及实现方法
CN107144265B (zh) * 2017-04-21 2018-04-10 深圳市数字城市工程研究中心 一种基于投影图解的三维产权体空间界址点的测量方法
CN107507167B (zh) * 2017-07-25 2020-03-10 上海交通大学 一种基于点云平面轮廓匹配的货物托盘检测方法及系统
CN108022307B (zh) * 2017-11-26 2021-06-18 中国人民解放军陆军装甲兵学院 基于增材再制造点云模型的自适应平面分层方法
WO2019104666A1 (zh) * 2017-11-30 2019-06-06 深圳市沃特沃德股份有限公司 扫地机器人及其进行区域划分的方法
CN108920765B (zh) * 2018-06-06 2021-11-30 东南大学 一种基于建筑物三维线段模型的假设平面拟合方法
CN109271654B (zh) * 2018-07-19 2023-04-18 平安科技（深圳）有限公司 模型轮廓的切分方法及装置、存储介质、终端
CN109187548B (zh) * 2018-07-24 2020-04-07 西南交通大学 一种岩体裂隙识别方法
CN109215064B (zh) * 2018-08-03 2021-06-08 华南理工大学 一种基于超像素向导的医学图像配准方法
CN109166149B (zh) * 2018-08-13 2021-04-02 武汉大学 一种融合双目相机与imu的定位与三维线框结构重建方法与系统
CN109344435B (zh) * 2018-08-28 2022-12-06 贵州数智联云工程科技有限公司 一种用于建筑排布的设计方法及其装置、存储介质
CN109934837B (zh) * 2018-12-26 2023-04-28 江苏名通信息科技有限公司 一种3d植物叶片轮廓的提取方法、装置及系统
CN111862179B (zh) * 2019-04-12 2022-02-18 北京城市网邻信息技术有限公司 三维对象建模方法与设备、图像处理装置及介质
CN110188616B (zh) * 2019-05-05 2023-02-28 上海盎维信息技术有限公司 基于2d和3d影像的空间建模方法及装置
CN112348887A (zh) * 2019-08-09 2021-02-09 华为技术有限公司 终端位姿确定方法以及相关装置
CN113739797A (zh) * 2020-05-31 2021-12-03 华为技术有限公司 视觉定位方法和装置
CN112001947A (zh) * 2020-07-30 2020-11-27 海尔优家智能科技（北京）有限公司 拍摄位置的确定方法及装置、存储介质、电子装置
CN112069585B (zh) * 2020-09-16 2021-06-22 中国测绘科学研究院 一种分区降维化处理的建筑遮挡高效识别方法
CN112346409B (zh) * 2020-11-17 2021-10-08 浙江大学 一种快速判断点云包含问题的方法
CN112465849B (zh) * 2020-11-27 2022-02-15 武汉大学 一种无人机激光点云与序列影像的配准方法
CN112288820B (zh) * 2020-12-29 2021-04-27 深圳市六合智能感知系统科技有限公司 规则物体的轮廓计算方法、装置、设备及存储介质
CN113298838B (zh) * 2021-07-27 2021-09-21 成都睿铂科技有限责任公司 一种物体轮廓线提取方法及系统
CN113724270B (zh) * 2021-08-25 2023-08-08 华南农业大学 一种牲畜表面点云智能分割方法及其系统
CN113469887B (zh) * 2021-09-02 2022-01-07 深圳市信润富联数字科技有限公司 物体数模的转换方法、装置、设备及存储介质
CN113971307A (zh) * 2021-10-27 2022-01-25 深圳须弥云图空间科技有限公司 一种关联关系生成方法、装置、存储介质及电子设备
CN114461104B (zh) * 2022-02-07 2024-04-16 深圳须弥云图空间科技有限公司 楼型拼合方法、装置、设备及存储介质
CN114429588B (zh) * 2022-04-06 2022-06-03 四川省冶勘设计集团有限公司 应用于矿山环境恢复治理的数据处理方法及系统
CN115082951B (zh) * 2022-07-19 2023-01-10 万翼科技有限公司 梁构件识别方法、装置、设备及存储介质
CN117808981A (zh) * 2023-12-25 2024-04-02 北京市市政工程设计研究总院有限公司 数字高程模型构建方法
CN117745955B (zh) * 2024-02-20 2024-05-07 北京飞渡科技股份有限公司 基于建筑基底矢量数据生成城市建筑场景的方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930540A (zh) * 2012-10-26 2013-02-13 中国地质大学（武汉） 城市建筑物轮廓检测的方法及系统
CN104484668A (zh) * 2015-01-19 2015-04-01 武汉大学 一种无人机多重叠遥感影像的建筑物轮廓线提取方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4115117B2 (ja) * 2001-10-31 2008-07-09 キヤノン株式会社 情報処理装置および方法
US20090232355A1 (en) * 2008-03-12 2009-09-17 Harris Corporation Registration of 3d point cloud data using eigenanalysis

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930540A (zh) * 2012-10-26 2013-02-13 中国地质大学（武汉） 城市建筑物轮廓检测的方法及系统
CN104484668A (zh) * 2015-01-19 2015-04-01 武汉大学 一种无人机多重叠遥感影像的建筑物轮廓线提取方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"A Comprehensive Automated 3D Approach for Building Extraction,Reconstruction,and Regularization from Airborne Laser Scanning Point Clouds";Peter Dorninger等;《sensors》;20081117;第8卷(第11期);第7323-7343页 *
"Fully Automatic Registration of 3D Point Clouds";Ameesh Makadia等;《CVPR ’06 Proceedings of the 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition》;20060622;第1卷;第1297-1304页 *
"从机载激光雷达点云提取建筑物屋顶边界的活动轮廓模型改进方法";孙颖等;《测绘学报》;20140630;第43卷(第6期);第620-626,636页 *
"点、线相似不变性的城区航空影像与机载激光雷达点云自动配准";张良等;《测绘学报》;20140430;第43卷(第4期);第372-379页 *

Also Published As

Publication number Publication date
CN104809689A (zh) 2015-07-29

Similar Documents

Publication Publication Date Title
CN104809689B (zh) 一种基于轮廓的建筑物点云模型底图配准方法
CN111815757B (zh) 基于图像序列的大型构件三维重建方法
CN110570428B (zh) 一种从大规模影像密集匹配点云分割建筑物屋顶面片的方法及系统
Lafarge et al. Creating large-scale city models from 3D-point clouds: a robust approach with hybrid representation
Zhao et al. Rectilinear parsing of architecture in urban environment
Cui et al. Automatic 3-D reconstruction of indoor environment with mobile laser scanning point clouds
CN109598290A (zh) 一种基于两级检测相结合的图像小目标检测方法
CN102804231B (zh) 三维场景的分段平面重建
CN107564062A (zh) 位姿异常检测方法及装置
CN105957059B (zh) 电子元件漏件检测方法和系统
CN110310320A (zh) 一种双目视觉匹配代价聚合优化方法
Wang et al. A robust approach for automatic registration of aerial images with untextured aerial lidar data
Li et al. Line segment matching and reconstruction via exploiting coplanar cues
CN108090435A (zh) 一种可停车区域识别方法、系统及介质
CN112163622B (zh) 全局与局部融合约束的航空宽基线立体像对线段匹配方法
Lafarge et al. Insertion of 3-D-primitives in mesh-based representations: Towards compact models preserving the details
CN103247045A (zh) 一种从多视图中得到人造场景主方向及图像边缘的方法
CN110400322A (zh) 基于颜色和三维几何信息的水果点云分割方法
CN107146239A (zh) 卫星视频运动目标检测方法及系统
CN104657986B (zh) 一种基于子空间融合和一致性约束的准稠密匹配扩展方法
CN108550166A (zh) 一种空间目标图像匹配方法
CN105335688B (zh) 一种基于视觉图像的飞机机型识别方法
CN110222638B (zh) 一种基于素描图候选框策略和Fast R-CNN的飞机目标检测方法
CN109187548A (zh) 一种岩体裂隙识别方法
Bartelsen et al. Orientation and dense reconstruction from unordered wide baseline image sets

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180330

Termination date: 20190515