周曉光，苗洪利，王云海，范陳清，崔廷偉，張 杰

1.中國海洋大學 信息科學與工程學院，山東 青島 266100；2.國家海洋局 第一海洋研究所，山東 青島 266061

1 引 言

衛(wèi)星繞地球運行一周的軌道地面軌跡可分為2個弧段（pass），分別為上升弧段和下降弧段，在一個軌道重復周期內，若干上升弧段和下降弧段所形成的交點稱為交叉點。交叉點是由于地球自轉所形成［1］。在交叉點處，上升弧段和下降弧段具有相同的經(jīng)緯度（即相同的地理位置）。

文獻［2］提出將衛(wèi)星的長期攝動和地球自轉結合，基于地球自轉模型的天文學計算得到了理論交叉點。由于重復周期的實際軌道并不嚴格重合，導致了天文學計算得到的理論交叉點位置與實際位置相差過大。為了解決此問題，文獻［3—4］提出了由理論交叉點位置出發(fā)迭代求解交叉點實際位置的方法。但是在將地球視為形狀規(guī)則、密度均勻的正球體的二體系統(tǒng)中，雖然衛(wèi)星運動軌跡大致可以由軌道傾角、升交點赤徑、軌道長半軸、軌道偏心率、近地點角距和平近點角6個參量確定，但由于高度計采用大地經(jīng)緯度，此時地球被視為密度不等的旋轉橢球體，而且衛(wèi)星的軌道受到地球引力、大氣阻力、太陽光壓力以及月球、太陽、行星等因素影響而產(chǎn)生攝動，所以很難通過數(shù)學計算得到交叉點的精確位置。為了適應區(qū)域平差的要求，一些學者將全球劃分為六大區(qū)域，分區(qū)域對其上升弧段和下降弧段進行二次項擬合，聯(lián)立求解交叉點位置［5］，如文獻［6］在研究多星聯(lián)合確定中國近海重力異常時，在中國近海及其鄰海區(qū)域內的T／P、ERS2和Geosat 3顆衛(wèi)星地面軌跡進行二次多項式擬合計算交叉點，但這種方法只能用于區(qū)域平差，在全球范圍內高度計的實際軌道與二次多項式曲線相差較大，如果對全段直接擬合所得出的二次多項式進行聯(lián)立，其得到的交叉點概略位置與精確位置會相差過大，這將導致從概略位置到精確位置的逼近效率變低；由于高緯度地區(qū)擬合曲線和實際地面軌跡偏離較大，使二次多項式聯(lián)立無法得到交叉點概略位置，導致高緯度區(qū)域交叉點數(shù)大量缺失。

交叉點精確位置對應的海表面信息是由上升弧段和下降弧段周邊實測點線性插值得到，因此交叉點精確位置確定直接影響交叉點測量不符值的精度。而高度計時標偏差計算［7］、衛(wèi)星數(shù)據(jù)定軌［8］、衛(wèi)星軌道誤差修正［9－11］、大地水準面模型的建立［12］、利用測高數(shù)據(jù)反演重力場［13］、海況偏差的計算［14－16］、潮汐的調和分析［17－18］以及 多星數(shù) 據(jù)聯(lián)合方面的研究［19－20］都是以交叉點數(shù)據(jù)為基礎的。天文學計算方法過于復雜，計算效率低。而適用于區(qū)域平差的擬合方法又無法應用到全球范圍。因此找到一種準確高效、適合于全球范圍的交叉點確定方法十分重要。

2 方法與數(shù)據(jù)處理

2.1 最優(yōu)分段

因為全球范圍內衛(wèi)星地面軌跡與二次多項式曲線差距較大，所以通過對全段軌跡的二次擬合得到交叉點概略位置的可行性不大。本文根據(jù)衛(wèi)星地面軌跡的特征，對其進行合理的分段，然后對分段后的軌跡進行二次擬合，并對升、降軌跡擬合二次多項式聯(lián)立得到所需交叉點概略位置。

由于衛(wèi)星的地面軌跡與其軌道傾角關系密切，因此不同衛(wèi)星地面軌跡特征差距較為明顯，本文以軌道傾角為66°的Jason－1衛(wèi)星高度計為例進行分析，圖1為Jason－1衛(wèi)星地面軌跡的示意圖。如圖可知衛(wèi)星地面軌跡為類正弦曲線，因此本文根據(jù)其曲線的彎曲特性，采用不等間隔的分段方法設計了11種分段類型，如表1所示。按照各種類型進行分段擬合并進行擬合效果分析，最終確定最優(yōu)分段類型。

提取Jason－1衛(wèi)星高度計第126個重復周期的第21弧段（上升）和64弧段（下降）的GDR（geophysical data record）2級數(shù)據(jù)產(chǎn)品中的經(jīng)緯度信息，按照表1中所列的分段方法進行分段，對每段經(jīng)度λi和緯度φi的實測值按照式（1）進行擬合，確定系數(shù)A、B、C

圖1 Jason－1衛(wèi)星地面軌跡示意圖Fig.1 The view of ground track of Jason－1altimeter

表1 各類型的分段方式Tab.1 Description of the region of latitude in different types

將GDR中的經(jīng)度值λi代入擬合方程（1）得到緯度的計算值φi，通過式（2）與實測值φi進行標準偏差（S）計算

式中，S為標準偏差；N表示測量點的數(shù)目。標準偏差表示了擬合值和實測值的偏離程度，其值越小則說明擬合值和實測值越貼近，擬合效果越好。在實際應用中，表2為各類型分段擬合后整個上升弧段或下降弧段的標準偏差。經(jīng)統(tǒng)計Jason－1衛(wèi)星高度計地面軌跡中兩個相鄰點緯度的差值均值為0.039 4，因此分段擬合之后整個弧段緯度的擬合標準偏差在小于0.039 4的前提下，擬合偏差最小的分段類型即認為是最優(yōu)分段類型。

分析表2可以發(fā)現(xiàn)，pass 21和pass 64的標準偏差均以Type－12為界，在此之前呈明顯下降趨勢而在此之后略有上升，值得注意的是Type－4的標準偏差與相鄰的Type－3和Type－5相比明顯要小，這是因為Jason－1高度計的軌道地面軌跡在－45°和45°緯度附近明顯曲折并不適宜進行二次多項式擬合，Type－4恰好以－45°和45°為分段點，而Type－3和Type－5則包含了這兩個拐點，導致擬合效果不如 Type－4好。Type－18和 Type－36盡管分段更細，不但運算效率會下降，而且平均分段忽略了衛(wèi)星地面軌跡的彎曲特征導致擬合效果下降。因此Type－12所述的12分段最優(yōu)。由于不同傾角衛(wèi)星地面軌跡會發(fā)生變化，分段節(jié)點需要根據(jù)軌跡曲線特征作適當調整。但本文分段設計方案及最優(yōu)分段選取方法對任何傾角的衛(wèi)星依然適用。

表2 各分段類型不同弧段的擬合標準偏差Tab.2 Standard deviation of different pass under different types

2.2 交叉點概略位置確定

按照Type－12的分段類型將升降弧段分別劃分為12個緯度區(qū)域，從－90°～－65°緯度區(qū)域開始，對每個緯度區(qū)域內的所有上升弧段分別與所有的下降弧段利用如下條件判斷二者是否存在交叉點：

（1）上升弧段在該緯度區(qū)域內第一點的經(jīng)度小于下降弧段在該緯度區(qū)域最后一點的經(jīng)度；

（2）上升弧段在該緯度區(qū)域內最后一點的經(jīng)度大于下降弧段在該緯度區(qū)域內第一點的經(jīng)度。

如果存在交叉點則在該緯度區(qū)域內對兩弧段利用GDR中逐點的經(jīng)緯度實測值分別進行擬合，得到二者的擬合二次多項式，聯(lián)立求解確定交叉點概略位置。

2.3 交叉點精確位置確定

確定交叉點的精確位置是為了獲取該位置對應的海況信息（如海表面高度、有效波高、風速等）。而恰好在交叉點處有測量值的概率很小。一般都是找到精確位置升、降弧段兩側各兩個數(shù)據(jù)點Qa1、Qa2和Qd1、Qd2（簡稱插值點），角標a表示升弧段，d表示降弧段，利用4個點的經(jīng)緯度值進行插值得到交叉點的精確位置。同時交叉點升弧段和降弧段的海況信息分別由各自附近兩個插值點的測量值線性插值得到。精確位置確定方法如下：

在升、降弧段上分別找到與概略位置相同緯度位置上下相鄰的兩個數(shù)據(jù)點，數(shù)據(jù)點位置分別表示為Pa1、Pa2、Pd1、Pd2（簡稱初始點），如果點Pa1和Pa2組成的線段與Pd1和Pd2組成的線段存在交點，則初始點Pa1、Pa2、Pd1、Pd2就是所求的插值點Qa1、Qa2、Qd1、Qd2，本文中此情況簡稱case 1，如圖2（a）所示。如果點Pa1和Pa2組成的線段與Pd1和Pd2組成的線段不存在交點，則必須擴展線段的長度，使其產(chǎn)生交點，以該交點為新的概略位置得到新的4個初始點，再次判斷新的4個初始點是否就是Qa1、Qa2、Qd1、Qd2，如果不是則重復上述步驟，直到得到滿足要求的4個插值點，此情況簡稱為case2，如果上述步驟重復n次則簡稱為case 2－n，如圖2（b）所示。

圖2 pass 21和pass 64以及pass 1和pass12所構成交叉點精確位置的確定Fig.2 Find the precise location of crossover formed by pass 21，pass64and pass1，pass12

3 結果分析

以Jason－1高度計第126重復周期的數(shù)據(jù)為例，針對表1中分段類型 Type－1、Type－9、Type－12和Type－36，找到的交叉點概略點的總數(shù)、符合case1和case2－n概略點的數(shù)量及占概率點總數(shù)的百分比，見表3。

表3 不同分段類型概略點的數(shù)量及占總數(shù)的百分比Tab.3 Number and percent of cursory location of crossover in different condition and different experiments

4個不同分段類型的結果表明：隨著分段區(qū)間的增加能夠得到的交叉點概略位置點的數(shù)量在增加，但是分段超過12段之后所增加的概略點數(shù)量已經(jīng)極少，這是由于Type－12的分段擬合效果已經(jīng)很好，基本接近交叉點概略位置數(shù)量極限。隨著分段區(qū)間的增加符合case 1的概略點所占比重逐漸增加，在Type－12和Type－36中大部分的概略點已經(jīng)不需要進行case2中的步驟即可找到精確位置周圍4點。由此說明，分段越細，找到精確位置的效率越高。綜合數(shù)量和效率兩方面考慮，12分段已是最優(yōu)分段。

使用分段類型Type－1、Type－12分別確定Jason－1衛(wèi)星高度計第126個重復周期的第113個弧段上的交叉點位置，結果見圖3。由于數(shù)據(jù)點過于密集，將一條上升軌跡分為3個緯度區(qū)間分別作圖。

分析圖3，可以發(fā)現(xiàn)未分段Type－1確定的交叉點基本都集中在南北緯40°之間的中低緯度地區(qū)，而且與Type－12確定的交叉點完全重合，說明使用Type－12分段類型與Type－1分段類型確定交叉點有著相同的精度，見圖3（b）。但在高緯度Type－1無法找到交叉點，而Type－12則能確定更多的有效交叉點，如圖3（a）、圖3（c）所示。

圖3 使用分段類型Type－1和Type－12確定的交叉點位置Fig.3 Location of crossover calculated using Type－1and Type－12

4 結 論

為了能夠高效、精確確定交叉點的位置，本文提出了分段擬合衛(wèi)星地面軌跡的方法，設計了11種分段類型，并應用于Jason－1高度計交叉點位置確定過程中。以整個弧段緯度的擬合標準偏差作為確定最優(yōu)分段類型的標準，其中擬合偏差最小的分段類型被認為是最優(yōu)分段類型。最終結果表明：不同的分段方法擬合的效果不同，能找到交叉點概略位置的數(shù)量隨著分段區(qū)間的增加而增加，但是超過12段之后概略位置點數(shù)量趨于穩(wěn)定。隨著分段區(qū)間數(shù)的增加，概略位置與精確位置更加接近，這樣能夠更快地確定交叉點的精確位置。綜合交叉點數(shù)量和效率兩方面考慮，采用12段的分段方法確定交叉點效果最優(yōu)。由于不同傾角衛(wèi)星地面軌跡會發(fā)生變化，分段節(jié)點需要根據(jù)軌跡曲線特征適當調整。本文分段設計方案及最優(yōu)分段選取方法對任何傾角的衛(wèi)星依然適用。

［1］ ZHAI Guojun，HUANG Motao，XIE Xijun，et al.The Theory and Method of Satellite Altimeter Data Processing［M］.Beijing：Surveying and Mapping Press，2000：106－107.（翟國君，黃謨濤，謝錫君，等.衛(wèi)星測高數(shù)據(jù)處理的理論與方法［M］.北京：測繪出版社，2000：106－107.）

［2］ HAGAR H.Engineering Memorandum 315－29：Computation of Circular Orbit Ground Trace Intersections［R］.［S.l.］：Jet Propulsion Laboratory，1977.

［3］ ROWLANDS D.The Adjustment of SEASAT Altimeter Data on a Global Basis for Geoid and Sea Surface Height Determination［R］.Columbus：Ohio State University，1981.

［4］ SHUM C K.Altimeter Methods in Satellite Geodesy［R］.Austin：University of Texas，1982.

［5］ WANG Guangyun，WANG Haiying，XU Guochang.The Theory of Satellite Altimeter［M］.Beijing：Science Press，1995：207－208.（王廣運，王海瑛，許國昌.衛(wèi)星測高原理［M］.北京：科學出版社，1995：207－208.）

［6］ LI Jiancheng，NING Jinsheng，CHEN Junyong，et al.Determination of Gravity Anomalies over the South China Sea by Combination of TOPEX／Poseidon，ERS2and Geosat Altimeter Data［J］.Acta Geodaetica et Cartographic Sinica，2001，30（3）：197－202.（李建成，寧津生，陳俊勇，等.聯(lián)合TOPEX／Poseidon，ERS2和Geosat衛(wèi)星測高資料確定中國 近 海 重 力 異 常 ［J］.測 繪 學 報，2001，30（3）：197－202.）

［7］ SCHUTZ B E，TAPLEY B D，SHUM C K.Evaluation of the SEASAT Altimeter Time Tag Bias［J］.Journal of Geophysical Research，1982，87（C5）：3239－3245.

［8］ ZHANG Rongzhi，ZHOU Fengqi，XIONG Jing.Using Crossover Data from Satellite Altimeter for Orbit Determination［J］.Acta Astronomica Sinica，2010，3（51）：317－323.（張榮之，周鳳歧，熊菁.星載測高儀交叉點數(shù)據(jù)定軌研究［J］.天文學報，2010，51（3）：317－323.）

［9］ CLOUTIER J R.A Technique for Reducing Low－frequency，Time－dependent Errors Present in Network－type Surveys［J］.Journal of Geophysical Research，1983，88（B1）：659－663.

［10］ KOZEL B J.Dual－satellite Altimeter Crossover Measurements for Precise Orbit Determination［D］.Austin：University of Texas，1995.

［11］ WANG Haiying.Satellite Altimeter Data Processing and Its Applications in China Seas and Vicinity［J］.Acta Geodaetica et Cartographic Sinica，2001，32（2）：185.（王海瑛.中國近海衛(wèi)星測高數(shù)據(jù)處理與應用研究［J］.測繪學報，2001，32（2）：185.）

［12］ KLOKOCNIK J，WABNER C A.A Test of GEM－T2from GEOSAT Crossovers Using Latitude Lumped Coefficients［J］.Journal of Geodesy，1994，68（2）：100－108.

［13］ HUANG Motao，WANG Rui，ZHAI Guojun，et al.Integrated Data Processing for Multi－satellite Missions and Recovery of Marine Gravity Field［J］.Geomatics and Information Science of Wuhan University，2007，11（32）：988－993.（黃謨濤，王瑞，翟國君，等.多代衛(wèi)星測高數(shù)據(jù)聯(lián)合平差及重力場反演［J］.武漢大學學報：信息科學版，2007，11（32）：988－993.）

［14］ GASPAR P，OGOR F，LE TRAON P Y，et al.Estimating the Sea State Bias of the TOPEX and POSEIDON Altimeters from Crossover Differences［J］.Journal of Geophysical Research，1994，99（C12）：24981－24994.

［15］ FLORENS J P，GASPAR P.Estimation of the Sea State Bias in Radar Altimeter Measurements of Sea Level：Results from a New Non Parametric Method［J］.Journal of Geophysical Research，1998，103（15）：803－814.

［16］ GASPAR P，LABROUE S，OGOR F，et al.Improving Nonparametric Estimates of the Sea State Bias in Radar Altimeter Measurements of Sea Level［J］.Journal of Atmospheric ＆ Oceanic Technology，2002，19（10）：1690－1707.

［17］ BAO Jingyang，CHAO Dingbo，LI Jiancheng.Tidal Harmonic Analysis Near Crossovers of TOPEX／POSEIDON Ground Track in South China Sea［J］.Acta Geodaetica et Cartographic Sinica，2001，29（1）：17－23.（暴景陽，晁定波，李建成.南中國海TOPEX／POSEIDON軌跡交叉點測高數(shù)據(jù)的潮汐調和分析［J］.測繪學報，2001，29（1）：17－23.）

［18］ MAO Qingwen，SHI Ping，QI Yiquan.Tide Separation from the Altimetry Data Using Harmonic Analysis Method［J］.The Ocean Engineering，2002，20（1）：41－45.（毛慶文，施平，齊義泉.運用調和分析方法分離衛(wèi)星高度計資料中的潮汐信息［J］.海洋工程，2002，20（1）：41－45.）

［19］ HUANG Motao，ZHAI Guojun，OUYANG Yongzhong，et al.Recovery of Marine Gravity Field Using Integrated Data from Multi－satellite Missions［J］.Science of Surveying and Mapping，2006，31（6）：37－39.（黃謨濤，翟國君，歐陽永忠，等.利用多代衛(wèi)星測高數(shù)據(jù)反演海洋重力場［J］.測繪科學，2006，6（31）：37－39.）

［20］ HU Jianguo，ZHANG Chuanyin，CHANG Xiaotao.Methods of Multi－satellite Altimetry Data Processing and Its Application［J］.Bulletin of Surveying and Mapping，2004（1）：1－4.（胡建國，章傳銀，常曉濤.近海多衛(wèi)星測高數(shù)據(jù)聯(lián)合處理的方法及應用［J］.測繪通報，2004（1）：1－4.）