机载LiDAR与地基SAR数据融合三维可视化方法研究.pdf

8月0引言位移监测是变形监测的主要监测内容之一，传统的位移监测主要是利用精密水准仪、全站仪、全球卫星导航系统等技术方法来实现。这些方法具有测量成本低、观测简单等优点，但是在一些存在安全隐患以及地势起伏较大的区域实施仍有较大困难；同时，上述技术手段采用的是对点进行测量，通过以点概面的方式来进行位移监测。因此，要满足获取目标区域的形变趋势和总体形变特征的需求，需布设大量的控制点和观测点，在实际工作中较难实现。激光雷达（Light Detection and Ranging,LiDAR）通过脉冲和光电成像技术可获取探测目标高分辨率的距离、坡度、粗糙度、反射率以及数字成像等信息，经过相应处理可获取区域性的三维点云。地基SAR采用主动式微波成像雷达技术，通过合成孔径雷达(SAR)技术和步进频率连续波(SF-CW）或线性调频连续波(FM-CW）技术的结合，从而获取目标区域的高分辨率的二维影像。之后，通过雷达干涉技术处理后即得到目标的形变值1。近年来，随着地基合成孔径雷达和机载LiDAR技术的发展，地基SAR与机载LiDAR技术在地表位移监测等方面的研究已有很多成果。刘玉溪等对地基SAR在矿山应急测量中的应用进行了研究2，摘要：本文基于机载LiDAR数据与地基合成孔径雷达数据（地基SAR）的数据，采取坐标转换结合影像映射的数据融合方法，将地基SAR获取的高精度二维位移图转换并映射至机载激光雷达获取的三维点云模型上，从而获得目标区域的三维干涉雷达点云模型。其结果不仅保留目标区域高精度的二维平面位移信息，也在一定程度上实现了雷达位移影像的三维可视化，在地质灾害、边坡位移等隐患体监测中具有直观快速识别和可视化监测的良好效果。关键词：地基SAR；机载LiDAR；数据融合；变形监测；位移监测中图分类号：P237文献标识码：A文章编号：2095-7319（2023）04-0091-05Research of 3D visualization methodon airborne LiDAR and Ground-Based SAR data fusionCHEN Tingjun(Geological Engineering Survey in Fujian Province，Fuzhou 350002，China)Abstract:Based on airborne LiDAR data and Ground-Based Synthetic-Aperture Radar data(Ground-Based SAR),this paper adoptsthe data fusion method of coordinate transformation and image mapping to transform and map the high-precision 2D displacement mapobtained from Ground-Based SAR to the 3D point cloud model obtained from airborne LiDAR,so as to obtain the 3D interferometric radar point cloud model of the target area.The result not only preserves high-precision 2D plane displacement information of the targetarea,but also achieves 3D visualization of radar displacement images to a certain extent.It has a good effect in intuitive and rapid identification and visual monitoring of hidden dangers such as geological disasters and slope displacement.Key words:Ground-Based SAR;airborne LiDAR;data fusion;deformation monitoring;displacement monitoring陈廷钧（福建省地质工程勘察院，福建 福州 350002）机载 LiDAR 与地基 SAR 数据融合三维可视化方法研究收稿日期:20230508作者简介:陈廷钧（1974），男，福建惠安人，本科，工程师，主要从事新型基础测绘、智能交通及 SAR 在地灾方面的应用研究。E-mail：陈廷钧：机载 LiDAR 与地基 SAR 数据融合三维可视化方法研究91经纬天地Survey World2023年第4期No.42023论证了地基SAR在沉降监测中的可行性。祁宁等分析了机载激光雷达技术在沉降监测应用中的有效性与点云模型的直观性及可视化3。以上研究主要集中在地基SAR与机载LiDAR监测精度以及可行性研究上，对两种数据的融合涉及较少。结合两种数据的特点，本文提出坐标转换结合影像映射的方法，将地基SAR获取的高精度二维位移图转换至机载激光雷达获取的三维点云模型上，获得目标区域的三维干涉雷达点云模型。通过此方法，保留了地基 SAR 雷达位移影像高精度的同时，在一定程度上实现了雷达位移影像的三维可视化，对地质灾害隐患体具有直观快速识别和实时可视化监测的良好效果。1技术原理1.1机载激光雷达测量原理激光雷达作为一种主动式观测手段，其数据形式为雷达点云。其工作方式可简要描述为：激光雷达向测量区域的目标点发射一束激光，目标点反射激光信号，测量反射信号到达发射机的时间。通过分析反射信号的强弱以及频率变化等参数，便可以获取激光雷达相对目标点的距离、运动速度以及方向角，从而获取目标的三维点云数据。相比传统测绘手段，激光雷达具有不受天气恶劣、光照不足等不利条件制约的优点，在复杂的观测环境中仍然可以得到高精度的空间信息。1.2地基SAR监测原理地基SAR是微波遥感技术的一种，通过步进频率连续波或线性调频连续波技术可以优化观测目标距离向和方位向分辨率，同时利用合成孔径雷达技术，可提高雷达影像的空间分辨率。基于此，通过步进频率连续波或线性调频连续波技术与合成孔径雷达技术的结合，地基SAR可实现目标区域空间二维雷达图像的获取。之后，利用干涉测量技术，即可融合步进频率连续波技术或线性调频连续波技术与合成孔径雷达获取目标物在天线方向上的高精度位移。地基SAR的监测示意如图1所示。地基SAR的距离向分辨率R与光的传播速度c及地基SAR频率带宽B有关，与雷达距待测目标物之间的距离无关。其计算式如式（1）所示：R=c2B（1）角度向分辨率与雷达线性滑动轨道长度L以及雷达信号的波长有关，其计算式如式（2）所示：=2L（2）弧长距离与雷达信号到待监测目标物体之间的距离D以及角度向分辨率有关，其计算式如式（3）所示：=D（3）通过弧长即可以求出单个像素点的面积。此外，由于地基SAR空间基线为0，因此干涉相位仅与地基合成孔径雷达的形变相位defo、大气延迟相位atmo以及噪声相位Noise有关，其计算式如式（4）所示：=defo+atmo+Noise+2K（4）1.3LiDAR点云与地基 SAR 数据融合方法本文通过坐标转换以及映射相结合的方式，实现了地基 SAR 和机载LiDAR点云的数据融合，生成了三维干涉雷达点云模型。模型结果不仅保留了机载LiDAR生成目标物表面的点云，同时还融入了地基 SAR获取的高精度位移影像。模型结果具有精度高以及可视化的特点，能够精确地反映目标物的变形情况，为后期预测构筑物变形的趋势以及进行相关决策提供了支撑。数据融合主要包括3个过程：1）数据获取：通过地基SAR设备与机载激光雷达设备获取目标区域的高精度二维位移影像以及点云模型；2）数据融合：通过坐标转换实现二维位移影像与点云模型坐标系的统一，之后将二维位移影像映射至点云模型；3）结果分析：通过对融合后的三维干涉雷达点云模型的形变分析，可以对相关构筑物的形变情况进行判断，对相关构筑物的形变趋势进行预测。本文采用二维四参数模型进行坐标转换5，转换模型公式如式（5）所示：图 1地基 SAR 监测示意图928月式中：（x1，y1)与（x2，y2)为两个坐标系下的坐标点；（x，y）为平移参数；为旋转参数；m 为尺度参数。2实验及结果分析2.1数据获取为验证本文方法在实际测量中的可行性，本文实验采用Metasensing公司研制的Fast-GBSAR设备和大疆M300+L1设备对某高速边坡进行测量。分别获取目标边坡的二维位移影像以及三维点云模型。该边坡位于福建省福银高速沿线，边坡总落差约40 m，边坡顶部高程101.45 m，底部高程71.08 m，边坡总长约100 m，如图2所示。边坡主要为裸岩结构，对地基 SAR 和激光雷达扫描的回波反射率较高，信号丢失较少，有助于获取高质量的观测数据，方便后续的实验分析。首先，对边坡及其周边区域环境进行机载雷达扫描，无人机航飞时间约20，航飞后经大疆智图处理获取该边坡的激光点云数据；之后采用地基SAR对该边坡进行观测。地基SAR设备布设于距离本边坡的正前方，距离约150 m。本次观测采样频率设定为每20 min获取1张二维雷达影像，实验观测总时间约8 h，最终共计获取了28幅二维位移影像。2.2数据融合本文采取的融合方法是通过坐标转换将地基SAR系统平面坐标系（像平面坐标系）转换至机载LiDAR系统的平面坐标系（投影坐标系），从而实现平面坐标系统的统一，之后将观测方向位移影像映射至三维点位模型，从而实现地基SAR数据与激光LiDAR数据的融合，得到三维干涉雷达点云模型。2.2.1坐标转换为实现二维位移影像与点云模型的统一，采用人工选点的方法，在地基SAR图像中选取6个目标点D1D6，具体如图3所示，测量相关点在地基SAR坐标系下的坐标。同时利用机载LiDAR点云模型获取其对应目标点D1D6在机载LiDAR坐标系下的坐标。利用本文算法把地基SAR的图像与三维点云进行匹配，将目标点地基SAR坐标系下的坐标转换为机载LiDAR坐标系下的对应坐标，对转换后的坐标误差进行评定。点D1D6在地基SAR平面坐标系以及机载LiDAR平面坐标系下的坐标值如表1所示。选取D1D5五个控制点，按照上述数据融合方法计算2个坐标系的转换参数，并利用点D6验证坐标转换的精度。利用D1D5计算得到的转换参数，计算后得到坐标转换参数如下：DX=50*4.8 m，DY=33*46.3 m，=-5.9694，m=-1.9468；点位中误差为0.5 m。计算 D6由地基 SAR坐标系下的二维坐标转换到机载 LiDAR 坐标系下的坐标，评估 D6坐标转换的误差，如表2所示，并判断是否可以进行映射。转换后目标点点位中误差为0.5 m2（1.0 m），其结果存在误差，但是仍在误差允许范围内且不影响可视化的实现，满足本次坐标转换要求。本次坐标转换存在误差，这是以下原因造成的：1）地基控制点D1D2D3D4D5D6地基SARx32.2-38.336.5-30.1-57.9-11.4y401.8398.1433.1426.3419.3425.9机载LiDARx50*4.250*0.950*0.350*3.950*9.550*5.8y33*65.533*69.833*36.833*42.333*49.233*43.6表 1控制点平面坐标单位：m图 2实验区域图图 3目标点点位图（5）陈廷钧：机载 LiDAR 与地基 SAR 数据融合三维可视化方法研究93经纬天地Survey World2023年第4期No.420232.2.2影像映射在完成地基SAR的XOY平面（像平面坐标系）到机载LiDAR的XOY平面（投影坐标系）的坐标转换之后，通过映射的方法，对地基SAR测量信息与机载LiDAR的测量信息进行融合。将地基SAR二维位移变化图映射到维点云模型表面，生成目标区域的三维干涉雷达点云模型。本次实验观测时间段为8 h，每2 h为一期进行观测，生成8幅影像图，共生成32幅边坡位移变化图。分析边坡的四期位移变化情况，如图4所示。2.3结果分析2.3.1相干分析在边坡上下层选取4个像素点D1D4作为变形分析目标点，位置图如图5所示，并用目标点之间的相干系数表明目标点是否稳定可靠5，计算式如式（6）所示：式中：M与N为两幅图像；*为共轭算子。利用式（6）对5个点进行相干性分析，根据相干分析的特性，相干系数g与干涉相位的可靠性正相关，即：相干系数越接近1，其干涉相位越可靠，则位移变化值也就越准确。根据分析结果，本次选取的4个特征点间的相干系数均超过了0.97，说明其相关性好，特征点的干涉相位稳定可靠。4个特征点的位移时间变化曲线如图6所示。本次观测时间段共8 m，在观测期间内，边坡目标点未发生明显的位移，最大的位移仅4.1 mm，且各目标点的位移变化量基本一致，说明该变化量可能是由测量误差引起的。这表明该边坡在测量时段内处于稳定状态。由图6变化曲线可以看出：采用地基SAR系统对该边坡进行观测时，地基SAR雷达影像的数据位移分辨率可达mm级。在区域变形观测中，区域的形变信息可通过选取的回波信息较好的目标点来表示。2.3.2三维变形分析利用两者数据融合之后建立的三维点云模型，选取目标点D2，并分析其在测量时间段内的位移变化。D2点的二维形变曲线如图7所示，D2点在三维雷达干涉点云的位移变化如图8所示。控制点D6转换前x50*5.8y33*43.6转换后x50*5.6y33*43.1差值x0.2y0.5点位误差P0.5单位：m表 2转换后坐标SAR系统的雷达影像分辨率与机载LiDAR数据的空间分辨率不一致，导致在手动选择控制点时存在人为误差；2）边坡表面材质、纹理趋于一致，在选取控制点时存在一定困难，实际应用中可考虑采用喷漆或布靶标等方式提高控制点的识别度，减小控制点的人工选取误差，进而提高整体转换精度。图 4三维干涉雷达点云模型位移/mm图 5变形特征点位置图位移/mm时间/h图 6各特征点的时间位移图（6）94（上接第68页）3结语在困难地区的大比例尺数字地形图的质量精度检验中，传统方法存在成本消耗过大、效率低、检测范围不全等问题。为解决上述问题，研究提出无人机 OP 方法，并进一步改进常规 OP 方案，使其更加符合地形困难地区的需求，最终建立八排山地区的实景三维模型，为大比例尺数字地形图的质检工作提供新方式。试验结果显示：在平面检测与高程检测中，新型质检方式在绝大部分指标结果中的效果优于传统质检方式的效果。在质量评分中，新型质检方式的评分为 89.8 分，较传统的质检方式高1.5 分。综上所述，研究提出基于倾斜摄影测量的实景三维模型与新型质检方案具有更高的精度与适用性。但研究仍存在不足，该方法仅为静态地获取区域地理状况，但地理空间信息是随时间动态变化的，故需在未来研究中寻找动态化使用层面的质检方法。参考文献：1张志艺,张华.测绘地理信息产业公共服务平台建设研究J.地理空间信息,2020(1):6-8,16.2祁芳,张训虎,王树东.基于网格化的基础测绘地形图困难类别确定方法及应用研究J.工程勘察,2020(9):57-60.3满亚洲,朱兰艳,吕文雅.低空无人机航测系统在高原复杂地形量测中的应用J.激光杂志,2019(8):22-25.4李国锋.免像控无人机航测技术在1 500地形图测绘中的应用J.经纬天地,2019(1):21-24.8月可以看出：利用地基 SAR观测的 D2点的二维形变特征与D2在三维雷达点云上的空间位移变化趋势基本相同。从结果分析上看，融合后模型保留了雷达干涉影像高精度的位移信息，在观测段时间内，边坡变形量很小，并无明显变形。3结语本文提出了一种地基 SAR 与机载激光 LiDAR测量数据融合的方法，利用二维四参数方法将目标区域点在地基 SAR测量下坐标与机载激光 LiDAR测量下坐标进行匹配，结合映射，初步实现两者得数据融合。利用该方法对某公路边坡进行实时位移监测，通过对数据融合精度的评估与分析，验证了本文算法的可行性和有效性。结果表明：地基SAR目标点经坐标转换后的坐标与机载 LiDAR获取的三维点云模型中对应点的坐标仍存在一定误差，但该误差在观测允许范围内，不影响可视化的实现，说明利用本方法进行两者数据融合是可行的。该方法不仅可应用于边坡高精度位移监测及位移三维可视化中，同时可作为地质灾害监测等场景应用的借鉴。参考文献：1刘文婷，龙四春，张宇轩，吴文豪，祝传广.地基SAR与地基RAR模式形变监测对比J.大地测量与地球动力学，2021，41(10)：1076-1079.2刘玉溪，秦宏楠，杨凤芸，侯杉山.地基SAR用于矿山应急救援地表沉陷监测研究J.中国安全生产科学技术,2022，18(S1):25-30.3祁宁，宋心怡，李孟.机载激光雷达技术在沉降监测中的应用J.陕西煤炭,2022,41(4):186-189.4曹雪娟,阳凡林,张龙平,等.不同区域范围的二维坐标系转换方法J.工程勘察，2013,40(12):58-63.5杨鸿.面向多应用的地基差分干涉雷达高相干点选择方法研究D.重庆：重庆三峡学院，2022.图8特征点D2四期位移变化空间显示图位移/mm图7特征点D2时间位移曲线图位移/mm时间/h陈廷钧：机载 LiDAR 与地基 SAR 数据融合三维可视化方法研究95