孫佳明 余俊鵬 李潔明 滿益云 沈剛

一種通用攝影測量畸變模型性能仿真分析

孫佳明1余俊鵬1李潔明1滿益云2沈剛1

（1 廣東工業(yè)大學土木與交通工程學院，廣州 510006）（2 錢學森空間技術實驗室，北京 100094）

攝影測量中，畸變模型的選擇對于保證畸變檢校效果十分重要。文章提出了一種基于二元傅里葉級數的通用型畸變模型，理論上可精確擬合各種形式的幾何畸變。對模擬畸變數據的擬合試驗結果表明，傅里葉模型相比一般多項式模型具有更優(yōu)的畸變擬合性能，擬合精度隨著階數的增加而穩(wěn)步提高，對Brown模型畸變擬合中誤差可達到0.3像元以內，對Ebner模型畸變擬合中誤差達到0.001像元的水平。二元傅里葉模型可在攝影測量數據處理和光學相機檢校中替代其他畸變模型，作為一種通用畸變模型使用。

攝影測量 傅里葉級數 畸變模型 仿真分析 航天測繪

0 引言

攝影測量畸變指由于攝影成像系統(tǒng)在視場范圍內的成像光束未嚴格滿足理想中心投影而產生的像點位置偏差，反映了鏡頭光學畸變、感光器件形變、光線折射等因素的綜合作用?；兊臋z校和補償對于保證攝影測量精度至關重要[1-4]。

畸變作為一種系統(tǒng)誤差，其誤差模型大致可分為物理模型和數學模型兩大類[5-6]。1965年，美國攝影測量學者Duane C．Brown通過大量試驗，對光學相機幾何畸變的物理因素深入分析后提出了Brown模型[7]。至今，Brown模型及其改進模型仍是近景和航空攝影測量中被普遍采用的“標準”模型[8-12]。1976年，Ebner提出一種數學模型——正交多項式模型，采用一組正交參數集對攝影測量畸變進行數學擬合，并通過自檢校空中三角網光束法平差驗證了其有效性[13]。與正交多項式模型類似的數學模型還有一般多項式模型[14]、球面調和函數模型等。

在攝影測量科研和生產實踐過程中，不同畸變模型往往根據某種特定情形而提出[15-19]，缺乏較好的通用性。在實際數據處理過程中，很難自動選擇最適合的模型進行畸變檢校和補償，從而導致處理結果的精度損失。針對這一問題，本文基于傅里葉級數對連續(xù)函數的逼近特性，提出一種可適應任意畸變的通用攝影測量畸變模型，通過仿真試驗表明，其畸變擬合性能隨著選用階數增加而提高，可望應用于航天航空遙感影像的高精度幾何處理。

1 目標模型

1.1 攝影測量畸變

式中1、2分別為、方向的畸變函數。

式（3）中函數關系可有多種選擇方案，以下分別給出Brown模型和Ebner模型兩種代表性模型。

1.2 Brown模型

Brown模型是一種模擬物理因素作用的經典畸變模型[24]，形式如下：

1.3 Ebner模型

Ebner模型為一種數學畸變模型。在傳統(tǒng)攝影測量的9個標準點位上，將畸變量表示為

式中1~12表示、方向的縮放、旋轉、錯切、彎曲等基本幾何形變項的系數；為標準點位距。各項具體形式可參考文獻[23]。在傳統(tǒng)攝影測量的9個標準點位上，各項參數之間具備正交性，有利于提高求解穩(wěn)定性和進行參數顯著性檢驗。然而，隨著各種數字傳感器的應用，像方不規(guī)則點位分布情況十分普遍，Ebner模型的適用性有所下降。

2 通用畸變模型

數學畸變模型不必考慮參數的物理意義，適合作為通用型畸變模型使用。一般多項式模型是目前較為常用的通用模型，其形式為：

式中為多項式階數；、分別為各項中、的冪次數；p,n，q,n為多項式各項對應的系數。

隨著一般多項式階數的增加，可以提升其對于特定目標函數的擬合精度，但提升幅度有限，無法實現任意程度逼近。此外實際應用中也不宜過大（一般不超過3），否則容易出現參數過多化而難以求解。因此，尋求一種能夠精確擬合各種幾何畸變并且便于求解的通用模型，仍然十分必要。

以上各項系數之間嚴格正交，有利于實現精確求解。以式（7）為基礎，得到用二元傅里葉級數逼近像方未知幾何畸變的形式，即以像點坐標表示的二元傅里葉畸變模型

二元傅里葉多項式的階數越高，對畸變的擬合精度越高。當=1時，傅里葉畸變模型簡化為常差模型：

當=2時，式（8）所確定的傅里葉畸變模型含16個待求參數，具體形式為

當=3時，畸變模型待求參數增加至36個，形式為

3 試驗及結果分析

3.1 試驗方法

基于上述理論，可通過仿真試驗方法對傅里葉畸變模型的性能進行驗證。試驗流程如圖1所示。首先，在10 000像元×10 000像元的影像范圍內，分別以Brown模型和Ebner模型作為目標模型，生成20×20=400個規(guī)則格網點位上的模擬畸變值；將各個格網點位坐標與其畸變相加得到模擬像點坐標觀測值；然后，以各個點位的模擬畸變值作為觀測值，將模擬像點坐標觀測值代入式（8），建立觀測方程組；按最小二乘原理構造法方程求解參數，統(tǒng)計各點位殘差，得到傅里葉模型對目標模型的擬合精度結果。

圖1 試驗方法流程

3.2 試驗結果

3.2.1 對Brown模型的擬合結果

對由Brown模型生成的不同大小畸變，分別用3~7階傅里葉模型進行擬合，結果見表1。

對不同大小的Brown模型畸變，分別用3階一般多項式模型和3階傅里葉模型進行擬合。結果如圖2所示。當邊緣畸變大小為20像元時，一般多項式擬合中誤差大于1個像元，而傅里葉模型擬合中誤差為0.8像元。

表1 對Brown模型擬合結果

Tab.1 The fitting results to the Brown model 單位：像元

圖2 對Brown模型畸變擬合精度對比

3.2.2 對Ebner模型的擬合結果

對由Ebner模型生成的不同大小畸變，用3~7階傅里葉模型分別擬合，結果見表2。

對由Ebner模型生成的邊緣畸變大小為20像元的情況，畸變擬合前后殘差分布如圖3所示?？梢钥闯觯洈M合后畸變的誤差系統(tǒng)性基本消除，擬合效果良好。

3.3 試驗分析

試驗結果表明，對于Brown模型畸變，當邊緣畸變大小分別為5、10、20像元時，采用3階傅里葉模型擬合，中誤差分別為0.2、0.4、0.8個像元，均優(yōu)于1個像元。采用7階傅里葉模型，擬合中誤差減小到0.1個像元以內，最大殘差0.89像元。對比傅里葉模型和一般多項式模型的擬合效果，兩者在采用相同階數時傅里葉模型的擬合精度更優(yōu)，且畸變較大時精度優(yōu)勢更加明顯。

對于Ebner模型畸變，當邊緣畸變大小分別為5、10、20像元時，3階傅里葉模型的擬合中誤差分別為0.07、0.14、0.23個像元，均優(yōu)于0.3個像元，且在行列方向上的殘差及中誤差大小基本一致。采用7階傅里葉模型，擬合中誤差已達到0.001個像元的水平。

從試驗結果可以看出，無論對于Brown模型還是Ebner模型，傅里葉模型的擬合精度都隨著階數增加而不斷提高，與其理論特性相符合。

表2 對Ebner模型擬合結果

Tab.2 Fitting results to the Ebner model 單位：像元

圖3 對Ebner模型殘差擬合前后對比

4 結束語

在攝影測量處理中，由于模擬物理因素作用的畸變模型形式趨于復雜，模型參數與攝影測量定向參數的強相關性導致其解算困難。本文提出一種基于二元傅里葉多項式的數學畸變模型，通過仿真擬合試驗表明，3階傅里葉模型對Brown模型的擬合精度小于1個像元，7階模型的擬合精度可小于0.1個像元；7階傅里葉模型對Ebner模型的擬合精度能夠達到千分之一像元的水平。在相同階數情況下，傅里葉模型對畸變擬合效果明顯優(yōu)于一般多項式通用模型。當選取足夠高的階數時，二元傅里葉函數理論上可以任意精度逼近多種畸變函數形式，是一種較為理想的通用模型。后續(xù)將開展傅里葉畸變模型在光學相機檢校、航空航天攝影測量數據處理等方面的實驗研究，進一步考察其實用性能。

[1] 汪佳麗, 王煒杰, 馬振玲. 相機畸變的混合模型迭代檢校法[J]. 測繪通報, 2019(4): 111-114.WANG Jiali,WANG Weijie, MA Zhenling. Hybrid-model Based Camera Distortion Iterative Calibration Method[J]. Bulletin of Surveying and Mapping, 2019(4): 111-114. (in Chinese)

[2] 崔璨璨, 汪佳麗, 馬振玲. 相機檢校的迭代處理方法[J]. 測繪通報, 2018(6): 50-54, 60.CUI Cancan, WANG Jiali, MA Zhenling. The Iterative Method of Camera Calibration[J]. Bulletin of Surveying and Mapping, 2018(6): 50-54, 60. (in Chinese)

[3] 徐文, 龍小祥, 李慶鵬, 等. “高分四號”衛(wèi)星影像輻射與幾何精度評價[J]. 航天返回與遙感, 2016, 37(4): 16-25.XU Wen, LONG Xiaoxiang, LI Qingpeng, et al. Image Radiometric and Geometric Accuracy Evaluation of GF-4 Satellite[J]. Spacecraft Recovery & Remote Sensing, 2016, 37(4): 16-25. (in Chinese)

[4] 朱瑤, 王治樂, 賀磊, 等. 線掃描相機標定及畸變矯正方法[J]. 航天返回與遙感, 2017, 38(6): 46-53.ZHU Yao, WANG Zhile, HE Lei, et al. Calibration and Distortion Correction Method of Line Scan Camera[J]. Spacecraft Recovery & Remote Sensing, 2017, 38(6): 46-53. (in Chinese)

[5] 黃留波. 不同畸變模型在航攝相機校準中的適用性研究[J]. 廣西水利水電, 2016(4): 24-27. HUANG Liubo. Study of Adaptability of Different Distortion Model in Aerial Photographic Camera Calibration[J]. Guangxi Water Resources & Hydropower Engineering, 2016(4): 24-27. (in Chinese)

[6] 陳國平. 攝影測量與遙感[M]. 武漢: 武漢大學出版社, 2011. CHEN Guoping. Photogrammetry and Remote Sensing[M]. Wuhan: Wuhan University Press, 2011. (in Chinese)

[7] BROWN D C. Close-range Camera Calibration[J]. Photogrammetric Engineering, 1971, 37(8): 855-866.

[8] 楊韞瀾, 韓玲, 胡海彥. Brown自檢校模型相關性分析與改進[J]. 測繪通報, 2014(5): 37-40.YANG Yunlan, HAN Ling, HU Haiyan. Correlation Analysis and Improvement of Brown Self-calibration Model[J]. Bulletin of Surveying and Mapping, 2014(5): 37-40. (in Chinese)

[9] 孫鵬. 大尺寸動態(tài)攝影測量關鍵技術研究[D]. 北京: 北京郵電大學, 2019. SUN Peng. Research on Key Techniques of Large Scale Dynamic Photogrammetry[D]. Beijing: Beijing University of Posts and Telecommunications, 2019. (in Chinese)

[10] 曹宗偉, 李星全, 杜沛遠. 數字攝影測量相機的畸變測量與糾正[J]. 測繪工程, 2017, 26(3): 60-64. CAO Zongwei, LI Xingquan, DU Peiyuan. Distortion Measurement and Correction of Digital Photogrammetry Camera[J]. Engineering of Surveying and Mapping , 2017, 26(3): 60-64. (in Chinese)

[11] 馮其強. 數字工業(yè)攝影測量技術研究與實踐[D]. 鄭州: 解放軍信息工程大學, 2010. FENG Qiqiang. Research and Practice of Digital Industrial Photogrammetry[D]. Zhengzhou: PLA Information Engineering University, 2010. (in Chinese)

[12] 馮文灝, 商浩亮, 侯文廣. 影像的數字畸變模型[J]. 武漢大學學報(信息科學版), 2006, 31(2): 99-103. FENG Wenhao, SHANG Haoliang, HOU Wenguang. A Digital Distortion Model for all Kinds of Imaging Systems[J]. Geomatics and Information Science of Wuhan University, 2006, 31(2): 99-103. (in Chinese)

[13] EBNER H. Self-calibrating Block Adjustment[J]. Bildmessung and Luftbildwesen, 1976, 44(4): 128-139.

[14] 楊韞瀾. 攝影測量自檢校技術研究[D]. 西安: 長安大學, 2015. YANG Yunlan. Study on the Self-calibration Technology in Photogrammetry[D]. Xi'an: Chang'an University, 2015. (in Chinese)

[15] 劉世杰, 唐新明, 王鳳香, 等. 基于遙感影像的衛(wèi)星顫振探測分析的關鍵技術與實現[J]. 地球信息科學學報, 2015, 17(6): 724-731. LIU Shijie, TANG Xinming, WANG Fengxiang, et al. Research on Key Technology and Implementation of Detection and Analysis of Satellite Attitude Jitter Using Remote Sensing Imagery[J]. Journal of Geo-Information Science, 2015, 17(6): 724-731. (in Chinese)

[16] 童小華, 葉真, 劉世杰. 高分辨率衛(wèi)星顫振探測補償的關鍵技術方法與應用[J]. 測繪學報, 2017, 46(10): 1500-1508. TONG Xiaohua, YE Zhen, LIU Shijie. Essential Technology and Application of Jitter Detection and Compensation for High Resolution Satellites[J].Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1500-1508. (in Chinese)

[17] 胡堃, 黃旭, 張永軍, 等. 基于高頻角位移數據的衛(wèi)星平臺顫振檢測與影像幾何質量補償[J]. 電子與信息學報, 2018, 40(7): 6-12. HU Kun, HUANG Xu, ZHANG Yongjun, et al. Satellite Platform Jitter Detection and Image Geometric Quality Compensation Based on High-frequency Angular Displacement Data[J]. Journal of Electronics and Information Technology, 2018, 40(7): 6-12. (in Chinese)

[18] 張宏偉, 張炳先, 侯作勛, 等. 顧及光行差改正的遙感衛(wèi)星成像模型及驗證[J]. 航天返回與遙感, 2019, 40(4): 112-121. ZHANG Hongwei, ZHANG Bingxian, HOU Zuoxun, et al. Construction and Verification of Imaging Model for Optical Remote Sensing Satellite Considering the Correction of Aberration of Light[J]. Spacecraft Recovery & Remote Sensing, 2019, 40(4): 112-121. (in Chinese)

[19] 王密, 朱映, 范城城. 高分辨率光學衛(wèi)星影像平臺震顫幾何精度影響分析與處理研究綜述[J]. 武漢大學學報(信息科學版), 2018, 43(12): 146-155. WANG Mi, ZHU Ying, FAN Chengcheng.Development of Platform Jitter Geometric Analysis and Processing for High-Resolution Optical Satellite Imagery[J]. Journal of Wuhan University (Information Science Edition), 2018, 43(12): 146-155. (in Chinese)

[20] 余俊鵬, 何熾榮. GPS-RTK輔助近景攝影測量光束法平差及精度分析[J]. 廣東工業(yè)大學學報, 2017, 34(6): 73-77. YU Junpeng, HE Chirong. An Accuracy Analysis of GPS-RTK Assisted Close-Range Photogrammetric Bundle Adjustment[J]. Journal of Guangdong University of Technology, 2017, 34(6): 73-77. (in Chinese)

[21] 姚娜, 林招榮, 任超鋒, 等. 一種適合攝像機非線性畸變校正的畸變模型[J]. 激光與光電子學進展, 2014, 51(2): 167-174. YAO Na, LIN Zhaorong, REN Chaofeng, et al. A Distortion Model Suitable for Nonlinear Distortion Correction of Digital Video Camera[J]. Laser & Optoelectronics Progress, 2014, 51(2): 167-174. (in Chinese)

[22] 蔣永華, 張過, 唐新明, 等. 資源三號測繪衛(wèi)星三線陣影像高精度幾何檢校[J]. 測繪學報, 2013, 42(4): 523-529, 553.JIANG Yonghua, ZHANG Guo, TANG Xinming, et al. High Accuracy Geometric Calibration of ZY-3 Three-line Image[J]. Acta Geodaetica et Cartographica Sinica, 2013, 42(4): 523-529,553. (in Chinese)

[23] 李德仁, 袁修孝. 誤差處理與可靠性理論[M]. 2版. 武漢: 武漢大學出版社, 2012. LI Deren, YUAN Xiuxiao. Error Processing and Reliability Theory[M]. 2nd ed. Wuhan: Wuhan University Press, 2012. (in Chinese)

[24] 高瞻宇. 采用簡化Brown模型及改進BFGS法的相機自標定[J]. 光學精密工程, 2017, 25(9): 2532-2540. GAO Zhanyu. Self-calibration Based on Simplified Brown Non-linear Camera Model and Modified BFGS Algorithm[J]. Optics and Precision Engineering, 2017, 25(9): 2532-2540. (in Chinese)

[25] TANG R, FRITSCH D, CRAMER M. New Rigorous and Flexible Fourier Self-calibration Models for Airborne Camera Calibration[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2012, 71: 76-85.

Performance Analysis of a Generic Photogrammetric Distortion Model

SUN Jiaming1YU Junpeng1LI Jieming1MAN Yiyun2SHEN Gang1

（1 School of Civil and Traffic Engineering, Guangdong University of Technology, Guangzhou 510006, China）（2 Qian Xuesen Laboratory of Space Technology, Beijing 100094, China）

The selection of distortion models in photogrammetry is very important to ensure distortion calibration effect. In this paper, a generic distortion model based on binary Fourier series is proposed, which can approximate various forms of geometric distortion with high precision theoretically. The fitting test results on simulated distortion data show that the Fourier model has better distortion fitting performance than the general polynomial model, the fitting accuracy steadily improves with the increase of the order, and the fitting error to the Brown distortion can reach within 0.3 pixel, the fitting error to the Ebner model reaches the level of 0.001 pixel. The binary Fourier model is expected to replace other distortion models in photogrammetric data processing and optical camera calibration as a generic distortion model.

photogrammetry; Fourier series; distortion model; simulation analysis; space mapping

P237

A

1009-8518(2020)05-0110-08

10.3969/j.issn.1009-8518.2020.05.013

2019-12-12

國家自然科學基金青年基金項目（41704019）

孫佳明, 余俊鵬, 李潔明, 等. 一種通用攝影測量畸變模型性能仿真分析[J]. 航天返回與遙感, 2020, 41(5): 110-117.

SUN Jiaming, YU Junpeng, LI Jieming, et al. Performance Analysis of a Generic Photogrammetric Distortion Model[J]. Spacecraft Recovery & Remote Sensing, 2020, 41(5): 110-117. (in Chinese)

孫佳明，男，1995年生，廣東工業(yè)大學測繪科學與技術專業(yè)在讀碩士研究生，主要研究方向為攝影測量與遙感。E-mail：1164783700@qq.com。

余俊鵬，男，1982年生，2009年獲武漢大學攝影測量與遙感專業(yè)博士學位，高級工程師。主要研究方向為攝影測量與遙感。E-mail：strength701@163.com。

（編輯：夏淑密）