一種簡(jiǎn)化的輸電線路點(diǎn)云電塔自動(dòng)定位方法

危雙豐， 尹思銘， 習(xí)曉環(huán)， 王 濮

(1.北京建筑大學(xué) 測(cè)繪與城市空間信息學(xué)院， 北京 100044； 2.北京建筑大學(xué) 北京未來城市設(shè)計(jì)高精尖創(chuàng)新中心， 北京 100044；3.北京建筑大學(xué) 代表性建筑與古建筑數(shù)據(jù)庫教育部工程研究中心， 北京 102616；4.中國科學(xué)院遙感與數(shù)字地球研究所， 北京 100094)

收稿日期：2019-05-14

基金項(xiàng)目：國家自然科學(xué)基金項(xiàng)目(41601409)； 北京市自然科學(xué)基金項(xiàng)目(8172016)

第一作者簡(jiǎn)介：危雙豐(1979—)，男，副教授，博士，研究方向：點(diǎn)云智能處理.

激光雷達(dá)技術(shù)在架空輸電線路中的應(yīng)用已經(jīng)得到國內(nèi)外電網(wǎng)公司的廣泛認(rèn)可，如線路規(guī)劃選線與施工、安全巡檢、危險(xiǎn)點(diǎn)預(yù)警與分析、數(shù)字電網(wǎng)建設(shè)等，特別是對(duì)于線路穿越地形復(fù)雜、人員難以到達(dá)的區(qū)域更顯示其突出優(yōu)勢(shì)[1-3].

激光雷達(dá)在輸電線路中的應(yīng)用主要包括輸電線路通道數(shù)據(jù)獲取、點(diǎn)云去噪、濾波分類、數(shù)字化建模與分析等，其中濾波分類是后續(xù)應(yīng)用中不可或缺的流程，直接影響后續(xù)建模的精度和分析應(yīng)用結(jié)果. 濾波分類的目的是將輸電通道點(diǎn)云數(shù)據(jù)分為地面點(diǎn)、植被點(diǎn)、建筑點(diǎn)、電力線點(diǎn)、電塔點(diǎn)云等，在此基礎(chǔ)上進(jìn)行各要素的建模與分析. 電塔作為輸電線路的重要組成部分，其安全穩(wěn)定性對(duì)線路穩(wěn)定運(yùn)行至關(guān)重要，特別是山區(qū)地形復(fù)雜、人員難以到達(dá)的區(qū)域，更需要對(duì)電塔進(jìn)行三維數(shù)字重建、定期/不定期進(jìn)行變化分析，如電塔位置的精細(xì)識(shí)別與動(dòng)態(tài)變化等. 國內(nèi)外學(xué)者開展了很多有關(guān)電塔定位的研究. 彭向陽等[3]3670-3677利用電力走廊區(qū)域點(diǎn)云的高差、密度和坡度特征提取電塔位置，但坡度特征易受到地表高大植被的影響而效果不佳，且其結(jié)果可用高差特征替換. 虢韜等[4]利用電塔的空間分布特征，首先粗提電塔點(diǎn)云，再分離電塔，最后利用最小二乘法精確擬合電塔位置. 該方法提取點(diǎn)云較為復(fù)雜，為了簡(jiǎn)化提取過程，本文利用電塔點(diǎn)云的特征直接提取其位置. 韓文軍[5]首先提取電力線點(diǎn)云，然后通過檢測(cè)電力線連接點(diǎn)定位電塔. 其方法過于依賴電力線，若電力線有異常情況(如異物、斷裂等)則難以取得良好的效果，限制了其在災(zāi)害監(jiān)測(cè)方面的應(yīng)用能力. 游安清等[6]將點(diǎn)云數(shù)據(jù)格網(wǎng)化后向水平方向投影降維成二維灰度圖像，提取其直方響應(yīng)最大的位置為電塔位置. 但進(jìn)行桿塔定位算法不夠精細(xì)，且沒有經(jīng)過較大規(guī)模數(shù)據(jù)的驗(yàn)證. 劉洋等[7]進(jìn)行了改進(jìn)，增加利用懸空特征點(diǎn)集來輔助識(shí)別電塔，但也使方法變得更加復(fù)雜. 葉嵐等[8]在提取電力線的基礎(chǔ)上提出利用數(shù)學(xué)形態(tài)學(xué)進(jìn)行電塔識(shí)別定位的方法. 其在電塔識(shí)別上有創(chuàng)新的方法，但主要針對(duì)電力線的提取和擬合. 賴廣陵等[9]提出基于車載激光雷達(dá)數(shù)據(jù)進(jìn)行電塔定位的方法. 經(jīng)過對(duì)比，機(jī)載激光點(diǎn)云地面噪聲點(diǎn)更少、電塔頂部信息更全面，更利于進(jìn)行電塔定位，且機(jī)載點(diǎn)云獲取方法更為方便，有更好的應(yīng)用前景. 林祥國等[10]則是直接基于已提供的電塔二維信息精確定位每座塔的位置，此方法需要先驗(yàn)數(shù)據(jù)的支持，因此適用度不高. 上述方法都為電塔定位的研究提供了良好的思路.

本文針對(duì)前人算法進(jìn)行改進(jìn)，針對(duì)電塔自身和機(jī)載激光雷達(dá)數(shù)據(jù)的特點(diǎn)，設(shè)計(jì)并整合了較為有效的提取、定位手段，簡(jiǎn)化了總體處理流程. 本算法利用電塔區(qū)域點(diǎn)云的大高差、高密度等特點(diǎn)，提出了一種對(duì)點(diǎn)云構(gòu)建二維格網(wǎng)，根據(jù)密度、歸一化數(shù)字表面模型(Normalized Digital Surface Model，nDSM)統(tǒng)計(jì)出興趣區(qū)(Area of Interest，AOI)，最后根據(jù)AOI的高程空間分布特征提取電塔坐標(biāo)的算法，用試驗(yàn)驗(yàn)證了算法的穩(wěn)健性和效率.

1 算法流程

在二維平面空間上，輸電線路通常呈折線段，其節(jié)點(diǎn)為電塔；在三維空間中，電塔通常具有特定的高度(與輸電電壓有關(guān))，具有一定的高程空間分布特征. 此外，由于機(jī)載激光雷達(dá)和電塔自身的特點(diǎn)，激光可以穿透電塔的鏤空結(jié)構(gòu)并返回得到其內(nèi)部構(gòu)造，使得電塔區(qū)域內(nèi)的點(diǎn)云投影到二維平面上時(shí)密度較高. 因此，本文根據(jù)電塔區(qū)域點(diǎn)云數(shù)據(jù)的nDSM特征和密度特征，對(duì)其進(jìn)行濾波、篩選，得到電塔的潛在AOI. 為了進(jìn)一步分離電塔點(diǎn)云，剔除周圍高大喬木等的影響，需要統(tǒng)計(jì)AOI的高程空間特征，從而排除干擾，提取出電塔區(qū)域并計(jì)算位置，實(shí)現(xiàn)電塔的自動(dòng)定位. 電力桿塔自動(dòng)定位流程如圖1所示.

1.1 預(yù)處理

由于機(jī)載激光雷達(dá)儀器存在噪聲，使得點(diǎn)云數(shù)據(jù)中存在一定量的孤立點(diǎn)和噪聲點(diǎn). 孤立點(diǎn)和噪聲點(diǎn)的存在會(huì)影響nDSM特征的判斷、提取，并降低電塔定位的最終精度，因此需要對(duì)原始數(shù)據(jù)進(jìn)行降噪處理. 本文使用的降噪方法為半徑異常值去除法. 其基本原理是設(shè)定一個(gè)半徑值，若待定點(diǎn)周圍半徑區(qū)域內(nèi)所含的其他點(diǎn)數(shù)小于給定的閾值(一般為經(jīng)驗(yàn)值)，則該點(diǎn)被視為噪聲點(diǎn)，該方法能夠有效去除離散點(diǎn).

此外，受限于無人機(jī)激光雷達(dá)數(shù)據(jù)的獲取方式，其掃描得到的電力走廊點(diǎn)云密度往往不夠均勻. 若點(diǎn)云在整個(gè)走廊上分布不均，則會(huì)導(dǎo)致統(tǒng)計(jì)的局部密度特征與地物的垂直分布特征失去相關(guān)性，從而無法準(zhǔn)確識(shí)別電塔，因此需要對(duì)電力走廊點(diǎn)云數(shù)據(jù)進(jìn)行密度均勻化處理. 本文采用了體素格網(wǎng)濾波器，其格網(wǎng)尺寸受提取的電塔尺寸和原數(shù)據(jù)的點(diǎn)云密度限制，而后者主要受無人機(jī)航高、航速和掃描角分辨率等因素共同影響.

為了進(jìn)一步提高處理效率，更有效提取點(diǎn)云特征，將經(jīng)過濾波之后的點(diǎn)云數(shù)據(jù)建立K- D樹索引，并進(jìn)行格網(wǎng)化，格網(wǎng)的尺寸大小與體素格網(wǎng)濾波器的格網(wǎng)尺寸、電力走廊區(qū)域的地形以及電塔的尺寸有關(guān). 若格網(wǎng)尺寸較大，則會(huì)影響定位精度；若格網(wǎng)區(qū)域過小，則會(huì)影響程序的處理效率且獲取不到足夠的統(tǒng)計(jì)數(shù)據(jù).

點(diǎn)云規(guī)則格網(wǎng)化的具體公式：

(1)

式中：(x,y)為待定點(diǎn)的二維平面坐標(biāo)，xmin、ymin分別為點(diǎn)云數(shù)據(jù)中x、y坐標(biāo)的最小值，g為格網(wǎng)尺寸(即格網(wǎng)邊長)，i、j分別為計(jì)算得到的待定點(diǎn)所屬區(qū)域的點(diǎn)云規(guī)則格網(wǎng)的行、列號(hào).

1.2 AOI提取方法

輸電走廊區(qū)域中含有電力設(shè)施(電力線、電塔等)的區(qū)域較其他區(qū)域有著較大的相對(duì)高度，因此可以根據(jù)nDSM快速得到電塔的范圍. 此外，電塔區(qū)域較其他區(qū)域又有著局部密度高的特點(diǎn). 因此，可以結(jié)合局部密度特征確定AOI的位置.

本文統(tǒng)計(jì)的具體點(diǎn)云特征包含點(diǎn)云局部密度特征，數(shù)字高程模型(Digital Eleration Model，DEM)特征和數(shù)字表面模型(Digital Surface Moder，DSM)特征[11]，并據(jù)此得到nDSM數(shù)據(jù).

點(diǎn)云格網(wǎng)的局部密度特征是定位電塔的重要依據(jù)，其計(jì)算方法是將點(diǎn)云降維，以一個(gè)格網(wǎng)為單位，統(tǒng)計(jì)落入格網(wǎng)中點(diǎn)的個(gè)數(shù)，以作為該格網(wǎng)區(qū)域內(nèi)的局部密度特征.

單個(gè)格網(wǎng)的DEM、DSM特征值的具體計(jì)算方法是統(tǒng)計(jì)其中最低點(diǎn)、最高點(diǎn)的高程值，并分別作為該格網(wǎng)區(qū)域的DEM、DSM特征值.

單個(gè)點(diǎn)云格網(wǎng)的nDSM特征數(shù)據(jù)是用其DSM特征值減去DEM特征值得到的，記錄了格網(wǎng)區(qū)域內(nèi)高于地面的地物相對(duì)于地面的最大高度信息. 該信息對(duì)于區(qū)分不同類型地物具有重要作用，能夠有效找到高大地物，如電塔、電力線、建筑物及高大喬木植被等[12].

計(jì)算出各個(gè)格網(wǎng)中的上述特征之后，可根據(jù)密度特征和nDSM特征分別設(shè)定相應(yīng)的閾值提取得到AOI. 電塔nDSM特征分離閾值H*的設(shè)立應(yīng)不高于數(shù)據(jù)區(qū)域中最低電塔的高度. 密度特征分離閾值D則受點(diǎn)云數(shù)據(jù)密度影響，其取值為經(jīng)驗(yàn)值. 根據(jù)上述H*和D可以提取AOI.

1.3 基于高程空間統(tǒng)計(jì)法的電塔存在區(qū)域提取

由于電力走廊區(qū)域環(huán)境復(fù)雜，不僅存在電塔，還存在電力線、復(fù)雜地形和高大植被. 這些都會(huì)對(duì)電塔的識(shí)別定位帶來影響. 例如，電力線具有較大nDSM值，而高大植被等兼具有nDSM值較大、密度較高的特點(diǎn). 為了在AOI中有效地提取出電塔，本文在保證便于計(jì)算的同時(shí)提取不同地物的高程空間特征并予以區(qū)分，將各個(gè)AOI的高程空間分為上下2個(gè)區(qū)域，如圖2所示. 各區(qū)域上下閾值以該區(qū)域nDSM值的特定倍數(shù)設(shè)定，目的是為避免復(fù)雜電力線噪聲或復(fù)雜植被噪聲等對(duì)電塔分類、定位的影響. 其具體過程如下：

第一，設(shè)置合理的閾值，將AOI的高程空間分為上部區(qū)域Cloudup和下部區(qū)域Clouddown，如圖2所示.

第二，分別統(tǒng)計(jì)上部區(qū)域和下部區(qū)域中的點(diǎn)云密度Dup、Ddown.

第三，計(jì)算各個(gè)AOI格網(wǎng)的上部區(qū)域和下部區(qū)域的點(diǎn)云密度之比r.

由于桿塔自身結(jié)構(gòu)特點(diǎn)，其形狀在高程空間內(nèi)分布較為均勻和連續(xù)，與其他高大地物(如喬木等)有一定區(qū)別. 因此，根據(jù)電塔點(diǎn)云數(shù)據(jù)在其高程空間內(nèi)的上、下部區(qū)域均有一定的密度、分布較為均勻的特征，可以在AOI內(nèi)的復(fù)雜環(huán)境中找出含有電塔的區(qū)域.

1.4 電塔擬合定位算法

由于經(jīng)過上述特征提取獲得的只是電塔存在的格網(wǎng)區(qū)域，而不是電塔實(shí)際的存在位置，因此需要對(duì)電塔存在區(qū)域進(jìn)行電塔位置的擬合. 為了獲得電塔的精確位置，本文使用了電塔位置粗定位和高程分層迭代精定位2個(gè)步驟，具體流程圖如圖3所示.

由于提取特征的格網(wǎng)尺寸小于電塔的尺寸，且同一電塔有著眾多不同的響應(yīng)區(qū)域，因此粗定位可以將經(jīng)過之前1.3提取得到的各個(gè)電塔存在區(qū)域的中心點(diǎn)取平均值，從而得到電塔大概位置，并將臨近的響應(yīng)區(qū)域識(shí)別為同一個(gè)電塔.

高程分層迭代精定位是在電塔位置粗定位基礎(chǔ)上進(jìn)行的. 當(dāng)粗定位得到每個(gè)電塔大概的位置坐標(biāo)后，以該坐標(biāo)為圓心，一定半徑R和一定搜索高度H來劃定圓柱形搜索區(qū)，在濾波后的點(diǎn)云數(shù)據(jù)中搜索得到在搜索區(qū)范圍內(nèi)的點(diǎn)，并擬合計(jì)算該點(diǎn)云質(zhì)心的x、y坐標(biāo)值作為電塔平面坐標(biāo).

其中，搜索高度H是根據(jù)固定系數(shù)乘以該區(qū)域nDSM值而確定的，因此并不能在劃分搜索高度時(shí)確定以H劃分的搜索區(qū)在電塔點(diǎn)云塔身的具體部位，搜索區(qū)中有可能包含有電力線點(diǎn)、絕緣子點(diǎn)、懸掛點(diǎn)等噪聲點(diǎn)，從而影響電塔擬合定位的精度. 為確定合適的擬合高度，本文采用了高程分層迭代的方法來選擇點(diǎn)云密度較低的搜索區(qū)進(jìn)行擬合. 當(dāng)搜索區(qū)點(diǎn)云密度大于設(shè)定的點(diǎn)云密度閾值Dn時(shí)，H則遞增Δh.

經(jīng)過高程分層迭代精定位后，可以獲得較高精度的電塔位置坐標(biāo)，從而完成基于無人機(jī)機(jī)載雷達(dá)數(shù)據(jù)的輸電線走廊電塔的識(shí)別與定位.

2 試驗(yàn)結(jié)果與分析

2.1 試驗(yàn)數(shù)據(jù)

所用試驗(yàn)數(shù)據(jù)包含了26 165 476個(gè)點(diǎn)，平均點(diǎn)密度47 pts/m2，電力走廊電壓等級(jí)為500 kV,電力走廊長度為4 728 m，地形復(fù)雜，包含平原和山地，植被茂盛. 此數(shù)據(jù)在提供時(shí)已經(jīng)經(jīng)過部分的預(yù)處理，包括配準(zhǔn)、裁剪等. 電力走廊部分有13座電塔，其中有1座電塔數(shù)據(jù)不完整. 電塔之間高度不統(tǒng)一，主要分布在45～100 m(圖4).

由于數(shù)據(jù)區(qū)內(nèi)存在交叉跨越塔架，且電塔數(shù)據(jù)不完整、掃描密度不均勻等，都會(huì)給電塔的自動(dòng)定位帶來了一定困難.

2.2 試驗(yàn)結(jié)果

本文使用C++語言在Windows平臺(tái)上利用Visual Studio 2015編譯器開發(fā)，集成了點(diǎn)云處理庫(Point Cloud Library,PCL)來實(shí)現(xiàn)算法的主要功能. 依據(jù)實(shí)驗(yàn)區(qū)的地形情況和電塔尺寸設(shè)置格網(wǎng)邊長為2 m. 在根據(jù)試驗(yàn)數(shù)據(jù)設(shè)置合適的密度分離閾值D和nDSM分離閾值H*后，便可根據(jù)H*、D得到試AOI. 本次試驗(yàn)中所設(shè)的密度分離閾值D為82.5 pts/m2，nDSM分離閾值H*設(shè)為25 m.

2.2.1 對(duì)數(shù)據(jù)進(jìn)行高程空間統(tǒng)計(jì)并提取電塔存在區(qū)

本文統(tǒng)計(jì)的高程空間點(diǎn)云分布特征為上、下區(qū)域點(diǎn)云密度特征和它們的比值，具體參見1.3.

其中，上、下部區(qū)域密度特征對(duì)電力桿塔較為敏感，可以在電塔位置AOI數(shù)據(jù)上提取出電塔所在的位置. 其中上部密度特征除了對(duì)電塔有著較強(qiáng)的響應(yīng)以外，還對(duì)電力線有著一定的相應(yīng). 而下部密度特征也會(huì)受到地表密集植被影響.

AOI上下密度比例特征則能較好地從AOI中提取電塔位置，其分離效果也較為可觀.

根據(jù)統(tǒng)計(jì)得到的AOI高程空間特征，分別設(shè)定上部區(qū)域和下部區(qū)域中的點(diǎn)云密度閾值Dup、Ddown及其比值r，從而精確提取得到電塔所在的區(qū)域. 本次試驗(yàn)中Dup、Ddown設(shè)為不小于10 pts/m2，半徑R則為0.5～0.8.

2.2.2 對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行電塔擬合定位

為了在各個(gè)電塔存在區(qū)域中精確定位得到電塔的位置，對(duì)電塔存在區(qū)域數(shù)據(jù)進(jìn)行了“先粗后細(xì)”的定位方式.

在進(jìn)行完粗定位后，可以在粗定位的成果之上進(jìn)行高程分層迭代精定位. 搜索高度H設(shè)為92% nDSM，搜索半徑R設(shè)為8 m. 經(jīng)過試驗(yàn)驗(yàn)證，進(jìn)行高程分層迭代會(huì)顯著提高電塔定位精度. 其最終定位結(jié)果如圖5所示.

2.3 討論與分析

經(jīng)過上述過程，可得到電塔平面坐標(biāo)，其定位精度分析結(jié)果見表1. 其中，人工真值X、Y為在試驗(yàn)點(diǎn)云上人工選取的各個(gè)電塔的中心坐標(biāo)，x、y則為本文算法所得到的定位坐標(biāo). 可以看到，本文算法能將試驗(yàn)數(shù)據(jù)中的電塔全部識(shí)別并進(jìn)行定位，證明了本算法在復(fù)雜環(huán)境中對(duì)非電塔目標(biāo)(如喬木等)有一定的抗干擾能力. 其中4號(hào)電塔由于原始掃描數(shù)據(jù)不全導(dǎo)致誤差較大(圖6)，其誤差達(dá)到了近2 m. 但這也反映出本算法對(duì)各種電塔識(shí)別的穩(wěn)健性，即使電塔數(shù)據(jù)不全也能進(jìn)行識(shí)別. 其余電塔中最大點(diǎn)位誤差為1.253 m,平均誤差為0.778 m. 由結(jié)果看出，本文算法精度達(dá)到米級(jí)，可以為電力走廊建模、點(diǎn)云精細(xì)分類以及災(zāi)害檢測(cè)等提供基礎(chǔ)初始數(shù)據(jù)，基本滿足電力巡檢數(shù)據(jù)處理、分析的要求.

本文為減少定位處理中所處理的數(shù)據(jù)量，提高算法運(yùn)行效率，在定位之前先提取AOI. 在進(jìn)行AOI的提取之后，可以排除原始數(shù)據(jù)中大量的無關(guān)點(diǎn)云(如植被點(diǎn)、地面點(diǎn)及部分電力線點(diǎn))，此時(shí)對(duì)提取后的數(shù)據(jù)進(jìn)行高程空間的統(tǒng)計(jì)則能極大減少處理的數(shù)據(jù)量，從而在保證定位精度的情況下提高算法的效率并減少資源消耗. 本文處理4.7 km數(shù)據(jù)共耗時(shí)136 s，大致處理速度達(dá)到了125 km/h.

為了提高效率，本文并沒有采用復(fù)雜的特征提取算法，如使用霍夫變換、懸鏈擬合、種子點(diǎn)生長算法以及計(jì)算坡度等復(fù)雜二維特征，但依舊獲得了較為良好的定位精度.

表1 試驗(yàn)數(shù)據(jù)定位精度分析

3 結(jié)論

作為一項(xiàng)自動(dòng)化程度較高的技術(shù)，無人機(jī)電力巡檢有著一定優(yōu)勢(shì). 本文針對(duì)電塔點(diǎn)云的nDSM、局部密度和特有高程空間分布特征，提出了一種由二維格網(wǎng)特征提取和高程空間特征分析相結(jié)合的電塔自動(dòng)定位方法. 本算法有如下特點(diǎn)：

1)降低了點(diǎn)云密度分布的均勻性對(duì)電塔定位精度的影響. 本文使用體素格網(wǎng)濾波器對(duì)數(shù)據(jù)進(jìn)行了處理，降低點(diǎn)云密度分布的不均勻性，并通過將迭代搜索的初始高程設(shè)置為較高的位置進(jìn)行改善. 根據(jù)無人機(jī)掃描的特點(diǎn)，在較高位置上其掃描分辨率較高，數(shù)據(jù)較近地面數(shù)據(jù)也更為全面和均勻. 同時(shí)，較高的位置也可以有效避免電塔位置擬合時(shí)受到地表植被的影響.

2)本文算法與其他算法相比，在結(jié)合有效提取手段的同時(shí)，簡(jiǎn)化了處理流程，并能為之后的點(diǎn)云精細(xì)分類、災(zāi)害檢測(cè)等提供基礎(chǔ)數(shù)據(jù).

3)本算法可針對(duì)不同類型的電塔設(shè)置特定的提取參數(shù)，有較好的普適性. 經(jīng)過試驗(yàn)驗(yàn)證，本方法能夠準(zhǔn)確提取并定位無人機(jī)機(jī)載激光雷達(dá)獲取的電力走廊點(diǎn)云數(shù)據(jù)中的電塔，擁有較好的穩(wěn)定性、適應(yīng)性和處理效率，對(duì)促進(jìn)和提高無人機(jī)巡檢的自動(dòng)化水平有著積極的作用，并為智能巡檢系統(tǒng)的開發(fā)帶來了一定的啟發(fā).