魏丹丹，甘甫平，尚　坤，梁樹能，何　嬌

(1.中國國土資源航空物探遙感中心，北京 100083；2.中國地質(zhì)大學(xué)(北京) 地球科學(xué)與資源學(xué)院，北京 100083)

0　引言

近年來由于遙感技術(shù)的不斷成熟和發(fā)展，通過遙感圖像可以獲得越來越多的信息[1]。例如，通過遙感圖像，可以測制地形圖[2]，為國土資源、農(nóng)業(yè)和林業(yè)等領(lǐng)域提供服務(wù)[3]。但由于各種因素，如平臺位置、運(yùn)動(dòng)狀態(tài)變化和地球的曲率等因素造成的誤差影響圖像的幾何定位精度[4]，造成圖像幾何質(zhì)量下降[5]。因此，對遙感圖像幾何定位精度的評價(jià)成為遙感領(lǐng)域的重要研究課題之一[6]。

地面接收單位接收遙感數(shù)據(jù)后首先對遙感平臺、地球傳感器的各種參數(shù)進(jìn)行處理和校正，但仍不能滿足用戶的需求，對于處理過的圖像還要進(jìn)行進(jìn)一步的幾何校正和輻射校正；即使對接收到的遙感圖像進(jìn)行了一系列校正，但在定位和定量精度方面仍存在一定的誤差，這就需要對該誤差進(jìn)行衡量、計(jì)算，這就使得定位精度成為衡量遙感圖像應(yīng)用價(jià)值的重要指標(biāo)之一。因此，為能夠給用戶提供更優(yōu)質(zhì)的遙感圖像，必須研究圖像的定位精度[7]和變形情況[8]。本文從用戶應(yīng)用角度出發(fā)結(jié)合在軌測試工作，對圖像幾何定位精度開展評價(jià)。

圖像外部幾何定位精度是指圖像在地理參考坐標(biāo)系中的絕對位置精度[9]，也指經(jīng)過系統(tǒng)幾何校正后的圖像上的地理位置和真實(shí)地理位置之間的差異[10]。圖像的內(nèi)部幾何變形包括長度、角度及放射變形等，圖像內(nèi)部幾何定位精度的評價(jià)包括圖像變形的絕對量評價(jià)[11]和整幅圖像變形的一致性評價(jià)[12]。

CEBRS-04星是中巴資源04號陸地觀測衛(wèi)星，其上搭載有5 m全色(簡稱P5)、10 m空間分辨率多光譜相機(jī)(P10)、20 m多光譜相機(jī)和紅外多光譜相機(jī)，寬視場成像儀5個(gè)載荷，本文主要針對5 m全色相機(jī)和10 m多光譜相機(jī)。

1　圖像定位精度評價(jià)

根據(jù)行業(yè)應(yīng)用要求，本文主要以遙感影像平面制作規(guī)范(GB/T 15968—2008)、數(shù)字地形圖系列和基本要求(GB/T 18315—2001)和基礎(chǔ)地理信息數(shù)字產(chǎn)品1∶10 000、1∶50 000數(shù)字影像圖(CH/T 1009—2001)行業(yè)標(biāo)準(zhǔn)或規(guī)范為參考。

圖像基本幾何質(zhì)量評價(jià)包括全色多光譜相對配準(zhǔn)精度[13]和圖像定位精度方面的評價(jià)[14]。具體來看主要包括：P5與P10影像范圍一致性、P5和P10影像間波段匹配精度評價(jià)以及無控制點(diǎn)的RPC處理圖像定位精度評價(jià)[15]。

根據(jù)遙感圖像像點(diǎn)位移規(guī)律[16-17]，通過比較控制點(diǎn)坐標(biāo)值(基準(zhǔn)圖幅匹配)和圖像上量算值的偏差，分別評價(jià)X方向和Y方向平均誤差，研究不同遙感圖像、不同地形地物類別的幾何精度。地形地物類別選取平地、丘陵和山地3種地區(qū)[18]。

采用5 m全色影像、已有的實(shí)測高精度幾何控制點(diǎn)以及Worldview-2正射校正影像為基礎(chǔ)，選取待校正影像和基礎(chǔ)底圖上均有的同名明顯特征地物點(diǎn)為控制點(diǎn)。對10 m多光譜影像進(jìn)行幾何校正。校正單元選取整個(gè)單相機(jī)影像，控制點(diǎn)均勻分布，根據(jù)地形需要，適當(dāng)確定控制點(diǎn)數(shù)量。

將波段匹配像元錯(cuò)位數(shù)、影像圖X、Y方向平均誤差、區(qū)域定位平均誤差等作為評價(jià)圖像定位精度優(yōu)劣的指標(biāo)。

2　測試評價(jià)分析

本文選取新疆哈密與云南東川區(qū)域作為評價(jià)示范區(qū)，選取該區(qū)域CEBRS-04星5 m全色和10 m多光譜影像數(shù)據(jù)作為評價(jià)所用數(shù)據(jù)，此外還有實(shí)測大比例尺幾何控制點(diǎn)、0.5 m分辨率Worldview-2正射影像和Google earth影像數(shù)據(jù)作為輔助數(shù)據(jù)。

2.1　P5與P10影像范圍一致性

選取新疆哈密、云南東川PMS數(shù)據(jù)，以P5影像為基準(zhǔn)對P10影像進(jìn)行配準(zhǔn)并套合，檢查配準(zhǔn)后P5與P10影像的范圍一致性。經(jīng)檢查，P5與P10影像的范圍基本一致，如圖1所示。

(a)　新疆哈密P5與P10影像范圍

(b)　云南測試區(qū)一P5與P10影像范圍

(c)　云南測試區(qū)二P5與P10影像范圍圖1　測試區(qū)影像范圍

2.2　P5和P10影像間波段匹配精度評價(jià)

從各測試區(qū)的影像來看，單相機(jī)多光譜數(shù)據(jù)各波段間匹配精度非常高，不存在波段不匹配情況，進(jìn)行彩色合成之后不存在錯(cuò)位、重影等現(xiàn)象。

為定量描述P5與P10影像各波段間配準(zhǔn)情況，對測試區(qū)部分?jǐn)?shù)據(jù)進(jìn)行了定量分析。針對尾號為8004的P10多光譜數(shù)據(jù)，以P5數(shù)據(jù)為基準(zhǔn)，對P5和P10手工選擇了11個(gè)檢查點(diǎn)，計(jì)算結(jié)果表明：所有檢查點(diǎn)X方向殘差平均值為0.000 623像元，Y方向殘差平均值為0.003 387像元，平面方向殘差平均值為0.007 707像元。在采用一次多項(xiàng)式計(jì)算模型下，P5和P10檢查點(diǎn)RMS為0.766 816像元，在二次多項(xiàng)式計(jì)算模型下，P5和P10檢查點(diǎn)RMS為0.015 735像元，波段匹配性非常好。波段匹配精度目視檢查如圖2所示。

圖2　新疆哈密測試區(qū)P5、P10波段匹配精度檢查結(jié)果

云南測試區(qū)一的P10多光譜數(shù)據(jù)中，以P5和P10為例，對P5和P10手工選擇了9個(gè)檢查點(diǎn)，計(jì)算結(jié)果表明：所有檢查點(diǎn)X方向殘差平均值為0.005 005像元，Y方向殘差平均值為0.000 03像元，平面方向殘差平均值為0.005像元。P5和P10檢查點(diǎn)RMS為0.378像元，波段匹配性非常好。波段匹配精度目視檢查如圖3所示。

圖3　云南測試區(qū)一P5與P10波段匹配精度目視結(jié)果

云南測試區(qū)二的P10多光譜數(shù)據(jù)中，以P5和P10為例，對P5和P10手工選擇了8個(gè)檢查點(diǎn)，計(jì)算結(jié)果表明：所有檢查點(diǎn)X方向殘差平均值為0.000 24像元，Y方向殘差平均值為0.004 58像元，平面方向殘差平均值為0.004 507像元，P5和P10檢查點(diǎn)RMS為0.288像元，波段匹配性非常好。波段匹配精度目視檢查如圖4所示。

圖4　云南測試區(qū)二P5與P10波段匹配精度目視結(jié)果

新疆哈密測試區(qū)P5與P10波段匹配精度檢查表如圖5所示，云南測試區(qū)一P5與P10波段匹配精度檢查表如圖6所示，云南測試區(qū)二P5與P10波段匹配精度檢查表如圖7所示。

圖5　新疆哈密測試區(qū)P5與P10波段匹配精度檢查表

圖6　云南測試區(qū)一P5與P10波段匹配精度檢查表

圖7　云南測試區(qū)二P5與P10波段匹配精度檢查表

2.3　無控制點(diǎn)的RPC處理圖像定位分析評價(jià)

主要是利用新疆哈密試驗(yàn)場示范區(qū)黑山口幅系統(tǒng)幾何校正后的2級產(chǎn)品P5全色和P10多光譜影像數(shù)據(jù)，并以野外實(shí)測的42個(gè)控制點(diǎn)和利用實(shí)測控制點(diǎn)生產(chǎn)的Worldview-2正射影像為基礎(chǔ)底圖，進(jìn)行CBERS-04衛(wèi)星P5全色和P10多光譜影像數(shù)據(jù)的幾何精度初步評價(jià)。具體的幾何比對如表1所示。

表1　影像坐標(biāo)與實(shí)測幾何標(biāo)志點(diǎn)坐標(biāo)對比　(s)

經(jīng)查詢相關(guān)資料，在新疆哈密試驗(yàn)場示范區(qū)范圍內(nèi)(42°10′～42°20′ N、94°45′～95°00″ E)，緯度1″≈30.8 m，經(jīng)度1″≈23.05 m，把表1中的單位轉(zhuǎn)換為m并取絕對值后，如表2所示。

由表1和表2可看出，CBERS-04星同名點(diǎn)影像坐標(biāo)與野外實(shí)測坐標(biāo)相比較：

表2　影像坐標(biāo)與實(shí)測幾何標(biāo)志點(diǎn)坐標(biāo)對比　(m)

① P5全色影像在緯度方向上相差0.038"～0.74"，約11.7～22.79 m，平均約17.56 m；在經(jīng)度方向上相差3.24"～4.44"，約74.68～102.34 m，平均約93.05 m；點(diǎn)位誤差77.01～104.59 m，平均為94.82 m，如圖8所示。

② P10多光譜影像在緯度方向上相差1.12"～1.68"，約34.5～54.74 m，平均約28.13 m；在經(jīng)度方向上相差在0.91"～1.79"，約20.98～41.26 m，平均約33.85 m；點(diǎn)位誤差在50.17～63.57 m，平均為57.55 m，如圖9所示。

圖8　P5全色影像幾何精度評價(jià)

圖9　P10多波段影像幾何精度評價(jià)

從圖8和圖9可以看出，對P5全色影像，其經(jīng)度方向上的誤差對整體的點(diǎn)位誤差影響較大；對P10多光譜影像，其緯度方向上的誤差對整體的點(diǎn)位誤差影響較大。

2.4　結(jié)果分析

在對CBERS-04星5 m全色和10 m多光譜影像進(jìn)行無控制點(diǎn)的RPC校正之后，結(jié)合已有的野外實(shí)測控制點(diǎn)位置和由其校正后得到的Worldview-2正射影像，進(jìn)行全色和多光譜影像數(shù)據(jù)幾何定位精度評價(jià)，得到以下結(jié)論：① P5與P10影像范圍一致；② 從各測試區(qū)的影像來看，單相機(jī)多光譜數(shù)據(jù)各波段間匹配精度非常高，不存在波段不匹配情況，進(jìn)行彩色合成之后不存在錯(cuò)位、重影等現(xiàn)象；③ 無控制點(diǎn)的RPC處理影像坐標(biāo)與實(shí)測幾何標(biāo)志點(diǎn)坐標(biāo)對比發(fā)現(xiàn)P5影像經(jīng)度方向上的誤差對整體點(diǎn)位誤差影響較大，P10影像緯度方向上的誤差對整體的點(diǎn)位誤差影響較大。

3　結(jié)束語

遙感圖像在基礎(chǔ)測繪、土地資源調(diào)查等應(yīng)用方面對地面定位、定量精度有較高的要求。本文從全色多光譜相對配準(zhǔn)精度、圖像定位精度方面展開評價(jià)，將波段匹配像元錯(cuò)位數(shù)、影像圖X、Y方向平均誤差、區(qū)域定位平均誤差等作為評價(jià)圖像定位精度優(yōu)劣的指標(biāo)，開展了CBERS-04星5 m全色與10 m多光譜影像范圍一致性、影像間波段匹配精度和無控制點(diǎn)的RPC處理圖像定位精度評價(jià)。5 m全色與10 m多光譜影像范圍一致性較好，10 m多光譜數(shù)據(jù)的波段之間匹配精度高；在圖像定位精度上，5 m全色影像經(jīng)度方向的誤差對整體點(diǎn)位誤差影響較大，10 m多光譜影像的整體點(diǎn)位誤差主要來源于緯度方向上的偏移。

隨著遙感衛(wèi)星的廣泛應(yīng)用，對衛(wèi)星的幾何定位精度將是更大的挑戰(zhàn)。這意味著不僅對衛(wèi)星平臺、飛行姿態(tài)、成像模式的控制程度和穩(wěn)定性提出了更高的要求，而且需要更加精確的影像像元與地面點(diǎn)位的轉(zhuǎn)換模型和地面參考點(diǎn)測量。

[1]張?jiān)鱿椋魹t，溫慶可，等.土地資源遙感應(yīng)用研究進(jìn)展[J].遙感學(xué)報(bào)，2016，20(5)：1243-1258.

[2]夏朝旭，何政偉，于歡，等.面向?qū)ο蟮耐恋馗脖蛔兓瘷z測研究[J].遙感技術(shù)與應(yīng)用，2014，29(1)：106-113.

[3]潘德爐，何賢強(qiáng).海洋水色水溫遙感應(yīng)用科學(xué)與技術(shù)[M].北京：海洋出版社，2014.

[4]陳藝蝦.遙感圖像幾何定位精度評價(jià)方法研究[D].南京：南京理工大學(xué)，2013.

[5]宋偉，朱岱寅，李勇.機(jī)載聚束SAR圖像定位精度研究[J].數(shù)據(jù)采集與處理，2014，29(4)：555-561.

[6]明洋，陳楚江，余紹淮，等.多條帶WorldView衛(wèi)星圖像幾何定位精度分析[J].測繪科學(xué)，2013，38(1)：160-162.

[7]李宗華，王新洲，彭明軍，等.高分辨率衛(wèi)星遙感影像在土地利用變更調(diào)查中的應(yīng)用[J].測繪信息與工程，2005，30(4)：13-16.

[8]白煥成.高分辨率衛(wèi)星遙感影像在基礎(chǔ)測繪中的應(yīng)用探討[J].測繪科學(xué)與工程，2008，28(2)：15-18.

[9]陳正超.中國DMC小衛(wèi)星在軌測試技術(shù)研究[D].北京：中國科學(xué)院遙感應(yīng)用研究所，2005.

[10] 王春媛.遙感圖像幾何校正及目標(biāo)識別技術(shù)研究 [D].哈爾濱：哈爾濱工業(yè)大學(xué)，2014.

[11] 章學(xué)靜.像素級圖像增強(qiáng)及配準(zhǔn)算法研究 [D].北京：北京理工大學(xué)，2014.

[12] 李暉，肖鵬峰，馮學(xué)智，等.結(jié)合光譜和尺度特征的高分辨率圖像邊緣檢測算法[J].紅外與毫米波學(xué)報(bào)，2012，31(5)：469-474.

[13] 朱琳.地理國情普查中遙感影像處理方法的研究[D].西安：西安科技大學(xué)，2015.

[14] 陳亞妮.遙感圖像質(zhì)量評價(jià)方法研究與軟件設(shè)計(jì)[D].杭州：浙江大學(xué)，2017.

[15] 薛萬勛.基于邊緣特征的遙感圖像質(zhì)量評價(jià)方法研究[D].北京：中國科學(xué)院大學(xué)，2016.

[16] 張繼賢，張永紅，林宗堅(jiān).SPOT影像像點(diǎn)位移的研究[J].測繪科學(xué)，2000，25(1)：19-22.

[17] 舒嘉，張林海.航攝像片像點(diǎn)位移的精度分析[J].測繪與空間地理信息，2008，31(3)：155-156.

[18] 戰(zhàn)鷹，史良樹，王金強(qiáng)，等.衛(wèi)星側(cè)視成像引起的像點(diǎn)位移誤差計(jì)算方法[J].河南理工大學(xué)學(xué)報(bào)(自然科學(xué)版)，2015，34(3)：370-373.