胡 舒，王樹(shù)根，王 越，李 欣

（1.武漢大學(xué) 遙感信息工程學(xué)院，湖北 武漢 430079；2.武漢市測(cè)繪研究院，湖北 武漢 430022）

隨著科技的發(fā)展，城市中私家車(chē)的數(shù)量急劇增加，違規(guī)占道、違章停車(chē)等情況層出不窮，導(dǎo)致道路安全形勢(shì)日趨嚴(yán)峻。夜間違停亂停導(dǎo)致的道路安全問(wèn)題越來(lái)越受到重視，公眾對(duì)公共停車(chē)場(chǎng)數(shù)量的需求也越來(lái)越大，因此改善當(dāng)前城市夜間道路環(huán)境已刻不容緩。為了獲取城市區(qū)域道路兩側(cè)的夜間停車(chē)信息，可利用夜間車(chē)流量相對(duì)較少、路面狀況較穩(wěn)定的特點(diǎn)開(kāi)展城市區(qū)域道路兩側(cè)停車(chē)調(diào)查工作，為建設(shè)公共停車(chē)場(chǎng)以及治理夜間違章停車(chē)提供依據(jù)。

車(chē)載激光掃描技術(shù)是一種動(dòng)態(tài)獲取城市三維空間數(shù)據(jù)的新技術(shù)[1]。該技術(shù)已在三維城市建設(shè)、桿狀地物檢測(cè)、建筑物立面檢測(cè)以及道路安全檢測(cè)等領(lǐng)域得到廣泛應(yīng)用。與機(jī)載LiDAR相比，車(chē)載LiDAR具有掃描頻率高、移動(dòng)速度快和測(cè)量距離近等特點(diǎn)，采集的點(diǎn)云數(shù)據(jù)量大且內(nèi)容豐富，包括道路沿線(xiàn)地物，如樹(shù)木、建筑物、行人、線(xiàn)桿、車(chē)輛以及其他目標(biāo)等。車(chē)載激光掃描技術(shù)不受光照因素的制約，可用來(lái)獲取夜間道路兩側(cè)的地物目標(biāo)點(diǎn)云數(shù)據(jù)。

針對(duì)基于三維激光點(diǎn)云的目標(biāo)檢測(cè)方法，國(guó)內(nèi)外已有不少研究，如譚賁[2]、PU S[3]等通過(guò)分析車(chē)載LiDAR點(diǎn)云中不同地物的掃描點(diǎn)特征來(lái)采集點(diǎn)云數(shù)據(jù)塊的各項(xiàng)特征，再進(jìn)行地物的自動(dòng)識(shí)別和分類(lèi)，但沒(méi)有討論車(chē)輛目標(biāo)的檢測(cè)問(wèn)題；YANG B S[4]等發(fā)現(xiàn)從點(diǎn)云數(shù)據(jù)中提取車(chē)輛目標(biāo)存在一定的困難，由于建筑物陰影、樹(shù)木等對(duì)車(chē)輛的遮擋，導(dǎo)致檢測(cè)準(zhǔn)確率降低。然而，基于車(chē)載LiDAR點(diǎn)云的車(chē)輛目標(biāo)自動(dòng)檢測(cè)問(wèn)題鮮有研究，因此本文根據(jù)車(chē)輛目標(biāo)的點(diǎn)云特征和自然空間特征等，提出了一種基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)方法，并通過(guò)實(shí)測(cè)數(shù)據(jù)驗(yàn)證了該算法的有效性。

1 基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)

1.1 點(diǎn)云預(yù)處理

車(chē)載激光掃描系統(tǒng)獲取的車(chē)載LiDAR點(diǎn)云數(shù)據(jù)量非常龐大，且存在數(shù)據(jù)冗余，為了便于后續(xù)處理，需對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行預(yù)處理。本文利用參考文獻(xiàn)[5]的方法分別在水平方向和垂直方向?qū)c(diǎn)云進(jìn)行處理，主要目的是去除地面點(diǎn)云。在水平方向上，根據(jù)道路寬度信息建立緩沖區(qū)，道路兩側(cè)以外的不屬于道路范圍的點(diǎn)云被全部去除；在垂直方向上，結(jié)合POS軌跡數(shù)據(jù)并按高度過(guò)濾，提取出地面到掃描中心高度范圍內(nèi)的點(diǎn)云。通過(guò)上述預(yù)處理，去除了地面點(diǎn)云，只保留了感興趣的道路范圍內(nèi)的點(diǎn)云數(shù)據(jù)，大大減少了點(diǎn)云數(shù)量，提高了點(diǎn)云的處理效率。

1.2 獨(dú)立地物分割

本文采用分割的方法對(duì)點(diǎn)云進(jìn)行聚類(lèi)分析，將點(diǎn)云分割為點(diǎn)云簇；再根據(jù)車(chē)輛目標(biāo)的特征，實(shí)現(xiàn)城市道路兩側(cè)車(chē)輛目標(biāo)的自動(dòng)檢測(cè)。聚類(lèi)算法主要包括層次聚類(lèi)方法、劃分方法、基于網(wǎng)格的方法、基于密度的方法以及基于模型的方法[6-8]。除此之外，還可將點(diǎn)云數(shù)據(jù)進(jìn)行投影，并對(duì)投影得到的圖像進(jìn)行二值化處理，再利用圖像的聚類(lèi)方法[9]進(jìn)行聚類(lèi)分析。

車(chē)載LiDAR點(diǎn)云剔除地面點(diǎn)后，路面目標(biāo)點(diǎn)云為獨(dú)立分布狀態(tài)，因此本文分別采用基于歐式距離分割的聚類(lèi)方法和基于密度的DBSCAN聚類(lèi)方法對(duì)車(chē)載LiDAR點(diǎn)云數(shù)據(jù)進(jìn)行聚類(lèi)，并比較兩種算法的性能，進(jìn)而選取最適合本文數(shù)據(jù)的方法進(jìn)行聚類(lèi)實(shí)驗(yàn)。

1.2.1 基于歐式距離分割的聚類(lèi)方法

該方法基于歐式距離度量點(diǎn)與點(diǎn)之間的遠(yuǎn)近[10]。在多維空間中，利用點(diǎn)與點(diǎn)之間的幾何距離來(lái)表示歐氏距離，如在三維空間中的任意兩點(diǎn)A(XA,YA,ZA)和B(XB,YB,ZB)，其歐式距離可表示為：

針對(duì)點(diǎn)云數(shù)據(jù)量龐大、冗余嚴(yán)重且處理效率低下的情況，在進(jìn)行基于歐式距離分割的聚類(lèi)之前，需先對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行體素柵格濾波，以實(shí)現(xiàn)體素網(wǎng)格化下采樣，該過(guò)程利用PCL庫(kù)[11]中的VoxelGrid濾波器實(shí)現(xiàn)。由于點(diǎn)云中不可避免的存在大量噪聲干擾，為了減少干擾，利用統(tǒng)計(jì)分析技術(shù)對(duì)點(diǎn)云進(jìn)行濾波，濾除噪聲點(diǎn)，可大大減少冗余的點(diǎn)云數(shù)據(jù)，使點(diǎn)云形狀特征得以保持，該過(guò)程利用PCL庫(kù)中的StatisticalOutlierRemoval濾波器實(shí)現(xiàn)。

在進(jìn)行基于歐式距離分割的聚類(lèi)過(guò)程中，利用KD Tree算法在多維空間內(nèi)對(duì)每個(gè)點(diǎn)的最鄰近點(diǎn)進(jìn)行搜索，小于設(shè)定的距離閾值D的點(diǎn)云數(shù)據(jù)被劃分到同一個(gè)點(diǎn)云簇。設(shè)定一個(gè)點(diǎn)云簇包含點(diǎn)云數(shù)量的最小值（Min_Size）和最大值（Max_Size），只有包含點(diǎn)云數(shù)量在這個(gè)范圍之間的點(diǎn)云簇才能被保存下來(lái)，其余點(diǎn)云簇則直接濾除，該閾值設(shè)定的目的是濾除分割得到的細(xì)小部件。

1.2.2 基于密度的DBSCAN聚類(lèi)方法

由于濾除地面點(diǎn)云后，點(diǎn)云密度特征明顯，因此本文同時(shí)采用基于密度的方法進(jìn)行聚類(lèi)。具有噪聲的基于密度的DBSCAN[12]聚類(lèi)方法是一個(gè)典型的基于密度的聚類(lèi)方法。其核心思想為：在點(diǎn)云數(shù)據(jù)集中選擇任意一個(gè)核心對(duì)象（即在一定范圍內(nèi)周?chē)c(diǎn)的數(shù)量大于一定閾值的點(diǎn)）作為種子點(diǎn)，然后向周?chē)鷶U(kuò)展，找到該種子點(diǎn)密度可達(dá)的所有樣本集合，聚類(lèi)為一個(gè)點(diǎn)云簇，重復(fù)上述過(guò)程，直到點(diǎn)云數(shù)據(jù)集內(nèi)所有的點(diǎn)都被聚類(lèi)或被判定為異常點(diǎn)為止。

在聚類(lèi)過(guò)程中，本文采用基于歐式距離分割的聚類(lèi)方法和基于密度的DBSCAN聚類(lèi)方法分別對(duì)同一組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行聚類(lèi)實(shí)驗(yàn)。聚類(lèi)結(jié)果如圖1所示，可以看出，圖1a中下采樣效果以及離群點(diǎn)濾除效果明顯，點(diǎn)云簇邊緣更加光滑；圖1b中點(diǎn)云簇邊緣噪聲點(diǎn)依然存在，且離群點(diǎn)去除效果較差?；跉W式距離分割的聚類(lèi)方法的聚類(lèi)結(jié)果明顯優(yōu)于基于密度的DBSCAN聚類(lèi)方法，如圖1中紅色矩形框所示，基于密度的DBSCAN聚類(lèi)方法分割嚴(yán)重不足，這是由于該點(diǎn)云數(shù)據(jù)集密度不均勻所導(dǎo)致的，因此本文選擇基于歐式距離分割的方法進(jìn)行聚類(lèi)。

圖1 兩種聚類(lèi)方法效果對(duì)比

1.3 車(chē)輛目標(biāo)的自動(dòng)檢測(cè)

1.3.1 特征提取

對(duì)分割后得到的獨(dú)立地物點(diǎn)云簇進(jìn)行點(diǎn)云特征的提取，提取的特征應(yīng)能區(qū)別車(chē)輛與其他點(diǎn)云數(shù)據(jù)中的地物。根據(jù)車(chē)輛的空間拓?fù)湮恢?、物理大小、?chē)輛與地面的關(guān)系以及車(chē)輛的幾何特征等，建立車(chē)輛目標(biāo)自動(dòng)檢測(cè)的知識(shí)規(guī)則，再根據(jù)該規(guī)則實(shí)現(xiàn)車(chē)輛的自動(dòng)檢測(cè)。知識(shí)規(guī)則涉及車(chē)輛目標(biāo)的大小、高度、水平投影面積、水平投影外接矩形邊界框的長(zhǎng)與寬等。

1）大小是區(qū)分不同地物目標(biāo)的重要特征之一。本文采用每個(gè)獨(dú)立地物包含的點(diǎn)云數(shù)量來(lái)表示地物大小。

2）高度是區(qū)分城市地物的重要特征之一。首先計(jì)算獨(dú)立地物的最大高度Zmax和最小高度Zmin，再根據(jù)式（2）計(jì)算得到獨(dú)立地物的高度H。車(chē)輛目標(biāo)的平均高度是相似的，約為1.5 m。

3）水平投影面積。車(chē)輛在三維、二維空間中的形態(tài)如圖2所示，城市中不同型號(hào)車(chē)輛在水平面上的投影在面積和形狀上是相似的，因此地物的水平投影面積S，即為物體點(diǎn)云簇在水平面上投影的邊界矩形框面積，計(jì)算公式為：

式中，xmax、xmin、ymax、ymin分別為x坐標(biāo)和y坐標(biāo)的最大值與最小值。

4）水平投影外接矩形邊界框的長(zhǎng)與寬，描述的是物體的形狀特性。車(chē)輛具有相似的形狀，其水平投影外接矩形邊界框的長(zhǎng)與寬也相似，邊長(zhǎng)過(guò)長(zhǎng)和過(guò)短的地物都不符合車(chē)輛的形狀特性。

圖2 車(chē)輛在三維、二維空間的示意圖

1.3.2 特征判定與車(chē)輛檢測(cè)

在圖像處理時(shí)，通常采用模板匹配的方法來(lái)檢測(cè)地物目標(biāo)，但在三維點(diǎn)云數(shù)據(jù)處理時(shí)，該方法不適用，因此本文基于上述知識(shí)規(guī)則和要素特征，分別計(jì)算各獨(dú)立地物點(diǎn)云的屬性，為進(jìn)行車(chē)輛目標(biāo)的檢測(cè)提供基礎(chǔ)知識(shí)。

根據(jù)知識(shí)規(guī)則定義計(jì)算獨(dú)立地物的特征向量Vi=[x,y,z,n,h,s,l]，其中x、y、z為空間位置，以獨(dú)立地物的質(zhì)心表示；n為大小，以獨(dú)立地物點(diǎn)云中所包含的點(diǎn)云數(shù)量來(lái)表示；h為高度；s為水平投影面積；l為水平投影外接矩形邊界框的長(zhǎng)與寬。結(jié)合車(chē)輛目標(biāo)的自然特征和空間幾何特征，與知識(shí)規(guī)則建立自然對(duì)應(yīng)關(guān)系，計(jì)算各獨(dú)立地物點(diǎn)云的特征，并得到特征向量；再與建立的車(chē)輛目標(biāo)自動(dòng)檢測(cè)知識(shí)規(guī)則進(jìn)行匹配，實(shí)現(xiàn)車(chē)輛目標(biāo)的自動(dòng)檢測(cè)。

基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)算法流程如圖3所示。

圖3 基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)算法流程圖

2 實(shí)驗(yàn)結(jié)果與分析

本文選取武漢市三環(huán)線(xiàn)內(nèi)某道路的車(chē)載LiDAR點(diǎn)云實(shí)測(cè)數(shù)據(jù)進(jìn)行車(chē)輛目標(biāo)檢測(cè)實(shí)驗(yàn)，以驗(yàn)證提出的基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)算法的有效性與可行性。

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

實(shí)驗(yàn)數(shù)據(jù)來(lái)源于武漢市三環(huán)線(xiàn)內(nèi)主要道路的夜間停車(chē)狀況調(diào)查項(xiàng)目。武漢市三環(huán)線(xiàn)內(nèi)某段道路預(yù)處理后的點(diǎn)云數(shù)據(jù)如圖4所示，呈長(zhǎng)條狀分布，全程長(zhǎng)約為1.4 km，共包含9 742 560個(gè)掃描點(diǎn)。

圖4 武漢市三環(huán)線(xiàn)內(nèi)某段道路的車(chē)載LiDAR點(diǎn)云數(shù)據(jù)

2.2 實(shí)驗(yàn)結(jié)果與分析

考慮到點(diǎn)云數(shù)據(jù)量龐大、處理效率低的情況，實(shí)驗(yàn)中將該道路點(diǎn)云分為8個(gè)小段，分別進(jìn)行車(chē)輛目標(biāo)的自動(dòng)檢測(cè)。首先對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行預(yù)處理、下采樣、濾波以及分割等處理；再利用基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)算法計(jì)算地物目標(biāo)的點(diǎn)云特征、幾何特征以及各項(xiàng)屬性；然后將其與車(chē)輛目標(biāo)知識(shí)規(guī)則庫(kù)進(jìn)行匹配，得到車(chē)輛目標(biāo)的檢測(cè)結(jié)果。某段道路的檢測(cè)結(jié)果如圖5所示，原始點(diǎn)云數(shù)據(jù)在FugroViewer中顯示，其余結(jié)果均在CloudCompare軟件中顯示。

由圖5b可知，基于歐式距離分割的聚類(lèi)方法能很好地去除離群點(diǎn)和異常點(diǎn)，并將地物目標(biāo)分割成獨(dú)立地物點(diǎn)云，地物目標(biāo)的分割效果較好；由圖5c可知，基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)算法能有效實(shí)現(xiàn)車(chē)輛目標(biāo)的自動(dòng)檢測(cè)。

為了驗(yàn)證本文算法的正確性，實(shí)驗(yàn)過(guò)程中對(duì)8個(gè)路段分別進(jìn)行了車(chē)輛目標(biāo)的自動(dòng)檢測(cè)，同時(shí)采用人工目視判讀的方法對(duì)點(diǎn)云數(shù)據(jù)中的車(chē)輛目標(biāo)進(jìn)行統(tǒng)計(jì)，并將本文算法得到的檢測(cè)結(jié)果與之進(jìn)行比較，結(jié)果如表1所示。

圖5 某段道路車(chē)輛目標(biāo)檢測(cè)結(jié)果

表1 樣本路段車(chē)輛檢測(cè)結(jié)果

由表1可知，基于知識(shí)規(guī)則的車(chē)輛目標(biāo)自動(dòng)檢測(cè)算法能有效檢測(cè)車(chē)輛目標(biāo)，自動(dòng)檢測(cè)準(zhǔn)確率可達(dá)88.10%，證明該方法是可行且有效的。

該路段某個(gè)位置造成的車(chē)輛漏檢（紅色框）和錯(cuò)檢（黃色框）現(xiàn)象如圖6所示。通過(guò)分析可以發(fā)現(xiàn)，影響自動(dòng)檢測(cè)準(zhǔn)確率的主要原因一方面是部分車(chē)輛目標(biāo)之間存在相互遮擋，導(dǎo)致距離過(guò)近、分割不足，或是車(chē)輛點(diǎn)云不連貫導(dǎo)致過(guò)分割、車(chē)輛信息不完整，進(jìn)而導(dǎo)致提取的特征無(wú)法與知識(shí)規(guī)則匹配；另一方面是由于極少數(shù)地物與車(chē)輛目標(biāo)有相似的特征，導(dǎo)致識(shí)別困難。

考慮到過(guò)分割與分割不足對(duì)車(chē)輛目標(biāo)檢測(cè)精度的影響，可通過(guò)優(yōu)化聚類(lèi)方法達(dá)到對(duì)點(diǎn)云數(shù)據(jù)的更優(yōu)分割，從而提高車(chē)輛目標(biāo)檢測(cè)的準(zhǔn)確率。除此之外，點(diǎn)云密度也是影響檢測(cè)精度的因素之一，因此在點(diǎn)云數(shù)據(jù)的獲取上也應(yīng)更加嚴(yán)謹(jǐn)。

圖6 車(chē)輛目標(biāo)漏檢結(jié)果示意圖

3 結(jié) 語(yǔ)

本文以武漢市三環(huán)線(xiàn)內(nèi)主要道路的夜間停車(chē)狀況調(diào)查項(xiàng)目為依托，采用武漢市三環(huán)線(xiàn)內(nèi)某段道路沿線(xiàn)點(diǎn)云數(shù)據(jù)進(jìn)行車(chē)輛目標(biāo)的自動(dòng)檢測(cè)。首先通過(guò)點(diǎn)云預(yù)處理移除地面點(diǎn)云，并對(duì)非地面點(diǎn)云進(jìn)行下采樣和濾波，以簡(jiǎn)化點(diǎn)云并去除離群點(diǎn)與異常點(diǎn)；然后采用基于歐式距離分割的聚類(lèi)方法將點(diǎn)云數(shù)據(jù)劃分為點(diǎn)云簇；再分析車(chē)輛目標(biāo)的點(diǎn)云特征、自然物理特征、空間幾何特征等相關(guān)特點(diǎn)，建立車(chē)輛目標(biāo)檢測(cè)的知識(shí)規(guī)則；最后根據(jù)知識(shí)規(guī)則實(shí)現(xiàn)車(chē)輛目標(biāo)的自動(dòng)檢測(cè)。實(shí)驗(yàn)表明，該方法的準(zhǔn)確率可達(dá)88.10%，減少了人工統(tǒng)計(jì)的工作量。在后續(xù)研究中，可進(jìn)一步改進(jìn)聚類(lèi)算法以及深入挖掘人工智能方法，以便提高車(chē)輛目標(biāo)檢測(cè)的精度和自動(dòng)化程度。