曹 斌，朱述龍，邱振戈，曹彬才

1. 上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306; 2. 信息工程大學(xué)五院，河南 鄭州 450001

曹斌，朱述龍，邱振戈，等.一種更嚴密的雙介質(zhì)立體攝影測量折射改正算法[J].測繪學(xué)報，2017,46(9)：1182-1192.

10.11947/j.AGCS.2017.20170119.

CAO Bin，ZHU Shulong，QIU Zhenge，et al.More Rigorous Correction of Refraction Effects in Two-media Stereophoto-grammetry[J]. Acta Geodaetica et Cartographica Sinica,2017,46(9):1182-1192. DOI:10.11947/j.AGCS.2017.20170119.

一種更嚴密的雙介質(zhì)立體攝影測量折射改正算法

曹 斌1，朱述龍2，邱振戈1，曹彬才2

1. 上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306; 2. 信息工程大學(xué)五院，河南 鄭州 450001

提出一種更嚴密的雙介質(zhì)立體攝影測量物方坐標折射改正算法。該算法用水下目標的空中同名直線光線公垂線段的中點作為攝影測量交會點的理論位置，解決了空中同名直線光線延長線不相交情況下攝影測量交會點不存在導(dǎo)致的點位關(guān)系不確定性問題，使攝影測量交會點到真實物點的坐標折射改正公式能夠嚴格推導(dǎo)出來。分析了空中同名直線光線不相交情況對物方大地坐標折射改正的不利影響，研究了攝影測量交會點與水下真實物點的相互位置關(guān)系，推導(dǎo)了水下目標點的水深和大地坐標計算公式(即折射改正公式)，通過WorldView-2立體影像淺海海底地形測量試驗對算法的正確性和測量精度進行驗證。研究表明，不論水下目標的空中同名直線光線延長線是否相交，該算法都是適用的，且能顯著改善水下目標的高程測量精度。

雙介質(zhì)立體攝影測量；折射改正；攝影測量交會點；真實物點；水下目標定位

在航空/航天光學(xué)遙感中，如果觀測對象為水下目標，則來自目標的電磁波要穿過水和空氣兩種介質(zhì)才能到達相機焦平面成像，此時可運用雙介質(zhì)立體攝影測量方法獲取水下目標點的三維大地坐標[1-5]。

成像光線在單介質(zhì)(水或空氣)中可視為一條直線，在雙介質(zhì)(水和空氣)中可視為一條折線。在雙介質(zhì)立體攝影測量中，如果能直接建立水下目標點與對應(yīng)像點間的折線關(guān)系模型，則利用立體影像上的同名像點，按照物像折線關(guān)系模型就可以直接交會得到水下目標點的三維大地坐標[6-11]。但在實際應(yīng)用中，通常僅有遙感圖像的單介質(zhì)空中物像直線關(guān)系或其近似模型，此時對于水下目標點，只能先利用立體影像上的相應(yīng)同名像點，按照空中物像直線關(guān)系或其近似模型，聯(lián)立求解出虛擬的水中點位(不妨也稱為“攝影測量交會點”)的三維大地坐標，然后再按照光線在空氣和水中的折線傳播規(guī)律，對攝影測量交會點的三維大地坐標進行折光改正，進而得到水下目標點的三維大地坐標[12-19]。本文主要研究后一種情況的折射改正問題。

現(xiàn)有的研究成果[12-19]幾乎都提到了“在雙介質(zhì)成像情況下水下目標的空中同名直線光線不相交情況的存在”。但是在推導(dǎo)折射改正公式時，卻都忽略了空中同名直線光線不相交情況對物方坐標折射改正的影響，因而得到的折射改正公式具有很大的近似性或不合理性，由此得到的水下目標點的三維大地坐標也必然存在不合理誤差。

針對上述問題，首先根據(jù)水下目標點的雙介質(zhì)立體成像幾何特點，分析空中同名直線光線不相交情況對物方大地坐標折射改正的不利影響，然后用水下目標的空中同名直線光線公垂線段的中點作為攝影測量交會點的理論位置，在此基礎(chǔ)上建立攝影測量交會點與水下真實物點的相互位置關(guān)系，推導(dǎo)出水下目標點的水深和大地坐標計算公式(即折射改正公式)。為了驗證所提算法的正確性和有效性，開發(fā)了相應(yīng)的雙介質(zhì)立體攝影測量軟件，并利用WorldView-2衛(wèi)星遙感影像及相關(guān)數(shù)據(jù)進行了算法的合理性試驗、理論誤差測算和精度改善情況試驗，取得了預(yù)期的結(jié)果。

1 空中同名光線延長線不相交對折射改正的影響

如圖1所示，P為水下目標點，P在立體像對左、右圖像上的像點分別為p1、p2；在雙介質(zhì)成像情況下，左圖像上像點p1的成像光線為折線S1P1P(其中像點p1在直線S1P1上，S1為獲取左圖像時的相機鏡頭中心位置，P1為成像光線與空氣和水分界面的交點)；右圖像上同名像點p2的成像光線為折線S2P2P(其中像點p2在直線S2P2上，S2為獲取右圖像時的相機鏡頭中心位置，P2為成像光線與空氣和水分界面的交點)。如果僅有遙感圖像的單介質(zhì)空中物像直線關(guān)系模型，則可先按照單介質(zhì)立體攝影測量原理，利用同名像點p1、p2的像坐標，按照直線S1P1和S2P2的數(shù)學(xué)模型，聯(lián)立求解出攝影測量交會點F的大地坐標(圖1中F為直線S1P1和S2P2延長線的交點)；然后再按照光在空氣和水中的折線傳播規(guī)律，對F的大地坐標進行折光改正，進而將F修正為目標點P。

圖1 雙介質(zhì)立體攝影測量原理平面示意圖Fig.1 Sketch of two-media stereophotogrammetry

實際上，在雙介質(zhì)成像情況下，三維空間中的直線S1P1和S2P2的延長線并不一定相交。如圖2所示，直線S1P1的延長線在折線S1P1P確定的垂直于水面的平面內(nèi)，直線S2P2的延長線在折線S2P2P確定的垂直于水面的平面內(nèi)。如果這兩個平面不重合且直線S1P1和S2P2關(guān)于L不對稱(L為兩個平面的交線，L經(jīng)過目標點P且垂直于水面)，則兩條空間異面直線S1P1和S2P2的延長線就不會相交。只有這兩個平面重合(即水下目標點P位于兩個攝站連線的正下方，參見圖3)或者成像光線關(guān)于L對稱(參見圖4)時，直線S1P1和S2P2的延長線才會相交，且交點一定位于L上。

圖2 空中同名光線延長線不相交情況Fig.2 The case of non-intersection of two aerial corresponding rays

圖3 空中同名光線延長線相交情形之一(空中同名光線所在平面重合)Fig.3 A case of intersection of two aerial corresponding rays (the two rays are in the same plane)

在雙介質(zhì)攝影測量中，直線S1P1和S2P2的延長線不相交情況是大量存在的，而相交情況只是特例。在直線S1P1和S2P2延長線不相交情況下，按照單介質(zhì)立體攝影測量原理求解出的攝影測量交會點F的理論位置，是決定物方折射改正的最主要依據(jù)。雖然文獻[12-19]都注意到了空中同名直線光線不相交情況的存在，但都沒有回答這種情況下攝影測量交會點F理論位置究竟在哪里？所以不可能從根本上解決不相交情況對折射改正的影響。例如，文獻[12—13]只是近似地按照直線S1P1和S2P2延長線相交時的幾何關(guān)系進行物方坐標折射改正(即按照圖3和圖4所示的幾何關(guān)系進行物方坐標折射改正)，文獻[14—19]僅按圖5所示的單像幾何關(guān)系進行物方坐標折射改正，其不合理性是顯而易見的，這樣做的結(jié)果必然導(dǎo)致明顯的折射改正誤差。為克服現(xiàn)有算法不足，本文從確定攝影測量交會點F的理論位置著手研究更嚴密的物方坐標折射改正算法。

圖4 空中同名光線延長線相交情形之二(空中同名光線關(guān)于L對稱)Fig.4 A case of intersection of two aerial corresponding rays (the two rays are symmetric with respect to L)

2 基于嚴密物像幾何關(guān)系的折光改正算法

不失一般性，假設(shè)L與水面的交點為M，直線S1P1的延長線與L的交點為F1，直線S2P2的延長線與L的交點為F2，按照直線S1P1和S2P2的數(shù)學(xué)模型聯(lián)立求解出的點位為F，各點的相對關(guān)系如圖6所示。理論上，F(xiàn)的最優(yōu)解應(yīng)為兩條空間異面直線S1P1和S2P2的公垂線段的中點。當S1P1和S2P2不相交時，F(xiàn)不在L上(注意P、F2、F1和M都在L上)。特殊情況下，當S1P1和S2P2的延長線相交時，F(xiàn)、F2和F1重合為一點且該點一定在L上。

圖5 單幅影像的物像幾何關(guān)系Fig.5 Geometry of two-media photogrammetry using a single image

圖6 雙介質(zhì)立體攝影測量主要點位之間的關(guān)系 Fig.6 Geometry of two-media stereophotogrammetry (including main photogrammetric points)

(1)

(2)

(3)

式中

(4)

(5)

式中，i1為左像點p1對應(yīng)的水下目標點P的輻射光線照射到水面點P1的入射角；r1為左像點p1對應(yīng)的水下目標點P的輻射光線穿過水面點P1的折射角；i2為右像點p2對應(yīng)的水下目標點P的輻射光線照射到水面點P2的入射角；r2為右像點p2對應(yīng)的水下目標點P的輻射光線穿過水面點P2的折射角；θ為空中同名直線光線S1P1和S2P2的夾角(即交會角)。

利用直線光線S1P1和直線光線S2P2的單位方向向量，可推導(dǎo)出θ的表達式如下

θ=cos-1(cosω1sinα1cosω2sinα2+

sinω1sinω2+cosω1cosα1cosω2cosα2)

(6)

式中，α1為左像點p1對應(yīng)的空中直線光線S1P1的前視角；ω1為左像點p1對應(yīng)的空中直線光線S1P1的側(cè)視角；α2為右像點p2對應(yīng)的空中直線光線S2P2的后視角；α2為右像點p2對應(yīng)的空中直線光線S2P2的側(cè)視角。

(7)

所以F到線段F1F2中點的距離δ為

(8)

式中，l1、l2分別為線段F1P1和F2P2的長度。

假設(shè)F1、F2、P的水深(即圖6中線段F1M、F2M、PM的長度)分別為h1、h2、hP，則由三角形幾何關(guān)系可得

(9)

(10)

(11)

(12)

在雙介質(zhì)光學(xué)遙感中，同名像點p1、p2對應(yīng)的r1、r2一般是已知的。由折射定律可知

(13)

(14)

式中,n是水相對空氣的折射率(在水溫為20°C時，n=1.332 99)。

利用WorldView-2衛(wèi)星遙感影像及其觀測參數(shù)實際驗算的結(jié)果表明：F到線段F1F2中點的距離δ并不大(參見表2)，相較于雙介質(zhì)攝影測量能力而言，可以近似地將攝影測量交會點F和線段F1F2中點視為同一點。由此可得：線段F1F2中點的大地坐標近似等于F的大地坐標，而F的水深hF可用F1F2中點的水深來代替，即

(15)

整理可得P的水深hP如下

(16)

hF=Z0-ZF

(17)

將式(17)代入式(16)可得

(18)

考慮到P與線段F1F2中點的平面坐標相同，可得P的大地坐標為

(19)

3 試驗與分析

下面通過兩個雙介質(zhì)立體攝影測量試驗來驗證本文算法的合理性以及對高程測量精度的影響。一是利用WorldView-2立體影像的觀測參數(shù)計算攝影測量交會點F到線段F1F2中點的距離δ(各符號含義參見圖6)，并根據(jù)計算結(jié)果來驗證本文算法的合理性以及由此引起的理論誤差大小；二是利用WorldView-2立體影像及其觀測參數(shù)進行淺海海底地形測量試驗，并將本文算法的結(jié)果與最有代表性的Murase算法[13]的結(jié)果進行比較，來驗證本文算法對高程測量精度的改善能力。

3.1 試驗衛(wèi)星影像及其特點

綜合考慮試驗需要和遙感數(shù)據(jù)保障條件，選擇海南省三沙市甘泉島和珊瑚島兩個地區(qū)的WorldView-2[20]多光譜立體影像及其觀測參數(shù)和定位模型參數(shù)進行試驗。圖7和圖8分別為甘泉島和珊瑚島地區(qū)WorldView-2多光譜立體影像的局部重疊區(qū)域，而整幅圖像對應(yīng)的航向視角(即前視角或后視角)、側(cè)視角、入視角等觀測參數(shù)如表1所示。實際試驗中，圖像上每個像素對應(yīng)的航向視角、側(cè)視角、入視角等觀測參數(shù)用雙線性內(nèi)插方法計算得到。

雖然WorldView-2提供了全色立體影像和多光譜立體影像，但由于多光譜影像的分辨率已經(jīng)很高(星下點分辨率已達1.85 m)，完全能夠滿足攝影測量需要，并且其中藍、綠波段影像的海底紋理特征比其他波段(包括全色)更清晰[21-22]，更有利于海底紋理特征的檢測和提取，這是本文選用WorldView-2多光譜立體影像的主要原因。

表1 WorldView-2多光譜立體影像的觀測角

說明：航向視角即為文中的α1或α2(正值表示前視，負值表示后視)，側(cè)視角即為文中的ω1或ω2(正值表示右視，負值表示左視)，入視角即為文中的r1或r2。

由于甘泉島、珊瑚島兩個地區(qū)的WorldView-2立體影像的左圖像都有嚴重的太陽耀斑，為了消除太陽耀光對水域影像質(zhì)量、海底紋理特征提取和水下同名像點匹配等后續(xù)工作的影響，運用文獻[23]的太陽耀斑消除的方法對甘泉島、珊瑚島兩個地區(qū)WorldView-2立體影像的左圖像進行了處理。

由于WorldView-2相機采用了圖9所示的交錯拼接型長線陣探測器[24]，所以WorldView-2影像(包括全色和多光譜)是由多個短線陣探測器推掃成像后經(jīng)過拼接處理得到的，已不是嚴格意義上的原始影像，所以無法建立嚴格的物像關(guān)系模型。為此，WorldView-2提供了有理函數(shù)模型RFM(rational function model)以及相應(yīng)的有理多項式系數(shù)RPC(rational polynomial coefficients)和坐標歸一化參數(shù)作為近似的目標定位模型[25]。該模型考慮了光學(xué)相機鏡頭畸變、光行差(由衛(wèi)星運動導(dǎo)致)、大氣折射的影響，其表達式如下

(20)

在實際應(yīng)用中，通常對式(20)進行線性化并整理得到如下誤差方程式

(21)

對于一對同名像點，依據(jù)式(21)可以列出4個方程式，聯(lián)立這4個方程式就可以求解出該同名像點對應(yīng)的地面點的大地坐標(X,Y,Z)。

3.2 算法的合理性試驗及理論誤差測算

本文算法的關(guān)鍵是利用水下目標的空中同名直線光線公垂線段的中點作為攝影測量交會點F的理論位置，并且近似地將攝影測量交會點F與線段F1F2中點視為同一點，在此基礎(chǔ)上嚴格推導(dǎo)出從攝影測量交會點到水下物點的高程坐標折射改正公式。因此，本文算法的合理性和理論誤差主要取決于攝影測量交會點F偏離線段F1F2中點的程度，即取決于攝影測量交會點F與線段F1F2中點之間的距離δ。

為此，利用WorldView-2立體影像的觀測參數(shù)和其他信息，實際計算攝影測量交會點F到線段F1F2中點的距離δ(符號含義見圖6)，并根據(jù)計算結(jié)果來分析“將攝影測量交會點F和線段F1F2中點視為同一點”的合理性以及由此造成的理論誤差大小。

為了計算同名像點的δ，需要已知空中同名光線的航向視角、側(cè)視角、入視角和水下目標點的水深值。雖然表1給出了整幅圖像航向視角、側(cè)視角、入視角的最大值、最小值和平均值，圖像上每個像素的航向視角、側(cè)視角、入視角等觀測參數(shù)也可以用線性內(nèi)插方法得到。但由于每個像素航向視角、側(cè)視角、入視角的差異都不大，由此引起的δ差異可以忽略不計，為使問題簡單化，此處使用整幅圖像航向視角、側(cè)視角、入視角的平均值作為圖像上每個像素的航向視角、側(cè)視角、入視角，并假設(shè)水下目標點的水深hP分別為5 m、10 m、15 m來計算δ，其結(jié)果(包括中間結(jié)果)如表2所示。

表2 hP=5、10、15 m時的δ、h1和h2

從表2可以看出：在θ、r1、r2一定的情況下，δ隨hP的增大而線性增大，但δ/hP基本不變。例如，當水深hP為5 m、10 m、15 m時，甘泉島的δ分別為0.123 m、0.251 m、0.377 m，δ/hP基本維持在2.5%左右；珊瑚島的δ分別為0.051 m、0.102 m、0.153 m，δ/hP基本維持在1%左右。雖然甘泉島的δ稍大一些，但δ/hP也只有2.5%左右，相對于雙介質(zhì)攝影測量的定位能力而言是不大的，可以忽略不計，因此將攝影測量交會點F與線段F1F2中點近似為同一點是合理的。

3.3 算法的測量精度改善情況試驗

下面通過不同方法計算得到的淺海海底地形數(shù)據(jù)的比較，來驗證本文算法對高程測量精度的改善能力。

首先分別按照本文方法和最有代表性的Murase算法[13]開發(fā)了相應(yīng)的雙介質(zhì)立體攝影測量軟件；然后分別運用每款軟件對甘泉島地區(qū)的WorldView-2多光譜綠色波段立體影像進行雙介質(zhì)立體攝影測量處理，共得到兩組甘泉島周圍淺海海底地形數(shù)據(jù)(如圖10(a)和圖10(b)所示)；接著再分別運用每款軟件對珊瑚島地區(qū)的WorldView-2多光譜綠色波段立體影像進行雙介質(zhì)立體攝影測量處理，共得到兩組珊瑚島周圍淺海海底地形數(shù)據(jù)(圖11(a)和圖11(b))。

本文方法的雙介質(zhì)立體攝影測量處理過程如下：

(1) 運用區(qū)域生長法提取遙感影像上的水域邊界，將水域和陸地區(qū)分開來，以便對水域和陸地區(qū)域?qū)嵤┎煌臄z影測量處理。

(2) 對陸地區(qū)域，按照特征點提取、同名像點匹配、目標點大地坐標立體攝影測量解算的流程，獲得陸地區(qū)域內(nèi)檢測到的所有特征點的大地坐標，然后以離散陸地特征點為基礎(chǔ)內(nèi)插得到陸地區(qū)域規(guī)則格網(wǎng)的大地坐標。

(3) 對水域，按照太陽耀斑消除、特征點提取、同名像點匹配、目標點大地坐標立體攝影測量解算、雙介質(zhì)折射改正的流程，盡可能多地獲得淺海區(qū)域海底特征點的大地坐標，然后以海底離散特征點為基礎(chǔ)內(nèi)插得到淺海區(qū)域規(guī)則格網(wǎng)的大地坐標。

Murase算法[13]的雙介質(zhì)立體攝影測量處理過程與本文方法完全一致，只是在物方坐標折射改正時所使用的公式不同。

圖7 甘泉島WorldView-2多光譜立體影像(MS1的2、3、4波段合成)Fig.7 WorldView-2 multispectral stereo images of Ganquan Island (compositing of MS1 bands2,3 and 4)

為了客觀進行精度比較，需要甘泉島、珊瑚島兩個地區(qū)的高精度淺海海底地形數(shù)據(jù)作為參考基準。但由于數(shù)據(jù)資源和研究經(jīng)費限制，僅購得甘泉島周圍淺海區(qū)域的高精度機載激光測深數(shù)據(jù)，如圖10(c)所示。該數(shù)據(jù)是采用加拿大Optech公司SHOALS-3000[26]測量得到的，測深精度為±0.3 m。由于僅有甘泉島周圍淺海區(qū)域的SHOALS-3000測深數(shù)據(jù)，所以精度比較試驗只能針對甘泉島地區(qū)。

圖8 珊瑚島WorldView-2多光譜立體影像(MS1的2、3、4波段合成)Fig.8 WorldView-2 multispectral stereo images of Shanhu Island (compositing of MS1 bands 2,3 and 4)

圖9 WorldView-2焦平面上的探測器Fig.9 Detectors on WorldView-2 focal plane

圖10 3種方法得到的甘泉島周圍淺海海底地形數(shù)據(jù)Fig.10 Three DEMs of Ganquan Island area derived from different methods

精度試驗方法是：以甘泉島周圍淺海區(qū)域為試驗場，用SHOALS-3000測深數(shù)據(jù)作為基準，將本文方法得到的結(jié)果、Murase方法得到的結(jié)果與SHOALS-3000測深數(shù)據(jù)進行疊置比較，得到各個格網(wǎng)點高程坐標的誤差，最后按式(22)統(tǒng)計計算所有格網(wǎng)點高程值的均方根誤差RMSE，結(jié)果如表3所示。

(22)

表3 甘泉島周圍淺海海底地形數(shù)據(jù)的RMSE

為了直觀反映每個點的雙介質(zhì)攝影測量值(水深值)與機載激光測深值的偏離情況，對每種方法，以機載激光測深值作為平面直角坐標系的橫軸、以雙介質(zhì)攝影測量值作為縱軸，將所有測量點繪制在該平面直角坐標系中，得到每種方法的雙介質(zhì)攝影測量值與機載激光測深值的偏離情況圖，結(jié)果如圖12所示。

圖12 兩種方法得到的甘泉島周圍淺海水深偏離情況圖Fig.12 The relationships between our or Murase’s results and SHOALS-3000’s results

從表3可以看出：對所有格網(wǎng)點，本文方法的均方根誤差為2.35 m，比Murase方法的3.13 m小0.78 m；對水深不低于4 m的格網(wǎng)點(即剔除水質(zhì)較差區(qū)域的淺水點)，本文方法的均方根誤差為1.89 m，仍比Murase方法的2.31 m小0.42 m。由于測量精度的差異主要是物方坐標折射改正算法不同導(dǎo)致的，因此試驗結(jié)果說明了本文方法對提高淺海海底的高程(或水深)測量精度有明顯的效果。

從圖12(a)可以看出：本文方法的結(jié)果中，水深大于4 m的海底點基本分布在吻合線上或附近兩側(cè)(所謂吻合線是指通過坐標原點且斜率為1的直線，該線上所有點的雙介質(zhì)攝影測量值等于機載激光測深值)。從圖12(b)可以看出：Murase方法的結(jié)果中，雖然水深大于4 m的海底測量點也基本分布在吻合線上或附近兩側(cè)，但偏離吻合線的幅度明顯大于本文方法的結(jié)果。圖12(a)與(b)的差異說明：本文方法獲取的淺海海底地形更接近基準數(shù)據(jù)，即說明了本文方法優(yōu)于Murase方法。

從圖12還可以看出：兩種方法在水深小于4 m的淺水區(qū)域都有明顯的系統(tǒng)性誤差，主要原因是淺水區(qū)水質(zhì)較差，水中雜質(zhì)較多，提取到的特征點可能是水中雜質(zhì)點而非海底特征點，從而導(dǎo)致水質(zhì)較差的淺水區(qū)雙介質(zhì)攝影測量精度不高。另外，SHOALS-3000的測深能力為2～50 m，水深小于2 m的機載激光測深數(shù)據(jù)精度也不高。綜合考慮這些因素，針對水深不低于4 m的測量點，統(tǒng)計得到的均方根誤差RMSE更有參考價值。

4 結(jié) 論

提出了一種更嚴密的雙介質(zhì)立體攝影測量物方坐標折射改正算法。該算法用水下目標的空中同名直線光線公垂線段的中點作為攝影測量交會點的理論位置，在此基礎(chǔ)上建立了攝影測量交會點與水下真實物點的相互位置關(guān)系，推導(dǎo)了水下目標點的水深和大地坐標計算公式(即折射改正公式)。理論研究和試驗結(jié)果表明：

(1) 提出的雙介質(zhì)立體攝影測量物方坐標折射改正算法，除了合理地將攝影測量交會點F近似為線段F1F2中點以外，其余公式都是在嚴密的成像幾何關(guān)系基礎(chǔ)上推導(dǎo)出來的，比現(xiàn)有的雙介質(zhì)折射改正算法更嚴密。

(2) 無論空中同名直線光線延長線是否相交，本文的折射改正算法都是適用的。因此，在使用本文算法時，不需要考慮物像相對位置關(guān)系(如空中同名光線所在平面是否重合？空中同名光線關(guān)于L是否對稱？)。

(3) 在空中同名直線光線的入視角r1和r2以及交會角θ一定的情況下，本文折射改正算法的理論誤差δ隨測量點的水深hP增大而線性增大，但δ/hP基本不變。另外，在交會角θ一定的情況下，選擇入視角r1和r2對稱性較好的立體像對進行雙介質(zhì)攝影測量，對降低本文折射改正算法的理論誤差、提高水下目標的測量精度是有利的。

(4) 本文算法能顯著提高水下目標的高程測量精度。以甘泉島周圍淺海海底地形測量為例，在水質(zhì)清澈情況下，能將高程測量精度從Murase方法的2.31 m提高到1.89 m。

應(yīng)當指出，利用航空航天遙感影像和雙介質(zhì)立體攝影測量方法進行水下目標定位，是建立在“空氣和水體均為單一介質(zhì)”、“空氣與水體的分界面為平面”、“水下目標在遙感影像上清晰可見”等條件之上的。該方法僅適用于水體清澈、無海浪的淺海區(qū)域測繪，對影像質(zhì)量和預(yù)處理技術(shù)(如太陽耀斑消除、水域影像增強等)也有更高要求。

[1] CHONG A K, SCHNEIDER K. Two-medium Photogrammetry for Bottlenose Dolphin Studies[J]. Photogrammetric Engineering and Remote Sensing, 2001, 67(5): 621-628.

[2] TELEM G, FILIN S. Photogrammetric Modeling of Underwater Environments[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2010, 65(5): 433-444.

[3] WIMMER M, MANDLBURGER G, RESSL C, et al. Estimation of High Resolution Shallow Water Bathymetry via Two-media-photogrammetry：A Case Study at the Pielach River[C]∥Geophysical Research Abstracts vol. 18, EGU General Assembly 2016. Vienna,Austria: EGU, 2016: 7197.

[4] WESTAWAY R M, LANE S N, HICKS D M. Remote Survey of Large-scale Braided, Gravel-bed Rivers Using Digital Photogrammetry and Image Analysis[J]. International Journal of Remote Sensing, 2003, 24(4): 795-815.

[5] LANE S N, WIDDISON P E, THOMAS R E, et al. Quantification of Braided River Channel Change Using Archival Digital Image Analysis[J]. Earth Surface Processes and Landforms, 2010, 35(8): 971-985.

[6] SHAN Jie. Relative Orientation for Two-media Photogrammetry[J]. The Photogrammetric Record, 1994, 14(84): 993-999.

[7] TELEM G, FILIN S. Photogrammetric Modeling of the Relative Orientation in Underwater Environments[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2013(86): 150-156.

[8] 王有年, 韓玲, 王云. 水下近景攝影測量試驗研究[J]. 測繪學(xué)報, 1988, 17(3): 217-224.

WANG Younian, HAN Ling, WANG Yun. Experimental Research of Underwater Close-range Photogrammetry[J]. Acta Geodaetica et Cartographica Sinica, 1988, 17(3): 217-224.

[9] 常本義. 雙介質(zhì)攝影測量基本公式[J]. 測繪學(xué)報, 1991, 20(4): 288-294.

CHANG Benyi. Basic Formulas of Two-media Photogrammetry[J]. Acta Geodaetica et Cartographica Sinica, 1991, 20(4): 288-294.

[10] 單杰. 雙介質(zhì)攝影測量的相對定向[J]. 測繪學(xué)院學(xué)報, 1993, 10(3): 38-44.

SHAN Jie. Relative Orientation for Two-media Photogrammetry[J]. Journal of the Institute of Surveying and Mapping, 1993, 10(3): 38-44.

[11] 周高偉, 李英成, 任延旭, 等. 低空無人機雙介質(zhì)水下礁盤深度測量試驗與分析[J]. 測繪學(xué)報, 2015, 44(5): 548-554. DOI:10.11947/j.AGCS.2015.20140259.

ZHOU Gaowei, LI Yingcheng, REN Yanxu, et al. Research of Two-media Underwater Reefs Depth Measurement Experiment Based on Low-altitude UAV[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(5): 548-554. DOI:10.11947/j.AGCS.2015.20140259.

[12] 曹彬才, 邱振戈, 朱述龍, 等. 高分辨率衛(wèi)星立體雙介質(zhì)淺水水深測量方法[J]. 測繪學(xué)報, 2016, 45(8): 952-963. DOI: 10.11947/j.AGCS.2016.20150583.

CAO Bincai, QIU Zhenge, ZHU Shulong, et al. Shallow Water Bathymetry through Two-medium Photogrammetry Using High Resolution Satellite Imagery[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(8): 952-963. DOI: 10.11947/j.AGCS.2016.20150583.

[13] MURASE T, TANAKA M, TANI T, et al. A Photogrammetric Correction Procedure for Light Refraction Effects at a Two-medium Boundary[J]. Photogrammetric Engineering & Remote Sensing, 2008, 74(9): 1129-1136.

[14] GEORGOPOULOS A, AGRAFIOTIS P. Documentation of a Submerged Monument Using Improved Two Media Techniques[C]∥Proceedings of the 18th International Conference on Virtual Systems and Multimedia. Milan, Italy: IEEE, 2012: 173-180.

[15] FRYER J G, KNIEST H T. Errors in Depth Determination Caused by Waves in Through-water Photogrammetry[J]. The Photogrammetric Record, 1985, 11(66): 745-753.

[16] WESTAWAY R M, LANE S N, HICKS D M. Remote Sensing of Clear-water, Shallow, Gravel-bed Rivers Using Digital Photogrammetry[J]. Photogrammetric Engineering and Remote Sensing, 2001, 67(11): 1271-1282.

[17] BUTLER J, LANE S, CHANDLER J, et al. Through-water Close Range Digital Photogrammetry in Flume and Field Environments[J]. The Photogrammetric Record, 2002, 17(99): 419-439.

[18] MAAS H G. A Modular Geometric Model for Underwater Photogrammetry[C]∥The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences：Volume XL-5/W5. Dresden, Germany: TSPRS, 2015: 139-141.

[19] DIETRICH J T. Bathymetric Structure-from-motion: Extracting Shallow Stream Bathymetry from Multi-view Stereo Photogrammetry[J]. Earth Surface Processes and Landforms, 2017, 42(2): 355-364.

[20] DigitalGlobe. WorldView-2 Overview[EB/OL]. [2017-03-10]. http:∥content.satimagingcorp.com.s3.amazonaws.com/static/satellite-sensor-specification/WorldView-2-PDF-Download.pdf.

[21] 曹斌, 邱振戈, 曹彬才. 四種遙感淺海水深反演算法的比較[J]. 測繪科學(xué)技術(shù)學(xué)報, 2016, 33(4): 388-393.

CAO Bin, QIU Zhenge, CAO Bincai. Comparison among Four Inverse Algorithms of Water Depth[J]. Journal of Geomatics Science and Technology, 2016, 33(4): 388-393.

[22] 曹斌, 邱振戈, 朱述龍, 等. BP神經(jīng)網(wǎng)絡(luò)遙感水深反演算法的改進[J]. 測繪通報, 2017(2): 40-44. DOI: 10.13474/j.cnki.11-2246.2017.0045.

CAO Bin, QIU Zhen’ge, ZHU Shulong, et al. Improvement of BPANN Based Algorithm for Estimating Water Depth from Satellite Imagery[J]. Bulletin of Surveying and Mapping, 2017(2): 40-44. DOI: 10.13474/j.cnki.11-2246.2017.0045.

[23] HEDLEY J D, HARBORNE A R, MUMBY P J. Technical Note: Simple and Robust Removal of Sun Glint for Mapping Shallow-water Benthos[J]. International Journal of Remote Sensing, 2005, 26(10): 2107-2112.

[24] DigitalGlobe. Radiometric Use of WorldView-2 Imagery[EB/OL]. [2017-03-10]. http:∥www.pancroma.com/downloads/Radiometric_Use_of_WorldView-2_Imagery.pdf.

[25] DigitalGlobe. Imagery Support Data (ISD) Documentation[EB/OL]. [2017-03-10]. http:∥www.docin.com/p-345895199.html.

[26] Optech. SHOALS-3000 Product Brochure[EB/OL]. [2017-03-10]. http:∥pdf.directindustry.com/pdf/optech/shoals-3000-product-brochure/25132-53146.html.

More Rigorous Correction of Refraction Effects in Two-media Stereophoto-grammetry

CAO Bin1，ZHU Shulong2，QIU Zhenge1，CAO Bincai2

1. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China; 2. The Fifth College of the Information Engineering University, Zhengzhou 450001, China

A more rigorous algorithm is presented for correction of refraction effects in two-media stereo photogrammetry. The mid-point of the shortest line segment joining two aerial corresponding rays of a point on an underwater object is used as a photogrammetric intersection point which doesn’t exist when the two rays are non-intersecting. As a result, the uncertainty of the intersection point is removed, the positional relationship between the intersection point and the true object point becomes definite, and the refraction correction formula from the intersection point to the true object point can be strictly derived. The bad effect on the refraction correction is firstly analyzed, which caused by that the two rays are non-intersecting. Then the positional relationship between the intersection point and the true object point is studied. After that, the formulas regarding water depth and geodetic coordinates of points on an underwater object are deduced, that is often known as correction of refraction effects. Finally, the algorithm is tested by two experiments using the data of WorldView-2. The results show that the algorithm is suitable for any case in which whether or not the two aerial corresponding rays of an underwater object point are intersecting, and it can significantly improve the measurement accuracy of underwater object’s elevation.

two-media stereophotogrammetry; correction of refraction effects; photogrammetric intersection point; true object point; underwater object positioning

The Scientific Research Program of Science and Technology Commission of Shanghai Municipality of China (No. 14590502200); The Foundation of the Key Laboratory for Aerial Remote Sensing Technology of National Administration of Surveying, Mapping and Geoinformation (NASG) of China (No. 2016B05)

CAO Bin (1992—), male, postgraduate, majors in marine remote sensing and surveying.

ZHU Shulong

P237

A

1001-1595(2017)09-1182-11

上海市科學(xué)技術(shù)委員會科研計劃(14590502200)；航空遙感技術(shù)國家測繪地理信息局重點實驗室基金(2016B05)

(責(zé)任編輯：叢樹平)

2017-03-13

修回日期： 2017-06-26

曹斌(1992—)，男，碩士生，研究方向為海洋遙感測繪。

E-mail： caobinalonzo@sina.com

朱述龍

E-mail： zhushulong668@sina.com