测绘人招聘、求职,首选“测绘人才网”
目前,轻小型航空遥感(无人飞行器低空航测系统)是不可缺少的补充技术手段。传统的UAV上装载单相机进行航测技术获取的数字地表模型(DSM),通常为2.5维,尤其在城市地区,很难区分地物点和地面点,获取DEM较为困难。但是UAV装载组合宽角相机拍摄的倾斜影像,利用4个不同视图方向DSM来弥补垂直DSM在立面及局部地面、建筑物屋顶面空间信息的不足,来获取全方位的DSM很容易检测和剔除建筑物等高出地面的地物,获得真实地形DEM。 倾斜影像获取的DSM具有几何约束特点,影像密集匹配生成全方位点云DSM。利用已知良好匹配点建立不规则三角网格(TIN),以实现纹理自适应的匹配传播,并利用新匹配点对TIN进行动态细化,有利于减少因纹理贫乏因素而导致的错误匹配,提高DSM的精度。这种点、三角网及纹理几何约束的关系,决定了DSM具有一定几何约束特点。 滤波是DSM获取高精度DEM的关键技术,必须从DSM检测和剔除建筑物等高出地面的地物,获得真实地形DEM现有滤波的方法很多,由于数据源的不同,滤波方法也大不相同。根据DSM的三角网格具有几何约束特点,需要在基于TIN模型上进行滤波处理。 以倾斜影像获取的点云DSM为主要研究对象,本文在基于TIN的分割方法的基础之上,提出基于法向量差值区域生长分割TIN的滤波算法。该算法首先以规则格网将TIN进行分块,进而在每个区域内将若干个地面点和地物点作为种子,用区域增长的方法进行生长,计算三角形的法向量,认为法向量变化到一定阈值时进行分割;其次利用分割段高程、法向量作为约束条件,利用分割后的点云实现了建筑物、植被和其他地物的识别分类,并在获得地面点后将其插值生成数字高程模型。
基于区域生长分割TIN的滤波算法
1 等格网DSM分块 DSM是点云构成的不规则三角网格,数据量十分庞大。通常在滤波之前,考虑将复杂的整体DSM进行分块,采用“先分块后融合”的思路构建子网格,最后融合子网格恢复原始实体的表面形状。这样就能把复杂的问题简单化,提高了提取DEM的效率。 等格网DSM分块的基本思路是:设定固定格网大小,对DSM进行规则格网划分,根据格网单元的范围确定每个区域中的三维点,即构建数据索引,从而实现对DSM数据的分块。如图1所示,实现了基于邻域关系的快速查找,在各自网格内进行区域生长,从而提高滤波算法的效率,便于区域生长分割TIN。
2 基于区域生长分割TIN方法 本文采用基于区域生长分割TIN方法。区域生长的基本思想是:首先将三角网内的三角形进行拓扑构建。采用自动选取种子的方式,在三角网平坦的地方,依据三角面片的几何原则,选择种子点;以一个种子点开始,根据生长原则选出可以进行生长的三角面片进行区域生长,把未在任何类中的三角面片元素放入种子类中,反复运行直到种子队列为空,停止分割,实现了三角网分割成段。 图2为区域生长过程示意图,图中颜色部分与非颜色部分分别代表区域生长过程中的已生长区域与未生长区域。该算法的流程如图3所示。
2.1 算法原理 (1) 种子选取:种子点的选择是进行区域生长的关键前提,对最后的分割结果影响很大。如果种子点选择不当,将导致整个分割过程的判断依据失真,最终得出错误的分割结果。因此,选择种子的几何原则如下: 几何原则是在大量三角片面中选择种子,考虑的是面片平坦的地方。D(Tt)可以用来控制一个三角形的法向量和它所有的相关三角形的法向量的差值距其期望值的距离;A_I(Tt)可以用来控制一个三角形的法向量和它所有相关的三角形的法向量的差值的变化幅度,即 A_I(Tt)<> 三角形Tt和其所有相关三角形中的一个三角形Ti法向量的差值为
三角形Tt和其所有相关三角形法向量的期望为E(Tt),方差为D(Tt)
保存一个三角形和其所有相关三角形法向量的最大值与最小值的差值记为
(2) 区域生长原理:DSM一般是由地面、建筑物、植被及其他地物组成。因此,TIN模型中只要存在两个面法向量差值较大的明显不是一类,根据DSM特点,确定相邻三角面片法向量差值作为区域生长的驱动信号。 三角形法向量Ni和Nj的差值数学描述为dis=|Ni-Nj|,其中Ni和Nj已经归一化,取值范围为[0,2],平行时为0和2,垂直时为1.414,可知相邻两个法向量的阈值为a、b (a 2.2 算法流程 (1) 构建点云数据三角网格的拓扑关系,使得点、线及三角形有一定的拓扑关系。 (2) 种子的选择采用候选种子的方法,考虑面平坦的地方,根据一定的几何原则,将符合要求的三角形压入一个候选种子队列中;然后在候选种子队列中从队列头开始查找,如果已经在类中则弹出,直到找到一个没有在任何类中的三角形,把它选为种子,把种子初始化为一个类。 (3) 从类的边界三角形队列中,根据生长原则,开始生长。 (4) 同样根据生长原则选出类中边界三角形的相关三角形进行生长,在生长的同时,记录下生长中类的新边界,并将它们压入到一个暂时队列中。当类边界三角形队列为空时,把类中边界三角形队列更新。 (5) 一个类边界三角形的队列为空时,如果类中三角形总体的数目同生长前没有发生变化,类就停止生长。 (6) 如果没有达到停止条件,转到步骤(2)。如果达到了停止的条件,生成了一个类,并把它压入存放类的队列中,并转到步骤(1)。 (7) 反复顺序执行步骤(1)—步骤(5)直到候选种子队列为空。把没有在任何类中的三角形按照期望从小到大的顺序压入种子候选队列二中,再执行步骤(1)—步骤(5),直到候选种子队列二的队列为空,结束分割。 3 基于聚类分析的DSM分类 由于基于法向量差值分割TIN具有局限性,还是可能会出现同一特征被划分为不同区域的过分割现象。在区域生长分割后,会由于过分割现象、模型误差而产生小区域,因此,根据相邻分割段法向量的朝向小于一定的阈值,进行合并处理。通过聚类分析的方法,将合并后的分割段,根据法向量、高度的约束条件,将TIN分类成地面、建筑物、植被及其他地物,方便以后的过滤。 3.1 建筑物判别 由于地面部分相对平坦,几乎没有突变,绝大部分建筑物的顶部是一个平面。根据分割段与段之间的几何拓扑,首先按照相邻分割段的平均法向量的夹角φ区分建筑物的墙面。类的平均法向量为P1(A1, B1, C1),相邻类的平均法向量为P2(A2, B2, C2),则两个法向量之间的夹角为
即φ>φ1或φ<> 屋顶和地面均相对平坦,很难区分。因此,需要根据墙面特征来判断屋顶,判断的条件为:参照图6为例,以一个墙面开始,找到该墙面的相邻分割段,从中找出高出该墙面平均高度的分割段,则该分割段为屋顶,低于墙面平均高度的分割段为地面点。 3.2 植被判别 植被通常表现为相邻点之间空间分布离散,且有较大的高差,因此植被通常会被分为破碎的分割段。这些分割段面积小、数量多、高差大,检测到这种分割段则认为是植被。 对于有一定坡度地面上的植被来说,由于通常滤波方法都是基于地物点明显高于地面点实现的,因此这种假设对于地形起伏较大时,地形点可能和植被点具有相同的高度从而造成滤波失败。当然,低矮的植被是低矮地物的一种,与低矮地物一样,仅仅依靠几何信息时,往往难以滤除。对于这种难以分辨植被的情况,采取人工选择地面种子点,进行区域生长,获取地面点。 3.3 其他地物判别 按照尺寸和几何特征检测出建筑物和植被,剩下的基本是地面点,但是也包含一些附在地面的低矮地物,如小汽车、人等。其他地物分割段之间的所有可能总结如下:①待定分割段比它周围的分割段都要高;②待定分割段面积小。 经过上述分类,将TIN分为地面点和地物点,将获取的地面点直接进行三角剖分,得到高精度DEM。
试验及分析
1 滤波试验 本文采用C++来实现基于分割TIN的滤波算法,获取高精度DEM,并用若干个复杂城市地形的倾斜影像测试算法。下面是其中一个试验区的试验结果及分析。 该试验区地形复杂,右下角是一个明显高低起伏的山地,整个区域覆盖了不同高度的植被,一个明显的居民区位于左部分,居民区内的建筑物大小不一。 试验数据为吉林省敦化市的倾斜影像匹配的DSM数据,大小为1000m*600m。点云数据如图4所示,TIN数据如图5所示。按照本文的算法,在全区范围内选取若干个种子点,取相邻三角形法向量差值进行区域生长分割三角网,根据分割段的平均法向量和平均高度将三角网分类成地面、建筑物和其他地物,图6为分类的晕渲图;根据分类的结果,过滤地物点,获取地面点,图7为全区滤波后的晕渲图,该晕渲图的颜色与图6一致,即两幅图中相同的高程具有相同的颜色;图7将获取的地面点生成DEM,图8为内插DEM结果晕渲图。
2 精度评定 DEM作为基础地理信息数据,其精度至关重要。DEM提取之后,需要对其质量进行精度评定,本试验采用检查点方法对DEM进行精度评定。 精度评定的目的是检验DEM的高程精度,常规检查点的方法是一个检查点对应图上一个点,这种一一对应的方法适用于数字测图的评定;本文是检验DEM的高程精度,因此需要一个检查点对应DEM局部位置,以进行评定。在已有DSM上选取5个GPS控制点作为检查点,高度为h;找到新生成的DEM所对应点的位置,在此附近选点(10~15个点),高程值为hi,逐一和检查点进行比较,计算高程的平方差值(用RMSE表示)。检查点的实际高程与提取结果的均方误差值见表1,该方法能够很好地对DEM精度进行评定。
经过精度评定,最大的RMSE为0.45m,能够有效地滤除地物点,获取DEM精度较高。 造成误差的原因有:①倾斜影像精度及匹配生成DSM方法的精度影响DEM的精度;②区域本身的众多地形因素如地形复杂度、地形坡度等影响DEM的精度;③滤波算法影响DEM的精度。 本文的算法适合倾斜影像,执行效率高、算法稳定,能有效地滤除不同尺寸的建筑物、植被和其他地物,地面点丢失较少,可获取高精度DEM,具有一定得可靠性与实用性。 文章来源:测绘通报 更多招聘信息请登录:www.3sjob.net测绘人才网
|