相控陣天氣雷達(dá)站環(huán)境參數(shù)智能分析系統(tǒng)研發(fā)

呂雪芹，何艷麗，敖振浪

(廣東省氣象探測數(shù)據(jù)中心，廣州 510080)

0 引言

天氣雷達(dá)是探測臺風(fēng)、暴雨、龍卷風(fēng)等強對流天氣的重要技術(shù)手段之一。雷達(dá)探測效果不僅要求雷達(dá)設(shè)備本身具有精準(zhǔn)的探測性能，更是要求雷達(dá)站周邊環(huán)境條件好，在綜合考慮工程建設(shè)基礎(chǔ)條件可行的前提下，要求阻擋角盡量越小越好。一般情況下，雷達(dá)站站址的拔海高度越高，四周阻擋角越小。但是拔海高度不是越高越好，因為天氣雷達(dá)一般不使用負(fù)仰角掃描，所以0°仰角以下的區(qū)域沒有掃描到，就會出現(xiàn)雷達(dá)掃描盲區(qū)。而在低空范圍里，往往又是強對流天氣經(jīng)常發(fā)生的關(guān)鍵區(qū)域。目前廣東省氣象部門快速推進(jìn)X波段雙極化相控陣天氣雷達(dá)建設(shè)，正是彌補原有S波段新一代天氣雷達(dá)探測網(wǎng)的低空探測盲區(qū)的不足。在原有雷達(dá)探測網(wǎng)的基礎(chǔ)上，如何更好的優(yōu)化布局X波段雙極化相控陣天氣雷達(dá)，站點選擇顯得尤為重要。按照以往的定性分析雷達(dá)站探測環(huán)境的常規(guī)做法已經(jīng)不能完全滿足精細(xì)化探測業(yè)務(wù)需要，因此如何更好地定量化計算、分析、判斷雷達(dá)站環(huán)境條件參數(shù)是非常重要的。近些年來有不少關(guān)于雷達(dá)站選址方面的文章，但主要是分析四周阻擋角和等射束高度圖的介紹，如伍洋等[1]借助GIS數(shù)據(jù)分析新一代天氣雷達(dá)站選址,并直觀展現(xiàn)站點的遮蔽角圖和等射束高度圖。景號然[2]等利用SRTM高程數(shù)據(jù)作為選址基礎(chǔ)數(shù)據(jù)，計算得到天氣雷達(dá)在0.5°，1.0°，2.4°仰角上地物遮擋情況；利用高程格點數(shù)據(jù)獲得3個仰角的地物剖面數(shù)據(jù)，獲得站點上空1 km,海拔3 km和6 km等射束高度圖，分析了四川省天氣雷達(dá)網(wǎng)探測環(huán)境，并給出了評估結(jié)果。李明鳳[3]等利用數(shù)字高程模型(DEM)對廣東省已建成的12部新一代天氣雷達(dá)分別進(jìn)行地形阻擋分析，主要是通過制作遮蔽角圖和等射束高度圖分析地形阻擋，根據(jù)阻擋分析探測效果優(yōu)劣。葉飛[4]通過獲取DEM高程數(shù)據(jù)、SRTM數(shù)據(jù)和雷達(dá)周邊高大建筑物的遮擋數(shù)據(jù)，分別繪制天氣雷達(dá)單站遮蔽角圖和內(nèi)蒙全區(qū)新一代天氣雷達(dá)網(wǎng)等射束高度拼圖，對內(nèi)蒙全區(qū)新一代天氣雷達(dá)探測環(huán)境進(jìn)行了評估分析，通過評估分析可以對雷達(dá)覆蓋率不夠的區(qū)域提出補償機制。高永芹、劉鋒、吳煥萍、楊磊等[5-8]介紹了基于組件式GIS技術(shù)的業(yè)務(wù)應(yīng)用系統(tǒng)建設(shè)的總體框架,并闡述了系統(tǒng)的建設(shè)內(nèi)容,主要包括電子地圖建設(shè),制圖產(chǎn)品輸出,實時與歷史資料檢索，評估以及參數(shù)化制圖等功能，對GIS應(yīng)用于決策服務(wù)系統(tǒng)的若干關(guān)鍵技術(shù)問題進(jìn)行了討論與分析。這些研究和系統(tǒng)，基本上是基于GIS地理信息獲取高程數(shù)據(jù)制作雷達(dá)站四周遮蔽角圖和等射束高度圖，再人工目視評估探測環(huán)境，費時費力，缺少定量化分析，評價完整性也存在局限性，不能很好地滿足精密探測業(yè)務(wù)需求。本系統(tǒng)重點解決定量化、自動化評估探測環(huán)境內(nèi)容。

1 架構(gòu)及界面設(shè)計

1.1 總體架構(gòu)

整個系統(tǒng)的架構(gòu)設(shè)計分為四個部分：第一個部分是站點基本信息的錄入，包括站點名稱、經(jīng)緯度、海拔高度、雷達(dá)天線計劃離地面高度等參數(shù)；第二部分是通過GIS獲取雷達(dá)站地理位置360°范圍內(nèi)距離站點一定半徑內(nèi)的高程數(shù)據(jù)[9-11]；第三部分是制作遮蔽角圖和等射束高度圖，進(jìn)而制作多個站點的等射束高度圖的拼圖；第四部分是自動綜合分析探測環(huán)境狀況，輸出綜合分析報告。流程框圖如圖1所示。

圖1 系統(tǒng)架構(gòu)圖

1.2 界面設(shè)計

采用C#編程語言構(gòu)建WINDOWS交互主界面，菜單欄分為站點基本參數(shù)、獲取高程數(shù)據(jù)、探測環(huán)境分析及顯示等4個部分，功能設(shè)計著重考慮了各個部分在實際使用中的關(guān)聯(lián)性和靈活性。使用廣東省連州雷達(dá)站作為背景圖，功能標(biāo)簽和按鍵采取圖形化凸凹立體輪廓，立體感強，對約定可選內(nèi)容采取可編輯文本框，方便輸入。界面力求美觀，交互友好便捷，如圖2所示。

圖2 主界面

在實際應(yīng)用中，對于某個擬選站點，可使用儀器測量它的經(jīng)緯度和海拔高度，根據(jù)雷達(dá)塔樓和雷達(dá)基座計算天線離地面的高度，輸入到對應(yīng)的“站點基本參數(shù)”菜單欄。

設(shè)計可編輯文本框的預(yù)設(shè)半徑范圍分別為50 km、75 km、100 km和150 km。對于X波段雙極化相控陣?yán)走_(dá)來說，有效探測范圍一般是50 km以內(nèi)，選取數(shù)字高程數(shù)據(jù)DEM的范圍一般選擇50 km即可，如果考慮為S波段雷達(dá)所用，則選擇100 km就足夠了，因為地球曲率原因，100 km以外只有超過8km以上的高山才可能產(chǎn)生阻擋[12-13]。讀取的數(shù)字高程數(shù)據(jù)存入對應(yīng)站點的TXT文件中，TXT文件還包含基本參數(shù)等信息，作為后續(xù)分析和評估所需基本數(shù)據(jù)。

考慮使用的靈活性，通過“獲取高程數(shù)據(jù)”菜單欄中一鍵讀取數(shù)據(jù)，系統(tǒng)自動搜索指定半徑范圍的高程數(shù)據(jù)，生成TXT文件保存于磁盤中用于任何時候使用。用戶可以在TXT文件中根據(jù)實際狀況增加或者修改高程數(shù)據(jù)，例如附近有建筑物阻擋而GIS里又沒有標(biāo)示的點。

設(shè)計“探測環(huán)境分析”菜單欄一次性完成阻擋角圖、等射束高度圖和多站拼圖的制作，生成對應(yīng)的JPG圖片，同時，自動分析不同仰角的連排阻擋物的方位角以及全方位占比，進(jìn)一步自動計算靜錐盲區(qū)、低空盲區(qū)、掃描面積和掃描立體空域完備度等，再自動分析區(qū)域內(nèi)的環(huán)境條件情況，輸出綜合分析結(jié)論。

2 高程數(shù)據(jù)讀取方法

一般來說，從網(wǎng)絡(luò)上商用的GIS信息系統(tǒng)在線讀取高程數(shù)據(jù)是有限制的，不方便且速度慢，不太適合雷站址選擇的應(yīng)用場景需求[14-15]。本文重點介紹離線讀取的方法，首先從地理信息系統(tǒng)下載相應(yīng)區(qū)域的TIF格式高程DEM文件，TIF格式高程文件很大，涵蓋一個省的矩形區(qū)域的文件大小至少2 GB，占用太多計算機資源，影響運行性能，所以選擇區(qū)域應(yīng)該適當(dāng)，夠用就行了。數(shù)字高程的分辨率當(dāng)然越高越好，本設(shè)計下載30米距離分辨率1米高度分辨率的數(shù)據(jù)，基本滿足業(yè)務(wù)需求。通過文件讀取的方法離線讀取數(shù)字高程數(shù)據(jù)，這樣速度非?？靃16]，而且使用起來方便，完全不依賴于互聯(lián)網(wǎng)即可隨時隨地使用，這對于野外現(xiàn)場選址的時候至關(guān)重要。下面是實現(xiàn)的主要代碼。在菜單點擊“獲取高程”按鍵，直接進(jìn)入事件代碼。

private void buttonScanning_Click(object sender, EventArgs e)

Gdal.AllRegister(); Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES"); //支持中文

Dataset ds = Gdal.Open("D:雷達(dá)站選址智能測評程序GIS基礎(chǔ)信息文件(不可刪除)GIS_DEM.tif", Access.GA_ReadOnly);

int rasterX = ds.RasterXSize; //影像寬度

int rasterY = ds.RasterYSize; //影像高度

int bandCount = ds.RasterCount; //波段數(shù)

double[] tmpD = new double[6];

ds.GetGeoTransform(tmpD); //影像坐標(biāo)變換參數(shù)

string proj = ds.GetProjection(); //影像坐標(biāo)系信息(WKT格式字符串)

backgroundWorker1=new System.ComponentModel.BackgroundWorker();

backgroundWorker1.DoWork+= ackgroundWorker1_DoWork;//注冊后臺操作

backgroundWorker1.RunWorkerCompleted+= backgroundWorker1_RunWorkerCompleted;

Control.CheckForIllegalCrossThreadCalls= false;//.NET禁止了跨線程調(diào)用控件 backgroundWorker1.RunWorkerAsync(); //啟動后臺操作

buttonScanning.Enabled = false;// 在獲取資料期間禁止再次按該按鈕

textBox7.ForeColor = Color.Black;

textBox7.Text = "正在讀取高程數(shù)據(jù)......";

button1.Enabled = true;

眾所周知，從地圖讀取數(shù)據(jù)是需要花費一定時間的。如果前臺長時間讀數(shù)據(jù)的時候不能進(jìn)行其他操作，則非常乏味且不可接受。采取的方法就是將前臺讀取數(shù)據(jù)的任務(wù)通過委托的方法放在后臺運行，前臺就能夠響應(yīng)用戶操作。使用委托語句backgroundWorker1.DoWork += backgroundWorker1_DoWork； 進(jìn)行注冊后臺操作，把讀取任務(wù)放在后臺，前臺無感知。后臺讀取事件由如下主要代碼實現(xiàn)。

private void backgroundWorker1_DoWork(object sender, DoWorkEventArgs e)

string filePathOne = "D:雷達(dá)站選址智能測評程序數(shù)據(jù)資料demo.TXT" StreamWriter Rad_fp = new StreamWriter(filePathOne);//實例化打開一個文本文件

//掃描四周的阻擋物的海拔高度值

for (loop = 0; loop < LOOP; loop++)

float AVG = 0, MaxAVG = 0;

Distant = 0;

NewDistant = 0;

MaxHight = 0;

for (i = 0; Distant < DataAskRange; i++)//尋找每條掃描線上的最高點

Distant = Distant + 50; //按照每50米距離間隔查詢一個點的海拔高度

double[] LatLon = new double[2];

LatLon = ComputerNewPointLatLon.computerThatLonLat(Lon0, Lat0, Distant, Angle);//計算新的坐標(biāo)的經(jīng)緯

AVG = ProcessGisTiffFile.GetMultifyElevation(ds, LatLon[0], LatLon[1]);

hight = (float)(0.5 * Distant * Distant / radiusEarthMetres);

hight = AVG - hight;// 地球曲率訂正

HigthTemp = (float)(MaxAVG - StationHeight - AntennaHeight);

BlockingAngle = (float)(Math.Asin(((HigthTemp-(0.5*NewDistant* NewDistant/radiusEarthMetres))/ NewDistant)) * 180 / Math.PI); //計算某方位角點的阻擋角

//保存數(shù)據(jù)到TXT文件

Rad_fp.Write((DataAskRange / 1000).ToString().PadLeft(5, ' ') + " ");//測評范圍

Rad_fp.Write(Angle.ToString().PadLeft(5, ' ') + " ");//方位角

Rad_fp.Write((NewDistant / 1000.0).ToString("F1").PadLeft(6, ' ') + " ");//最高點的距離

Rad_fp.Write((MaxAVG).ToString("F1").PadLeft(6, ' ') + " ");//最高點的高度

Rad_fp.Write((MaxHight).ToString("F1").PadLeft(6, ' ') + " ");//訂正后的高度

Rad_fp.Write(BlockingAngle.ToString("F2").PadLeft(6, ' ') + " ");//阻擋角

Angle = Angle + AngleInterval;//增加一個間隔的方位角

if (backgroundWorker1.CancellationPending) //檢查是否有中止后臺繼續(xù)執(zhí)行的按鈕操作，有則退出

e.Cancel = true;

break;

Rad_fp.Flush();//清除內(nèi)存釋放資源

Rad_fp.Close();//關(guān)閉資料文件

為了便于后續(xù)處理，掃描到的高程數(shù)據(jù)存儲到TXT文件的格式如表1所示。

表1 存儲格式

需要注意的是，以站點為中心掃描四周的阻擋物，究竟以多少方位角間隔和距離間隔進(jìn)行步進(jìn)掃描，應(yīng)該由實際需求決定[17]。本設(shè)計以0.25°間隔和50 m距離間隔在掃描范圍內(nèi)掃描地形的海拔高度值，從實際效果來看還是比較理想的。

3 遮蔽角和等射束高度圖及拼圖

3.1 遮蔽角和等射束高度圖的制作

評價雷達(dá)探測能力的三種圖[18]的制作流程如圖3所示。

圖3 阻擋圖等制作流程圖

首先打開某一個單站的txt文件，逐行讀取全部數(shù)據(jù)，分別對每一個方位對應(yīng)的原始高度值和距離值進(jìn)行提取，作進(jìn)一步的距離訂正。利用三角函數(shù)計算方位角對應(yīng)點的阻擋角，HigthTemp=hight-StationHeight-AntennaHeight;BlockingAngle=(float)(Math.Atan(HigthTemp/(NewDistant*1000))*180/Math.PI);其中HigthTemp是雷達(dá)天線訂正后的高度。NewDistant是阻擋物與雷達(dá)站之間的距離。

在極坐標(biāo)系繪制阻擋圖和在直角坐標(biāo)系上繪制阻擋圖。從圖4可以看出，制作的阻擋角圖既具有同圓極坐標(biāo)形式，也具有對應(yīng)的直角坐標(biāo)形式，比傳統(tǒng)只有極坐標(biāo)形式更加形象直觀，容易理解，如圖4所示。

圖4 阻擋角圖

天氣雷達(dá)探測能力不僅受雷達(dá)性能、各種衰減、電磁波折射和降水云性質(zhì)等因素的影響，也往往會受到雷達(dá)站四周的高大建筑物、山體的阻擋，有的影響相當(dāng)嚴(yán)重[19-20]。為了具體、定量地分析掌握雷達(dá)站在各個方向的探測能力，需要制作表征探測能力的等射束高度圖。

等射束高度圖是指在一定折射條件下，測站四周由于地物阻擋，繪制的各個方向上、各種斜距下，雷達(dá)波束中心軸線能夠到達(dá)的最低高度等值線圖。等射束高度圖考慮了雷達(dá)天線高度的情況，地球球面及大氣折射的影響．特別是它考慮了周圍地物、地形對雷達(dá)探測能力的影響[21]，這對分析雷達(dá)探測環(huán)境條件是非常有用的工具。應(yīng)用測高公式：

H=h+1000Rsinδ+0.058 9R2

(1)

式(1)中,h為天線高度(m)，δ為遮蔽角(°)，H為觀測高度(m) ，R是最大探測距離以km為單位。根據(jù)遮蔽角δ，業(yè)務(wù)上常常分別設(shè)置H為1 km、3 km和6 km，可以分別計算雷達(dá)最大觀測距離R。

以0.25°方位間隔的各個方向上，應(yīng)用測高公式可以分別求出各個方向上1 km高度的雷達(dá)可以達(dá)到的最大探測距離，采用極坐標(biāo)的形式將探測距離依次標(biāo)出，順序連線形成等射束高度圖，如圖5 所示。

圖5 等射束高度圖

3.2 多站拼圖的制作

多站點拼圖實際上是將每個雷達(dá)站點的位置分別在所處的經(jīng)緯度上繪制等射束高度圖，單站使用不同顏色繪制，形成一幅多站拼圖。必須非常注意：低緯度地區(qū)應(yīng)該使用麥卡托投影坐標(biāo)轉(zhuǎn)換變成真實的地理坐標(biāo)[22]！以減少位置誤差。等高度通常取1 km、3 km和6 km。由于參與拼圖的站點相距的距離可能不一樣，有些相距很近，有些相距很遠(yuǎn)，相距太近繪出來拼圖容易密集重疊不好看，如果相距太遠(yuǎn)又會超出屏幕范圍。因此，在屏幕固定區(qū)域范圍里繪制多個站點的拼圖，為了盡可能的利用屏幕可顯示范圍，增強拼圖的清晰度，采用了自動適應(yīng)屏幕大小的技術(shù)。這里通過計算最左、最右、最上、最下的站點經(jīng)緯度，計算出拼圖的中心位置，確定向四邊延伸或者縮小一定比例，進(jìn)行歸一化處理自適應(yīng)屏幕大小。海拔高度1 km、單站半徑45 km的拼圖如圖6所示。

圖6 拼圖

4 自動分析功能實現(xiàn)

打開單站txt文件，逐行讀取各個參數(shù)。統(tǒng)計和分析三個設(shè)定門限的阻擋面，即阻擋角0.5°以上、1.0°以上和1.5°以上。假定如果阻擋角0.5°以上連續(xù)有10°方位阻擋，那么就視作連排阻擋。同理，阻擋角1.0°以上連續(xù)7°以上方位阻擋，就視作連排阻擋，阻擋角1.5°以上連續(xù)5°以上方位阻擋，就視作連排阻擋。根據(jù)阻擋嚴(yán)重程度分為優(yōu)、良、一般或者差等4個等次，作為阻擋優(yōu)劣的判據(jù)在輸出報告中給出建議。下面是主要代碼實現(xiàn)。

TempSaveFileName ="客觀分析報告.txt";

StreamWriter fp1 = new StreamWriter(TempSaveFileName);//實例化打開一個文本文件

fp1.Write("{0}雷達(dá)站探測環(huán)境凈空條件客觀分析報告： ", StationName);

float sum = 0, percent;

int flagCnt = 0;

bool Startflag = false;

int More10BlockingAngle = 0;

for (loop = 0; loop < LOOP; loop++)

if (Angle_BlockingAngle[loop, 1] >= 0.5)//阻擋角>=0.5的占比統(tǒng)計

{ ++sum;

if (Startflag == true)

{ flagCnt++;

} else

{ Startflag = true; flagCnt = 1;

else

{ Startflag = false;

if (flagCnt >= 20)//步進(jìn)0.5°，方位連排0.5X20=10°以上

{ fp1.Write("方位從{0}度到{1}度的區(qū)段，共有{2}度范圍，并且擋角高于0.5°的連排(連續(xù)多于10°方位角)山體阻擋。 ", Angle_BlockingAngle[loop - flagCnt, 0], Angle_BlockingAngle[loop - 1, 0], flagCnt/2);

flagCnt = 0;

More10BlockingAngle++;

percent = sum * 100 / LOOP;

fp1.Write("本站四周高于0.5°阻擋角一共有{0}個點(方位0.5°間隔一個點),", sum);

fp1.Write("即是高于0.5°的阻擋面占比{0}%。 ", percent.ToString("#0.0"));

探測凈空面積的計算方法是采用作者提出的細(xì)分逼近法，限于篇幅，細(xì)分逼近算法不在這里介紹。使用這種算法計算出X波動雙極化相控陣?yán)走_(dá)45 km掃描半徑范圍內(nèi)，1 km高度探測凈空投影面積、凈空占總投影面積比例、投影阻擋面積比例、有效探測空域及其占比等。對于其他類型天氣雷達(dá)，不同掃描范圍和3 km、6 km高度探測凈空面積的計算方法類同。通過自動分析功能，給出對應(yīng)的綜合客觀評價。下面是獨樹崗雷達(dá)站探測環(huán)境凈空條件客觀分析報告實例。

方位從82.25°～98°的區(qū)段，共有32°范圍，并且擋角高于0.5°的連排(連續(xù)多于10°方位角)山體阻擋。

本站四周高于0.5°阻擋角一共有515個點(方位0.25°間隔一個點),即是高于0.5°的阻擋面占比35.8%。

方位從86.75～98°的區(qū)段，共有23°范圍，并且擋角高于1.0°的連排(連續(xù)7度方位角以上)山體阻擋。

本站四周高于1.0°阻擋角一共有118個點(方位0.25°間隔一個點),即是高于1.0°的阻擋面占比8.2%。

本站四周有1個方向存在高于1.0°阻擋的連排(連續(xù)多于7°方位角)區(qū)域，反映探測環(huán)境條件不太理想。

方位從89.5°～97°的區(qū)段，共有15°范圍，并且擋角高于1.5°的連排(連續(xù)5°方位角以上)山體阻擋。

本站四周高于1.5°阻擋角一共有42個點(方位0.25°間隔一個點),即是高于1.5°的阻擋面占比2.9%。

本站四周有1個方向存在高于1.5°阻擋的連排(連續(xù)多于5°方位角)區(qū)域，說明探測環(huán)境條件不理想。

本站45 km掃描半徑、饋源上空25 km以下范圍內(nèi)，有效探測立體空域占總空域體積的百分比(雷達(dá)站45 km探測半徑、饋源上空25 km以下范圍內(nèi)):80.64%。

本站45 km掃描半徑、饋源上空25 km以下范圍內(nèi)，有效探測立體空域占默認(rèn)掃描空域的百分比(雷達(dá)站0～37°仰角、雷達(dá)饋源上空0～25 km高度范圍內(nèi)):98.85%。

本站45 km掃描半徑、饋源上空25 km以下范圍內(nèi)，靜椎區(qū)空域體積占總空域體積的百分比:18.05%。

本站45 km掃描半徑、饋源上空25 km以下范圍內(nèi)，低空掃描盲區(qū)占總空域體積的百分比(天線饋源以下、仰角0°以下):0.38%。

5 結(jié)束語

以往雷達(dá)選址的通常做法，是使用經(jīng)緯儀在擬選點實地測量四周的遮蔽角，這種人工方法相當(dāng)費時費力和粗糙，非常容易受地理位置和天氣狀況影響，較遠(yuǎn)地方阻擋物看不見而忽視了其阻擋影響。也有采用人工在地圖上作業(yè)，尋找各個方位角上最大海拔高度的點，再使用尺子測量與雷達(dá)站之間的距離，然后計算阻擋角。所有這些方法很難保證客觀地反映實際探測環(huán)境狀況，而且效率非常緩慢，往往需要幾天時間才能完成，難以滿足現(xiàn)代業(yè)務(wù)需要。

本文基于離線的GIS地理信息系統(tǒng)數(shù)字高程模型DEM獲取雷達(dá)站探測范圍內(nèi)的山脈、森林、高樓等阻擋情況，非常適合現(xiàn)場初步尋找天氣雷達(dá)最佳候選站址，靈活性便利性不言而喻。自動畫出雷達(dá)探測范圍內(nèi)的阻擋圖、測站上空1 km和海拔3 km、6 km高度的等高度射束圖、多站拼圖，自動計算探測凈空投影面積、有效探測空域等覆蓋率指標(biāo)，并進(jìn)行智能分析，給出最優(yōu)合理化評估報告。系統(tǒng)能夠在數(shù)十秒內(nèi)完成整個計算、畫圖及分析過程，其客觀性和效率凸顯。引入有效探測空域的分析和計算功能是傳統(tǒng)雷達(dá)選址方法不具有的，定量化的分析極大地提高了天氣雷達(dá)站選址決策的科學(xué)性，提高全網(wǎng)探測能力，為雷達(dá)監(jiān)測精密提供良好基礎(chǔ)，該系統(tǒng)在廣東省26部X波段雙極化相控陣?yán)走_(dá)選址中得到很好的應(yīng)用。

為了便于讀者參考借鑒，給出有關(guān)的實現(xiàn)代碼以期對讀者有用。