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

1、 基于GIS的天氣雷達(dá)Z數(shù)據(jù)補(bǔ)償 林輝 張鐸 楊建斌Summary：天氣雷達(dá)是監(jiān)測和預(yù)警強(qiáng)對流天氣的主要工具，通過發(fā)射一系列脈沖電磁波，利用云霧、雨、雪等降水粒子對電磁波的散射和吸收效應(yīng)，來探測降水的空間分布和鉛直結(jié)構(gòu)。雷達(dá)探測波束在低仰角經(jīng)常會(huì)出現(xiàn)被地物遮擋的情況，導(dǎo)致某些方位探測到的回波會(huì)偏弱。多雷達(dá)多重覆蓋拼圖大多數(shù)時(shí)候可以很好解決這個(gè)問題，但是每個(gè)關(guān)注的區(qū)域都做到多雷達(dá)覆蓋在實(shí)際中很難實(shí)現(xiàn)，也可能存在某個(gè)區(qū)域多部雷達(dá)波束均受到遮擋的情況。Key：GIS;天氣雷達(dá);數(shù)據(jù)引言：利用GIS高程地形信息，演算出不同仰角下雷達(dá)波束在每個(gè)方位的遮擋情況，計(jì)算出波束遮擋比例以及強(qiáng)度修正值，在幾乎不追

2、加任何成本的前提下就可以對雷達(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á)探測波束是細(xì)長的

3、針狀，C波段探測有效距離一般是200km。雷達(dá)掃描有水平與俯仰2個(gè)自由度，采用極坐標(biāo)（方位，距離）方式記錄數(shù)據(jù)。雷達(dá)體積掃描就是在不同仰角下，雷達(dá)天線水平轉(zhuǎn)動(dòng)360度，發(fā)出電磁波探測波束，接收天氣目標(biāo)散射回來的信號(hào)，數(shù)字化采樣后并保存到體積掃描的原始數(shù)據(jù)文件里。原始數(shù)據(jù)文件一般由文件標(biāo)識(shí) 、文件頭和徑向數(shù)據(jù)記錄塊組成。文件頭記載雷達(dá)站名、站址、雷達(dá)型號(hào)、主要參數(shù)、觀測時(shí)間、掃描類型、工作狀況等內(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è)庫組成，探測有效距離200km就

4、分成400個(gè)庫來量化，一個(gè)庫代表了500m內(nèi)采樣的回波的均值。實(shí)際的庫長，徑向數(shù)據(jù)庫數(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ù)文件，并截取出所需要的那部分高程文

5、件。高程文件拼接時(shí)，會(huì)出現(xiàn)重疊的邊界，要切除這部分邊界數(shù)據(jù)。（三）高程數(shù)據(jù)高斯投影GDEMV2高程數(shù)據(jù)是采用經(jīng)緯度來表示位置，應(yīng)用上需要通過投影算法轉(zhuǎn)化成我們平時(shí)用的空間直角坐標(biāo)系，此處采用了高斯投影算法。（四）遮蔽角計(jì)算雷達(dá)波束中心高度與目標(biāo)物實(shí)際距離關(guān)系如下公式所示，為雷達(dá)高度，為仰角，h為雷達(dá)波束中心高度，R為目標(biāo)物到雷達(dá)的實(shí)際距離。那么距離為R，高度為h的地物，對海拔高度為的雷達(dá)形成的遮擋角為以雷達(dá)站位置為原點(diǎn)，計(jì)算高程數(shù)據(jù)每個(gè)格點(diǎn)所處雷達(dá)方位角以及距離，以及格點(diǎn)高度形成的遮蔽角。由于雷達(dá)探測波速是接近直線傳播，因此同個(gè)方位角每個(gè)距離庫范圍內(nèi)波束的遮擋角是由該距離庫與雷達(dá)位置之間所有距

6、離庫的最大遮擋角決定。假設(shè)探測雷達(dá)波束從原點(diǎn)射出，地物B形成遮蔽角于小于地物A的遮蔽角，因此雷達(dá)波束上沿（雷達(dá)波束有波瓣寬度，一般是1）要高于才能探測到地物B上空以及更遠(yuǎn)處的天氣目標(biāo)。地物C雖然高度很低，但是之前地物形成遮蔽角，且大于，所以雷達(dá)波束上沿要高于才能探測到地物C上空以及的天氣目標(biāo)。因此對每個(gè)方位的地物形成的遮擋角由近到遠(yuǎn)進(jìn)行最大值保持處理，也就是只要當(dāng)前地物形成遮蔽角小于或等于之前最大遮蔽角，當(dāng)前地物遮蔽角就等效于之前最大值。（五）波束遮擋比例演算與補(bǔ)償值計(jì)算從雷達(dá)原點(diǎn)發(fā)射出雷達(dá)探測波束，假設(shè)仰角為，波束寬度為。在地物A以及之前，波束下沿角度（-0.5）高于地物，遮擋比例為0，探測

7、的回波不需要修正。地物B處波束上沿角度（+0.5）高于地物，下沿低于地物，發(fā)生部分遮擋，地物B形成遮蔽角為，則遮擋比例，也就是雷達(dá)波束上沿角與地物遮蔽角的差值占總波束寬度的比例，如下邊公式所示。在地物C處，雷達(dá)波束上沿低于地物，波束被完全遮擋，該處即更遠(yuǎn)處回波不需要修正。將遮擋比例換算為對數(shù)，并對Z數(shù)據(jù)進(jìn)行修正：是原始Z數(shù)據(jù)，是修正后的Z數(shù)據(jù)。三、驗(yàn)證多雷達(dá)拼圖通?？梢垣@得很好的地物遮擋補(bǔ)償效果。以汕頭氣象局、梅州氣象局兩部雷達(dá)的拼圖作為參照物。生成汕頭空管站天氣雷達(dá)體掃基數(shù)據(jù)補(bǔ)償前后的最強(qiáng)平面投影圖。通過對補(bǔ)償前后的投影圖與作為參照物的拼圖進(jìn)行對比分析來驗(yàn)證該算法。選取了2019年12月05日北京時(shí)間7點(diǎn)06分的三部雷達(dá)的數(shù)據(jù)進(jìn)行分析。采用本算法對汕頭空管站天氣雷達(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ǔ)償前只探測到依稀少量強(qiáng)回波，補(bǔ)償后，強(qiáng)回波數(shù)量有所增加，接近帶狀分布，東北方向強(qiáng)回波增加較為明顯，接近拼圖所展示的探測結(jié)果?？梢?，該算法對地物遮蔽具有一定補(bǔ)償作用。四