五、畸變矯正—讓世界不在扭曲

? ? ? ? ? 這篇博文所要講述的內(nèi)容，是標(biāo)定的主要用途之一：矯正攝像機的畸變。對于圖像畸變矯正的方法，張正友教授也在其大作“A Flexible New Technique forCamera Calibration”中給出。

? ? ? ? ?玉米在這里先為大家介紹一下，攝像機畸變主要有哪幾種以及這些討厭的畸變從何而來。

攝像機畸變主要包含：鏡像畸變、切向畸變。玉米在這里分別講述一下，并引用一些圖片，使大家對他們有一個直觀的認(rèn)識。

? ? ? ? ?1．徑向畸變：產(chǎn)生原因是光線在遠(yuǎn)離透鏡中心的地方比靠近中心的地方更加彎曲徑向畸變主要包含桶形畸變和枕形畸變兩種。下面兩幅圖是這兩種畸變的示意：

? ? ? ? ? 它們在真實照片中是這樣的：

? ? ? ? ? ?2.切向畸變：產(chǎn)生的原因透鏡不完全平行于圖像平面，這種現(xiàn)象發(fā)生于成像儀被粘貼在攝像機的時候。下面圖片來自于《學(xué)習(xí)opencv》p413。

? ? ? ? ? 畸變矯正，首先應(yīng)該知道畸變系數(shù)，然后做與畸變相反的變換，消除畸變。

? ? ? ? ? 張氏標(biāo)定法中只關(guān)注徑向畸變。我們是按照張氏標(biāo)定，計算畸變系數(shù)的。那么，下面讓我們來看一下徑向畸變的數(shù)學(xué)表達(dá)。

? ? ? ? ? 因為在實際情況下，徑向畸變較小，所以其可以用主點(principle point)周圍的泰勒級數(shù)展開的前幾項進行描述。張氏標(biāo)定法，利用前兩項來確定徑向畸變的畸變系數(shù)。數(shù)學(xué)表達(dá)式如下：

其中，代表理想無畸變的像素坐標(biāo)，代表實際徑像畸變的情況下的像素坐標(biāo)，代表主點，代表理想無畸變時的連續(xù)圖像坐標(biāo)，代表實際徑像畸變的情況下的連續(xù)圖像坐標(biāo)。K₁,k₂代表前兩階的畸變參數(shù)。。

???????? 那么對于圖像上的任意一點，我們有兩個個等式。化成矩陣形式：

? ? ? ? 通過前面兩篇博文所介紹的內(nèi)容，求得的攝像機模型，我們可以計算出（可通過攝像機模型直接解出）、（與是否畸變無關(guān)，直接在已求得的內(nèi)參陣中得到）、（即直接讀取的有畸變的像素坐標(biāo)），（由攝像機模型可以通過物體的世界坐標(biāo)點解出）。那么，兩個方程兩個未知數(shù)，我們用一點就可以求的徑向畸變了。因為我們有n張圖片，每張圖片上有m個點，所以我們可以得到2mn個等式。運用最小二乘法對結(jié)果進行優(yōu)化。可用下式解徑向畸變k=[k₁,k₂]。

其中，D等式左邊的方程的系數(shù)矩陣，d是等式右邊的有畸變的像素坐標(biāo)與無像素坐標(biāo)之差構(gòu)成的矩陣。有上式，一步到位，計算出畸變系數(shù)。

?????? 但是，但是，大家先別急。我們的參數(shù)求解之路還沒有結(jié)束。為了使標(biāo)定得到的估算實際的、存在徑向畸變的攝像機參數(shù)。我們需要把，剛剛求解得到的畸變參數(shù)，連同前面得到的理想無畸變條件下的內(nèi)外參數(shù)一起，進行極大似然估計。

?????? 這里的極大似然估計方法與上一篇博文《極大似然參數(shù)估計》相同，玉米在這里就不在贅述。

依然是以最小化下列函數(shù)為目標(biāo)，只不過此次參數(shù)估計中又多加了k1,k2：

? ? ? ? 然后，仍然用Levenberg-Marquardt算法進行計算。最終就得到了，我們想要在考慮畸變情況下，攝像機的參數(shù)。

張氏標(biāo)定，大功告成！

? ? ? ?下面讓我們矯正畸變：

? ? ? ? 就是這樣一個簡單的公式，就完成了畸變矯正。代表矯正畸變后的像素坐標(biāo)，代表實際徑像畸變的情況下的圖像的像素坐標(biāo)。

? ? ? ? 最后給大家展示兩幅圖，讓大家看一下畸變矯正的成果：

? ? ? ? ? 這就是張氏標(biāo)定帶給我們的第一個直觀的好處，后面還有兩個大好處，玉米將為您介紹。

? ? ? ? ? 玉米才疏學(xué)淺，懇請請大家批評指正！

總結(jié)

以上是生活随笔為你收集整理的五、畸变矫正—让世界不在扭曲的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。