機載三線陣ccd傳感器幾何定位處理

機載三線陣ccd傳感器幾何定位處理

sd數(shù)字航空攝影系統(tǒng)的開發(fā)是當前攝影測量的重要趨勢之一?；?d矩陣概念的數(shù)字航空攝影機在行業(yè)中引起了越來越多的關注。但從實際應用來看,直接獲取的外方位元素精度還難以滿足高精度定位的要求,因此,有必要在少量地面控制的支持下,通過區(qū)域網(wǎng)平差來獲取加密點坐標和更高精度的定向數(shù)據(jù)。光束法區(qū)域網(wǎng)平差是理論上最直接、最嚴密、精度最高的空中三角測量區(qū)域網(wǎng)平差方法,但在解析空中三角測量的生產和科學實踐中,其優(yōu)勢卻往往難以體現(xiàn)?；诟郊訁?shù)的自檢校光束法區(qū)域網(wǎng)平差是最有效的系統(tǒng)誤差補償方法,將自檢校光束法區(qū)域網(wǎng)平差技術應用于傳統(tǒng)膠片式模擬相機影像,攝影測量學者進行過大量研究和實驗,但目前針對新型機載三線陣CCD傳感器影像的研究相對較少。本文將基于附加參數(shù)的自檢校技術應用于機載三線陣CCD影像的精確定位處理,建立了自檢校光束法區(qū)域網(wǎng)平差模型,并利用登封地區(qū)ADS40全套航攝飛行數(shù)據(jù),對不同平差模型、不同控制點數(shù)量對影像定位精度的影響進行了評估。1網(wǎng)絡沙面模型的自適應模型1.1帶權觀測值的修正自檢校光束法平差的基本公式為:x+Δx=-fa1(X-Xs)+b1(Y-Ys)+c1(Ζ-Ζs)a3(X-Xs)+b3(Y-Ys)+c3(Ζ-Ζs)y+Δy=-fa2(X-Xs)+b2(Y-Ys)+c2(Ζ-Ζs)a3(X-Xs)+b3(Y-Ys)+c3(Ζ-Ζs)(1)x+Δx=?fa1(X?Xs)+b1(Y?Ys)+c1(Z?Zs)a3(X?Xs)+b3(Y?Ys)+c3(Z?Zs)y+Δy=?fa2(X?Xs)+b2(Y?Ys)+c2(Z?Zs)a3(X?Xs)+b3(Y?Ys)+c3(Z?Zs)(1)式中,Δx、Δy表示該像點處引入的附加參數(shù)函數(shù)。如果將控制點也處理成帶權觀測值的話,則平差的基本誤差方程式為:Vx=AX+BXg+CXs-Lx?ΡxVg=EgXg-Lg?ΡgVs=EsXs-Ls?Ρs(2)式中,X為外方位元素和坐標未知數(shù)改正數(shù)矢量;A為相應的誤差方程式系數(shù)矩陣;Lx為像點坐標的觀測值矢量;Px為像點坐標的權矩陣;Xg為控制點坐標的改正數(shù)矢量;B為相應的誤差方程式系數(shù)矩陣;Lg為控制點坐標改正數(shù)的觀測值矢量;Pg為控制點坐標的權矩陣;Xs為自檢校附加參數(shù)矢量;C為相應的誤差方程式系數(shù)矩陣;Ls為附加參數(shù)的觀測值矢量;Pc為附加參數(shù)的權矩陣。1.2參數(shù)模型在自適應學校1.2.1偏心畸變公式1)鏡頭光學畸變差是指相機物鏡系統(tǒng)設計、制作和裝配引起的像點偏離其理想位置的點位誤差。鏡頭的光學畸變是非線性的,主要包括徑向畸變和偏心畸變。徑向畸變可表示為:Δx=k1r3+k2r5+k3r7+?(3)分別將其分解到像平面坐標系的x軸和y軸上:{Δxr=k1ˉxr2+k2ˉxr4+k3ˉxr6+?Δyr=k1ˉyr2+k2ˉyr4+k3ˉyr6+?(4)其中,ˉx=(x-xp);ˉy=(x-yp);r2=ˉx2+ˉy2;(xp,yp)為像主點坐標。偏心畸變主要是由光學系統(tǒng)光心與幾何中心不一致造成的,即鏡頭器件的光學中心不能嚴格共線,它使像點既產生徑向偏差又產生切向偏差,表達式為:p(r)=√p21+p22?r2(5)將其分解到像平面坐標系的x軸和y軸上,則有:{Δxd=p1(r2+2ˉx2)+2p2ˉxˉyΔyd=2p1ˉxˉy+p2(r2+2ˉy2)(6)式中,p1、p2為偏心畸變系數(shù)。偏心畸變在數(shù)量上要比徑向畸變小得多。2)傳感器鏡頭主點坐標(xp,yp)的偏移,該偏移量在x、y方向上的偏移量用常量Δxp、Δyp表示,同時,設傳感器焦距的變化量為Δf,該變化量在像點上的影響可以模擬為:dxf=-Δff(x-xp)dyf=-Δff(y-yp)(7)則某像點p,鏡頭畸變造成的像點誤差模型為:{Δx′=Δxp-Δffˉx+(k1r2+k2r4+k3r6)ˉx+p1(r2+2ˉx2)+2p2ˉxˉyΔy′=Δyp-Δffˉy+(k1r2+k2r4+k3r6)ˉy+2p1ˉxˉy+p2(r2+2ˉy2)(8)1.2.2y軸方向的移動影響1)像元尺寸的變化主要影響成像比例尺。如圖1所示,建立線陣CCD掃描線坐標系,定義x軸為飛行方向,y軸為CCD掃描方向,Np為線陣CCD像元數(shù)。設單個像元原尺寸為(px,py),變化率為(dxp,dyp),沿x和y軸方向總變化為(dpx,dpy),由于沿y軸方向有Np個像元,x軸方向僅有1個像元,顯然有:{dpx=dxpdpy=Νp?dyp(9)可以只考慮沿y軸方向的變形影響,忽略x軸方向。2)CCD在焦平面內的平移將使得像主點偏離原定位置,沿x和y軸方向的移動影響可用常量dxc和dyc表示,如圖2所示。3)設線性CCD在焦平面內旋轉了角度θ,dxθ和dyθ分別為由于CCD旋轉而造成的在飛行方向和掃描方向的像坐標誤差,如圖3所示,則有:{dxθ=ˉysinθdyθ=ˉy-ˉycosθ=ˉy(1-cosθ)(10)一般情況下dyθ很小,常不予考慮,只進行x軸方向的改正。因此,對單線陣CCD,在考慮沿y軸方向的比例尺變化影響、線陣CCD在焦平面內的平移影響、沿x軸方向的旋轉影響后像點誤差模型為:{Δx′=(Δxp+dxc)-Δffˉx+(k1?r2+k2?r4+k3?r6)ˉx+p1(r2+2ˉx2)+2p2ˉxˉy+ˉysinθΔy′=(Δyp+dyc)-Δffˉy+(k1?r2+k2?r4+k3?r6)ˉy+2p1ˉxˉy+p2(r2+2ˉy2)+ˉysy(11)其中,dxc和dyc可合并到主點坐標偏移量Δxp、Δyp中;θ為線性CCD在焦平面內旋轉的角度;sy為比例因子。ADS40采用的是單鏡頭多線陣CCD傳感器,各線陣CCD安置在同一焦平面上,共用一套光學系統(tǒng),因此可采用同一組光學畸變系數(shù)。對線陣CCDj(j=1,2,3),將dxcj和dycj合并到主點坐標偏移量Δxpj、Δypj中,旋轉角度為θj,比例因子為syj,則有:{Δx′j=Δxpj-Δffˉx+(k1?r2+k2?r4+k3?r6)ˉx+p1(r2+2ˉx2)+2p2ˉxˉy+ˉysinθjΔy′j=Δypj-Δffˉy+(k1?r2+k2?r4+k3?r6)ˉy+2p1ˉxˉy+p2(r2+2ˉy2)+ˉysyj(12)1.3外方位元素變化模型建立如何在三線陣影像的自檢校光束法平差中引入POS數(shù)據(jù)是關鍵問題之一。三線陣傳感器定向參數(shù)眾多,且每個采樣周期都對應一套外方位元素,解算過程中不可能、也沒有必要一一求解,因此,選擇合適的描述外方位元素變化的數(shù)學模型就十分關鍵。本文采用3種模型進行研究,分別為低階多項式模型(LPM)、分段多項式模型(PPM)和定向片內插模型(OIM)。1測量誤差類型化考慮POS系統(tǒng)直接測定的外方位元素仍存在一定測量誤差,且呈現(xiàn)較強的系統(tǒng)性,因此,采用LPM來描述POS數(shù)據(jù)的測量誤差,其實質是對POS數(shù)據(jù)進行漂移誤差改正,此時,瞬時時刻t的外方位元素可表示為:{XS=XGΡS+aX+bX?(t-t0)YS=YGΡS+aY+bY?(t-t0)ΖS=ΖGΡS+aΖ+bΖ?(t-t0)ω=ωΙΜU+aω+bω?(t-t0)φ=φΙΜU+aφ+bφ?(t-t0)κ=κΙΜU+aκ+bκ?(t-t0)(13)2外方位元素約束條件將飛行軌道分成若干段,每一段采用一個時間t的LPM來描述外方位元素,并在軌道分段處考慮外方位元素變化的連續(xù)和光滑。對于第i個軌道分段內的時刻t,PPM模型可表示為:{XS=XGΡS+Xi0+Xi1?ˉt+Xi2?ˉt2YS=YGΡS+Yi0+Yi1?ˉt+Yi2?ˉt2ΖS=ΖGΡS+Ζi0+Ζi1?ˉt+Ζi2?ˉt2ω=ωΙΜU+ωi0+ωi1?ˉt+ωi2?ˉt2φ=φΙΜU+φi0+φi1?ˉt+φi2?ˉt2κ=κΙΜU+κi0+κi1?ˉt+κi2?ˉt2(14)其中,ˉt=t-t0,在分段邊界處,由相鄰分段多項式i和i+1計算出的外方位元素應滿足相等的約束條件。此外,如果考慮到軌道光滑,也可附加一階導數(shù)相等的條件。3lagrange廣義內插在飛行軌道上按照一定時間間隔抽取若干離散的掃描周期,稱為定向片,并以定向片時刻的外方位元素作為平差未知數(shù),其他采樣時刻的外方位元素利用定向片時刻的外方位元素通過Lagrange多項式內插得到。設下視像點ΡN成像于掃描行j,其位于定向片Κ和Κ+1之間,當采用三次Lagrange多項式內插時,則第j掃描行的外方位元素(XjS,YjS,ΖjS,ωj,φj,κj)需利用相鄰的4個定向片內插得出:Ρ(tj)=Κ+2∑i=Κ-1(Ρ(ti)?Κ+2∏k=Κ-1k≠it-tkti-tk)(15)其中,Ρ(t)表示t時刻的某一外方位元素分量。2實驗與分析2.1兩組平差模型的建立實驗采用2009年8月獲取的河南登封地區(qū)ADS40影像數(shù)據(jù),選取航高為600m和1000m的數(shù)據(jù)進行測試,分別記為數(shù)據(jù)A和數(shù)據(jù)B。數(shù)據(jù)A相對于地面的平均航高約600m,相應地面分辨率為6cm(GSD),選擇了其中的4條東西航線和2條南北構架航線,影像覆蓋區(qū)內共選取控制點65個。數(shù)據(jù)B相對于平均地面的航高約1000m,平均地面分辨率為10cm,包括12條東西向平行航線,測區(qū)兩端加飛4條南北向構架航線,選擇其中的4條東西航線和2條構架航線,影像覆蓋區(qū)內共選取控制點43個。采用自行研制的ADS40影像平差軟件,對每一組ADS40數(shù)據(jù)分別進行常規(guī)光束法平差和自檢校光束法區(qū)域網(wǎng)平差,平差模型分別采用LPM、PPM和OIM三種實驗結果見表1、2。為驗證控制點數(shù)量對平差精度的影響,實驗中分別設置了5種不同數(shù)量的控制點配置方案。2.2控制平差精度分析通過實驗可以看出,本文設計的自檢校光束法區(qū)域網(wǎng)平差模型可以有效消除或削弱系統(tǒng)誤差的影響,顯著提升定位精度。①在無地面控制的條件下,常規(guī)區(qū)域網(wǎng)平差精度較差,如果以GSD計,采用LPM模型,數(shù)據(jù)A的X、Y、Z方向精度分別為3.9、3.6和4.6個GSD,數(shù)據(jù)B則分別為2、2.7和2.9個GSD,采用PPM和OIM模型情況稍好,但總體看定位精度不夠理想;采用自檢校平差后則能有效提高定位精度,數(shù)據(jù)A提高到3.0、2.0和2.2個GSD,數(shù)據(jù)B則分別提高到1.8、2.1和1.4個GSD。②采用不同控制點數(shù)量的實驗表明,無論是常規(guī)區(qū)域網(wǎng)平差,還是自檢校平差,增加控制點數(shù)量均有助于改善定位精度,但程度有限,當控制點數(shù)量到一定程度后,平差精度趨于平穩(wěn),再增加控制點意義不大。③采用自檢校光束法平差后,同種控制條件下定位精度有明顯改善,以PPM模型為例,在8個控制點條件下數(shù)據(jù)A的X、Y、Z方向精度分別由0.059m、0.172m和0.151m,提高到0.043m、0.046m和0.092m,數(shù)據(jù)B的精度分別由0.106m、0.179m和0.209m,提高到0.087m、0.126m和0.205m。④在3種平差模型中,LPM模型最為簡單方便,解算速度最快,但定位精度較差;PPM模型由于在每個軌道分段內引入一組漂移改正參數(shù),相比LPM模型精度有所提升,但幅度不大,且由于引入過多模型參數(shù),導致解算不夠穩(wěn)定;OIM模型解算最為穩(wěn)定,精度最優(yōu),對三線陣影像區(qū)域網(wǎng)平差有良好的適用性,但解算速度稍慢,有待進一步改進。圖4顯示的是數(shù)據(jù)A和B采用OIM模型進行常規(guī)區(qū)域網(wǎng)平差和自檢校平差后的精度對比。其中,OIMX、OIMY、OIMZ分別表示常規(guī)區(qū)域網(wǎng)平差后X、Y、Z三個方向的均方根誤差;SOIMX、SOIMY、SOIMZ分別表示自檢校區(qū)域網(wǎng)平差后X、Y、Z三個方向的均方根誤差?？梢钥闯?隨著控制點數(shù)量的增加,兩種方式的平差在精度上均逐步提升,當控制點多于8個以后,精度變化趨于平穩(wěn);在相同控制條件下,自檢校平差的平面及高程精度均優(yōu)于常規(guī)區(qū)域網(wǎng)平差。其中,平面Y方向精度改善最為明顯,在采用12個控制點解算趨于穩(wěn)定時,數(shù)據(jù)A由0.165m提高到0.034m,數(shù)據(jù)B由0.119m提高到0.045m,提升幅度為62%;Z方向次之;平面X方向精度改善幅度不大。進一步分析發(fā)現(xiàn),在常規(guī)區(qū)域網(wǎng)平差情況下,Y方向定位精度要明顯低于X方向,這表明沿線陣CCD掃描方向可能存在較大的系統(tǒng)性誤差,采用自檢校平差后對Y方