荐言献策

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

三维坐标转换在震后建筑物LiDAR点云中的应用

杨凡 温超 范志伟 刘晓丹 王晓山 李志强 李晓丽

杨凡,温超,范志伟,等. 三维坐标转换在震后建筑物LiDAR点云中的应用[J]. 华北地震科学,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005
引用本文: 杨凡,温超,范志伟,等. 三维坐标转换在震后建筑物LiDAR点云中的应用[J]. 华北地震科学,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005
Yang Fan,Wen Chao,Fan Zhiwei,et al. Application of Three-dimensional Coordinate Transformation to LiDAR Point Clouds of Post-earthquake Buildings[J]. North China Earthquake Sciences,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005
Citation: Yang Fan,Wen Chao,Fan Zhiwei,et al. Application of Three-dimensional Coordinate Transformation to LiDAR Point Clouds of Post-earthquake Buildings[J]. North China Earthquake Sciences,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005

三维坐标转换在震后建筑物LiDAR点云中的应用

doi: 10.3969/j.issn.1003-1375.2023.03.005
基金项目: 河北省地震科技星火计划项目(DZ2021120300001)
详细信息
    作者简介:

    杨凡(1986—),男,高级工程师,主要研究方向为地震应急、遥感应用、地震灾害. E-mail:yangfan1182@126.com

    通讯作者: 温超(1981—),男,高级工程师,主要研究方向为地震灾害. E-mail:wenchao810101@126.com
  • 中图分类号: P315-39

Application of Three-dimensional Coordinate Transformation to LiDAR Point Clouds of Post-earthquake Buildings

  • 摘要: 为了解决建筑物LiDAR点云数据震害信息定量分析中坐标系转换的问题,对经过滤波、去噪、分割提取的建筑物独立墙面LiDAR点云数据进行了三维坐标转换,并与仪器配套软件的转换结果进行对比验证,坐标系转换算法的结果与x轴方向水平线吻合,优于软件转换结果。该方法的建立降低了数据处理过程中人为误差因素,提高了地基LiDAR点云数据坐标系转换的精准性和准确性,有助于提高震后建筑物LiDAR点云数据处理分析的效率,为后期建筑物震害信息定量分析判定奠定理论基础和技术支撑。
  • 图  1  采样点分布图

    Figure  1.  The distribution of sampling point

    图  2  建筑物原始点云数据可视化

    Figure  2.  Visualization of the original point cloud data of the building

    图  3  坐标系及转换示意图

    Figure  3.  Coordinate system and conversion schematic

    图  4  原始点云坐标转换对比

    Figure  4.  Comparison of coordinates transformation of original point cloud

    图  5  不同坐标转换方法结果对比

    Figure  5.  Comparison of the results of different coordinate transformation methods

    表  2  采样点数量及烈度分布明细表

    Table  2.   Number of sampling points and intensity distribution

    穿斗木结构合计共计
    宏观烈度基本完好轻微破坏中等破坏严重破坏毁坏31359
    6
    75
    82111
    94431
    砌体结构
    宏观烈度基本完好轻微破坏中等破坏严重破坏毁坏299
    686
    728431321
    8193423173
    9202329228
    框架结构
    宏观烈度基本完好轻微破坏中等破坏严重破坏毁坏29
    6
    71
    88113
    91221
    下载: 导出CSV

    表  3  坐标转换后的数据表结构(x,y,z单位m)

    Table  3.   Data table structure after coordinates transformation (x,y,z units m)

    12345678910
    xyzScalarRedGreenBlueNxNyNz
    −0.066 588−0.016 0443.657 68921151145163−0.41−0.86−0.28
    −0.070 478−0.011 9923.638 101202292382530.01−1.000.04
    −0.066 784−0.015 853.616 868172552552550.01−1.000.04
    −0.065 418−0.017 2783.596 783182482492550.01−1.000.01
    −0.065 515−0.017 1823.576 198152482492550.01−1.00−0.02
    −0.065 561−0.017 1383.555 983192442452540.01−1.000.01
    −0.068 492−0.014 0883.536 02714250248255−0.02−1.000.01
    −0.066 379−0.016 2943.515 47214249250255−0.02−1.000.01
    −0.067 311−0.015 3273.494 6615244244255−0.02−1.000.01
    −0.065 444−0.017 2793.474 535172422482550.01−1.00−0.02
    −0.065 635−0.017 0833.454 364112492502550.01−1.000.01
    −0.064 461−0.018 3113.433 614162492502550.01−1.000.01
    −0.066 353−0.016 3443.413 711182492502550.01−1.000.01
    −0.066 214−0.016 4933.393 50215248251255−0.02−1.00−0.02
    −0.066 389−0.016 3143.372 96916247249255−0.02−1.00−0.02
    −0.071 332−0.011 1663.353 73415248249255−0.02−1.00−0.02
    −0.067 71−0.014 9473.333 01816255255255−0.02−1.000.04
    −0.065 861−0.016 8793.312 57518248249255−0.05−1.000.04
    −0.062 491−0.020 3963.291 54815255255255−0.05−1.000.01
    −0.067 97−0.014 6883.272 17015254246255−0.02−1.00−0.02
    下载: 导出CSV
  • [1] 张健雄, 赵国强, 吴新辉. 基于非线性最小二乘算法的点云坐标转换研究[J]. 河南理工大学学报:自然科学版, 2012, 31(4): 429-433.
    [2] 李瑞雪. 基于地面LiDAR数据的建筑物立面识别及提取研究[D]. 赣州: 江西理工大学, 2019.
    [3] 郭迎钢, 李宗春, 何华, 等. 三维坐标转换公共点最优权值的单纯形搜索算法[J]. 测绘学报, 2020, 49(8): 1004-1013. doi:  10.11947/j.AGCS.2020.20190409
    [4] 方兴, 曾文宪, 刘经南, 等. 三维坐标转换的通用整体最小二乘算法[J]. 测绘学报, 2014, 43(11): 1139-1143. doi:  10.13485/j.cnki.11-2089.2014.0193
    [5] 徐湘宇, 王继刚, 汪玲, 等. 一种简化的小区域三维坐标系转换模型[J]. 地矿测绘, 2019, 35(1): 1-3. doi:  10.3969/j.issn.1007-9394.2019.01.001
    [6] Vasilakos C, Chatzistamatis S, Roussou O, et al. Terrestrial photogrammetry vs laser scanning for rapid earthquake damage assessment[C]//International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. Istanbul: ISPRS, 2018: 527-533.
    [7] Olsen M J, Kuester F, Chang B J, et al. Terrestrial laser scanning-based structural damage assessment[J]. Journal of Computing in Civil Engineering, 2010, 24(3): 264-272. doi:  10.1061/(ASCE)CP.1943-5487.0000028
    [8] Ziółkowski P. Remote sensing in laboratory diagnostics of reinforced concrete elements – current development and vision for the future[C]//Proceedings of TheIRES 8th International Conference. London, 2015.
    [9] Li D S, Liu J P, Feng L, et al. Terrestrial laser scanning assisted flatness quality assessment for two different types of concrete surfaces[J]. Measurement, 2020, 154: 107436. doi:  10.1016/j.measurement.2019.107436
    [10] 刘兴奇. 地面三维激光扫描仪在建筑物检测中的应用研究[D]. 北京: 北京建筑大学, 2021.
    [11] 赵兴友. 三维激光扫描仪在建筑立面测绘中的应用[J]. 测绘与空间地理信息, 2021, 44(S1): 206-208,212.
    [12] Jiao Q S, Jiang H B, Li Q. Building earthquake damage analysis using terrestrial laser scanning data[J]. Advances in Civil Engineering, 2019, 2019: 8308104.
    [13] Jiang H B, Li Q, Jiao Q S, et al. Extraction of wall cracks on earthquake-damaged buildings based on TLS point clouds[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11(9): 3088-3096. doi:  10.1109/JSTARS.2018.2857564
    [14] 崔驿宁. 基于3D深度学习的LiDAR建筑物震害信息提取[D]. 北京: 中国地震局预测研究所, 2021.
    [15] Zhang J M, Lu C Q, Wang J, et al. Concrete cracks detection based on FCN with dilated convolution[J]. Applied Sciences, 2019, 9(13): 2686. doi:  10.3390/app9132686
    [16] Liao Y J, Wood R L, Mohammadi M E, et al. Damage assessment of a sixteen story building following the 2017 central Mexico earthquake[C]//Proceedings of the 12th Canadian Conference on Earthquake Engineering. Quebec City: Canadian Association for Earthquake Engineering, 2019.
    [17] Turkan Y, Hong J, Laflamme S, et al. Adaptive wavelet neural network for terrestrial laser scanner-based crack detection[J]. Automation in Construction, 2018, 94: 191-202. doi:  10.1016/j.autcon.2018.06.017
    [18] Rastiveis H, Khodaverdi Zahraee N, Jouybari A. Object-oriented classification of lidar data for post-earthquake damage detection[C]//International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. Istanbul: ISPRS, 2018: 421-427.
    [19] Puente I, Lindenbergh R, Van Natijne A, et al. Monitoring of progressive damage in buildings using laser scan data[C]//International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. Riva del Garda: ISPRS, 2018: 923-929.
    [20] Yang F, Wen X T, Wang X S, et al. A model study of building seismic damage information extraction and analysis on ground-based lidar data[J]. Advances in Civil Engineering, 2021, 2021: 5542012.
    [21] Janalipour M, Mohammadzadeh A. Evaluation of effectiveness of three fuzzy systems and three texture extraction methods for building damage detection from post-event LiDAR data[J]. International Journal of Digital Earth, 2018, 11(12): 1241-1268. doi:  10.1080/17538947.2017.1387818
    [22] Yang F, Fan Z W, Wen C, et al. Three-dimensional point cloud analysis for building seismic damage information[J]. Photogrammetric Engineering & Remote Sensing, 2022, 88(2): 103-111.
    [23] 陈运泰, 杨智娴, 张勇, 等. 从汶川地震到芦山地震[J]. 中国科学:地球科学, 2013, 43(6): 1064-1072.
    [24] 刘杰, 易桂喜, 张致伟, 等. 2013年4月20日四川芦山M7.0级地震介绍[J]. 地球物理学报, 2013, 56(4): 1404-1407. doi:  10.6038/cjg20130434
    [25] 陈立春, 冉勇康, 王虎, 等. 芦山地震与龙门山断裂带南段活动性[J]. 科学通报, 2013, 58(20): 1925-1932.
    [26] 柯金朴. 空间三维坐标转换原理及实现[J]. 江西测绘, 2008(3): 16-17,20.
    [27] 杨卫东, 刘玉树. 一种不同坐标系之间的变换矩阵的转换方法[J]. 计算机辅助设计与图形学学报, 2000, 12(1): 53-56. doi:  10.3321/j.issn:1003-9775.2000.01.011
    [28] 胡风华, 方庆法, 张晶, 等. 浅谈三维坐标系的转换[J]. 煤炭技术, 2003, 22(11): 82-83. doi:  10.3969/j.issn.1008-8725.2003.11.044
  • [1] 孙鸿博, 孙艺, 赵楠, 解惠婷, 陆立, 荆涛.  基于“三网一员”的震后灾情速报系统设计与应用 . 华北地震科学, 2023, 41(2): 1-7. doi: 10.3969/j.issn.1003-1375.2023.02.001
    [2] 于晓虹, 叶晶, 洪赢政, 楼文高.  地震死亡人数评估的投影寻踪回归建模研究 . 华北地震科学, 2022, 40(4): 19-27. doi: 10.3969/j.issn.1003-1375.2022.04.004
    [3] 董博, 纪春玲, 张合, 李姜, 吕国军.  唐山区域概率地震危险性研究 . 华北地震科学, 2020, 38(1): 38-48. doi: 10.3969/j.issn.1003-1375.2020.01.006
    [4] 谢瑞杰, 管斌, 樊冬, 徐强, 洪阳, 汪树理.  典型矿业城市建筑物抗震安全性能调查与分析—以铜陵市为例 . 华北地震科学, 2020, 38(1): 77-83. doi: 10.3969/j.issn.1003-1375.2020.01.012
    [5] 宫玥, 申文庄, 张守洁.  地震应急信息产品的效能评估研究 . 华北地震科学, 2018, 36(4): 55-62. doi: 10.3969/j.issn.1003-1375.2018.04.008
    [6] 文翔, 毕熙荣, 向魏.  基于SVM多源遥感影像的面向对象倒塌建筑物提取研究 . 华北地震科学, 2015, 33(3): 11-17. doi: 10.3969/j.issn.1003-1375.2015.03.003
    [7] 郭斌, 王斌, 张月辉.  苏北-南黄海盆地地震构造基本特征 . 华北地震科学, 2013, 31(1): 38-44. doi: 10.3969/j.issn.1003-1375.2013.01.008
    [8] 曲利, 董晓娜, 胡旭辉, 蔡伟光, 李永红.  浅谈山东“十五”地震前兆仪器的运行维护 . 华北地震科学, 2012, 30(2): 56-59.
    [9] 冯志泽, 何钧, 任锡泰, 杜宪宋.  大中城市建筑物震害类比预测研究 . 华北地震科学, 2003, 21(3): 32-36.
    [10] 李天生, 范文.  重力变化与小汤山4.0级地震 . 华北地震科学, 1992, 10(2): 87-90.
  • 加载中
图(5) / 表 (2)
计量
  • 文章访问数:  1435
  • HTML全文浏览量:  677
  • PDF下载量:  7
出版历程
  • 收稿日期:  2022-10-20
  • 网络出版日期:  2023-07-31
  • 刊出日期:  2023-07-30

三维坐标转换在震后建筑物LiDAR点云中的应用

doi: 10.3969/j.issn.1003-1375.2023.03.005
    基金项目:  河北省地震科技星火计划项目(DZ2021120300001)
    作者简介:

    杨凡(1986—),男,高级工程师,主要研究方向为地震应急、遥感应用、地震灾害. E-mail:yangfan1182@126.com

    通讯作者: 温超(1981—),男,高级工程师,主要研究方向为地震灾害. E-mail:wenchao810101@126.com
  • 中图分类号: P315-39

摘要: 为了解决建筑物LiDAR点云数据震害信息定量分析中坐标系转换的问题,对经过滤波、去噪、分割提取的建筑物独立墙面LiDAR点云数据进行了三维坐标转换,并与仪器配套软件的转换结果进行对比验证,坐标系转换算法的结果与x轴方向水平线吻合,优于软件转换结果。该方法的建立降低了数据处理过程中人为误差因素,提高了地基LiDAR点云数据坐标系转换的精准性和准确性,有助于提高震后建筑物LiDAR点云数据处理分析的效率,为后期建筑物震害信息定量分析判定奠定理论基础和技术支撑。

English Abstract

杨凡,温超,范志伟,等. 三维坐标转换在震后建筑物LiDAR点云中的应用[J]. 华北地震科学,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005
引用本文: 杨凡,温超,范志伟,等. 三维坐标转换在震后建筑物LiDAR点云中的应用[J]. 华北地震科学,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005
Yang Fan,Wen Chao,Fan Zhiwei,et al. Application of Three-dimensional Coordinate Transformation to LiDAR Point Clouds of Post-earthquake Buildings[J]. North China Earthquake Sciences,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005
Citation: Yang Fan,Wen Chao,Fan Zhiwei,et al. Application of Three-dimensional Coordinate Transformation to LiDAR Point Clouds of Post-earthquake Buildings[J]. North China Earthquake Sciences,2023, 41(3):28-35. doi:10.3969/j.issn.1003−1375.2023.03.005
    • 地基三维激光扫描仪(Terrestrial Laser Scaning,TLS)作为新兴的主动遥感技术,具有效率高、精度高等特点,逐渐被应用到测量等各个领域,是获取地物目标高精度LiDAR点云数据的主要途径。TLS是一种利用激光反射信号进行测量的新型传感技术,可以通过发射和接收激光光束测量目标与扫描仪之间的距离,并通过激光光束发射时扫描仪的姿态角计算目标物表面所有采样点的三维坐标,同时可以获取目标物的RGB值、反射强度值及点云法向量,但是TLS获取的LiDAR点云数据坐标系是相对于仪器系统内部坐标系随机存在的[1]

      近年来,地面激光扫描技术在三维建模、文物保护、变形监测、森林结构调查、建筑物几何质量评估等不同领域得到广泛应用[2]。郭迎钢等基于最优化算法提出了一种稳健的公共点加权坐标转换算法,提高了三维坐标转换参数求解质量[3];方兴等提出适用于任意旋转角度及一般性权矩阵的三维坐标转换模型,该模型基于通用加权整体最小二乘算法,综合考虑了结构性系数矩阵、随机和非随机元素的系数矩阵等情况[4];张健雄等基于非线性最小二乘算法实现对点云坐标系的转换[1];徐湘宇等提出了基于小区域简化算法的三维坐标系转换模型[5];地面摄影测量被证明是一种准确的地震损害快速评估方法[6]。利用TLS对建筑物结构构件损坏的持续监测与分析实验,证明了其在结构建模和分析应用方面的有效性[7-8],利用TLS获取的高精度数据结合最小二乘法可应用于建筑物平面规整度的质量评价[9]。刘兴奇基于地面三维激光扫描技术进行逆向建模,对墙体垂直度和整体倾斜度进行检测[10];赵兴友通过三维激光扫描仪在建筑物立面测绘精度实例分析,验证了三维激光扫描仪在建筑物立面检测中的可行性[11]。为了从LiDAR数据中快速提取建筑物震害数据,判定建筑物的震害等级,有必要了解LiDAR数据特征和分类器的有效性。在震后建筑物的损失分析方面,学者提出了基于地基LiDAR数据的建筑物形状分析模型,有效解决了建筑物等高多边形序列提取、形状离散参数提取、不规则建筑物区块分割与震害分析等问题[12-13]。在震后建筑物墙体裂缝提取方面,学者采用平面三角剖分建模方法构建三角形不规则网络数据集,基于裂缝宽度的反距离加权点云光栅化方法生成栅格曲面,根据裂缝的形状特征提取相关信息[13]。崔驿宁通过3D点云深度学习的方法建立单体建筑物震害分类识别和震后点云场景分割模型,进行建筑物震害信息提取[14];Jianming Zhang提出了基于扩张卷积的全卷积网络用于混凝土裂缝检测[15];Yijun Liao等基于地面激光雷达扫描数据结合有限元模型分析预测建筑物的抗震性能[16];C. Vasilakos等将地面摄影测量数据与激光扫描仪的点云数据进行对比分析,验证了地面摄影测量在震后快速损失评估中的可行性[6];Yelda Turkan等提出了一种基于自适应小波神经网络的方法将TLS数据压缩成低分辨率和高分辨率的组合,实现自动检测混凝土裂缝和其他形式的损坏[17];H. Rastiveis等提出一种基于对象的方法,利用LiDAR数据绘制地震后的破坏性建筑[18];Puente等提出了一种基于拟合平面与点云距离的方法,用于监测建筑物LiDAR数据的渐进式破坏[19];Yang等利用地基LiDAR数据的高精度特点,结合Alpha shapes算法创建了三角网矢量模型(TIN-shape modle),解决了震后建筑物破坏特征提取偏差较小、肉眼无法识别的问题[20];Milad Janalipour等研究了基于LiDAR数据的三种纹理提取方法和三种模糊系统用于建立建筑物破坏图[21]。TLS高精度数据的获取为提取偏差较小、肉眼无法识别的破坏特征提供了技术帮助,利用三维激光扫描仪对建(构)筑物进行扫描,不仅可以获得建(构)筑物高精度平面数据信息,同时还能获取建(构)筑物垂直于平面的变形信息,为震后建(构)筑物的震害等级定量分析提供了高精度数据支撑,提高了建(构)筑物中肉眼难辨的倾斜和变形识别认知度及震害信息的判别精度,尤其是地震现场科考过程中对基本完好、轻微破坏建筑物识别判定模糊的问题和震后建筑物的安全隐患问题[22]。利用TLS扫描震后建筑物获取三维点云数据,快速提取建筑物震害信息能够为震后的应急救援和损失评估提供技术支撑和科学保障。

      为了实现地基LiDAR建筑物震害信息自动化提取,解决点云数据依赖仪器配套软件转换坐标系手动操作复杂繁琐的问题,本文以2013年四川芦山地震后采集的震区建筑物地基LiDAR点云数据为例开展相关研究。

    • 2013年4月20日四川省芦山县(30.3°N, 103.0°E)发生MS7.0地震,震源深度13 km,四川省芦山“4·20”7.0级强烈地震共计造成196人死亡, 21人失踪,13 486人受伤,室外避难人数约63万人,大量房屋受损破坏[23]。截至28日8时,四川省芦山“4·20”7.0级强烈地震共记录到余震5 531次,其中3.0级以上余震113次,包括5.0~5.9级4次,4.0~4.9级21次,3.0~3.9级88次[24]。芦山县位于龙门山前缘构造带南段,龙门山断裂带位于青藏高原东缘,这是2008年汶川8.0级特大地震发生5年之后,龙门山断裂带发生的又一次破坏性地震[25]

      研究团队于2013年8月进入震区调查,使用Trimble GX200三维激光扫描仪(技术指标见表1)在Ⅶ、Ⅷ、Ⅸ区开展建筑物LiDAR数据抢救性采集(图1),共计359站。其中木构架建筑31站,砌体结构(包含砖木结构、砖混结构预制板顶、砖混结构现浇顶)299站,框架结构29站(表2)。

      图  1  采样点分布图

      Figure 1.  The distribution of sampling point

      表 2  采样点数量及烈度分布明细表

      Table 2.  Number of sampling points and intensity distribution

      穿斗木结构合计共计
      宏观烈度基本完好轻微破坏中等破坏严重破坏毁坏31359
      6
      75
      82111
      94431
      砌体结构
      宏观烈度基本完好轻微破坏中等破坏严重破坏毁坏299
      686
      728431321
      8193423173
      9202329228
      框架结构
      宏观烈度基本完好轻微破坏中等破坏严重破坏毁坏29
      6
      71
      88113
      91221

      建筑物震害信息提取的基础是分离出LiDAR数据中的建筑物独立墙面点云,并进行坐标系转换。将地基三维激光扫描仪(TLS)扫描的建筑物三维LiDAR点云数据进行分割、滤波、去噪,提取建筑物独立墙面点云,建立基于独立墙面的点云坐标系统,是规范化和自动化震害信息分析的重要基础。在数据采集过程中,为提高工作效率,每个建筑物采用夹角式扫描,即同时采集建筑的正面和侧面,为降低坐标误差,采用相对坐标的形式,即三维激光扫描仪所在点为控制点坐标。采集建筑物原始点云数据如图2所示。从可视化和数据处理分析角度,相对坐标不便于后期震害信息点云数据的提取。针对此问题,本文提出了一种建筑物单体独立墙面原点坐标转换方法,转换后的点云坐标系统即定义墙体点云数据的左下角为原点坐标x/y/z的起点(0,0,0),对处理后的数据按照单位为米的规范进行转存。其中:x为横坐标(建筑物墙面长度),y为纵坐标(建筑物墙面高度),z为变形值(建筑物垂直于墙面的起伏变形量)。

      图  2  建筑物原始点云数据可视化

      Figure 2.  Visualization of the original point cloud data of the building

    • 在仪器坐标向原点坐标(建筑物平面实际坐标)转换的过程中,空间三维坐标系转换主要有3个步骤:①旋转矩阵和旋转向量;②欧拉角(旋转角);③四元数[1]。旋转欧拉角按照相对坐标(仪器坐标)转换到原点坐标(建筑物平面实际坐标)的过程,先按x轴旋转,之后按y轴旋转,最后按z轴旋转,最终得到垂直于视觉屏幕的原点坐标系[5, 26-28]。其中,旋转角度分别为αβθ,因此仪器坐标系到原点坐标系的旋转矩阵按如下定义:

      $$ {x'}=\left[\begin{array}{ccc}1& 0& 0\\ 0& \mathrm{cos}\alpha & \mathrm{sin}\alpha \\ 0& -\mathrm{sin}\alpha & \mathrm{cos}\alpha \end{array}\right] $$ (1)
      $$ y'=\left[\begin{array}{ccc}\mathrm{cos}\beta & 0& -\mathrm{sin}\beta \\ 0& 1& 0\\ \mathrm{sin}\beta & 0& \mathrm{cos}\beta \end{array}\right] $$ (2)
      $$ z'=\left[\begin{array}{ccc}\mathrm{cos}\theta & \mathrm{sin}\theta & 0\\ -\mathrm{sin}\theta & \mathrm{cos}\theta & 0\\ 0& 0& 1\end{array}\right] $$ (3)

      仪器坐标系到原点坐标系的旋转矩阵为:

      $$\begin{split} x'y'z'=&\left[\begin{array}{ccc}1& 0& 0\\ 0& \mathrm{cos}\alpha & \mathrm{sin}\alpha \\ 0& -\mathrm{sin}\alpha & \mathrm{cos}\alpha \end{array}\right]\left[\begin{array}{ccc}\mathrm{cos}\beta & 0& -\mathrm{sin}\beta \\ 0& 1& 0\\ \mathrm{sin}\beta & 0& \mathrm{cos}\beta \end{array}\right]\left[\begin{array}{ccc}\mathrm{cos}\theta & \mathrm{sin}\theta & 0\\ -\mathrm{sin}\theta & \mathrm{cos}\theta & 0\\ 0& 0& 1\end{array}\right] \\ =&\left[\begin{array}{ccc}1& 0& 0\\ 0& \mathrm{cos}\alpha & \mathrm{sin}\alpha \\ 0& -\mathrm{sin}\alpha & \mathrm{cos}\alpha \end{array}\right]\left[\begin{array}{ccc}\mathrm{cos}\beta \mathrm{cos}\theta & \mathrm{cos}\theta \mathrm{sin}\beta & -\mathrm{sin}\beta \\ -\mathrm{sin}\theta & \mathrm{cos}\theta & 0\\ \mathrm{sin}\beta \mathrm{cos}\theta & \mathrm{sin}\beta \mathrm{sin}\theta & \mathrm{cos}\beta \end{array}\right]\\ =&\left[\begin{array}{ccc}\mathrm{cos}\beta \mathrm{cos}\theta & \mathrm{cos}\theta \mathrm{sin}\beta & -\mathrm{sin}\beta \\ -\mathrm{cos}\alpha \mathrm{sin}\theta +\mathrm{sin}\alpha \mathrm{sin}\beta \mathrm{cos}\theta & \mathrm{cos}\alpha \mathrm{cos}\theta +\mathrm{sin}\alpha \mathrm{sin}\beta \mathrm{sin}\theta & \mathrm{sin}\alpha \mathrm{cos}\beta \\ \mathrm{sin}\alpha \mathrm{sin}\theta +\mathrm{cos}\alpha \mathrm{sin}\beta \mathrm{cos}\theta & -\mathrm{sin}\alpha \mathrm{cos}\theta +\mathrm{cos}\alpha \mathrm{sin}\beta \mathrm{sin}\theta & \mathrm{cos}\alpha \mathrm{cos}\beta \end{array}\right]\end{split}$$

      为了使转换完的坐标系原点o能够与建筑物平面实际坐标系oʹ完全吻合,需要引入平移矩阵,即仪器坐标系在建筑物平面实际坐标系下的坐标δ,其变换方程如下:

      $$ {X}_{\text{建}}={X}_{\text{仪}}{x'}{y'}{z'}+{\delta }_{\text{仪}\to \text{建}} $$ (4)

      式中:X为建筑物平面实际三维坐标;X为仪器采集三维坐标。

    • 建筑物地基LiDAR点云数据经过滤波、去噪、分割等步骤后,将建筑物独立墙面LiDAR数据提取保存为满足本研究需要的文件格式,最关键的环节是LiDAR数据的坐标系转换。地基三维激光扫描仪采集获取的建筑物LiDAR数据是仪器自定义坐标系的相对坐标,为了便于建筑物震害信息提取并达到数据最佳的可视化效果,需将数据采集时随机获取的三维坐标转换为墙体垂直于视觉屏幕的原点(0,0,0)的坐标(图3)。

      图  3  坐标系及转换示意图

      Figure 3.  Coordinate system and conversion schematic

      本文依据空间三维坐标系转换算法,结合本研究仪器采集的数据结构特点,建立了适用于本文的坐标系转换方法。

      仪器采集的建筑物点云数据,x为场景横坐标,y为建筑物高度坐标,z为垂直于xoy面的坐标。对建筑物三维点云数据分割后,提取建筑物独立墙面并建立独立坐标系,用于后期震害信息分析。其中,x′为墙面的长度,y′为墙面的高度,z′为墙面的倾斜变形值,此时构建的oxyz′坐标系即为原点坐标系(建筑物平面实际坐标系)。本文坐标系的转换方法如下:

      1)确定原点坐标系原点o′,即平移到o′为零点坐标,选取原始数据的xz值,利用最小二乘法拟合一元一次直线方程,即:

      $$ z=Fx+\delta $$ (5)

      式中:Fxz值的拟合直线斜率,该斜率值的反正切值(arctan(F))即为欧拉角(墙面旋转角度)θ

      在原始数据中选取xy值的最小值xminymin,并计算求得xmin在公式(5)中的zz1,即o’点坐标为(xminyminz1)。

      通过平移向量方程即可将所有仪器坐标点平移到原点坐标系。

      $$ \left.\begin{array}{c}{\mathrm{x}}'={x-x}_{min}\\ {y}'={y-y}_{min}\\ z'={z-z}_{1}\end{array}\right\} $$ (6)

      2)在坐标系旋转过程中,遵循右手法则,本文的数据结构特点是yy′平行,在旋转过程中y值不变,将墙面点云数据投影到xoz平面(图3),xz为仪器坐标系,xz′为建筑物平面实际坐标系。其本质就是将xzy为轴转到xz′,其旋转角度为上一步骤计算得到的θ角,即$ {x′}=x\mathrm{cos}\theta -z\mathrm{sin}\theta $$ {z′}=x\mathrm{sin}\theta + z\mathrm{cos}\theta $。引用坐标旋转矩阵,转换后的xyz′坐标系坐标为:

      $$ (x'y'z')=\left(xyz\right)\left\{\begin{array}{ccc}\mathrm{cos}\theta & \mathrm{sin}\theta & 0\\ -\mathrm{sin}\theta & \mathrm{cos}\theta & 0\\ 0& 0& 1\end{array}\right\}\left\{\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right\} $$ (7)

      式中:θ=arctan(F);矩阵$ \left\{\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right\} $用于调整坐标轴方向(正负值)。

      结合上述算法分析,对建筑物独立墙面原始点云数据进行坐标系转换,其可视化结果与扫描仪配套软件Realworks的转换结果基本吻合(图4)。转换后的数据表结构如表3,其中,x为墙面的长度,y为墙面的高度,z为垂直于墙面的起伏变形量,Scalar为仪器信号反射强度,Red/Green/Blue为颜色色彩值,Nx/Ny/Nz为点云法向量信息。当z值为正,则该点为垂直于墙面向内凸出;当z值为负,则该点为垂直于墙面向外凸出。

      图  4  原始点云坐标转换对比

      Figure 4.  Comparison of coordinates transformation of original point cloud

      表 3  坐标转换后的数据表结构(x,y,z单位m)

      Table 3.  Data table structure after coordinates transformation (x,y,z units m)

      12345678910
      xyzScalarRedGreenBlueNxNyNz
      −0.066 588−0.016 0443.657 68921151145163−0.41−0.86−0.28
      −0.070 478−0.011 9923.638 101202292382530.01−1.000.04
      −0.066 784−0.015 853.616 868172552552550.01−1.000.04
      −0.065 418−0.017 2783.596 783182482492550.01−1.000.01
      −0.065 515−0.017 1823.576 198152482492550.01−1.00−0.02
      −0.065 561−0.017 1383.555 983192442452540.01−1.000.01
      −0.068 492−0.014 0883.536 02714250248255−0.02−1.000.01
      −0.066 379−0.016 2943.515 47214249250255−0.02−1.000.01
      −0.067 311−0.015 3273.494 6615244244255−0.02−1.000.01
      −0.065 444−0.017 2793.474 535172422482550.01−1.00−0.02
      −0.065 635−0.017 0833.454 364112492502550.01−1.000.01
      −0.064 461−0.018 3113.433 614162492502550.01−1.000.01
      −0.066 353−0.016 3443.413 711182492502550.01−1.000.01
      −0.066 214−0.016 4933.393 50215248251255−0.02−1.00−0.02
      −0.066 389−0.016 3143.372 96916247249255−0.02−1.00−0.02
      −0.071 332−0.011 1663.353 73415248249255−0.02−1.00−0.02
      −0.067 71−0.014 9473.333 01816255255255−0.02−1.000.04
      −0.065 861−0.016 8793.312 57518248249255−0.05−1.000.04
      −0.062 491−0.020 3963.291 54815255255255−0.05−1.000.01
      −0.067 97−0.014 6883.272 17015254246255−0.02−1.00−0.02

      为了进一步验证坐标转换后的点云数据质量,将墙体点云变形值z投影到地面,用以计算x方向z值线性趋势线的水平度是否与x轴平行(图5)。其中,绿色线为通过Realworks软件手动转换坐标系的结果,转换后代表墙面起伏变化的点云z值线性趋势线R2为0.032 4,与x轴方向水平线存在一定误差;蓝色线为本文算法自动转换坐标系的结果,转换后代表墙面起伏变化的点云z值线性趋势线R2为0.001,与x轴方向水平线吻合。经过对比分析发现,软件自带坐标系转换方法在手动选点确定坐标轴方向时,受墙面点云变形值起伏变化的影响,不易确定位于同一水平面的2个点。坐标系转换算法的建立,避免了这一误差的出现,提高了地基LiDAR点云数据坐标系转换的精准性和准确性,为后期建筑物震害信息定量分析判定奠定理论基础和技术支撑。

      图  5  不同坐标转换方法结果对比

      Figure 5.  Comparison of the results of different coordinate transformation methods

    • 本文对经过滤波、去噪、分割提取的建筑物独立墙面LiDAR点云数据开展了三维坐标系转换,并与软件转换的数值结果进行了对比。坐标系转换算法的结果与x轴方向水平线吻合,优于软件转换结果。该算法的建立,降低了数据处理过程中人为误差因素,提高了地基LiDAR点云数据坐标系转换的精准性和准确性,为后期建筑物震害信息定量分析判定奠定理论基础和技术支撑。

参考文献 (28)

目录

    /

    返回文章
    返回