梯度下降法在雙目視覺導(dǎo)航中的一個運用

廣東工業(yè)大學(xué) 袁文波 陳慶瑜 高洪兵 王欽若

1 引言

本文將討論的是一種基于雙目視覺且路線固定的視覺導(dǎo)航方式，如無人搬運機器人，無人快遞小車的導(dǎo)航。其導(dǎo)航特點是路線已知且相對固定，機器只需沿著設(shè)定好的路線前進就能達到目的地。

2 無人快遞車導(dǎo)航原理概述

要實現(xiàn)自動導(dǎo)航，那么這就要求機器必須知道自己的實時位姿與正確位姿的相對偏差。由于快遞小車的運行路線固定，我們可以事先人為控制讓安裝有雙目攝像頭的小車沿著正確的道路行走一遍，雙目攝像頭相對小車位置固定，并對前進路線環(huán)境進行拍攝，我們每一幀將有兩張照片，左攝像頭及右攝像頭所拍攝的照片，我們對兩張照片進行特征點匹配，本文選用的是SIFT特征點匹配方式，并將匹配上的特征點的特征向量及其在小車坐標(biāo)系（一般選小車的中心點或者小車的某個頂點作為小車坐標(biāo)系的原點，小車的正前基線、水平基線、豎直方向基線分別與小車坐標(biāo)系的三個坐標(biāo)軸平行）的坐標(biāo)按照片先后順序進行存儲。雙目視覺確定匹配上的特征點在小車坐標(biāo)系的坐標(biāo)方法已經(jīng)很成熟在這里將不再贅述。當(dāng)小車在正常運作時，不斷拍攝前方環(huán)境，并在左右?guī)掌羞M行特征點匹配，將匹配上的特征點在與之前保存的特征點匹配，選取匹配度最高的一幀作為此次調(diào)整的正確位姿。并根據(jù)匹配上的特征點的坐標(biāo)計算出小車實時位姿與正確位姿的偏差，從而達到導(dǎo)航的目的。本文將要詳細闡述的是在小車正常運行時如何檢測小車實時位姿與正確位姿的偏差。

圖1 SIFT特征點匹配效果

3 實時位姿與正確位姿偏差計算

圖2 實時位姿與正確位姿偏差示意圖

當(dāng)小車正常工作時，進行左右兩張照片特征點的匹配，及計算出匹配點相對小車坐標(biāo)系的坐標(biāo)。將匹配上的特征點與系統(tǒng)里事先掃描路線預(yù)存的特征點以幀為跨度進行匹配，選取匹配點數(shù)最高的一幀，接下來就是要根據(jù)這些匹配點在其各自坐標(biāo)系的坐標(biāo)求出這兩個坐標(biāo)系的變換矩陣，有了變換矩陣也就知道了小車實際位姿與正確位姿的偏差。有了這個偏差便能對小車位姿進行調(diào)控。其實這就轉(zhuǎn)化成了一個數(shù)學(xué)問題，已知同一組點在兩個三維坐標(biāo)系的坐標(biāo)求兩三維坐標(biāo)系之間的變換矩陣。首先我們來觀察一下，兩個空間三維直角坐標(biāo)系能夠通過怎樣的變換使得兩個坐標(biāo)系重合，很明顯我們可以先進行平移使得兩個坐標(biāo)系原點重合，重合后分別先后繞三個坐標(biāo)軸進行三次旋轉(zhuǎn)，一定能使兩坐標(biāo)系重合。以下四個變換矩陣分別對應(yīng)以上的四個步驟，平移，繞X軸旋轉(zhuǎn)，繞Y軸旋轉(zhuǎn)，繞Z軸旋轉(zhuǎn)。

最終變換矩陣R計算如下：

隨機選取五個點，在坐標(biāo)坐標(biāo)系1的坐標(biāo)構(gòu)成矩陣Y1則Y1矩陣如下,每列代表一個點的坐標(biāo)x,y,z，添1為了計算方便。

則這五點在第二個坐標(biāo)系的坐標(biāo)計算公式如下：

實際情況是我們有了五個點在兩個坐標(biāo)系的坐標(biāo)，要求坐標(biāo)系的變換矩陣,即求那六個參數(shù)，如果這些坐標(biāo)都是精確的，那就是一個解方程的問題，但實際運用時測量都是有誤差的，本文介紹一種在利用梯度下降的方法在坐標(biāo)有誤差的情況下求出精度符合要求的變換矩陣。我們假設(shè)xt=50，yt=-80zt=200，xr=3.1416/2，yr=3.1416/3，zr=3.1416/5帶入各變換矩陣，利用公式(7)求得該五點在坐標(biāo)系二的坐標(biāo)如下：

再以xt,yt,zt,xr,yr,zr為未知數(shù)通過公式(9)計算，

則在坐標(biāo)系二的五點的坐標(biāo)為包含xt,yt,zt,xr,yr,zr的表達式，再通過公式（10）計算Loss：

此時Loss為包含xt，yt，zt，xr，yr，zr的表達式，將Loss分別對xt,yt,zt,xr,yr,zr求偏導(dǎo)，利用梯度下降法，不斷迭代，使Loss下降到一個可以接受的范圍，則便可得到各個參數(shù)較為準(zhǔn)確的值。之后再在Y2中加入隨機的正態(tài)分布誤差來模仿實際測量帶來的誤差，用同樣的方法得到各個參數(shù)較為準(zhǔn)確的值。如下兩矩陣分別是以u=0，σ=1，u=0，σ=2，隨機生成的干擾數(shù)據(jù)。

noise將以如下形式加載到數(shù)據(jù)中：

表1所示反應(yīng)了在不加噪聲及加不同噪聲時梯度下降法的計算精度。圖2為Loss，及其中一個參數(shù)xr在無噪聲情況下，隨著迭代次數(shù)增加而收斂的圖像，其余參數(shù)及在有噪聲情況下的收斂趨勢與下圖類似。

表1 梯度下降法求的結(jié)果

圖2 Loss及xr收斂圖