顧及陸湖反射差異的衛星測高監測湖泊水位的波形分析與重定

田山川¹ , 郝衛峰¹ , 李斐^1,2,3 , 羅天文²     

1. 武漢大學中國南極測繪研究中心, 湖北 武漢 430079; 
2. 武漢大學測繪遙感信息工程國家重點實驗室, 湖北 武漢 430079; 
3. 武漢大學國家領土 *** 與海洋權益協同創新中心, 湖北 武漢 430079

基金項目：湖北省自然科學基金（2015CFB500）

之一作者簡介：田山川(1994-), 男, 工程碩士, 研究方向為衛星測高。E-mail: shanchuantian@163.com

微信號：MeetyXiao
添加微信好友, 獲取更多信息
復制微信號

通信作者：郝衛峰, E-mail: haowf@whu.edu.cn

摘要：以Jason-2衛星雷達測高數據監測洪澤湖水位為例，分析了湖泊水位沿衛星軌跡呈"V"字形分布的原因。分析發現：由于湖泊面積相對較小，受陸地影響，測高波形中陸地反射信號占主導地位，水面反射產生的上升前緣非常小，導致波形跟蹤出錯。解釋了波形形成的機制、造成水位計算錯誤的原因，結合遙感影像證明了解釋的正確性，并給出了相應的波形重定算法。結果表明，該重定算法優于常用的波形重定算法，剔除異常值后水位標準差為0.09 m，為湖泊區域的測高數據處理提供了一種參考。

關鍵詞：波形分析    反射差異    波形重定    湖泊水位    

Waveform Analysis and Retracking of Lake Level Monitoring by Satellite Altimeter Considering the Difference Between Land and Lake Reflection

TIAN Shanchuan¹ , HAO Weifeng¹ , LI Fei^1,2,3 , LUO Tianwen²     

Abstract: The phenomenon that the distribution of water level along the satellite track is like a "V" is *** yzed.To show the problem, the water level of Hongze Lake is calculated by Jason-2 satellite altimeter data.The *** ysis shows that the signal of ground reflection dominates the shape of waveform due to the *** all size of the lake and influence of land.In this case, the leading edge produced by water is not obvious, so it is easy to retrack the wrong middle point.The mechani *** of waveform formation and the causes of water level calculation errors is explained and proved with remote sensing images.And the corresponding waveform retracking algorithm is proposed.The results indicate that the proposed algorithm performs better than common waveform retracking algorithm.After removing the abnormal value, the standard deviation of the water level reaches 0.09m, which means the *** ysis method could be a reference for altimeter data processing in lake region.

Key words: waveform *** ysis     reflection differences     waveform retracking     lake water level   

湖泊對氣候變化和人類活動較為敏感，它的任何變化都會對自然生態系統及社會經濟系統產生影響。水位是反映湖泊變化的一項重要指標，傳統水位測量主要由水文站完成，其測量精度高，可以連續獲取湖泊水位，但不足之處也很明顯，水文站維護成本高，偏遠欠發達地區的水文站還不能實現全覆蓋，各個地區數據難以共享，給相關科學研究帶來了一定的障礙^[1]。近年來，衛星測高技術成為監測湖泊水位變化的一種新手段，是湖泊水位持續監測的潛在重要數據源。

由于衛星測高技術最初是為海洋任務設計的^[2]，多采用雷達高度計^[3]，故其波束較寬，地表覆蓋面積大，所測高度為衛星至該區域的平均高度。海面相對平坦，采用雷達高度計是合適的，并且雷達高度計可提供額外的信息，計算海面有效波高、后向散射系數、風場等^[4, 5]。但是，將其應用于內陸湖泊時可能會產生一些問題。內陸湖泊面積普遍較小，雷達波可能會觸及湖泊以外的地物，接收到的雷達回波會受到污染，導致距離跟蹤錯誤。在計算湖泊水位時有這樣一種現象，湖泊水位沿衛星軌跡呈如下分布：從湖邊到湖中心降低，再從湖中心到湖邊升高，更低點和更高點高程之差更大可超過10 m，水位沿衛星軌跡呈“V”字形分布。由于水位并非在某一高程上下波動，數據編輯、粗差剔除等策略無法使用^[6-8]，直接利用現有波形重定算法也無法得到較好結果^[9-11]。有文獻對這種情況進行了分析，認為這種現象是由衛星指向角偏差造成的，并利用拋物線擬合的方式對水位進行改正^[12, 13]，雖然其做法能夠較好地對數據進行改正，但并未從波形形成的機制解釋這種現象形成的原因。為了能從實質上解決該問題，需要從波形形成機制出發，探究水位沿衛星軌跡呈“V”字形分布的原因及解決方案。

為了突出“V”字形分布特征，本文以洪澤湖為研究對象，利用Jason-2雷達測高數據計算其水位。首先對洪澤湖的Jason-2測高數據進行處理，展示水位沿衛星軌跡呈“V”字形分布的現象，然后解釋陸地對測高波形的影響及其導致測高波形跟蹤錯誤的原因，指出由水面和陸地反射分別造成的上升前緣在波形中的位置，并結合Landsat影像從波形形成機制證明解釋的正確性。由此設計波形重定算法，從復雜的波形中提取出正確的波形前緣中點，得到正確的湖泊水位。

1 湖泊水位計算

洪澤湖(32°45′N—33°45′N，117°45′E—119°00′E)為中國第四大淡水湖，位于江蘇省西部淮河中下游，湖面面積約2000 km²。Jason-2衛星于2008年6月發射，工作至今，其搭載的雷達高度計的軌跡正好穿越該湖泊(圖 1)，軌跡長度為12 km。Jason-2標稱軌道高為1336 km，軌道傾角66°，軌道重復周期約10 d(文獻[14])。Jason-2數據來自于法國海洋衛星歸檔、驗證及解釋數據中心(Archivage，Validation et Interprétation des données des Satellites Océanographiques，AVISO)發布的SGDR數據集。

圖 1 Jason-2經過洪澤湖的地面軌跡Fig. 1 The track of Jason-2 over Hongze lake

由測高數據計算湖泊水位高的公式如下

 (1)

式中，H為湖面相對于大地水準面的高程，即水位；Alt為衛星軌道的參考橢球高(WGS-84)；R為衛星至湖面的距離(測距值)；G為大地水準面相對于參考橢球高，由重力場模型EGM96得到；V為各項改正值之和，包括電離層改正、對流層干分量和濕分量改正、固體潮改正以及極潮改正等。Jason-2搭載Poseidon-3雙頻高度計，兩種頻率分別為Ku波段(13.575 GHz)和C波段(5.3 GHz)，通常Ku波段測距精度高于C波段，故Ku波段測距值使用較多^[15]。SGDR數據集的測距值有兩種采樣頻率，分別為1 Hz和20 Hz，1 Hz的測距值由20 Hz的測距值線性回歸得到，本文使用20 Hz的Ku波段測距值計算湖泊水位。

衛星接收雷達回波后對波形進行跟蹤獲得測距值，使用不同的跟蹤算法得到的測距值可能會有差異。Jason-2在獲取波形后使用默認跟蹤算法即可解算出測距值，該測距值在數據中對應的字段名為tracker_20hz_ku。衛星將數據傳輸至地面站后對波形進行重新處理，該過程稱為波形重定，用戶獲取衛星SGDR數據后也可進行波形重定以獲取精度更高的測距值。SGDR數據中提供了3種不同波形重定算法得到的測距值，字段名稱分別為：range_20hz_ku、ice_range_20hz_ku、range_20hz_ku_mle3。利用以上不同測距值進行計算將得到不同的湖泊水位，本文計算了cycle0-300的沿衛星軌跡水位，但為了便于觀察，只顯示疊加cycle0-100的結果，圖 2展示了不同測距值得到的洪澤湖水位沿衛星軌跡的分布。可以看出，不同測距值計算得到的水位各有不同，但都有同樣一個現象：大多數沿衛星軌跡水位從低緯度到高緯度(從湖邊到湖中心再到湖邊)先下降再上升，水位沿衛星軌跡呈“V”字形分布，更低點和更高點的差值可以超過10 m。實際上，同一條沿衛星軌跡的湖泊水位應該基本一致，而此處的水位變化明顯與常理不符。需要說明的是，緯度范圍兩端數據的波形中水面和陸地反射信號夾雜在一起無法分離，在實際的計算中會被舍棄，但為了完整地展示水位分布特征，此處給出了整個范圍內的水位分布。

圖 2 不同跟蹤 *** 的測距值計算得到的洪澤湖水位沿衛星軌跡分布
Fig. 2 Water level of Hongze Lake along the satellite track by different tracking methods

2 湖泊沿軌水位呈“V”字形原因分析

衛星自帶的波形跟蹤器是針對海洋波形設計的，當雷達波照射區域為海洋或者大型水域時，測高回波波形類似海洋波形，跟蹤器能夠較好地跟蹤波形。由于湖泊面積較小，測高波形易受地形和地物的影響，回波波形不再與海洋波形類似，原有跟蹤器無法準確跟蹤，從而導致水位計算錯誤。盡管地面處理中心給出了幾種不同波形重定算法跟蹤得到的測距值，當實際回波類型與算法對應的回波類型不符時仍無法有效跟蹤波形。

為了說明波形跟蹤錯誤的成因，以Jason-2 cycle31在洪澤湖上的數據為例，用ice_range_20hz_ku測距值進行水位的計算，其水位沿衛星軌跡分布如圖 3(a)所示。此時水位有明顯的“V”字形現象，其對應的回波波形(圖 3(b))與海洋波形(圖 3(d))明顯不同，海洋波形只有一個明顯的上升前緣，且波形后緣下降緩慢。將水位更低點的波形單獨取出(圖 3(c))，可以看出隨著采樣門的增加，該波形先產生一個小的上升前緣，隨后產生一個大的上升前緣，波形后緣迅速下降。基于實測數據比較，可以看出先產生的小上升前緣由水面反射導致，大上升前緣主要由陸地反射導致，波形重定時跟蹤小上升前緣才能得到正確的測距值。需要說明的是，小上升前緣極易被當作熱噪聲剔除，傳統算法錯誤地跟蹤大上升前緣將使跟蹤到的波形前緣中點靠后，因此計算的測距值偏大，導致水位偏低。

圖 3 Jason-2 cycle31沿衛星軌跡水位及波形
Fig. 3 Water level and waveform along satellite track of Jason-2 at cycle31

雷達高度計波束寬度為1°左右，測量時發射一束脈沖寬度為3.125 ns的脈沖^[16]。脈沖發出后向下傳播首先接觸水面(圖 4(a))，對應圖 4(b)中波形的小上升前緣(灰 *** 域)，隨后脈沖在水面均勻擴散，波形基本平穩。隨著脈沖的繼續擴散，接觸到陸地(圖 4(c))，對應圖 4(d)中波形的大上升前緣(灰 *** 域)，隨后脈沖繼續擴散，接收到的能量減弱，波形后緣呈下降趨勢。

圖 4 小前緣波形形成機制Fig. 4 Mechani *** of *** all front edge waveform formation

上述波形形成過程有一個特點，水面返回能量導致的波形上升量明顯小于夾雜陸地信號的波形上升量。為解釋此現象，圖 5給出完全由水面和完全由陸地反射的波形，其中，圖 5(a)為一段海洋波形，完全由水面反射；圖 5(b)為洪澤湖附近完全由陸地反射的波形。可以看出，海洋波形的更大回波強度不超過300，多數在150~300之間。洪澤湖附近的陸地因地物復雜，波形回波能量峰值波動較大，但多數波形回波能量峰值接近1400，遠大于海洋波形。正因如此，在水陸信號混合的波形中，水面反射形成的上升前緣峰值小，陸地反射形成的上升前緣峰值大，最終產生了具有小前緣的波形。

圖 5 海洋和陸地波形Fig. 5 Wave of ocean and land

如果將波形中小上升前緣當作熱噪聲去除，則極易跟蹤陸地反射形成的上升前緣，由此計算水位的方式如圖 6所示。衛星距離水面的實際距離為R，波形跟蹤時將陸地反射信號誤以為水面跟蹤得到距離R′，將其反算星下點水位時會產生大小為Δh的偏差，Δh=R′－R。可以看出，當星下點離陸地越遠時，產生的水位偏差越大，計算得到的水位越低，與前述湖泊水位沿衛星軌跡呈“V”字形分布現象相符。

圖 6 水位計算錯誤示意圖Fig. 6 Schematic diagram of wrong calculation of water level

為驗證本文分析的正確性，計算衛星星下點與陸地反射點的距離，并利用Landsat影像做驗證。有小上升前緣的波形中，認為小上升前緣由星下點水面反射，隨后的大上升前緣主要由陸地反射，那么跟蹤兩處上升前緣得到的前緣中點做差，即得到脈沖接觸水面和陸地的時間差Δt，由此可計算Δh=Δt·c，其中c為光速，再由式(2)計算星下點與陸地反射點的距離r

 (2)

Landsat影像由美國地質勘探局(United States Geological Survey，USGS)免費提供，下載一幅與Jason-2 cycle31時間最接近的影像，時間相差5 d，以影像為背景，以衛星星下點為圓心，星下點與陸地反射點的距離r為半徑作圓，結果如圖 7所示。

圖 7 Landsat影像驗證結果Fig. 7 Results of Landsat image validation

Jason-2 cycle31衛星軌跡方向為從北到南，剛進入湖泊時星下點離陸地較近，此時的波形如圖 7(a)所示，虛線為跟蹤兩處上升前緣得到的前緣中點位置，可以看出小上升前緣與大上升前緣相隔較近，對應的圓半徑也較小。圖 7(b)中角點圓圈代表了可能引起大上升前緣的地方，可以知道當圓圈還未到達陸地時(也即脈沖未到達陸地)波形沒有上升趨勢，當圓圈正好與陸地相接觸時，波形產生了一個大的上升前緣，故可以認為大的上升前緣由陸地反射造成。當衛星向南移動時，波形中的大、小上升前緣間隔增大(圖 7(c))，對應的圓圈半徑也增大(圖 7(d))，但圓圈始終與陸地相接。當衛星移動至湖泊中心時，波形的大、小上升前緣間隔達到更大(圖 7(e))，對應角點圓圈半徑也達到更大(圖 7(f))，此時圓圈與南北兩岸相接，說明波形中大上升前緣可能受南北兩岸陸地的共同影響。衛星越過湖泊中心繼續向南移動時，波形中大上升前緣主要受南面陸地反射的影響，波形中大、小上升前緣間隔減小(圖 7(g))，對應的圓圈半徑也減小，只與南面陸地邊界相接(圖 7(h))。衛星不斷向南移動的過程中，波形中大、小上升前緣間隔不斷減小(圖 7(i))，圓圈半徑也不斷減小(圖 7(j))，但圓圈邊界始終與陸地相接。

由上述分析可以看出，洪澤湖的測高波形受陸地影響，波形復雜，與海洋波形相差較大。由于陸地反射信號較強，使得波形中陸地反射造成的上升前緣明顯，而水面反射造成的上升前緣不明顯。由于水面反射的信號仍然存在于波形中，設計對應的波形重定算法，即可跟蹤正確的波形前緣中點，得到正確水位。

3 波形重定

現有波形重定算法主要有兩類：①考慮整體波形的重定算法，如重心偏移法(OCOG)^[17]、閾值法^[18]、β參數法^[19]等，由于陸地反射的信號在波形中占主導地位，重定后的波形前緣中點偏離真值較遠；②考慮局部波形的重定算法，如子波形重定算法^[20]，由于不考慮整體波形，只顧及子波形，算法會檢測所有有明顯上升的前緣，故有可能將水面反射形成的上升前緣檢測到，但無法說明提取的上升前緣是否由水面反射產生，只能通過先驗信息加以判斷，如模型值、已有測量值，此 *** 有很大的局限性。綜上分析，無法直接利用現有波形重定算法進行跟蹤，需要根據波形特征設計新的波形重定算法。

受陸地影響的波形有這樣一個特征：波形中先有一個小的上升前緣，隨后有一個大的上升前緣，小的由水面反射，大的受陸地影響較大。設計波形重定算法時，先采用子波形^[20]的思想，提取小的上升前緣，再使用閾值法^[18]確定子波形前緣中點。使用閾值法時，根據文獻[18]的建議，并在試驗中多次嘗試，閾值取50%為更優。由于離岸太近時水面和陸地反射信號夾雜在一起無法分離，故此處將數據范圍縮小一定程度，以有效跟蹤波形前緣中點。利用上述思想計算得到洪澤湖水位，如果水位結果中存在異常值則需要進行異常值剔除。常用異常值剔除 *** 有剔除偏離均值3倍中誤差的異常值^[21]、將水位可視化后人工剔除^[22]、使用四分位數間距剔除^[23]等，本文使用3倍中誤差的 *** 進行異常值剔除。將cycle0-100沿衛星軌跡水位疊加，結果如圖 8所示，可以清晰地看到，水位沿衛星軌跡呈“V”字形分布現象消失。

圖 8 經本文波形重定后的水位疊加圖(cycle0-100)
Fig. 8 The overlay of water level after waveform retracking presented by this paper (cycle0-100)

為進一步驗證波形重定結果的質量，選取了部分洪澤湖水文站水位的觀測數據，該數據為2009年每個月中Jason-2通過湖面時的瞬時水位，每月取一個與衛星軌跡同步的水位觀測值。水文站水位與波形重定后水位對比如圖 9(a)所示。由于兩種水位的高程基準不同，故存在一定的系統偏差，兩者的相關系數為0.99，為強相關。逐月作差得到的水位變化可以消除基準不一致的影響，如圖 9(b)所示，通過兩種手段得到的逐月水位變化的互差的絕對值更大為9 cm，最小為1 cm，標準差為5 cm，說明波形重定算法效果較好，可以用于一般的水位監測。

圖 9 水文站水位與波形重定后水位對比
Fig. 9 Water level contrast of hydrologic station and waveform retracking
 

4 試驗分析與討論

目前波形重定算法繁多，定量分析這些波形重定算法能否處理本文所提取的波形數據，同時將本文提出的波形重定算法與目前常用的波形重定算法進行比較。常見的波形重定算法有β參數法^[19]、重心偏移法(OCOG)^[17]、閾值法(threshold)^[18]、改進后的閾值法^[24, 25]、子波形閾值法^[20]等，Jason-2的SGDR數據集中自帶了Ice、Ocean、MLE3幾種重定算法重定后的測距值^[15]。將幾種常用波形重定算法分別應用于處理洪澤湖區的測高數據，并用IMP(the improvement percentage)^[20]描述波形重定后的改善度，IMP的定義如下

 (3)

式中，σ_raw、σ_re分別為波形重定前、后水位的標準差；IMP越高表明波形重定算法的改善效果越好。對洪澤湖Jason-2cycle0-300的測高數據進行波形重定，各個波形重定算法結果如表 1所示。可以看出，本文設計的波形重定算法表現更好，剔除異常值后水位平均標準差為0.09 m，IMP為95%。閾值法、改進后的閾值法、重心偏移法等優于未重定(tracker)的結果，可見其對本文的波形數據有一定的效果，但改善效果不如本文的波形重定算法好。SGDR數據中自帶的3種波形重定算法Ice、Ocean、MLE3的結果比未重定的結果差，說明其 *** 不適用于本文的波形數據。值得注意的是文獻[20]的波形重定算法結果最不樂觀，雖本文的波形重定思想與其相似，表現卻差異很大，主要原因在于，如果上升前緣中有一個點突然下降，那么這個前緣基本不會被文獻[20]的波形重定算法跟蹤到，而本文的數據恰在較多的上升前緣中有一個下降的點，所以該算法表現差。

表 1 各波形重定算法結果對比Tab. 1 Comparison of different waveform retracking algorithms

5 結論

本文以Jason-2雷達測高數據計算洪澤湖水位為例，探討了湖泊水位沿衛星軌跡呈“V”字形分布的原因。該現象對應的測高波形先有一小上升前緣隨后有一大上升前緣，波形中小上升前緣由水面反射，大上升前緣主要由陸地反射，結合Landsat影像從波形形成機制證實了解釋的正確性。正是由于水面反射形成的上升前緣很小，數據自帶以及常用的波形重定算法很難有效跟蹤正確的波形前緣中點，導致水位計算錯誤。針對此種波形設計了重定算法，結果顯示本文的波形重定算法優于數據自帶以及常用的波形重定算法，水位標準差達0.09 m，可用于一般的水位監測。

本文指出的水位沿衛星軌跡呈“V”字形分布現象不局限于Jason-2和洪澤湖，搭載雷達高度計的衛星獲取的數據如果有相似特征，則可以使用本文的分析 *** 和波形重定算法。特別是衛星離岸較近時，如果陸地反射信號較強，則此區域測高波形中陸地反射信號占主導地位，常用的波形重定算法不適用，如不對其進行分析和波形重定將得到錯誤的水位。

需要說明的是，并不是所有的湖面測高波形都有小前緣并且小前緣是水面反射的，這取決于波形是否受陸地影響以及陸地反射信號的強弱、地形分布等諸多因素，需要綜合分析。本文的波形重定算法表現良好，前提是波形中小上升前緣由水面反射形成，如果波形有小上升前緣但卻是噪聲，則此算法不再適用。
 

【引文格式】田山川, 郝衛峰, 李斐, 等. 顧及陸湖反射差異的衛星測高監測湖泊水位的波形分析與重定[J]. 測繪學報，2018，47(4)：498-507. DOI: 10.11947/j.AGCS.2018.20170348