一種簡單而精確的魚眼圖像校正算法研究

舒旭

摘 要： 針對魚眼圖像的校正提出了一種有效區域提取算法，并在魚眼圖像的球面物投影平面展開時與目標半立方體的校正平面建立線性映射。實驗結果表明，該算法能有效地提取魚眼圖像輪廓的有效區域，同時半立方體校正平面極大地減少了校正圖像需插值的像素點數量，提高了校正速度。

關鍵詞： 魚眼圖像； 等距投影； 畸變校正； 輪廓提取

中圖分類號：TP751 文獻標志碼：A 文章編號：1006-8228（2017）06-51-04

A simple and accurate fisheye image correction algorithm

Shu Xu

（XiangNan Preschool Education College， Chenzhou， Hunan 423000， China）

Abstract： This paper presents an effective region extraction algorithm for fish eye image correction， and establishes a linear mapping between semi-spherical surface projection plane and semi-cube projection plane in the process of fisheye image correction. The experimental results show that the algorithm can effectively extract the effective contour region of the fish eye image， and at the same time， the semi-cube correction plane can greatly reduce the number of pixels that need to be interpolated， which improve the correction speed.

Key words： fisheye image； equidistant projection； distortion correction； contour extraction

0 引言

人們對“看的更廣”的要求日益提高，低成本的、具有廣角或魚眼鏡頭的大視場攝像系統備受市場青睞。然而采用大視場攝像系統采集信息時，因鏡頭成像面大小固定，采集過程中為獲得更多的拍攝信息使得采集的圖像變形較嚴重，所以必須對畸變的圖像進行校正以適合人眼觀看。

目前，具有代表性的魚眼圖像校正算法有以下幾類：經緯映射法、多項式坐標變換法、球面透視投影約束法。楊玲和成運使用經緯映射來校正魚眼圖像[1]。該算法認為魚眼圖像中扭曲的場景可以用地球儀經度線來表示，利用同一經度上的像素點在校正后的圖像中應具有相同縱坐標的特點，把經緯度網格點上的像素值映射到校正平面上與之對應的位置。這種直接根據比例建立相應投影關系來實現圖像校正的算法雖較容易實現，但校正后的圖像存在一定的失真，校正圖像質量有待提高。廖士中和高培煥等提出的多項式坐標變換算法[2]是根據魚眼圖像的像素點位置通過最小二乘發擬合誤差平方最小以計算校正后圖像對應點，該方法對于魚眼鏡頭產生的徑向畸變非線性失真，實驗證明往往需要使用5階二元多項式，需要N個控制點，計算量非常巨大，很難應用于實際操作應用中。英向華和胡占義提出的球面透視投影約束校正算法[3]思想是在球面成像投影模型下，利用空間直線的球面透視投影為大圓這一約束來恢復魚眼變形校正，該算法在選取理想的采樣點后效果較好，但建立目標函數時，為權衡精度和速度同樣需要計算5階二元多項式，且采樣點的合理選取也是一個很難把握的地方，因該算法不能自適應自動選取采樣點，所有采樣點需手動選取，選取的結果直接關系到目標函數的擬合計算次數和圖像校正的效果，不具有普遍性。

通過以上分析，本文旨在尋求一種不需要對魚眼鏡頭進行參數標定，算法結構簡單易于實現、具有較大的普適性，且校正效果較好的畸變圖像校正算法。提出以等距投影相關理論為基礎，在精確提取魚眼圖像輪廓的前提下，運用球面等距投影校正算法能夠兼顧算法速度和圖像校正效果。

1 魚眼圖像輪廓的提取及圓心的確定

目前，提取魚眼圖像輪廓的常用算法有兩種。

第一種是面積統計法，該算法的原理是事先確定一個閾值T，求出圖像中像素值大于閾值的像素點個數和sum，并同時統計出像素值大于或等于閾值T的橫坐標和∑X、縱坐標和∑Y。計算圓心（x0，y0）和半徑R的公式如下：

⑴

⑵

⑶

該算法的缺點是只適應于魚眼圖像圓形區域內像素值較小的像素點個數少的情況，當圓形區域內有較多黑像素點時，該算法提取出來的圖像輪廓誤差較大。

第二種是逐行逐列掃描法[4-5]。該算法的原理是首先確定一個閾值T，分別從四周開始向中心逐行逐列掃描。如果逐行掃描時存在像素點的亮度值大于閾值T，則該像素點的橫坐標為有效圖像區域的邊界，同理確定出縱坐標邊界。該算法的缺點是只適應于圖像圓形有效區域的亮度值遠大于圓形區域之外的情況，如圓形區域外存在一些亮度值較大的噪聲干擾點，則該算法提取的圖像輪廓誤差較大。

基于以上分析，本文進行了一些改進。提出一種自適應掃描線逼近法，具體算法步驟如下。

⑴ 為突出邊緣信息，將采集的魚眼圖像做銳化處理，并計算出平均閾值T。

⑵ 將銳化后的圖像從上到下按水平方向掃描，如果像素點（i，j）值小于T，則繼續掃描。

⑶ 如果像素點（i，j）的值大于或等于閾值T，則判斷以改點為中心的3×3領域內其他像素點的值是否存在全為零的情況。

⑷ 如果存在，則改點為干擾點，故繼續掃描；否則y=j為圖像的上水平切線top。

⑸ 同理，重復步驟⑵至步驟⑷，分別確定圖像的下水平切線bottom，以及左垂直切線left、右垂直切線right。

⑹ 計算半徑R和圓心坐標（x0，y0）。

⑷

2 球面等距投影校正及展開模型

球面等距投影[6-8]是一種非線性投影，其描述的是球面物與平面像的映射關系。其投影原理是，球面物上的像素點如果具有相同的視場角，那么其在平面像上對應相等的徑向距離r'。即用公式表示為：

⑸

設XOY平面為魚眼圖像的成像平面，那么將該像平面上的每一個像素點P1校正為空間上的點P3需分成三步：第一步，首先將像平面點非線性的逆投影到球面上，而這一過程必須滿足等距投影逆變化的約束；第二步，每一個球面點P2被映射為連接P2與球心O的射線OP2；第三步，將射線OP2線性的映射到校正平面上。

具體算法步驟如下。

⑴ 計算魚眼圖像上的點到圖像圓心的距離。設點P1為魚眼圖像的像素點，其坐標為P1=（x，y），求出點P1到圓心的距離r并歸一化。

⑹

⑵ 計算球平面任一點P2的方位角和仰角。根據等距投影成像原理，那么等距投影校正模型就是其成像原理的逆推過程。

P2為魚眼像平面的點P1映射到半球面上的點，其坐標用表示（θ為方位角，為俯仰角），其計算公式如下：

⑺

⑻

⑶ 將球面上的任意像素點，線性地投影到校正平面上。每一個球面上的點P2被映射為連接圓心0與P2的射線OP2交立方體面于點P3，如圖1所示。

取Z=1作為校正目標平面時，那么球面上像素點投影到目標平面上點的坐標步驟如下：

首先，計算點P3到目標平面中心的距離 。

其次，計算的坐標。；。

當仰角較大時，較大，特別是近似90?時，趨于∞，即校正目標平面趨于無窮大。為改進這一不足，基于投影物球面建立一個以球心為中心的正切立方體，將球面上的像素點投影到立方體的上、前、后、左、右五個面上，如圖2所示。

本文取=45?作為臨界點。仰角∈[0，45?]、方位角θ∈[-180?，180?]區間內的球面像素點投影到Z=1的平面上；仰角∈[45?，90?]、方位角θ∈[-180?，180?]區間內的球面像素點分別投影到立方體四周的投影面上。

3 實驗結果展示及結論

本文在matlab平臺下，采集180?左右格式為.gif的魚眼圖像（如圖3所示）后首先提取其有效區域，接著對有效區域做歸一化處理，然后通過等距投影逆推原理將有效區域的魚眼原始圖像投影到單位為1的半球面，最后將半球面的像素點一一映射到校正平面上，完成圖像校正。圖4為提取的魚眼輪廓圖像。圖5為采集該圖像0?～90?部分并對其校正后的圖像，圖6為該圖像0?～180?的校正展開圖像。

從實驗結果看，第一，本文提出的掃描算法對魚眼圖像有效區域的提取有較高的抗干擾能力；第二，魚眼圖像0?～90?部分校正效果較好，隨著視場的增大，當視場接近180?時，略微出現畸變情況。

相對于多項式魚眼圖像畸變校正算法，本文的算法優勢是不需要多項式計算，算法復雜度簡單，運行速度快且效果較好，比較適應于在硬件平臺轉化為工程算法。在接下來的研究中，如果進一步改善圖像視覺質量和提高算法在視場接近180?時曲線校正精度是工作重點之一。

參考文獻（References）：

[1] 楊玲，成運.應用經緯映射的魚眼圖像校正設計方法[J].工程

圖學學報，2010.31（6）：19-22

[2] 廖士中，高培煥，蘇藝，王大鵬.一種光學鏡頭攝像機圖像幾何

畸變的修正方法[J].中國圖像圖形學報，2000.5（7）：593-596

[3] 英向華，胡占義.一種基于球面透視投影約束的魚眼鏡頭校

正方法[J].計算機學報，2003.26（12）：1702-1708

[4] 王大宇，崔漢國，陳軍.魚眼圖像輪廓提取及校正研究[J].計算

機工程與設計，2007.28（12）：2878-2880

[5] 張軍，王志舟，楊正瓴.單幅圓形魚眼圖像校正[J]. 計算機應

用，2015.35（5）：1444-1448

[6] 楊秋慧，玉昊昕，李陶深，玉振明.一種大視場鏡頭的畸變校正

方法[J].計算機應用與軟件，2015.32（09）：136-140

[7] 孟靈.等距投影的魚眼圖像校正算法設計[J].數字技術與應

用，2016.7（6）：125-126

[8] 吳桂平，吳魏，王成，畢昆.基于雙線性插值魚眼圖像校正算法[J].

計算機應用與軟件，2012.29（2）：122-124