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

    基于光流算法的粒子圖像測(cè)速技術(shù)研究

    2015-06-22 14:46:36王宏偉
    實(shí)驗(yàn)流體力學(xué) 2015年3期
    關(guān)鍵詞:光流流場(chǎng)灰度

    王宏偉, 黃 湛

    (中國(guó)航天空氣動(dòng)力技術(shù)研究院, 北京 100074)

    基于光流算法的粒子圖像測(cè)速技術(shù)研究

    王宏偉, 黃 湛*

    (中國(guó)航天空氣動(dòng)力技術(shù)研究院, 北京 100074)

    光流測(cè)量技術(shù)作為一種新的空氣動(dòng)力學(xué)實(shí)驗(yàn)技術(shù),以其像素級(jí)分辨率的矢量場(chǎng)測(cè)量?jī)?yōu)勢(shì)獲得廣泛的應(yīng)用。光流測(cè)量技術(shù)使用光流約束方程,配合平滑限定條件,可以進(jìn)行速度場(chǎng)測(cè)量,獲得高分辨率的全局矢量場(chǎng)。首先通過研究積分最小化光流測(cè)速理論和算法,采用C++編寫光流速度測(cè)量程序,然后通過3種典型人工位移圖像對(duì)光流計(jì)算程序進(jìn)行驗(yàn)證,并將結(jié)果和標(biāo)準(zhǔn)位移分布進(jìn)行比對(duì)分析,以指導(dǎo)如何在實(shí)際應(yīng)用中獲得高精度光流速度場(chǎng),最后進(jìn)行小型風(fēng)洞后向臺(tái)階實(shí)驗(yàn),利用高速相機(jī)拍攝示蹤粒子圖像,使用光流計(jì)算程序獲得速度矢量場(chǎng),同采用互相關(guān)算法的粒子圖像測(cè)速計(jì)算結(jié)果進(jìn)行比較,體現(xiàn)出光流計(jì)算方法像素級(jí)分辨率的矢量場(chǎng)測(cè)量?jī)?yōu)勢(shì)。

    光流;速度測(cè)量;粒子圖像;運(yùn)動(dòng)估計(jì);積分最小化

    0 引 言

    在空氣動(dòng)力學(xué)發(fā)展過程中,應(yīng)用實(shí)驗(yàn)手段獲得流體的空間運(yùn)動(dòng)結(jié)構(gòu),獲取全場(chǎng)流動(dòng)信息,是了解空氣動(dòng)力學(xué)諸多流動(dòng)現(xiàn)象及機(jī)理的關(guān)鍵,同時(shí)這一需求也不斷刺激空氣動(dòng)力學(xué)實(shí)驗(yàn)測(cè)試技術(shù)水平的提升。傳統(tǒng)的流動(dòng)測(cè)量方法大多采用單點(diǎn)測(cè)量或介入式測(cè)量手段,不能滿足復(fù)雜流動(dòng)情況以及全場(chǎng)速度測(cè)量的需求。

    近年來,粒子圖像測(cè)速技術(shù)利用拍攝跟隨流場(chǎng)運(yùn)動(dòng)的粒子圖像并通過相關(guān)計(jì)算獲得速度場(chǎng),并以其非介入式、全場(chǎng)流動(dòng)速度測(cè)量、測(cè)量精度高及測(cè)速范圍廣等優(yōu)點(diǎn)在全局速度測(cè)量方面具有應(yīng)用優(yōu)勢(shì)。不過,粒子圖像測(cè)速得到的速度矢量為其判讀區(qū)域內(nèi)粒群的平均速度,忽略了此判讀區(qū)域內(nèi)速度場(chǎng)的變化,這使得粒子圖像測(cè)速方法的空間分辨率受到一定的限制,尤其對(duì)于速度梯度大、速度變化劇烈的區(qū)域,其速度場(chǎng)測(cè)量精度是受到制約的。

    相機(jī)所記錄的物體的圖像,本質(zhì)上是物體所反射的光的分布。當(dāng)相機(jī)連續(xù)拍攝運(yùn)動(dòng)物體時(shí),物體所反射的光就會(huì)在相機(jī)底片上形成一系列連續(xù)變化的圖像,由這些圖像所組成的連續(xù)變化的信息不斷“流過”相機(jī)底片,如同光在流動(dòng),這是光流定義的物理基礎(chǔ)。物體在光的照射下,其表面灰度會(huì)呈現(xiàn)出一定的分布,而所謂光流,就是指圖像中灰度模式的運(yùn)動(dòng)速度。粒子圖像測(cè)速技術(shù)中拍攝到的粒子圖像顯示的是粒子濃度的分布,具有一定濃度的運(yùn)動(dòng)粒子群通過其散射光在相機(jī)底片上成像,所以粒子圖像測(cè)速技術(shù)中采用相關(guān)方法得到的速度場(chǎng)實(shí)際上是光流速度場(chǎng)(Optical Flow),光流速度場(chǎng)能否代表真實(shí)速度場(chǎng)要取決于示蹤物質(zhì)能否充分地跟隨流體運(yùn)動(dòng)。

    Horn和Schunck早在20世紀(jì)80年代就開創(chuàng)性地進(jìn)行了光流研究和計(jì)算[3],Horn等認(rèn)為同一運(yùn)動(dòng)物體引起的光流場(chǎng)應(yīng)該是連續(xù)的、平滑的,即同一物體上相鄰點(diǎn)的速度是相似的,那么其投影到圖像上的光流變化也應(yīng)該是平滑的,并提出了一種利用加在光流場(chǎng)上的附加速度平滑約束,即全局平滑約束將光流場(chǎng)的計(jì)算問題轉(zhuǎn)化為一個(gè)變分問題,給出像素級(jí)分辨率的運(yùn)動(dòng)估計(jì)。隨后,光流技術(shù)在發(fā)展過程中吸引了各個(gè)領(lǐng)域的科研人員進(jìn)行研究[4-6],也產(chǎn)生了很多其它解算方法,如圖像相關(guān)測(cè)速法、基于能量的方法[7]、基于相位的方法[8]和神經(jīng)動(dòng)力學(xué)方法等,針對(duì)Horn等提出的微分法在計(jì)算較大像素位移時(shí)遇到的困難,Tianshu Liu等在其工作中提到[9],Corpetti等人提出在時(shí)間步上采用積分形式的光流方程代替微分形式的變分方程[10],Ruhnau等使用高斯濾鏡生成二元圖像金字塔的多分辨率策略[11]。這些算法有各自的優(yōu)勢(shì),不過由于具有一定的適用范圍和算法復(fù)雜度使其相對(duì)于微分法無法獲得廣泛應(yīng)用。

    本文引入光流計(jì)算方法,編寫光流速度測(cè)量程序,通過人工位移圖像對(duì)該程序進(jìn)行驗(yàn)證和分析,以合理的測(cè)速范圍和最優(yōu)的加權(quán)值指導(dǎo)如何在實(shí)際應(yīng)用中獲得高精度光流速度場(chǎng);進(jìn)行小型風(fēng)洞實(shí)驗(yàn),以粒子圖像為速度測(cè)量載體,使用光流計(jì)算方法獲得像素級(jí)分辨率的速度矢量場(chǎng)。旨在通過這些研究,發(fā)展以光流技術(shù)為基礎(chǔ)的全場(chǎng)流動(dòng)測(cè)量技術(shù),期望在提高空間分辨率和避免速度梯度影響、了解復(fù)雜流動(dòng)的空間結(jié)構(gòu)、獲取更多的動(dòng)力學(xué)信息方面獲得進(jìn)步,以此豐富實(shí)驗(yàn)流體力學(xué)測(cè)量技術(shù)。

    1 基于粒子圖像的光流速度場(chǎng)求解

    Horn的光流場(chǎng)解算方法是微分法的基礎(chǔ),它以灰度守恒方程為基礎(chǔ),附加速度平滑約束條件,利用連續(xù)變化的圖像上各點(diǎn)灰度的時(shí)間、空間梯度來計(jì)算每一點(diǎn)的速度矢量。

    光流的灰度約束條件認(rèn)為,圖像上各點(diǎn)灰度隨時(shí)間和空間變化滿足如下關(guān)系:

    式中:I(x,y,t)是t時(shí)刻(x,y)位置的圖像灰度,I(x+Δx,y+Δy,t+Δt)為原(x,y)位置的物體運(yùn)動(dòng)到(x+Δx,y+Δy,t+Δt)位置的灰度,令u=Δx/Δt,v=Δy/Δt,將上式進(jìn)行泰勒展開并忽略高階項(xiàng)可以獲得光流約束方程:

    式中:Ix,Iy,It是灰度的二維空間和時(shí)間梯度。顯然該方程是求解u,v2個(gè)未知數(shù)的不定方程,需要引入輔助約束方程進(jìn)行求解。

    速度梯度平方和判定量為

    實(shí)際拍攝的圖像灰度將會(huì)受到硬件條件和噪聲的影響而偏離基本的光流約束條件,所以不能肯定地認(rèn)為ξI恒為0,因此有必要找到合適的加權(quán)因子β,使得2個(gè)判定誤差在二維空間上的積分為最小,即對(duì)

    ξ2=?

    進(jìn)行最小化。通過變分法建立每個(gè)因變量的歐拉特征方程:

    2 光流速度場(chǎng)求解算法驗(yàn)證

    基于微分法編寫了光流速度場(chǎng)計(jì)算程序,擬用于風(fēng)洞試驗(yàn)速度場(chǎng)測(cè)量。為了驗(yàn)證光流算法的精度,本文選取灰度圖像,通過人工位移獲得標(biāo)準(zhǔn)速度場(chǎng),然后采用所編寫的光流計(jì)算程序?qū)υ搱D像對(duì)進(jìn)行處理,將光流計(jì)算結(jié)果與標(biāo)準(zhǔn)速度場(chǎng)進(jìn)行比較。實(shí)際計(jì)算中,將采用1個(gè)單位的時(shí)間梯度,所得的速度值將直接反映圖像像素位移值,同時(shí),也可考察2張圖像的像素位移對(duì)光流計(jì)算精度的影響。

    首先考察第1組圖像(見圖1),圖像為600pixel×600pixel大小的風(fēng)暴云圖,通過處理將圖像向右平移一個(gè)像素(Δx=1pixel),通過光流速度場(chǎng)計(jì)算程序得到的速度矢量分布(Skip 16pixel×16pixel,下同)和流線如圖2所示,可以看出,在1像素位移條件下,光流計(jì)算方法所得到的速度場(chǎng)分布十分均勻。

    圖1 灰度原圖與平移一個(gè)像素后圖像

    圖2 采用光流計(jì)算方法得到的速度矢量圖和流線圖

    圖3 橫向速度分布和縱向速度分布分析

    在圖1中所示的縱向考察線和橫向考察線(下同)上對(duì)光流速度場(chǎng)和標(biāo)準(zhǔn)速度進(jìn)行分析如圖3所示,可以看出, 在1像素位移條件下光流算法獲得的速度值與標(biāo)準(zhǔn)速度幾乎沒有偏差。

    第2組圖像(見圖4)采取了Δx=1 pixel和Δy=1 pixel的斜向移動(dòng),處理得到的速度矢量分布和流線分布如圖5所示,速度矢量分布很均勻,流場(chǎng)平滑性很好,不過已開始出現(xiàn)小范圍的速度偏差。

    通過縱向考察線和橫向考察線將所計(jì)算的速度場(chǎng)與標(biāo)準(zhǔn)速度進(jìn)行比較發(fā)現(xiàn)(見圖6), 在Δx=1 pixel和Δy=1 pixel的斜向運(yùn)動(dòng)條件下,所計(jì)算的速度值能夠保證較高的準(zhǔn)度分布,不過已開始出現(xiàn)一定量的抖動(dòng)偏差。這種偏差很有可能與2個(gè)方向速度偏差的合成有關(guān);同時(shí)也暗示隨著連續(xù)圖像間運(yùn)動(dòng)位移的增大,光流算法的計(jì)算誤差也會(huì)隨之增加。

    圖4 灰度原圖與x,y各平移一個(gè)像素后圖像

    圖5 采用光流計(jì)算方法得到的速度矢量圖和流線圖

    圖6 橫向速度分布和縱向速度分布與標(biāo)準(zhǔn)速度比較

    第3組圖像(見圖7)為中心旋轉(zhuǎn)1°前后的2張圖像,通過計(jì)算獲得的速度矢量和流線分布圖如圖8所示,可以看出,圖像中心的速度矢量分布較為均勻和平滑,邊緣處的速度偏差較大,平滑性較差。

    通過在前述2個(gè)方向上將光流計(jì)算的速度矢量和標(biāo)準(zhǔn)速度進(jìn)行比較(見圖9),可以明顯發(fā)現(xiàn):圖像位移量不超過2個(gè)像素部分的計(jì)算結(jié)果與標(biāo)準(zhǔn)速度結(jié)果十分吻合,偏差很??;圖像位移量超過2個(gè)像素部分的計(jì)算結(jié)果則出現(xiàn)了明顯的偏差。這很好地說明了所得到的速度場(chǎng)分布中為何中心區(qū)域速度分布均勻,平滑性好,而邊緣區(qū)域速度偏差較大,平滑性較差,同時(shí)也在很大程度上對(duì)如何在實(shí)際測(cè)量應(yīng)用中保證光流計(jì)算獲得高精度的結(jié)果進(jìn)行了指導(dǎo)。

    圖7 灰度原圖與旋轉(zhuǎn)1度后的圖像

    圖8 采用光流計(jì)算方法得到的速度矢量圖和流線圖

    圖9 橫向速度分布和縱向速度分布與標(biāo)準(zhǔn)速度比較

    3 加權(quán)因子最優(yōu)值估計(jì)

    β2作為加權(quán)因子表征了實(shí)驗(yàn)和數(shù)值離散噪音的相對(duì)大小。在實(shí)際計(jì)算中β2通過調(diào)整速度梯度平方和比重來影響計(jì)算結(jié)果的平滑性,即通過選取不同的加權(quán)值獲得不同平滑程度的流場(chǎng)速度分布。下面將以旋轉(zhuǎn)灰度圖像為計(jì)算基礎(chǔ),在不同β2值下對(duì)光流計(jì)算結(jié)果中縱向考察線速度分布與標(biāo)準(zhǔn)速度分布進(jìn)行比較(見圖10~14)。

    由圖中數(shù)據(jù)可以看出,隨β2增大,縱向考察線上速度場(chǎng)的分布能夠更好地吻合標(biāo)準(zhǔn)速度分布。在實(shí)際應(yīng)用粒子圖像進(jìn)行光流速度場(chǎng)計(jì)算時(shí),β2相對(duì)較大,說明噪聲較大,這與粒子圖像不同于激光誘導(dǎo)熒光這類圖像有關(guān),粒子圖像的灰度分布存在一定的間斷性,需要大的加權(quán)值來限定誤差。

    圖10 β2=1時(shí)縱向考察線速度分布與標(biāo)準(zhǔn)速度分布比較

    圖11 β2=4時(shí)縱向考察線速度分布與標(biāo)準(zhǔn)速度分布比較

    圖12 β2=7時(shí)縱向考察線速度分布與標(biāo)準(zhǔn)速度分布比較

    圖13 β2=10時(shí)縱向考察線速度分布與標(biāo)準(zhǔn)速度分布比較

    圖14 β2=13時(shí)縱向考察線速度分布與標(biāo)準(zhǔn)速度分布比較

    為了進(jìn)一步分析β2的選取對(duì)速度場(chǎng)特性影響,根據(jù)不同β2值解算出的光流速度場(chǎng)與標(biāo)準(zhǔn)速度分布之間的相關(guān)系數(shù)來表現(xiàn)β2對(duì)解吻合的影響(見圖15)。根據(jù)圖中數(shù)據(jù)結(jié)果可以看出,解的準(zhǔn)確性在β2增長(zhǎng)初期會(huì)迅速增加,在10~18之間達(dá)到平穩(wěn)階段并取得極值,超過20以后,反而會(huì)下降,因此選取適當(dāng)?shù)摩?值將對(duì)光流算法獲得的速度場(chǎng)的準(zhǔn)確性有很大幫助。不過這里β2是與灰度的時(shí)間和空間梯度間隔相關(guān)的常數(shù),實(shí)際應(yīng)用中應(yīng)結(jié)合具體參數(shù)進(jìn)行選取計(jì)算。

    圖15 相關(guān)系數(shù)隨β2變化曲線

    4 光流速度場(chǎng)測(cè)量試驗(yàn)與結(jié)果

    本文針對(duì)光流測(cè)速算法對(duì)流場(chǎng)的要求,進(jìn)行了后向臺(tái)階內(nèi)流場(chǎng)渦結(jié)構(gòu)觀測(cè)實(shí)驗(yàn)。

    4.1 光流速度場(chǎng)測(cè)量試驗(yàn)

    實(shí)驗(yàn)的難點(diǎn)在于如何保證實(shí)驗(yàn)條件滿足光流流場(chǎng)計(jì)算的約束方程。根據(jù)前述分析,應(yīng)用如上微分法對(duì)光流流場(chǎng)進(jìn)行解算時(shí),要求連續(xù)拍攝的每張圖像具有相同的光照條件,相鄰圖像對(duì)應(yīng)像素點(diǎn)的位移不至太大,最好保證在2個(gè)像素以內(nèi),從而進(jìn)行有效的光流速度場(chǎng)計(jì)算。

    實(shí)驗(yàn)在實(shí)驗(yàn)室自行搭建的小型低速風(fēng)洞系統(tǒng)內(nèi)進(jìn)行,整個(gè)系統(tǒng)由風(fēng)洞主體、流動(dòng)控制系統(tǒng)、示蹤粒子發(fā)生和播發(fā)系統(tǒng)、激光照明及附屬光路系統(tǒng)、圖像采集系統(tǒng)和計(jì)算機(jī)處理系統(tǒng)組成。實(shí)驗(yàn)布局如圖16所示。

    圖16 光流測(cè)速實(shí)驗(yàn)布局

    風(fēng)洞主體通過其前段可控轉(zhuǎn)速的渦輪風(fēng)扇電機(jī)控制流動(dòng)速度,通過調(diào)節(jié)變頻器可以將流經(jīng)中部實(shí)驗(yàn)段的流動(dòng)速度控制在1~20m/s的范圍內(nèi),風(fēng)洞主體通過其尾部的蜂窩整流器吸入混有示蹤粒子的氣流并進(jìn)行整流,保證來流的穩(wěn)定條件,并與3個(gè)示蹤粒子播撒管共同工作使示蹤粒子與吸入氣流混合均勻。為保證恒定的光照條件,實(shí)驗(yàn)采用氬離子激光器輸出的連續(xù)激光通過擴(kuò)束鏡和片光源為流場(chǎng)提供連續(xù)的片光照明;為保證相鄰圖像間位移不至太大,實(shí)驗(yàn)采用高速相機(jī)進(jìn)行圖像采集,觀測(cè)流場(chǎng)前先對(duì)標(biāo)尺進(jìn)行拍攝用以標(biāo)定空間尺度,采集結(jié)束后傳輸?shù)接?jì)算機(jī)上存儲(chǔ)并處理,現(xiàn)場(chǎng)設(shè)備及布置照片如圖17~20所示。

    圖17 實(shí)驗(yàn)光路及實(shí)驗(yàn)段

    圖18 粒子播撒架

    圖19 粒子發(fā)生器

    4.2 粒子圖像光流場(chǎng)分析與圖像均一化

    高速相機(jī)拍攝的后向臺(tái)階內(nèi)流場(chǎng)粒子圖像如圖21所示,大量連續(xù)拍攝的粒子圖像記錄了連續(xù)變化的光流流場(chǎng)信息,虛線框內(nèi)為計(jì)算區(qū)域。為了對(duì)實(shí)際拍攝示蹤粒子圖像的光流場(chǎng)特性進(jìn)行驗(yàn)證,將實(shí)際拍攝的光流場(chǎng)圖像進(jìn)行旋轉(zhuǎn)獲得標(biāo)準(zhǔn)位移場(chǎng),如圖22所示;采用與前述光流解算方法驗(yàn)證相同的比較方式對(duì)計(jì)算獲得的速度場(chǎng)進(jìn)行考察,如圖23和24所示;并分析了速度分布與標(biāo)準(zhǔn)速度分布間的相對(duì)誤差和絕對(duì)誤差,如圖25所示。

    圖20 實(shí)驗(yàn)現(xiàn)場(chǎng)

    圖21 后向臺(tái)階流場(chǎng)粒子圖像及流場(chǎng)計(jì)算區(qū)域

    圖22 使用示蹤粒子圖像旋轉(zhuǎn)1°前后圖像及驗(yàn)證計(jì)算區(qū)域

    圖23 用光流算法對(duì)旋轉(zhuǎn)前后粒子圖像進(jìn)行計(jì)算的速度矢量與流線圖

    分析可以看出,在示蹤粒子圖像上,當(dāng)像素位移小于一定值時(shí),計(jì)算獲得的速度分布能夠很好的與標(biāo)準(zhǔn)值吻合;橫向考察線和縱向考察線上速度分布絕對(duì)誤差基本處于5%線以下,不過由于靠近0位移位置時(shí),速度絕對(duì)值較小,故而相對(duì)誤差較大;在±1像素位置絕對(duì)誤差和相對(duì)誤差最小,都達(dá)到了1%左右,該位置±0.5像素鄰域是速度相對(duì)精度和絕對(duì)精度都比較高的區(qū)域。

    圖24 橫向速度分布和縱向速度分布與標(biāo)準(zhǔn)速度比較

    圖25 考察線上速度分布的絕對(duì)誤差和相對(duì)誤差

    由于激光經(jīng)過柱透鏡制造的片光存在一定的擴(kuò)張角度,在激光片光傳播過程中,片光的寬度始終在變大,所以激光片光在展開區(qū)域的能量分布隨時(shí)空變化也會(huì)有一定的差異。此外由于實(shí)驗(yàn)段的有機(jī)玻璃制造工藝以及存在一定的磨損,其中存在許多細(xì)微的氣泡或表面有輕微劃痕,當(dāng)激光入射進(jìn)有機(jī)玻璃遇到這些氣泡或劃痕時(shí),激光就要被削減一部分能量,導(dǎo)致射入流場(chǎng)里的激光片光含有大量的條紋,即激光片光的能量分布是不連續(xù)的,自然所得的灰度場(chǎng)里也含有大量的條紋。

    如果還要考慮噪聲的話,可以在實(shí)驗(yàn)前或?qū)嶒?yàn)后獲得一個(gè)噪聲灰度場(chǎng)。獲得噪聲灰度場(chǎng)的一種方法是,在不加示蹤物質(zhì)的情況下,將攝像機(jī)對(duì)準(zhǔn)拍攝區(qū)域,連續(xù)拍攝一段時(shí)間,將得到的圖像進(jìn)行平均,就可以得到噪聲灰度場(chǎng)INoise。這時(shí)相對(duì)灰度場(chǎng)的表達(dá)如式(9)所示。

    4.3 試驗(yàn)結(jié)果與分析

    經(jīng)過光流算法計(jì)算得到的真實(shí)流場(chǎng)速度分布如圖26所示,結(jié)果顯示最大位移為2.2像素,表明區(qū)域內(nèi)速度計(jì)算結(jié)果基本位于可信區(qū)間內(nèi)。

    由于所觀測(cè)流場(chǎng)的絕對(duì)速度較低(<1m/s),后向臺(tái)階內(nèi)流場(chǎng)結(jié)構(gòu)不具有典型性,所以將計(jì)算區(qū)域內(nèi)的光流處理結(jié)果同采用互相關(guān)算法的PIV結(jié)果(見圖27)進(jìn)行比較。圖28為PIV顯示的流場(chǎng)結(jié)構(gòu),同圖29的光流算法流場(chǎng)結(jié)構(gòu)比較可以發(fā)現(xiàn),在同樣像素分辨率下,光流算法的結(jié)果在速度大小、方向和等速度線分布上基本與PIV計(jì)算結(jié)果基本一致,大尺度渦結(jié)構(gòu)和位置也與其基本無差異,但光流計(jì)算結(jié)果中可以觀測(cè)到更小尺度的渦結(jié)構(gòu),等速度線更加平滑。

    圖26 光流計(jì)算結(jié)果(Skip 8pixel×8pixel)

    圖27 PIV計(jì)算結(jié)果(8pixel×8pixel判讀區(qū))

    比較兩者流線圖(圖28和29)可以看出,在同樣像素分辨率下,PIV算法的計(jì)算結(jié)果中較大尺度的渦結(jié)構(gòu)可以顯示出來,但是對(duì)于比較小尺度的渦結(jié)構(gòu)顯示卻不明顯,而光流算法的計(jì)算結(jié)果中小尺度渦結(jié)構(gòu)被明顯地顯示出來。

    為進(jìn)一步分析,將兩者小尺度渦結(jié)構(gòu)區(qū)域圖像放大(圖30和31),可以看出,光流計(jì)算結(jié)果在每個(gè)像素上都獲得一個(gè)速度矢量,可以顯示更加精細(xì)的流場(chǎng)結(jié)構(gòu),而PIV算法由于采用了8pixel×8pixel的判讀小區(qū),其流場(chǎng)結(jié)構(gòu)的精細(xì)顯示受到限制。

    圖28 PIV計(jì)算獲得的流場(chǎng)結(jié)構(gòu)

    圖29 光流計(jì)算獲得的流場(chǎng)結(jié)構(gòu)

    圖30 PIV結(jié)果中小尺度結(jié)構(gòu)區(qū)域放大

    圖31 光流結(jié)果中小尺度結(jié)構(gòu)區(qū)域放大

    5 結(jié)論與展望

    通過對(duì)光流算法研究與編程實(shí)現(xiàn),成功搭建了可用于光流測(cè)速的實(shí)驗(yàn)平臺(tái),進(jìn)行了小型風(fēng)洞內(nèi)的光流速度場(chǎng)測(cè)量實(shí)驗(yàn),比較光流算法和PIV對(duì)拍攝的粒子圖像處理得到的速度場(chǎng),可以得出結(jié)論:在本實(shí)驗(yàn)流動(dòng)條件下,光流算法可以在像素級(jí)分辨率情況下獲得優(yōu)于PIV獲得的速度矢量場(chǎng),適合用于復(fù)雜流動(dòng)的速度場(chǎng)測(cè)量。不過基于微分法的光流解算方法在計(jì)算較大位移區(qū)間時(shí)存在缺陷,未來將進(jìn)一步通過算法研究和實(shí)驗(yàn)手段進(jìn)行改進(jìn)。

    [1] Lang D B. Laser doppler velocity and vorticity measure- ments in turbulent shear layers[D]. California:Caltech, Pasadena, 1985.

    [2] Adrian R J. Particle-imaging techniques for experimental mechanics[J]. Annual Review of Fluid Mechanics, 1991, 23(1): 261-304.

    [3] Horn B K P, Schunck B G. Determining optical flow[J]. Artificial Intelligence, 1981, 17: 185-203.

    [4] Liu T, Shen L. Determination of velocity and skin friction fields from images by solving projected motion equations[C]. 22nd International Congress on Instrumenta- tion in Aerospace Simulation Facilities (ICIASF), International Congress on Instrumentation in Aerospace Simulation Facilities, Pacific Grove, CA, June 2007.

    [5] Liu T, Sullivan J P. Pressure and temperature sensitive paints experimental fluid mechanics[M]. Berlin: Springer-Verlag, 2004.

    [6] Tikhonov A N, Arsenin V Y. Solutions of ill-posed problems[M]. Winston, Washington D C, 1977.

    [7] Heeger D J. Model for the extraction of image flow[J]. Optical Society of America, 1987, A4: 1455-1471.

    [8] Fleet D J, Jespon A D. Computation of component image velocity from local phase information[J]. International Journal of Computer Vision, 1990, 5: 77- 104.

    [9] Liu Tianshu, Shen Lixin. Fluid flow and optical flow[J]. J. Fluid Mech, 2008, 614: 253-291.

    [10] Corpetti T, Heitz D, Arroyo G, et al. Fluid experimental flow estimation based on an optical flow scheme[J]. Exps Fluids, 2006, 40(1): 80-97.

    [11] Ruhnau P, Kohlberger T, Schnorr C, et al. Variational optical flow estimationfor particle image velocimetry[J]. Exps Fluids, 2005, 38(1): 21-32.

    [12] Zhan Huang, Hong Wei Wang, et al. The research of particle image velocimetry based on optical flow[C]. 16th International Symposium on Flow Visualization, Okinawa, Japan, 2004.

    (編輯:楊 娟)

    Research on particle image velocimetry based on optical flow

    Wang Hongwei, Huang Zhan*
    (China Academy of Aerospace Aerodynamics, Beijing 100074, China)

    As a new aerodynamics experimental technique, optical flow test technique gradually attracts more and more attention due to its advantages in vector field measurement with pixel scale resolution and strong smoothness ability. By means of scalar constraint equation combined with smoothness constraint condition, optical flow test technique can measure global velocity vector field with high space resolution. In this paper, the integration minimization optical flow velocimetry theory and algorithm were firstly studied and programmed by C++; then, in order to verify the accuracy of the optical flow algorithm, three groups of grayscale images shifted by given displacements were used for processing by the optical flow algorithm program and the result can be used to guide how to obtain high-precision optical flow calculation result in the actual measurement applications; at last, backward step verification experiment, in which tracer particle images were acquired by high speed camera and velocity vector field was calculated by optical flow algorithm, was completed in a small wind tunnel. The result calculated by optical flow algorithm was compared with the result calculated by PIV′s cross-correlation algorithm, which shows that optical flow test technique possesses the advantages of pixel scale resolution.

    optical flow;velcocity measurement;particle image;motion estimation;integration minimization

    1672-9897(2015)03-0068-08

    10.11729/syltlx20140115

    2014-10-11;

    2014-12-02

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2014CB744801)

    WangHW,HuangZ.Researchonparticleimagevelocimetrybasedonopticalflow.JournalofExperimentsinFluidMechanics, 2015, 29(3): 68-75. 王宏偉, 黃 湛. 基于光流算法的粒子圖像測(cè)速技術(shù)研究. 實(shí)驗(yàn)流體力學(xué), 2015, 29(3): 68-75.

    O355

    A

    王宏偉(1987-),男,黑龍江訥河人,助理工程師。研究方向:實(shí)驗(yàn)流體力學(xué)。通信地址:北京市豐臺(tái)區(qū)云崗西路17號(hào)(100074)。E-mail: www_whw_cn@163.com

    *通信作者 E-mail: xunfang05@sohu.com

    猜你喜歡
    光流流場(chǎng)灰度
    利用掩膜和單應(yīng)矩陣提高LK光流追蹤效果
    采用改進(jìn)導(dǎo)重法的拓?fù)浣Y(jié)構(gòu)灰度單元過濾技術(shù)
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場(chǎng)計(jì)算
    基于灰度拉伸的圖像水位識(shí)別方法研究
    基于物理學(xué)的改善粒子圖像測(cè)速穩(wěn)健光流方法研究
    轉(zhuǎn)杯紡排雜區(qū)流場(chǎng)與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場(chǎng)統(tǒng)計(jì)分析
    基于最大加權(quán)投影求解的彩色圖像灰度化對(duì)比度保留算法
    基于灰度線性建模的亞像素圖像抖動(dòng)量計(jì)算
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    极品教师在线免费播放| 亚洲国产精品合色在线| 又黄又爽又免费观看的视频| 久久人人精品亚洲av| 欧美黄色淫秽网站| 午夜精品一区二区三区免费看| 久久伊人香网站| 亚洲av成人精品一区久久| 香蕉久久夜色| 高潮久久久久久久久久久不卡| 欧美乱码精品一区二区三区| 国产精品98久久久久久宅男小说| 亚洲欧美日韩东京热| 女人爽到高潮嗷嗷叫在线视频| 在线播放国产精品三级| 可以在线观看毛片的网站| 国产成人精品久久二区二区免费| 亚洲精品在线美女| 午夜福利18| 欧美黑人欧美精品刺激| 麻豆av在线久日| 中文字幕精品亚洲无线码一区| 老司机福利观看| 欧美另类亚洲清纯唯美| 丁香六月欧美| av福利片在线观看| 欧美又色又爽又黄视频| 露出奶头的视频| 国模一区二区三区四区视频 | 精品国产美女av久久久久小说| 久久久久久国产a免费观看| av在线天堂中文字幕| 给我免费播放毛片高清在线观看| 精品午夜福利视频在线观看一区| 黑人操中国人逼视频| 好男人电影高清在线观看| 国产蜜桃级精品一区二区三区| 两个人视频免费观看高清| 后天国语完整版免费观看| 成人手机av| 国产精品 欧美亚洲| 中文资源天堂在线| 成人三级黄色视频| 国内毛片毛片毛片毛片毛片| 18禁裸乳无遮挡免费网站照片| 久久久国产成人精品二区| 99久久99久久久精品蜜桃| 精品国产亚洲在线| 两个人免费观看高清视频| 久久精品国产亚洲av高清一级| 欧美人与性动交α欧美精品济南到| 99热这里只有是精品50| 亚洲人成网站在线播放欧美日韩| 97碰自拍视频| 香蕉丝袜av| 国产熟女午夜一区二区三区| 熟女少妇亚洲综合色aaa.| 精品久久久久久久人妻蜜臀av| 久久国产精品影院| 国产精品久久久人人做人人爽| 亚洲人成77777在线视频| 日韩高清综合在线| 精品久久蜜臀av无| 亚洲va日本ⅴa欧美va伊人久久| 国产成人av激情在线播放| 777久久人妻少妇嫩草av网站| 无遮挡黄片免费观看| 黄色视频不卡| 久9热在线精品视频| 欧美国产日韩亚洲一区| 国产精品野战在线观看| xxxwww97欧美| av天堂在线播放| 亚洲一区二区三区色噜噜| 少妇粗大呻吟视频| 日韩成人在线观看一区二区三区| 一进一出好大好爽视频| 亚洲精品美女久久久久99蜜臀| 老汉色av国产亚洲站长工具| 久久久久久大精品| 又大又爽又粗| 久久精品成人免费网站| 女警被强在线播放| 一进一出抽搐动态| 亚洲欧美精品综合久久99| 后天国语完整版免费观看| 亚洲人与动物交配视频| 国产成人啪精品午夜网站| 色在线成人网| 五月伊人婷婷丁香| 国产高清视频在线观看网站| 五月玫瑰六月丁香| 欧美zozozo另类| 久久久国产成人免费| 免费在线观看黄色视频的| 国产v大片淫在线免费观看| 亚洲avbb在线观看| 亚洲av片天天在线观看| 国产精品一区二区三区四区免费观看 | 狂野欧美激情性xxxx| 色精品久久人妻99蜜桃| 90打野战视频偷拍视频| 天天躁狠狠躁夜夜躁狠狠躁| 欧美在线黄色| 国产不卡一卡二| 亚洲激情在线av| 亚洲精品久久成人aⅴ小说| 美女大奶头视频| 国产真人三级小视频在线观看| 麻豆成人av在线观看| 麻豆久久精品国产亚洲av| 老司机靠b影院| 老汉色∧v一级毛片| 中文亚洲av片在线观看爽| 老司机靠b影院| 中文字幕人成人乱码亚洲影| 国产不卡一卡二| 亚洲成人中文字幕在线播放| 精品人妻1区二区| 人人妻,人人澡人人爽秒播| 亚洲专区中文字幕在线| 九色国产91popny在线| 久99久视频精品免费| 夜夜夜夜夜久久久久| svipshipincom国产片| 日日爽夜夜爽网站| 日本熟妇午夜| 精品国产亚洲在线| 欧美成人午夜精品| 国产精品亚洲美女久久久| 99久久国产精品久久久| 国产精品日韩av在线免费观看| 国产爱豆传媒在线观看 | 90打野战视频偷拍视频| 国产精品av久久久久免费| 亚洲七黄色美女视频| 精品一区二区三区视频在线观看免费| 在线观看午夜福利视频| 男人舔女人的私密视频| 国产69精品久久久久777片 | 天堂av国产一区二区熟女人妻 | 亚洲七黄色美女视频| 一边摸一边做爽爽视频免费| 黄色片一级片一级黄色片| 亚洲国产精品sss在线观看| 国产av一区在线观看免费| 国产精品久久久久久精品电影| 久久这里只有精品19| 国产精品综合久久久久久久免费| 国产麻豆成人av免费视频| 久久久久久人人人人人| 国产视频内射| 此物有八面人人有两片| av片东京热男人的天堂| 午夜福利18| 国产高清videossex| 1024香蕉在线观看| 夜夜躁狠狠躁天天躁| 琪琪午夜伦伦电影理论片6080| 精品国产乱子伦一区二区三区| а√天堂www在线а√下载| 亚洲中文字幕日韩| 国产高清videossex| 99精品欧美一区二区三区四区| 啪啪无遮挡十八禁网站| 不卡一级毛片| 精品熟女少妇八av免费久了| 草草在线视频免费看| 国产成人精品久久二区二区免费| 国产精品久久视频播放| 女人高潮潮喷娇喘18禁视频| www.熟女人妻精品国产| 国产精品精品国产色婷婷| 午夜激情av网站| 国产视频内射| 亚洲av成人av| 色老头精品视频在线观看| 在线观看免费午夜福利视频| 亚洲乱码一区二区免费版| 国产高清videossex| 天天添夜夜摸| 亚洲欧美精品综合一区二区三区| 成人av在线播放网站| 国产主播在线观看一区二区| 久久久久精品国产欧美久久久| 国产欧美日韩精品亚洲av| 色老头精品视频在线观看| 亚洲欧美激情综合另类| 亚洲国产欧洲综合997久久,| 亚洲精品一卡2卡三卡4卡5卡| 黄色片一级片一级黄色片| 久久久久久久久久黄片| 久久精品国产综合久久久| 欧美色视频一区免费| 国产高清激情床上av| 欧美乱色亚洲激情| 精品无人区乱码1区二区| 深夜精品福利| 国产乱人伦免费视频| 亚洲第一欧美日韩一区二区三区| 丰满的人妻完整版| 窝窝影院91人妻| 国产一区二区三区视频了| 成人亚洲精品av一区二区| 欧美性猛交黑人性爽| 午夜两性在线视频| 亚洲第一电影网av| 国产精华一区二区三区| 美女高潮喷水抽搐中文字幕| 国产精品亚洲一级av第二区| 成熟少妇高潮喷水视频| 成人av在线播放网站| 国产成人影院久久av| 亚洲人与动物交配视频| 啪啪无遮挡十八禁网站| 久久久久久久久久黄片| 波多野结衣高清无吗| 99久久久亚洲精品蜜臀av| 亚洲真实伦在线观看| 欧美性长视频在线观看| 日日干狠狠操夜夜爽| 亚洲美女视频黄频| 99国产精品一区二区三区| 午夜成年电影在线免费观看| 国产成人精品无人区| 巨乳人妻的诱惑在线观看| а√天堂www在线а√下载| 嫁个100分男人电影在线观看| 欧美成人一区二区免费高清观看 | 99久久99久久久精品蜜桃| 亚洲成av人片在线播放无| 国产探花在线观看一区二区| 手机成人av网站| 国内精品久久久久精免费| 麻豆成人午夜福利视频| 在线观看美女被高潮喷水网站 | 变态另类丝袜制服| 午夜视频精品福利| 久久伊人香网站| 99久久精品热视频| 一级黄色大片毛片| 午夜成年电影在线免费观看| 舔av片在线| 曰老女人黄片| 国产熟女午夜一区二区三区| 亚洲,欧美精品.| 国内少妇人妻偷人精品xxx网站 | 97人妻精品一区二区三区麻豆| 99久久99久久久精品蜜桃| 日本撒尿小便嘘嘘汇集6| 18禁美女被吸乳视频| 男女下面进入的视频免费午夜| av福利片在线观看| 黄色丝袜av网址大全| 毛片女人毛片| 久久精品影院6| 人人妻人人澡欧美一区二区| 国内精品一区二区在线观看| 日韩大码丰满熟妇| 亚洲中文字幕一区二区三区有码在线看 | 亚洲熟妇熟女久久| 九九热线精品视视频播放| 国产日本99.免费观看| av福利片在线观看| 精品久久久久久久毛片微露脸| 成人国语在线视频| 动漫黄色视频在线观看| 制服丝袜大香蕉在线| 国产人伦9x9x在线观看| 日韩 欧美 亚洲 中文字幕| 后天国语完整版免费观看| 欧美日韩瑟瑟在线播放| 亚洲人成77777在线视频| 亚洲天堂国产精品一区在线| 亚洲国产精品sss在线观看| 一本久久中文字幕| 免费无遮挡裸体视频| 亚洲熟妇熟女久久| 久久久精品大字幕| 香蕉av资源在线| 国产97色在线日韩免费| 亚洲 国产 在线| 午夜免费激情av| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲五月天丁香| 亚洲午夜理论影院| 男女做爰动态图高潮gif福利片| 可以在线观看的亚洲视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| 女人被狂操c到高潮| 成人av一区二区三区在线看| 制服丝袜大香蕉在线| 午夜福利成人在线免费观看| 在线a可以看的网站| 色综合亚洲欧美另类图片| 国产欧美日韩精品亚洲av| 亚洲国产精品999在线| 美女免费视频网站| 九九热线精品视视频播放| 亚洲免费av在线视频| 日日摸夜夜添夜夜添小说| 久久精品国产99精品国产亚洲性色| cao死你这个sao货| 中国美女看黄片| 狂野欧美激情性xxxx| 亚洲最大成人中文| 亚洲中文字幕日韩| 国产成人啪精品午夜网站| 最近最新中文字幕大全电影3| av免费在线观看网站| 婷婷精品国产亚洲av在线| 国产一区二区三区在线臀色熟女| 久久久精品国产亚洲av高清涩受| 在线国产一区二区在线| 国产亚洲精品第一综合不卡| 女生性感内裤真人,穿戴方法视频| 午夜老司机福利片| 国产伦一二天堂av在线观看| 人人妻,人人澡人人爽秒播| 久久精品影院6| 好男人在线观看高清免费视频| 久久 成人 亚洲| 午夜福利18| 久久亚洲真实| 久久久久久久久免费视频了| 国产高清视频在线播放一区| 国产黄色小视频在线观看| 国产成人一区二区三区免费视频网站| 韩国av一区二区三区四区| 亚洲av成人av| 脱女人内裤的视频| 在线免费观看的www视频| a在线观看视频网站| 久久精品国产综合久久久| 五月伊人婷婷丁香| 久久这里只有精品19| 每晚都被弄得嗷嗷叫到高潮| 99在线视频只有这里精品首页| 在线十欧美十亚洲十日本专区| 久久久久久大精品| 久久久久久久久中文| 久久香蕉国产精品| 日日摸夜夜添夜夜添小说| 久久精品影院6| 亚洲人成网站高清观看| 国产精品久久久久久精品电影| 色av中文字幕| 搡老岳熟女国产| 99热6这里只有精品| 免费在线观看影片大全网站| 久9热在线精品视频| 精品国产亚洲在线| 国产欧美日韩一区二区精品| 国产午夜福利久久久久久| 国产精品久久视频播放| 国产精品98久久久久久宅男小说| 日本 欧美在线| 亚洲va日本ⅴa欧美va伊人久久| 久久中文字幕人妻熟女| 欧美日韩国产亚洲二区| 欧美性猛交╳xxx乱大交人| 成人精品一区二区免费| 久久久久久国产a免费观看| 欧美一级a爱片免费观看看 | 成人三级做爰电影| 国产在线观看jvid| 欧美中文日本在线观看视频| 美女大奶头视频| 免费在线观看成人毛片| 村上凉子中文字幕在线| 色综合欧美亚洲国产小说| 亚洲国产精品999在线| 狂野欧美激情性xxxx| 久久国产乱子伦精品免费另类| 国产视频内射| 国产亚洲精品第一综合不卡| 精品午夜福利视频在线观看一区| 97人妻精品一区二区三区麻豆| 国产区一区二久久| 国产野战对白在线观看| 成年免费大片在线观看| 欧美日本亚洲视频在线播放| 首页视频小说图片口味搜索| 精品久久久久久,| 性欧美人与动物交配| 久久国产乱子伦精品免费另类| 精品无人区乱码1区二区| 身体一侧抽搐| 国产欧美日韩一区二区精品| 欧美日韩黄片免| 午夜福利视频1000在线观看| 麻豆久久精品国产亚洲av| 日本一二三区视频观看| 不卡av一区二区三区| 五月玫瑰六月丁香| 高潮久久久久久久久久久不卡| 婷婷六月久久综合丁香| 国产成人一区二区三区免费视频网站| 三级国产精品欧美在线观看 | √禁漫天堂资源中文www| 一级a爱片免费观看的视频| 成人亚洲精品av一区二区| 国产在线观看jvid| 国产精品永久免费网站| 香蕉av资源在线| 91老司机精品| 亚洲熟妇熟女久久| 91在线观看av| 免费观看人在逋| 搞女人的毛片| 国产精品久久久久久人妻精品电影| 国产成+人综合+亚洲专区| 国产精品av久久久久免费| 天天一区二区日本电影三级| 人成视频在线观看免费观看| 怎么达到女性高潮| 久久久久久久久久黄片| 色综合欧美亚洲国产小说| avwww免费| 午夜福利成人在线免费观看| 他把我摸到了高潮在线观看| 国产成人av激情在线播放| 免费在线观看完整版高清| 一本精品99久久精品77| 亚洲18禁久久av| 中文字幕久久专区| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区免费欧美| www.999成人在线观看| 久久久久久亚洲精品国产蜜桃av| 色噜噜av男人的天堂激情| 亚洲欧洲精品一区二区精品久久久| 久久婷婷人人爽人人干人人爱| 亚洲电影在线观看av| 亚洲五月天丁香| 一本大道久久a久久精品| 亚洲成人国产一区在线观看| 午夜免费激情av| 91麻豆av在线| 亚洲全国av大片| cao死你这个sao货| 91九色精品人成在线观看| 在线观看午夜福利视频| 亚洲国产欧美网| 亚洲熟妇中文字幕五十中出| 欧美成人免费av一区二区三区| 精品少妇一区二区三区视频日本电影| 在线免费观看的www视频| 久久伊人香网站| 午夜福利18| 日韩有码中文字幕| 一级黄色大片毛片| 日日干狠狠操夜夜爽| 一本综合久久免费| 精品乱码久久久久久99久播| 一级毛片女人18水好多| 久久久久久国产a免费观看| 日韩大尺度精品在线看网址| 精品久久久久久久毛片微露脸| 欧美绝顶高潮抽搐喷水| av超薄肉色丝袜交足视频| av福利片在线观看| 制服诱惑二区| 女人高潮潮喷娇喘18禁视频| 精品欧美一区二区三区在线| 巨乳人妻的诱惑在线观看| 一本精品99久久精品77| 91成年电影在线观看| 国产亚洲精品一区二区www| 欧美黑人巨大hd| 国产精品99久久99久久久不卡| 黄片大片在线免费观看| 亚洲 欧美 日韩 在线 免费| 美女 人体艺术 gogo| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲第一电影网av| 免费高清视频大片| 女人爽到高潮嗷嗷叫在线视频| 亚洲aⅴ乱码一区二区在线播放 | 最近最新中文字幕大全免费视频| 国产亚洲av高清不卡| 精品欧美国产一区二区三| 久久久久亚洲av毛片大全| 精品日产1卡2卡| 国产私拍福利视频在线观看| 可以在线观看毛片的网站| 亚洲狠狠婷婷综合久久图片| 欧美日韩瑟瑟在线播放| 后天国语完整版免费观看| 婷婷精品国产亚洲av在线| 亚洲国产中文字幕在线视频| 国产精品久久久久久精品电影| 国产99久久九九免费精品| 亚洲人成77777在线视频| 亚洲国产精品合色在线| 又大又爽又粗| 国产高清视频在线播放一区| 久久精品国产99精品国产亚洲性色| 欧美午夜高清在线| 国产99久久九九免费精品| 国产伦人伦偷精品视频| 在线a可以看的网站| 亚洲人成电影免费在线| 此物有八面人人有两片| 欧美日韩亚洲国产一区二区在线观看| 少妇被粗大的猛进出69影院| 国产成人精品久久二区二区免费| 久久精品国产亚洲av高清一级| 在线观看www视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| 巨乳人妻的诱惑在线观看| 级片在线观看| 国产三级黄色录像| 好看av亚洲va欧美ⅴa在| 午夜福利在线观看吧| 国产一级毛片七仙女欲春2| 欧美黄色淫秽网站| 他把我摸到了高潮在线观看| 中文字幕av在线有码专区| 亚洲精品中文字幕一二三四区| 全区人妻精品视频| 最新在线观看一区二区三区| 日韩欧美一区二区三区在线观看| 丁香六月欧美| 久久精品国产亚洲av高清一级| 亚洲人成网站在线播放欧美日韩| 很黄的视频免费| 久久亚洲精品不卡| 国产99久久九九免费精品| 国产高清有码在线观看视频 | 欧美在线一区亚洲| 欧美日韩黄片免| 久久草成人影院| 宅男免费午夜| 好看av亚洲va欧美ⅴa在| a级毛片在线看网站| 悠悠久久av| 国产成人欧美在线观看| 十八禁网站免费在线| 国产亚洲精品一区二区www| 在线观看免费日韩欧美大片| 神马国产精品三级电影在线观看 | 久久香蕉国产精品| 国产精品一区二区三区四区免费观看 | 99久久综合精品五月天人人| 女人高潮潮喷娇喘18禁视频| 99久久国产精品久久久| 日本 av在线| 啦啦啦观看免费观看视频高清| 19禁男女啪啪无遮挡网站| 一级毛片精品| 国产伦人伦偷精品视频| 我的老师免费观看完整版| 首页视频小说图片口味搜索| 一二三四社区在线视频社区8| 国产免费男女视频| 一本一本综合久久| 久久国产乱子伦精品免费另类| 久久国产精品人妻蜜桃| 听说在线观看完整版免费高清| 变态另类成人亚洲欧美熟女| aaaaa片日本免费| 亚洲av美国av| 不卡av一区二区三区| 成人av在线播放网站| 男人舔女人下体高潮全视频| 啦啦啦韩国在线观看视频| 99国产精品一区二区蜜桃av| 免费在线观看影片大全网站| 欧美另类亚洲清纯唯美| 19禁男女啪啪无遮挡网站| 白带黄色成豆腐渣| 一本一本综合久久| 一区福利在线观看| 午夜精品一区二区三区免费看| 大型av网站在线播放| 亚洲专区中文字幕在线| 色av中文字幕| 久久久水蜜桃国产精品网| 国内精品一区二区在线观看| 久99久视频精品免费| 草草在线视频免费看| 国产亚洲精品综合一区在线观看 | 国产男靠女视频免费网站| 观看免费一级毛片| 在线观看美女被高潮喷水网站 | 午夜福利在线观看吧| av天堂在线播放| 中文在线观看免费www的网站 | 久久久久久免费高清国产稀缺| 中出人妻视频一区二区| 岛国视频午夜一区免费看| www日本在线高清视频| 可以在线观看的亚洲视频| 久久精品国产清高在天天线| 搡老熟女国产l中国老女人| 久9热在线精品视频| 国产主播在线观看一区二区| netflix在线观看网站| 禁无遮挡网站| 国产精品久久久av美女十八| 怎么达到女性高潮| а√天堂www在线а√下载| 午夜福利高清视频| 蜜桃久久精品国产亚洲av| 久久精品综合一区二区三区| 国产亚洲av高清不卡| 色综合欧美亚洲国产小说| 亚洲av日韩精品久久久久久密| 欧美成人午夜精品| 少妇的丰满在线观看| 国产一区二区激情短视频| 后天国语完整版免费观看| av在线播放免费不卡|