• <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)識別方法研究
    亚洲精品国产色婷婷电影| 日日撸夜夜添| 美女高潮的动态| 免费大片18禁| 在线观看一区二区三区| 欧美区成人在线视频| 久久精品久久久久久久性| 卡戴珊不雅视频在线播放| 汤姆久久久久久久影院中文字幕| 亚洲久久久国产精品| 精品一区二区三区视频在线| 国产精品人妻久久久影院| 国内少妇人妻偷人精品xxx网站| 丝袜喷水一区| 国产成人aa在线观看| 国产精品女同一区二区软件| 有码 亚洲区| 国产一区亚洲一区在线观看| 亚州av有码| 成人国产麻豆网| 精品久久久久久久久亚洲| 日本午夜av视频| 欧美精品一区二区免费开放| 永久免费av网站大全| 1000部很黄的大片| 在线天堂最新版资源| 99视频精品全部免费 在线| 亚洲成色77777| h日本视频在线播放| 成人高潮视频无遮挡免费网站| 国产成人freesex在线| 午夜视频国产福利| 国产成人精品一,二区| 国产色婷婷99| 日韩视频在线欧美| 亚洲av不卡在线观看| 99热这里只有是精品在线观看| 亚洲欧美日韩另类电影网站 | 欧美日韩精品成人综合77777| 日韩一本色道免费dvd| 中文字幕av成人在线电影| 中文天堂在线官网| 汤姆久久久久久久影院中文字幕| 中文天堂在线官网| 久久久亚洲精品成人影院| 男人舔奶头视频| 国产乱人视频| 国产精品国产三级专区第一集| 亚洲美女黄色视频免费看| 欧美亚洲 丝袜 人妻 在线| 自拍偷自拍亚洲精品老妇| 日产精品乱码卡一卡2卡三| 免费高清在线观看视频在线观看| 男女边摸边吃奶| 国产爱豆传媒在线观看| 日本黄大片高清| 国产精品久久久久久久电影| 亚洲欧美一区二区三区黑人 | 性色avwww在线观看| 丰满少妇做爰视频| 在线观看免费高清a一片| 亚洲精品国产av成人精品| 热99国产精品久久久久久7| 黄色视频在线播放观看不卡| 下体分泌物呈黄色| 少妇人妻精品综合一区二区| 国产一区有黄有色的免费视频| 色婷婷av一区二区三区视频| 日日啪夜夜撸| 亚洲精品久久久久久婷婷小说| 1000部很黄的大片| 精品人妻偷拍中文字幕| 日韩人妻高清精品专区| 免费高清在线观看视频在线观看| 亚洲欧美日韩无卡精品| 精品久久久久久久久亚洲| 日韩伦理黄色片| 十分钟在线观看高清视频www | 欧美日韩综合久久久久久| 大香蕉久久网| 亚洲中文av在线| 成年美女黄网站色视频大全免费 | 亚洲高清免费不卡视频| 人人妻人人添人人爽欧美一区卜 | 水蜜桃什么品种好| 日本黄色片子视频| 天天躁日日操中文字幕| 国产成人精品福利久久| 啦啦啦中文免费视频观看日本| 久久国产精品男人的天堂亚洲 | 黄色日韩在线| 日韩国内少妇激情av| 亚洲色图av天堂| 夜夜看夜夜爽夜夜摸| 边亲边吃奶的免费视频| 国产精品偷伦视频观看了| 3wmmmm亚洲av在线观看| 亚洲国产欧美在线一区| 久久久精品94久久精品| 天堂8中文在线网| 丰满迷人的少妇在线观看| 嫩草影院入口| av女优亚洲男人天堂| 日韩一区二区三区影片| 亚洲美女搞黄在线观看| av在线老鸭窝| 男人添女人高潮全过程视频| 97超碰精品成人国产| www.av在线官网国产| 国产免费一级a男人的天堂| 亚洲精品久久久久久婷婷小说| 婷婷色综合大香蕉| 亚洲av.av天堂| 久久精品久久久久久噜噜老黄| 日本wwww免费看| av又黄又爽大尺度在线免费看| 黑人猛操日本美女一级片| 一级毛片久久久久久久久女| 国产成人a∨麻豆精品| 国内精品宾馆在线| 久久久久国产精品人妻一区二区| 如何舔出高潮| 亚洲人成网站高清观看| 亚洲精品一二三| 少妇高潮的动态图| 99视频精品全部免费 在线| 男女啪啪激烈高潮av片| 黄片无遮挡物在线观看| 我的女老师完整版在线观看| 国产精品三级大全| 久久精品久久精品一区二区三区| 好男人视频免费观看在线| 久久久久精品性色| 美女视频免费永久观看网站| 又爽又黄a免费视频| a级一级毛片免费在线观看| 亚洲国产精品999| 欧美最新免费一区二区三区| 麻豆成人午夜福利视频| 国精品久久久久久国模美| 日韩免费高清中文字幕av| 国内精品宾馆在线| 99热这里只有是精品在线观看| 精品99又大又爽又粗少妇毛片| 大又大粗又爽又黄少妇毛片口| 伦精品一区二区三区| 中文资源天堂在线| 男女免费视频国产| 日韩一本色道免费dvd| av国产免费在线观看| 久久久久精品久久久久真实原创| 一级毛片 在线播放| 欧美国产精品一级二级三级 | 亚洲av男天堂| 日本vs欧美在线观看视频 | 网址你懂的国产日韩在线| 久久鲁丝午夜福利片| 在线看a的网站| 精华霜和精华液先用哪个| 丰满人妻一区二区三区视频av| 男女国产视频网站| 最后的刺客免费高清国语| 少妇熟女欧美另类| 欧美极品一区二区三区四区| 日韩一区二区视频免费看| 五月开心婷婷网| 亚洲精品aⅴ在线观看| 日韩精品有码人妻一区| 日韩三级伦理在线观看| 精品少妇久久久久久888优播| 在现免费观看毛片| 国产一级毛片在线| 涩涩av久久男人的天堂| 亚洲欧美成人精品一区二区| 在线天堂最新版资源| 2022亚洲国产成人精品| 精品国产三级普通话版| 一级毛片电影观看| 久久 成人 亚洲| 人妻少妇偷人精品九色| 中文字幕免费在线视频6| 联通29元200g的流量卡| 国产成人午夜福利电影在线观看| 夫妻性生交免费视频一级片| 在线看a的网站| 免费观看a级毛片全部| 午夜激情久久久久久久| 国精品久久久久久国模美| 女性被躁到高潮视频| 免费看日本二区| 联通29元200g的流量卡| 麻豆成人av视频| 插逼视频在线观看| 免费人妻精品一区二区三区视频| 亚洲欧美成人综合另类久久久| 成人漫画全彩无遮挡| 亚洲四区av| 欧美成人精品欧美一级黄| 大又大粗又爽又黄少妇毛片口| 看免费成人av毛片| 中文字幕av成人在线电影| 国产精品秋霞免费鲁丝片| 免费久久久久久久精品成人欧美视频 | 一级a做视频免费观看| 观看美女的网站| 亚洲综合精品二区| 欧美丝袜亚洲另类| av在线app专区| 午夜免费男女啪啪视频观看| 人人妻人人澡人人爽人人夜夜| 亚洲av国产av综合av卡| 日韩伦理黄色片| 亚洲国产精品一区三区| 免费观看的影片在线观看| 国产一级毛片在线| 男人添女人高潮全过程视频| 伊人久久精品亚洲午夜| 亚洲精品自拍成人| 久久久久久久久大av| 少妇裸体淫交视频免费看高清| 成年人午夜在线观看视频| 久久精品国产鲁丝片午夜精品| 十八禁网站网址无遮挡 | 99re6热这里在线精品视频| 日本-黄色视频高清免费观看| 精品国产乱码久久久久久小说| 少妇的逼好多水| 国产精品久久久久成人av| av线在线观看网站| 欧美高清成人免费视频www| 成年美女黄网站色视频大全免费 | 亚洲精品中文字幕在线视频 | 在线 av 中文字幕| 亚洲av不卡在线观看| 中国美白少妇内射xxxbb| 交换朋友夫妻互换小说| 少妇人妻 视频| 天堂8中文在线网| 少妇人妻精品综合一区二区| 少妇精品久久久久久久| 色综合色国产| 婷婷色麻豆天堂久久| 亚洲一级一片aⅴ在线观看| 天堂8中文在线网| 国产精品嫩草影院av在线观看| av线在线观看网站| 97超碰精品成人国产| 永久网站在线| 最黄视频免费看| 国产视频首页在线观看| 中文字幕精品免费在线观看视频 | 国产精品国产三级国产专区5o| 成人美女网站在线观看视频| 亚洲中文av在线| 国产永久视频网站| 精品久久久久久电影网| 一本色道久久久久久精品综合| 亚洲aⅴ乱码一区二区在线播放| 久久99热6这里只有精品| 成人综合一区亚洲| 精品国产一区二区三区久久久樱花 | 大陆偷拍与自拍| 免费观看a级毛片全部| 伦理电影免费视频| 久久久色成人| 涩涩av久久男人的天堂| 日本黄色日本黄色录像| 亚洲精品亚洲一区二区| 啦啦啦在线观看免费高清www| 最后的刺客免费高清国语| 中文字幕免费在线视频6| 久久av网站| 边亲边吃奶的免费视频| 中国三级夫妇交换| 成人无遮挡网站| 欧美丝袜亚洲另类| 欧美变态另类bdsm刘玥| 欧美精品人与动牲交sv欧美| 欧美国产精品一级二级三级 | 精品久久久久久久末码| 香蕉精品网在线| 大香蕉久久网| av国产免费在线观看| 久久韩国三级中文字幕| 丝瓜视频免费看黄片| 色综合色国产| 亚洲精品日韩av片在线观看| 在线观看美女被高潮喷水网站| 婷婷色av中文字幕| 麻豆成人av视频| 少妇人妻一区二区三区视频| 久久99热这里只频精品6学生| 久久久久久人妻| 国产高清不卡午夜福利| 亚洲av男天堂| 在线 av 中文字幕| 国产av一区二区精品久久 | 久久久久人妻精品一区果冻| 91精品国产九色| 在线天堂最新版资源| 国产精品三级大全| 午夜激情福利司机影院| 久久这里有精品视频免费| 亚洲精品自拍成人| 青春草亚洲视频在线观看| 久久精品久久久久久久性| 欧美人与善性xxx| 久久久久精品久久久久真实原创| 高清不卡的av网站| 国产欧美日韩一区二区三区在线 | 在线观看美女被高潮喷水网站| 极品教师在线视频| 不卡视频在线观看欧美| 国产精品无大码| 一级毛片久久久久久久久女| 女的被弄到高潮叫床怎么办| 97超碰精品成人国产| 精品人妻视频免费看| 色视频在线一区二区三区| 超碰97精品在线观看| 成人高潮视频无遮挡免费网站| 亚洲av成人精品一二三区| 亚洲精品久久午夜乱码| 欧美成人精品欧美一级黄| www.av在线官网国产| 五月伊人婷婷丁香| 啦啦啦中文免费视频观看日本| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久久电影| 九色成人免费人妻av| 日韩 亚洲 欧美在线| 赤兔流量卡办理| 2022亚洲国产成人精品| 免费大片18禁| 婷婷色av中文字幕| 亚洲久久久国产精品| 欧美bdsm另类| 美女内射精品一级片tv| 又爽又黄a免费视频| 精品国产露脸久久av麻豆| 99久久中文字幕三级久久日本| 午夜福利影视在线免费观看| 老师上课跳d突然被开到最大视频| 99视频精品全部免费 在线| 亚洲精品乱久久久久久| 女性被躁到高潮视频| 日韩欧美一区视频在线观看 | 欧美人与善性xxx| 日本色播在线视频| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式| 国产有黄有色有爽视频| 成人国产麻豆网| 插阴视频在线观看视频| 又黄又爽又刺激的免费视频.| 午夜福利网站1000一区二区三区| 丰满乱子伦码专区| 色视频www国产| 成人高潮视频无遮挡免费网站| 日韩不卡一区二区三区视频在线| 少妇猛男粗大的猛烈进出视频| 高清毛片免费看| 国产熟女欧美一区二区| 久久久久久久久大av| 欧美日韩国产mv在线观看视频 | 久久99热这里只有精品18| 简卡轻食公司| 国产精品不卡视频一区二区| 国产精品成人在线| 高清不卡的av网站| av黄色大香蕉| 三级国产精品欧美在线观看| 啦啦啦在线观看免费高清www| 久久精品夜色国产| 91狼人影院| 中文欧美无线码| 六月丁香七月| 国产深夜福利视频在线观看| 欧美最新免费一区二区三区| 亚洲美女视频黄频| 我要看黄色一级片免费的| 国产成人精品一,二区| 国产成人精品久久久久久| 国产欧美日韩一区二区三区在线 | 免费观看无遮挡的男女| 精品亚洲成a人片在线观看 | 插阴视频在线观看视频| 久热这里只有精品99| 亚洲美女黄色视频免费看| 好男人视频免费观看在线| 免费人成在线观看视频色| 日本与韩国留学比较| 97热精品久久久久久| 一级毛片aaaaaa免费看小| 在线观看免费日韩欧美大片 | 人妻一区二区av| 视频区图区小说| 日韩亚洲欧美综合| 国产成人午夜福利电影在线观看| 成人国产麻豆网| 高清欧美精品videossex| av专区在线播放| 一区二区三区精品91| 毛片女人毛片| 另类亚洲欧美激情| 久久精品国产自在天天线| 五月天丁香电影| 国产成人freesex在线| 18禁裸乳无遮挡动漫免费视频| 韩国高清视频一区二区三区| 女性生殖器流出的白浆| 一边亲一边摸免费视频| 美女视频免费永久观看网站| 欧美日韩综合久久久久久| 午夜激情久久久久久久| 国产精品一区www在线观看| 色综合色国产| 有码 亚洲区| 内地一区二区视频在线| 成人毛片60女人毛片免费| 一级二级三级毛片免费看| 国产一区有黄有色的免费视频| 91午夜精品亚洲一区二区三区| av国产精品久久久久影院| 老女人水多毛片| 在线观看免费日韩欧美大片 | 91狼人影院| 七月丁香在线播放| 人妻系列 视频| 国产精品久久久久久久电影| 成人漫画全彩无遮挡| 国产久久久一区二区三区| 午夜福利高清视频| av卡一久久| 精品视频人人做人人爽| 最黄视频免费看| 亚洲精品日本国产第一区| 国产精品一区二区三区四区免费观看| 最近最新中文字幕免费大全7| av女优亚洲男人天堂| 国产一区二区三区综合在线观看 | 久久久久精品久久久久真实原创| 性色avwww在线观看| 国产成人午夜福利电影在线观看| 国产精品国产三级专区第一集| 亚洲婷婷狠狠爱综合网| 中文字幕久久专区| 免费观看a级毛片全部| 在现免费观看毛片| 日本av手机在线免费观看| 在线观看免费高清a一片| 国产日韩欧美亚洲二区| 精品人妻偷拍中文字幕| 美女福利国产在线 | 黑丝袜美女国产一区| 国产视频首页在线观看| 精品人妻一区二区三区麻豆| 国产黄色免费在线视频| 大片免费播放器 马上看| 欧美高清成人免费视频www| 亚洲av.av天堂| 亚洲av综合色区一区| 人体艺术视频欧美日本| 久久久精品94久久精品| 涩涩av久久男人的天堂| 亚洲精品,欧美精品| 毛片女人毛片| 内射极品少妇av片p| 久久精品国产自在天天线| 亚洲精品久久午夜乱码| 人妻系列 视频| 国产色爽女视频免费观看| 久久久久久久久久成人| 熟女电影av网| av免费观看日本| 一级片'在线观看视频| 成人特级av手机在线观看| 免费看av在线观看网站| 亚洲欧美精品自产自拍| 中国国产av一级| 精品久久久噜噜| 最新中文字幕久久久久| 国产成人aa在线观看| 欧美区成人在线视频| 日本爱情动作片www.在线观看| 午夜视频国产福利| 国产成人免费观看mmmm| 免费看av在线观看网站| 国内少妇人妻偷人精品xxx网站| 久久午夜福利片| 亚洲在久久综合| av一本久久久久| 国产成人a区在线观看| 亚洲欧洲国产日韩| 美女高潮的动态| 午夜日本视频在线| av播播在线观看一区| 丰满迷人的少妇在线观看| 亚洲第一av免费看| 看非洲黑人一级黄片| 人妻少妇偷人精品九色| 你懂的网址亚洲精品在线观看| 黄色怎么调成土黄色| 成人午夜精彩视频在线观看| 99久久中文字幕三级久久日本| 欧美 日韩 精品 国产| 久久精品人妻少妇| 中文字幕精品免费在线观看视频 | 国产熟女欧美一区二区| 一二三四中文在线观看免费高清| 日日摸夜夜添夜夜添av毛片| 啦啦啦在线观看免费高清www| 亚洲国产精品国产精品| 99热国产这里只有精品6| 妹子高潮喷水视频| 这个男人来自地球电影免费观看 | 久久99热6这里只有精品| 男女免费视频国产| 成人无遮挡网站| 三级经典国产精品| 亚洲,一卡二卡三卡| 亚洲国产精品999| 国产亚洲最大av| 欧美一级a爱片免费观看看| 少妇人妻 视频| 蜜桃亚洲精品一区二区三区| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品电影小说 | 亚洲av中文av极速乱| 国产 精品1| 纵有疾风起免费观看全集完整版| 毛片女人毛片| 精品国产一区二区三区久久久樱花 | 国产精品一区二区在线不卡| 成人午夜精彩视频在线观看| 欧美少妇被猛烈插入视频| 欧美97在线视频| 中国三级夫妇交换| 又黄又爽又刺激的免费视频.| 高清日韩中文字幕在线| 免费黄网站久久成人精品| 中文欧美无线码| 国产无遮挡羞羞视频在线观看| 欧美精品一区二区免费开放| 青春草国产在线视频| freevideosex欧美| 日韩 亚洲 欧美在线| 欧美精品人与动牲交sv欧美| 亚洲av免费高清在线观看| 夫妻性生交免费视频一级片| 亚洲婷婷狠狠爱综合网| av福利片在线观看| 在线观看免费高清a一片| 黄色配什么色好看| 国产av国产精品国产| www.av在线官网国产| 国产成人免费无遮挡视频| 成人高潮视频无遮挡免费网站| 亚洲色图综合在线观看| 五月开心婷婷网| 国产成人a∨麻豆精品| 亚洲经典国产精华液单| 国产免费视频播放在线视频| 五月开心婷婷网| 王馨瑶露胸无遮挡在线观看| 全区人妻精品视频| av不卡在线播放| 寂寞人妻少妇视频99o| 中文乱码字字幕精品一区二区三区| 日本-黄色视频高清免费观看| 一本一本综合久久| 大话2 男鬼变身卡| 国产人妻一区二区三区在| 国产爱豆传媒在线观看| 亚洲国产色片| 国内少妇人妻偷人精品xxx网站| 午夜激情福利司机影院| 亚洲精品久久久久久婷婷小说| 免费在线观看成人毛片| 中文欧美无线码| 韩国高清视频一区二区三区| 免费观看a级毛片全部| 亚洲精品第二区| 久久99精品国语久久久| 久久99热6这里只有精品| 自拍偷自拍亚洲精品老妇| 欧美丝袜亚洲另类| 精华霜和精华液先用哪个| 国产深夜福利视频在线观看| 色吧在线观看| 联通29元200g的流量卡| 精品久久国产蜜桃| 伦精品一区二区三区| 高清不卡的av网站| 国产男女超爽视频在线观看| 一本一本综合久久| 亚洲av欧美aⅴ国产| 久久精品国产亚洲网站| 黄片无遮挡物在线观看| 亚洲精品国产av蜜桃| 欧美日韩在线观看h| 蜜桃久久精品国产亚洲av| 一级毛片 在线播放| 国产色婷婷99| 成人特级av手机在线观看| 色5月婷婷丁香| 亚洲色图综合在线观看| av天堂中文字幕网| 中文乱码字字幕精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 99热全是精品| 看非洲黑人一级黄片| 高清在线视频一区二区三区| 大码成人一级视频| 亚洲丝袜综合中文字幕|