第43卷 第2期
华北理工大学学报(自然科学版)V o l .43 N o .22021年04月J o u r n a l o fN o r t hC h i n aU n i v e r s i t y o f S c i e n c e a n dT e c h n o l o g y (N a t u r a l S c i e n c eE d i t i o n )A p
r .2021 收稿日期:2020-09-23 修回日期:2021-03-29第一作者：李京，男，硕士研究生㊂研究方向：地理信息技术与应用㊂ 通讯作者：刘亚静，女，博士，教授㊂研究方向：地理信息技术与应用㊂E -m a i l :l y j 2206@163.c o m D O I :10.3969/j
.i s s n .2095-2716.2021.02.004文章编号:2095-2716(2021)02-0024-06
李京，刘亚静，刘明月
(华北理工大学矿业工程学院，河北唐山063210
)关键词：点云数据；震后建筑物；北川遗址；三维建模；三维激光扫描
摘 要：三维激光扫描技术具有高效率㊁高精度和非接触式主动测量的特点，为震后建筑物三
维可视化建模提供一种全新手段㊂以北川旧县城遗址的单体不规则建筑为例，利用三维激光
扫描仪所获取的激光点云数据，经过配准拼接㊁数据去冗㊁去噪滤波等预处理操作，基于
G e o m a g i cS t u d i o 平台完成震后建筑物三维可视化建模，并对构建模型进行偏差分析以及精度验证㊂结果表明，该方法可为震后建筑损毁程度及灾情评估提供有效的技术支持㊂ 中图分类号:P 627 文献标识码:A
0引言
进行震后建筑的三维建模对于充分了解震后建筑损毁程度以及进行灾情评估具有重要意义㊂传统的建筑物震害调查以实地获取震害照片为主，但由于数据采集能力的限制，很难快速获取高精度的建筑物三维数据[1]㊂随着现代测量技术不断更新与发展，无人机摄影测量㊁遥感技术以及三维激光扫描仪等手段不断应用于震后建筑震害评估[ 1-3]㊂其中三维激光扫描仪可高精度快速地获取目标物体表面的三维点云数据并且突破传统测量方法单点测量局限性的优势，目前已广泛应用于林业调查㊁地形测绘㊁变形监测㊁考古与文物保护
及城市建设等多个领域[4-7]㊂基于激光点云建立震害建筑物的三维可视化模型，并将震后建筑物的真
实表面以数字的形式长期完整地保存下来，有利于建筑物的灾后重建，并且进行震害信息提取从而进行联合分析,
有利于满足对于地震遥感应急评估的需求，并为之提供可用的数据源评估㊂如李强等以北川县城扫描获取的地面激光雷达数据为实验数据，将面向对象方法引入地面L i D A R 点云建筑物提取中,
实现了震害建筑物自动提取[3]；焦其松等采集汶川地震震后彭州市白鹿中学的精细点云数据，对建筑物的典型破坏特征及地表破裂情况进行了研究，再现了白鹿中学破坏场景，实现了震害信息的永久保存[8]；沈大勇等以震后北川县城为例，探讨了震后三维场景构建过程中，通过几何修正和纹理映射实现损毁道路三维可视化[9];D o m i n i k B r u n n e r 等提出利用震前高分辨率光学影像模拟震后S A R 影像，与震后S A R 影像进行变化检测来提取震
害建筑物信息的方法[10]㊂R e h o r 等人利用震后机载L I D A R 数据提取了震害信息[1
1]㊂该项研究选取北川县城遗址内单体建筑物为研究对象，通过布置方案全面扫取研究区域的震害现场，获取完整的震后建筑物表面三维激光点云数据，然后选取部分震后受损的建筑点云数据并进行配准拼接㊁点云分割㊁去噪滤波等预处理操作㊂之后基于三维激光点云数据实现震后建筑物的三维建模，实现三维实体重建，并给出受损建模具体过程以及精度分析㊂
1数据与方法
1.1 点云数据的获取为了获取完整的三维激光点云数据，应遵从以下布设原则：首先，通过对现场进行实地踏勘，根据测区现场地形特征和地物特征布设作业范围，并布设扫描站点，合理安排扫描仪器的位置确保相邻站点数据具有重叠部分㊂同时尽可能减少测站数，以便减少数据量和测站之前的拼接误差㊂其次，确保遮挡较少，避开其他人员干扰，从而使数据包含完整的建筑信息㊂
该数据采集利用R i g e lV Z -400三维激光扫描仪对建筑物进行扫描，扫描速度为125000p t s /s ，扫描仪的点间距<1mm ；距离精度为ʃ5mm ，扫描视场范围为360ʎˑ100ʎ(水平ˑ垂直)，设置此次的扫描模式为：长距模式㊂在不同的站点获取点云数据，使建筑物完全被扫描，确保目标不被漏点㊂
1.2 点云数据预处理
1.2.1数据拼接
为了获取目标建筑物点云数据全貌，将在不同扫描站所获取的点云数据进行点云拼接㊂点云拼接是通过2个扫描站间点云数据的对应关系，将2个坐标系下的点云数据转换到同一个坐标系中，也称之为点云配准[12]㊂点云拼接分为粗配准和精配准，当粗配准精度不满足要求时，利用精配准进一步提高精度㊂粗配准，点云数据的拼接使用R i s c a nP r o 软件对点云数据进行初步处理,
以选取同名点的配准方法，相邻两站依次进行配准㊂在选取同名点时，尽可能精准选取两站重合的同
名点，以便减小误差㊂拼接完成之后可为每站的点云赋予不同的颜，从而观察拼接之后的效果㊂常采用的精配准算法是I C P 迭代法(I t e r a t i v eC l o s e s t P o i n t )也叫迭代最近点法，该算法是由B e s l 和M c k a y 在1992年提出[12]，通过利用最小二乘迭代法使平移参数和旋转参数最优，并使误差函数最小㊂该项研究通过点云拼接之后的精度可达到4.6mm ㊂1.2.2数据去冗
为了提高数据处理效率和对电脑软硬件的配置要求㊂需要数据去冗，数据去冗分为粗删和精删㊂粗删主要是对震后建筑物研究对象以外无关的事物进行删除，而对于建筑物倒塌的部分进行保留，以便还原㊁保留更真实的震后情况㊂精删主要是对局部具体细节，如通过透过门窗等扫描的点云在建筑物内部进行剔除㊂
1.2.3去噪滤波
为避免数据量过大，致使在R i s a nP r o 软件中无法生成p o l y d a t a ，引起数据导出失败，在数据预处理过程中需要对数据进行多次滤波和抽稀㊂在不影响数据质量的情况下尽可能减少数据量，以便在软件中导出数
据，并提高点云计算速率从而提高建模效率㊂去噪滤波是通过对激光扫描回波信号强度的辨别来对点云数据中错误点及粗差点进行剔除㊂通过中值滤波剔除奇异点，利用曲面拟合去除前端遮挡物，当回波信号强度
低于阈值时，距离信号值为无效[13,14]㊂以上数据拼接数㊁数据去冗以及去噪滤波主要通过R i s a nP r o 软件进
行，在完成数据预处理之后，以文件形式保存提取出的目标建筑物并用于三维建模㊂
1.3 基于点云的震后建筑物三维建模方法
G e o m a g i cS t u d i o 具有强大的逆向建模功能，在我国已得到广泛的应用㊂该软件遵循点阶段 多边形阶段 曲面阶段的三阶段作业流程，可以轻易地从点云创建出完美的多边形模型和样条四边形网格，并可自动
转换为N U R B S 曲面，建模效率高[15]㊂在通过R i s a nP r o 软件完成点云数据预处理基础上，基于G e o m a g
i c S t u d i o 软件实现震害建筑物三维可视化建模,
通过数据的转换导入㊁分割以及处理完成目标建筑物的模型的构建，在处理过程中一般分为点㊁多边形和曲面3个阶段㊂图1所示为震后建模流程㊂5
2 第2期 李京,
等：基于激光点云数据的震后建筑物三维可视化建模
图1 震后建模流程1.3.1点阶段
在点云预处理阶段已进行拼接㊁冗余剔除以及去噪滤波㊂点阶段继续进行点云阶段的处理㊂首先利用
C l o u dC o m p a r e 转换为p l y 格式，导入G e o m a g i cS t u d i o 软件㊂为使得建筑物特征明显㊁清晰以及提高计算速率和建模效率，在点阶段进行数据分割㊂点云分割就是将同一片曲面上获得的点云数据赋予相同标志的
过程[16]㊂点云分割的方法主要有2种：手动分割和自动分割，其中手动分割适用于数据特征明显，分割次数少的情况㊂目前R a i n D r o p 公司推出的软件G e o m a g i cS t u d i o 未对点云数据分割进行深入研究和应用㊂所以该项目运用软件进行手动分割，在分割过程中要尽可能减少分割次数，使分割后的数据凸显各个形状特征，确保模型建模贴近真实情况，从而使模型在后期封装的精度更高㊁特征更加明显㊁可视化效果更好㊂
1.3.2多边形阶段
首先进行封装，所谓封装就是通过点云数据创建多边形的过程㊂点云的封装主要有曲面封装和体积封装㊂曲面建模从局部进行三角网格化，然后不断延伸直至实现全局网格化，曲面封装建模速率较快，并且节省计算机内存，主要用于数据均匀或完整的曲面㊂体积封装主要用于数据不完整不均匀的数据，需要时间较
多，计算速率较慢[17]㊂该项研究采用曲面封装的方法，利用封装工具对点云进行封装，完成基本模型的初步建立㊂在完成封装后，再对一些之前漏掉的杂点封装成的点云集合进行去除，对于与所要建模的建筑物无关
的数据进行删除或者去除特征处理㊂
由于所获得的点云的质量问题而导致封装后所建立的基本模型缺少点云数据而产生孔洞㊂需要对孔洞进行修补，修补处应尽量与周围特征面保持一致㊂对于模型曲面上的凸形部分，可以选择去除特征或是删除钉状物功能加以去除，并且可以使用快速平滑和砂纸功能使模型表面更加平滑㊂在分块的模型合并的前后，为防止由于封装成型的三角形数量过多而导致在显示时形成点云状，需要对模型进行简化，在简化时选取
70%~90%之间㊂选取百分比的太小要适中,
太少会使简化效果不好，太多会导致丢失掉许多细节㊂模型合并主要有2种方法：一种是直接对点云进行拼合，另一种是各个部分构造几何形体最后再拼接几何形体[
17]㊂该项研究选取的是先将各个部分的点云构造几何体处理好后再进行模型的合并㊂把最初分割后处理好的模型通过导入的方式导入到一起，由于模型的坐标系是相同的，所以在分割后再次合并的时候可
以相互吻合㊂选取2个模型通过合并工具进行合并，使之成为一个完整的模型㊂1.3.3曲面阶段
曲面重构就是根据点云数据构造出一个用精确的数学函数描述的曲面㊂一般采用过插值和近似这2种方法进行曲面的构建，其中近似法只需要逼近原始数据即可㊂
三维曲面模型构建原理:
f (x ,y ,z )=a 1x 2+a 2y 2+a 3z 2+a 4x y +a 5y z +a 6y z +a 7x +a 8y +a 9
z =0(1)其中x ㊁y ㊁z 分别为点云数据的三维坐标,a i 为拟合系数㊂基于曲面的模型重建时直接对测量数据进行曲面片的拟合，然后再将获得的曲面片经过过渡㊁混合㊁连
62 华北理工大学学报(自然科学版) 第43卷
接形成最终曲面模型的过程㊂
1.4 偏差分析方法
1.4.1仪器测量误差
点云数据误差主要来源于测距和测角的误差㊂这两方面对精度的影响为: m 1=m 2d +m 2β/ρ2S 2(2)其中,m 1为点云数据总误差,m d 为仪器测距误差,m βm β为仪器测角误差,S 为测站与目标间的水平距离㊂点云的位置信息对距离误差影响较大㊂距离误差不仅取决于所测目标的距离，还受到测量次数的影响㊂距离越长，误差越大㊂测量次数越少，误差越大㊂
R i g e lV Z -400三维激光扫描仪的测距精度为3mm+2p p m ，测角误差最大不超过0.5ᵡ，测站距离目标在100m 以内㊂代入(2)式可得,m 1=3.2mm ㊂1.4.2三维建模误差
根据建模应用分析之后可得建模误差m 3分别为0.412mm 和0.030mm ㊂根据误差传播定律可得，点位误差的理论值为:
m =m 21+m 22+m 23(3
)将以上数据代入(3)中，得m 分别为5.62mm 和5.60mm ㊂2建模结果及其应用偏差分析
2.1 建模结果
基于由R i g e lV Z -400扫描获取的激光点云数据，采用G e o m a g
i cS t u d i o 实现对震害受损建筑物进行三维可视化建模㊂建模结果如图2所示㊂图2 震害受损建筑物三维可视化建模结果
2.2 模型的应用分析建模完成之后对模型进行偏差分析，最大偏差分别为
3.814mm ㊁-3.829mm ㊁1.425mm ㊁-2.076mm ㊂模
型与点云数据的标准偏差分别为为0.412mm ,0.030mm ㊂通过偏差计算得出,2个受损建筑精度分别为5.62mm ，和5.60mm ㊂对震后建筑物的重建不仅是形状上的还原同时还对还原的精度进行了计算和说明,7
2 第2期 李京,
等：基于激光点云数据的震后建筑物三维可视化建模
82华北理工大学学报(自然科学版)第43卷经偏差分析表明，该模型可获得较好的精度，并能够很好地保存实体的灾害信息，从而实施三维可视化的浏览和定位量测㊂如图3㊁图4所示为偏差分析㊂
图3偏差分析a
图4偏差分析b
3结论
(1)以北川旧县城遗址的部分震后建筑物为例，利用R i g e lV Z-400三维激光扫描仪获取点云数据，先后进行了点云数据拼接㊁特征点提取和数据转换等预处理，利用G e o m a g i cS t u d i o完成了震后建筑物三维可视化建模，并进行了偏差分析㊂
(2)三维激光扫描技术作为新兴的技术手段，它在一定程度上可提高震害建筑物提取精度，同时还能较好地应用于地震应急与震害评估中，不仅可为之后的震后灾情研究提供有效而又丰富的实体信息，并可为地震灾害精准评估及受损建筑物复原提供科学参考，还可为遥感应急救援与评估增添新的可用数据源㊂
参考文献:
[1]薄景山，李孝波，段玉石，等.九寨沟7.0级地震学校建筑震害调查与分析[J].地震工程与工程振动,2018,38(05):120-128.
[2]袁小祥，王晓青，丁香，等.基于无人机影像的九寨沟地震建筑物震害定量评估[J].中国地震,2017,33(04):582-589.
[3]王晓青，王龙，章熙海，等.汶川8.0级地震震害遥感定量化初步研究 以都江堰城区破坏为例[J].地震,2009,29(01):174-181.
[4]李强，焦其松，张景发.基于地面激光雷达技术点云的北川县城震害建筑物提取研究[J].科学技术与工程,2016,16(19):244-249.
[5]朱赞，甘淑，袁希平，等.三维激光扫描技术对岩溶区公路滑坡要素提取[J/O L].测绘科学,2019(11):
1-9[2019-11-02].h t t p://k n s.c n k i.
n e t/k c m s/d e t a i l/11.4415.P.20190805.0917.006.h t m l.
[6]李波.机载L i D A R点云数据在城市遥感信息提取中的应用探索[J].测绘通报,2018(s1):66-69.
[7]刘芳，冯仲科，杨立岩，等.基于三维激光点云数据的树冠体积估算研究[J].农业机械学报,2016,47(03):328-334.
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议