基于數字地圖的電波傳播特性預測方法

侯松高

(煙臺特種設備檢驗研究院，山東 煙臺 264001)

基于數字地圖的電波傳播特性預測方法

侯松高

(煙臺特種設備檢驗研究院，山東 煙臺 264001)

針對實際地理環境中的電波傳播特性預測問題，提出了一種基于數字地圖的電波傳播特性預測方法。該方法基于拋物方程建立地形條件下的電波傳播預測模型，研究從DEM數字地圖中提取任意二維地形剖面的方法，最后形成基于數字地圖的電波傳播預測流程。仿真結果表明，該方法能夠基于數字地圖有效預測實際地理環境中的電波傳播特性，具有較大的工程應用價值。

地形；數字地圖；拋物方程；電波傳播

0 引言

隨著電子與信息技術的不斷進步，無線通信、導航定位、雷達與遙感探測等系統得到迅速發展并廣泛應用[1-2]，而這些系統中信息的傳輸都以電磁波為載體。在實際環境中，電磁波的傳播過程非常復雜，會受到山地、高原、丘陵、平原和盆地等不規則地形和大氣環境影響，引起反射、折射、繞射及散射等現象，會出現多徑效應。拋物方程(Parabolic Equation,PE)法能是一種前向全波法，能同時處理不規則地形和復雜大氣結構對電波傳播的影響，在對流層電波傳播問題中得到廣泛應用[3-4]。

基于PE的地形條件下電波傳播特性預測，國內外學者開展了大量研究，但多側重于理論方法研究[5-6]，而與真實的地理環境、地表信息結合較少，在實際工程中的應用研究有待進一步深入。而數字地圖能夠較為準確地描述地球表面的地理特征，其中，數字高程模型(Digital Elevation Model，DEM)是地形起伏的數字描述。為此，本文研究基于數字地圖的電波傳播特性預測方法，首先利用PE法建立地形條件下電波傳播預測模型，然后研究從DEM數字地圖中提取任意二維地形剖面的方法，形成基于數字地圖的電波傳播預測流程，最后通過仿真算例證明該方法的有效性。

1 電波傳播預測模型

PE模型是從Maxwell方程組出發，由二維標量波動方程作前向傳播近似，把橢圓型微分方程簡化為拋物型微分方程形式。當阻抗邊界為不規則地形時，對不規則地形邊界進行寬角平移變換，并進行Feit-Fleck近似，可得寬角PE模型[7]：

(1)

式中，u為二維標量場，k(=2π/λ)為自由空間波數，n為傳播媒介的折射率，T為描述地形起伏的函數，T″為地形函數的二階導數，x、z分別為直角坐標系下的傳播距離與高度。

采用SSFT算法可求得式(1)的步進解：

(2)

式中，F、F-1分別為傅里葉變換與逆變換。m(x,z)=n(x,z)-zT″(x)是平移變換后的修正折射率，不僅包含大氣折射因素，還反映了地形邊界對電波傳播的影響[8]。若將真實的大氣環境與地形數據信息應用到PE模型中，則可提高該模型的實際應用價值。

2 地理信息獲取

2.1 DEM數據

DEM作為地形的數字描述，是地學分析計算的基礎數據，被廣泛應用于土木工程、城市規劃以及軍事模擬等諸多領域。描述地形表面形態屬性的數字地形模型(Digital Terrain Model，DTM)，帶有空間位置和地形屬性兩類特征信息，若地形屬性取高程信息，即為DEM。從數學角度，DEM可看作表示地形起伏的高程屬性z關于二維地理空間平面坐標x、y的一種離散函數。

SRTM 3和ASTER GEDEM是目前可以免費獲取的最完整、可靠的2種全球DEM[9-10]。其中SRTM[11-12]是由美國、德國和意大利航天局共同實施的，采用雷達測圖技術對北緯60°～南緯56°之間地球上80%的陸地進行測繪，分辨率有1″(約30 m)和3″(SRTM 3)2種，其中，分辨率為1″的SRTM 1僅供美國本土使用，分辨率為3″的SRTM 3在全球均可免費獲取。ASTER GDEM是由NASA與日本經濟產業省(METI)共同推出的項目，屬于美國地球觀測系統(EOS)的一部分，該數據覆蓋北緯83°～南緯83°間的所有陸地，占地球陸地表面的99%，分辨率為1″(約30 m)。這2類數據的高程基準皆為EGM96，平面基準皆為WGS84。本文主要采用GeoTiff格式存儲的SRTM 3DEM數據。

圖1中，x表示沿經度遞增的方向，y表示沿緯度遞減的方向。該存儲方式是以左上角的格網點為起點，以一定采樣間隔，按行從上往下、列從左到右的順序存儲各網格交叉點的高程值。讀取GeoTiff格式DEM數據時，計算出待求交叉點與起始點位置之間的經緯度偏移量就可獲取該點的高程值。

圖1 GeoTiff格式的DEM數據

2.2 地形剖面的提取

在實際工程應用中，地理環境下的電波傳播預測需要從DEM中提取任意兩點間的二維地形剖面。由2.1節可知，直接讀取DEM數據，只能得到網格交叉點處的高程值，而電波傳播計算中采樣點往往是隨機的，為了獲得更加準確的地形剖面，需要采用內插算法獲得傳播路徑上任意點的高程值。

DEM內插就是根據若干個相鄰參考點的高程值計算出待定點的高程值。根據內差點的分布范圍，內插方法分為整體內插、分塊內插和逐點內插3類[13]，每一類方法中又有多種插值方式。其中分塊內插中的規則格網雙線性多項式內插法是目前DEM分析中最常用的方法，本文在地形剖面抽取中采用該方法。

雙線性多項式內插法[14]是使用與采樣點最近的4個已知節點確定一個雙線性多項式，來計算采樣點的高程值。設采樣點與4個已知點的經緯度和高程值滿足函數關系：

z=a0+a1x+a2y+a3xy，

(3)

式中，a0、a1、a2、a3是待求參數。將4個已知點的經緯度和高程值代入式(3)可得一個四元一次方程組，由此可求出4個參數，再將采樣點的經緯度坐標代入式(3)中，便可求得該點處的高程值。若DEM數據呈正方形規則格網分別，則可直接采用下式：

(4)

式中，A、B、C、D為正方形的4個格網點，l為格網邊長，P為采樣點。

3 電波傳播特性預測仿真

3.1 預測流程

根據2節可知，由DEM數字地圖可獲得網格節點的高程信息，通過雙線性多項式內插法可以計算非節點的高程值，從而可以實現從DEM 中提取任意兩點之間的二維地形剖面，將該地形信息應用于第1節的電波傳播預測模型中，可實現地理環境下的電波傳播預測，具體流程如下：

① 根據電波傳播預測需求，確定發射點與接收點所在位置的經緯度，并由DEM獲取相應位置的高程值；

② 根據收發位置的經緯度和DEM的分辨率，計算出二維地形的水平距離以及所有節點的高程值；

③ 根據PE計算的步進數確定地形采樣點數，通過雙線性內插法計算出采樣點的高程值，從而確定收發點之間的二維地形剖面；

④ 根據地形剖面信息，確定單個步進上PE計算的下邊界以及修正折射率，通過式(2)可確定該步進上的場分布；

⑤ 重復步驟④，經過若干次步進迭代可實現二維地形剖面整個計算區域的場分布，從而實現二維地理環境下的電波傳播預測。

3.2 仿真算例

發射點的經緯度坐標為(E112.05°，N34.29°)，接收點的經緯度坐標為(E112.69°，N33.15°)，采用上述方法從SRTM 3 DEM中提取二維地形剖面，如圖2所示。發射天線高50 m，工作頻率1 000 MHz，全向天線，水平極化；標準大氣環境，地表的相對介電常數為20，電導率為10 mS/m。采用上述地形條件下寬角PE模型預測標準大氣條件下圖2地形剖面上的電波傳播特性，如圖3所示。

圖2 二維地形剖面

圖3 傳播損耗空間分布偽彩圖

分析圖3可知，電波傳播受復雜地形的影響，在山峰背后場強發生衰減，尤其在最高峰之后場強明顯變小。其中在近距離范圍內，由于較大高度上對應的計算仰角很大，超出了Feit-Fleck型PE的計算仰角，圖3中未給出較大高度點上傳播損耗的有效計算結果。

圖4為海拔高度hf=1 600 m的傳播損耗值，其中，左坐標表示地形高程值，右坐標表示電波傳播損耗值。從圖中可以看出，在第3個山峰(次高峰)48.7 km之前，傳播損耗隨距離增大緩慢增長，且曲線發生上下震蕩。這是由于在該高度處地形對電波傳播的影響很小，近于自由空間傳播；且由于山峰的存在引起電波干涉致使曲線震蕩。在最高峰和次高峰之間，受地形影響，傳播損耗值急劇變化。在最高峰之后，傳播損耗值陡增。隨著距離增大，地形高度降低，在該高度處電波傳播受地形的影響變得很小，傳播損耗值增長緩慢。上述結論與文獻[15]采用Longley-Rice模型對地形影響的分析相一致。

圖4 1 600 m高度的傳播損耗隨距離的變化

圖5顯示了不同距離處的傳播損耗隨高度的變化。距離70 km處，位于最高峰之后，在山峰高度以上區域，傳播損耗隨著高度變化較小，且有震蕩現象；山峰高度以下，電波繞射波能夠到達的區域，隨著高度的下降傳播損耗變化相對平緩，而隨著陰影區域的深入，傳播損耗急劇增大。距離120 km處，地形變化相對平緩，地形對電波傳播影響較小，且隨著高度的增大而減小。

圖5 不同距離處傳播損耗隨高度的變化

4 結束語

針對實際地理環境中的電波傳播特性預測問題，從DEM數字地圖中提取真實地理環境下的地形數據，并采用雙線性多項式內插法獲得任意兩點間的地形剖面，利用PE模型預測了該地形剖面上的電波傳播損耗，并通過仿真算例證明了該方法的有效性。將電波傳播特性預測模型與具體的地理環境相結合，提高了電波傳播預測研究的工程實用性，尤其是在通信、電子對抗等領域的實際工程中具有較大的應用價值。

[1] Levy M F.Parabolic Equation Methods for Electromagnetic Wave Propagation[M]. London:IEE Press,2000:33-34,287-290.

[2] 張青洪,廖成,盛楠,等.拋物方程方法的非均勻網格技術研究[J].電波科學學報,2013,28(8):635-640.

[3] 劉曉娣,周新力,侯松高.電波傳播拋物方程模型在航空通信中的應用[J].電訊技術,2016,56(6):624-628.

[4] 侯松高,劉曉娣,周新力.一種低空突防中的雷達探測盲區預測方[J].無線電工程,2016,46(9):45-47.

[5] 常以濤,胡繪斌,姜永金,等.寬角拋物方程中不規則地形的寬角移位變換處理[J].空軍雷達學院學報,2011,25(2):101-104.

[6] 劉曉娣,肖金光,周新力.基于PE模型的地形條件下電波傳播損耗預測[J].計算機仿真,2015,32(8):221-234.

[7] Kuttler J R,Dockery G D.Propagation Modeling Over Terrain Using the Parabolic Wave Equation[J].IEEE Transactions on Antennas and Propagation,2000,48(6):260-277.

[8] 劉曉娣,周新力,肖金光,等.一種地形條件下電波傳播入射余角估計算法[J].電波科學學報,2015,30(6):1211-1217.

[9] 張朝忙,劉慶生,劉高煥,等.SRTM 3與ASTER GDEM數據處理及應用進展[J].地理與地理信息科學,2012,28(5):29-34.

[10]趙國松,杜耘,凌峰,等. ASTER GDEM與SRTM 3高程差異影響因素分析[J]. 測繪科學,2012,37(4):167-170.

[11]Bamler R.The SRTM Mission:A World-wide 30 m Resolution DEM from SAR Interferometry in 11 days[C]∥Wichmann Verlag Heidelberg,Photogrammetric Week,1999:145-153.

[12]Rodriguez E,Morris C S,Belz J E. A Globel Assessment of the SRTM Performance[J].Photogrammetric Engineering & Remote Sensing,2006,72(3):249-260.

[13]盧霞,郭淑艷,劉付程.典型鹽沼植被區數字地形模型的構建[J].淮海工學院學報(自然科學版),2011,20(4):60-64.

[14]李志林,朱慶.數字高程模型[M].武漢:武漢測繪科技大學出版社,2000.

[15]張文波,曹耀欽.電磁環境仿真中電波傳播模型研究及仿真分析[J].電波科學學報,2012,27(3):538-542.

Method of Radio Propagation Characteristics Prediction Based on Digital Map

HOU Song-gao

(Yantai Special Equipment Inspection Institute,Yantai Shandong 264001,China)

A new method of radio propagation characteristics prediction is proposed to improve the practicability of calculation in geographical environment. Radio propagation model on terrain based on parabolic equation is used. The method of extracting the two-dimensional terrain elevation profile from DEM digital map is developed. Then the routine for predicting radio propagation characteristics based on digital map is built. Simulation results show that this method could effectively predict radio propagation characteristics in practical geographical environment.

terrain;digital map;parabolic equation;radio propagation

10. 3969/j.issn. 1003-3114. 2017.04.19

侯松高.基于數字地圖的電波傳播特性預測方法[J].無線電通信技術,2017,43(4):80-82,95.

[HOU Songgao. Method of Radio Propagation Characteristics Prediction Based on Digital Map [J]. Radio Communications Technology,2017,43(4):80-82,95.]

2017-03-24

國家自然科學基金項目(41476089)

侯松高(1981—)，男，助理工程師,主要研究方向：電波傳播、應用電子技術、特種設備檢驗檢測。

TN011

A

1003-3114(2017)04-80-3