賈新強,高關根,張亞崇
(1 海軍裝備部,西安 710065;2 西安飛行自動控制研究所,西安 710065)
機載地形輔助導航作為強電磁對抗環(huán)境下一種有效的自主航空導航手段,日益受到國際導航界的關注和重視?,F(xiàn)階段,機載地形輔助導航除了能完成傳統(tǒng)的定位功能之外,還具有地形跟隨/地形回避、合理路由選擇、近地告警、目標截獲和精確武器投放、導航系統(tǒng)水平通道/垂直通道完好性監(jiān)控等極具軍事應用價值的多項功能[1]。如英國BAE公司開發(fā)的地形剖面匹配系統(tǒng)TERPROM?(terrain profile matching,TERPROM)是目前世界上應用最廣的一種地形輔助導航系統(tǒng),可提供地形匹配、地形防撞、先進地形回避提示、障礙物/電線告警和提示、空地測距、地形感知顯示等功能,能夠有效增強機組的態(tài)勢感知能力、減輕機組的工作負擔及防止可控飛行撞地,其中地形匹配是其核心功能[2]。
地形匹配技術產生于20世紀40年代末50年代初,美、英等發(fā)達國家相繼提出了多種地形匹配方案,并陸續(xù)付諸實施。典型應用方案包括地形輪廓匹配(terrain contour matching, TERCOM)、桑迪亞慣性地形輔助導航(Sandia terrain aided navigation, SITAN)、TERPROM?和精密地形輔助導航系統(tǒng)(precision terrain aided navigation system, PTAN)等。表1是對國外典型地形匹配系統(tǒng)應用情況的統(tǒng)計分析。
表1 國外典型地形匹配系統(tǒng)對比
隨著全球導航衛(wèi)星系統(tǒng)(global navigation satellite system, GNSS)的廣泛應用,航空領域地形匹配系統(tǒng)的作用在逐漸減弱,例如在戰(zhàn)斧Block IV上已用全球定位系統(tǒng)(global positioning system,GPS)替代了地形匹配系統(tǒng),地形匹配系統(tǒng)僅作為備份系統(tǒng)。然而,由于GNSS信號的脆弱性,特別是水下載體和山區(qū)機載導航應用時,GNSS信號存在拒止、有意或無意的干擾或欺騙,地形匹配依然是可靠的自主定位導航信息源。針對水下載體自主導航定位需求,已經形成了地形輪廓匹配、多波束匹配和三維匹配等水下地形匹配系統(tǒng)和相關的理論方法[3-5];針對機載應用,當前地形匹配的主要工作聚焦在如何適應機載大范圍機動飛行和系統(tǒng)可靠性方面[6]。隨著機載激光掃描系統(tǒng)的發(fā)展,通過激光點云與獲取的影像進行匹配,可生成高精度數(shù)字高程數(shù)據,可有效提升地形匹配系統(tǒng)的性能[7]。此外,在GPS拒止環(huán)境下,結合地形特征的視覺導航和地磁異常輔助導航等也可為載機提供有效的定位導航信息[8]。
由于機載應用環(huán)境的復雜性、不確定性以及機載任務自身的復雜性,有必要對機載地形匹配的關鍵技術特征及典型應用場景進行較為全面的辨析,以便于充分理解機載地形匹配系統(tǒng)的技術應用特征及其適用場景,從而為該技術后續(xù)的發(fā)展提供必要的參考和借鑒。
地形匹配,是源自合成孔徑雷達之父Carl Wiley提出的地形匹配導航概念(map matching),在加入古德依爾飛機公司(Goodyear aircraft Co.)后將概念變?yōu)楝F(xiàn)實,并成為該公司的拳頭產品——自動地形識別和導航地圖匹配雷達系統(tǒng)(automatic terrain recognition and navigation map-matching radar system,ATRAN)[9]。
利用地形特征進行飛機導航是人們所熟知的古老導航技術。自從飛機出現(xiàn),飛行員就通過目視地形、地物進行導航。然而,現(xiàn)代機載地形匹配與傳統(tǒng)的地形導航技術截然不同,它是利用飛機正下方地形高程的起伏特性進行定位的一種系統(tǒng),即通過機載高度測量設備(如無線電高度表/激光雷達、大氣數(shù)據系統(tǒng)等)測量出飛機下方地形的一序列高程數(shù)據(測量地形剖面),利用模式識別技術將測量的地形剖面與預先存儲的基準地形剖面進行相關處理,找出最佳的匹配地形剖面并計算得到飛機的位置。本質上,機載地形匹配系統(tǒng)是通過綜合處理地形特征傳感器(如大氣數(shù)據系統(tǒng)和無線電高度表的組合)的測量數(shù)據和載機存儲的數(shù)字地形高程數(shù)據來估計飛機的精確位置,然后再用這個精確位置對機載慣導進行修正,如此不斷循環(huán),就能連續(xù)獲得飛機的精確位置。圖1是機械地形,匹配系統(tǒng)工作過程的簡要原理示意圖。
圖1 機載地形匹配系統(tǒng)示意圖
一般情況下,機載地形匹配系統(tǒng)在地形起伏特征明顯的地區(qū)定位精度很高,而在地形平坦地區(qū)或水面上,僅能提供高度信息。
機載地形匹配系統(tǒng)是在實測地形剖面數(shù)據和基準地形高程數(shù)據兩者均存在一定程度的噪聲條件下尋求最佳導航定位結果的一種技術,系統(tǒng)組成如圖2所示,主要包括[10]:
1) 地形特征傳感設備(無線電高度表/激光雷達、大氣數(shù)據系統(tǒng)等);
2) 主導航設備(慣導或其他推算導航設備);
3) 匹配處理器(地形高程數(shù)據庫、捕獲濾波器、跟蹤濾波器等)。
圖2是從系統(tǒng)軟、硬件及工作邏輯兩方面給出的機載地形匹配系統(tǒng)組成示意圖。
圖2 機載地形匹配系統(tǒng)組成示意圖
地形特征傳感設備包含無線電高度表/激光雷達、大氣數(shù)據系統(tǒng)等。其中大氣數(shù)據系統(tǒng)通過對氣流動壓、靜壓和溫度等的測量,經計算與修正可自主獲得氣壓高度和真空速等參數(shù)[11]。由于機載地形匹配系統(tǒng)不敏感因大氣數(shù)據系統(tǒng)工作機理、季節(jié)、氣候等原因造成的氣壓高度的常值偏差,所以氣壓高度的常值偏差不影響正確匹配概率和定位精度,只有其隨機誤差會對正確匹配概率和定位精度有一定的影響。
無線電高度表用于自主測量飛機與地面的相對距離,其波束示意圖如圖3所示。目前無線電高度表的硬件測量精度很高(一般優(yōu)于1 m),但是飛機正下方離地高度的測量精度只有高度值的1%~3%,這主要是由于無線電高度表有一定的波束角,飛行高度越高,波束覆蓋范圍越大,其波束示意圖如圖3所示。覆蓋范圍內的地形變化越復雜,因此其測高誤差越大[12]。從原理上分析,無線電高度表的測量值對地形高度有平滑效果。因此,在機載地形匹配系統(tǒng)實現(xiàn)中,為了尋求定位精度與數(shù)據存儲量之間的最佳折衷,數(shù)字地形高程數(shù)據庫的分辨率應和無線電高度表的測量精度相匹配。目前,國外大多數(shù)軍用地形匹配系統(tǒng)使用的數(shù)字地形高程數(shù)據庫的網格間距在100 m左右,能滿足地形匹配要求。
圖3 無線電高度表波束示意圖
激光雷達波束覆蓋范圍很小,具有更精確的相對高度測量值。但其測量距離較短、測量精度受天氣的影響大,且易跟蹤地面樹木、建筑物等獨物體,上述情況均會對匹配性能產生不利影響。
數(shù)字地形高程數(shù)據庫是通過對地形高程的離散采樣并量化后得到的,其性能一般由地圖大小、水平和垂直參考坐標系、格網尺寸、圓誤差(circular error,CEP)和線誤差(linear error,LEP)等指標決定[13],圖4是美國軍用1~5級數(shù)字地形高程數(shù)據庫的具體情況。按照國內制圖規(guī)范或要求,1∶50 000數(shù)字高程模型對鄰近高程控制點的高程中誤差一般為平地≤3.0 m,丘陵≤4.5 m,山地≤6.0 m,高山地≤10.0 m。
圖4 美國軍用1~5級數(shù)字地形高程數(shù)據庫示意
由上述分析可見,通過氣壓高度和無線電高度可自主測量地形高程,極少受外界干擾;地表的地形起伏基本不受時間的影響,人為的地形地貌的改變也不常發(fā)生。因此,機載地形匹配是一種自主、抗干擾和全天候的理想定位系統(tǒng)。
主導航設備指慣導或其他推算導航設備。慣導用于提供一個粗略的初始位置和誤差范圍,機載地形匹配系統(tǒng)據此確定在地形高程數(shù)據庫中的搜索范圍。通過限制搜索范圍一是降低計算量,二是減少多重匹配情況。此外,僅通過一個位置的高程測量值尚不足以獲得明確的定位結果,需要通過慣導的速度信息將多個位置的測量參數(shù)結合形成一個定位特征(剖面),與數(shù)字地形高程數(shù)據庫中的特征進行匹配,如圖5所示;同時,慣導的速度用于控制測量數(shù)據的采樣,慣導的位置和速度還用于描述地形剖面的走向。在設計匹配定位準則時,需利用慣導短時精度高、輸出連續(xù)等特點,對最優(yōu)匹配濾波器(與先驗垂直通道模型匹配程度最好的子濾波器)是否就是所要尋找的正確濾波器(水平位置估值最接近真實位置的子濾波器)進行判決。因此,慣導性能對于機載地形匹配系統(tǒng)的精度和可靠性有著重要的影響。
圖5 采用測量剖面的模式匹配示意
地形匹配系統(tǒng)的核心基本上是一種算法功能,這個功能可在飛機上已有的任何計算機上運行。而要達到所預期的定位精度和正確匹配概率等技術指標要求,主要取決于適度精度的高度傳感器、適度精度的數(shù)字高程模型數(shù)據以及所采用的地形匹配算法。
目前,典型的地形匹配算法包括序列相關算法、單點迭代算法以及組合匹配算法。
序列相關算法是當采集到的高程序列長度達到設定的數(shù)量之后,就進行一次匹配運算,即僅依靠地形高程剖面進行相關處理。匹配完成之后,將修正信息提供給慣導。這類算法主要包括TERCOM算法和等值線匹配(iterative closest contour point,ICCP)算法。
由于需要對地形高程序列進行相關分析,TERCOM一般僅適用于航跡相對固定的飛行器。TERCOM常采用交叉相關(cross correlation,COR)算法和均方差(mean square difference,MSD)算法。MSD和COR都屬于統(tǒng)計決策法中的最基本匹配算法,其中MSD屬于距離決策,COR算法屬于形狀決策。在具體實現(xiàn)中可采用改進的MSD算法和COR算法消除氣壓高度表和無線電高度表測量中的系統(tǒng)偏差,并采用MSD算法和COR算法融合表決機制以及多數(shù)表決決策機制作為提高定位概率和匹配概率的技術措施,詳細分析參見文獻[10]。
ICCP算法源于圖像配準問題中的ICP(iterative closest point,ICP)算法,采用慣導航跡與高度序列之間的調優(yōu),即通過慣導指示航跡進行剛性旋轉和平移變換來實現(xiàn)迭代配準[14],匹配誤差隨慣導誤差的累積而增加,一般常用于水下地形匹配、水下重力匹配導航和水下地磁匹配等應用領域的研究。
單點迭代算法是同時利用地形高程和從地形高程數(shù)據導出的地形斜率對每一個高程采樣點都進行匹配,并將修正信息提供給慣導。由于單點迭代算法不需要事先規(guī)劃航路,可較好滿足機載實時應用的需求,SITAN算法為典型代表[15],其主要模塊對應的任務及核心技術見表2。
表2 SITAN算法的主要模塊
公開可見的典型應用包括F-l6飛機的地形輔助導航系統(tǒng)和面向直升飛機應用的HELI/SITAN。原理上,單點迭代算法更適合飛行高度高、高度測量誤差大、慣導初始積累誤差大、飛行方式靈活等機載應用場合。以桑迪亞實驗室為直升機應用而研制的HELI/SITAN為例,在滿足定位準則時,就能給出離散點上的定位信息,而不需要進行事先航路規(guī)劃,解決了機載環(huán)境適應性問題,較好地滿足了機載應用需求。
類似其他工程問題,機載地形匹配系統(tǒng)同樣也面臨著存在很大位置不確定性的情況下如何精確捕獲初始位置的困難和復雜性。在HELI/SITAN中采用基于密集并行濾波器的捕獲模式設計技術,以數(shù)字高程地圖的網格間距為子濾波器的搜索間距,在搜索期間濾波器整體上一直保持規(guī)則的網格結構的并行濾波器陣列,對搜索區(qū)內每一個可能的路徑進行遍歷,記錄并跟蹤每一采樣時刻的最優(yōu)匹配濾波器,同時利用慣導緩變、連續(xù)的誤差漂移特性建立相應的定位準則,實時判斷并找到當前的正確濾波器,以正確濾波器的當前位置作為飛機的真實位置的估值,就能獲得定位誤差小于1個網格間距的搜索定位精度。獲得精確的初始位置后,在跟蹤模式中通過相關處理實現(xiàn)信噪比估算、匹配性能的實時監(jiān)控以及提供定位信息和定位信息品質等信息,從而可獲得優(yōu)于2個網格間距的定位精度。在跟蹤模式中需要重點關注衡量匹配程度變量的值以及跟蹤濾波器誤差協(xié)方差估算值的精度,這些值不僅代表了匹配的可靠性及跟蹤模式的性能,而且是進行模式轉換控制邏輯轉換的主要依據。圖6是機載地形匹配系統(tǒng)工作過程示意圖。
圖6 系統(tǒng)工作過程示意圖
具體實現(xiàn)中,采用變化緩慢的一維垂直通道偏差作為系統(tǒng)狀態(tài),該偏差由氣壓高度表偏差、無線電高度表偏差以及數(shù)字高程圖偏差引起。在地形匹配系統(tǒng)運行過程中應注意:上述各個偏差沒有一個能單獨觀測到,但由各個偏差綜合而成的偏差能夠觀測到,需要采用綜合處理的方式進行垂直通道偏差的衡量。系統(tǒng)處理過程為:
1)系統(tǒng)狀態(tài)方程為:
xn=xn-1+wn
(1)
2)系統(tǒng)量測方程為:
zn=zn-1+vn
(2)
3)匹配程度的衡量
當某個濾波器最接近飛機真實位置時,根據該濾波器位置上的地形高程數(shù)據所計算出的地形高程值和載機測量到的地形高程值之差幾乎等于偏差,因此該濾波器最符合先驗偏差模型。
采用平滑加權殘差平方(smoothed weighted residual squared,SWRS)一方面能夠提供單個濾波器與先驗垂直通道偏差模型匹配程度的量度:SWRS值越小,匹配程度越好;另一方面,SWRS的計算過程也容易引入指數(shù)衰減的權,方便計入慣導的漂移,因此可通過SWRS來增強較新匹配信息的作用強度,通過控制衰減因子適應不同精度等級的慣導應用。SWRS的計算過程為:
SWRSjn=aWRSjn+(1-a)WRSj(n-1)
(3)