基于GIS的天氣雷達(dá)Z數(shù)據(jù)補(bǔ)償

林輝 張鐸 楊建斌

摘要：天氣雷達(dá)是監(jiān)測(cè)和預(yù)警強(qiáng)對(duì)流天氣的主要工具，通過(guò)發(fā)射一系列脈沖電磁波，利用云霧、雨、雪等降水粒子對(duì)電磁波的散射和吸收效應(yīng)，來(lái)探測(cè)降水的空間分布和鉛直結(jié)構(gòu)。雷達(dá)探測(cè)波束在低仰角經(jīng)常會(huì)出現(xiàn)被地物遮擋的情況，導(dǎo)致某些方位探測(cè)到的回波會(huì)偏弱。多雷達(dá)多重覆蓋拼圖大多數(shù)時(shí)候可以很好解決這個(gè)問(wèn)題，但是每個(gè)關(guān)注的區(qū)域都做到多雷達(dá)覆蓋在實(shí)際中很難實(shí)現(xiàn)，也可能存在某個(gè)區(qū)域多部雷達(dá)波束均受到遮擋的情況。

關(guān)鍵詞：GIS;天氣雷達(dá);數(shù)據(jù)

引言：

利用GIS高程地形信息，演算出不同仰角下雷達(dá)波束在每個(gè)方位的遮擋情況，計(jì)算出波束遮擋比例以及強(qiáng)度修正值，在幾乎不追加任何成本的前提下就可以對(duì)雷達(dá)回波Z數(shù)據(jù)，即反射率因子進(jìn)行補(bǔ)償。

一、算法流程

算法流程如下：

解析多普勒天氣雷達(dá)的體掃原始數(shù)據(jù)，除了提取Z數(shù)據(jù)以外，還需獲取雷達(dá)體積掃描各仰角角度，雷達(dá)海拔高度，雷達(dá)經(jīng)緯度等信息，在計(jì)算遮蔽信息時(shí)需要用到。根據(jù)雷達(dá)經(jīng)緯度，截取以雷達(dá)為中心，東西跟南北跨度為1度（接近111km）的高程數(shù)據(jù)。由于高程數(shù)據(jù)采用經(jīng)緯度表示方位，因此要用高斯投影算法投影到平面直角坐標(biāo)。 根據(jù)雷達(dá)高度，結(jié)合大氣折射率跟地球曲率，計(jì)算雷達(dá)每個(gè)方位距離的遮蔽角。演算不同仰角下不同方位距離的波束遮擋比例，計(jì)算出強(qiáng)度（Z數(shù)據(jù)）的修正值。

二、算法實(shí)現(xiàn)

（一） 民航多普勒天氣雷達(dá)格式

天氣雷達(dá)探測(cè)波束是細(xì)長(zhǎng)的針狀，C波段探測(cè)有效距離一般是200km。雷達(dá)掃描有水平與俯仰2個(gè)自由度，采用極坐標(biāo)（方位，距離）方式記錄數(shù)據(jù)。雷達(dá)體積掃描就是在不同仰角下，雷達(dá)天線水平轉(zhuǎn)動(dòng)360度，發(fā)出電磁波探測(cè)波束，接收天氣目標(biāo)散射回來(lái)的信號(hào)，數(shù)字化采樣后并保存到體積掃描的原始數(shù)據(jù)文件里。

原始數(shù)據(jù)文件一般由文件標(biāo)識(shí) 、文件頭和徑向數(shù)據(jù)記錄塊組成。文件頭記載雷達(dá)站名、站址、雷達(dá)型號(hào)、主要參數(shù)、觀測(cè)時(shí)間、掃描類(lèi)型、工作狀況等內(nèi)容。? 徑向數(shù)據(jù)記錄以極坐標(biāo)方式順序排列。雷達(dá)一般掃描1到15個(gè)仰角層，在每個(gè)仰角層雷達(dá)天線水平轉(zhuǎn)動(dòng)360度，一般保存為360個(gè)方位的徑向數(shù)據(jù)，一般每個(gè)徑向數(shù)據(jù)是由400個(gè)庫(kù)組成，探測(cè)有效距離200km就分成400個(gè)庫(kù)來(lái)量化，一個(gè)庫(kù)代表了500m內(nèi)采樣的回波的均值。實(shí)際的庫(kù)長(zhǎng)，徑向數(shù)據(jù)庫(kù)數(shù)，各層徑向數(shù)據(jù)記錄數(shù)，體掃層數(shù)與角度都在文件頭內(nèi)有定義。

（二）高程數(shù)據(jù)解析與截取

高程數(shù)據(jù)采用的是格點(diǎn)分辨率約為30米的GDEMV2，數(shù)據(jù)以TIFF格式保存。TIFF文件主要由圖像文件頭，圖像文件目錄，目錄入口，圖像數(shù)據(jù)組成。一個(gè)TIF文件表示經(jīng)度緯度范圍為1°的高程數(shù)據(jù)，分辨率為3601*3601。每個(gè)像素點(diǎn)代表經(jīng)度緯度為1”的高程數(shù)據(jù)，高度值用2個(gè)字節(jié)表示，可表示高度范圍為0到65536米。

根據(jù)雷達(dá)數(shù)據(jù)文件定義的雷達(dá)經(jīng)緯度，計(jì)算雷達(dá)為中心，經(jīng)緯度跨度為1度所需要用到的高程數(shù)據(jù)文件，并截取出所需要的那部分高程文件。高程文件拼接時(shí)，會(huì)出現(xiàn)重疊的邊界，要切除這部分邊界數(shù)據(jù)。

（三）高程數(shù)據(jù)高斯投影

GDEMV2高程數(shù)據(jù)是采用經(jīng)緯度來(lái)表示位置，應(yīng)用上需要通過(guò)投影算法轉(zhuǎn)化成我們平時(shí)用的空間直角坐標(biāo)系，此處采用了高斯投影算法。

（四）遮蔽角計(jì)算

雷達(dá)波束中心高度與目標(biāo)物實(shí)際距離關(guān)系如下公式所示，

為雷達(dá)高度，θ為仰角，h為雷達(dá)波束中心高度，R為目標(biāo)物到雷達(dá)的實(shí)際距離。那么距離為R，高度為h的地物，對(duì)海拔高度為的雷達(dá)形成的遮擋角為

以雷達(dá)站位置為原點(diǎn)，計(jì)算高程數(shù)據(jù)每個(gè)格點(diǎn)所處雷達(dá)方位角以及距離，以及格點(diǎn)高度形成的遮蔽角。由于雷達(dá)探測(cè)波速是接近直線傳播，因此同個(gè)方位角每個(gè)距離庫(kù)范圍內(nèi)波束的遮擋角是由該距離庫(kù)與雷達(dá)位置之間所有距離庫(kù)的最大遮擋角決定。假設(shè)探測(cè)雷達(dá)波束從原點(diǎn)射出，地物B形成遮蔽角于小于地物A的遮蔽角θ，因此雷達(dá)波束上沿（雷達(dá)波束有波瓣寬度，一般是1°）要高于θ才能探測(cè)到地物B上空以及更遠(yuǎn)處的天氣目標(biāo)。地物C雖然高度很低，但是之前地物形成遮蔽角β，且β大于θ，所以雷達(dá)波束上沿要高于β才能探測(cè)到地物C上空以及的天氣目標(biāo)。因此對(duì)每個(gè)方位的地物形成的遮擋角由近到遠(yuǎn)進(jìn)行最大值保持處理，也就是只要當(dāng)前地物形成遮蔽角小于或等于之前最大遮蔽角，當(dāng)前地物遮蔽角就等效于之前最大值。

（五）波束遮擋比例演算與補(bǔ)償值計(jì)算

從雷達(dá)原點(diǎn)發(fā)射出雷達(dá)探測(cè)波束，假設(shè)仰角為β，波束寬度為α。在地物A以及之前，波束下沿角度（β-0.5α）高于地物，遮擋比例為0，探測(cè)的回波不需要修正。地物B處波束上沿角度（β+0.5α）高于地物，下沿低于地物，發(fā)生部分遮擋，地物B形成遮蔽角為θ，則遮擋比例，也就是雷達(dá)波束上沿角與地物遮蔽角的差值占總波束寬度的比例，如下邊公式所示。

在地物C處，雷達(dá)波束上沿低于地物，波束被完全遮擋，該處即更遠(yuǎn)處回波不需要修正。

將遮擋比例換算為對(duì)數(shù)，并對(duì)Z數(shù)據(jù)進(jìn)行修正：

是原始Z數(shù)據(jù)，是修正后的Z數(shù)據(jù)。

三、驗(yàn)證

多雷達(dá)拼圖通常可以獲得很好的地物遮擋補(bǔ)償效果。以汕頭氣象局、梅州氣象局兩部雷達(dá)的拼圖作為參照物。生成汕頭空管站天氣雷達(dá)體掃基數(shù)據(jù)補(bǔ)償前后的最強(qiáng)平面投影圖。通過(guò)對(duì)補(bǔ)償前后的投影圖與作為參照物的拼圖進(jìn)行對(duì)比分析來(lái)驗(yàn)證該算法。

選取了2019年12月05日北京時(shí)間7點(diǎn)06分的三部雷達(dá)的數(shù)據(jù)進(jìn)行分析。采用本算法對(duì)汕頭空管站天氣雷達(dá)體掃基數(shù)據(jù)每個(gè)仰角層Z數(shù)據(jù)進(jìn)行補(bǔ)償，生成補(bǔ)償前后的最強(qiáng)平面投影圖，與汕頭、梅州氣象局兩部雷達(dá)拼圖進(jìn)行比較，拼圖中心點(diǎn)已經(jīng)設(shè)置為汕頭空管天氣雷達(dá)站所在位置。從拼圖看，在空管雷達(dá)站西北跟東北方向，有呈帶狀分布的較強(qiáng)回波，色標(biāo)為淡黃色。而空管站天氣雷達(dá)在強(qiáng)度補(bǔ)償前只探測(cè)到依稀少量強(qiáng)回波，補(bǔ)償后，強(qiáng)回波數(shù)量有所增加，接近帶狀分布，東北方向強(qiáng)回波增加較為明顯，接近拼圖所展示的探測(cè)結(jié)果。可見(jiàn)，該算法對(duì)地物遮蔽具有一定補(bǔ)償作用。

四、結(jié)語(yǔ)

本文闡述了利用GIS高程地物信息，對(duì)天氣雷達(dá)強(qiáng)度Z數(shù)據(jù)進(jìn)行補(bǔ)償?shù)姆椒āＲ陨穷^空管站天氣雷達(dá)體掃數(shù)據(jù)為樣本，參照汕頭、梅州氣象局天氣雷達(dá)拼圖，對(duì)算法有效性進(jìn)行驗(yàn)證，結(jié)果表明該算法從一定程度上修正因?yàn)榈匚镎趽醵鴮?dǎo)致回波強(qiáng)度偏弱情況。