一种移动二维激光雷达的点云强度改正方法技术

一种地面移动二维激光的点云强度改正方法，该方法从影响强度差异的因素入手，根据同类地物反射强度相同这一原则，利用二维激光雷达提取标准漫反射板在不同扫描距离、角度下的强度数据，分别建立基于数据驱动的距离和入射角改正模型，并利用最小二乘法解算改正参数；然后采用移动二维激光扫描系统对实际场景进行移动扫描，获取实际场景的点云距离、扫描角度、强度等信息，再根据已建立的改正模型消除距离和入射角对强度值的影响，使改正后的同质区域强度趋于统一。

A Moving Two-Dimensional Lidar Point Cloud Intensity Correction Method

【技术实现步骤摘要】
一种移动二维激光雷达的点云强度改正方法
本专利技术涉及激光点云强度改正方法，尤其是一种地面移动二维激光雷达的点云强度改正方法。
技术介绍
激光强度表征为被扫描目标对发射激光光束后向散射回波的光功率。由于受到扫描仪特性、大气传输、探测器与放大电路噪声和扫描物体几何形状等因素响，激光强度值不能作为反射特性被直接提取，而是需要对其值进行改正，以消除各种因素的干扰。通过激光强度改正，扫描目标的属性信息被充分挖掘，可以开展利用目标特性进行分类的研究工作。目前研究表明，改正后的激光在森林和地表资源勘探、历史文物修复、路面和隧道破损检测、冰川表面分类等领域有着重要的价值。激光强度信息常用来辅助提高地物点云分类精度。针对不同的扫描方式，现有的激光强度改正方法，主要分为以下两种：(1)地面三维激光扫描(Terrestriallaserscanning，TLS)强度改正。地面三维激光扫描除了能获取高精度的三维空间信息外，还能记录目标的激光强度。由于地面三维激光雷达扫描距离相对较近，大气等环境因素可以忽略不计，激光强度改正主要考虑扫描距离和入射角因素。通过对实测数据的拟合，可以得到激光强度值与距离、入射角之间的关系。同时，通过反求运算，得到目标反射率。(2)机载激光扫描(Airbornelaserscanning，ALS)强度改正。机载激光扫描是依托全球定位系统、惯性测量装置、激光扫描仪获得高精度的三维点云数据，从激光强度在传播过程中的衰减机理出发，对激光传播距离、大气衰减、入射角因素进行改正，消除干扰因素对激光强度值的影响。ALS和TLS都采用成本较高的三维激光雷达，移动二维激光雷达成本低、数据冗余度小，可根据测量点的帧序号和帧内编号直接索引点云。目前还没有方法来消除距离和入射角对移动二维激光雷达强度的影响，因此，针对移动二维激光雷达强度改正方法有待深入研究。
技术实现思路
本专利技术的目的是针对移动二维激光雷达获得的同质区域点云强度差异，提出了一种地面移动二维激光的点云强度改正方法。本专利技术的技术方案是：本专利技术提供一种移动二维激光雷达的点云强度改正方法，该方法包括以下步骤：S1、利用二维激光雷达提取标准漫反射板在不同扫描距离、入射角度下的强度数据，分别建立扫描距离和强度值的函数，以及入射角度和强度值的函数；S2、采用最小二乘法解算改正参数，根据拟合优度和均方根误差，确定模型阶数；S3、采用移动二维激光扫描系统对实际场景进行移动扫描，获取实际场景的点云距离ρ、扫描角度θ，计算激光落脚点到激光原点的入射角α；S4、建立基于扫描距离和入射角度对激光强度数据的改正模型，得到改正后的激光强度IS。进一步地，步骤S1具体为：S1-1、设置入射角度、激光雷达的距离范围[Rmin，Rmax]以及激光雷达的距离间隔xm；采用激光雷达以同样的入射角度获取标准漫反射板在前述距离范围内，按照固定间隔调整距离所接收的各距离下的激光强度数据，建立扫描距离和强度值的拟合函数f2(R)；其中：R表示激光雷达到标准漫反射板的距离；Rsp表示拟合函数分段点，选取前述各距离下的激光强度最大值处的距离值，n和m分别表示f21(R)和f22(R)的多项式阶数，[a0，a1，…an]和[b0，b1，…bm]分别表示多项式的系数；S1-2、设置激光雷达的距离、激光雷达的入射角度范围[0°，90°]以及标准漫反射板的旋转间隔；采用激光雷达以同样的扫描距离获取标准漫反射板在前述入射角度范围内，按照固定旋转间隔调整角度所接收的各入射角度下的激光强度数据，建立入射角度和强度值的拟合函数f3(cosα)：f3(cosα)＝c0+c1cosα+…+ck(cosα)k其中：α表示激光雷达的入射角度；[c0，c1，…ck]表示f3(cosα)多项式系数，k表示多项式阶数。进一步地，步骤S1-1中：对于激光雷达获取的激光强度数据，取一帧数据的中间点作为当前距离下的强度值；步骤S1-2中：对于激光雷达获取的激光强度数据，取一帧数据的中间点作为当前入射角下的强度值。进一步地，步骤S1-2中：标准漫反射板的旋转间隔为10°。进一步地，步骤2具体为：S2-1、根据步骤1所建立的扫描距离和强度值的函数，以及入射角度和强度值的函数，通过最小二乘法计算，获取前述两个函数中的多项式系数[a0，a1，a2，…，an]T、[b0，b1，b2，…，bm]T、[c0，c1，c2，…，ck]T；S2-2、采用以下公式计算均方根误差(RMSE)和拟合优度(R-square)，获取前述两个函数中的模型阶数：其中：RMSE1表示参与拟合的距离-强度数据点(R(i)，I2(i))与拟合函数f2(R(i))的均方根误差；i表示距离-强度数据点的编号；M表示距离-强度数据点的总数；f2(R(i))表示编号为i的数据点的扫描距离和强度函数的拟合值，I2(i)表示编号为i的数据点对应距离下的激光雷达记录的原始强度值；RMSE2表示参与拟合的入射角-强度的数据点(αj，I3(j))与拟合函数f3(cosα(j))的均方根误差；j表示入射角-强度数据点的编号，N表示入射角-强度数据点的总数；α表示激光入射角；f3(cosα(j))表示编号为j的数据点的入射角和强度函数的拟合值，I3(j)表示编号为j的数据点对应入射角下的激光雷达记录的原始强度值；R-square1表示f2(R(i))的拟合优度，R-square2表示f3(cosα(j))的拟合优度；μ2表示原始强度I2(i)均值，μ3表示原始强度I3(j)均值；设定多项式阶数n、m和k的范围，n、m、k∈[0，7]，计算预设范围内所有系数对应的均方根误差(RMSE)和拟合优度(R-square)，选择拟合优度接近1且均方根误差最小的阶数。进一步地，步骤3具体为：步骤31、以激光雷达初始位置为坐标原点O，建立直角坐标系O-xyz，x轴方向为激光雷达在滑台上的运动方向，y轴方向为激光雷达扫描深度方向，z轴方向为被扫描目标与地面垂直的高度方向，利用二维激光雷达采集第j帧内yz平面的极坐标数据其中N为激光雷达扫描点数，坐标变换公式：其中，ρ(i′，j′)和θ(i′，j′)分别为激光雷达第j′帧内第i′个激光点的距离和扫描角度，x(i′，j′)表示第j′帧内第i′个激光点在x方向上的坐标，y(i′，j′)表示第j′帧内第i′个激光点在深度方向的坐标，z(i′，j′)表示第j′帧内第i′个激光落脚点在高度方向的坐标，Δt表示激光雷达的扫描周期，v表示激光雷达移动速度v；步骤32、计算点云的法向量，对于点云中每个点P(i′，j′)，通过最邻近算法建立邻域集，即获取最相近的k个相邻点，然后为这些点拟合一个最小二乘意义上的局部平面，即对于邻域中每个点qr，建立协方差矩阵M：其中：k为最邻近点的个数，r为最邻近点的编号，是经过最邻近点的三维质心；采用下述公式计算法向量；其中：L为特征值的编号，取值1、2、3，λL为第L个特征值，为第L个特征值对应的特征向量，；通过主成分分析，即对协方差矩阵M进行特征值分解，将M矩阵中最小特征值λmin对应的特征向量作为拟合平面的法向量，这里记作步骤33、采用下式得到激光落脚点P(i′，j′)((x(i′，j′)，y(i′，j′)，z(i′，j′))到激光原点O(本文档来自技高网...

【技术保护点】
1.一种移动二维激光雷达的点云强度改正方法，其特征在于，该方法包括以下步骤：S1、利用二维激光雷达提取标准漫反射板在不同扫描距离、入射角度下的强度数据，分别建立扫描距离和强度值的函数，以及入射角度和强度值的函数；S2、采用最小二乘法解算改正参数，根据拟合优度和均方根误差，确定模型阶数；S3、采用移动二维激光扫描系统对实际场景进行移动扫描，获取实际场景的点云距离ρ、扫描角度θ，计算激光落脚点到激光原点的入射角α；S4、建立基于扫描距离和入射角度对激光强度数据的改正模型，得到改正后的激光强度IS。

【技术特征摘要】
1.一种移动二维激光雷达的点云强度改正方法，其特征在于，该方法包括以下步骤：S1、利用二维激光雷达提取标准漫反射板在不同扫描距离、入射角度下的强度数据，分别建立扫描距离和强度值的函数，以及入射角度和强度值的函数；S2、采用最小二乘法解算改正参数，根据拟合优度和均方根误差，确定模型阶数；S3、采用移动二维激光扫描系统对实际场景进行移动扫描，获取实际场景的点云距离ρ、扫描角度θ，计算激光落脚点到激光原点的入射角α；S4、建立基于扫描距离和入射角度对激光强度数据的改正模型，得到改正后的激光强度IS。2.根据权利要求1所述的移动二维激光雷达的点云强度改正方法，其特征在于，步骤S1具体为：S1-1、设置入射角度、激光雷达的距离范围[Rmin，Rmax]以及激光雷达的距离间隔xm；采用激光雷达以同样的入射角度获取标准漫反射板在前述距离范围内，按照固定间隔调整距离所接收的各距离下的激光强度数据，建立扫描距离和强度值的拟合函数f2(R)；其中：R表示激光雷达到标准漫反射板的距离；Rsp表示拟合函数分段点,选取前述各距离下的激光强度最大值处的距离值，n和m分别表示f21(R)和f22(R的多项式阶数，[a0，a1，…an]和[b0，b1，…bm]分别表示多项式的系数；S1-2、设置激光雷达的距离、激光雷达的入射角度范围[0°，90°]以及标准漫反射板的旋转间隔；采用激光雷达以同样的扫描距离获取标准漫反射板在前述入射角度范围内，按照固定旋转间隔调整角度所接收的各入射角度下的激光强度数据，建立入射角度和强度值的拟合函数f3(cosα)：f3(cosα)＝c0+c1cosα+…+ck(cosα)k其中：α表示激光雷达的入射角度；[c0，c1，…ck]表示f3(cosα)多项式系数，k表示多项式阶数。3.根据权利要求2所述的移动二维激光雷达的点云强度改正方法，其特征在于，步骤S1-1中：对于激光雷达获取的激光强度数据，取一帧数据的中间点作为当前距离下的强度值；步骤S1-2中：对于激光雷达获取的激光强度数据，取一帧数据的中间点作为当前入射角下的强度值。4.根据权利要求2所述的移动二维激光雷达的点云强度改正方法，其特征在于，步骤S1-2中：标准漫反射板的旋转间隔为10°。5.根据权利要求2所述的移动二维激光雷达的点云强度改正方法，其特征在于，步骤2具体为：S2-1、根据步骤1所建立的扫描距离和强度值的函数，以及入射角度和强度值的函数，通过最小二乘法计算，获取前述两个函数中的多项式系数[a0，a1，a2，…，an]T、[b0，b1，b2，…，bm]T、[c0，c1，c2，…，ck]T；S2-2、采用以下公式计算均方根误差(RMSE)和拟合优度(R-square)，获取前述两个函数中的模型阶数：其中：RMSE1表示参与拟合的距离-强度数据点(R(i)，I2(i))与拟合函数f2(R(i))的均方根误差；i表示距离-强度数据点的编号；M表示距离-强度数据点的总数；f2(R(i))表示编号为i的数据点的扫描距离和强度函数的...

【专利技术属性】
技术研发人员：陶冉，李秋洁，袁鹏成，顾洲，刘旭，周宏平，
申请(专利权)人：南京林业大学，
类型：发明
国别省市：江苏,32