基于雙線性內(nèi)插值的 CD雷達(dá)垂直剖面算法

熊 偉,周明飛,潘徐艷

(1.貴州省氣象臺(tái),貴州 貴陽(yáng) 550002;2.貴州省氣候中心,貴州 貴陽(yáng) 550002)

基于雙線性內(nèi)插值的 CD雷達(dá)垂直剖面算法

熊 偉1,周明飛1,潘徐艷2

(1.貴州省氣象臺(tái),貴州 貴陽(yáng) 550002;2.貴州省氣候中心,貴州 貴陽(yáng) 550002)

該文針對(duì)貴州 CD型號(hào)新一代多普勒天氣雷達(dá)反射率因子,采用 C#編程,經(jīng)坐標(biāo)轉(zhuǎn)換后處理成三維網(wǎng)格數(shù)據(jù),基于雙線性內(nèi)插值法,實(shí)現(xiàn)任意方向雷達(dá)反射率因子垂直剖面制作,結(jié)合貴陽(yáng)雷達(dá)分析,結(jié)果表明雙線性內(nèi)插值法放大后的圖像質(zhì)量較高,解決了因垂直剖面圖像較小不利于顯示及放大的問題。

CD雷達(dá);垂直剖面;雙線性內(nèi)插值;圖像

1 引言

目前,新一代多普勒天氣雷達(dá)在天氣預(yù)報(bào)業(yè)務(wù)中得到了廣泛的應(yīng)用,主要針對(duì)災(zāi)害性天氣,特別是短時(shí)強(qiáng)降水、雷雨大風(fēng)、冰雹等災(zāi)害性天氣的監(jiān)測(cè)和預(yù)警。為了清楚地看出風(fēng)暴的回波結(jié)構(gòu),需分析雷達(dá)回波特征及對(duì)風(fēng)暴相應(yīng)位置作垂直剖面分析。貴州所有的雷達(dá)均為成都 784廠 CD型號(hào)新一代多普勒天氣雷達(dá) (以下簡(jiǎn)稱 CD雷達(dá)),目前沒有成熟的基于雷達(dá)反射率因子的垂直剖面分析軟件,而業(yè)務(wù)工作又急需實(shí)現(xiàn)貴州 CD雷達(dá)反射率因子垂直剖面分析及高質(zhì)量圖像的顯示應(yīng)用。

本文針對(duì)貴州 CD雷達(dá)反射率因子,采用 C#編程,經(jīng)坐標(biāo)轉(zhuǎn)換后處理成三維網(wǎng)格數(shù)據(jù),實(shí)現(xiàn)任意方向垂直剖面分析,利用雙線性內(nèi)插值法解決垂直剖面圖像顯示放大且不失真的現(xiàn)象,結(jié)合貴陽(yáng)雷達(dá)分析雙線性內(nèi)插值法在雷達(dá)垂直剖面圖像中的應(yīng)用。

2 CD雷達(dá)反射率因子分析處理

2.1 極坐標(biāo)下 CD雷達(dá)反射率因子轉(zhuǎn)換為三維網(wǎng)格數(shù)據(jù)

CD雷達(dá)理論探測(cè)距離為 250Km,實(shí)際有效探測(cè)距離在 160Km左右,徑向庫(kù)長(zhǎng)為 250m。采用VCP21模式 (6min完成 9個(gè)仰角的掃描：0.5°,1.45°,2.4°,3.4°,4.3°,6.0°,9.9°,14.6°,19.5°)進(jìn)行掃描,如圖 1所示。

本方法在垂直方向取 20Km,分辨率為 0.5Km,水平方向取 300Km,分辨率為 1Km的三維數(shù)據(jù)矩陣Data(40,300,300),這樣可得到 40層不同高度平面上的 CAPPI數(shù)據(jù),此 CAPPI數(shù)據(jù),實(shí)際上是由不同仰角不同距離上的雷達(dá)回波觀測(cè)值通過內(nèi)插得到的。利用同一高度上的 CAPPI數(shù)據(jù),可以較方便地分析雷達(dá)回波信息在某高度上的水平分布,便于與鄰近該高度的天氣圖分析相結(jié)合;而用不同高度上的 CAPPI數(shù)據(jù)還可以了解風(fēng)暴的三維結(jié)構(gòu)。

不同高度上的 CAPPI數(shù)據(jù)可通過網(wǎng)格反插得到,網(wǎng)格反查是以雷達(dá)站為網(wǎng)格中點(diǎn)。利用數(shù)學(xué)公式得到網(wǎng)格中任意網(wǎng)格點(diǎn) D到雷達(dá)站的距離 R及與正北方向的夾角β。這樣,D網(wǎng)格點(diǎn)的分析值就近似等于方位角為β,徑向距離為 R的點(diǎn)的觀測(cè)值。

對(duì)于三維網(wǎng)格點(diǎn),如圖 2,D網(wǎng)格點(diǎn)到雷達(dá)站高度為 H、距離為 R、與雷達(dá)站水平面的夾角為θ、方位角為β,D網(wǎng)格點(diǎn)的雷達(dá)回波分析值按照最鄰近插值法及內(nèi)插平均有以下 3種情況：

①θ ＜0.5°或 θ＞19.5°(第 1個(gè)和最后 1個(gè)仰角)

③α＜θ＜α′(α′為α仰角的下一掃描仰角 )

式中 f(H,R,θ,β)代表 D網(wǎng)格點(diǎn)的雷達(dá)回波分析值,f′(R,αβ)代表仰角為α,方位角為β,徑向距離為 R處的雷達(dá)回波觀測(cè)值。α為 VCP21掃描模式中的體掃仰角。經(jīng)過循環(huán)計(jì)算可求出三維網(wǎng)格數(shù)據(jù) Data(40,300,300)。

2.2 實(shí)例對(duì)比分析

圖 3、圖 4分別是 2008-07-28T 18：18貴陽(yáng)雷達(dá) 2Km及 5Km高度雷達(dá)反射率因子 CAPPI圖像,圖 (a)為本方法得出,圖 (b)是 784廠雷達(dá)顯示軟件計(jì)算的,從圖中可以看出,雷達(dá)反射率因子的強(qiáng)度及位置比較一致,本文還經(jīng)過大量的對(duì)比分析,表明本方法得到的結(jié)果是準(zhǔn)確的,算法是可行的。

2.3 反射率因子垂直剖面計(jì)算

通過對(duì)雷達(dá)反射率因子處理后得到三維網(wǎng)格數(shù)據(jù) Data(40,300,300),網(wǎng)格中心的經(jīng)緯度為雷達(dá)站經(jīng)緯度,如圖 5,在屏幕上從 A到 B點(diǎn)畫任意一直線,通過屏幕坐標(biāo)轉(zhuǎn)換得到 A和 B點(diǎn)的經(jīng)緯度,而網(wǎng)格點(diǎn)的邊界經(jīng)緯度是已知的,通過計(jì)算,A、B點(diǎn)的雷達(dá)反射率因子大小為 Data(40,x1,y1)、Data(40,x2,y2),進(jìn)而可求 AB直線上任意 C點(diǎn)的雷達(dá)反射率因子 Data(40,x,y)。

x1,y1,x2,y2分別為 A、B點(diǎn)與網(wǎng)格點(diǎn)邊界的網(wǎng)格距離,dis為 A、B兩點(diǎn)之間的網(wǎng)格距離,通過相關(guān)數(shù)學(xué)計(jì)算可得沿 AB直線切面的網(wǎng)格數(shù)據(jù) PouMianDbz[40,dis]。

3 雙線性內(nèi)插值法

最簡(jiǎn)單的插值算法是最鄰近插值法,該方法計(jì)算量小,速度快,在多數(shù)情況下能得到令人滿意的結(jié)果,但是當(dāng)對(duì)圖像進(jìn)行多倍放大處理時(shí),最鄰近算法往往滿足不了放大的要求。雙線性內(nèi)插值算法計(jì)算量大,但計(jì)算結(jié)果更符合實(shí)際情況,成像后圖像質(zhì)量高,不會(huì)出現(xiàn)像素值不連續(xù)的的情況。

對(duì)于一個(gè)新網(wǎng)格點(diǎn) NewGridData(i+a,j+b),其中 i、j均為非負(fù)整數(shù),a、b為 [0,1]區(qū)間的浮點(diǎn)數(shù),則這個(gè)網(wǎng)格點(diǎn)的值NewGridData(i+a,j+b)可由原網(wǎng)格中 OldGridData(i,j)、OldGridData(i+l,j)、Old-GridData(i,j+1)、OldGridData(i+l,j+1)4個(gè)網(wǎng)格點(diǎn)的值決定,即：NewGridData(i+a,j+b)=(1-a)×(1-b)×OldGridData(i,j)+(1-a)×b×Old-GridData(i,j+1)+a ×(1-b)×OldGridData(i+l,j)+a ×b ×OldGridData(i+1,j+1),依次類推 ,求出所有新網(wǎng)格點(diǎn)數(shù)值,這就是雙線性內(nèi)插值法。

3.1 主要 C#代碼

由于篇幅原因,在這里省略了部分代碼：

{

//zoomW為 X軸方向的放大倍數(shù),zoomH為 Y軸方向的放大倍數(shù)

通過以上計(jì)算得到新的網(wǎng)格數(shù)據(jù) NewGridData,利用此網(wǎng)格數(shù)據(jù)即可進(jìn)行成像處理。

3.2 經(jīng)雙線性內(nèi)插值法與未經(jīng)插值得到的圖像對(duì)比

圖 6是 2008-07-28 T 18：18貴陽(yáng)雷達(dá)組合反射率因子圖像,圖 7是圖 6中 A點(diǎn)向 B點(diǎn)作垂直剖面得到的圖像,其中 (a)為未經(jīng)插值處理得到的圖像,圖像小,視覺效果差,(c)是 (a)通過看圖軟件放大后的圖像,出現(xiàn)模糊現(xiàn)象,(b)為經(jīng)雙線性內(nèi)插值后得到的圖像,清晰度高,視覺效果好。

4 小結(jié)

本文介紹了基于 CD型號(hào)雷達(dá)反射率因子垂直剖面的實(shí)現(xiàn)方法以及利用雙線性內(nèi)插值法應(yīng)用到垂直剖面圖像向中。通過比較說明本方法中將雷達(dá)原始數(shù)據(jù)轉(zhuǎn)換成三維網(wǎng)格數(shù)據(jù)、任意方向作垂直剖面分析方法是可行的。最后利用雙線內(nèi)插值法實(shí)現(xiàn)圖像縮放的算法,結(jié)果表明經(jīng)雙線性內(nèi)插值法放大后的圖像質(zhì)量較高,解決了因垂直剖面圖像較小不利于顯示及放大的問題。

[1] 俞小鼎,姚秀萍,等 .多普勒天氣雷達(dá)原理與業(yè)務(wù)應(yīng)用[M].北京：氣象出版社,2006.

[2] 黃云仙,張英 .多普勒天氣雷達(dá)數(shù)據(jù)插值方法比較研究[J].遙感信息 ,2008,2：39-45.

[3] 楊毅,等 .多普勒雷達(dá)資料格點(diǎn)化方案的比較研究[J].干旱氣象,2004,22(2)：6-10.

[4] 肖艷姣,劉黎平 .新一代天氣雷達(dá)網(wǎng)資料的三維格點(diǎn)化及拼圖方法研究[J].氣象學(xué)報(bào),2006,64(5)：648-656.

[5] 羅斌,越飛,等 .VisualC#2005編程實(shí)例精粹[M].北京：中國(guó)水利水電出版社,2006.

Vertical profile and dual-linear interpolation’s application of CD radar in its image

XiongWei1,ZhouMing-fei1,Pan Xu-yan2

(1.GuizhouMeteorologicalObservatory,Guiyang 550002;2.The Climate centre of Guizhou,Guiyang 550002)

The New Generation DopplerWeatherRadarReflectivity factorof Guizhou is transfor med to the coordinate transfor mation into three-dimensional grid data post-processing byC#programming.Based on bilinear interpolation method,Any direction vertical section of radar reflectivity factor is achieved.Combining Guiyang radar analysis,the results showed that the bilinear interpolation method of high quality image amplifier and the method solve for vertical section is s maller and enlarge images.

CD radar;vertical section;quadratic interpolation;image

TP31

A

1003-6598(2010)03-0007-04

2010-01-11

熊偉 (1978-),男,工程師,主要從事短時(shí)臨近天氣預(yù)報(bào)工作。

貴州省科技廳黔科合重大專項(xiàng)字 (2006)6034號(hào)。