一种杆塔线路点云数据校核方法与流程

1.本发明涉及户外点云数据的校核技术领域，具体涉及一种杆塔线路点云数据校核方法。


背景技术：

2.无人机对高压输电杆塔进行精细化巡检，需要对线路采集三维点云数据，以规划三维航线实现自动巡检，但在三维激光扫描仪实际采集点云数据由于受到各种条件的影响，导致数据精度不高，所以往往需要进行数据校核。传统的控制点采集方式有全站仪、手持rtk采集等，由于高压输电杆塔较多架设在山体上，地理环境复杂，传统的控制点采集方式存在效率低、视野受限、采集困难等缺点，因此难以对杆塔的点云数据进行准确地校核。


技术实现要素：

3.本发明的目的在于克服现有技术中的缺点与不足，提供一种杆塔线路点云数据校核方法，可以对杆塔的点云数据进行准确的校核，从而提高杆塔所谓区域的点云数据的准确性。
4.本发明的一个实施例提供一种杆塔线路点云数据校核方法，包括以下步骤：
5.通过无人机获取杆塔上预设的多个选取点的经纬度数据；
6.通过所述无人机获取各个所述选取点的高程数据；
7.根据所述经纬度数据和高程数据生成各个所述选取点的三维空间坐标；
8.获取所述杆塔的点云数据，根据所述点云数据获取各个所述选取点的点云坐标；
9.获取各个所述选取点的三维空间坐标与对应的所述点云坐标的第一平均差值；
10.若所述第一平均差值小于或等于所述误差值，确定所述点云数据为有效。
11.相对于现有技术，本发明的杆塔线路点云数据校核方法通过在所述杆塔上选择多个选取点，然后通过无人机获取各个选取点的三位空间坐标，再将获得的三维空间坐标与杆塔的点云数据进行比对，以判断杆塔的点云数据准确性，从而判断杆塔的点云数据是否有效或是否需要校核。
12.进一步，若所述第一平均差值大于预设的误差值，还包括以下步骤：
13.获取所述各个所述选取点的三维空间坐标对应的所述点云坐标的平均偏差参数，根据所述平均偏差参数校正所述点云数据，获得校正后的点云数据；
14.根据所述校正后的点云数据获取各个所述选取点校正后的点云坐标；
15.获取各个所述选取点的三维空间坐标与对应的所述校正后的点云坐标的第二平均差值；
16.若所述第二平均差值小于或等于所述误差值，确定所述校正后的点云数据为有效。当所述第一平均差值大于预设的误差值时，根据各个所述选取点的三维空间坐标对应的所述点云坐标的平均偏差参数对所述点云数据进行校正，从而获得校正后的点云数据，并将校正后的点云数据与所述三位空间坐标进行比对，以判断校正后的点云数据准确性。
17.进一步，当设置在所述无人机下方的镜头位于所述选取点的竖直上方，且所述镜头的中心处对准所述选取点时，获取所述镜头的中心处的经纬度坐标；
18.将所述镜头的中心处的经纬度坐标确定为对应的所述选取点的经纬度数据。通过镜头对准所述选取点，再获取镜头的经纬度坐标，从而准确地得到所述选取点的经纬度数据。
19.进一步，当所述镜头位于所述选取点的水平侧面，且所述镜头的中心处对准所述选取点时，获取所述镜头的中心处的高程；
20.将所述镜头的中心处的高程确定为对应的所述选取点的高程数据。通过镜头对准所述选取点，再获取镜头的高程，从而准确地得到所述选取点的高程数据。
21.进一步，所述获取镜头的中心处的经纬度坐标，包括以下步骤：
22.获取所述无人机的中心处的经纬度、所述无人机中心处与所述镜头的中心处的水平偏差以及镜头朝向方位角；
23.blng＝(d
·
sin(a
·
π/180)/(rj (rc-rj)
·
(90-lat)/90
·
cos(lat
·
π/180)) lng
·
π/180)
·
180/π；
24.blat＝(d
·
cos(a
·
π/180)/(rj (rc-rj)
·
(90-lat)/90) lat
·
π/180)
·
180/π；
25.其中，blng表示所述镜头的中心处的纬度，blat示所述镜头的中心处的经度，rc表示地球的椭球长半轴；rj表示地球的椭球短半轴；lat表示所述无人机的中心处的经度；lng示所述无人机的中心处的纬度；a表示所述镜头朝向方位角，d表示所述水平偏差。可以根据无人机的经纬度坐标获取镜头的中心处的经纬度坐标。
26.进一步，所述获取所述镜头的中心处的高程，包括以下步骤：
27.获取所述无人机的中心处的高程，以及所述无人机中心处与所述镜头的中心处的高程偏差；
28.balt＝alt-h；
29.其中，balt表示所述镜头的中心处的高程，alt表示所述无人机的中心处的高程，h表示所述高程偏差。可以根据无人机的高程获取镜头的中心处的高程。
30.进一步，所述无人机与rtk网络信号连接；
31.所述获取所述无人机的中心处的经纬度，包括以下步骤；
32.当所述rtk网络的解算结果为固定解时，获取所述无人机的中心处的经纬度；
33.所述获取所述无人机的中心处的高程，包括以下步骤；
34.当所述rtk网络的解算结果为固定解时，获取所述无人机的中心处的高程。降低获取的经纬度和高程的误差，从而提高获取的数据的准确性。
35.进一步，所述获取各个所述选取点的三维空间坐标与对应的所述点云坐标的第一平均差值包括以下步骤：
36.x1＝(lng1-blng)
·
π
·
rc
·
cos((lat1-blat)/2
·
π/180))/180；
37.y1＝(lat1-blat)
·
π
·
rc/180；
[0038][0039]
h1＝alt1-balt；
[0040]
[0041][0042]
其中，lng1表示所述点云坐标的经度；lat1表示所述点云坐标的纬度；alt1表示所述点云坐标的高程；d1表示对应的所述选取点的三维空间坐标与点云坐标的水平差值；h1表示对应的所述选取点的三维空间坐标与点云坐标的高程差值；s表示对应的所述选取点的三维空间坐标与点云坐标的直线差值；n表示所述选取点的个数；k表示所述第一平均差值。通过上述数据计算出所述第一平均差值。
[0043]
进一步，所述获取各个所述选取点的三维空间坐标与对应的所述校正后的点云坐标的第二平均差值，包括以下步骤：
[0044]
x2＝(lng2-blng)
·
π
·
rc
·
cos((lat2-blat)/2
·
π/180))/180；
[0045]
y2＝(lat2-blat)
·
π
·
rc/180；
[0046][0047]
h2＝alt2-balt；
[0048][0049][0050]
其中，lng2表示所述校正后的点云坐标的经度；lat2表示所述校正后的点云坐标的纬度；alt2表示所述校正后的点云坐标的高程；s'表示对应的所述选取点的三维空间坐标与校正后的点云坐标的直线差值；k'表示所述第二平均差值。通过上述数据计算出所述第二平均差值。
[0051]
为了能更清晰的理解本发明，以下将结合附图说明阐述本发明的具体实施方式。
附图说明
[0052]
图1为本发明一个实施例的杆塔线路点云数据校核方法的步骤s1-s6的流程图。
[0053]
图2为本发明一个实施例的杆塔线路点云数据校核方法的步骤s71-s74的流程图。
[0054]
图3为本发明一个实施例的杆塔线路点云数据校核方法的步骤s11-s12的流程图。
[0055]
图4为本发明一个实施例的杆塔线路点云数据校核方法的步骤s21-s22的流程图。
具体实施方式
[0056]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0057]
请参阅图1，其是本发明一个实施例的，所述杆塔线路点云数据校核方法包括以下步骤：
[0058]
s1：通过无人机获取杆塔上预设的多个选取点的经纬度数据；
[0059]
在步骤s1中，当所述无人机下方的镜头位于所述选取点的竖直上方，且所述镜头的中心处对准所述选取点时，采集所述选取点的经纬度数据。
[0060]
优选地，所述选取点的数量为至少3个，且各个选取点之间的距离至少为2m。
[0061]
s2：通过所述无人机获取各个所述选取点的高程数据；
[0062]
步骤s1中，当所述镜头位于所述选取点的水平侧面，且所述镜头的中心处对准所述选取点时，采集所述选取点的高程数据。
[0063]
在所述步骤s1和s2中，所述无人机可以通过在底部设置可以转动的镜头，以便于所述镜头在选取点的竖直上方对准所述选取点时，所述镜头处于竖直朝向所述选取点，而所述镜头在所述选取点的水平侧面，且所述镜头的中心处对准所述选取点时，所述镜头水平朝向所述选取点。而在其他实施例中，本领域技术人员还可以通过在所述无人机上分别设置朝向水平方向的镜头和竖直向下的镜头以代替上述内容说的可以转动的镜头。
[0064]
其中，所述步骤s1和s2是独立进行的步骤，因此所述步骤s1和s2之间并不存在步骤顺序的限制。
[0065]
s3：根据所述经纬度数据和高程数据生成各个所述选取点的三维空间坐标；
[0066]
s4：获取所述杆塔的点云数据，根据所述点云数据获取各个所述选取点的点云坐标；
[0067]
其中，所述杆塔的点云数据可以激光雷达采集得到。
[0068]
s5：获取各个所述选取点的三维空间坐标与对应的所述点云坐标的第一平均差值；
[0069]
s6：若所述第一平均差值小于或等于所述误差值，确定所述点云数据为有效。
[0070]
其中，所述误差值小于或等于10cm。
[0071]
在本实施例中，本发明的杆塔线路点云数据校核方法通过在所述杆塔上选择多个选取点，然后通过无人机获取各个选取点的三位空间坐标，再将获得的三维空间坐标与杆塔的点云数据进行比对，以判断杆塔的点云数据准确性，从而判断杆塔的点云数据是否有效或是否需要校核。
[0072]
请参阅图2，在一个可行的实施例中，若所述第一平均差值大于预设的误差值，还包括以下步骤：
[0073]
s71：获取所述各个所述选取点的三维空间坐标对应的所述点云坐标的平均偏差参数，根据所述平均偏差参数校正所述点云数据，获得校正后的点云数据；
[0074]
s72：根据所述校正后的点云数据获取各个所述选取点校正后的点云坐标；
[0075]
s73：获取各个所述选取点的三维空间坐标与对应的所述校正后的点云坐标的第二平均差值；
[0076]
s74：若所述第二平均差值小于或等于所述误差值，确定所述校正后的点云数据为有效。
[0077]
在本实施例中，当所述第一平均差值大于预设的误差值时，根据各个所述选取点的三维空间坐标对应的所述点云坐标的平均偏差参数对所述点云数据进行校正，从而获得校正后的点云数据，并将校正后的点云数据与所述三位空间坐标进行比对，以判断校正后的点云数据准确性。
[0078]
请参阅图3，在一个可行的实施例中，所述步骤s1包括以下步骤：
[0079]
s11：当设置在所述无人机下方的镜头位于所述选取点的竖直上方，且所述镜头的中心处对准所述选取点时，获取所述镜头的中心处的经纬度坐标；
[0080]
s12：将所述镜头的中心处的经纬度坐标确定为对应的所述选取点的经纬度数据。
[0081]
在本实施例中，通过镜头对准所述选取点，再获取镜头的经纬度坐标，从而准确地得到所述选取点的经纬度数据，避免了镜头中心点与所述无人机中心点的距离造成获取的所述选取点的经纬度数据出现误差。
[0082]
优选地，所述获取镜头的中心处的经纬度坐标，包括以下步骤：
[0083]
获取所述无人机的中心处的经纬度、所述无人机中心处与所述镜头的中心处的水平偏差以及镜头朝向方位角；
[0084]
blng＝(d
·
sin(a
·
π/180)/(rj (rc-rj)
·
(90-lat)/90
·
cos(lat
·
π/180)) lng
·
π/180)
·
180/π；
[0085]
blat＝(d
·
cos(a
·
π/180)/(rj (rc-rj)
·
(90-lat)/90) lat
·
π/180)
·
180/π；
[0086]
其中，blng表示所述镜头的中心处的纬度，blat示所述镜头的中心处的经度，rc表示地球的椭球长半轴；rj表示地球的椭球短半轴；lat表示所述无人机的中心处的经度；lng示所述无人机的中心处的纬度；a表示所述镜头朝向方位角，d表示所述水平偏差。可以根据无人机的经纬度坐标获取镜头的中心处的经纬度坐标。
[0087]
请参阅图4，在一个可行的实施例中，所述步骤s2包括以下步骤：
[0088]
s21：当所述镜头位于所述选取点的水平侧面，且所述镜头的中心处对准所述选取点时，获取所述镜头的中心处的高程；
[0089]
s22：将所述镜头的中心处的高程确定为对应的所述选取点的高程数据。
[0090]
在本实施例中，通过镜头对准所述选取点，再获取镜头的高程，从而准确地得到所述选取点的高程数据。避免了镜头中心点与所述无人机中心点的距离造成获取的所述选取点的高程数据出现误差。
[0091]
优选地，所述获取所述镜头的中心处的高程，包括以下步骤：
[0092]
获取所述无人机的中心处的高程，以及所述无人机中心处与所述镜头的中心处的高程偏差；
[0093]
balt＝alt-h；
[0094]
其中，balt表示所述镜头的中心处的高程，alt表示所述无人机的中心处的高程，h表示所述高程偏差。可以根据无人机的高程获取镜头的中心处的高程。
[0095]
在一个可行的实施例中，所述无人机与rtk网络信号连接；
[0096]
所述获取所述无人机的中心处的经纬度，包括以下步骤；
[0097]
当所述rtk网络的解算结果为固定解时，获取所述无人机的中心处的经纬度；
[0098]
所述获取所述无人机的中心处的高程，包括以下步骤；
[0099]
当所述rtk网络的解算结果为固定解时，获取所述无人机的中心处的高程。
[0100]
其中，所述rtk网络是由基准站网，数据处理中心和数据通信线路组成的。基准站上应配备双频全波长gps接收机，该接收机最好能同时提供精确的双频伪距观测值。基准站的站坐标应精确已知，其坐标可采用长时间gps静态相对定位等方法来确定。此外，这些站还应配备数据通信设备及气象仪器等。基准站应按规定的采样率进行连续观测，并通过数据通信链实时将观测资料传送给数据处理中心。数据处理中心根据流动站送来的近似坐标(可据伪距法单点定位求得)判断出该站位于由哪三个基准站所组成的三角形内。然后根据这三个基准站的观测资料求出流动站处所受到的系统误差,并播发给流动用户来进行修正
以获得精确的结果。而当所述rtk网络的解算结果为固定解时，其定位精度在1厘米以内，因此降低获取的经纬度和高程的误差，从而提高获取的数据的准确性。
[0101]
在一个可行的实施例中，所述获取各个所述选取点的三维空间坐标与对应的所述点云坐标的第一平均差值包括以下步骤：
[0102]
x1＝(lng1-blng)
·
π
·
rc
·
cos((lat1-blat)/2
·
π/180))/180；
[0103]
y1＝(lat1-blat)
·
π
·
rc/180；
[0104][0105]
h1＝alt1-balt；
[0106][0107][0108]
其中，lng1表示所述点云坐标的经度；lat1表示所述点云坐标的纬度；alt1表示所述点云坐标的高程；d1表示对应的所述选取点的三维空间坐标与点云坐标的水平差值；h1表示对应的所述选取点的三维空间坐标与点云坐标的高程差值；s表示对应的所述选取点的三维空间坐标与点云坐标的直线差值；n表示所述选取点的个数；k表示所述第一平均差值。通过上述数据计算出所述第一平均差值。
[0109]
在一个可行的实施例中，所述获取各个所述选取点的三维空间坐标与对应的所述校正后的点云坐标的第二平均差值，包括以下步骤：
[0110]
x2＝(lng2-blng)
·
π
·
rc
·
cos((lat2-blat)/2
·
π/180))/180；
[0111]
y2＝(lat2-blat)
·
π
·
rc/180；
[0112][0113]
h2＝alt2-balt；
[0114][0115][0116]
其中，lng2表示所述校正后的点云坐标的经度；lat2表示所述校正后的点云坐标的纬度；alt2表示所述校正后的点云坐标的高程；s'表示对应的所述选取点的三维空间坐标与校正后的点云坐标的直线差值；k'表示所述第二平均差值。通过上述数据计算出所述第二平均差值。
[0117]
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
[0118]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护
范围。