• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于改進(jìn)光流法的雷達(dá)圖像運(yùn)動估計(jì)

    2017-12-20 10:06:22王志斌肖艷姣
    關(guān)鍵詞:光流法光流流場

    王志斌,肖艷姣,吳 濤

    (1.中國氣象局 武漢暴雨研究所,湖北 武漢 430205;2.武漢中心氣象臺,湖北 武漢 430074)

    基于改進(jìn)光流法的雷達(dá)圖像運(yùn)動估計(jì)

    王志斌1,肖艷姣1,吳 濤2

    (1.中國氣象局 武漢暴雨研究所,湖北 武漢 430205;2.武漢中心氣象臺,湖北 武漢 430074)

    介紹了一種改進(jìn)的光流方法對天氣雷達(dá)圖像運(yùn)動進(jìn)行估計(jì),計(jì)算天氣雷達(dá)圖像運(yùn)動的光流變化,從而得到矢量場。用傳統(tǒng)的Horn & Schunck方法中的全局平滑化算子估計(jì)時,其對數(shù)據(jù)的噪聲敏感,且難以達(dá)到全場最優(yōu)條件。通過加入改進(jìn)后的高價平滑算子,可以減少靠向二級最小值的可能。另外,通過加入Lucas & Kanade方法計(jì)算局部光流得到小范圍的光流場,在小的區(qū)域內(nèi)光流場容易滿足最優(yōu)條件。適當(dāng)運(yùn)用平滑處理得到整場的矢量場,并以它作為改進(jìn)Horn & Schunck方法的初猜場進(jìn)行計(jì)算,且應(yīng)用二維連續(xù)無輻散質(zhì)量方程約束,同時對得到的風(fēng)場進(jìn)行風(fēng)速以及風(fēng)向的質(zhì)量控制,得到比較完整的雷達(dá)運(yùn)動矢量場。通過實(shí)驗(yàn)分析發(fā)現(xiàn),改進(jìn)的光流法優(yōu)于氣象上傳統(tǒng)的交叉相關(guān)法,且優(yōu)于單一使用某種光流方法的結(jié)果。

    雷達(dá)圖像;運(yùn)動估計(jì);光流法;約束;改進(jìn)

    1 概 述

    多普勒天氣雷達(dá)圖像中包含了豐富的天氣現(xiàn)象,通過圖像可以對大風(fēng)、冰雹、雷電、暴雨等強(qiáng)天氣進(jìn)行監(jiān)測,也可通過多幅圖像的變化對它們進(jìn)行預(yù)報(bào)[1]。目前氣象預(yù)報(bào)業(yè)務(wù)上使用自動外推預(yù)報(bào)技術(shù),包括兩類:交叉相關(guān)法[2-3]和單體質(zhì)心法[4-7]。單體質(zhì)心法是將雷達(dá)的三維數(shù)據(jù)進(jìn)行分析;交叉相關(guān)法利用不同時次天氣雷達(dá)圖像求最優(yōu)空間的相關(guān),選取不同時次相關(guān)最好的空間匹配塊,來確定圖像的移動矢量特征,并進(jìn)行圖像外推預(yù)報(bào)。交叉相關(guān)算法計(jì)算簡單,可以求不同時刻的相關(guān)系數(shù),也可以求歐氏距離。但對形狀變化較快的天氣雷達(dá)圖像,通過交叉相關(guān)法計(jì)算出的運(yùn)動矢量場質(zhì)量會降低[8]。

    為了克服這些缺點(diǎn),引入了在計(jì)算機(jī)領(lǐng)域運(yùn)用較多的光流方法。采用改進(jìn)的光流法計(jì)算得到天氣雷達(dá)圖像的光流場來代替交叉相關(guān)法得到的風(fēng)場(后面提到的光流場和風(fēng)場及矢量場是等價的)。光流法最早由Gibson[9]提出,Horn & Schunck[10]提出了一種基于全局約束的求解方法(簡稱HS方法)。Lucas-Kanade[11]提出了一種局部的約束求解方法(簡稱LK方法)?;诠饬鞯幕炯s束方程,還產(chǎn)生了其他光流方法,但大多數(shù)可以用變分模型來描述,即最小化某個能量泛函的過程。這個泛函可以由兩項(xiàng)組成,即平滑項(xiàng)和數(shù)據(jù)項(xiàng)。數(shù)據(jù)項(xiàng)和圖像本身的數(shù)據(jù)有關(guān),如雷達(dá)回波數(shù)據(jù)、圖像的灰度等,數(shù)據(jù)項(xiàng)決定了運(yùn)動的真實(shí)程度。平滑項(xiàng)是泛函的附加項(xiàng),進(jìn)行條件約束,主要解決方程的不適定性。確定了泛函方程后,可以運(yùn)用歐拉-拉格朗日方程進(jìn)行求解,其數(shù)值解一般采用松弛迭代或最優(yōu)下降算法求出。光流法在計(jì)算機(jī)圖像處理、模式識別、交通、醫(yī)學(xué)等領(lǐng)域[12-14]應(yīng)用廣泛。韓雷[8]運(yùn)用HS方法進(jìn)行了外推預(yù)報(bào),其不足之處是全場不容易滿足最優(yōu)條件,容易陷入局部最小化,特別是在天氣雷達(dá)圖像的場能量不集中時。曹春燕[15]使用LK方法,運(yùn)用局部最優(yōu)進(jìn)行求解計(jì)算,對局部容易滿足最優(yōu)條件[16-18],但對整場不容易滿足最優(yōu)條件。

    文中將兩者結(jié)合起來,先運(yùn)用LK求局部最優(yōu);再以局部最優(yōu)結(jié)果作為初始風(fēng)場,利用改進(jìn)HS方法計(jì)算全場光流值,得到相對合理的運(yùn)動矢量場;最后通過二維連續(xù)質(zhì)量方程進(jìn)行約束,得到更加完善的運(yùn)動矢量場。與傳統(tǒng)的交叉相關(guān)法和單一的光流法相比較,結(jié)果表明該方法更合理有效,在氣象的短時預(yù)報(bào)業(yè)務(wù)中得到了較好的應(yīng)用。

    2 光流方法改進(jìn)及天氣雷達(dá)數(shù)據(jù)的質(zhì)量控制

    光流方法的原理是,光流是空間運(yùn)動物體在平面上投影時運(yùn)動的速度[19-20],也是研究圖像灰度在時間上的變化,因此可以將光流矢量定義為投影平面特定坐標(biāo)點(diǎn)上的灰度瞬時變化率,代表二維矢量場。二維光流場實(shí)際上是實(shí)際物體的三維運(yùn)動場的投影[21-22]。光流法的實(shí)質(zhì)是由二維光流場重構(gòu)三維運(yùn)動場。下面介紹兩種經(jīng)典的光流方法:HS和LK。

    設(shè)在平面坐標(biāo)上有一點(diǎn)M(x,y),它代表了三維空間中一點(diǎn)V(x,y,z)在平面坐標(biāo)中的投影。該點(diǎn)在t時刻的灰度為I(x,y,t)。另假設(shè)該點(diǎn)在t+Δt時刻運(yùn)動到(x+Δx,y+Δy),在一定時間間隔Δt內(nèi)灰度值保存不變。即

    I(x,y,t)=I(x+Δx,y+Δy,t+Δt)

    (1)

    將式(1)按照泰勒公式展開:

    (2)

    其中,Rn包含了Δx,Δy,Δt的二次以上的無窮階小項(xiàng)。

    變換式(2)即為:

    (3)

    式(3)除以Δt,并且取Δt→0的極限,忽略掉高階項(xiàng)Rn可得:

    (4)

    可以簡寫為:

    Ixu+Iyv+It=0

    (5)

    式(5)就是光流方程。Ix,Iy為圖像灰度的空間梯度;It為圖像灰度的時間變化率;u,v為光流場。在方程中,Ix,Iy,It可以通過相鄰的兩幅圖像的灰度值計(jì)算出來,但式中有兩個變量,而只有一個方程,如果不做任何假設(shè),方程式無法求解。因此,需要引入更多的約束條件,才能求解u,v。

    2.1 基于全局的HS算法及其改進(jìn)

    如果相隔時間很短的兩幅圖像且灰度變化也很小,在光照不變的條件下,Horn與Schunk等導(dǎo)出了光流的基本方程。巧妙地將二維運(yùn)動和圖像灰度相關(guān)聯(lián)。這個方程有兩個假設(shè),灰度守恒和平滑約束,并且提出了全局約束條件方程。通過變分方法可以求解該方程。全局約束方程如下:

    J=J0+α2JHS

    (6)

    其中,α為平滑系數(shù);J0為灰度守恒項(xiàng);JHS為平滑約束項(xiàng)。

    J0=?(Ixu+Iyv+II)2dxdy

    (7)

    (8)

    根據(jù)HS約束方程,用歐拉-拉格朗日方程求解,可得:

    Ix(Ixu+Iyv+It)-α2Δu=0

    (9)

    Iy(Ixu+Iyv+It)-α2Δv=0

    (10)

    其中,Δ是二階的拉普拉斯算子。

    (11)

    用雅克比迭代法可以求解式(9)和式(10)。

    (12)

    (13)

    用上述方法求解光流場,由于采用全局優(yōu)化方案,很容易陷入局部最小值,而且對圖像噪聲敏感。因此改平滑項(xiàng)為其他方法,即改用Wahba and Wendelberger[23]方法(簡稱Jww)。這可以減少靠向二級最小值的可能。即JHS改為Jww:

    (14)

    同樣可以用高階歐拉-拉格朗日方程求解。

    Ix(Ixu+Iyv+It)-α2Δ2u=0

    (15)

    Iy(Ixu+Iyv+It)-α2Δ2v=0

    (16)

    其中,Δ2u,Δ2v為:

    (17)

    (18)

    同樣可以用雅克比迭代法求解。方程如下:

    (19)

    (20)

    圖1中描述了各格點(diǎn)。其中f0表示為ui,j,i,j為對應(yīng)的格點(diǎn)坐標(biāo)。

    i-2i-1ii+1i+2j+2f12j+1f8f4f5jf11f3f0f1f9j-1f7f2f6j-2f10

    圖1 格點(diǎn)表示

    2.2 LK方法

    在HS方法中,采用Jww來改進(jìn)數(shù)據(jù)平滑項(xiàng),但仍然有可能不是全局最優(yōu),特別是數(shù)據(jù)噪聲大時。因此LK提出了局部最優(yōu)化的方法。基本思路是把圖像分為多個小區(qū)域,因?yàn)樵谛》秶鷥?nèi)容易滿足光流條件,抗噪聲也強(qiáng)。

    假設(shè)一個小區(qū)域范圍內(nèi)的光流恒定,根據(jù)光流方程有:

    (21)

    式中有兩個未知變量u,v,但方程多余兩個,是一個超定方程,需要用最小二乘法求解。

    將式(15)記為如下行列式:

    (22)

    (23)

    (24)

    得到u,v的解為:

    (25)

    如果在選擇的小區(qū)域范圍內(nèi)使用一個窗口權(quán)重函數(shù),使得靠近中心區(qū)域的影響比其他范圍大,得到的速度場為:

    (26)

    其中,W為權(quán)重函數(shù),用最小二乘法求解可得:

    (27)

    2.3 HS-LW混合方法求解

    LW方法對一個小區(qū)域內(nèi)進(jìn)行光流計(jì)算,HS是對整場進(jìn)行最優(yōu)計(jì)算,將兩者結(jié)合起來,先用LW計(jì)算出全場的光流場,再用這個光流場作為HS方法的初始值并用HS方法計(jì)算出全場的光流。計(jì)算出的光流場相對合理。

    利用HS和LW方法進(jìn)行灰度導(dǎo)數(shù)計(jì)算時,用Sobel算子進(jìn)行計(jì)算。Sobel算子有對噪聲處理的魯棒性。LK方法中的窗口函數(shù)為楊輝三角系數(shù),一般取n=9。

    2.4 無輻散質(zhì)量方程約束

    最后,用二維連續(xù)無輻散質(zhì)量方程對HS-LW方法求出的矢量場進(jìn)行約束。

    二維連續(xù)無輻散質(zhì)量方程[24]為:

    (28)

    同樣用變分的方法,對HS-LW求出的矢量場u0,v0和實(shí)際的矢量場u,v用一個代價函數(shù)表示,并最終求解出實(shí)際矢量場。

    代價函數(shù)為:

    (29)

    根據(jù)歐拉-拉格朗日方程求解得:

    (30)

    (31)

    邊界條件為λ(x,y)=0。

    將式(30)和式(31)代入式(28)得到:

    (32)

    式(32)可以通過有限差分方法計(jì)算出λ。同理也可以通過式(30)和式(31)計(jì)算出u,v。

    2.5 數(shù)據(jù)處理質(zhì)量控制及程序計(jì)算流程

    文中使用的資料為湖北省內(nèi)的6部多普勒天氣雷達(dá)資料。在進(jìn)行拼圖前,對單部天氣雷達(dá)圖像數(shù)據(jù)進(jìn)行了質(zhì)量控制,方法采用文獻(xiàn)[25]提出的算法進(jìn)行。首先針對不同雷達(dá)型號對數(shù)據(jù)進(jìn)行規(guī)整和缺測填補(bǔ),對原始數(shù)據(jù)進(jìn)行有效性檢查,對孤立回波進(jìn)行過濾。然后選取了五個回波特征量作為候選因子,對超折射回波進(jìn)行抑制,計(jì)算出超折射回波和正?;夭ǖ母怕史植紙D,找出兩者的差別,建立各因子的隸屬函數(shù)。最后通過模糊邏輯算法,把這些因子進(jìn)行綜合運(yùn)算,得到控制值,大于某值則去掉該回波,并用周圍正常回波進(jìn)行替代,否則該值正常不進(jìn)行處理。

    進(jìn)行單站質(zhì)量控制后,把極坐標(biāo)的單站雷達(dá)資料轉(zhuǎn)換成三維的直角坐標(biāo)格式,采用文獻(xiàn)[26]提出的方法進(jìn)行。取3 km高度的CAPPI拼圖作為光流場的輸入灰度值。在進(jìn)行了單站質(zhì)量控制的基礎(chǔ)上,為了進(jìn)一步減少數(shù)據(jù)噪聲,仍用9點(diǎn)平滑方法對3 km的CAPPI場進(jìn)行質(zhì)量控制。

    通過上述各種光流方法計(jì)算出的光流場,有些格點(diǎn)的值仍有不合理的地方,需要進(jìn)一步加上附加條件進(jìn)行約束和控制。首先對計(jì)算出的光流場某些格點(diǎn)值大于50 m/s的值進(jìn)行剔除。將某一格點(diǎn)和其周圍格點(diǎn)(一般是9*9格點(diǎn)范圍)的平均值進(jìn)行比較。當(dāng)超過平均速度6 m/s,以及夾角超過30°,這個值也剔除,并用9*9格點(diǎn)范圍內(nèi)的平均值進(jìn)行替換,同時這個值和整場的平均值偏差大于15 m/s時給予剔除。

    3 實(shí) 驗(yàn)

    為了評價各種光流方法的效果,把各種光流算法和TREC計(jì)算出的風(fēng)場對天氣雷達(dá)圖像運(yùn)用歐拉后向外推算法進(jìn)行外推預(yù)報(bào),并用氣象上三個通用指標(biāo)進(jìn)行判斷。三個指標(biāo)分別為擊中率(POD)、虛警率(FAR)和臨界成功指數(shù)(CSI)。

    (33)

    (34)

    (35)

    將預(yù)報(bào)結(jié)果和實(shí)況位置的雷達(dá)圖像數(shù)據(jù)進(jìn)行比較,預(yù)報(bào)出現(xiàn)而實(shí)況也出現(xiàn)為成功;實(shí)況出現(xiàn)但預(yù)報(bào)沒有出現(xiàn)為失?。活A(yù)報(bào)出現(xiàn)但實(shí)況沒有出現(xiàn)為虛假。外推時間為30 min和60 min。把雷達(dá)圖像數(shù)據(jù)數(shù)值分為5個等級(5~15 dbz,15~25 dbz,25~35 dbz,35~45 dbz,45 dbz)進(jìn)行評分。選取3個不同的影響湖北省的降水天氣過程,分別是2016年4月15日到16日的層狀云降水,2016年6月30日到7月1日以及2016年7月3日到7月4日對流性降水過程。

    從表1和表2可以看出,HS+LW+二維約束的效果略好,無論是擊中率還是成功指數(shù)都比其他方法好,虛警率也較低,其外推30 min的成功指數(shù)比TREC高5個百分點(diǎn)。其他光流方法和TREC的效果相當(dāng),指標(biāo)表現(xiàn)上各有千秋,外推60 min的效果明顯不如外推30 min的,說明預(yù)報(bào)時間越長,效果越差。

    表1 2016年4月15日-16日30 min外推結(jié)果平均

    表2 2016年4月15日-16日60 min外推結(jié)果平均

    從表3和表4可以看出,HS+LW+二維約束在擊中率、成功指數(shù)以及虛警率方面比其他方法好,其30 min擊中率比TREC高9個百分點(diǎn),成功指數(shù)高6個百分點(diǎn)。除HS方法和TREC相當(dāng)外,其他光流方法比TREC效果略好。

    表3 2016年6月30日-7月1日30 min外推結(jié)果平均

    表4 2016年6月30日-7月1日60 min外推結(jié)果平均

    從表5和表6可以看出,HS+LW+二維約束在擊中率、成功指數(shù)以及虛警率方面都比其他方法略好,其30 min擊中率比HS高8個百分點(diǎn),成功指數(shù)高7個百分點(diǎn)。除HS方法和TREC相當(dāng)外,其他光流方法比TREC效果略好 。

    表5 2016年7月3日-7月4日30 min外推結(jié)果平均

    表6 2016年7月3日-7月4日60 min外推結(jié)果平均

    天氣雷達(dá)圖像光流估計(jì)處理流程見圖2。

    圖2 雷達(dá)圖像光流計(jì)算處理流程

    4 結(jié)束語

    介紹了光流方法的基本原理,提出了一種改進(jìn)的HS方案。結(jié)合HS及LW各自優(yōu)點(diǎn)聯(lián)合計(jì)算出光流場,通過二維連續(xù)方程進(jìn)行約束。通過天氣雷達(dá)圖像的變化計(jì)算出實(shí)用的風(fēng)場,并利用外推方法對各種結(jié)果進(jìn)行了比較。

    改進(jìn)方法在對流性降水時優(yōu)于TREC,而在層狀云降水時差別不明顯,只是HS+LW+二維約束方法略好于其他方法。通過改進(jìn)的HS方法,進(jìn)一步減少了計(jì)算結(jié)果趨向局部最小的可能性,結(jié)果優(yōu)于HS。通過使用LW方法產(chǎn)生的初始光流的結(jié)果,作為改進(jìn)的HS方法的初始場,使得原來使用單一的LW方法的結(jié)果得到了進(jìn)一步優(yōu)化。通過結(jié)合改進(jìn)的HS、LW和二維連續(xù)方程約束等方法,得到了效果較好的光流結(jié)果。

    與傳統(tǒng)的TREC方法相比,光流法側(cè)重于變化,因?yàn)樵谟?jì)算光流場時不僅考慮了天氣雷達(dá)圖像在空間上的變化,同時也考慮了時間上的變化。該中和方法計(jì)算耗時小,在一般微機(jī)上幾秒即可完成,因此有較大的應(yīng)用推廣價值。

    [1] Wong W K.Merging radar nowcast with convection-permitting rapidly-updated NWP model for significant convection forecast[C]//4th international symposium on nowcasting and very-short-range forecast.[s.l.]:[s.n.],2016.

    [2] Rinehart R E,Garvey E T.Three-dimensional storm motion detection by conventional weather radar[J].Nature,1978,273(5660):287-289.

    [3] Wilson J W,Crook N A,Mueller C K,et al.Nowcasting thunderstomes:a status report[J].Bulletin of the American Meteorological Society,1998,79(10):2079-2100.

    [4] Crane R K.Automatic cell detection and tracking[J].IEEE Transactions on Geoscience Electronics,1979,17(4):250-262.

    [5] Dixon M, Wiener G. TITAN: Thunderstorm identification, tracking,analysis,and nowcasting:a radar-based methodology[J].Journal of Atmospheric & Oceanic Technology,1993,10(6):785-797.

    [6] Rasmussen R,Dixon M,Hage F,et al.Weather support to deicing decision making (WSDDM):a winter weather nowcasting system[J].Bulletin of the American Meteorological Society,2001,82(4):579-595.

    [7] Mueller C,Saxen T,RobertsR,et al.NCAR auto-nowcast system[J].Weather & Forecasting,2003,18(4):545-561.

    [8] 韓 雷,王洪慶,林隱靜.光流法在強(qiáng)對流天氣臨近預(yù)報(bào)中的應(yīng)用[J].北京大學(xué)學(xué)報(bào):自然科學(xué)版,2008,44(5):751-755.

    [9] Gibson J J.The perception of the visual world[M].[s.l.]:[s.n.],1950.

    [10] Horn B K P,Schunck B G.Determining optical flow[J].Artificial Intelligence,1981,17(1-3):185-204.

    [11] Lucas B D,Kanade T.An iterative image registration technique with an application to stereo vision[C]//Proceedings of seventh international joint conference on artificial intelligence.Vancouver:IEEE,1981.

    [12] 史金龍,白素琴,施金宛,等.基于CLG光流算法的云的運(yùn)動分析與研究[J].計(jì)算機(jī)技術(shù)與發(fā)展,2011,21(12):135-138.

    [13] 馬 龍,王魯平,李 飚,等.自然場景圖像的光流場估計(jì)[J].系統(tǒng)工程與電子技術(shù),2012,34(6):1278-1282.

    [14] 柳士俊,張 蕾.光流法及其在氣象領(lǐng)域里的應(yīng)用[J].氣象科技進(jìn)展,2015,5(4):16-21.

    [15] 曹春燕,陳元昭,劉東華,等.光流法及其在臨近預(yù)報(bào)中的應(yīng)用[J].氣象學(xué)報(bào),2015,73(3):471-480.

    [16] 楊國亮,王志良,牟世堂,等.一種改進(jìn)的光流算法[J].計(jì)算機(jī)工程,2006,32(15):187-188.

    [17] 陳曉靜,敬忠良,張 軍.基于衛(wèi)星圖像數(shù)據(jù)的大范圍風(fēng)場重建[J].計(jì)算機(jī)工程,2013,39(12):233-236.

    [18] 魏國劍,侯志強(qiáng),李 武,等.融合光流檢測與模板匹配的目標(biāo)跟蹤算法[J].計(jì)算機(jī)應(yīng)用研究,2014,31(11):3498-3501.

    [19] 鄭來芳,孫 煒,歐陽明華,等.結(jié)合光流和人工勢場的風(fēng)管機(jī)器人避障方法[J].計(jì)算機(jī)工程與應(yīng)用,2016,52(9):243-247.

    [20] 吳振杰,毛曉波.一種改進(jìn)的幀間差光流場算法[J].計(jì)算機(jī)工程與應(yīng)用,2013,49(18):200-203.

    [21] 陳 震,高滿屯,曲仕茹,等.基于直線光流場的三維運(yùn)動和結(jié)構(gòu)重建[J].計(jì)算機(jī)學(xué)報(bào),2004,27(8):1046-1055.

    [22] 宋 爽,楊 健,王涌天.全局光流場估計(jì)技術(shù)及展望[J].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào),2014,26(5):841-850.

    [23] Wahba G,Wendelberger J.Some new mathematical methods for variational objective analysis using splines and cross validation[J].Monthly Weather Review,1980,108:1122-1143.

    [24] Li L,Schmid W,Joss J.Nowcasting of motion and growth of precipitation with radar over a complex orography[J].Journal of Applied Meteorology,1995,34:1286-1300.

    [25] 吳 濤,萬玉發(fā),沃偉鋒,等.SWAN系統(tǒng)中雷達(dá)反射率因子質(zhì)量控制算法及其應(yīng)用[J].氣象科技,2013,41(5):809-817.

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

    MotionEstimationforRadarImageBasedonImprovedOpticalFlowMethod

    WANG Zhi-bin1,XIAO Yan-jiao1,WU Tao2

    (1.Institute of Wuhan Heavy Rain,China Meteorological Administration,Wuhan 430205,China;2.Wuhan Central Meteorological Observatory,Wuhan 430074,China)

    An improved optical flow method is introduced to estimate the radar images motion and calculate the optical flow change of motion for radar image,obtaining the vector field.It is sensitive to data noise and difficult to meet the optimal condition when global smoothing operator from traditional Horn & Schunck method is used to estimate.After adding the improved smoothing operator with high level,it can reduce the possibility for the secondary minimum.In addition,through Lucas & Kanade method for calculation of the local optical flow,a small range of optical flow field is obtained,which can meet the optimal requirements.Smooth processing is applied to get the entire vector field appropriately,with it as first guess filed of improved Horn & Schunck method for calculation.At the same time,through the two-dimensional continuous quality equation constraints without radiation,controlling the wind speed and direction of wind field acquired,the more complete motion vector field for radar is obtained.Experimental analysis shows that the improved optical flow method is better than traditional cross-correlation method and a single optical flow method.

    radar images;motion estimation;optical flow method;constraints;improvement

    TP39

    A

    1673-629X(2017)12-0170-06

    10.3969/j.issn.1673-629X.2017.12.037

    2016-10-20

    2017-02-23 < class="emphasis_bold">網(wǎng)絡(luò)出版時間

    時間:2017-08-01

    國家自然科學(xué)基金資助項(xiàng)目(41275106);科技部行業(yè)專項(xiàng)公益性行業(yè)(氣象)科研專項(xiàng)(GYHY201306008);湖北省雷電專項(xiàng)( FL-Z-201401)

    王志斌(1965-),男,正研高工,研究方向?yàn)橛?jì)算機(jī)圖像處理、數(shù)據(jù)庫開發(fā)。

    http://kns.cnki.net/kcms/detail/61.1450.TP.20170801.1550.020.html

    猜你喜歡
    光流法光流流場
    利用掩膜和單應(yīng)矩陣提高LK光流追蹤效果
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場計(jì)算
    基于物理學(xué)的改善粒子圖像測速穩(wěn)健光流方法研究
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計(jì)分析
    基于背景分類的監(jiān)控視頻中的運(yùn)動目標(biāo)檢測算法綜述
    Matlab下視頻處理系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)
    權(quán)重系數(shù)自適應(yīng)光流法運(yùn)動目標(biāo)檢測
    基于瞬態(tài)流場計(jì)算的滑動軸承靜平衡位置求解
    初始段安控圖像目標(biāo)識別方法研究
    国产毛片a区久久久久| 欧美极品一区二区三区四区| aaaaa片日本免费| 久久人妻av系列| 麻豆一二三区av精品| 激情在线观看视频在线高清| 黄色日韩在线| 在线a可以看的网站| 成年女人看的毛片在线观看| 欧美一区二区精品小视频在线| 又紧又爽又黄一区二区| 搡老熟女国产l中国老女人| av专区在线播放| 国产激情偷乱视频一区二区| 小蜜桃在线观看免费完整版高清| 国产激情偷乱视频一区二区| 91狼人影院| 搞女人的毛片| 757午夜福利合集在线观看| 欧美成人性av电影在线观看| 少妇人妻精品综合一区二区 | 97人妻精品一区二区三区麻豆| 国产一区二区亚洲精品在线观看| 午夜激情福利司机影院| 一个人看视频在线观看www免费| 国产乱人视频| 最近视频中文字幕2019在线8| 简卡轻食公司| 午夜视频国产福利| 精品人妻1区二区| 女生性感内裤真人,穿戴方法视频| a在线观看视频网站| 亚洲不卡免费看| 黄色配什么色好看| 国产伦精品一区二区三区视频9| 日日干狠狠操夜夜爽| 成年免费大片在线观看| 亚洲国产精品合色在线| 国产主播在线观看一区二区| a在线观看视频网站| 午夜久久久久精精品| 三级男女做爰猛烈吃奶摸视频| 久久精品久久久久久噜噜老黄 | 精品久久久久久成人av| 亚洲人成网站在线播放欧美日韩| 国产伦精品一区二区三区视频9| 夜夜爽天天搞| 直男gayav资源| 国产精品1区2区在线观看.| 国产不卡一卡二| av视频在线观看入口| 婷婷丁香在线五月| 在线a可以看的网站| 国产精品久久久久久精品电影| 亚洲精品影视一区二区三区av| 成人亚洲精品av一区二区| 亚洲最大成人中文| 丰满的人妻完整版| 一级作爱视频免费观看| 在线播放无遮挡| 此物有八面人人有两片| 观看美女的网站| 女人被狂操c到高潮| 欧美乱色亚洲激情| 久久精品久久久久久噜噜老黄 | 国产 一区 欧美 日韩| 在线观看66精品国产| 国产精品久久久久久精品电影| 国产精品久久视频播放| 免费观看精品视频网站| 亚洲精品一区av在线观看| 最近视频中文字幕2019在线8| 两人在一起打扑克的视频| 9191精品国产免费久久| avwww免费| 51午夜福利影视在线观看| 免费看日本二区| 伊人久久精品亚洲午夜| 国产伦精品一区二区三区视频9| 欧美成人一区二区免费高清观看| 免费观看精品视频网站| 午夜免费成人在线视频| 99国产精品一区二区蜜桃av| 国产色爽女视频免费观看| 夜夜看夜夜爽夜夜摸| 亚洲18禁久久av| 精品一区二区三区视频在线| 搞女人的毛片| 久久精品国产清高在天天线| 一进一出抽搐动态| 国产成+人综合+亚洲专区| 精品人妻视频免费看| 69人妻影院| 亚洲男人的天堂狠狠| 国产毛片a区久久久久| 日本 av在线| 国产精品久久视频播放| 国产av麻豆久久久久久久| 日本免费一区二区三区高清不卡| 露出奶头的视频| 桃色一区二区三区在线观看| 国产乱人伦免费视频| 久久久久久久久大av| 最后的刺客免费高清国语| 国产在线精品亚洲第一网站| 男人舔奶头视频| 99国产综合亚洲精品| 美女高潮的动态| 少妇高潮的动态图| 亚州av有码| 麻豆成人午夜福利视频| 国产亚洲精品综合一区在线观看| 黄色女人牲交| 国产欧美日韩一区二区精品| 久久久久免费精品人妻一区二区| 午夜福利欧美成人| 久久久色成人| 日韩国内少妇激情av| 3wmmmm亚洲av在线观看| 久久久久九九精品影院| 国产在线男女| 亚洲成人中文字幕在线播放| 午夜精品在线福利| 国产熟女xx| 色噜噜av男人的天堂激情| 亚洲精品久久国产高清桃花| 啦啦啦观看免费观看视频高清| 日韩精品青青久久久久久| 免费大片18禁| 一本精品99久久精品77| av天堂在线播放| 性色av乱码一区二区三区2| 精品99又大又爽又粗少妇毛片 | 亚洲av免费在线观看| 国产成人aa在线观看| 在线播放无遮挡| 少妇被粗大猛烈的视频| 久久久国产成人免费| 日韩欧美在线乱码| 成人三级黄色视频| 毛片女人毛片| 91九色精品人成在线观看| 色噜噜av男人的天堂激情| 午夜免费激情av| 国产伦精品一区二区三区四那| 无遮挡黄片免费观看| 91久久精品国产一区二区成人| 国产av一区在线观看免费| 国产午夜福利久久久久久| 亚洲无线在线观看| 欧美高清性xxxxhd video| 99久久久亚洲精品蜜臀av| 亚洲av电影在线进入| 国产成人欧美在线观看| 亚洲一区二区三区色噜噜| 欧美日韩中文字幕国产精品一区二区三区| 亚洲美女搞黄在线观看 | 精品福利观看| 九九热线精品视视频播放| 色综合欧美亚洲国产小说| 国产老妇女一区| 变态另类丝袜制服| 日韩欧美在线二视频| 日本黄大片高清| 欧美性猛交╳xxx乱大交人| 日韩欧美在线二视频| 免费在线观看影片大全网站| 免费大片18禁| 亚洲在线观看片| 日本与韩国留学比较| 男女视频在线观看网站免费| 99精品在免费线老司机午夜| www.色视频.com| 麻豆成人午夜福利视频| 一个人看的www免费观看视频| 欧美黄色片欧美黄色片| 色av中文字幕| 999久久久精品免费观看国产| av在线天堂中文字幕| 蜜桃久久精品国产亚洲av| 国内精品美女久久久久久| 嫩草影院精品99| www.熟女人妻精品国产| 亚洲美女黄片视频| 亚洲成av人片免费观看| 757午夜福利合集在线观看| 亚洲18禁久久av| 最新在线观看一区二区三区| 身体一侧抽搐| 啪啪无遮挡十八禁网站| 久久久久免费精品人妻一区二区| 国产欧美日韩精品亚洲av| 夜夜夜夜夜久久久久| 亚洲精品一卡2卡三卡4卡5卡| 亚洲内射少妇av| 免费看日本二区| 国产极品精品免费视频能看的| 黄色视频,在线免费观看| 国产成人啪精品午夜网站| 亚洲电影在线观看av| 国内毛片毛片毛片毛片毛片| 国产精品野战在线观看| 男人舔女人下体高潮全视频| 能在线免费观看的黄片| 国产视频一区二区在线看| 欧美日韩中文字幕国产精品一区二区三区| 在线播放国产精品三级| 日本五十路高清| 一区二区三区免费毛片| 亚洲国产日韩欧美精品在线观看| 精品午夜福利在线看| 欧美中文日本在线观看视频| 国产爱豆传媒在线观看| 中亚洲国语对白在线视频| 国产精品野战在线观看| 亚洲,欧美精品.| 首页视频小说图片口味搜索| 欧美+亚洲+日韩+国产| 久久天躁狠狠躁夜夜2o2o| 中文字幕av在线有码专区| 精品一区二区免费观看| 久久久国产成人精品二区| 动漫黄色视频在线观看| 日韩av在线大香蕉| 精品久久久久久久人妻蜜臀av| 男女做爰动态图高潮gif福利片| 国产人妻一区二区三区在| 婷婷亚洲欧美| 99国产综合亚洲精品| 亚洲人成网站高清观看| 国产蜜桃级精品一区二区三区| 99久久精品国产亚洲精品| 三级男女做爰猛烈吃奶摸视频| 国产精品自产拍在线观看55亚洲| 国产色婷婷99| 少妇丰满av| 国产精品美女特级片免费视频播放器| 九九在线视频观看精品| 亚洲欧美日韩卡通动漫| 人人妻人人看人人澡| 久久久成人免费电影| 精品无人区乱码1区二区| 国产老妇女一区| 天美传媒精品一区二区| netflix在线观看网站| 色吧在线观看| 日本在线视频免费播放| 日韩中字成人| 首页视频小说图片口味搜索| 国产精品久久久久久亚洲av鲁大| 国产成人影院久久av| 老司机深夜福利视频在线观看| 日本精品一区二区三区蜜桃| 又爽又黄a免费视频| 国产真实伦视频高清在线观看 | 天堂av国产一区二区熟女人妻| 亚洲片人在线观看| 国产单亲对白刺激| 国产男靠女视频免费网站| 精品欧美国产一区二区三| 国产毛片a区久久久久| 综合色av麻豆| 看免费av毛片| 婷婷六月久久综合丁香| 人妻久久中文字幕网| 亚洲专区中文字幕在线| 久久久国产成人免费| 又紧又爽又黄一区二区| 日本 欧美在线| 国产成+人综合+亚洲专区| АⅤ资源中文在线天堂| 九色国产91popny在线| 精品人妻视频免费看| 国产麻豆成人av免费视频| 亚洲中文字幕日韩| 最近最新免费中文字幕在线| 国产白丝娇喘喷水9色精品| 我要看日韩黄色一级片| 国产av麻豆久久久久久久| 国产精品一及| 中文字幕人妻熟人妻熟丝袜美| 精品人妻视频免费看| bbb黄色大片| 亚洲电影在线观看av| 毛片一级片免费看久久久久 | 激情在线观看视频在线高清| 国产精品,欧美在线| 亚洲精品久久国产高清桃花| 国产淫片久久久久久久久 | 国产真实乱freesex| 国产熟女xx| 色噜噜av男人的天堂激情| 日韩高清综合在线| 色尼玛亚洲综合影院| 免费无遮挡裸体视频| 色吧在线观看| 国产人妻一区二区三区在| 麻豆一二三区av精品| 深夜a级毛片| 欧美zozozo另类| 久久久久精品国产欧美久久久| 国产极品精品免费视频能看的| 女生性感内裤真人,穿戴方法视频| 国产日本99.免费观看| 亚洲va日本ⅴa欧美va伊人久久| 美女大奶头视频| 毛片一级片免费看久久久久 | 91麻豆av在线| .国产精品久久| 亚洲内射少妇av| 男人狂女人下面高潮的视频| 亚洲美女黄片视频| 久久久国产成人精品二区| 国产精品久久久久久亚洲av鲁大| 国产中年淑女户外野战色| 日韩中文字幕欧美一区二区| 日本一二三区视频观看| 亚洲国产高清在线一区二区三| 日本 欧美在线| 国产成+人综合+亚洲专区| 午夜福利免费观看在线| 99久久精品一区二区三区| 男女床上黄色一级片免费看| 99热这里只有是精品50| 国产精品一区二区三区四区久久| 国内精品久久久久精免费| 久久久久久久久大av| 国产日本99.免费观看| av天堂中文字幕网| 国产精品影院久久| 2021天堂中文幕一二区在线观| 一区二区三区免费毛片| netflix在线观看网站| 久久久成人免费电影| 午夜激情福利司机影院| 欧美不卡视频在线免费观看| 国产精品久久久久久精品电影| 精品无人区乱码1区二区| 看十八女毛片水多多多| 久久久久久国产a免费观看| 老司机福利观看| 国产精品一区二区性色av| 一级毛片久久久久久久久女| 国产 一区 欧美 日韩| 午夜免费成人在线视频| 国产精品久久视频播放| 国产精品电影一区二区三区| 亚洲欧美日韩无卡精品| 久久久久国内视频| 看黄色毛片网站| 精品国产三级普通话版| 一级作爱视频免费观看| 日日摸夜夜添夜夜添av毛片 | 日韩精品青青久久久久久| 久久6这里有精品| 观看免费一级毛片| 午夜免费成人在线视频| 亚洲精华国产精华精| 国产欧美日韩精品一区二区| 久久精品夜夜夜夜夜久久蜜豆| 美女cb高潮喷水在线观看| 亚洲国产日韩欧美精品在线观看| 国产精品,欧美在线| 又爽又黄无遮挡网站| 日本黄色视频三级网站网址| 日韩欧美免费精品| 尤物成人国产欧美一区二区三区| 人妻夜夜爽99麻豆av| 偷拍熟女少妇极品色| 国产大屁股一区二区在线视频| 亚洲精品一区av在线观看| 俺也久久电影网| 欧美+日韩+精品| 欧美黑人欧美精品刺激| 女人被狂操c到高潮| 国产精品久久久久久久电影| 乱人视频在线观看| 亚洲第一电影网av| 亚洲美女黄片视频| 日韩亚洲欧美综合| 成年版毛片免费区| 国产一区二区三区在线臀色熟女| 级片在线观看| 亚洲国产精品999在线| 国产欧美日韩一区二区精品| 老司机深夜福利视频在线观看| 久久精品影院6| 免费黄网站久久成人精品 | 亚洲人成电影免费在线| 老熟妇乱子伦视频在线观看| 亚洲av电影在线进入| 男女那种视频在线观看| 深夜a级毛片| 国产人妻一区二区三区在| 激情在线观看视频在线高清| 欧美性猛交黑人性爽| 亚洲第一电影网av| 亚洲美女黄片视频| 精品一区二区三区视频在线| 成年版毛片免费区| 欧美日韩中文字幕国产精品一区二区三区| 国产av一区在线观看免费| 日韩欧美精品v在线| 日韩有码中文字幕| 变态另类成人亚洲欧美熟女| 99国产综合亚洲精品| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 久久久色成人| 亚洲av二区三区四区| 国产精品女同一区二区软件 | 在线播放无遮挡| 亚洲精品色激情综合| 女人被狂操c到高潮| 国产高清视频在线播放一区| 色精品久久人妻99蜜桃| 91在线观看av| 精品午夜福利在线看| 精品99又大又爽又粗少妇毛片 | 久久精品夜夜夜夜夜久久蜜豆| 日本 欧美在线| 国产精品永久免费网站| 成人三级黄色视频| 男女视频在线观看网站免费| av黄色大香蕉| 丁香欧美五月| 99热这里只有是精品在线观看 | 能在线免费观看的黄片| 亚洲电影在线观看av| 最近中文字幕高清免费大全6 | 嫩草影院入口| 特级一级黄色大片| 亚洲国产欧美人成| 真实男女啪啪啪动态图| 美女大奶头视频| 国内精品久久久久久久电影| 亚洲一区二区三区色噜噜| 精品久久久久久久末码| 免费看光身美女| 国产三级黄色录像| 真实男女啪啪啪动态图| 亚洲片人在线观看| 亚洲一区高清亚洲精品| 91久久精品电影网| 久久久久国内视频| 亚洲国产精品合色在线| 国产精品自产拍在线观看55亚洲| 高潮久久久久久久久久久不卡| 香蕉av资源在线| 亚洲avbb在线观看| 久久人人精品亚洲av| 99精品久久久久人妻精品| 在线天堂最新版资源| 制服丝袜大香蕉在线| 国产亚洲精品av在线| 国产亚洲欧美在线一区二区| 免费在线观看日本一区| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 一卡2卡三卡四卡精品乱码亚洲| 国产黄片美女视频| 黄色视频,在线免费观看| 丁香六月欧美| 亚洲18禁久久av| 级片在线观看| 久久精品久久久久久噜噜老黄 | 亚洲国产精品成人综合色| 搡老熟女国产l中国老女人| 少妇的逼好多水| 久久人妻av系列| 九九久久精品国产亚洲av麻豆| 亚洲av不卡在线观看| 亚洲国产精品999在线| 久久99热6这里只有精品| 波多野结衣高清无吗| 精品人妻熟女av久视频| 亚洲一区二区三区色噜噜| 18美女黄网站色大片免费观看| 国产毛片a区久久久久| 69人妻影院| 欧美乱色亚洲激情| 免费大片18禁| 亚洲午夜理论影院| 国产国拍精品亚洲av在线观看| 99热这里只有是精品50| 亚洲激情在线av| 亚洲精品乱码久久久v下载方式| 3wmmmm亚洲av在线观看| 国产精品伦人一区二区| 亚洲精品一区av在线观看| 国产麻豆成人av免费视频| 欧美一区二区精品小视频在线| 久久久国产成人免费| 国产一级毛片七仙女欲春2| 一本久久中文字幕| 婷婷精品国产亚洲av| 免费搜索国产男女视频| 又黄又爽又免费观看的视频| 国产伦人伦偷精品视频| 亚洲一区二区三区不卡视频| 欧美性感艳星| 午夜免费男女啪啪视频观看 | 午夜福利在线观看免费完整高清在 | 国产成人福利小说| 夜夜躁狠狠躁天天躁| 国产精品,欧美在线| 一本精品99久久精品77| 91麻豆精品激情在线观看国产| 国产精品永久免费网站| 婷婷六月久久综合丁香| 99热这里只有精品一区| 久久亚洲精品不卡| 亚洲精品乱码久久久v下载方式| 黄色女人牲交| 精品人妻一区二区三区麻豆 | 直男gayav资源| 久久这里只有精品中国| 亚洲欧美精品综合久久99| 精华霜和精华液先用哪个| 午夜福利成人在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲美女搞黄在线观看 | 91狼人影院| 偷拍熟女少妇极品色| 特级一级黄色大片| 搡老妇女老女人老熟妇| 天天一区二区日本电影三级| 国产高清有码在线观看视频| 欧美精品啪啪一区二区三区| 在线观看美女被高潮喷水网站 | 深夜a级毛片| 91麻豆av在线| avwww免费| 高清在线国产一区| 亚洲av.av天堂| 国产在视频线在精品| 亚洲av.av天堂| 男女那种视频在线观看| 深夜精品福利| 99精品在免费线老司机午夜| 亚洲av日韩精品久久久久久密| 免费av观看视频| av在线天堂中文字幕| 亚洲最大成人中文| 亚洲aⅴ乱码一区二区在线播放| 亚洲av二区三区四区| 在线观看午夜福利视频| 熟女电影av网| 给我免费播放毛片高清在线观看| 日本黄色视频三级网站网址| 亚洲欧美日韩高清专用| 国产高清视频在线播放一区| 青草久久国产| 中国美女看黄片| 精品福利观看| 性色av乱码一区二区三区2| 老熟妇乱子伦视频在线观看| 国产精华一区二区三区| 在线观看av片永久免费下载| 又爽又黄无遮挡网站| 日韩欧美一区二区三区在线观看| 亚洲av成人精品一区久久| 国产成年人精品一区二区| 丰满人妻熟妇乱又伦精品不卡| av欧美777| 蜜桃久久精品国产亚洲av| www.色视频.com| 精品久久久久久成人av| 午夜福利欧美成人| 午夜久久久久精精品| 亚洲熟妇熟女久久| 18+在线观看网站| 狠狠狠狠99中文字幕| 欧美区成人在线视频| 国产淫片久久久久久久久 | 最新中文字幕久久久久| 国产亚洲欧美98| 日本五十路高清| 国产精品亚洲av一区麻豆| 日韩中字成人| 一级黄片播放器| 久久精品久久久久久噜噜老黄 | 少妇人妻精品综合一区二区 | 久久99热这里只有精品18| 中文在线观看免费www的网站| 少妇的逼好多水| 小蜜桃在线观看免费完整版高清| 国产高清有码在线观看视频| 女同久久另类99精品国产91| 国产精华一区二区三区| 一进一出抽搐gif免费好疼| 亚洲国产精品成人综合色| 国内精品久久久久久久电影| 日韩人妻高清精品专区| 欧美3d第一页| 欧美在线黄色| av在线老鸭窝| 久久久久久久午夜电影| 少妇熟女aⅴ在线视频| 午夜免费激情av| 好男人电影高清在线观看| 成年女人毛片免费观看观看9| 香蕉av资源在线| 真人做人爱边吃奶动态| 嫩草影院精品99| 网址你懂的国产日韩在线| 国产精华一区二区三区| 久久久久久九九精品二区国产| 美女高潮的动态| 日韩欧美一区二区三区在线观看| 啦啦啦韩国在线观看视频| 乱人视频在线观看| 看十八女毛片水多多多| 真人一进一出gif抽搐免费| 国产男靠女视频免费网站|