應(yīng)用于無(wú)人駕駛車的激光雷達(dá)雪天去噪方法研究

鐘科娣,劉元盛,張 軍,路 銘

(1.北京聯(lián)合大學(xué) 北京市信息服務(wù)工程重點(diǎn)實(shí)驗(yàn)室,北京 100101；2.北京聯(lián)合大學(xué)機(jī)器人學(xué)院,北京 100101；3.北京聯(lián)合大學(xué)應(yīng)用科技學(xué)院,北京 100101)

1 引 言

近年來(lái),無(wú)人駕駛技術(shù)的快速發(fā)展促進(jìn)了傳感器技術(shù)的發(fā)展。激光雷達(dá)傳感器以其高精度和高可靠性的特點(diǎn)而廣泛應(yīng)用于無(wú)人駕駛車輛中,其工作原理是通過(guò)旋轉(zhuǎn)發(fā)射激光束,掃描周邊三維空間范圍內(nèi)的物體實(shí)際位置,并將反射點(diǎn)以點(diǎn)云的信息存儲(chǔ)起來(lái)。通過(guò)分析激光雷達(dá)的點(diǎn)云數(shù)據(jù),可實(shí)現(xiàn)障礙物檢測(cè)、位姿估計(jì)、點(diǎn)云地圖建立等功能,為無(wú)人駕駛車輛提供必要的感知信息。由于激光雷達(dá)自身的缺陷或外部環(huán)境因素的影響,激光雷達(dá)的點(diǎn)云數(shù)據(jù)會(huì)包含一些不可忽視的噪聲,從而影響到激光雷達(dá)的檢測(cè)性能[1-3]。為了消除噪聲造成的不利影響,研究人員通常采用濾波算法對(duì)激光雷達(dá)數(shù)據(jù)進(jìn)行預(yù)處理,達(dá)到濾除噪聲的目的。針對(duì)點(diǎn)云噪聲的分布特點(diǎn),前人已經(jīng)研究出許多有效的去噪算法。這些去噪算法按空間分布可分為二維空間去噪算法和三維空間去噪算法。

二維空間去噪算法的應(yīng)用范圍為相機(jī)圖像的噪聲去噪,處理目標(biāo)為二維空間數(shù)據(jù),如用于處理邊緣細(xì)節(jié)和噪聲的濾波方法[4-5],以及通過(guò)自適應(yīng)形態(tài)學(xué)濾波抑制噪聲和雜波信號(hào)的檢測(cè)方法[6],而激光雷達(dá)點(diǎn)云是三維點(diǎn)云,在使用二維去噪算法進(jìn)行去噪前需將三維點(diǎn)云轉(zhuǎn)換到二維空間后再進(jìn)行去噪處理。常規(guī)的轉(zhuǎn)換方法是將一幀激光雷達(dá)點(diǎn)云投影到圖像上,其中激光掃描點(diǎn)代表圖像像素,點(diǎn)到傳感器的歐氏距離代表像素值,將激光雷達(dá)點(diǎn)云數(shù)據(jù)圖像化后再通過(guò)二維空間去噪方法進(jìn)行降噪處理[7-8]。

三維空間去噪算法可直接處理激光雷達(dá)點(diǎn)云,進(jìn)行點(diǎn)云去噪處理。Jenke等人提出了基于貝葉斯曲面重建的方法[9],通過(guò)貝葉斯規(guī)則創(chuàng)建一個(gè)分布模型去估計(jì)點(diǎn)云的分布和重建,在降噪的同時(shí)保持了點(diǎn)云的分布和結(jié)構(gòu)。敖建鋒等人有針對(duì)性的應(yīng)用多種濾波方法和格網(wǎng)插值方法提取信息,獲取高質(zhì)量的基礎(chǔ)數(shù)據(jù)[10]。張芳菲等人針對(duì)孤立噪聲點(diǎn),使用k-d樹(shù)處理點(diǎn)云數(shù)據(jù)得到有序點(diǎn)云,統(tǒng)計(jì)噪聲點(diǎn)與近鄰點(diǎn)的正態(tài)分布特性,自適應(yīng)計(jì)算去噪的閾值,對(duì)點(diǎn)云中孤立的噪聲點(diǎn)具有較好的去噪效果[11]。

常規(guī)的點(diǎn)云噪聲都有相對(duì)應(yīng)的降噪算法處理并且都具有較好的去噪效果,但對(duì)于一些特殊場(chǎng)景中的噪聲,常規(guī)的去噪算法的去噪效果較差。如在降雪場(chǎng)景中,雪花飄落時(shí)充斥在整個(gè)三維空間內(nèi),其反射的點(diǎn)云數(shù)據(jù)極易造成誤檢,故激光雷達(dá)點(diǎn)云必須要經(jīng)過(guò)去噪處理才能可靠的使用。雪花噪聲在激光雷達(dá)的近距離處密度大,而隨著距離的增加,雪花點(diǎn)的密度也逐漸降低,這些特性使得雪花具有不同于傳統(tǒng)點(diǎn)云噪聲的獨(dú)特噪聲分布,需要專門的去噪方法進(jìn)行濾除。

為了更好地解決雪花噪聲影響的問(wèn)題,本文采用了點(diǎn)云反射強(qiáng)度和鄰域搜索組合的去噪方法。反射強(qiáng)度是激光雷達(dá)掃描到物體時(shí)返回的脈沖回波強(qiáng)度的測(cè)量指標(biāo),并且在常規(guī)場(chǎng)景中不同的物體具有不同的反射強(qiáng)度,可根據(jù)反射強(qiáng)度實(shí)現(xiàn)物體的識(shí)別分類。但僅通過(guò)反射強(qiáng)度對(duì)物體進(jìn)行識(shí)別分類則會(huì)出現(xiàn)漏檢誤檢的情況,比如在目標(biāo)物體與其他物體具有相同的反射強(qiáng)度時(shí)不能精確識(shí)別分類?；卩徲蛩阉鞯姆椒赏ㄟ^(guò)定義鄰域,驗(yàn)證鄰域內(nèi)點(diǎn)云的分布屬性,達(dá)到識(shí)別分類的目的。但鄰域搜索的方法在雪天去噪的時(shí)不能保證環(huán)境特征的質(zhì)量。為此提出了基于反射強(qiáng)度和鄰域搜索的組合方法,在提高去噪效果的同時(shí),減少環(huán)境特征的損傷,為無(wú)人駕駛車輛的安全行駛提供優(yōu)質(zhì)的環(huán)境特征。

2 點(diǎn)云數(shù)據(jù)識(shí)別分類方法描述

2.1 基于激光雷達(dá)反射強(qiáng)度的識(shí)別分類

激光雷達(dá)數(shù)據(jù)包含了物體在三維空間的位置信息,也包含了目標(biāo)物體的反射強(qiáng)度等信息[12]。反射強(qiáng)度是反映目標(biāo)物體的重要的物理量,根據(jù)不同材質(zhì)物體有不同的反射強(qiáng)度,激光雷達(dá)的反射強(qiáng)度可用于區(qū)域劃分和物體識(shí)別分類[13]。激光雷達(dá)掃描到物體時(shí)返回的反射強(qiáng)度范圍為1～100。在實(shí)驗(yàn)中為了獲得雪花反射強(qiáng)度和環(huán)境特征反射強(qiáng)度的分布模型,對(duì)單幀點(diǎn)云數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。點(diǎn)云反射強(qiáng)度的分布模型可分為兩個(gè)步驟:第一步是對(duì)純雪花點(diǎn)云反射強(qiáng)度分析,通過(guò)提取雪花點(diǎn)云數(shù)據(jù)統(tǒng)計(jì)分析,得到雪花點(diǎn)云的反射強(qiáng)度分布模型。第二步是對(duì)環(huán)境特征的反射強(qiáng)度分析,通過(guò)晴天的激光雷達(dá)數(shù)據(jù)統(tǒng)計(jì)分析,得到環(huán)境特征的反射強(qiáng)度分布模型。在降雪場(chǎng)景中,積雪覆蓋會(huì)對(duì)環(huán)境特征的反射強(qiáng)度分布模型造成影響,所得到的強(qiáng)度分布模型不能準(zhǔn)確反映環(huán)境特征物的強(qiáng)度分布。故在統(tǒng)計(jì)環(huán)境特征物體的強(qiáng)度時(shí)采用晴天采集的激光雷達(dá)數(shù)據(jù)。針對(duì)晴天和降雪天氣場(chǎng)景,從采集的數(shù)據(jù)中各隨機(jī)抽取20幀數(shù)據(jù)統(tǒng)計(jì)分析,得到的激光雷達(dá)數(shù)據(jù)的強(qiáng)度分布,如圖1所示。通過(guò)驗(yàn)證可知雪花的反射強(qiáng)度值在1～8的范圍內(nèi),環(huán)境特征物體的反射強(qiáng)度為2～50,環(huán)境特征的反射強(qiáng)度在低強(qiáng)度處分布較少,當(dāng)強(qiáng)度大于45時(shí),統(tǒng)計(jì)到的激光點(diǎn)數(shù)較少,可根據(jù)得到的反射強(qiáng)度分布對(duì)物體進(jìn)行粗識(shí)別分類。若在降雪場(chǎng)景中僅通過(guò)反射強(qiáng)度對(duì)物體識(shí)別分類,會(huì)受到積雪覆蓋或部分與雪花具有相同反射強(qiáng)度環(huán)境特征的影響,在識(shí)別分類上具有單一屬性不可靠的局限性。

圖1 點(diǎn)云反射強(qiáng)度分布Fig.1 Reflection intensity distribution of point cloud

2.2 基于鄰域搜索算法的識(shí)別分類

鄰域搜索算法的本質(zhì)是通過(guò)迭代搜索目標(biāo)點(diǎn)的鄰域,獲得局部最優(yōu)解。鄰域搜索的關(guān)鍵是鄰域的確定,即如何定義鄰域。鄰域搜索中的動(dòng)態(tài)半徑離群點(diǎn)濾波器(DROR,Dynamic Radius Outlier Removal Filter)主要基于空間分布的規(guī)律進(jìn)行識(shí)別分類,通過(guò)確定搜索半徑來(lái)定義鄰域[14]。搜索半徑通過(guò)聯(lián)合目標(biāo)點(diǎn)到傳感器的歐式距離與激光雷達(dá)的水平分辨率的關(guān)系獲取。以目標(biāo)點(diǎn)為圓心,結(jié)合搜索半徑共同定義搜索鄰域,并根據(jù)搜索鄰域內(nèi)點(diǎn)云數(shù)對(duì)目標(biāo)點(diǎn)識(shí)別分類。當(dāng)激光點(diǎn)到傳感器的距離小于設(shè)定的最小的搜索半徑(SRmin),則以SRmin為搜索半徑進(jìn)行鄰域搜索。DROR的優(yōu)點(diǎn)是通過(guò)動(dòng)態(tài)變化的搜索半徑定義搜索鄰域,適應(yīng)了激光雷達(dá)掃描點(diǎn)隨著距離增加而減少的特性。DROR示意圖如圖2所示,聯(lián)合計(jì)算點(diǎn)到激光雷達(dá)的歐式距離和激光雷達(dá)水平分辨率獲得搜索半徑,從而確定以點(diǎn)為圓心,半徑為SR的圓形區(qū)域,散亂分布的點(diǎn)為激光雷達(dá)的掃描點(diǎn),l為激光束在k時(shí)間內(nèi)的變化位置,l1、l2為相鄰時(shí)刻的激光束,α為激光雷達(dá)的水平分辨率,η平面為t時(shí)刻內(nèi)激光束的掃描平面。DROR的局限性在于不能保持較好的去噪效果的同時(shí)保留較好的環(huán)境特征質(zhì)量。

圖2 DROR示意圖Fig.2 DROR diagrammatic sketch

3 基于反射強(qiáng)度與鄰域搜索組合的濾波算法

通過(guò)第二章對(duì)反射強(qiáng)度和鄰域搜索的介紹,可知單獨(dú)使用這兩種方法在對(duì)雪花噪聲去噪時(shí)都存在自身的局限性,對(duì)雪花噪聲的去噪效果較差。為降低雪花噪聲對(duì)激光雷達(dá)檢測(cè)產(chǎn)生的影響,可將兩種方法組合互補(bǔ),從而獲得更加精準(zhǔn)的分類識(shí)別點(diǎn)云,提高算法的去噪性能和環(huán)境特征質(zhì)量保留。

本文提出基于激光雷達(dá)反射強(qiáng)度與鄰域搜索組合的濾波算法,主要分為兩個(gè)步驟,第一步是對(duì)激光雷達(dá)點(diǎn)云反射強(qiáng)度統(tǒng)計(jì)分析,將帶噪聲的點(diǎn)云進(jìn)行粗識(shí)別分類。其中,對(duì)經(jīng)過(guò)統(tǒng)計(jì)得到的激光雷達(dá)點(diǎn)云反射強(qiáng)度,通過(guò)最小類內(nèi)方差法粗識(shí)別分類,得到非雪花點(diǎn)云和雪花點(diǎn)云,并將非雪花點(diǎn)云排除,不參與后續(xù)的點(diǎn)云處理,為下一步鄰域搜索的識(shí)別分類提供一個(gè)較好的初始點(diǎn)云。第二步是通過(guò)鄰域搜索的方法進(jìn)行精準(zhǔn)識(shí)別分類。通過(guò)搜索鄰域內(nèi)點(diǎn)云的數(shù)目,完成點(diǎn)云識(shí)別分類,將雪花精確識(shí)別分類并濾除?；诜瓷鋸?qiáng)度和鄰域搜索算法的組合算法,針對(duì)降雪場(chǎng)景中雪花噪聲分布的去噪問(wèn)題,具有比較好的效果。基于激光雷達(dá)反射強(qiáng)度和鄰域搜索的組合濾波算法去噪示意圖如圖3所示。

圖3 雪花噪聲去噪示意圖Fig.3 Snow noise denoising diagram

基于激光雷達(dá)反射強(qiáng)度的粗識(shí)別分類的動(dòng)態(tài)邊界閾值通過(guò)最小類內(nèi)方差法獲取。最小類內(nèi)方差法是圖像中一種常用的閾值分割方法,也適用于激光雷達(dá)數(shù)據(jù)反射強(qiáng)度邊界閾值分割。將激光雷達(dá)點(diǎn)云數(shù)據(jù)有序化,獲得數(shù)據(jù)點(diǎn)個(gè)數(shù)N,統(tǒng)計(jì)數(shù)據(jù)點(diǎn)云反射強(qiáng)度的分布概率,如表1所示。

表1 激光點(diǎn)反射強(qiáng)度的概率分布Tab.1 Probability distribution of laser pointreflection intensity

式(1)表示所有強(qiáng)度等級(jí)的概率和為1:

(1)

激光雷達(dá)數(shù)據(jù)的總個(gè)數(shù)為N,反射強(qiáng)度為i的點(diǎn)個(gè)數(shù)為ni,則反射強(qiáng)度為i的點(diǎn)概率pi為:

(2)

通過(guò)最小類內(nèi)方差法統(tǒng)計(jì)反射強(qiáng)度時(shí),若在等級(jí)強(qiáng)度為Tk時(shí)取得動(dòng)態(tài)邊界閾值,可根據(jù)動(dòng)態(tài)邊界閾值Tk劃分點(diǎn)云為兩個(gè)點(diǎn)云A、B,若A∈{1,2,…,Tk},B∈{Tk+1,Tk+2,…,Tm},點(diǎn)云A標(biāo)記為雪花點(diǎn)云,B標(biāo)記為非雪花點(diǎn)云,則目標(biāo)點(diǎn)屬于雪花點(diǎn)云A或?qū)儆诜茄┗c(diǎn)云B的概率分別為η1、η2:

(3)

(4)

式(5)、(6)表示雪花點(diǎn)云A的反射強(qiáng)度等級(jí)的均值μ1和方差σ1分別為:

(5)

(6)

式(7)為計(jì)算雪花點(diǎn)云A的反射強(qiáng)度等級(jí)為i時(shí)點(diǎn)的概率P(i|A):

(7)

式(8)、(9)表示雪花點(diǎn)云B的反射強(qiáng)度等級(jí)的均值μ2和方差σ2分別為:

(8)

(9)

式(10)計(jì)算雪花點(diǎn)云B的反射強(qiáng)度等級(jí)為i時(shí)點(diǎn)的概率P(i|B):

(10)

式(11)表示點(diǎn)云A、B類內(nèi)方差公式:

(11)

通過(guò)迭代的方式使得方差函數(shù)F(Tk)取得最小,即滿足min|F(Tk)|,1≤Tk≤Tm時(shí),可獲得基于反射強(qiáng)度的動(dòng)態(tài)邊界閾值Tk。動(dòng)態(tài)邊界閾值Tk對(duì)激光雷達(dá)原始點(diǎn)云進(jìn)行粗識(shí)別分類,當(dāng)動(dòng)態(tài)邊界閾值大于激光雷達(dá)點(diǎn)云數(shù)據(jù)強(qiáng)度pj時(shí),即Tk>pj,則將該點(diǎn)標(biāo)記為雪花點(diǎn),否則標(biāo)記為非雪花點(diǎn)。根據(jù)反射強(qiáng)度分割的點(diǎn)云是一個(gè)粗分類的過(guò)程,為了提高精確的點(diǎn)云分類,將經(jīng)過(guò)反射強(qiáng)度處理后的非雪花云作為鄰域搜索的輸入數(shù)據(jù),對(duì)該點(diǎn)進(jìn)行鄰域搜索。當(dāng)鄰域搜索的點(diǎn)云數(shù)的值大于設(shè)定的最小點(diǎn)云數(shù),即k>kmin時(shí),標(biāo)記該點(diǎn)為非雪花點(diǎn),否則標(biāo)記為雪花點(diǎn)。最后將標(biāo)記為雪花的點(diǎn)云濾除,獲得環(huán)境特征點(diǎn)云,實(shí)現(xiàn)雪花噪聲的去噪?；诜瓷鋸?qiáng)度和鄰域搜索組合算法的流程如圖4所示。

圖4 基于反射強(qiáng)度和鄰域搜索組合算法流程Fig.4 Algorithm flow based on reflection intensityand neighborhood search

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

本文算法通過(guò)RS-Bpearl 32線激光雷達(dá)采集的數(shù)據(jù)進(jìn)行驗(yàn)證,硬件設(shè)備為Inteli 7處理器,16 GB內(nèi)存電腦,系統(tǒng)環(huán)境為L(zhǎng)inux的機(jī)器人操作系統(tǒng)(ROS,Robot Operating System)。在實(shí)驗(yàn)驗(yàn)證中所用到的激光雷達(dá)數(shù)據(jù)源自北京聯(lián)合大學(xué)校區(qū)采集的實(shí)際激光雷達(dá)點(diǎn)云數(shù)據(jù),分為降雪和晴天兩種場(chǎng)景,總幀數(shù)均為3000幀,采集數(shù)據(jù)的路線均為同一環(huán)境,如圖5(e)所示。對(duì)以上場(chǎng)景的數(shù)據(jù)通過(guò)驗(yàn)證統(tǒng)計(jì)離群點(diǎn)濾波器、半徑離群點(diǎn)濾波器和動(dòng)態(tài)半徑離群點(diǎn)濾波器等不同的點(diǎn)云濾波器在降雪場(chǎng)景的雪花去噪效果,并與本文的方法相對(duì)比。

圖5 原始點(diǎn)云Fig.5 Original point cloud

如圖5中,(a)為實(shí)際道路場(chǎng)景和數(shù)據(jù)采集路線圖,且道路兩邊有建筑物、植物和車輛等環(huán)境特征。(b)為雪天場(chǎng)景原始點(diǎn)云俯視圖,(c)為前視圖；(d)為晴天場(chǎng)景原始點(diǎn)云俯視圖,(e)為前示圖；對(duì)比圖4的(b)與(d)、(c)與(e)可知標(biāo)記的區(qū)域H是激光雷達(dá)檢測(cè)到的雪花噪聲。在該段時(shí)間內(nèi),激光雷達(dá)檢測(cè)到單幀數(shù)據(jù)中降雪量的變化,如圖6所示。該變化趨勢(shì)可作為降雪量變化的參考量,并且可根據(jù)降雪量的變化驗(yàn)證算法在降雪場(chǎng)景的去噪效果和穩(wěn)定性。

圖6 降雪量變化趨勢(shì)Fig.6 Variation trend of snowfall

通過(guò)實(shí)驗(yàn)結(jié)果驗(yàn)證,本文提出的反射強(qiáng)度和鄰域搜索算法的組合算法是可行有效的。實(shí)驗(yàn)結(jié)果驗(yàn)證分為兩步分進(jìn)行,第一步是對(duì)環(huán)境特征質(zhì)量保留的驗(yàn)證,驗(yàn)證算法對(duì)環(huán)境特征質(zhì)量的保留效果。第二步對(duì)雪花噪聲分析,驗(yàn)證本文算法對(duì)雪花噪聲的去噪效果。通過(guò)實(shí)驗(yàn)和數(shù)據(jù)分析,定性與定量的分析本文算法與常規(guī)濾波器的去噪性能。

4.1 實(shí)驗(yàn)結(jié)果

統(tǒng)計(jì)離群點(diǎn)濾波器:統(tǒng)計(jì)離群點(diǎn)濾波器首先通過(guò)(SOR,Statistical Outlier Removal Filter)迭代點(diǎn)云中的每個(gè)點(diǎn),并計(jì)算到其k個(gè)最近鄰點(diǎn)的距離的平均值和標(biāo)準(zhǔn)差。該濾波器遍歷所有點(diǎn),并刪除到其k個(gè)鄰居的平均距離大于閾值的點(diǎn)。統(tǒng)計(jì)離群點(diǎn)濾波器應(yīng)用于降雪場(chǎng)景時(shí),可以達(dá)到將離散的雪花噪聲濾除,但是對(duì)于密度大的雪花噪聲的濾除效果較差,該濾波器的濾除效果如圖7所示。

圖7 SOR濾波器Fig.7 SOR filter

半徑離群點(diǎn)濾波器:半徑離群點(diǎn)濾波器(ROR,Radius Outlier Removal Filter)遍歷點(diǎn)云中的每個(gè)點(diǎn),并計(jì)算在設(shè)定搜索半徑R內(nèi)找到的點(diǎn)云數(shù)[16]。若搜索到的點(diǎn)云數(shù)小于設(shè)定的最小點(diǎn)云數(shù),則該點(diǎn)被移除。該濾波器的搜索半徑是固定不變的,當(dāng)搜索半徑足夠小時(shí)可以達(dá)到較好的濾除效果,但是,在去噪的同時(shí)會(huì)嚴(yán)重?fù)p傷環(huán)境特征的質(zhì)量。當(dāng)搜索半徑R=0.05時(shí),可以濾除較多的雪花噪聲,但密度大的一簇雪花噪聲沒(méi)有濾除,去噪的效果如圖8所示。

圖8 ROR濾波器Fig.8 ROR filter

動(dòng)態(tài)半徑離群點(diǎn)濾波器:該濾波器考慮了激光雷達(dá)掃描點(diǎn)距離越遠(yuǎn)越稀疏的特性,為每個(gè)掃描點(diǎn)定義了一個(gè)動(dòng)態(tài)變化的搜索半徑,進(jìn)一步解決了遠(yuǎn)距離損害環(huán)境特征質(zhì)量的難點(diǎn)[15]。去噪性能與半徑濾波器相比,該濾波器在總體效果上有了一定的提升,但是密度大的一簇雪花噪聲仍然存在,沒(méi)有濾除,其去噪效果如圖9所示。

圖9 DROR濾波器Fig.9 DROR filter

本文方法:本文方法是基于反射強(qiáng)度與鄰域搜索組合的算法,首先通過(guò)反射強(qiáng)度的粗濾波后為鄰域搜索提供一個(gè)較好的初始點(diǎn)云,再通過(guò)鄰域搜索方法將雪花噪聲識(shí)別濾除,去噪的效果如圖10所示。

圖10 本文方法濾除效果Fig.10 the filtering effect of this method

濾波器參數(shù)如表2所示。

表2 濾波器參數(shù)Tab.2 Filter parameters

4.2 實(shí)驗(yàn)分析

實(shí)驗(yàn)驗(yàn)證的數(shù)據(jù)采集于校園道路,故只對(duì)距離激光雷達(dá)20 m范圍內(nèi)的數(shù)據(jù)展開(kāi)分析,檢測(cè)范圍可滿足無(wú)人車的安全行駛。針對(duì)單幀激光雷達(dá)點(diǎn)云數(shù)據(jù)的去噪性能和保留環(huán)境特征性能進(jìn)行分析,圖11展示了本文算法與常規(guī)濾波器對(duì)環(huán)境特征質(zhì)量保留的結(jié)果。從整體來(lái)看,本文方法在20 m范圍內(nèi)的環(huán)境特征保留效果僅次于SOR濾波器。雖然SOR濾波器保留環(huán)境特征的效果最好,但是雪花噪聲去噪性能最差,只將局部離群的雪花噪聲濾除,不能濾除密度大的雪花噪聲,如圖12所示。

圖11 環(huán)境特征保留Fig.11 Preservation of environmental features

本文算法對(duì)雪花噪聲去噪的效果如圖12所示。從3000幀數(shù)據(jù)中均衡抽樣100幀數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,本文算法在該時(shí)間段內(nèi)降噪效果處于96 %左右,DROR算法的去噪效果為88.5 %,提高了約8 %。本文方法濾除率的波動(dòng)幅度不超過(guò)1 %,而DROR在降雪量發(fā)生變化時(shí)濾除效果波動(dòng)比較大,波動(dòng)范圍為3 %。ROR濾波器的去噪效果比DROR的效果差,而SOR濾波器降噪效果在20 %左右,去噪效果最差。常規(guī)的濾波器對(duì)雪花噪聲去噪時(shí)具有局限性,不能達(dá)到去噪的同時(shí)保留環(huán)境特征質(zhì)量的目的。

圖12 平均濾除率Fig.12 Average filtration rate

本文方法結(jié)合了激光雷達(dá)掃描點(diǎn)的強(qiáng)度信息,在雪花噪聲去噪和環(huán)境特征保留上都具有很大的改善。與DROR濾波器相比,本文方法在進(jìn)行粗識(shí)別分類時(shí)分割出局部的環(huán)境特征,所以在環(huán)境特征保留上具有更好的效果,同時(shí)通過(guò)改變給定點(diǎn)云數(shù)提高對(duì)雪花噪聲的濾除。

5 結(jié) 語(yǔ)

針對(duì)降雪場(chǎng)景中雪花噪聲造成激光雷達(dá)的誤檢的問(wèn)題,提出了一種基于激光雷達(dá)反射強(qiáng)度與鄰域搜索組合濾波算法。該方法充分利用了點(diǎn)云反射強(qiáng)度和噪聲點(diǎn)云空間分布的特點(diǎn),實(shí)現(xiàn)了降雪場(chǎng)景中雪花噪聲的濾除,有效的解決了雪花噪聲引起的影響的問(wèn)題。在降雪場(chǎng)景采集的真實(shí)數(shù)據(jù)驗(yàn)證下實(shí)驗(yàn)結(jié)果表明:與常規(guī)的去噪算法相比,其雪花噪聲去噪率和環(huán)境特征質(zhì)量皆得到了提升,雪花噪聲去噪率由88.5 %提升至96 %,且具有較好的穩(wěn)定性；在達(dá)到去噪目的的前提下,環(huán)境特征質(zhì)量的保留效果皆優(yōu)于其余方法。本文方法能夠降低降雪場(chǎng)景中激光雷達(dá)的誤檢率,滿足無(wú)人車所需的高精度、高可靠性的環(huán)境感知信息。下一步的工作重點(diǎn)則是繼續(xù)對(duì)算法進(jìn)行優(yōu)化,提高實(shí)時(shí)性以及去噪效果,同時(shí)對(duì)降雨等天氣進(jìn)行去噪效果評(píng)估。