基于升降軌數據融合的月球微波輻射亮溫制圖

陳興梅，連 懿，曾曉明

（天津師范大學 地理與環境科學學院，天津300387）

在新一輪月球探測活動中，嫦娥一號（CE-1）和嫦娥二號（CE-2）搭載的多波段微波輻射計（microwave radiometer，MRM）具有良好的空間覆蓋度，為分析月表微波亮溫特征、研究月球淺表層物質的物理化學和地理分布特性提供了良好的數據支撐[1-5].基于嫦娥微波輻射計數據的全月微波亮溫制圖是分析月球微波輻射亮溫時空變化特征的關鍵數據基礎.

近年來，全月微波亮溫的制圖工作取得了大量的研究成果.Fa 等[6]基于嫦娥一號微波輻射計數據，利用觀測中太陽入射角的最近鄰插值法構建了月球白天和夜間的全月微波亮溫分布圖，并討論了月球表面亮溫與緯度、FeO+TiO2含量等的關系.王振占等[7]對嫦娥一號數據進行在軌兩點定標處理，獲得全月的亮溫分布圖并對其特點進行分析，但該研究在制圖過程中沒有考慮天文時間和光照條件等因素. 針對該問題，Chan 等[8]提出利用時角來校正亮溫數據，并基于嫦娥一號微波輻射亮溫數據從宏觀角度分析了月球表層的熱輻射特性.Zheng 等[9]基于時角概念進一步構建亮溫日變化模型，并獲取了第一張不同時刻嫦娥一號微波亮溫數據的全月分布圖，校正后的數據比以往數據表現出更豐富的信息.在此基礎上，連懿等[10]利用嫦娥二號數據，采用時角計算和亮溫日變化模型反演月壤介電常數，并制作出全月微波亮溫分布圖.朱永超等[11]基于嫦娥二號數據，采用反距離差值法繪制形成12幅月球亮溫圖，但由于部分數據異常以及無數據等原因導致部分亮溫圖數據丟失.

在微波輻射計數據采集過程中，不同軌道高度造成的信號衰減差異和不同儀器姿態造成的差異均會造成升降軌數據不一致，進而可能會影響微波輻射亮溫制圖的精度，從而對月球微波輻射亮溫時空變化特征研究產生影響. 現有研究均未考慮嫦娥衛星升軌（一側軌道從南向北飛）與降軌（一側軌道從北向南）的差異性.因此，本研究基于嫦娥一號和嫦娥二號搭載微波輻射計數據特征和時角校正，對嫦娥一號和嫦娥二號數據質量及其數據覆蓋度的差異進行對比分析，并利用亮溫日變化模型和四次線性擬合模型對嫦娥二號衛星所得亮溫數據進行不同緯度的擬合校正及亮溫圖制作，進而提高嫦娥二號升降軌數據的一致性，并在一定程度上減少數據誤差，為后續月表研究提供更加準確的數據基礎.

1 數據來源與處理

1.1 嫦娥微波輻射計

多波段微波輻射計是嫦娥一號和嫦娥二號衛星搭載的有效載荷之一，是利用被動接收的各個高度的溫度輻射微波信號判斷溫度曲線，從而定量測量月球表面溫度的高靈敏度接收裝置.嫦娥一號衛星攜帶的微波輻射計（CE-1 MRM）具有4個通道，其頻率分別為3.00、7.80、19.35 和37.00 GHz，波長分別為10.00、3.84、1.55 和0.81 cm[12]，4個通道空間分辨率分別約為50.00、35.00、35.00 和35.00 km[13]，觀察角和測量精度均為0°和0.50 K[10].CE-1 衛星（軌道高度為200 km）在2007 年11 月—2008 年7 月共獲得1307 軌微波輻射計亮溫數據[13-14].

嫦娥二號衛星攜帶的微波輻射計（CE-2 MRM）具備的4個通道頻率和波長與嫦娥一號衛星相同，4個通道空間分辨率分別約為25、15、15 和15 km，輻射分辨率分別約為0.16 K、0.16 K、0.12 K 和0.14 K，測量精度均為0.5 K[14-15]. CE-2 衛星（軌道高度為100 km）在2010 年10 月—2011 年5 月共獲得2394 軌微波輻射計亮溫數據[16-17]，軌道周期為118 min，在一個軌道周期內，CE-2 衛星升軌（衛星一側軌道從南向北）獲取觀測數據1825個，降軌（衛星另一側從北向南）獲取相同數量的觀測數據[11].

本文采取探月工程地面應用系統發布的CE-1 MRM 和CE-2 MRM 2C 級3.0 GHz 微波亮溫數據.2C級軌道數據包括了數據采樣時間、4個頻率通道的亮溫、星下點太陽入射角和方位角、月球表面經緯度以及軌道高度等信息，其中星下點太陽入射角和方位角以及月球表面緯度可以用于計算月球時角，即月球當地時間[18].

1.2 數據處理

月表物理溫度是月表微波輻射亮溫的重要影響因素，而嫦娥衛星的月球采樣時間在很大程度影響著采樣時刻的物理溫度，從而對數據的一致性產生影響.為了消除月球時間對亮溫的影響，也為了便于后續的數據評價及亮溫校正，本研究引入時角的概念，對嫦娥微波輻射計數據進行時間校正處理[8-9].

時角是天文學名詞，指天子午圈與天體的赤經圈在北極所成的球面角，或在天赤道上所夾的弧度.其中1 h 表示15°，根據太陽的入射角和方位角以及經緯度等幾何定位數據可以計算出月球時間：

式（1）中：i 為太陽入射角，取值范圍為[0，π]；a 為太陽方位角，取值范圍為[0，2π]，可以從2C 級數據中直接獲取；φ 為負極角[-π，π]；λ 為緯度，取值范圍為[-π/2，π/2].

本研究利用式（1）計算獲得亮溫數據的月球時角.隨后在以20°緯度間隔進行區域劃分的基礎上[9]，對嫦娥亮溫數據按每5°范圍設置采樣緯度間隔，將全月球表面分為18個區域，對每個緯度區域建立亮溫日變化模型.模型采用7 次線性公式擬合參數，即

式（2）中：h 為時角；TB（h）為時角為h 時的亮溫值；ai（i =0，1，…，7）為7 次線性擬合模型參數變量.當h=0°時，即為正午時刻；當h=180°時，為午夜時刻.

基于該模型將不同時間的采樣數據擬合到特定的觀測時刻，獲得特定時刻的亮溫數據，降低了由采樣時刻不同造成的月表溫度差異對嫦娥數據質量評價和升降軌數據校正的影響.

2 嫦娥微波輻射計數據評價

CE-1 和CE-2 衛星在軌期間對月探測獲取了大量的微波輻射計數據，本研究從其數據覆蓋度和數據質量2 方面對CE-1 和CE-2 微波輻射計數據間的差異性進行定量評價.

2.1 數據覆蓋度評價

本研究以0.5 K 的分辨率將月表劃分為720×360的規則格網，對CE-1 和CE-2 的亮溫數據在全月面的空間覆蓋狀況進行統計，圖1 為嫦娥衛星獲取微波輻射計數據的空間分布情況.

圖1 CE-1 MRM 和CE-2 MRM 的數據分布Fig.1 Data distribution of CE-1 and CE-2 microwave radiometer

由圖1 可以看出，CE-1 和CE-2 的覆蓋密度隨著緯度的升高而增加，在兩極最高，這可能是由嫦娥衛星的軌道特性決定的；且CE-2 的數據覆蓋度明顯優于CE-1；該分辨率下網格內CE-2 的重復采樣數量超過10個和超過20個的區域較CE-1 更多，其數據的空間覆蓋度也大于CE-1 衛星.嫦娥衛星的采樣數據統計結果如表1 所示.由表1 可知，CE-1 衛星的數據覆蓋率為99.76%（＞0個），CE-2 衛星的數據覆蓋率達到99.90%（＞0個）.其中CE-2 MRM 數據覆蓋率（＞20個）相比CE-1 衛星增加了44.52%，對應的亮溫數據量增加了115404個（＞20個）.綜合分析統計數據可知，CE-1 網格內重復采樣數據的個數主要集中在20個以下，而CE-2 衛星每個網格的采樣數據個數主要集中在20個以上，覆蓋率約為CE-1 衛星的2 倍. 因此，相比CE-1 衛星，CE-2 衛星獲取的亮溫數據具有更高的時間覆蓋率和空間覆蓋度. 同時由于CE-2 衛星相比CE-1 衛星的飛行高度低100 km，因此微波探測器天線波束在月面的覆蓋就會縮小，其空間分辨率也優于CE-1號.

表1 CE-1 MRM 和CE-2 MRM 的采樣數據統計Tab.1 Sampling data statistics of CE-1 and CE-2 microwave radiometer

2.2 數據質量評價

由于衛星飛行的方向和軌道高度不同，相同采樣點在同一月球時刻的亮溫值也可能存在較大差異.本研究從其升降軌數據角度對數據質量進行評價.首先利用式（2）的亮溫日變化模型實現不同時刻亮溫數據的校正，然后將亮溫數據全部歸一化到正午時刻（時角為0°），再對擬合后的亮溫數據進行統計分析，得到正午時刻全球亮溫升降軌數據的2D 散點圖，結果如圖2 所示.

圖2 CE-1 MRM 和CE-2 MRM 的數據質量對比Fig.2 Data quality comparison of CE-1 and CE-2 microwave radiometer

由圖2 可以看出，CE-1 的升降軌數據變化存在分帶性且較為離散，而CE-2 升降軌數的一致性相對較好，基本為線性分布.本研究利用相關系數r 評價比較兩者間差異性的大小，r 的絕對值越大表明相關性越強.CE-1 的相關系數r 為0.9356，CE-2 的相關系數r 為0.9523，高于CE-1.這可能是因為與CE-1 衛星軌道相比，CE-2 衛星飛行軌道距離月球表面更近，因此受電磁波輻射和光照幾何等因素的影響較小，數據質量更好.

綜合數據覆蓋評價和數據質量評價結果可知，CE-2 衛星數據在空間覆蓋度和數據質量方面均優于CE-1 衛星.因此，本研究中采用CE-2 MRM 數據進行后續的微波輻射亮溫數據校正和微波亮溫制圖工作.

3 CE-2 MRM 升降軌數據校正與亮溫制圖

3.1 不同緯度擬合的升降軌校正

基于統計學思想，通過實驗選取了四次線性擬合模型，用于消除不同軌道間的亮溫差異.由于亮溫受到月表物理溫度的影響，具有緯度地帶性，因此針對不同緯度分別建立擬合模型，提高數據擬合的精度，該方法能夠在一定程度上降低升降軌數據校正過程中產生的過擬合現象，并最大程度地還原數據本身.

將CE-2 衛星的升軌和降軌數據按高緯度（0°～60°和300°～360°）和低緯度（60°～300°）劃分，對不同緯度的亮溫數據建立四次線性擬合模型，模型的形式為

式（3）中：M 為多項式中的最高次冪；xj為x 的j 次冪；wj為xj的系數.

利用該模型可以實現亮溫數據升軌到降軌的數據轉換，同時利用相關系數評價方法對轉換后的升降軌數據進行數據質量評價.結果表明CE-2 升降軌數據擬合校正后相關系數r=0.9812，而未校正CE-2 升降軌數據的相關系數r=0.9714，這說明針對緯度亮溫數據建立的四次線性擬合模型可以在一定程度上減少升降軌數據的差異.圖3 為分緯度模型校正后升降軌數據的2D 散點圖.由圖3 可以看出，分緯度校正后升降軌數據的分布趨于一條直線，表明數據質量有所提升.

3.2 全月的亮溫制圖

基于校正后的CE-2 MRM 亮溫數據，以0.5°的分辨率將月表劃分為360 行×720 列的規則網格進行全球亮溫投影，將月球表面亮溫全部歸一化到正午時刻（時角為0°），結合亮溫日變化線性擬合模型[11]，實現了升降軌數據融合的全月微波亮溫制圖，結果如圖4 所示.

圖3 CE-2 分緯度校正后亮溫分布Fig.3 CE-2 TB distribution after correction of different latitudes

圖4 校正前后正午時刻的亮溫分布Fig.4 Brightness temperature distribution before and after calibration at noon

由圖4 可以看出，校正前后的亮溫數據均具有緯度地帶性，這是因為緯度越高月球表面物理溫度越低，緯度越低月球表面物理溫度越高.此外，由于3 GHz電磁波穿透深度大，地形對亮溫的影響并不明顯，只有模糊的輪廓.受溫度的緯度地帶性影響，極地地區亮溫分布較為明顯，從極地地區向中間地帶的顏色過渡較大.

對比圖4（a）和圖4（b）可知，校正前30°W 和120°E附近有幾條比較明顯的條帶，校正后條帶明顯減少，升降軌融合的方式有效改善了極地地區的亮溫數據質量，拓展了亮溫數據的分析范圍.為了進一步比較不同緯度校正前后的亮溫，本研究提取cabeus 撞擊坑前后范圍（115°S～155°S，70°W～90°W）的數據進行放大對比.對于所選取cabeus 撞擊坑范圍內的地形，由于亮溫誤差較大，在校正前的局部圖中該區域內的地形不能明顯的體現，地形信息幾乎被掩蓋.經過不同緯度擬合校正后，數據間的差異明顯減少，制作所得亮溫圖中亮溫色彩一致性更好，投影所呈現出的地形更接近cabeus 撞擊坑的特征，這有利于地形的判讀和后續的月表分析.

4 結論

（1）本研究從數據分布和數據質量兩方面對CE-1和CE-2 微波輻射計2C 級亮溫數據進行差異性對比分析，結果表明CE-2 的數據質量和數據覆蓋度均優于CE-1.

（2）在時角校正的前提下，基于CE-2 號微波輻射計數據采用不同緯度分別擬合的方式對升降軌數據進行校正.校正后CE-2 號升降軌數據間的相關性有所提升.結合亮溫日變化模型分別制作出校正前后全月亮溫圖和局部區域亮溫圖.結果表明校正后的全月亮溫圖中右側及中間區域的條帶明顯減少，局部地區地形信息明顯，亮溫色彩一致性較好.

（3）基于統計學方法對升降軌數據進行不同緯度擬合校正，所得升降軌數據的誤差較小，提高了全月微波亮溫圖的制圖精度，為后續月球微波遙感數據處理以及通過微波遙感反演月壤厚度、成分等的研究提供了數據基礎.

本研究僅對嫦娥二號亮溫數據進行分緯度擬合，擬合時是否考慮按時角劃分并對不同時角建立不同擬合方程進行擬合仍需進行進一步討論.