利用ADS80影像測制西部地區(qū)大比例尺地形圖的快速流程

陳建斌，王明孝，張 萍，蔡曉麗

(1.68029部隊，甘肅蘭州730020;2.甘肅林學(xué)院，甘肅天水741000)

利用ADS80影像測制西部地區(qū)大比例尺地形圖的快速流程

陳建斌1，王明孝1，張 萍1，蔡曉麗2

(1.68029部隊，甘肅蘭州730020;2.甘肅林學(xué)院，甘肅天水741000)

利用ADS80影像對西部某區(qū)域進(jìn)行無控空中三角測量的精度作對比。研究考慮地形高差影響的高程異常擬合方法，提出利用ADS80影像快速測制西部大比例尺地形圖的作業(yè)流程，并將該流程成功應(yīng)用于某風(fēng)電場基礎(chǔ)測繪項目中。

ADS80;西部地區(qū);快速流程;高程擬合;數(shù)據(jù)重構(gòu)

一、引 言

西部地區(qū)以高原高寒和沙漠戈壁為主，地形起伏大，航攝影像投影誤差大，外業(yè)工作開展困難。航測成圖如果采用傳統(tǒng)模擬航片或面陣數(shù)字影像，由于其連接強(qiáng)度弱，接邊工作量大，地面依賴性高，勢必造成作業(yè)速度慢、周期長、耗資巨大的結(jié)果。

ADS80是一種改進(jìn)的多線陣推掃式成像傳感器。它集成了POS系統(tǒng)，能直接獲取影像的外方位元素;它拍攝的長條帶影像提高了后期處理的效率;大基高比能滿足大比例尺測圖要求;100%三度重疊影像能提供多種立體模型。

在某風(fēng)電場1300 km2基礎(chǔ)測繪項目中，采用ADS80相機(jī)，共規(guī)劃拍攝33個條帶，條帶長度從17 257 m到43561 m不等，像元大小6.5 μm，地面分辨率GSD為18 cm。

二、GNSS/IMU輔助空中三角測量精度分析

由于ADS采取的是連續(xù)推掃的成像方式，因此飛機(jī)的飛行姿態(tài)變化將直接影響成像的效果。雖然ADS配置的PAV陀螺穩(wěn)定平臺具有自動調(diào)整姿態(tài)的功能，在進(jìn)行中小比例尺航空攝影時，由于飛機(jī)高度相對較高，氣流比較平穩(wěn)，而且使用的運8、安30等較大型航測飛機(jī)受氣流影響也相對較小，俯仰和傾斜變化通常不超過5°。但在大比例尺航空攝影時，一般使用運5、運12等飛機(jī)，該類飛機(jī)自重小，受氣流的影響很大，而低空氣流變化本身就較劇烈，俯仰和傾斜變化經(jīng)常超過5°。另一方面，PAV陀螺穩(wěn)定平臺對姿態(tài)變化的修正需要一定的響應(yīng)時間，對瞬間大姿態(tài)變化的修正往往反應(yīng)滯后，因此，空中飛機(jī)姿態(tài)的變化情況也能在連續(xù)推掃記錄的影像中完全體現(xiàn)出來，并直接影響影像的后續(xù)處理質(zhì)量。

GNSS/IMU觀測數(shù)據(jù)的引入，大量減少了外業(yè)控制測量的工作量。一直以來，許多學(xué)者都對GNSS/IMU輔助空三所需的控制點數(shù)量進(jìn)行了研究［1］。在該項目中，筆者分別用無控、位于區(qū)域中心1個控制點、位于區(qū)域四周4個控制點、分布均勻的9個控制點等4種方案對測區(qū)進(jìn)行了對比試驗。空三解算使用ORIMA軟件。結(jié)果如表1所示，其中，精度用檢查點中誤差表示。

表1 4種加密方案精度比較

試驗結(jié)論如下:

1)在沒有地面控制的情況下，將POS(IPAS)數(shù)據(jù)作為帶權(quán)觀測值引入到區(qū)域網(wǎng)平差系統(tǒng)中，平差結(jié)果存在一定的系統(tǒng)誤差。

2)當(dāng)引入1個地面控制點時，精度得到了較大提高，平面達(dá)到1個像素，高程達(dá)到1.3個像素，高程精度提高明顯。這表明POS(IPAS)數(shù)據(jù)存在一定的系統(tǒng)誤差，高程分量尤為明顯。

3)當(dāng)引入4個地面控制點時，精度得到了進(jìn)一步提高，平面、高程精度都達(dá)到1個像素，但提高的幅度已不明顯。

4)當(dāng)引入分布均勻的9個控制點時，精度幾乎沒有提高。

試驗表明，POS(IPAS)數(shù)據(jù)的引入，對地面控制點的要求降到了最低，只需1個位于區(qū)域中心的控制點，就能很好地抵償POS數(shù)據(jù)的系統(tǒng)誤差。如果要解決坐標(biāo)系統(tǒng)的轉(zhuǎn)換問題，最好使用區(qū)域四周4個控制點的方案。

航測內(nèi)業(yè)規(guī)范對空三加密的精度的規(guī)定如表2所示［2］。

表2 空三加密內(nèi)業(yè)規(guī)范 m

限差放寬后，無控加密的最大誤差為0.40 m，小于規(guī)范中平地的高程中誤差限差±0.42 m，符合規(guī)范的要求。因此無地面控制的GNSS/IMU輔助空中三角測量的精度能滿足西部困難地區(qū)1∶2000成圖。

三、快速作業(yè)流程

1.ADS80通用作業(yè)流程

ADS80的推廣應(yīng)用，使得航測成圖效率大大提高。這得益于ADS80本身的特點:

1)集成POS系統(tǒng)，可以直接獲取影像外方位元素，減少對地面控制的依賴。

2)長條帶影像，提高了影像處理效率，航帶內(nèi)無需接邊，強(qiáng)度高。

3)100%三度重疊影像具有影像冗余，能以多種組合提供立體模型。

如圖1所示，對同一個測量任務(wù)，如果采用框幅式像片或面陣影像，外業(yè)控制測量、空三加密工作量將非常大。

當(dāng)采用ADS80通用流程時，空三加密工序被外業(yè)控制分割為前后兩部分:第1部分包括POS數(shù)據(jù)解算、連接點選取、L1級影像糾正、控制點概略選點、控制片制作;第2部分包括控制點轉(zhuǎn)刺量測、空三解算。兩部分之間是外業(yè)控制測量工序。這種流程不適用于應(yīng)急測繪項目。因為空三加密、外業(yè)控制測量耗時太長。進(jìn)行外業(yè)控制測量時，內(nèi)業(yè)人員只能等待，在西部地區(qū)，外業(yè)展開困難，影響外業(yè)進(jìn)度的因素太多，經(jīng)常會延誤。這種流程難以滿足緊急任務(wù)(如應(yīng)急保障、搶險救災(zāi))的需求。

2.ADS80快速作業(yè)流程

在總結(jié)了上述兩種流程的基礎(chǔ)上，本項目提出ADS80快速作業(yè)流程，如圖1所示。

快速流程的特點如下:

1)快速流程采用無控空三加密(無控空三加密的精度下文有具體分析)，因為不再依賴于地面控制，空三加密還可以分區(qū)域進(jìn)行，子區(qū)域空三加密成果立刻可以提交下工序，這樣減少了下工序的等待時間。

2)外業(yè)控制測量工序可以與4D數(shù)據(jù)采編工序同時進(jìn)行，避免了內(nèi)業(yè)等待造成的窩工。

3)外業(yè)控制測量不再是提供控制點的必須工序，而是為矢量重構(gòu)提供轉(zhuǎn)換參數(shù)和高程異常的可選工序，如果有測區(qū)精確的轉(zhuǎn)換參數(shù)和高程異常，則可以省去外業(yè)控制測量工作。

4)流程的最后是對矢量數(shù)據(jù)進(jìn)行重構(gòu)，包括坐標(biāo)系統(tǒng)轉(zhuǎn)換、高程轉(zhuǎn)換。

圖1 作業(yè)流程對比

四、快速獲取更精確的水準(zhǔn)面精化文件

ADS80快速作業(yè)流程的最后是對矢量數(shù)據(jù)進(jìn)行重構(gòu)，主要解決高程基準(zhǔn)問題，高程重構(gòu)的精度依賴于水準(zhǔn)面精化文件或高程異常模型。

我國似大地水準(zhǔn)面分為3個等級［3］，如表3所示。

表3 我國似大地水準(zhǔn)面等級

在西部廣大地區(qū)，缺少省級似大地水準(zhǔn)面成果，國家級成果精度只有±0.3 m(平地)，無法滿足1∶2000地形圖需要。因此，必須要用合適的擬合方法獲得更精確的似大地水準(zhǔn)面。

航測成圖時，一般不考慮地形高差變化對高程異常的影響。在西部地區(qū)，當(dāng)海拔為2000 m時，由地形因素引起的高程異常變化量最大可以達(dá)到0.2 m，已經(jīng)接近空中三角測量檢查點的中誤差限差，因此進(jìn)行高程異常擬合時，必須考慮地面高差引起的高程異常變化量。

精確地確定大地水準(zhǔn)面需要地球重力場模型及DEM高程模型，其公式為

對于小區(qū)域而言，前兩項NGM、NΔg變化相對穩(wěn)定，只需考慮Nh。如果將高程異常量作為研究對象，有如下公式

式中，ζ0為高程異常的長波項;ζTC為短波項。式(1)中的Nh，是由地形起伏引起的。

地形改正可以看做是重力場的噪聲。如果能求解ζTC，并在ζ中將ζTC扣除，則可近似認(rèn)為ζ0=ζ－ζTC為一個光滑的幾何曲面，這樣就可以利用測區(qū)中已知GPS水準(zhǔn)點的正常高和大地高確定一個多項式曲面函數(shù)，即小區(qū)域高程異常模型，再利用該模型確定其他點的ζ0值。然后，利用ζ=ζ0+ζTC，可求得ζ。

式中，Tc為地形起伏對地面擾動位的影響;r為積分元到積分點p的距離;hi為參考高程;G為萬有引力常數(shù);ρ為地球質(zhì)量密度;hi－h(huán)r是正常高差(小區(qū)域內(nèi)可近似用大地高差Hi－Hr代替)。

詳細(xì)計算步驟為:

1)引入全球格網(wǎng)間距為90 m的SRTM(shuttle radar topography mission)，SRTM數(shù)據(jù)主要是由美國太空總署NASA和國防部國家測繪局NIMA聯(lián)合測量的全球DEM，用已知的地面控制點對區(qū)域SRTM進(jìn)行多項式糾正，得到最佳附合于已知點的區(qū)域DEM。再按式(3)求解DEM格網(wǎng)點的ζTC(格網(wǎng)點)，結(jié)果如圖2所示。

2)利用多項式擬合的方法(此時，用二階多項式即可)，確定一個二次曲面，將DEM格網(wǎng)點帶入此二次曲面函數(shù)，求得每一格網(wǎng)點的ζ0(格網(wǎng)點)，結(jié)果如圖3所示。

3)利用式(2)，求得每一格網(wǎng)點的ζ(格網(wǎng)點)，結(jié)果如圖4所示。

圖2 ζTC(格網(wǎng)點)

圖3 ζ0(格網(wǎng)點)

圖4 考慮地形高差因素生成的高程異常模型

五、矢量數(shù)據(jù)重構(gòu)

矢量數(shù)據(jù)重構(gòu)解決兩個問題:平面坐標(biāo)轉(zhuǎn)換、高程基準(zhǔn)轉(zhuǎn)換。

1.平面坐標(biāo)轉(zhuǎn)換

平面坐標(biāo)轉(zhuǎn)換通過七參數(shù)模型來實現(xiàn)，七參數(shù)模型如式(4)、圖5所示。涉及3個線元素、3個角元素、1個尺度元素［4］。

圖5 七參數(shù)轉(zhuǎn)換示意圖

2.高程基準(zhǔn)轉(zhuǎn)換

矢量重構(gòu)之高程基準(zhǔn)轉(zhuǎn)換較復(fù)雜。因為測區(qū)內(nèi)的高程異常值不是常數(shù)，不能用統(tǒng)一增加常數(shù)的方法解決。要先按本文前述的方法擬合出一個高程異常格網(wǎng)模型，再與矢量數(shù)據(jù)套合疊加，生成新的矢量數(shù)據(jù)。這種方法適用于對水系、植被、高程點等要素。

因為等高線有其特殊性:一根等高線上的高程值唯一。而用這種方法重構(gòu)等高線，生成的曲線上的高程值不再唯一，得到的將是“非等高線”。

解決方法:先將等高線構(gòu)建TIN、生成DEM，對DEM進(jìn)行重構(gòu)建，最后再反生成等高線，如圖6所示。

圖6 等高線重構(gòu)流程

在構(gòu)建TIN時，如果遇到長舌狀的等高線，很容易會產(chǎn)生平三角，即單根等高線自構(gòu)。

GeoTIN軟件采用自動增加特征點的方法來解決此問題。如圖7所示，在平三角處自動增加特征點，進(jìn)行二次構(gòu)網(wǎng)。但仍不可避免會有自構(gòu)問題。

圖7 自動增加特征點并重新構(gòu)網(wǎng)

本項目采用了另外的構(gòu)網(wǎng)思路，如圖8所示，通過設(shè)定適當(dāng)?shù)乃阉靼霃剑沟雀呔€在一定范圍內(nèi)只能與相鄰等高線構(gòu)網(wǎng)，避免“自構(gòu)”，減少了等高線重構(gòu)產(chǎn)生的錯誤。

圖8 強(qiáng)制與相鄰等高線構(gòu)網(wǎng)

圖9為在WGS-84坐標(biāo)系、大地高基準(zhǔn)上采集的矢量數(shù)據(jù)。圖10為重構(gòu)后的矢量數(shù)據(jù)。

圖9 基于WGS-84坐標(biāo)系、大地高的矢量數(shù)據(jù)

圖10 重構(gòu)后的矢量數(shù)據(jù)

六、應(yīng) 用

該技術(shù)在某風(fēng)電場基礎(chǔ)測繪任務(wù)中得到了實際應(yīng)用。測區(qū)覆蓋面積1300 km2，海拔高度在830～1420 m之間，地貌上表現(xiàn)為戈壁平原，以山前沖積為主，地勢開闊，植被稀少。利用該技術(shù)在2個月時間內(nèi)完成了風(fēng)電場航攝、外業(yè)控制、1∶2000地形圖內(nèi)業(yè)成圖任務(wù)，極大地節(jié)約了人力物力。

［1］ 劉軍.GPS/IMU輔助機(jī)載線陣CCD影像定位技術(shù)研究［D］.鄭州:信息工程大學(xué)，2007:66-82.

［2］ 國家測繪局.GB/T 23236—2009數(shù)字航空攝影空中三角測量規(guī)范［S］.北京:中國標(biāo)準(zhǔn)出版社，2009:2-3.

［3］ 國家測繪局.GB/T 23709—2009區(qū)域似大地水準(zhǔn)面精化基本技術(shù)規(guī)定［S］.北京:中國標(biāo)準(zhǔn)出版社，2009: 2-3.

［4］ 熊介.橢球大地測量學(xué)［M］.北京:解放軍出版社，1985:311-312.

A Quick Working Flow of Large Scale Mapping in Western Region Based on ADS80 Images

CHEN Jianbin，WANG Mingxiao，ZHANG Ping，CAI Xiaoli

0494-0911(2012)08-0043-04

P23

B

2012-01-13

陳建斌(1976—)，男，甘肅隴西人，碩士，工程師，主要研究方向為攝影測量與遙感。