高賢文 金濤勇,2 黃海蘭,2

1 武漢大學(xué)測(cè)繪學(xué)院，武漢市珞喻路129號(hào)，430079 2 武漢大學(xué)地球空間環(huán)境與大地測(cè)量教育部重點(diǎn)實(shí)驗(yàn)室，武漢市珞喻路129號(hào)，430079

衛(wèi)星雷達(dá)測(cè)高已被廣泛應(yīng)用于海面高及其變化、海洋環(huán)流、海洋重力場(chǎng)、海底地形、冰蓋消融、內(nèi)陸水域變化等研究，在大地測(cè)量學(xué)、地球物理學(xué)、海洋學(xué)、水文學(xué)等領(lǐng)域發(fā)揮著重要作用。脈沖有限雙頻雷達(dá)高度計(jì)是衛(wèi)星雷達(dá)測(cè)高的主要模式，但受限于沿軌跡空間分辨率和測(cè)距精度低，難以滿足近岸、河流、海冰等區(qū)域以及中小尺度海洋現(xiàn)象監(jiān)測(cè)的需求。為此，Ka波段雷達(dá)高度計(jì)、合成孔徑雷達(dá)高度計(jì)、寬刈幅合成孔徑雷達(dá)干涉高度計(jì)以及多星編隊(duì)測(cè)高等模式被相繼提出[1-2]。合成孔徑雷達(dá)(synthetic aperture radar, SAR)測(cè)量模式相對(duì)于傳統(tǒng)脈沖有限雷達(dá)測(cè)量模式具有較高的沿軌跡空間分辨率和測(cè)距精度[3]，該技術(shù)已在Cryosat-2、Sentinel-3A、Sentinel-3B衛(wèi)星中實(shí)現(xiàn)，并將用于Topex/Poseidon系列的后續(xù)衛(wèi)星Jason-CS/Sentinel-6[4]，為今后衛(wèi)星測(cè)高的主要模式。

合成孔徑雷達(dá)高度計(jì)具有較高的沿軌空間分辨率，但在近岸區(qū)域仍會(huì)受復(fù)雜地形海況影響，需要進(jìn)行波形重跟蹤處理。高度計(jì)對(duì)回波波束處理模式不同會(huì)使其回波波形存在較大差別，傳統(tǒng)的脈沖有限雷達(dá)高度計(jì)回波后緣緩慢下降，而合成孔徑雷達(dá)高度計(jì)回波后緣迅速下降，使得傳統(tǒng)的布朗海恩模型無法適用，由此導(dǎo)致OCEAN算法、ALES算法[5]、BP算法[6]等基于布朗海恩模型的波形重跟蹤算法無法用于合成孔徑雷達(dá)高度計(jì)的數(shù)據(jù)處理。部分學(xué)者針對(duì)合成孔徑雷達(dá)高度計(jì)回波問題提出一些基于函數(shù)模型的擬合算法，包括基于延遲多普勒原理[7]的SAMOSA(SAR altimetry mode studies and applications)重跟蹤算法[8]；利用一個(gè)五段函數(shù)模型擬合回波波形的IceSheet重跟蹤算法，每段函數(shù)對(duì)應(yīng)合成孔徑雷達(dá)回波的一個(gè)區(qū)域[9]。此外，傳統(tǒng)數(shù)學(xué)統(tǒng)計(jì)經(jīng)驗(yàn)算法仍適用，包括OCOG算法[10]、閾值法[11]、ICE1算法[12]、多閾值法[13-14]等。

作為未來主要的衛(wèi)星測(cè)高模式，合成孔徑雷達(dá)高度計(jì)的波形重跟蹤算法精度分析對(duì)其應(yīng)用至關(guān)重要。目前基于模型的SAR重跟蹤算法較少，SAMOSA和IceSheet為Sentinel-3A衛(wèi)星官方采用的算法；而傳統(tǒng)經(jīng)驗(yàn)算法中，ICE1算法在獲取內(nèi)陸江河湖泊水面高時(shí)被證實(shí)具有較好的準(zhǔn)確性[15]，多閾值算法在近岸地區(qū)精度較好[13, 16-17]。本文將選取ICE1、多閾值、SAMOSA和IceSheet四種算法，利用全球驗(yàn)潮站觀測(cè)數(shù)據(jù)分析其在近岸20 km范圍內(nèi)的重跟蹤精度，為后續(xù)合成孔徑雷達(dá)高度計(jì)觀測(cè)數(shù)據(jù)處理提供參考。

1 數(shù)據(jù)與方法

1.1 采用數(shù)據(jù)

本文選取Sentinel-3A衛(wèi)星海洋數(shù)據(jù)的Level-2級(jí)非實(shí)時(shí)產(chǎn)品，該數(shù)據(jù)可提供儀器校正及各類地球物理和環(huán)境誤差改正項(xiàng)。Sentinel-3A衛(wèi)星海洋數(shù)據(jù)提供3種下載途徑，其1 a內(nèi)的觀測(cè)數(shù)據(jù)由哥白尼在線數(shù)據(jù)訪問平臺(tái)(https://coda.eumetsat.int/)提供下載；長期歷史數(shù)據(jù)存檔可以在哥白尼數(shù)據(jù)中心(https://www.eumetsat.int/website/home/Data/DataDelivery/EUMETSATDataCentre/index.html)提交訂單后獲取；同時(shí)，考慮到后期算法升級(jí)，哥白尼在線數(shù)據(jù)再處理平臺(tái)(https://codarep.eumetsat.int/)也提供再處理后的數(shù)據(jù)。Sentinel-3A衛(wèi)星軌道高度為814 km，軌道傾角為98.56°，運(yùn)行周期為27 d，每個(gè)周期共770個(gè)弧段，赤道處軌道間距約100 km，其星下點(diǎn)20 Hz采樣間距約為330 m。

本文采用Sentinel-3A衛(wèi)星海洋數(shù)據(jù)的SAR模式20 Hz多視平均回波數(shù)據(jù)，每個(gè)20 Hz波形包含128個(gè)采樣門，每個(gè)距離門對(duì)應(yīng)的時(shí)間延遲為3.125 ns，即46.84 cm的海面高改正值，其典型近岸波形如圖1所示。由于近岸微波輻射計(jì)對(duì)流層濕分量、雙頻電離層改正等依賴星載儀器的改正項(xiàng)精度相對(duì)較差，為便于比較算法的精度，在計(jì)算海面高時(shí)，需進(jìn)行基于ECMWF模型的對(duì)流層干分量和濕分量改正、基于GIM模型的電離層改正、顧及高頻大氣響應(yīng)的逆氣壓改正、基于GOT4.10c模型的海洋潮汐和海洋潮汐負(fù)荷改正、固體潮改正和極潮改正。海況偏差改正一般采用基于有效波高和風(fēng)速的參數(shù)估計(jì)方法進(jìn)行計(jì)算，或?qū)ζ渌姓`差改正后的殘差進(jìn)行非參數(shù)估計(jì)得到；而近岸波形受污染嚴(yán)重，導(dǎo)致有效波高計(jì)算誤差較大，從而使近岸海況偏差改正精度較差?？紤]到選取的驗(yàn)潮站附近測(cè)高弧段僅20 km，在小范圍內(nèi)可認(rèn)為海況偏差改正為常數(shù)，且不同的重跟蹤算法采用的海況偏差改正一致，為避免引入額外誤差影響重跟蹤算法精度，因此未添加海況偏差改正[6]。

圖1 Sentinel-3A衛(wèi)星在近岸地區(qū)典型回波示意圖Fig.1 Typical waveform of Sentinel-3A satellite in coastal area

為驗(yàn)證測(cè)高數(shù)據(jù)的精度，選取美國夏威夷大學(xué)海平面研究中心發(fā)布的1 h采樣。經(jīng)過質(zhì)量控制的驗(yàn)潮站觀測(cè)數(shù)據(jù)(https://uhslc.soest.hawaii.edu/opendap/rqds/contents.html)用于外部驗(yàn)證，并進(jìn)行逆氣壓改正達(dá)到與測(cè)高數(shù)據(jù)一致。首先，以Sentinel-3A衛(wèi)星地面軌跡為參考，篩選出距離地面軌跡觀測(cè)點(diǎn)小于20 km且重疊時(shí)間大于2.5 a的驗(yàn)潮站，并剔除最近星下點(diǎn)軌跡位于內(nèi)陸或島嶼的驗(yàn)潮站；然后，在驗(yàn)潮站附近選取Sentinel-3A衛(wèi)星觀測(cè)弧度上距離海岸線1～20 km、以1 km為采樣間隔的20個(gè)觀測(cè)序列，分別計(jì)算與相鄰驗(yàn)潮站海面高序列的相關(guān)性，剔除平均相關(guān)性小于0.5的驗(yàn)潮站及相應(yīng)的測(cè)高弧段；最后，共選取27個(gè)驗(yàn)潮站及相應(yīng)的Sentinel-3A有效弧段，站點(diǎn)分布如圖2所示，由于有效弧段長度僅20 km，在圖中顯示不明顯，為此給出有效弧段從開始到結(jié)束的整個(gè)部分。此外，對(duì)驗(yàn)潮站數(shù)據(jù)進(jìn)行滑動(dòng)平均濾波去除高頻信號(hào)和噪聲，并利用線性插值方法將驗(yàn)潮站數(shù)據(jù)內(nèi)插到測(cè)高時(shí)間歷元。由于驗(yàn)潮站處的海平面數(shù)據(jù)是相對(duì)于站點(diǎn)基準(zhǔn)面或其他自定義參考面，而衛(wèi)星測(cè)高海面高的參考面為參考橢球面，為將兩者進(jìn)行比較，降低海面高粗差值的影響，分別在兩種時(shí)間序列中移去各自的中值，進(jìn)而消除基準(zhǔn)差，同時(shí)重跟蹤算法的系統(tǒng)偏差也會(huì)被消除。通過中值對(duì)齊，還可消除測(cè)高數(shù)據(jù)中未考慮海況偏差改正的影響。

圖2 選取的27個(gè)驗(yàn)潮站站點(diǎn)以及相應(yīng)的Sentinel-3A有效弧段分布Fig.2 Distribution of 27 selected tide gauges and corresponding Sentinel-3A valid passes

1.2 波形重跟蹤處理與精度評(píng)估方法

在近岸區(qū)域，傳統(tǒng)的多閾值算法利用沿軌開放海域方向最近的重跟蹤海面高作為參考高度，但必須保證從開放海域向沿岸方向搜索，然后將最接近該參考高度的子波形前緣作為有效前緣[16]。而衛(wèi)星測(cè)高在開闊海域沿軌方向20 Hz的海面高噪聲水平較高，如Jason-2數(shù)據(jù)可達(dá)到7.5 cm[18]，該噪聲在近岸區(qū)域更高。因此，如果以鄰近的海面高作為參考，可能會(huì)引入較大不確定性，而通常沿軌星下點(diǎn)的海平面相對(duì)大地水準(zhǔn)面的正高序列相對(duì)穩(wěn)定，故本文采用正高序列的中值作為參考，以抑制沿軌噪聲引入的不確定性。因此，對(duì)于SAR模式的多閾值算法，本文對(duì)其實(shí)施以下步驟進(jìn)行修改：1)對(duì)所有波形采用Hwang等[16]的方法進(jìn)行子波形搜索，并采用閾值重跟蹤得到相應(yīng)海面高；2)扣除所有子波形海面高所對(duì)應(yīng)的EGM2008大地水準(zhǔn)面高得到正高序列；3)當(dāng)波形只有單個(gè)子波形前緣時(shí)，表明其受污染較小，可選擇此類波形對(duì)應(yīng)的正高序列的中值作為參考值；對(duì)于其他多前緣的子波形，選取相應(yīng)正高值距離該參考值最近的子波形前緣所對(duì)應(yīng)的海面高作為該波形的重跟蹤海面高。

針對(duì)4種重跟蹤算法，SAMOSA算法和IceSheet算法的重跟蹤距離采用Sentinel-3A衛(wèi)星數(shù)據(jù)產(chǎn)品提供的結(jié)果，ICE1算法和多閾值算法經(jīng)測(cè)試后選擇最優(yōu)閾值用于獲取重跟蹤距離，然后對(duì)各算法采用相同的地球物理改正。利用4種重跟蹤算法對(duì)近岸20 km內(nèi)的觀測(cè)數(shù)據(jù)進(jìn)行處理得到所有海面高；將海面高以1 km為間隔進(jìn)行平均，得到20個(gè)等間隔的平均海面高時(shí)間序列；利用EGM2008模型校正各點(diǎn)的大地水準(zhǔn)面梯度；最后計(jì)算各點(diǎn)海面高序列與最鄰近驗(yàn)潮站海面高序列的相關(guān)性和差值均方根誤差。此外，本文對(duì)4種算法重跟蹤后保留的有效波形個(gè)數(shù)相對(duì)原始波形的比例進(jìn)行比較，并分析重跟蹤后海面高序列相對(duì)EGM2008模型大地水準(zhǔn)面高差值的標(biāo)準(zhǔn)差。

2 結(jié)果分析

2.1 多閾值和ICE1算法的最優(yōu)閾值選取

由于ICE1算法和多閾值算法涉及閾值的選取，為獲得最佳閾值和最好的重跟蹤效果，對(duì)閾值范圍為0.3～1.0的兩種算法的重跟蹤海面高進(jìn)行分析，并計(jì)算其與所選取的所有驗(yàn)潮站站點(diǎn)海面高序列的相關(guān)性和差值均方根誤差(圖3)，為更好地展示最優(yōu)閾值的效果，未給出范圍內(nèi)所有閾值的重跟蹤結(jié)果。從圖3可以看出，當(dāng)ICE1閾值取0.85、多閾值算法閾值取0.75時(shí)，在近岸1～20 km范圍內(nèi)，兩種海面高序列的相關(guān)性和差值均方根誤差最優(yōu)。由于回波波形存在差異，傳統(tǒng)脈沖有限雷達(dá)高度計(jì)的星下點(diǎn)跟蹤位置大致位于波形前緣中點(diǎn)，故一般閾值較??；而合成孔徑雷達(dá)高度計(jì)的星下點(diǎn)跟蹤位置更接近于峰值，故其閾值較大(圖4)。

圖3 不同閾值情況下驗(yàn)潮站海面高與重跟蹤海面高序列的平均相關(guān)性和差值RMSFig.3 The mean correlation and RMS of difference between tide gauge sea levels and retracked sea levels under different thresholds

圖4 有效波高4 m時(shí)脈沖有限雷達(dá)回波波形和合成孔徑雷達(dá)多視回波波形Fig.4 The waveforms of pulse limited radar altimeter and synthetic aperture radar altimeter under 4 m significant wave height

2.2 重跟蹤后保留的有效波形比例

對(duì)4種算法重跟蹤后的20個(gè)海面高時(shí)間序列進(jìn)行分析，在每1 km處統(tǒng)計(jì)原始觀測(cè)的回波波形個(gè)數(shù)，然后計(jì)算不同距離處重跟蹤成功波形個(gè)數(shù)與原始波形個(gè)數(shù)之比，以分析不同離岸距離時(shí)各算法的重跟蹤成功率[17]，將明顯偏離正常區(qū)間的重跟蹤結(jié)果認(rèn)定為重跟蹤失敗，本文采用3倍標(biāo)準(zhǔn)差準(zhǔn)則進(jìn)行判定。表1為各算法統(tǒng)計(jì)結(jié)果，由表可見，除IceSheet算法在1 km處外，4種算法在所有距離段的重跟蹤率均大于50%；多閾值算法在近岸9 km范圍內(nèi)成功率最高，且在3 km 時(shí)已超過90%。在近岸9 km范圍內(nèi)，IceSheet算法的成功率始終低于其他算法，ICE1算法的成功率略高于SAMOSA算法；在近岸9～20 km范圍內(nèi)，SAMOSA算法成功率最高，表明在開闊海域該算法成功率較高。

表1 不同離岸距離情況下4種重跟蹤算法的有效波形比例Tab.1 The rates of effective waveform of four retrackers in different distances to the coast

2.3 相關(guān)性與差值均方根誤差分析

分別計(jì)算4種算法重跟蹤后的20個(gè)海面高時(shí)間序列與最鄰近驗(yàn)潮站海面高時(shí)間序列的相關(guān)性和差值均方根誤差，并統(tǒng)計(jì)不同距離情況下4種算法的平均相關(guān)性和平均差值均方根誤差(圖5)。由圖可見，在近岸20 km范圍內(nèi)，SAMOSA算法相關(guān)性最低且均方根誤差最大；在6 km范圍內(nèi)，多閾值算法平均相關(guān)性最高且差值均方根誤差最?。辉? km范圍外，多閾值和ICE1算法基本相當(dāng)；IceSheet算法始終處于中等水平。

圖5 不同距離情況下4種算法重跟蹤海面高與驗(yàn)潮站海面高時(shí)間序列的平均相關(guān)性和差值RMSFig.5 The mean correlation and RMS of difference between retracked sea levels and tide gauge sea levels of four retrackers in different distances to the coast

2.4 相對(duì)大地水準(zhǔn)面模型的差值穩(wěn)定性分析

由于平均海面高模型在近岸精度較差，且其主要采用傳統(tǒng)衛(wèi)星測(cè)高數(shù)據(jù)，而大地水準(zhǔn)面模型在近岸具有較好的穩(wěn)定性，故可通過計(jì)算重跟蹤海面高序列與大地水準(zhǔn)面的差值標(biāo)準(zhǔn)差來衡量算法自身的變異性[19]。將近岸20 km區(qū)域劃分為若干個(gè)距離區(qū)間，分別統(tǒng)計(jì)每個(gè)驗(yàn)潮站在各距離區(qū)間的重跟蹤海面高序列與EGM2008大地水準(zhǔn)面模型的差值標(biāo)準(zhǔn)差，最后取區(qū)間內(nèi)所有驗(yàn)潮站的平均標(biāo)準(zhǔn)差作為該區(qū)間的內(nèi)部可靠性評(píng)價(jià)。為保證足夠的沿軌跡采樣數(shù)量，以5 km為間隔進(jìn)行區(qū)間計(jì)算，結(jié)果見表2。從表中可以看出，在近岸0～5 km范圍內(nèi)，IceSheet算法略優(yōu)于其他算法；在近岸5～20 km范圍內(nèi)，SAMOSA算法顯著優(yōu)于其他算法，IceSheet、ICE1和多閾值算法效果相當(dāng)，表明SAMOSA算法具有較高的穩(wěn)定性(特別是在開闊海域)，而在近岸近距離區(qū)域回波波形偏離標(biāo)準(zhǔn)海洋波形，導(dǎo)致基于回波模型的重跟蹤算法出現(xiàn)偏差，與驗(yàn)潮站實(shí)測(cè)海面高的相關(guān)性和均方根誤差偏差較大。此外，在離岸距離15～20 km范圍內(nèi)的均值略大于5～15 km區(qū)間的均值，其原因?yàn)閭€(gè)別弧段的大地水準(zhǔn)面起伏異常，但其影響對(duì)各重跟蹤算法基本一致，因此對(duì)各算法結(jié)果的比較不會(huì)造成影響。

表2 不同離岸距離下4種算法重跟蹤海面高與EGM2008大地水準(zhǔn)面高差異的沿軌標(biāo)準(zhǔn)差Tab.2 The along-track STD of the difference between retracked sea surface height and EGM2008 geoid heights of four retrackers in different distances to the coast

3 結(jié) 語

基于Sentinel-3A衛(wèi)星合成孔徑雷達(dá)高度計(jì)觀測(cè)數(shù)據(jù)，利用全球范圍內(nèi)27個(gè)驗(yàn)潮站海面高數(shù)據(jù)，采用4種波形重跟蹤算法對(duì)近岸20 km范圍內(nèi)鄰近測(cè)高波形數(shù)據(jù)進(jìn)行重跟蹤，對(duì)各算法的精度進(jìn)行分析。結(jié)果表明，ICE1和多閾值算法分別選取閾值0.85和0.75時(shí)精度最優(yōu)，在此條件下多閾值算法在近岸6 km范圍內(nèi)可保留最多的有效波形個(gè)數(shù)，且與驗(yàn)潮站海面高序列的相關(guān)性最高和均方根誤差最??；IceSheet算法具有較低的有效波形比例，但在近岸5 km范圍內(nèi)具有最好的穩(wěn)定性；SAMOSA算法在近岸20 km范圍內(nèi)的重跟蹤海面高與驗(yàn)潮站海面高相差最大，但在離岸距離大于5 km后算法穩(wěn)定性最高，適用于開闊海域?？傮w而言，在近岸近距離海域，基于數(shù)學(xué)統(tǒng)計(jì)的多閾值重跟蹤算法精度較好；在開闊海域，基于物理模型的SAMOSA重跟蹤算法精度較好。

雖然現(xiàn)在已積累大量合成孔徑雷達(dá)高度計(jì)實(shí)測(cè)數(shù)據(jù)，但相對(duì)脈沖有限雷達(dá)高度計(jì)，在波形數(shù)據(jù)處理及應(yīng)用上的研究較少。新的高度計(jì)體系在理論上具有更高的沿軌分辨率和測(cè)距精度，但仍需進(jìn)行更深入的驗(yàn)證和分析，以推進(jìn)合成孔徑雷達(dá)測(cè)高的應(yīng)用，彌補(bǔ)脈沖有效雷達(dá)高度計(jì)的不足。