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

    粒子運(yùn)動(dòng)軌跡的圖像處理及流場(chǎng)重構(gòu)算法研究

    2020-01-10 01:53:20蔡小舒
    實(shí)驗(yàn)流體力學(xué) 2019年4期
    關(guān)鍵詞:插值骨架流場(chǎng)

    吳 凡, 周 騖, 蔡小舒

    (上海理工大學(xué) 顆粒與兩相流測(cè)量研究所,上海 200093)

    0 引 言

    流動(dòng)廣泛存在于化工和能源等行業(yè),復(fù)雜流動(dòng)現(xiàn)象中流動(dòng)參數(shù)的獲取通常十分重要[1- 2]。 目前有很多測(cè)量顆粒速度的方法,其中非侵入式的光學(xué)測(cè)量按照原理有兩大類,即基于多普勒原理的測(cè)量方法和基于TOF(Time of Flight)的測(cè)量方法[3]。 前者利用激光照射顆粒產(chǎn)生的多普勒頻移測(cè)速,后者則根據(jù)單位時(shí)間內(nèi)顆粒移動(dòng)的距離來計(jì)算速度[4]。

    典型的基于TOF的測(cè)量方法包括PIV(Particle Image Velocimetry)和PTV(Particle Tracking Veloci- metry)。PIV是20世紀(jì)70年代末發(fā)展起來的一種非接觸式流體力學(xué)測(cè)速方法,是目前應(yīng)用最為廣泛的全場(chǎng)無干擾測(cè)速手段[5]。PTV技術(shù)本質(zhì)上是PIV技術(shù)的衍生,兩者的不同之處在于:PIV獲取的速度來自于觀測(cè)區(qū)域內(nèi)多個(gè)粒子位移的統(tǒng)計(jì)平均值,而PTV則是通過追蹤單個(gè)粒子的位移來獲得速度。PTV由于其追蹤顆粒運(yùn)動(dòng)的特點(diǎn),所得結(jié)果具有直觀性和準(zhǔn)確性,不過同時(shí)也具有粒子濃度不宜過高的局限性[6]。早期PTV技術(shù)由于CCD攝像系統(tǒng)水平的

    限制,其時(shí)間和空間分辨率均不高,精度也不如PIV[7]。隨著計(jì)算機(jī)和CCD相機(jī)等硬件水平的不斷提高,在輸移現(xiàn)象和傳送現(xiàn)象等尤其關(guān)注顆粒的拉格朗日運(yùn)動(dòng)過程的情況下,PTV的作用逐漸顯現(xiàn)出來,因而也得到了廣泛的應(yīng)用[8]。

    根據(jù)圖像記錄的方式不同,PTV可以分為3種類型:?jiǎn)螏瑔蚊}沖、單幀多脈沖和多幀單脈沖[9]。近年來PTV領(lǐng)域的研究大多集中在多幀單脈沖形式[10]。多幀單脈沖形式具有更高的精度且能很好地解決速度二義性問題,因而得到廣泛應(yīng)用。不過當(dāng)涉及湍流這種十分復(fù)雜、瞬息萬變的流動(dòng)時(shí), 通常需要頻率非常高的激光脈沖, 才能實(shí)現(xiàn)對(duì)復(fù)雜流動(dòng)顆粒的精確追蹤, 但同時(shí)也會(huì)造成成本的急劇增加。

    單幀單曝光圖像法,即PTV的單幀單脈沖形式,具有系統(tǒng)簡(jiǎn)潔、圖像直觀等優(yōu)點(diǎn)。由于只有一次曝光,因此不需要功率非常大的激光器即可完成對(duì)示蹤粒子的追蹤。不過單幀單曝光圖像法的圖像處理過程相對(duì)復(fù)雜,不同于多幀單脈沖形式的PTV,單幀單曝光雖然沒有粒子匹配的問題,但由于其處理對(duì)象是軌跡而非單個(gè)點(diǎn),復(fù)雜度也有所提升。

    近年來對(duì)PTV的研究多集中在多幀單脈沖形式,對(duì)單幀單曝光圖像法的處理算法研究相對(duì)較少。其中,劉春嶸等[11]通過形心主軸的概念來計(jì)算運(yùn)動(dòng)軌跡的速度; 陳晶麗等[12]提出針對(duì)顆粒運(yùn)動(dòng)軌跡采用外接矩形的方法來提取顆粒速度與粒度;Jing等[13]則提出了通過檢測(cè)軌跡邊緣8個(gè)特征點(diǎn)來計(jì)算速度與粒度的方法。上述3種方法各有優(yōu)勢(shì)和特點(diǎn),不過三者有一個(gè)共同點(diǎn),那就是針對(duì)一條軌跡處理出一個(gè)速度矢量,導(dǎo)致信息的利用率較低。同時(shí),在處理彎曲軌跡時(shí),這些方法的精度均有所下降,而復(fù)雜流場(chǎng)中的彎曲軌跡是無法避免的。因此,本文提出一種方法,不僅能同時(shí)提取直線段軌跡和彎曲段軌跡的速度矢量,同時(shí)也提高軌跡上運(yùn)動(dòng)信息的利用率。

    1 單幀單曝光圖像法測(cè)速原理

    在利用單幀單曝光圖像法測(cè)量示蹤粒子的速度時(shí),一般CCD相機(jī)固定不動(dòng),通過視窗記錄測(cè)量區(qū)的圖像,此時(shí),離散相顆粒(如水滴、煤粉、氣泡等)與CCD相機(jī)發(fā)生相對(duì)位移, 其成像在曝光時(shí)間內(nèi)的疊加就形成了模糊的“拖影”,如圖1所示[14]。

    單幀單曝光運(yùn)動(dòng)模糊圖像包含了顆粒與相機(jī)的相對(duì)位移信息。因?yàn)槠毓鈺r(shí)間已知,在物體勻速直線運(yùn)動(dòng)假設(shè)的前提下,可以計(jì)算出運(yùn)動(dòng)速度,如式(1)所示。

    圖1 運(yùn)動(dòng)模糊圖像及其成像模型

    (1)

    式中,T為曝光時(shí)間且已知,故速度v的計(jì)算關(guān)鍵就在于運(yùn)動(dòng)長(zhǎng)度L的計(jì)算?,F(xiàn)有的很多算法都將軌跡近似視為一段直線,且最終得出一個(gè)速度矢量。以外接矩形法為例,該方法求出軌跡的最小外接矩形,提取矩形的長(zhǎng)和寬,長(zhǎng)度減去寬度便得到了顆粒的運(yùn)動(dòng)距離。其他幾種算法各有不同,但都具有一定的相似性。上述算法均需滿足顆粒運(yùn)動(dòng)軌跡接近直線的前提,在顆粒軌跡為曲線的情況下會(huì)產(chǎn)生一定誤差。所以需要提出另一種計(jì)算速度的方法,能同時(shí)適用于直線軌跡和彎曲軌跡的情況。

    2 基于骨架提取的顆粒軌跡處理算法

    骨架提取,也被稱為二值圖像細(xì)化。這種算法能將一個(gè)連通區(qū)域細(xì)化成一個(gè)像素的寬度,同時(shí)保持圖像的形態(tài)學(xué)特征,主要用于特征提取和目標(biāo)拓?fù)浔硎尽膱D2可以很直觀地看出骨架提取的特點(diǎn)。

    圖2 骨架提取示意圖

    結(jié)合骨架提取的特點(diǎn)可以將顆粒的軌跡進(jìn)行細(xì)化,得到單像素寬度的軌跡圖,也就是像素精度的跡線。在此基礎(chǔ)上進(jìn)行軌跡長(zhǎng)度計(jì)算要準(zhǔn)確得多,且對(duì)于彎曲軌跡也更能區(qū)分各像素點(diǎn)處速度方向的細(xì)微區(qū)別。為得到效果更好的軌跡骨架圖,需要對(duì)原圖依次進(jìn)行去噪、自適應(yīng)閾值分割、閉運(yùn)算、去除小顆粒及骨架提取等操作,主要的處理流程如圖3所示[15]。

    2.1 基于軌跡骨架圖的速度提取

    在原理介紹部分已經(jīng)提到,速度計(jì)算的前提是計(jì)算運(yùn)動(dòng)長(zhǎng)度。在獲得軌跡的單像素寬度跡線之后,通過連通區(qū)域分割函數(shù)可以識(shí)別出每條軌跡上的所有像素點(diǎn),傳統(tǒng)算法采用跡線起點(diǎn)與終點(diǎn)之間的距離作為長(zhǎng)度,這樣在處理彎曲軌跡時(shí)就會(huì)產(chǎn)生較大誤差。本文提出類似于求線積分的求長(zhǎng)度算法,首先將每2個(gè)像素點(diǎn)之間的距離求出來,然后迭加得到總的軌跡長(zhǎng)度,2種算法的對(duì)比如圖4所示。

    圖3 圖像處理主要流程[15]

    圖4 新舊算法求軌跡長(zhǎng)度的對(duì)比

    Fig.4Thedifferenceontrajectorylengthsbetweenthetraditionalandthenewalgorithm

    對(duì)于速度的方向,以往的計(jì)算方式是采取軌跡首尾相連的一條直線作為方向向量,如圖4中的舊算法。但這種方法在處理彎曲軌跡的時(shí)候會(huì)產(chǎn)生很大誤差,且一條軌跡上所有的點(diǎn)共享一個(gè)速度矢量,十分籠統(tǒng)。因此,提出了求跡線切線的方法,將像素點(diǎn)前后若干個(gè)像素點(diǎn)之間的割線近似視為切線,并作為軌跡上該點(diǎn)的速度方向,不僅對(duì)于彎曲軌跡能很好地提高精度,而且同一條軌跡上各個(gè)像素點(diǎn)的速度矢量也有區(qū)分,能很好地體現(xiàn)顆粒的實(shí)際運(yùn)動(dòng)過程。圖5顯示了新舊算法求速度方向的區(qū)別,新算法(左側(cè))各個(gè)像素點(diǎn)處的速度由其切線決定,能更好地顯示粒子的流動(dòng)過程,舊算法(右側(cè))則過于籠統(tǒng),所有像素點(diǎn)的速度矢量完全相同,不能體現(xiàn)實(shí)際的流動(dòng)狀態(tài)。此外,若將軌跡上各個(gè)點(diǎn)處的速度矢量逐幀播放,也能更清晰生動(dòng)地顯示顆粒的動(dòng)態(tài)流動(dòng)過程。

    圖5 新舊算法求速度方向的對(duì)比

    Fig.5Thedifferenceonvelocitydirectionbetweenthetraditionalandthenewalgorithm

    圖6為新舊算法在計(jì)算不同彎曲程度軌跡時(shí)的精度對(duì)比。從仿真圓形軌跡上截取不同弧長(zhǎng)的軌跡,然后分別采用2種算法求長(zhǎng),并與弧長(zhǎng)的理論值進(jìn)行對(duì)比。圖中,橫坐標(biāo)為軌跡的弧長(zhǎng),也代表了軌跡的彎曲程度,縱坐標(biāo)為軌跡長(zhǎng)度。紅色圓點(diǎn)為舊算法計(jì)算的長(zhǎng)度,藍(lán)色三角點(diǎn)為新算法計(jì)算的長(zhǎng)度,黑色方點(diǎn)為實(shí)際弧長(zhǎng)。由圖6可見,在弧長(zhǎng)較低,即軌跡接近直線的情況下,2種算法的結(jié)果均十分接近真值;隨著弧長(zhǎng)的增加(即軌跡彎曲程度增加),新算法則逐漸體現(xiàn)出了優(yōu)勢(shì)。

    圖6 2種速度提取算法精度比較

    Fig.6Accuracycontrastbetweenthetraditionalandthenewalgorithm

    2.2 速度二義性的判斷

    在前文的敘述中,速度的計(jì)算忽略了一個(gè)很重要的問題,就是速度二義性的判斷。對(duì)一條軌跡來說,一般有2個(gè)端點(diǎn),那么軌跡是從端點(diǎn)A運(yùn)動(dòng)到B還是從B運(yùn)動(dòng)到A是2種完全不同的情況,所以需要對(duì)該問題進(jìn)行判斷。下面提供2種方法:一種是利用流體的角速度進(jìn)行判斷,另一種則利用多幀圖像之間的運(yùn)動(dòng)關(guān)系進(jìn)行判斷。

    2.2.1 角速度判別法

    角速度判別法是一種利用右手螺旋定則來確定軌跡運(yùn)動(dòng)方向的方法。如圖7所示,有一組圍繞同一中心點(diǎn)進(jìn)行旋轉(zhuǎn)的4條軌跡,此時(shí)執(zhí)行前述的速度提取流程后,上下2條軌跡的速度方向會(huì)保持一致(水平向右)。但實(shí)際上,由于這是一組旋轉(zhuǎn)軌跡,上下2條軌跡的方向應(yīng)該相反。為使同一旋轉(zhuǎn)中心的一組軌跡旋轉(zhuǎn)方向保持一致,使用角速度判別法進(jìn)行判斷。圖7中,通過右手螺旋定則可以知道上方軌跡的角速度垂直紙面朝里,而下方的角速度垂直紙面向外,二者互相矛盾,通過改變其中一條軌跡的速度矢量即可解決速度二義性問題。在實(shí)際執(zhí)行過程中可預(yù)先定義軌跡的旋轉(zhuǎn)方向(順時(shí)針或逆時(shí)針),從而統(tǒng)一各個(gè)軌跡的速度方向,消除二義性。

    圖7 角速度判據(jù)

    對(duì)于直線運(yùn)動(dòng)狀態(tài),可以將旋轉(zhuǎn)中心點(diǎn)選在垂直于軌跡的無窮遠(yuǎn)處。此外,對(duì)于復(fù)雜的流動(dòng),可以通過選取多個(gè)旋轉(zhuǎn)中心分別進(jìn)行角速度判別,由于過程比較復(fù)雜,故方法具有一定的局限性。

    2.2.2 多幀連續(xù)曝光圖像匹配

    單幀圖像很難解決二義性問題,而多幀圖像則可以。因此,在相機(jī)進(jìn)行第一段曝光之后,立刻進(jìn)行第二段曝光,這樣得到2幀在時(shí)間序列上幾乎銜接的圖片,由于2幀圖片的時(shí)間先后順序已知,就能很方便地求出顆粒的速度方向。

    在多幀圖像匹配方法中,關(guān)鍵在于2幀圖片中同一條軌跡的匹配識(shí)別[16]。由于2次曝光之間的時(shí)間間隔非常短暫,所以2段軌跡理論上是能首尾相連的。實(shí)際處理中,對(duì)于第一幀圖片里的所有軌跡,于2個(gè)端點(diǎn)附近搜索其在第二幀圖片中是否有軌跡,搜索半徑由2段曝光之間的間距時(shí)間以及流場(chǎng)的速度決定。為了保證在搜索半徑中只有該軌跡的下一幀軌跡出現(xiàn),因此顆粒密度不宜太大。此外,考慮到邊緣部分的軌跡可能由于進(jìn)出視場(chǎng)而導(dǎo)致無法匹配,故應(yīng)提前將邊緣部分的軌跡剔除。圖8展示了2幀連續(xù)曝光圖片(曝光時(shí)間為7ms,間隔時(shí)間為0.077ms)匹配并疊加的過程,其軌跡匹配率為84.5%。

    圖8 2幀圖片的匹配與疊加過程

    在完成2幀圖片之間的軌跡匹配之后,速度二義性就解決了。舉例來說:假設(shè)AB是第一幀圖片的某條軌跡,端點(diǎn)分別是A和B,而CD是該軌跡在第二幀圖片上的位置,端點(diǎn)分別是C和D;若B和C是2條軌跡的連接點(diǎn),則軌跡的方向就依次是A—B—C—D。

    多幀圖像匹配之后,不僅解決了速度二義性問題,同時(shí),由于得到了不同時(shí)間序列上的速度值,也能因此獲得加速度。

    2.3 二值化與骨架提取的誤差分析

    基于骨架提取的軌跡處理算法可以分成2部分:圖像處理部分和速度計(jì)算部分。速度計(jì)算部分的誤差已經(jīng)在2.1節(jié)進(jìn)行了描述,此處對(duì)二值化和骨架提取的誤差進(jìn)行分析。

    2.3.1 二值化過程誤差分析

    單閾值分割通常采用Ostu (最大類間方差)閾值,其通過圖像灰度將圖像分為背景和目標(biāo),通過最大化背景和目標(biāo)之間的類間方差來找尋最優(yōu)閾值??紤]到單閾值分割在灰度分布不均的情況下并不適用,故使用了自適應(yīng)閾值分割。這是一種多閾值分割方法,通過分析像素周圍局部范圍的灰度特性來決定該像素點(diǎn)的閾值,能夠使閾值分割后的圖片保留更多信息。

    二值化過程的誤差分析實(shí)驗(yàn)采用了仿真的運(yùn)動(dòng)模糊圖片,如圖9所示。根據(jù)實(shí)際射流實(shí)驗(yàn)圖片,將仿真顆粒圖像的背景灰度設(shè)為40,目標(biāo)灰度設(shè)為200,粒徑為5pixel,將顆粒圖像進(jìn)行卷積并添加方差為0.008的高斯加性噪聲即得到仿真運(yùn)動(dòng)模糊圖像,卷積核則由運(yùn)動(dòng)長(zhǎng)度決定。

    圖9 多種運(yùn)動(dòng)的仿真運(yùn)動(dòng)模糊圖片

    仿真實(shí)驗(yàn)的結(jié)果如圖10所示,可見隨著運(yùn)動(dòng)長(zhǎng)度的增加,2種閾值分割的誤差均顯著下降,且自適應(yīng)閾值分割的誤差要明顯低于Ostu 閾值分割。

    圖10 不同二值化過程的誤差

    2.3.2 骨架提取過程誤差分析

    骨架提取過程將二值圖轉(zhuǎn)變成骨架圖,故誤差分析過程采用仿真的軌跡二值圖(由一系列圓形疊加形成的類似于橢圓的圖形)進(jìn)行一系列仿真實(shí)驗(yàn)。初步認(rèn)為骨架提取的結(jié)果受軌跡大小、軌跡長(zhǎng)短軸比、軌跡方向(位置)三方面的影響。實(shí)驗(yàn)過程中分別計(jì)算二值圖與骨架圖中軌跡對(duì)應(yīng)的運(yùn)動(dòng)長(zhǎng)度。在二值圖中,運(yùn)動(dòng)長(zhǎng)度為軌跡的長(zhǎng)度減去寬度;在骨架圖中,軌跡長(zhǎng)度即為像素距離迭代求和的長(zhǎng)度。仿真實(shí)驗(yàn)結(jié)果如圖11~13所示。

    從圖11可以看出骨架提取結(jié)果與軌跡的方向無關(guān);從圖12和13看出,骨架提取操作對(duì)運(yùn)動(dòng)長(zhǎng)度的影響很小。進(jìn)一步考察原始數(shù)據(jù)發(fā)現(xiàn),在所有情況下,骨架圖中的運(yùn)動(dòng)長(zhǎng)度都比二值圖的運(yùn)動(dòng)長(zhǎng)度大一

    個(gè)固定值1,這是由于長(zhǎng)度計(jì)量的起點(diǎn)不同所導(dǎo)致的。舉例來說,假設(shè)有一個(gè)僅包含一個(gè)像素點(diǎn)的軌跡,其長(zhǎng)度應(yīng)該為1,但長(zhǎng)度減去寬度卻等于0,從而導(dǎo)致了二者相差1的現(xiàn)象,這種情況在實(shí)際應(yīng)用過程中統(tǒng)一即可。

    圖11 骨架提取結(jié)果與軌跡方向的關(guān)系

    圖12 不同軌跡尺寸下二值圖與骨架圖的運(yùn)動(dòng)長(zhǎng)度

    Fig.12Motionlengthofbinaryandskeletonimageondifferenttrajectorysizes

    圖13 不同長(zhǎng)寬比(φ)下二值圖與骨架圖的運(yùn)動(dòng)長(zhǎng)度

    Fig.13Motionlengthofbinaryandskeletonimagewithdifferentwidth/lengthrate(φ)

    3 流場(chǎng)重構(gòu)算法

    在執(zhí)行了上述軌跡處理算法之后,可以得到每條軌跡的速度信息。然而多數(shù)情況下實(shí)際需要的是整個(gè)流場(chǎng)的速度信息,因此需要進(jìn)行流場(chǎng)重構(gòu)工作。

    3.1 RBF插值法重構(gòu)流體速度場(chǎng)

    從隨機(jī)的數(shù)據(jù)擬合出整個(gè)平面場(chǎng)的信息屬于數(shù)學(xué)上插值的范疇。選用徑向基函數(shù)(Radial Basis Function,RBF)插值的方法得出流場(chǎng)的速度場(chǎng)。該插值方法的基本假設(shè)是插值點(diǎn)的函數(shù)值受插值源函數(shù)值影響,且影響程度與插值點(diǎn)和插值源之間的距離有關(guān),其數(shù)學(xué)表達(dá)式如下:

    (2)

    根據(jù)不同應(yīng)用實(shí)例,可以選取不同的距離系數(shù)函數(shù)Φ以達(dá)到最佳的插值效果。比較多種函數(shù)的插值效果,從中選取名為“multiquadric” 的距離系數(shù),其距離系數(shù)函數(shù)如下:

    (3)

    考慮到流體速度場(chǎng)的特點(diǎn),任意一點(diǎn)的速度受來流方向速度值的影響大于下游方向速度值的影響。需要在RBF插值的基礎(chǔ)上進(jìn)行類似于“迎風(fēng)插值”的優(yōu)化。優(yōu)化的基本思想是在插值點(diǎn)的上游增加對(duì)插值點(diǎn)的影響;在插值點(diǎn)的下游減弱對(duì)插值點(diǎn)的影響。因此,可以構(gòu)建迎風(fēng)修正系數(shù)函數(shù)ψ,ψ的值由插值點(diǎn)與已知點(diǎn)之間速度向量的夾角決定,當(dāng)夾角為銳角時(shí),ψ大于0,對(duì)結(jié)果進(jìn)行增強(qiáng);當(dāng)夾角為鈍角時(shí),ψ小于0,對(duì)結(jié)果削弱。其表達(dá)式如下:

    (4)

    修正后的插值表達(dá)式如下:

    (5)

    3.2 RBF插值的誤差分析

    對(duì)插值的誤差進(jìn)行分析,采用仿真輸出的流場(chǎng)結(jié)果,隨機(jī)選取一定的數(shù)據(jù)點(diǎn)之后進(jìn)行插值,分析插值結(jié)果與原流場(chǎng)的差別即可。原流場(chǎng)與插值流場(chǎng)的速度云圖如圖14所示。

    為了量化2幅圖之間的相似度,以便更直觀地獲知不同插值方法的優(yōu)劣,定義2幅圖像之間的平均距離δ,如式(6)所示。δ越大,表示兩者之間的差別越大,δ越小,表示兩者越為接近。

    (6)

    表1顯示了選用不同插值函數(shù)并采用“迎風(fēng)插值”優(yōu)化之后的插值結(jié)果的平均距離。從中可以看出基于 “multiquadric” 并采用“迎風(fēng)插值”優(yōu)化之后得到插值效果是最好的,與肉眼直觀相符。

    圖14 仿真流場(chǎng)(a)與插值流場(chǎng)(b)對(duì)比

    Fig.14Contrastbetweensimulationflowfield(a)andinterpolationflowfield(b)

    表1 插值結(jié)果誤差分析Table 1 Error analysis on interpolation

    4 處理結(jié)果

    為驗(yàn)證本文算法在實(shí)際情況下的適用性,采用射流實(shí)驗(yàn)中拍攝的圖片作為處理對(duì)象,如圖15所示。實(shí)驗(yàn)中射流出口流速為0.15m/s,拍攝部位為射流口上部的卷吸區(qū)域。射流口位于圖片右下方,視場(chǎng)大小為2.88mm×3.86mm。圖16展示了軌跡骨架圖,圖17和18分別是最終速度場(chǎng)的云圖和矢量圖。

    圖15 原始圖像

    在速度云圖中,每個(gè)像素點(diǎn)處的顏色代表了該處的速度大小,因而能直觀地看出速度場(chǎng)在空間上的能量分布;在速度矢量圖中,箭頭方向?yàn)樗俣确较?,箭頭長(zhǎng)度表征了速度大小,從而可以看出流場(chǎng)整體運(yùn)動(dòng)趨勢(shì)。

    圖16 軌跡骨架圖

    圖17 速度云圖

    圖18 速度矢量圖

    為驗(yàn)證實(shí)驗(yàn)結(jié)果的可信度,利用插值后的速度場(chǎng)來計(jì)算速度的散度場(chǎng)??紤]到此處速度場(chǎng)并非連續(xù),因此采用變分方式來計(jì)算散度,同時(shí)假定流場(chǎng)為二維流動(dòng),散度計(jì)算公式如下:

    (7)

    采用云圖方式展示計(jì)算的散度場(chǎng),如圖19所示,由于常溫常壓下水流可視為不可壓流動(dòng),速度的散度應(yīng)處處為零,圖19中大多數(shù)散度值接近為零,少數(shù)幾處出現(xiàn)稍大的散度值(正值或負(fù)值),不過總體上結(jié)果符合事實(shí)。

    圖19 速度的散度云圖

    5 結(jié) 論

    為獲取顆粒的速度分布,針對(duì)單幀單曝光圖像的處理提出一種顆粒軌跡處理算法。通過骨架提取獲得軌跡的跡線,然后對(duì)跡線的像素坐標(biāo)進(jìn)行迭代求和來計(jì)算軌跡長(zhǎng)度,同時(shí)通過求跡線的切線來計(jì)算速度方向。相比傳統(tǒng)算法,該算法在速度大小和方向上均有更高精度,且不同于傳統(tǒng)算法一條軌跡一個(gè)速度矢量的模式,該算法的一條軌跡可以獲得軌跡上任意點(diǎn)處的速度矢量,大大提高了信息利用率。

    在得到原始的離散速度場(chǎng)之后,利用迎風(fēng)優(yōu)化后的RBF 插值算法,可以重構(gòu)一個(gè)相對(duì)連續(xù)且平滑的全分布流場(chǎng),可在此基礎(chǔ)上計(jì)算加速度、散度等參數(shù),方便后續(xù)研究;也可以很方便地與CFD 等仿真模擬結(jié)果進(jìn)行比較。

    本算法處理結(jié)果的精度上限取決于閾值分割過程,因而需要在拍攝過程中調(diào)整曝光時(shí)間、光源、顆粒濃度等因素,盡可能提高所拍攝圖像中軌跡的清晰度,以期達(dá)到好的閾值分割效果,減少圖片中的干擾軌跡,從而獲得更高精度的結(jié)果。

    猜你喜歡
    插值骨架流場(chǎng)
    淺談管狀骨架噴涂方法
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場(chǎng)計(jì)算
    骨架密度對(duì)炭/炭多孔骨架壓力浸滲銅的影響
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    轉(zhuǎn)杯紡排雜區(qū)流場(chǎng)與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場(chǎng)統(tǒng)計(jì)分析
    一種改進(jìn)FFT多譜線插值諧波分析方法
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    內(nèi)支撐骨架封抽技術(shù)在突出煤層瓦斯抽采中的應(yīng)用
    中國煤層氣(2014年3期)2014-08-07 03:07:45
    免费观看的影片在线观看| 中文资源天堂在线| 美女大奶头视频| 免费黄频网站在线观看国产| 激情五月婷婷亚洲| 人人妻人人看人人澡| 99九九线精品视频在线观看视频| 麻豆久久精品国产亚洲av| 日韩不卡一区二区三区视频在线| 午夜免费观看性视频| 婷婷色综合www| 日日摸夜夜添夜夜爱| av在线天堂中文字幕| 欧美日韩精品成人综合77777| 最近最新中文字幕大全电影3| 欧美成人一区二区免费高清观看| 男女下面进入的视频免费午夜| 看黄色毛片网站| 亚洲欧美一区二区三区国产| 国语对白做爰xxxⅹ性视频网站| 深夜a级毛片| 免费无遮挡裸体视频| av专区在线播放| 99热这里只有是精品在线观看| 一本一本综合久久| 亚洲av免费在线观看| 欧美日韩在线观看h| 两个人的视频大全免费| 国产男女超爽视频在线观看| 亚洲精华国产精华液的使用体验| 成人二区视频| 亚洲综合精品二区| 国产欧美另类精品又又久久亚洲欧美| 日本免费在线观看一区| 欧美激情久久久久久爽电影| 91精品一卡2卡3卡4卡| 成人亚洲精品av一区二区| 亚洲国产色片| 亚洲精品日韩av片在线观看| 国产有黄有色有爽视频| 中文字幕亚洲精品专区| 免费观看无遮挡的男女| 亚洲精品自拍成人| 最近视频中文字幕2019在线8| 国产午夜精品一二区理论片| 最近手机中文字幕大全| 国产一级毛片七仙女欲春2| 国产黄色视频一区二区在线观看| 午夜福利视频精品| 日韩中字成人| 精品国产露脸久久av麻豆 | 高清在线视频一区二区三区| 国产一区亚洲一区在线观看| 91午夜精品亚洲一区二区三区| 免费av不卡在线播放| 日韩一区二区视频免费看| 日日啪夜夜爽| 国产精品嫩草影院av在线观看| 色综合站精品国产| 91av网一区二区| 一级毛片黄色毛片免费观看视频| 99视频精品全部免费 在线| 久久久久免费精品人妻一区二区| 别揉我奶头 嗯啊视频| 亚洲国产最新在线播放| 婷婷六月久久综合丁香| 精品久久久久久久久久久久久| 亚洲精品国产av成人精品| 亚州av有码| 久久99蜜桃精品久久| 大话2 男鬼变身卡| 午夜视频国产福利| videos熟女内射| 国产v大片淫在线免费观看| 久久久久精品性色| 国产免费视频播放在线视频 | 国产精品久久久久久av不卡| 国产一区亚洲一区在线观看| 国产 亚洲一区二区三区 | 一区二区三区高清视频在线| 老司机影院毛片| 精品国产一区二区三区久久久樱花 | 一边亲一边摸免费视频| 免费播放大片免费观看视频在线观看| 亚洲国产高清在线一区二区三| 国产午夜精品论理片| 国产v大片淫在线免费观看| 丝袜美腿在线中文| 成人漫画全彩无遮挡| 日韩强制内射视频| av在线播放精品| 成人特级av手机在线观看| 久久99热这里只有精品18| 亚洲怡红院男人天堂| 国产 一区精品| 男女啪啪激烈高潮av片| 一级毛片黄色毛片免费观看视频| 国产高清不卡午夜福利| 国产伦一二天堂av在线观看| 波多野结衣巨乳人妻| 日本一本二区三区精品| 国产在线一区二区三区精| 晚上一个人看的免费电影| 久久久久性生活片| 啦啦啦啦在线视频资源| 99热这里只有是精品在线观看| 免费黄色在线免费观看| 麻豆成人午夜福利视频| 春色校园在线视频观看| 偷拍熟女少妇极品色| 国产综合精华液| 嫩草影院新地址| 亚洲国产精品成人综合色| 校园人妻丝袜中文字幕| 一级二级三级毛片免费看| 国产成人精品婷婷| 我的女老师完整版在线观看| 久久精品久久久久久噜噜老黄| 五月伊人婷婷丁香| 国产黄频视频在线观看| 精华霜和精华液先用哪个| 亚洲激情五月婷婷啪啪| 日韩人妻高清精品专区| 中文天堂在线官网| 熟女电影av网| 你懂的网址亚洲精品在线观看| 成人欧美大片| 99久久精品一区二区三区| 联通29元200g的流量卡| 国产成人a区在线观看| 在线观看人妻少妇| 久久久久久久久久久丰满| 中文字幕亚洲精品专区| 亚洲精品456在线播放app| 亚洲欧美日韩卡通动漫| 亚洲精品,欧美精品| 欧美日韩在线观看h| 热99在线观看视频| 日韩精品有码人妻一区| 日日撸夜夜添| 床上黄色一级片| 国产黄a三级三级三级人| 一级a做视频免费观看| 久久久国产一区二区| 日韩av免费高清视频| 一边亲一边摸免费视频| 听说在线观看完整版免费高清| 亚洲精品国产成人久久av| 亚洲国产欧美在线一区| 国产成人精品婷婷| 欧美另类一区| 噜噜噜噜噜久久久久久91| 日本爱情动作片www.在线观看| 天天躁日日操中文字幕| 久久久色成人| 三级国产精品片| 日韩一区二区三区影片| 成人二区视频| 国产在视频线在精品| 免费黄网站久久成人精品| 内地一区二区视频在线| 日本一本二区三区精品| 欧美日本视频| 国产成人freesex在线| 精品国产一区二区三区久久久樱花 | 久久鲁丝午夜福利片| 熟妇人妻久久中文字幕3abv| 插逼视频在线观看| ponron亚洲| 婷婷色综合大香蕉| 爱豆传媒免费全集在线观看| 男女边吃奶边做爰视频| 搡老妇女老女人老熟妇| 性插视频无遮挡在线免费观看| 男人和女人高潮做爰伦理| 亚洲精品成人av观看孕妇| 国产精品久久久久久久久免| 国产色爽女视频免费观看| 欧美高清成人免费视频www| 伊人久久精品亚洲午夜| 久久久久国产网址| 日韩,欧美,国产一区二区三区| 国产男女超爽视频在线观看| 日韩欧美 国产精品| 亚洲精品视频女| 亚洲经典国产精华液单| 天堂av国产一区二区熟女人妻| 91狼人影院| 人体艺术视频欧美日本| 国产永久视频网站| 黄片无遮挡物在线观看| 99久久精品热视频| 亚州av有码| 激情 狠狠 欧美| 丰满乱子伦码专区| 成人午夜高清在线视频| 欧美成人一区二区免费高清观看| 美女内射精品一级片tv| 国产成人aa在线观看| 熟妇人妻不卡中文字幕| 国产综合懂色| 人妻制服诱惑在线中文字幕| 精品国产一区二区三区久久久樱花 | 日韩成人伦理影院| 精品少妇黑人巨大在线播放| 波野结衣二区三区在线| 亚洲av日韩在线播放| 汤姆久久久久久久影院中文字幕 | 床上黄色一级片| 亚洲国产精品成人综合色| 中文乱码字字幕精品一区二区三区 | 精品亚洲乱码少妇综合久久| 久久这里有精品视频免费| 久久亚洲国产成人精品v| 我的老师免费观看完整版| eeuss影院久久| 国产精品国产三级国产av玫瑰| 日本wwww免费看| 中国国产av一级| 久久久久九九精品影院| 国产精品美女特级片免费视频播放器| 欧美 日韩 精品 国产| 99久久人妻综合| 联通29元200g的流量卡| 国产男女超爽视频在线观看| 免费看a级黄色片| 又粗又硬又长又爽又黄的视频| 欧美性猛交╳xxx乱大交人| 亚洲最大成人手机在线| 中文字幕亚洲精品专区| 国产伦一二天堂av在线观看| 亚洲国产高清在线一区二区三| 国产亚洲5aaaaa淫片| 色吧在线观看| 国产高清有码在线观看视频| 欧美xxxx性猛交bbbb| 久久久久久久亚洲中文字幕| 黄色一级大片看看| 亚洲精品亚洲一区二区| 久久99热这里只有精品18| 淫秽高清视频在线观看| 偷拍熟女少妇极品色| 一区二区三区免费毛片| 女人被狂操c到高潮| 精品人妻一区二区三区麻豆| 精品久久久久久久人妻蜜臀av| 精品国内亚洲2022精品成人| 国产又色又爽无遮挡免| 欧美最新免费一区二区三区| 2021天堂中文幕一二区在线观| 看非洲黑人一级黄片| 精品久久久久久久久av| 18禁裸乳无遮挡免费网站照片| 成年免费大片在线观看| 欧美3d第一页| 成人一区二区视频在线观看| 91狼人影院| 22中文网久久字幕| 日本wwww免费看| 18禁在线播放成人免费| 午夜爱爱视频在线播放| 啦啦啦啦在线视频资源| 日韩一区二区三区影片| 亚洲一区高清亚洲精品| 啦啦啦韩国在线观看视频| 国产精品99久久久久久久久| 午夜福利高清视频| 国产 亚洲一区二区三区 | 国国产精品蜜臀av免费| 亚洲欧美清纯卡通| 国产一级毛片在线| 夜夜爽夜夜爽视频| 欧美激情久久久久久爽电影| 国产三级在线视频| 国产视频首页在线观看| 国产综合精华液| 美女cb高潮喷水在线观看| 亚洲精品一区蜜桃| 久久久久久久久久成人| 神马国产精品三级电影在线观看| 99久久精品热视频| 一区二区三区四区激情视频| 久久精品夜色国产| 免费黄频网站在线观看国产| 岛国毛片在线播放| 亚洲精品日韩在线中文字幕| 97超视频在线观看视频| 亚洲高清免费不卡视频| 午夜激情福利司机影院| 成人一区二区视频在线观看| 欧美日韩综合久久久久久| 人人妻人人澡欧美一区二区| 欧美精品国产亚洲| 日本一二三区视频观看| 久久久精品免费免费高清| 久久久久免费精品人妻一区二区| 久久久久久久久中文| 亚洲第一区二区三区不卡| 久久久亚洲精品成人影院| 亚洲av成人精品一二三区| 最近中文字幕高清免费大全6| 最近的中文字幕免费完整| 青春草视频在线免费观看| 欧美三级亚洲精品| 精品久久久久久电影网| 丝袜美腿在线中文| 又爽又黄a免费视频| 免费看不卡的av| 成人毛片a级毛片在线播放| 国产黄a三级三级三级人| 好男人视频免费观看在线| 晚上一个人看的免费电影| 亚洲av中文字字幕乱码综合| 国产免费福利视频在线观看| 99久久精品热视频| 一区二区三区免费毛片| 成人鲁丝片一二三区免费| 91精品一卡2卡3卡4卡| 2021天堂中文幕一二区在线观| 激情五月婷婷亚洲| 插阴视频在线观看视频| 国产亚洲av片在线观看秒播厂 | 在线观看av片永久免费下载| 一级毛片电影观看| 非洲黑人性xxxx精品又粗又长| 秋霞伦理黄片| 国产成人精品久久久久久| 国内揄拍国产精品人妻在线| 一区二区三区四区激情视频| 如何舔出高潮| 国产亚洲精品av在线| 国产永久视频网站| 国产成年人精品一区二区| 色播亚洲综合网| 女的被弄到高潮叫床怎么办| 欧美激情久久久久久爽电影| 日韩精品有码人妻一区| 美女xxoo啪啪120秒动态图| 国产69精品久久久久777片| 国产精品不卡视频一区二区| 91午夜精品亚洲一区二区三区| 一级片'在线观看视频| 2021少妇久久久久久久久久久| 久久久久久久亚洲中文字幕| 成人毛片60女人毛片免费| 国产成人精品婷婷| 久久精品夜色国产| 尾随美女入室| 婷婷色综合大香蕉| 国产成人a∨麻豆精品| 麻豆av噜噜一区二区三区| 天堂av国产一区二区熟女人妻| 久久久久久九九精品二区国产| 夫妻午夜视频| 能在线免费观看的黄片| 午夜精品一区二区三区免费看| 国产老妇女一区| 亚洲av日韩在线播放| 日韩欧美三级三区| 欧美xxxx黑人xx丫x性爽| 成年人午夜在线观看视频 | 欧美变态另类bdsm刘玥| 亚洲精品视频女| 欧美性感艳星| 简卡轻食公司| 精品久久久噜噜| 在线观看av片永久免费下载| 丰满乱子伦码专区| 亚洲精品日本国产第一区| 色尼玛亚洲综合影院| 免费无遮挡裸体视频| 亚洲美女视频黄频| 欧美zozozo另类| 亚洲av在线观看美女高潮| 国产有黄有色有爽视频| 精品久久久久久久久av| 久久精品国产亚洲网站| 99久久精品一区二区三区| 麻豆久久精品国产亚洲av| 久久久精品欧美日韩精品| 日本黄色片子视频| 高清视频免费观看一区二区 | 毛片一级片免费看久久久久| 熟妇人妻不卡中文字幕| 黄色日韩在线| 亚洲一级一片aⅴ在线观看| 在线免费十八禁| 欧美高清成人免费视频www| 91aial.com中文字幕在线观看| 久久6这里有精品| 久久久a久久爽久久v久久| 国产高潮美女av| 淫秽高清视频在线观看| 天堂俺去俺来也www色官网 | 久久久久久九九精品二区国产| 91久久精品国产一区二区三区| 午夜福利视频1000在线观看| 久久精品熟女亚洲av麻豆精品 | 国产精品一及| 成人av在线播放网站| 精品久久久久久久人妻蜜臀av| 少妇人妻一区二区三区视频| 天天躁日日操中文字幕| 国产伦精品一区二区三区四那| 麻豆久久精品国产亚洲av| 久久久精品欧美日韩精品| 如何舔出高潮| 插逼视频在线观看| 麻豆国产97在线/欧美| 一级av片app| 在线天堂最新版资源| 我要看日韩黄色一级片| 精品亚洲乱码少妇综合久久| 国产精品一区二区性色av| or卡值多少钱| 特大巨黑吊av在线直播| 九九在线视频观看精品| 中国国产av一级| av线在线观看网站| 国产黄色免费在线视频| 色播亚洲综合网| 中文字幕久久专区| 国产亚洲5aaaaa淫片| 久久草成人影院| 日韩成人伦理影院| 久久精品国产自在天天线| 性色avwww在线观看| 80岁老熟妇乱子伦牲交| 哪个播放器可以免费观看大片| 一区二区三区乱码不卡18| 成人一区二区视频在线观看| 91aial.com中文字幕在线观看| 免费观看的影片在线观看| 国产伦理片在线播放av一区| 亚洲国产精品专区欧美| 免费黄色在线免费观看| 国产白丝娇喘喷水9色精品| 国产国拍精品亚洲av在线观看| 99热这里只有是精品在线观看| 国产精品一区二区在线观看99 | 六月丁香七月| 日韩一区二区视频免费看| 久久精品夜夜夜夜夜久久蜜豆| 在线 av 中文字幕| 不卡视频在线观看欧美| 亚洲18禁久久av| 免费看日本二区| 国产精品综合久久久久久久免费| 亚洲精品久久午夜乱码| 国产成人精品福利久久| 国内精品一区二区在线观看| av免费在线看不卡| 亚洲欧美中文字幕日韩二区| 岛国毛片在线播放| 久久久久久国产a免费观看| av线在线观看网站| 日产精品乱码卡一卡2卡三| 一本久久精品| av在线亚洲专区| 久久久成人免费电影| 最近2019中文字幕mv第一页| 老司机影院毛片| av一本久久久久| 男女那种视频在线观看| 日本色播在线视频| 日韩亚洲欧美综合| 免费无遮挡裸体视频| 精品一区在线观看国产| 国产综合懂色| 只有这里有精品99| 一本久久精品| 亚洲欧美日韩东京热| 一级a做视频免费观看| 久久精品国产亚洲av天美| 午夜福利网站1000一区二区三区| 嫩草影院新地址| 80岁老熟妇乱子伦牲交| 丰满人妻一区二区三区视频av| 色视频www国产| 国内精品美女久久久久久| 两个人视频免费观看高清| 秋霞伦理黄片| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的 | 成人综合一区亚洲| 又大又黄又爽视频免费| 天堂俺去俺来也www色官网 | www.色视频.com| 啦啦啦啦在线视频资源| 91久久精品电影网| 国产片特级美女逼逼视频| 精品熟女少妇av免费看| 男女视频在线观看网站免费| 成人午夜精彩视频在线观看| 我的女老师完整版在线观看| 极品教师在线视频| 91午夜精品亚洲一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产精品99久久久久久久久| 一级片'在线观看视频| 亚洲图色成人| 日韩伦理黄色片| 建设人人有责人人尽责人人享有的 | 亚洲最大成人中文| 男女国产视频网站| 国产成人a区在线观看| 看非洲黑人一级黄片| 少妇熟女欧美另类| 三级国产精品欧美在线观看| 欧美成人a在线观看| 亚洲性久久影院| 97在线视频观看| 免费少妇av软件| 精品久久久久久电影网| 国产精品久久久久久久久免| 男女边摸边吃奶| 五月天丁香电影| 插逼视频在线观看| 一个人观看的视频www高清免费观看| 欧美日本视频| 亚洲内射少妇av| 丝袜喷水一区| 精品国产三级普通话版| 男人狂女人下面高潮的视频| 神马国产精品三级电影在线观看| 五月天丁香电影| 免费观看av网站的网址| 午夜免费激情av| 欧美不卡视频在线免费观看| 亚洲av一区综合| 亚洲精华国产精华液的使用体验| 色哟哟·www| 亚洲三级黄色毛片| 亚州av有码| 伊人久久国产一区二区| 欧美bdsm另类| 久久草成人影院| 亚洲成人精品中文字幕电影| 欧美xxxx性猛交bbbb| 成人毛片a级毛片在线播放| 男插女下体视频免费在线播放| 精品不卡国产一区二区三区| 日韩欧美精品免费久久| 人人妻人人澡人人爽人人夜夜 | 亚洲精品日本国产第一区| 在线天堂最新版资源| 嫩草影院入口| 看十八女毛片水多多多| 卡戴珊不雅视频在线播放| 日韩,欧美,国产一区二区三区| 一夜夜www| 亚洲在久久综合| 国产av国产精品国产| 国产一区二区三区综合在线观看 | 亚洲四区av| 高清午夜精品一区二区三区| 欧美日韩国产mv在线观看视频 | 亚洲精品国产成人久久av| 国产又色又爽无遮挡免| 国产精品爽爽va在线观看网站| 午夜免费男女啪啪视频观看| a级毛片免费高清观看在线播放| 免费人成在线观看视频色| 天堂√8在线中文| 久久精品熟女亚洲av麻豆精品 | 性插视频无遮挡在线免费观看| 精品酒店卫生间| 色播亚洲综合网| 免费看光身美女| 国产av国产精品国产| 三级经典国产精品| 天天躁日日操中文字幕| 国产精品1区2区在线观看.| 欧美成人一区二区免费高清观看| 老司机影院毛片| 男插女下体视频免费在线播放| 国内精品美女久久久久久| 大香蕉97超碰在线| 午夜福利高清视频| 久久人人爽人人片av| 深夜a级毛片| 亚洲国产成人一精品久久久| 精品国产三级普通话版| 伊人久久国产一区二区| 久久久久久伊人网av| 亚洲在线观看片| 一级a做视频免费观看| 天堂av国产一区二区熟女人妻| 夫妻午夜视频| 中国美白少妇内射xxxbb| 乱人视频在线观看| 精品少妇黑人巨大在线播放| 日本免费a在线| 国产成人精品福利久久| 97热精品久久久久久| 国产午夜福利久久久久久| 国产精品久久久久久精品电影小说 | 高清日韩中文字幕在线| 十八禁国产超污无遮挡网站| 久久久国产一区二区| 午夜视频国产福利| 乱人视频在线观看| 日韩,欧美,国产一区二区三区| 欧美zozozo另类| 最近中文字幕2019免费版| 亚洲最大成人手机在线| 日韩av在线免费看完整版不卡| 97精品久久久久久久久久精品| 亚洲成人精品中文字幕电影| 高清毛片免费看| 亚洲精品乱久久久久久| 日韩电影二区| 欧美性感艳星| 女人久久www免费人成看片| 丰满少妇做爰视频|