杜 芳 任明武
(南京理工大學(xué)計(jì)算機(jī)科學(xué)與工程學(xué)院 南京 210094)
近年來(lái),隨著科技發(fā)展,人們對(duì)智能交通的關(guān)注度持續(xù)上升,關(guān)于移動(dòng)機(jī)器人、自主導(dǎo)航車以及高級(jí)輔助駕駛系統(tǒng)(Advanced Driver Assistant Systems,ADAS)的研究與應(yīng)用日益增多,道路環(huán)境的感知成為其研究的重要組成部分[1]。激光雷達(dá)以其較高的測(cè)距精度及角度分辨率,在道路環(huán)境感知中發(fā)揮了重要的作用。目前激光雷達(dá)在智能交通環(huán)境理解領(lǐng)域的研究主要集中在障礙物檢測(cè)、道邊檢測(cè)以及可通行區(qū)域信息提取等方面[2]。
根據(jù)文獻(xiàn)[3~4]可知,采用激光雷達(dá)檢測(cè)障礙物主要根據(jù)障礙物的高度信息進(jìn)行判定。本文結(jié)合多層激光雷達(dá)的特點(diǎn),使用二維信息,根據(jù)雷達(dá)數(shù)據(jù)點(diǎn)在物體表面呈密集分布的特征,提出了一種基于數(shù)據(jù)凹凸性分布的改進(jìn)多密度DBSCAN算法的快速障礙檢測(cè)方法(Convexity and Multi Density-Based Spatial Clustering of Applications with Noise,CM-DBSCAN)以及一種基于幀間慣導(dǎo)信息匹配的動(dòng)態(tài)障礙物檢測(cè)方法,本文方法能夠很好地提取道路環(huán)境中的障礙物信息,實(shí)現(xiàn)障礙物檢測(cè)。
與傳統(tǒng)攝像頭相比,在環(huán)境感知應(yīng)用上,激光雷達(dá)受環(huán)境光影響很小且穩(wěn)定性高[5]。本文采用德國(guó)生產(chǎn)的IBEO-LUX-2010型號(hào)四線激光雷達(dá),它是一種三維多層激光雷達(dá),與單線激光雷達(dá)相比,其數(shù)據(jù)量大、檢測(cè)精度高、探測(cè)范圍廣,可以很好地滿足實(shí)際應(yīng)用中的實(shí)時(shí)性以及有效性的要求[6]。
從實(shí)際應(yīng)用角度出發(fā),為保證激光雷達(dá)對(duì)前方道路觀察的全面和準(zhǔn)確,將四線激光雷達(dá)安裝于汽車正前方的中間位置。如圖1所示,車體坐標(biāo)系以車輛的慣導(dǎo)中心為坐標(biāo)原點(diǎn),用Ov表示,激光雷達(dá)坐標(biāo)系以激光雷達(dá)的中心為坐標(biāo)原點(diǎn),用Ol表示。雷達(dá)掃描方式?jīng)Q定了掃描距離ρ和掃描角度唯一確定雷達(dá)點(diǎn)P(x ,y,z) 的位置[7],故本文使用一種基于掃描角度劃分的分層式雷達(dá)坐標(biāo)轉(zhuǎn)換方法。其中,Vres是其垂直掃描精度為固定值0.8°,θv表示各層雷達(dá)數(shù)據(jù)點(diǎn)垂直方向的掃描角度,θh表示水平方向的掃描角度,AngleTics為固定值11520,c為常量,根據(jù)雷達(dá)屬性確定c值。關(guān)系式為
式(3)中 i=0,1,2,3分別表示對(duì)應(yīng)的四層數(shù)據(jù)。由于安裝固定后的雷達(dá)與車體的姿態(tài)一致,因此利用車輛的位姿傳感器測(cè)量車輛姿態(tài),由車體運(yùn)動(dòng)的航向角,俯仰角和翻滾角得到旋轉(zhuǎn)變換矩陣R[8]。 x′、y′、z′表示從激光雷達(dá)坐標(biāo)系到車體坐標(biāo)系的偏移。結(jié)合式(1)、式(2),將雷達(dá)采集到的數(shù)據(jù)轉(zhuǎn)換到車輛坐標(biāo)系下,并進(jìn)行姿態(tài)校正,可以克服道路結(jié)構(gòu)及車體運(yùn)動(dòng)導(dǎo)致的數(shù)據(jù)漂移情況。由于本文中的實(shí)驗(yàn)只需二維信息,故將得到的三維點(diǎn)云數(shù)據(jù)投影至二維平面內(nèi)。
圖1 坐標(biāo)系關(guān)系示意圖
本實(shí)驗(yàn)只考慮車前30m內(nèi)的障礙物,故除去距離大于30m的雷達(dá)點(diǎn)。對(duì)于同一層雷達(dá)線,雷達(dá)數(shù)據(jù)掃描時(shí)的障礙返回點(diǎn)與其他點(diǎn)相比,整個(gè)數(shù)據(jù)分布呈凹凸性。根據(jù)DBSCAN原始聚類方法特點(diǎn)可知,縮小核心點(diǎn)的選取范圍是降低時(shí)間復(fù)雜度、改進(jìn)搜索速度的有效途徑之一[9]。因此本實(shí)驗(yàn)根據(jù)數(shù)據(jù)分布的凹凸特性,計(jì)算同一掃描角度θh出發(fā)的第i層雷達(dá)點(diǎn)相對(duì)于其它三層(用 j表示)雷達(dá)點(diǎn)的距離表示:
本方法只需知道同一角度掃描的點(diǎn)距雷達(dá)的最近距離distmin,對(duì)于同一角度出發(fā)的點(diǎn),只要distmin小于特定閾值則為障礙候選點(diǎn)。得出中所有接近distmin的對(duì)應(yīng)點(diǎn),認(rèn)為其為初步濾取的障礙點(diǎn)。
DBSCAN是經(jīng)典的基于密度的聚類算法,密度相連點(diǎn)的最大集合即為簇,該方法將具有足夠高密度的區(qū)域劃分為簇,可發(fā)現(xiàn)任意形狀的聚類[10~11]。DBSCAN需要兩個(gè)參數(shù):掃描半徑Eps和最小包含點(diǎn)數(shù)MinPts,這兩個(gè)參數(shù)值在大多數(shù)情況下是人為給定的,其中參數(shù)Eps的值在很大程度上影響最終的聚類結(jié)果[12]。因此確定參數(shù)Eps的值是DBSCAN算法獲得較好聚類結(jié)果的關(guān)鍵。
因此,本文基于雷達(dá)數(shù)據(jù)特性,在3.1節(jié)數(shù)據(jù)預(yù)處理的基礎(chǔ)上針對(duì)參數(shù)敏感且無(wú)法適用于多密度數(shù)據(jù)集聚類的缺點(diǎn)[13],提出了一種快速的基于多密度雷達(dá)數(shù)據(jù)的動(dòng)態(tài)參數(shù)DBSCAN聚類方法CM-DBSCAN。
由于該雷達(dá)為分層式對(duì)稱掃描,故Eps參數(shù)確定如下:
其中Ce為聚類半徑系數(shù),θhres為水平方向的掃描角分辨率。完整的基于凹凸性分布的多密度聚類算法流程如下:
1)同一角度掃描出發(fā)的點(diǎn) pi,i表示為第i層的點(diǎn),判斷點(diǎn) pi與其它層的同一角度出發(fā)的點(diǎn)的最小距離,若最小距離小于閾值ε,則認(rèn)為其是候選點(diǎn)point;
2)設(shè)定MinPts:=4,根據(jù)掃描角度和距離確定Eps的值;
3)判斷輸入點(diǎn)point是否為核心對(duì)象;
4)找出核心對(duì)象的E鄰域(即Eps范圍)中的所有直接密度可達(dá)點(diǎn);
5)Until所有輸入點(diǎn)都判斷完畢;
6)Repeat;
7)針對(duì)所有核心對(duì)象的E鄰域內(nèi)所有直接密度可達(dá)點(diǎn)找到最大密度相連對(duì)象集合,中間涉及到一些密度可達(dá)對(duì)象的合并;
8)Until所有核心對(duì)象的E領(lǐng)域都遍歷完畢。
在無(wú)人駕駛車輛行駛過(guò)程中,動(dòng)態(tài)障礙物如車輛、行人帶來(lái)的安全隱患非常大[14]。根據(jù)文獻(xiàn)[15]可知,目前研究的障礙物檢測(cè)方法大多都是基于單幀數(shù)據(jù)來(lái)進(jìn)行處理的,無(wú)法判斷障礙物的狀態(tài)。文獻(xiàn)[16~17]中,Aftatah M、陸興華等解釋了慣性導(dǎo)航系統(tǒng)(簡(jiǎn)稱慣導(dǎo))是一種不依賴于外部信息、以牛頓力學(xué)定律為基礎(chǔ)的自主式導(dǎo)航系統(tǒng),其通過(guò)測(cè)量載體在慣性參考系的加速度,將它對(duì)時(shí)間進(jìn)行積分,且把它變換到導(dǎo)航坐標(biāo)系中,就能夠得到在導(dǎo)航坐標(biāo)系中的速度、偏航角和位置等信息。因此本文根據(jù)車體運(yùn)動(dòng)信息,提出了一種基于慣導(dǎo)信息的動(dòng)態(tài)障礙物檢測(cè)方法。該方法在實(shí)車試驗(yàn)中得到了驗(yàn)證。根據(jù)第3章實(shí)驗(yàn)方法可以得出單幀數(shù)據(jù)中的障礙物檢測(cè)結(jié)果,對(duì)相鄰幀的檢測(cè)結(jié)果進(jìn)行障礙物匹配,得出匹配結(jié)果。具體的匹配方法如下:
1)計(jì)算候選障礙物質(zhì)心的位置:通過(guò)雷達(dá)數(shù)據(jù)點(diǎn)的距離平均值和角度平均值計(jì)算。
2)對(duì)檢測(cè)出的障礙物結(jié)果進(jìn)行標(biāo)號(hào),根據(jù)車體運(yùn)動(dòng)信息即車輛的位移矢量可知,若滿足相鄰幀的障礙物位移矢量與車體位移矢量之差η0<Δd<η1,且相鄰幀的障礙物點(diǎn)數(shù)之差 Δn<η2、面積之差Δs<η3則為動(dòng)態(tài)障礙物。
本文中的實(shí)驗(yàn)平臺(tái)是南京理工大學(xué)無(wú)人駕駛車平臺(tái)。根據(jù)獲得的實(shí)驗(yàn)數(shù)據(jù),利用上文中所提出的激光雷達(dá)數(shù)據(jù)預(yù)處理方法、障礙檢測(cè)方法以及動(dòng)態(tài)障礙檢測(cè)方法進(jìn)行了實(shí)驗(yàn),且不考慮道路邊沿為障礙。實(shí)驗(yàn)結(jié)果表明在結(jié)構(gòu)化道路環(huán)境中,單幀數(shù)據(jù)中的障礙物以及幀間的動(dòng)態(tài)障礙物的檢測(cè)結(jié)果能夠很好地滿足準(zhǔn)確性和實(shí)時(shí)性要求。
本文的實(shí)驗(yàn)場(chǎng)景如圖2所示。表1為算法改進(jìn)前后的實(shí)驗(yàn)場(chǎng)景中對(duì)應(yīng)的連續(xù)20幀數(shù)據(jù)的搜索時(shí)間統(tǒng)計(jì)表。單幀雷達(dá)數(shù)據(jù)的障礙檢測(cè)結(jié)果如圖3、圖4所示,其中圖3為DBSCAN算法改進(jìn)前的檢測(cè)結(jié)果,圖4表示改進(jìn)算法后的檢測(cè)結(jié)果,白色矩形表示單幀數(shù)據(jù)障礙物。圖5為基于幀間數(shù)據(jù)的靜態(tài)、動(dòng)態(tài)障礙物檢測(cè)實(shí)驗(yàn)結(jié)果,淺色矩形表示靜態(tài)障礙物,深色矩形表示動(dòng)態(tài)障礙物。
圖2 實(shí)驗(yàn)場(chǎng)景
圖3 改進(jìn)前原始DBSCAN算法檢測(cè)結(jié)果
圖4 改進(jìn)后CM-DBSCAN算法檢測(cè)結(jié)果
表1 算法改進(jìn)前后搜索時(shí)間(單位:ms)比較
圖5 幀間動(dòng)態(tài)障礙物檢測(cè)結(jié)果
由5.1實(shí)驗(yàn)結(jié)果可知,基于凹凸性分布的多密度障礙檢測(cè)方法以及基于慣導(dǎo)信息匹配的動(dòng)態(tài)障礙物檢測(cè)方法得到了很好的驗(yàn)證。實(shí)驗(yàn)表明,利用數(shù)據(jù)的凹凸性分布進(jìn)行數(shù)據(jù)預(yù)處理可以極大地提高后續(xù)聚類算法的搜索速度。由表1可知,改進(jìn)后的CM-DBSCAN算法的搜索時(shí)間相對(duì)于改進(jìn)前縮短了45.23%。障礙檢測(cè)準(zhǔn)確率方面,根據(jù)圖3、圖4和圖5可以看出基于CM-DBSCAN算法的障礙物檢測(cè)數(shù)目更加符合實(shí)際場(chǎng)景,且檢出率更高?;趹T導(dǎo)信息匹配檢測(cè)算法可成功檢測(cè)出場(chǎng)景中的動(dòng)態(tài)障礙物,多次實(shí)驗(yàn)得出單幀數(shù)據(jù)的障礙物檢測(cè)成功率可達(dá)92.32%,動(dòng)態(tài)障礙物檢測(cè)成功率為89.79%,由此可知,CM-DBSCAN算法具有較高的檢測(cè)準(zhǔn)確率。
本文基于四線激光雷達(dá)來(lái)檢測(cè)道路環(huán)境中的障礙物,提出了一種基于數(shù)據(jù)凹凸性分布多密度DBSCAN的快速障礙檢測(cè)方法(CM-DBSCAN)。該算法通過(guò)雷達(dá)數(shù)據(jù)的凹凸性分布得出障礙候選點(diǎn),縮小搜索范圍,從而加快了搜索速度。本文根據(jù)雷達(dá)點(diǎn)的密度分布,提出了一種多密度DBSCAN聚類方法,該算法能夠根據(jù)數(shù)據(jù)特征自適應(yīng)地調(diào)整參數(shù)Eps,解決了原始DBSCAN聚類方法的單一參數(shù)問(wèn)題以及只能對(duì)單一密度區(qū)域進(jìn)行聚類的問(wèn)題。另外,本文還提出了一種基于慣導(dǎo)信息的動(dòng)態(tài)障礙物檢測(cè)方法。該方法根據(jù)相鄰幀的障礙物檢測(cè)結(jié)果,結(jié)合車體運(yùn)動(dòng)信息進(jìn)行障礙物匹配,得出障礙檢測(cè)結(jié)果。根據(jù)得到的障礙檢測(cè)和道邊檢測(cè)的結(jié)果建立局部地圖,可作為未來(lái)的研究方向之一。