[摘 要]安全擋墻作為排土場安全措施之一，用來保障排土場正常作業(yè)，避免排巖車輛未保持安全卸載距離導致車輛滑落。文章提出基于無人機點云數(shù)據(jù)的排土場安全擋墻提取方法。以遼寧省鞍山市齊大山鐵礦排土場點云數(shù)據(jù)作為數(shù)據(jù)源，利用坡度濾波算法提取排土平臺點云數(shù)據(jù)，再對其進行alphashape邊緣提取算法來提取出安全擋墻一側(cè)的邊緣點云數(shù)據(jù)，再根據(jù)高程約束與向量叉乘位置判別來提取安全擋墻點云數(shù)據(jù)。通過對安全擋墻進行三維可視化建模實現(xiàn)安全擋墻的監(jiān)測。

[關(guān)鍵詞]無人機點云；安全擋墻；坡度濾波算法；alphashape邊緣提取算法

[中圖分類號]P231 [文獻標志碼]A [文章編號]2095–6487（2024）06–0168–03

Research on the Extraction Method of Safety Retaining Walls in Waste Disposal Sites Based on Point Cloud Data

YANG Wei，GUO Hangshen，WU Yixuan，HUAN Liu，MAO Yachun

[Abstract]As one of the safety measures for the dumping site， the safety retaining wall is used to ensure the normal operation of the dumping site and avoid the sliding of vehicles due to the failure of rock dumping vehicles to maintain a safe unloading distance. This article proposes a method for extracting safety retaining walls in waste disposal sites based on unmanned aerial vehicle point cloud data. Using the point cloud data of Qidashan Iron Mine dumping site in Anshan City， Liaoning Province， China as the data source， the slope filtering algorithm is first used to extract the point cloud data of the dumping platform. Then， the alphashape edge extraction algorithm is applied to extract the edge point cloud data on one side of the safety retaining wall. Finally， the safety retaining wall point cloud data is extracted based on height constraints and vector cross product position discrimination. Monitoring of safety barriers is achieved through three-dimensional visualization modeling of safety barriers.

[Keywords]drone point cloud； safety retaining wall； slope filtering algorithm； alphashape edge extraction algorithm

隨著礦山的不斷開采，排土場的高度持續(xù)增高，出現(xiàn)了大量高陡排土場邊坡。為了保障排巖車輛在進行卸載時保持安全距離，避免在排土場邊緣發(fā)生滑落，有必要設立與監(jiān)測排土場安全擋墻。

傳統(tǒng)排土場安全擋墻的安全巡檢多以人工到現(xiàn)場進行目視檢查為主。目前，已有不少國內(nèi)外學者對激光點云數(shù)據(jù)的濾波方法進行了研究。在機載方面，文獻[1]提出了自適應TIN三角網(wǎng)加密方法處理城市區(qū)域機載數(shù)據(jù)。該方法根據(jù)少量的種子點建立一個稀疏的TIN，然后通過迭代方法依據(jù)參數(shù)閾值對TIN進行不斷加密。文獻[2]提出了一種基于坡度變化的濾波算法，算法的關(guān)鍵在于選擇合適的坡度閾值，認為應當根據(jù)實驗區(qū)域地形的先驗知識選擇坡度閾值，由于需要包含所有地面形態(tài)的樣本，增大了算法的難度。文獻[3]提出利用不同尺寸窗口進行多次數(shù)學形態(tài)學濾波的方法，從最小窗口尺寸開始，對每個判斷為地面點的數(shù)據(jù)根據(jù)窗口尺寸的大小賦予權(quán)值，逐步增加窗口尺寸，迭代處理，窗口越大，權(quán)值越大，最終根據(jù)所有判斷出的地面點的權(quán)重計算地形表面。

安全擋墻的特征不明顯，導致常規(guī)地面濾波方法無法實現(xiàn)安全擋墻的提取，為此文章通過提取安全擋墻與排土場平臺之間的界線，再通過高程約束與向量叉乘位置判別來提取安全擋墻點云數(shù)據(jù)。

1 數(shù)據(jù)獲取及預處理

1.1 實驗數(shù)據(jù)

選擇遼寧省鞍山市齊大山鐵礦排土場，并以大疆M300 RTK無人機搭載激光雷達禪思L1獲取的排土場點云數(shù)據(jù)作為數(shù)據(jù)源開展研究。排土場位置無人機正射影像如圖1所示。

1.2 數(shù)據(jù)預處理

排土場點云數(shù)據(jù)存在航帶重疊區(qū)域，導致重疊區(qū)域的點云密度遠大于非重疊區(qū)域，為了后續(xù)算法能獲得良好效果，利用格網(wǎng)下采樣算法對排土場點云進行預處理，以此實現(xiàn)點云數(shù)據(jù)的下采樣，改善點云數(shù)據(jù)密度不均的問題。

2 研究方法

2.1 排土平臺點云提取

排土場卸載區(qū)的平臺區(qū)域地形平坦，該區(qū)域點云點與點之間的坡度變化較小，而安全擋墻近似梯形或三角形，使得安全擋墻點云點與點之間坡度變化較大。因此，需使用坡度濾波算法實現(xiàn)對排土場平臺與邊坡點云的提取。具體步驟如下。

（1）使用坡度濾波來提取出高程梯度較小的點云。

（2）給每個點云賦予對應的格網(wǎng)索引。

（3）確定劃分的每個格網(wǎng)中z坐標最小的點，格網(wǎng)中的其他點分別與該點進行相對高程及坡度計算。

（4）根據(jù)設定的高程閾值H及坡度閾值G，將小于高程閾值H和坡度閾值G的點提取出來，實現(xiàn)平坦區(qū)域點云提取，將大于高程閾值H和坡度閾值G的點提取出來，實現(xiàn)陡峭區(qū)域點云提取。

2.2 排土界線特征點提取

排土場平臺與安全擋墻的分界點提取是間接提取安全擋墻的重要步驟。為了成功提取安全擋墻，需要實現(xiàn)排土場平臺靠近安全擋墻那一側(cè)的邊緣點云的提取，為此在alpha shape邊緣檢測算法的基礎(chǔ)上，通過歐式聚類[4]提取出來的安全擋墻頂部點云數(shù)據(jù)來對邊緣點云進行最鄰近約束，實現(xiàn)排土界線的提取。

首先找到邊界點云數(shù)據(jù)的質(zhì)心點，計算每個點與該點連線到x軸正方向的夾角，排序則為排序結(jié)果為逆時針，夾角從大到小的順序排序，則排序結(jié)果為順時針。根據(jù)夾角的大小進行排序。按照夾角從小到大的順序再通過遍歷邊坡點云對有序邊界點進行最近鄰搜索，在邊坡點云的約束下實現(xiàn)有序排土界線特征點的提取如圖2所示。

2.3 排土場安全擋墻提取

安全擋墻完全在上述線狀點云的一側(cè)，為了避免提取安全擋墻時，提取出線狀點云另一側(cè)的非安全擋墻點云，需要判定安全擋墻點云與線狀點云的位置關(guān)系，來實現(xiàn)安全擋墻的提取。

（1）在工程實際應用中，NURBS（非均勻有理B樣條）算法可以精確描述空間曲線的姿態(tài)，為此本研究引入NURBS算法來擬合優(yōu)化臺階特征線。該方法主要通過調(diào)整和控制頂點來實現(xiàn)線條擬合。

（2）文章基于條件約束的安全擋墻提取算法實現(xiàn)安全擋墻的提取。

其中條件約束主要從兩部分來約束，具體如下。

（1）基于半徑鄰域搜索的高程約束部分，通過遍歷提取出來的排土場平臺與安全擋墻交界處的線狀點云中的每個點，以每個點為中心點，對整個排土場點云數(shù)據(jù)進行半徑搜索，并對搜索出來的點與對應的中心點的高程進行比較，通過設定的高程約束條件實現(xiàn)安全擋墻的初步提取。

（2）基于矢量叉乘來判定安全擋墻點與非安全擋墻點，實現(xiàn)安全擋墻點云數(shù)據(jù)的細化提取。通過遍歷上述從小到大排好序的線狀點云數(shù)據(jù)，計算兩兩相鄰點XOY平面的向量，并以線狀點云為中心，對初步提取的安全擋墻點云數(shù)據(jù)進行半徑搜索，計算搜索中心點與搜索點的向量，進行叉乘。若上述排序為順時針，只提取叉乘結(jié)果大于0的點，即位于線狀點云左側(cè)的安全擋墻點。若上述排序為逆時針，只提取叉乘結(jié)果小于0的點，即位于線狀點云右側(cè)的安全擋墻點。至此，實現(xiàn)了安全擋墻的提?。▓D3）。

3 安全擋墻提取精度評價

為了評估安全擋墻點集提取的精度，參考國際攝影測量和遙感學會提出的統(tǒng)計學誤差評判標準[5]，對安全擋墻提取結(jié)果的誤分率分為3種：第Ⅰ類誤差、第Ⅱ類誤差及總誤差。第Ⅰ類誤差、第Ⅱ類誤差反映了方法的適應性，總誤差反映了方法的可行性，總誤差越小則表示安全擋墻點集提取結(jié)果越精確。第Ⅰ類誤差、第Ⅱ類誤差、總誤差的定義分別為：

利用上述評估方法對本研究排土場點云安全擋墻提取結(jié)果進行誤差評定，具體結(jié)果見表1。

4 結(jié)束語

使用坡度濾波及alpha-shape等一系列算法提取排土界線，進而提取安全擋墻，該提取方法解決了傳統(tǒng)濾波方法提取安全擋墻出現(xiàn)誤提、漏提的情況。研究結(jié)果表明，基于排土界線進行范圍約束的安全擋墻提取算法能夠準確、自動地提取出排土場安全擋墻，可為安全擋墻健康評估提供技術(shù)支持，對露天礦智慧礦山建設具有重要的現(xiàn)實意義。

參考文獻

[1] 朱曉強.基于LiDAR點云和航空影像的城市三維重建[D].合肥：合肥工業(yè)大學，2009.

[2] 趙明君，劉超，高翔，等.基于多尺度網(wǎng)格的點云自適應坡度濾波算法[J].合肥工業(yè)大學學報（自然科學版），2021，44（2）：266-272.

[3] 張瑩，張凱義.數(shù)學形態(tài)學改進濾波方法在車載激光點云濾波中的應用[J].測繪標準化，2022，38（4）：81-85.

[4] 沈長表，夏永華，張坤祥，等.基于歐式聚類的車載點云道路提取[J].城市勘測，2023（3）：76-81.

[5] 毛亞純，伏雨文，曹旺，等.基于無人機點云數(shù)據(jù)的露天采場礦車提取方法[J].東北大學學報（自然科學版），2021，42（6）：842-848，863.