HY-2衛星掃描微波輻射計多通道分辨率匹配技術研究

黃 磊，周 武，李延民

（1.國家衛星海洋應用中心，北京 100081；2.中國航天科工集團公司第五研究院西安分院，西安 710100）

HY-2衛星掃描微波輻射計多通道分辨率匹配技術研究

黃 磊1，周 武1，李延民2

（1.國家衛星海洋應用中心，北京 100081；2.中國航天科工集團公司第五研究院西安分院，西安 710100）

星載微波輻射計多通道海洋參數反演需要基于同一位置、同一分辨率的多通道觀測值，但由于星載微波輻射計設計的局限性，多通道觀測地面像元位置和分辨率均不同，因此運用多通道分辨率匹配技術統一觀測面元是掃描微波輻射計數據處理中的關鍵技術，基于海洋二號（HY-2）衛星掃描微波輻射計的成像原理和天線方向圖，模擬各通道天線方向圖在地面的投影，用Backus-Gilbert（BG）算法將高頻觀測亮溫重采樣到最低的低頻6.6 GHz通道的觀測位置和像元分辨率。結果表明，BG降低分辨率重采樣算法能得到很好的擬合效果，并且不另引入噪聲。

掃描微波輻射計；HY-2；分辨率匹配；BG算法

1 前言

被動微波輻射計能夠監測海洋大范圍的變化，相對于主動微波遙感，其傳感器體積和功耗小，在沒有暴雨的情況下，儀器幾乎能全天候透過云層觀測海洋表面，是獲取大氣、海面和陸地等多領域環境信息的重要手段。但微波輻射計受自身系統的限制，其空間分辨率依賴于天線的尺寸，并與頻率相關，受衛星有效載荷體積和重量的限制，星載微波輻射計無法配裝大型天線，不同的頻率的饋源共用一個反射面天線，導致各通道觀測的入射角度不同，對應的地面觀測位置、覆蓋范圍和分辨率也不相同。在進行微波輻射計海洋參數反演時，需要用到多通道亮溫進行綜合計算，這要求多通道地面觀測位置和面元要一致，即將高分辨率的觀測結果降低到低分辨率的尺度，或者將低分辨率的觀測結果提高到高分辨率的尺度，前者會導致高分辨率通道的細節信息丟失，后者會引入大量的噪聲，降低觀測精度。由于海面溫度和海面風場的觀測主要受低頻影響，且為了保證觀測的準確性，所以星載微波輻射計多采用降低分辨率的多通道分辨率匹配技術進行地面像元重構，將同一時刻所有通道的觀測值統一到同一位置和分辨率。

目前星載微波傳感器空間匹配算法主要有Backus和Gilbert提出的BG算法[1]、Sethmann的圖像反卷積技術等。其中BG算法廣泛應用于國外的星載多通道掃描微波輻射計，如SSM/I、AMSR-E等微波輻射計[2，3]。在國內，也有眾多將BG算法應用于微波輻射計的研究[4，5]，并將此算法應用于FY-3衛星MWRI上[4～6]。BG算法通過模擬不同頻率天線增益到地面的投影，通過重采樣算法，利用鄰近觀測像元重疊部分，重構指定位置的觀測結果，并盡可能地接近真實的天線增益方向圖。

HY-2衛星上搭載的掃描微波輻射計是一臺設計有9個通道的輻射計，可以監測海面溫度、海面風場、海冰、降雨等海洋大氣參數，其頻率設置為6.6 GHz、10.7 GHz、18.7 GHz、23.8 GHz、37.0 GHz，除23.8 GHz頻率僅有垂直極化通道，其余頻率均包含水平極化通道和垂直極化通道，輻射計采用圓錐掃描的方式，天線的俯仰角為40°，幅寬1600 km，地面分辨率優于100 km、70 km、40 km、35 km、25 km。根據HY-2掃描微波輻射計的觀測特性，本文應用BG算法進行了多通道分辨率匹配研究，其中涉及觀測幾何定位、天線方向圖地面投影、目標方向圖重構，并對結果進行了分析，根據評價結果給出了HY-2掃描微波輻射計多通道分辨率匹配的方案。

2 BG反演算法

BG算法是由Stogryn將Backus-Gilbert矩陣反演方法引入星載微波輻射計系統提出的。根據Stogryn[7]，用G（ρA，ρ）表示地球坐標系下中心位置在ρA的增益方向圖ρ位置的有效天線增益，是采樣積分時間τ和瞬時傳感器方向圖Gl的函數

式（1）中，s^0（t）表示t時刻的天線視線方向；s^（t）表示傳感器到位置 ρ的單位矢量。 ρA對應的有效天線指向為

則第i次觀測亮溫則是地球表面亮溫在天線增益上的積分

BG反演算法是尋找相鄰觀測值TBi一組最優的線性組合，重構出一個最接近于真實測量值TBC

根據最小二乘原理，使得重構亮溫與參考亮溫最為接近必須滿足3個條件，首先是G為歸一化，即ai的總和為1，第2個約束條件Stogryn將其稱為分辨率價值函數

式（5）中，J為補償函數，用來產生理想的天線方向圖，通常將J設置為1。第3個約束條件是重構觀測亮溫的噪聲最小，實際觀測亮溫中疊加有隨機噪聲信號，噪聲的方差為（ΔTi）2，則重構亮溫的方

式（6）中，a是元素ai的列向量；E是誤差協方差矩陣。

為了在噪聲和重構分辨率之間尋求折中，Stogryn建立了分辨率噪聲權衡函數

式（8）中，

式（9）中，I為單位矩陣，G為N×N 對稱矩陣，表示參與采樣足印面元的重疊度。

3 BG算法在HY-2掃描微波輻射計數據上的實現

針對于HY-2衛星掃描微波輻射計各頻率觀測面元不同，運用BG算法來實現對各頻率數據的分辨率匹配計算。大致的過程包括微波輻射計定位、各通道天線方向圖地面投影、BG匹配算法重采樣系數生成和面元重構結果評價分析。

3.1 HY-2掃描微波輻射計觀測幾何

通過HY-2掃描微波輻射計儀器的姿態和衛星的軌道位置，可以推算出各頻率在地面的足印位置，主要包括由傳感器坐標系-平面大地坐標系-地心坐標系的坐標系轉換和俯仰、航偏、滾動等姿態角度的旋轉，結合各頻率的天線方向圖，即可模擬各通道在地面的投影區域。

HY-2掃描微波輻射計上各個饋源通過反射面接收輻射能量，由于多個饋源接受信號，各通道觀測角度不同，所以同一時刻，各頻率在地面的投影位置也不相同，如圖1所示，37 GHz地面觀測中心點位置與6.6 GHz觀測面元中心位置并不重合。此外，由于HY-2掃描微波輻射計低頻6.6 GHz，10.7 GHz與高頻18.7 GHz、23.8 GHz、37 GHz的饋源分別位于飛行方向的兩側，高頻通道與低頻通道觀測的覆蓋范圍也不一樣，如圖2所示，它們只有中間143個觀測面元是重疊覆蓋并包含所有觀測通道的，每個掃描行的初始和結尾部分，只有高頻觀測或者低頻觀測，這部分觀測區域無法反演全部的海洋大氣參數。中間重合的部分即是需要進行空間分辨率匹配的區域。

圖1 6.6 GHz和37 GHz觀測面元地面投影位置分布Fig.1 The observation area center arrangement of 6.6 GHz and 37 GHz

圖2 6.6 GHz和37 GHz覆蓋范圍重疊圖Fig.2 The overlap of observation coverage between 6.6 GHz and 37 GHz

3.2 天線方向圖地面投影

在獲得準確的定位算法后，按各通道的觀測幾何，將天線方向圖投影到地面，方向圖地面投影的采樣分辨率為1 km，選取沿軌道方向31行，掃描向150列，天線方向圖選取增益下降－10 db的區間，生成的各通道天線方向圖地面投影的結果如圖3所示。

圖3 地面足印Fig.3 Simulated scan antenna footprint

3.3 BG匹配算法重采樣系數計算

由生成的天線方向圖地面足印，用BG算法進行通道分辨率降低匹配計算，由于微波輻射計6.6 GHz頻段為海面溫度反演的關鍵通道，而且在多通道綜合反演算法中，大氣和海洋參數的反演均需要用到低頻部分，所以本文將6.6 GHz的天線方向圖作為參考方向圖。用18.7 GHz的垂直極化（V）天線方向圖擬合得出6.6 GHz V極化天線方向圖，用37 GHz的V極化天線方向圖擬合得出6.6 GHz V極化天線方向圖，如圖4所示。根據天線方向圖地面高低頻足印的重疊范圍，本文采用15×15的重采樣窗口。由圖4可知，對分辨率降低的計算，BG算法能很好地用高分辨率天線方向圖擬合低分辨率天線方向圖。

圖4 極化天線方向圖Fig.4 Polarization antenna pattern

BG算法可同時進行通道分辨率降低和增強兩種運算，但進行分辨率增強會引入噪聲，根據楊虎的研究[5]，針對掃描微波輻射計，利用低頻數據進行重采樣，不能完全擬合出高頻分辨率的天線方向圖，本文對此就不再進行計算分析。

3.4 擬合系數應用及評估

為檢驗BG算法，本文將重采樣系數應用于HY-2掃描微波輻射計L1B亮溫數據中。由于HY-2掃描微波輻射計有高頻和低頻兩組饋源指向，造成了高頻數據和低頻數據在地面的投影僅有中間的134列重合，另在掃描方向上，各掃描點的位置不同，觀測角度有所變化，則相鄰面元的相互關系不同，因此需要對掃描行中的各點分別選取不同的擬合系數。圖5表示的是37 GHz垂直極化通道觀測亮溫匹配到6.6 GHz垂直極化通道分辨率的情況。從圖5可看出，經過圖像降低分辨率的處理后，圖像的細節信息有所弱化，但是仍然保留了其通道原有的觀測特點，能反映出海上的水汽分布。

圖5 37 GHz V極化匹配到6.6 GHz V極化分辨率的結果Fig.5 37 GHz V matched to 6.6 GHz V resolution image

4 結語

本文詳細闡述了BG算法的原理，并針對HY-2掃描微波輻射計各觀測通道地面觀測位置和足印覆蓋范圍不匹配的情況，以天線方向圖為基礎，利用BG算法進行了降分辨率足印匹配研究，結果表明，通過合理的參數選取，能很好地將高頻通道的觀測結果擬合到低頻通道的觀測，且不引入噪聲，這對后期掃描微波輻射計海洋大氣參數的反演將起到十分重要的作用。

[1]Backus G，Gibert F.Numerical applications of a formalism for geophysical inverse problem[J].Geophysical Journal Royal Astronomical Society，1967，13：247-276.

[2]Farrar M R，Smith E A.Spatial resolution enhancement of terrestrial features using deconvolved SSM/I microwave brightness temperatures[J].IEEE Transactions on Geoscience and Remote Sensing 1992，30（2）：349-355.

[3]Robinson W D，Kummerow C，Olson W S.A technique for enhanceing and matching the resolution of microwave measurements from the SSM/I instrument[J].IEEE Transactions on Geoscience and Remote Sensing，1992，30（3）：419-429.

[4]李 靖，王 濤，張升偉，等.利用BG算法提高微波輻射計空間分辨率[J].遙感學報，2004，8（5）：409-413.

[5]楊 虎，商 建，呂利清，等.星載微波輻射計通道分辨率匹配技術及其在FY-3衛星MWRI中的應用[J].上海航天，2012（1）：23-28.

[6]尹紅剛，張德海.星載微波輻射計空間分辨率增強算法研究[J].海洋科學進展，2004，22（增）：192-197.

[7]Stogryn A.Estimates of brightness temperatures from scanning radiometer data[J].IEEE Transaction Antennas and Propagation，1978，26（5）：720-726.

Study on channel resolution matching algorithm for HY-2 satellite scanning microwave radiometer

Huang Lei1，Zhou Wu1，Li Yanmin2

(1.National Satellite Ocean Application Service，Beijing 100081，China；2.The Fifth Institute Xi’an Branch of China Aerospace Science and Industry Corporation，Xi’an 710100，China)

Accuracy satellite scanning microwave radiometer ocean parameters retrieval needs brightness temperature of multiple channels from the same area and identical resolution.Due to the limitation of spaceborne radiometer antenna design and feedhorn arrangement，each channel has different resolutions and observation positions，which will increase the ocean parameter retrieval error.Based on the HY-2 scanning microwave radiometer antenna pattern and imaging geometry，we simulate the ground footprint form different channels，use Backus-Gilbert(BG)algorithms to produce a single composite sample at a particular location and with a particular spatial weighting pattern from combination of the adjacent measurements；the low frequency 6.6 GHz is set as the reference channel.The result shows that the matching with resolution decreasing could simulate the real instrument observation without introducing noise.

scanning microwave radiometer；HY-2；resolution matching；Backus-Gilbert algorithm

V443

A

1009-1742（2014）06-0065-05

2014-04-10

國家高技術研究發展計劃（“863”計劃）項目“海洋動力環境微波遙感信息提取技術與應用”（2013AA09A505）

黃 磊，1982年出生，男，湖北潛江市人，助理研究員，研究方向為微波遙感；E-mail：huanglei@mail.nsoas.gov.cn