李慧娜，吉笑笑，韓 杰，關(guān)小果

(1.許昌學(xué)院電氣與機(jī)械工程學(xué)院，461000，河南，許昌；2.許昌學(xué)院城市與環(huán)境學(xué)院，461000，河南，許昌)

0 引言

衛(wèi)星影像幾何定位精度是決定基于衛(wèi)星遙感數(shù)據(jù)開展高精度反演信息提取的重要因素，只有在準(zhǔn)確掌握其幾何定位特性后，才能有針對性的建立相關(guān)校正模型，消除其定位誤差，進(jìn)而提高空間信息提取精度。目前，眾多學(xué)者針對國產(chǎn)高分辨率衛(wèi)星傳感器特點(diǎn)，開展了衛(wèi)星影像幾何定位精度評價(jià)研究。如劉其悅等人建立了基于地基數(shù)據(jù)和參考影像數(shù)據(jù)的衛(wèi)星影像幾何產(chǎn)品真實(shí)性檢驗(yàn)流程[1]。韓杰和謝勇以ZY-3衛(wèi)星傳感器校正產(chǎn)品為參考，對GF-1衛(wèi)星WFV傳感器四相機(jī)影像幾何定位精度進(jìn)行長時(shí)間序列評價(jià)，并開展了相關(guān)的定位誤差補(bǔ)償理論探討[2]。黃世存等人以1:5萬比例尺數(shù)字正射影像為參考，對GF-1衛(wèi)星PMS影像在平原和山區(qū)的幾何定位精度進(jìn)行評價(jià)[3]。樊文鋒等人在對正射糾正后的GF-2衛(wèi)星影像幾何定位精度進(jìn)行評價(jià)后，驗(yàn)證了在少量控制點(diǎn)的輔助下，GF-2衛(wèi)星影像可以滿足全地形1:1萬和山地1:5萬比例尺地形圖的更新要求[4]。徐文等人針對靜止軌道GF-4衛(wèi)星面陣成像特點(diǎn)，在構(gòu)建其幾何成像模型后，評價(jià)了該衛(wèi)星影像數(shù)據(jù)的幾何定位精度[5]。同時(shí)，徐文等人通過對ZY-3衛(wèi)星三線陣TLC影像進(jìn)行幾何精度評價(jià)，證明其在無控的狀態(tài)下，定位精度優(yōu)于6個(gè)像素[6]。

高分六號(GF-6)衛(wèi)星自2018年6月2日于酒泉衛(wèi)星發(fā)射中心成功發(fā)射后，其搭載的寬幅蓋WFV傳感器憑借單相機(jī)超大幅寬成像(優(yōu)于800 km)和新增特征波段(紅邊波段1：0.69～0.73 μm，紅邊波段2：0.73～0.77 μm，紫波段：0.40～0.45 μm和黃波段：0.59～0.63 μm)的特點(diǎn)，已經(jīng)成為環(huán)境保護(hù)、作物估產(chǎn)、農(nóng)業(yè)資源調(diào)查、森林火災(zāi)監(jiān)測等領(lǐng)域研究的重要數(shù)據(jù)源[7-9]。由于WFV傳感器采用單相機(jī)超大幅寬成像，其單景影像大小超過13 GB，這會(huì)對數(shù)據(jù)傳輸及處理帶來不便。因此，WFV單景影像采用分塊存儲(chǔ)方式，即分為左、中、右3塊子影像，其中中間子影像覆蓋范圍及影像大小約為兩側(cè)子影像的2倍[10]。同樣在利用該數(shù)據(jù)進(jìn)行上述研究和應(yīng)用前，需要對該數(shù)據(jù)幾何定位特性進(jìn)行綜合分析。雖然相關(guān)學(xué)者已經(jīng)對GF-1、GF-2、GF-4、ZY-3等國產(chǎn)衛(wèi)星影像進(jìn)行幾何定位精度評價(jià)，但是GF-6衛(wèi)星WFV影像超大幅寬成像及分塊分發(fā)的特點(diǎn)是否對其幾何定位精度有影響，長時(shí)間序列影像幾何定位特性如何，仍需要開展進(jìn)一步的分析與研究。

本文將華北地區(qū)作為實(shí)驗(yàn)區(qū)，以谷歌地球影像為參考，通過比對谷歌地球影像和WFV影像中同名點(diǎn)坐標(biāo)，對WFV單景分塊影像及長時(shí)間序列影像幾何定位精度進(jìn)行評價(jià)。本文研究結(jié)果將為后續(xù)幾何定位誤差源分析和幾何校正模型構(gòu)建提供重要參考信息。

1 數(shù)據(jù)介紹

以華北地區(qū)為研究區(qū)，通過中國資源衛(wèi)星應(yīng)用中心陸地觀測衛(wèi)星數(shù)據(jù)服務(wù)平臺，檢索自衛(wèi)星發(fā)射至2021年4月間云量較少的GF-6衛(wèi)星WFV影像，共檢索得到可用WFV影像47景。因?yàn)榍捌谟嘘P(guān)谷歌地球影像幾何定位精度的研究結(jié)果表明，其影像定位精度能夠滿足WFV影像(空間分辨率16 m)幾何定位精度評價(jià)要求[11-12]。因此，本文將谷歌地球影像作為WFV影像幾何定位精度評價(jià)的參考數(shù)據(jù)。在谷歌地球影像中，選取能夠在WFV影像中可識別的幾何檢查點(diǎn)，共計(jì)79個(gè)，如圖1所示。雖然可用WFV影像集中影像覆蓋范圍有所差異，但是本文在進(jìn)行幾何定位精度檢驗(yàn)時(shí)，保證了每景影像至少均勻覆蓋25個(gè)以上的檢查點(diǎn)。

圖1 檢查點(diǎn)分布位置

2 研究方法

本研究將谷歌地球影像作為幾何參考，將GF-6衛(wèi)星WFV傳感器所獲得的研究區(qū)影像作為研究對象，對其單景影像以及長時(shí)間序列影像的定位精度進(jìn)行分析。本文所采取的研究方法主要包括以下4個(gè)部分：1)利用ENVI軟件分別打開左、中、右3塊子影像，編輯影像頭文件信息，使其關(guān)聯(lián)各自有理多項(xiàng)式參數(shù)(Rational Polynomial Coefficients, RPCs)；2)基于SRTM DEM數(shù)據(jù)和WFV影像RPCs參數(shù)，對WFV影像進(jìn)行正射校正[13-14]，其中SRTM DEM高程精度引起的平面最大誤差約為9 m，相當(dāng)于WFV影像的0.56個(gè)像素，不會(huì)對正射校正精度產(chǎn)生影響[15]；3)將谷歌地球影像和WFV影像中同名點(diǎn)坐標(biāo)轉(zhuǎn)換成通用墨卡托投影坐標(biāo)[16]；4)比對同名點(diǎn)通用墨卡托投影坐標(biāo)差異，統(tǒng)計(jì)誤差平均值和標(biāo)準(zhǔn)差，繪制誤差分布圖，評價(jià)影像幾何定位精度。本文技術(shù)流程圖如圖2所示。

圖2 本文技術(shù)流程圖

3 實(shí)驗(yàn)結(jié)果與分析

3.1 單景影像幾何定位誤差分析

根據(jù)上述幾何定位精度評價(jià)方法，對單景WFV影像幾何定位精度分析。以景號為312708、338701、382955和406694的影像為例，對評價(jià)結(jié)果進(jìn)行展示，其成像時(shí)間分別為2020年3月17日、6月3日、11月10日和2021年1月31日。由于單景影像覆蓋范圍會(huì)涉及到跨帶現(xiàn)象，本文以左側(cè)子影像帶號為基準(zhǔn)，對中間和右側(cè)子影像坐標(biāo)進(jìn)行換算。

從圖3所展示的定位誤差分布圖(箭頭指向表示誤差方向，箭頭長度表示誤差大小，下同)可以看出，單景影像均出現(xiàn)了明顯的系統(tǒng)誤差，其幾何定位誤差方向大致均朝向西側(cè)，即跨軌方向的誤差明顯大于沿軌方向的誤差。從表1可看出，4景WFV影像平均平面定位誤差約為659.87 m，左、中、右3塊子影像平均平面誤差約為244.96 m、570.61 m和1 318.92 m，即定位誤差大小逐漸增加。與此同時(shí)，從圖3和表1還可以看出，沿軌方向的誤差在左側(cè)子影像中呈現(xiàn)負(fù)值，而在中間和右側(cè)子影像中呈現(xiàn)正值，且該值逐漸增加。上述單景影像幾何定位誤差分布特性，將為后續(xù)該影像幾何定位誤差源分析及校正模型構(gòu)建提供重要的參考依據(jù)。

表1 單景影像定位誤差統(tǒng)計(jì)結(jié)果

(a)景號312708 (b)景號338701

3.2 長時(shí)間序列影像幾何定位誤差分析

由于本文對其它成像時(shí)刻的WFV影像進(jìn)行幾何定位精度評價(jià)后，發(fā)現(xiàn)這些影像也具有相似的幾何定位特性。因此，在開展長時(shí)間序列影像幾何定位誤差分析時(shí)，選取每景影像平面平均誤差作為評價(jià)指標(biāo)進(jìn)行分析。

從圖4和表2中可以得出以下結(jié)論：1)無論是左、中、右子影像還是整景影像，跨軌方向的誤差均比沿軌方向的誤差大，且波動(dòng)更為劇烈，即跨軌方向的平均誤差和標(biāo)準(zhǔn)差更大；2)左、中、右側(cè)子影像中沿軌方向的誤差平均值分別為-40.75 m、40.49 m和326.88 m，即沿軌方向的誤差由負(fù)值變?yōu)檎?，且右?cè)子影像的沿軌方向誤差明顯大于其它2塊子影像；3)左、中、右側(cè)子影像平均平面誤差分別為307.89 m、563.81 m和1 188.21 m，即平均平面誤差逐漸增加；4)中間子影像無論在沿軌、跨軌還是平面方向的定位誤差波動(dòng)性最小，左側(cè)子影像次之，右側(cè)子影像波動(dòng)最大，中間子影像的沿軌、跨軌和平面方向定位誤差的標(biāo)準(zhǔn)差分別為29.91 m、80.95 m和83.15 m。

(a)左側(cè)子影像 (b)中間子影像

表2 長時(shí)間序列影像平均定位誤差統(tǒng)計(jì)結(jié)果

4 結(jié)束語

1)本文針對GF-6衛(wèi)星WFV傳感器超大幅寬成像和分塊分發(fā)的特點(diǎn)，以谷歌地球影像為參考，對華北地區(qū)該衛(wèi)星單景影像和長時(shí)間序列影像的幾何定位誤差進(jìn)行分析評價(jià)。研究結(jié)果表明，GF-6衛(wèi)星單景WFV影像呈現(xiàn)出明顯的系統(tǒng)誤差，其誤差方向均朝向西側(cè)，且左、中、右側(cè)子影像平均平面誤差逐漸增加。同時(shí)，長時(shí)間序列影像幾何定位誤差研究結(jié)果還表明對于整景WFV影像而言，具有較好的定位穩(wěn)定性，其平面幾何定位誤差標(biāo)準(zhǔn)差82.74 m(約5個(gè)像素)。

2)由于本文采用的是人工刺點(diǎn)的方式完成同名點(diǎn)信息提取，這可能會(huì)對評價(jià)精度產(chǎn)生一定的影響，后續(xù)將考慮選取與WFV影像空間分辨率相近的Landsat8衛(wèi)星OLI傳感器全色影像(空間分辨率15 m)為參考，采用高精度影像匹配方法實(shí)現(xiàn)同名點(diǎn)信息提取，降低人為因素的影響。

3)由于本文僅對WFV影像幾何定位誤差特性進(jìn)行了系統(tǒng)分析，并未進(jìn)一步探討是何種因素導(dǎo)致了該誤差特性。因此，后續(xù)將在分析誤差源的基礎(chǔ)上，構(gòu)建適用于WFV影像幾何精校正模型。

致謝：感謝中國資源衛(wèi)星應(yīng)用中心為本文研究提供GF-6衛(wèi)星WFV影像數(shù)據(jù)，谷歌地球提供參考影像！