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

    基于側(cè)窗濾波改進(jìn)的神經(jīng)網(wǎng)絡(luò)非均勻性校正算法

    2023-11-18 12:13:08黎明清王宇慶孫海江
    液晶與顯示 2023年11期
    關(guān)鍵詞:側(cè)窗粗糙度校正

    黎明清,王宇慶,孫海江

    (1.中國科學(xué)院 長春光學(xué)精密機(jī)械與物理研究所,吉林 長春 130033;2.中國科學(xué)院大學(xué),北京 100049)

    1 引 言

    紅外焦平面陣列探測器(IRFPA)具有靈敏度高、幀率高、結(jié)構(gòu)簡單等優(yōu)點(diǎn),已經(jīng)成為當(dāng)前紅外成像系統(tǒng)中使用的主流器件。該器件從最初的軍事偵察領(lǐng)域逐漸覆蓋醫(yī)療成像、工業(yè)檢測、消防監(jiān)控、自動駕駛等多種應(yīng)用場景[1],發(fā)揮著不可替代的作用。然而基于IRFPA 的紅外成像系統(tǒng)無法避免非均勻性噪聲問題。產(chǎn)生這一問題的主要原因包括各像元自身的響應(yīng)不均勻性以及信號讀出電路及放大器增益的非一致性[2]。目前的制造工藝水平短時間內(nèi)無法從硬件層面上完全克服非均勻性噪聲,因此在工程應(yīng)用中為了提升紅外成像系統(tǒng)圖像質(zhì)量必須進(jìn)行非均勻性校正。

    目前,非均勻性校正算法主要包括兩個不同的技術(shù)路線,分別是基于定標(biāo)實(shí)驗(yàn)的校正和場景變化的校正。定標(biāo)實(shí)驗(yàn)算法代表方案有兩點(diǎn)校正以及多點(diǎn)校正算法[3-6]。Zhang T 等人[7]提出了基于自適應(yīng)多點(diǎn)法的sCMOS 實(shí)時非均勻性校正,通過閾值比較的方式分別確定最優(yōu)分段點(diǎn)的位置以及最佳分段數(shù)量實(shí)現(xiàn)自適應(yīng)校正。該類算法需要預(yù)先做黑體輻射標(biāo)定實(shí)驗(yàn),當(dāng)環(huán)境溫度改變時器件參數(shù)會發(fā)生漂移,連續(xù)工作一段時間后需要重新標(biāo)定?;趫鼍白兓乃惴ㄖ饕〞r域高通濾波算法[8-9]、恒定統(tǒng)計算法[10-11]與神經(jīng)網(wǎng)絡(luò)算法[12-16]等。Qian 等人[8]將引導(dǎo)濾波引入時域高通濾波算法,有效改進(jìn)了“鬼影”問題與邊緣退化問題。Liu 等人[10]提出一種基于分段兩點(diǎn)校正法和局部恒定統(tǒng)計的聯(lián)合校正方法,有效增強(qiáng)了恒定統(tǒng)計法的有效性和魯棒性。Zhang L 等人[16]將雙邊濾波引入神經(jīng)網(wǎng)絡(luò)法,進(jìn)一步提升了該算法對處理后圖像的細(xì)節(jié)保護(hù)程度。此類算法基于場景圖像的變化更新校正參數(shù),具有無需標(biāo)定實(shí)驗(yàn)、校正過程不需停機(jī)的優(yōu)點(diǎn)。在基于場景變化的非均勻性校正算法中,主要考核校正后圖像殘留的非均勻性噪聲強(qiáng)度、校正后圖像有無“鬼影”現(xiàn)象以及算法的計算復(fù)雜度等方面的性能指標(biāo)。

    近年來,伴隨著視頻圖像處理設(shè)備計算能力的提升,為了解決校正參數(shù)失效導(dǎo)致的反復(fù)標(biāo)定實(shí)驗(yàn)問題,大量學(xué)者將研究方向轉(zhuǎn)移到基于場景的非均勻性校正方式,提出了基于最優(yōu)化理論計算的非均勻性校正算法。此類算法主要包括小波法、變分法、深度學(xué)習(xí)法等。2021 年,Zhang T等人[17]提出了基于小波主成分分析法,該算法從信號分解的角度尋求最優(yōu)化圖像信息主要成分達(dá)到去除非均勻性噪聲的目的,其最大的優(yōu)勢在于計算速度快,適用于尋求處理效率與處理結(jié)果平衡的應(yīng)用場景。2020 年,Yang 等人[18]利用各向異性變分法和低秩正則項相結(jié)合的方法提取條紋噪聲,實(shí)現(xiàn)了噪聲圖像分離的非均勻性校正。2022 年,Mou 等人[19]提出了基于全卷積網(wǎng)絡(luò)的紅外圖像非均勻校正算法,該網(wǎng)絡(luò)包含非均勻性等級估計子網(wǎng)絡(luò)和校正主網(wǎng)絡(luò)兩部分,有效避免了校正過程中出現(xiàn)過擬合問題。該類算法與濾波算法、統(tǒng)計法、傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法相比通常有更好的校正效果,但是此類方法也存在明顯的不足之處。小波分解主成分分析法、變分法數(shù)學(xué)計算復(fù)雜度較高、最優(yōu)參數(shù)求解復(fù)雜。深度學(xué)習(xí)法依賴專用設(shè)備且需要大量數(shù)據(jù)集標(biāo)定處理以及長時間的訓(xùn)練過程,修改網(wǎng)絡(luò)結(jié)構(gòu)后需要反復(fù)重新訓(xùn)練,不利于算法硬件實(shí)現(xiàn)。

    神經(jīng)網(wǎng)絡(luò)法作為一種經(jīng)典的基于場景變化的非均勻性校正算法,校正效果穩(wěn)定、計算復(fù)雜度低,便于工程中硬件實(shí)現(xiàn)。但是,此類算法多數(shù)利用幀間多次迭代計算校正圖像,容易出現(xiàn)圖像模糊退化以及“鬼影”現(xiàn)象,影響紅外探測器成像質(zhì)量。本文針對這一現(xiàn)象,提出一種將側(cè)窗濾波[20]與神經(jīng)網(wǎng)絡(luò)法結(jié)合的校正算法。該算法在去除非均勻性噪聲的同時較同類算法能更好地抑制“鬼影”現(xiàn)象從而獲得更好的圖像質(zhì)量。該算法基于傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)法改進(jìn),濾波算子采用側(cè)窗濾波算法,計算過程與均值濾波相似,計算復(fù)雜度遠(yuǎn)低于雙邊濾波改進(jìn)算法以及最優(yōu)化算法,為小算力、低功耗移動平臺實(shí)現(xiàn)實(shí)時非均勻性校正提供了新的研究思路。

    2 神經(jīng)網(wǎng)絡(luò)校正模型分析

    2.1 神經(jīng)網(wǎng)絡(luò)法工作原理

    圖1 為神經(jīng)網(wǎng)絡(luò)算法的結(jié)構(gòu)圖[12],它的關(guān)鍵步驟是校正層與隱含層的設(shè)計。

    輸入圖像序列首先要經(jīng)過輸入層提取單幀圖像數(shù)據(jù)送往校正層,校正層讀取增益校正參數(shù)與偏置校正參數(shù)完成圖像序列的輸出結(jié)果更新,即:

    式中:xi,j表示原始圖像位于(i,j)位置的灰度值,Gi,j表示相應(yīng)位置的增益校正參數(shù),Oi,j表示相應(yīng)位置的偏置校正參數(shù),yi,j表示相應(yīng)位置校正后的輸出結(jié)果。

    校正參數(shù)的計算需要使用隱含層。隱含層的作用是產(chǎn)生期望圖像,傳統(tǒng)算法計算過程如式(2)所示:

    式中,fi,j表示四鄰域均值濾波結(jié)果。神經(jīng)網(wǎng)絡(luò)法完成非均勻性校正的關(guān)鍵步驟是尋求使損失函數(shù)ei,j最小的增益參數(shù)和偏置參數(shù),損失函數(shù)的定義為:

    在計算損失函數(shù)最小值時需要利用梯度下降法獲得增益參數(shù)和偏置參數(shù),即:

    式中,α為迭代步長。當(dāng)前幀校正后的圖像與期望圖像之間的損失函數(shù)數(shù)值增大時,網(wǎng)絡(luò)結(jié)構(gòu)立即向減小損失函數(shù)方向傳播校正參數(shù)G、O的迭代結(jié)果,用于下一幀校正圖像計算,實(shí)現(xiàn)基于場景的自適應(yīng)校正。

    2.2 圖像模糊與偽影現(xiàn)象

    傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)法利用均值濾波獲取待校正像素點(diǎn)的期望值。均值濾波可以有效去除場景圖像中的高頻非均勻性噪聲,同樣也會去除圖像中的邊緣細(xì)節(jié)。在圖像的非邊緣區(qū)域,待校正像素鄰域內(nèi)的灰度差值較小,主要由噪聲信號引起;而在圖像的邊緣細(xì)節(jié)位置,鄰域內(nèi)的灰度差值較大,均值濾波會降低邊緣灰度梯度,引起圖像模糊退化。為了使期望圖像更加準(zhǔn)確,有學(xué)者提出了使用雙邊濾波[16]、引導(dǎo)濾波[21]以及非局部均值濾波[22]等結(jié)果作為期望圖像。上述算法在一定程度上能夠達(dá)到更好的校正效果,但是上述保邊濾波算法都是根據(jù)鄰域內(nèi)的像素灰度值相似程度賦予濾波權(quán)重,并沒有考慮真實(shí)邊緣的空間分布特征,無法從噪聲分布的結(jié)構(gòu)特征上對圖像進(jìn)行保邊濾波處理,并且隨著保邊效果的提升,計算復(fù)雜度也大幅度提升,難以在小算力、低功耗移動平臺實(shí)現(xiàn)實(shí)時校正。

    同時,對單幅圖像完成非均勻性校正后需要將期望圖像與校正圖像做差后利用梯度下降算法獲得校正參數(shù),反向傳遞至下一幀用于校正圖像。如果圖像序列相鄰幀之間元素移動過快,沿著物體運(yùn)動方向會殘留上一幀未移動圖像校正參數(shù),此時校正區(qū)域已經(jīng)沒有待校正目標(biāo),就會出現(xiàn)明顯的拖影現(xiàn)象。

    通過上文分析可知,紅外圖像非均勻性校正出現(xiàn)模糊和偽影與兩個因素有關(guān),即期望圖像獲取和幀間校正參數(shù)誤差傳遞。本文將從這兩個方向出發(fā),對傳統(tǒng)算法進(jìn)行改進(jìn)。

    3 神經(jīng)網(wǎng)絡(luò)校正模型改進(jìn)

    3.1 基于側(cè)窗濾波改進(jìn)的期望圖像獲取

    側(cè)窗濾波[20]算法由殷慧、龔元浩等在2019 年提出。該算法可以有效解決傳統(tǒng)濾波算法造成圖像邊緣模糊問題。通常為了簡化計算,在離散的情況下只定義8 個側(cè)窗,如圖2 所示。

    圖2 3×3 側(cè)窗濾波算子Fig.2 3×3 side window filter operator

    圖2 中NW 表示左上側(cè)窗,L 表示左側(cè)窗,U表示上側(cè)窗,NE 表示右上側(cè)窗,SW 表示左下側(cè)窗,D 表示下側(cè)窗,R 表示右側(cè)窗,SE 表示右下側(cè)窗。每個側(cè)窗中綠色區(qū)域表示側(cè)窗濾波有效區(qū)域,紅色區(qū)域表示側(cè)窗濾波無效區(qū)域。

    側(cè)窗濾波計算過程可以表示為:

    式中:wij表示像素i的側(cè)窗內(nèi)像素j的影響權(quán)重,q表示像素點(diǎn)的灰度值,Nn表示側(cè)窗內(nèi)權(quán)重之和,S表示側(cè)窗覆蓋的范圍,In表示第n個側(cè)窗的輸出結(jié)果。8 個側(cè)窗在離散的角度覆蓋了圖像在該點(diǎn)邊緣的所有可能幾何形狀,從中取出與原始圖像邊緣最相近的側(cè)窗輸出結(jié)果作為校正輸出結(jié)果,就可以實(shí)現(xiàn)根據(jù)邊緣分布的結(jié)構(gòu)特征對圖像進(jìn)行保邊濾波。在數(shù)值計算中,取側(cè)窗輸出與原始圖像在該點(diǎn)的灰度值之差的最小結(jié)果表示正確邊緣形狀,即:

    式中,Im表示側(cè)窗濾波輸出結(jié)果。經(jīng)典的側(cè)窗濾波技術(shù)的邊緣保持性能良好,但正是其十分出色的保邊效果導(dǎo)致殘余條紋噪聲明顯,因此在非均勻性校正問題中直接使用該算法也有著明顯的缺點(diǎn)。

    本文在側(cè)窗濾波的基礎(chǔ)上設(shè)計去噪強(qiáng)度參數(shù)θ。傳統(tǒng)側(cè)窗濾波在有效區(qū)域內(nèi)權(quán)重均衡,在無效區(qū)域內(nèi)權(quán)重為0。引入去噪?yún)?shù)θ表示無效區(qū)域像素權(quán)重與有效區(qū)域灰度權(quán)重的比值,θ越大去噪效果越明顯,但邊緣保持效果會下降。當(dāng)θ=1 時,改進(jìn)的側(cè)窗濾波算法等價于對圖像做均值濾波;當(dāng)θ=0 時,改進(jìn)側(cè)窗濾波算法等價于經(jīng)典側(cè)窗濾波算法。圖3 所示為θ=0.5 時改進(jìn)后的NW 側(cè)窗權(quán)重參數(shù)。通過改進(jìn)可以在保留側(cè)窗濾波邊緣保持優(yōu)勢情況下,將邊緣周圍圖像信息適度引入邊緣計算,達(dá)到去噪效果與保邊效果的平衡。

    圖3 改進(jìn)的3×3 側(cè)窗濾波算子Fig.3 Improved 3×3 side window filter operator

    3.2 局部發(fā)散校正參數(shù)的飽和非線性抑制

    在對校正參數(shù)反向傳遞之前,首先必須要明確校正參數(shù)是否發(fā)散,無效的校正參數(shù)傳遞會影響校正后圖像質(zhì)量。

    如式(8)所示,分析目標(biāo)所占像素點(diǎn)區(qū)域大小與濾波器寬度關(guān)系可知,使用直徑為dFLTR的濾波器能保護(hù)的最小邊緣如圖2 中最小的綠色區(qū)所示。定義該區(qū)域最小有效邊緣直徑為dValid,若紅外圖像中目標(biāo)小于最小邊緣dValid,校正參數(shù)幀間迭代過程中除了包含非均勻性噪聲外,還包含邊緣信息,導(dǎo)致該點(diǎn)參數(shù)值過大而發(fā)散,如圖4 所示。無效的局部矯正參數(shù)會對圖像產(chǎn)生過度校正,在圖像序列中的表現(xiàn)就是移動拖影現(xiàn)象。

    圖4 小目標(biāo)校正參數(shù)的發(fā)散(圖中尖峰位置)Fig.4 Divergence of small target correction parameters(peak position in the figure)

    針對圖中微小細(xì)節(jié)邊緣濾波后導(dǎo)致的校正參數(shù)發(fā)散現(xiàn)象,本文使用一種飽和非線性特性函數(shù)抑制校正參數(shù)上限,算法原理如式(9)所示:式 中:Ei,j(f)為當(dāng)前幀中位 于(i,j)位置的校正參數(shù)輸入信號;Xi,j(f)為當(dāng)前幀相應(yīng)位置的校正參數(shù)飽和特性輸出信號;Emax(f)表示當(dāng)前幀飽和非線性特性的輸出閾值,該值的確定需要統(tǒng)計當(dāng)前輸入圖像全局最大相鄰列像素點(diǎn)灰度均值的差值,非均勻性校正目的在于減小各列響應(yīng)不均勻性,因此該參數(shù)不應(yīng)當(dāng)超過原始誤差圖像中相鄰非均勻性列之間的最大灰度差值,從而實(shí)現(xiàn)視覺上的一致性。

    校正參數(shù)在Emax(f)的范圍內(nèi)是線性區(qū),取線性增益為k。當(dāng)超過線性區(qū)進(jìn)入飽和區(qū),輸入的校正參數(shù)繼續(xù)增大時,其輸出量趨于0,直到當(dāng)前位置像素點(diǎn)不再作為圖像弱小邊緣區(qū)域,校正參數(shù)會重新進(jìn)入線性區(qū),繼續(xù)經(jīng)過迭代更新適應(yīng)新的場景變化。改進(jìn)后的神經(jīng)網(wǎng)絡(luò)校正模型如圖5 所示,期望圖像采用含去噪強(qiáng)度參數(shù)的側(cè)窗濾波器獲取,校正參數(shù)的更新引入飽和非線性特性閾值函數(shù)修正。

    圖5 改進(jìn)后的神經(jīng)網(wǎng)絡(luò)模型流程框圖Fig.5 Flow chart of improved neural network correction algorithm

    4 實(shí)驗(yàn)與結(jié)果分析

    4.1 實(shí)驗(yàn)設(shè)備

    本文所有實(shí)驗(yàn)運(yùn)行計算機(jī)系統(tǒng)為Win10(CPU Intel Core i5-7300HQ 主 頻:2.30G,內(nèi)存8G),所有的算法均利用Matlab2018a 編程實(shí)現(xiàn)。

    4.2 算法對比及質(zhì)量評價

    考慮算法的硬件易實(shí)現(xiàn)性,將本文所提算法與經(jīng)典神經(jīng)網(wǎng)絡(luò)算法[12]、Zhang L 等人[16]提出的基于雙邊濾波改進(jìn)算法以及Zhang T 等人[17]提出的基于小波主成分分析法這3 種校正算法進(jìn)行比較。在實(shí)驗(yàn)中共采用3 組不同環(huán)境下的紅外相機(jī)拍攝圖像,分別為序列1:天空背景弱小飛機(jī)跟蹤目標(biāo)序列[23],共599 幀;序列2:地面復(fù)雜背景弱小飛機(jī)目標(biāo)圖像序列[23],共400 幀;序列3:道路行人檢測紅外圖像序列[24],共400 幀。

    圖6 展示的是序列1 原始圖像,圖像分辨率為256 像素×256 像素。該序列背景純凈單一,靶機(jī)目標(biāo)小,運(yùn)動速度快,圖像信噪比較低,校正過程中應(yīng)主要防止高速運(yùn)動的小目標(biāo)拖影現(xiàn)象。

    圖6 序列1 原始圖像Fig.6 Sequence 1 original images

    圖7 展示的是序列1 的非均勻性校正實(shí)驗(yàn)結(jié)果。從視覺效果上看,傳統(tǒng)的基于均值濾波的校正算法具有很強(qiáng)的平滑作用,但是在目標(biāo)邊緣產(chǎn)生大量鬼影,如圖7(a)所示,在靶機(jī)附近及右側(cè)出現(xiàn)的黑色陰影即為校正過程中引入的“鬼影”,該現(xiàn)象在目標(biāo)高速移動時尤為明顯?;陔p邊濾波的校正算法在邊緣保持以及鬼影抑制效果都有明顯提升,但如圖7(b)所示,通過肉眼仍可以觀察到有殘留的非均勻性噪聲,在靶機(jī)的右側(cè)也存在微弱的“鬼影”現(xiàn)象,并且靶機(jī)的機(jī)翼與中間部位的對比度相比原圖有所下降?;谛〔ㄖ鞒煞址治龅男U惴ㄊ轻槍螏瑘D像進(jìn)行校正,不會產(chǎn)生幀間迭代“鬼影”現(xiàn)象,但由于低信噪比圖像噪聲主成分仍包含較多噪聲信息導(dǎo)致該算法殘留非均勻性噪聲較多,如圖7(c)所示,圖中飛機(jī)局部放大圖像中存在肉眼可見非均勻性噪聲。圖7(d)為本文提出的改進(jìn)校正算法的校正結(jié)果。對比輸入圖像序列,本文算法的非均勻性噪聲去除效果好,校正后圖像沒有明顯的縱向條紋噪聲。對比文獻(xiàn)[12]、文獻(xiàn)[16]所用校正算法,本文算法的鬼影現(xiàn)象有明顯改善,細(xì)節(jié)邊緣保持程度也達(dá)到雙邊濾波效果,特別是在靶機(jī)機(jī)翼位置與中間位置的對比度更接近原始圖像序列,明顯優(yōu)于文獻(xiàn)[16]所用算法。對比文獻(xiàn)[17]所用算法,本文算法的殘余非均勻性噪聲明顯降低,圖像細(xì)節(jié)保護(hù)程度無明顯差別。

    圖7 不同算法對序列1 圖像的校正效果對比Fig.7 Comparison of correction effects of different algorithms on sequence 1

    為了進(jìn)一步評估3 種校正算法,本文引入列均值[25]與圖像粗糙度指標(biāo)ρ[26]進(jìn)行定量評估。圖像粗糙度ρ的定義為:

    式中:h1=[-1,1]為圖像在水平方向的卷積核模板,h2=[-1,1]T為圖像豎直方向的卷積核模板,I表示輸入圖像,‖ ‖表示取1 范數(shù)運(yùn)算,*表示做卷積操作。列均值指標(biāo)反映圖像在各通道增益及偏置的不均勻性,該指標(biāo)的變化曲線越平坦說明圖像均勻性越好。圖像粗糙度反映圖像整體灰度值的均勻性程度,該指標(biāo)越低說明圖像均勻性越好。

    從圖8(a)中可以看出,原始圖像在90~120 列存在明顯的平均值波動,表現(xiàn)在圖像上就是明暗相間的條紋噪聲。經(jīng)過文獻(xiàn)[12]算法處理后,各列均值波動幅度明顯下降。文獻(xiàn)[16]算法因?yàn)槭褂秒p邊濾波獲取期望圖像,處理后各列均值波動幅度較傳統(tǒng)算法略大,整體波動幅度較原始圖像序列有明顯改善。文獻(xiàn)[17]算法使用小波分解后的主成分分析分離非均勻性噪聲,在針對低信噪比圖像的處理過程中殘留有較大的非均勻性,導(dǎo)致各列均值波動幅度降低程度有限。使用本文提出的改進(jìn)算法校正后,圖像的各列均值過渡平坦,沒有明顯波動情況,處于另外兩種算法幅度變化的包絡(luò)曲線之中,說明相比經(jīng)典神經(jīng)網(wǎng)絡(luò)校正算法、基于雙邊濾波改進(jìn)的校正算法以及基于小波主成分分析算法所獲得的列均值更穩(wěn)定。從圖8(b)中可以觀察到,經(jīng)過文獻(xiàn)[17]算法處理后,圖像的整體粗糙度最大,但是該算法由于不涉及幀間迭代過程,所以從第一幀開始算法效果就趨于穩(wěn)定,而其他3 種算法在經(jīng)過約80 幀圖像數(shù)據(jù)迭代過程后才趨于穩(wěn)定。經(jīng)過文獻(xiàn)[16]算法處理后,圖像的整體粗糙度明顯降低。經(jīng)過本文算法校正后,粗糙度指標(biāo)明顯低于基于雙邊濾波改進(jìn)的校正算法,與文獻(xiàn)[12]所提經(jīng)典校正算法的粗糙度指標(biāo)相似,說明相比經(jīng)典神經(jīng)網(wǎng)絡(luò)校正算法、基于雙邊濾波改進(jìn)的校正算法以及基于小波主成分分析算法,本文方法的圖像整體均勻性最好。

    圖8 (a)序列1 列均值曲線;(b)序列1 粗糙度曲線。Fig.8 (a)Sequence 1 mean curves;(b)Sequence 1 roughness curves.

    由上述分析可知,本文提出的基于側(cè)窗濾波改進(jìn)校正算法在迭代效果趨于穩(wěn)定后的圖像視覺效果、列均值指標(biāo)、粗糙度指標(biāo)均優(yōu)于另外3 種對比算法,綜合考慮,本文算法針對圖像序列1 的校正效果是4 種方法中最好的。

    圖9 展示的是序列2 原始圖像,圖像分辨率為256 像素×256 像素。該環(huán)境下背景復(fù)雜、圖像對比度低、靶機(jī)目標(biāo)小且運(yùn)動快,圖像信噪比較低,校正過程中應(yīng)主要避免背景輪廓過度退化以及“鬼影”現(xiàn)象。

    圖9 序列2 原始圖像Fig.9 Sequence 2 original images

    圖10展示的是序列2非均勻性校正實(shí)驗(yàn)結(jié)果。從圖10(a)中可以看出,使用經(jīng)典校正算法處理后的圖像序列中非均勻性噪聲有明顯改善,但是沿著飛機(jī)移動方向存在明顯拖影現(xiàn)象,并且由于鏡頭移動,復(fù)雜場景中房屋位置周圍也產(chǎn)生了明顯的“鬼影”。使用基于雙邊濾波改進(jìn)的校正算法無明顯“鬼影”現(xiàn)象,但是圖像整體存在微弱的殘余噪聲,如圖10(b)所示?;谛〔ㄖ鞒煞址治龅男U惴▽儆趩螏?,不會產(chǎn)生幀間迭代“鬼影”現(xiàn)象,但校正后圖像仍殘留明顯非均勻性噪聲,如圖10(c)所示,圖中房屋處局部放大圖像中存在肉眼可見的非均勻性噪聲。圖10(d)為使用本文提出的校正算法的校正結(jié)果,對比輸入圖像,本文算法的非均勻性噪聲去除效果較好;對比文獻(xiàn)[12]算法,“鬼影”現(xiàn)象有明顯的改善;對比文獻(xiàn)[16]、文獻(xiàn)[17]算法,殘余非均勻性噪聲有明顯改善,在邊緣保持方面達(dá)到雙邊濾波改進(jìn)的水平。

    圖10 不同算法對序列2 圖像的校正效果對比Fig.10 Comparison of correction effects of different algorithms on sequence 2

    從圖11(a)中可以觀察到,原始圖像在90~120 列存在明顯的平均值波動。經(jīng)過文獻(xiàn)[12]算法處理后,各列均值波動幅度下降明顯。經(jīng)過文獻(xiàn)[16]算法處理后,均值波動幅度略大于傳統(tǒng)校正算法。經(jīng)過文獻(xiàn)[17]算法處理后,均值波動幅度較原始圖像數(shù)據(jù)有所降低,但仍存在較大的列間差異性。使用本文改進(jìn)算法校正后,圖像在各列均值性能均優(yōu)于另外3 種校正算法,列均值指標(biāo)均位于經(jīng)典神經(jīng)網(wǎng)絡(luò)校正算法、基于雙邊濾波改進(jìn)的校正算法以及基于小波主成分分析校正算法的相應(yīng)指標(biāo)包絡(luò)線內(nèi),校正結(jié)果穩(wěn)定。從圖11(b)中可以看出,復(fù)雜背景下場景快速運(yùn)動序列的圖像粗糙度在時間軸上呈現(xiàn)波動趨勢。場景的劇烈變化引起校正參數(shù)波動,導(dǎo)致傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法、基于雙邊濾波改進(jìn)算法和本文算法3 種校正算法校正后的粗糙度相近,基于小波主成分分析的校正算法同樣因?yàn)閳D片信噪比低而呈現(xiàn)較高粗糙度。分析曲線可知,在整體上本文算法校正后的圖像粗糙度低于文獻(xiàn)[16]、文獻(xiàn)[17]所用算法,略大于文獻(xiàn)[12]所用算法。

    圖11 (a)序列2 列均值曲線;(b)序列2 粗糙度曲線。Fig.11 (a)Sequence 2 mean curves;(b)Sequence 2 roughness curves.

    由上述分析可知,本文提出的基于側(cè)窗濾波改進(jìn)的校正算法在視覺效果上優(yōu)于文獻(xiàn)[12]、文獻(xiàn)[17]算法,與文獻(xiàn)[16]算法無明顯區(qū)別。列均值指標(biāo)優(yōu)于其余3 種對比算法。圖像粗糙度指標(biāo)優(yōu)于文獻(xiàn)[16]、文獻(xiàn)[17]算法,與文獻(xiàn)[12]算法無明顯區(qū)別。綜合考慮,本文算法針對圖像序列2 的校正效果是4 種方法中最好的。

    圖12 展示的是序列3 的原始圖像,圖像分辨率為1 024 像素×512 像素。該環(huán)境下,目標(biāo)邊緣細(xì)節(jié)復(fù)雜,圖像信噪比較高,場景運(yùn)動緩慢,校正過程中應(yīng)主要避免目標(biāo)復(fù)雜輪廓過度退化。

    圖12 序列3 原始圖像Fig.12 Sequence 3 original images

    圖13展示的是序列3非均勻性校正實(shí)驗(yàn)結(jié)果。從圖13(a)中可以看出,使用經(jīng)典校正算法處理后的圖像序列的非均勻性噪聲去除效果明顯,如放大圖中路面停止線處所示。圖像在整體上有明顯的模糊退化現(xiàn)象,局部存在細(xì)節(jié)丟失,對比度下降問題,如圖中放大后的市政護(hù)欄已經(jīng)無法分辨出內(nèi)部欄桿的邊緣。從圖13(b)中可以看出,使用基于雙邊濾波改進(jìn)的校正算法后,圖像清晰度較傳統(tǒng)算法有明顯改善,圖中放大后的市政護(hù)欄的內(nèi)部欄桿邊緣清晰可見,與原始圖像序列無明顯差異,但該算法也存在肉眼可見的殘余非均勻性噪聲問題,如放大圖中路面停止線處所示。從圖13(c)中可以看出,使用基于小波主成分分析算法校正后,圖像非均勻噪聲去除效果明顯,圖像邊緣細(xì)節(jié)保持性強(qiáng),局部放大圖中欄桿細(xì)節(jié)保存完好,可見該算法在處理高信噪比圖像時具有良好性能。但是由于小波分解后,最優(yōu)主成分替換圖像原始信息,在原圖中部分無噪聲區(qū)域形成少量明暗斑紋。圖13(d)為本文提出的改進(jìn)校正算法的校正結(jié)果。對比輸入圖像序列,本文算法的非均勻性噪聲去除效果較好,校正后圖像沒有明顯縱向條紋噪聲。對比文獻(xiàn)[12]所用算法,圖像清晰度有明顯改善。對比文獻(xiàn)[16]所用算法,有效減少了殘余非均勻性噪聲現(xiàn)象,在邊緣保持與細(xì)節(jié)保護(hù)方面達(dá)到基于雙邊濾波改進(jìn)的算法校正效果。對比文獻(xiàn)[17]所用算法,在邊緣保持性能上略有降低,圖像整體一致性高,沒有引入新的明暗斑紋。

    圖13 不同算法對序列3 圖像校正效果對比Fig.13 Comparison of correction effects of different algorithms on sequence 3

    從圖14(a)中可以看出,原始圖像在595~605列存在明顯的平均值波動。經(jīng)過文獻(xiàn)[12]算法處理后,各列均值波動幅度明顯下降。經(jīng)過文獻(xiàn)[16]算法處理后,均值波動幅度略大于傳統(tǒng)校正算法。使用文獻(xiàn)[17]算法處理后,均值波動幅度整體較原始序列有所下降,且沒有迭代過程,經(jīng)過第一張輸入圖像校正后即達(dá)到穩(wěn)定狀態(tài),但整體仍存在較大波動幅度。使用本文改進(jìn)算法校正后,圖像在各列均值性能優(yōu)于另外3 種校正算法,列均值指標(biāo)均位于經(jīng)典校正算法、基于雙邊濾波校正算法和基于小波主成分分析算法相應(yīng)指標(biāo)包絡(luò)線內(nèi),校正結(jié)果穩(wěn)定。從圖14(b)中可以觀察到,經(jīng)過文獻(xiàn)[17]算法處理后,圖像的整體粗糙度最大;經(jīng)過文獻(xiàn)[16]算法處理后,圖像的整體粗糙度略高于本文算法;經(jīng)過本文算法校正后,粗糙度指標(biāo)略高于文獻(xiàn)[12]所提經(jīng)典校正算法過度平滑后的粗糙度指標(biāo)。

    圖14 (a)序列3 列均值曲線;(b)序列3 粗糙度曲線。Fig.14 (a)Sequence 3 mean curves;(b)Sequence 3 roughness curves.

    由上述分析可知,本文提出的基于側(cè)窗濾波改進(jìn)的校正算法在列均值指標(biāo)性能優(yōu)于另外3 種對比算法;視覺效果明顯優(yōu)于傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法與基于雙邊濾波改進(jìn)算法;邊緣保持程度略低與文獻(xiàn)[17]小波主成分分析算法;圖像粗糙度指標(biāo)優(yōu)于文獻(xiàn)[16]、文獻(xiàn)[17]算法,低于文獻(xiàn)[12]算法過度平滑后的結(jié)果。綜合考慮,本文算法針對圖像序列3 是4 種方法中效果最好的。

    表1 統(tǒng)計了3 組圖像序列中第80~400 幀圖像粗糙度ρ的平均值,該指標(biāo)越低表明圖像一致性越好。從表1 可以看出,基于小波主成分分析算法使用單幀圖像校正,校正后的低信噪比圖像粗糙度指標(biāo)下降幅度小,高信噪比圖像粗糙度指標(biāo)下降明顯,但該算法在粗糙度指標(biāo)整體表現(xiàn)不足,各項數(shù)據(jù)均高于對比算法。本文提出的基于側(cè)窗濾波改進(jìn)校正算法在3 組測試序列中均表現(xiàn)出良好性能,其中在測試序列1、測試序列2 中粗糙度指標(biāo)為4 種校正算法中最低,在測試序列3中該指標(biāo)明顯低于雙邊濾波改進(jìn)的神經(jīng)網(wǎng)絡(luò)算法與基于小波主成分分析算法,高于傳統(tǒng)神經(jīng)網(wǎng)絡(luò)校正算法,但是傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法存在嚴(yán)重的模糊退化,過度的平滑性能導(dǎo)致該算法產(chǎn)生大量的“鬼影”,不能適用于真實(shí)檢測環(huán)境。

    表1 不同算法處理后的圖像序列平均粗糙度指標(biāo)對比Tab.1 Comparison of average roughness indexes of image sequences processed by different algorithms

    表2 統(tǒng)計了3 組圖像序列中處理400 幀圖像的耗時。該指標(biāo)涉及兩個影響參數(shù),即算法計算復(fù)雜度和輸入圖形序列像素尺寸。分析表2 中的數(shù)據(jù),橫向?qū)Ρ瓤芍獙τ谙嗤斎雸D像序列,本文所提算法處理相同幀圖像耗時略長于傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法,明顯優(yōu)于基于雙邊濾波改進(jìn)的非均勻性校正算法,與基于小波主成分分析法相比也有一定優(yōu)勢??v向?qū)Ρ? 種算法處理速度相對于不同尺寸輸入圖形序列可知,傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法、基于小波主成分分析算法和本文所提算法的耗時與圖形尺寸接近線性關(guān)系,而基于雙邊濾波改進(jìn)的神經(jīng)網(wǎng)絡(luò)算法在處理高分辨率圖像時耗時增加嚴(yán)重,原因是雙邊濾波器計算復(fù)雜度是輸入圖像空域尺寸平方項函數(shù)關(guān)系,該特性限制了算法的硬件實(shí)現(xiàn)。綜上可知,在處理速度方面,傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法優(yōu)于本文算法;基于小波主成分分析的校正算法處理時間與本文提出算法較為接近,表明該算法優(yōu)勢在于平衡處理速度和處理后圖像效果;基于雙邊濾波算法時耗性能最差。

    為滿足紅外成像器件非均勻性校正實(shí)時性要求,需要處理速度大于視頻圖像輸入速度。通常監(jiān)控錄像機(jī)系統(tǒng)視頻輸入接口幀率為25 幀/s,單幀耗時不超過40 ms。針對256×256 分辨率圖像序列,傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法處理一幅圖像耗時約為7.5 ms,雙邊濾波改進(jìn)算法耗時約為275 ms,小波主成分分析算法耗時約為35 ms,本文所提算法耗時約為12.5 ms;針對1 024×512 分辨率圖像序列,傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法處理一幅圖像耗時約為25 ms,雙邊濾波改進(jìn)算法耗時約為1 900 ms,小波主成分分析算法耗時約為175 ms,本文所提算法耗時約為95 ms??梢娫诨贑PU 處理器運(yùn)算平臺上,傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法、小波主成分分析算法、本文所提算法滿足小分辨率圖像的實(shí)時校正,大分辨率圖像僅有傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法滿足要求。為進(jìn)一步提升算法的處理效率,實(shí)現(xiàn)小算力、低功耗移動平臺的實(shí)時非均勻性校正,可采用現(xiàn)場可編程門陣列(FPGA)流水線處理實(shí)現(xiàn)并行加速,而小波主成分分解算法中涉及最優(yōu)化主成分計算過程,F(xiàn)PGA 硬件實(shí)現(xiàn)較為復(fù)雜,不利于實(shí)現(xiàn)實(shí)時非均勻性校正。綜合表1 和表2 可以看出,本文算法在非均勻性校正效果和算法運(yùn)行效率方面都具有明顯優(yōu)勢。

    5 結(jié) 論

    本文提出了一種基于側(cè)窗濾波技術(shù)改進(jìn)的神經(jīng)網(wǎng)絡(luò)非均勻性校正算法。該算法通過改進(jìn)的側(cè)窗濾波技術(shù)獲得期望校正圖像,通過幀間飽和非線性函數(shù)抑制校正參數(shù)局部發(fā)散,優(yōu)化了傳統(tǒng)神經(jīng)網(wǎng)絡(luò)算法的保邊效果,能夠有效避免校正后圖像出現(xiàn)“鬼影”問題,降低了基于雙邊濾波改進(jìn)算法的計算復(fù)雜度,較基于小波主成分分析算法更適用于不同信噪比的圖像序列且處理速度更快。實(shí)驗(yàn)結(jié)果表明,本文提出的算法在測試圖像序列1 上的平均粗糙度為0.030 3,連續(xù)處理400 幀圖像耗時4.938 5 s;在測試圖像序列2 上的平均粗糙度為0.058 6,連續(xù)處理400 幀圖像耗時4.916 2 s;在測試圖像序列3 上的平均粗糙度為0.079 2,連續(xù)處 理400 幀圖像耗時37.417 0 s。該算法在3 組測試圖像序列中平均圖像粗糙度降低了30.17%,平衡處理時間與處理后圖像質(zhì)量性能優(yōu)于對比算法,在運(yùn)算時間上有顯著提升。目前所提算法面臨的最大問題是如何進(jìn)一步提高運(yùn)算速度,實(shí)現(xiàn)高分辨率紅外圖像實(shí)時非均勻性校正。而在該算法中使用的關(guān)鍵步驟側(cè)窗濾波本質(zhì)為多次加權(quán)均值濾波,便于FPGA 并行加速處理,因此本文也為小算力、低功耗移動平臺實(shí)現(xiàn)實(shí)時非均勻性校正提供了新的研究思路。

    猜你喜歡
    側(cè)窗粗糙度校正
    劉光第《南旋記》校正
    國學(xué)(2020年1期)2020-06-29 15:15:30
    基于無人機(jī)影像的巖體結(jié)構(gòu)面粗糙度獲取
    甘肅科技(2020年20期)2020-04-13 00:30:18
    DJ4電力機(jī)車司機(jī)室活動側(cè)窗鎖閉器的故障原因分析及處理
    冷沖模磨削表面粗糙度的加工試驗(yàn)與應(yīng)用
    模具制造(2019年4期)2019-06-24 03:36:48
    一類具有校正隔離率隨機(jī)SIQS模型的絕滅性與分布
    機(jī)內(nèi)校正
    城軌車輛側(cè)窗粘接工藝分析
    基于BP神經(jīng)網(wǎng)絡(luò)的面齒輪齒面粗糙度研究
    鋼材銹蝕率與表面三維粗糙度參數(shù)的關(guān)系
    轎車乘員艙內(nèi)氣流對氣動阻力影響的數(shù)值模擬
    日韩av在线免费看完整版不卡| 少妇裸体淫交视频免费看高清| 丝袜脚勾引网站| 亚洲美女搞黄在线观看| 黄色一级大片看看| 99热国产这里只有精品6| 91久久精品国产一区二区三区| 亚洲av免费高清在线观看| 久久久久久人妻| 久久久欧美国产精品| 国产精品一区二区性色av| 久久久久久久久大av| 一二三四中文在线观看免费高清| 热re99久久精品国产66热6| 成人午夜精彩视频在线观看| 人妻夜夜爽99麻豆av| 国产精品秋霞免费鲁丝片| 精品少妇久久久久久888优播| 麻豆精品久久久久久蜜桃| 日韩欧美一区视频在线观看 | 大香蕉97超碰在线| 日日啪夜夜爽| 免费观看av网站的网址| 97超视频在线观看视频| 亚洲高清免费不卡视频| 精品人妻一区二区三区麻豆| av.在线天堂| 99久久精品热视频| 能在线免费看毛片的网站| 黑人巨大精品欧美一区二区蜜桃 | 另类亚洲欧美激情| 美女大奶头黄色视频| 爱豆传媒免费全集在线观看| 男人爽女人下面视频在线观看| 国产亚洲av片在线观看秒播厂| 日韩熟女老妇一区二区性免费视频| 亚洲美女黄色视频免费看| 亚洲欧洲精品一区二区精品久久久 | 人人澡人人妻人| 亚洲精品第二区| av视频免费观看在线观看| 人妻夜夜爽99麻豆av| 自拍欧美九色日韩亚洲蝌蚪91 | 日日啪夜夜撸| 久久久久人妻精品一区果冻| 亚洲精品中文字幕在线视频 | 亚洲内射少妇av| 国产精品熟女久久久久浪| 观看免费一级毛片| 国产女主播在线喷水免费视频网站| 一本色道久久久久久精品综合| 成人无遮挡网站| 99久久人妻综合| 久久久国产精品麻豆| 日韩伦理黄色片| 久久6这里有精品| 国产精品蜜桃在线观看| 三级经典国产精品| 色婷婷av一区二区三区视频| 观看美女的网站| 国产精品久久久久成人av| 97精品久久久久久久久久精品| 一级av片app| 亚洲精品456在线播放app| 欧美激情国产日韩精品一区| av福利片在线| 卡戴珊不雅视频在线播放| 国产精品伦人一区二区| 亚洲欧美日韩东京热| 欧美日韩视频精品一区| 日日爽夜夜爽网站| 亚洲精品国产av蜜桃| 国产在线免费精品| 亚洲av男天堂| 精品一区二区三区视频在线| 久久av网站| 久久99蜜桃精品久久| 99九九线精品视频在线观看视频| 18禁动态无遮挡网站| 黄色配什么色好看| av免费在线看不卡| 一本大道久久a久久精品| 老女人水多毛片| 日本免费在线观看一区| 亚洲欧洲日产国产| 免费看av在线观看网站| 久久久久久久久久久免费av| 日韩三级伦理在线观看| 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 日日撸夜夜添| 国产精品一区www在线观看| 人妻 亚洲 视频| 久久久亚洲精品成人影院| a级毛片免费高清观看在线播放| 秋霞在线观看毛片| 国产伦在线观看视频一区| 如日韩欧美国产精品一区二区三区 | 日本vs欧美在线观看视频 | 男人舔奶头视频| 成人影院久久| 熟妇人妻不卡中文字幕| 熟女电影av网| a级毛片免费高清观看在线播放| 国产精品一二三区在线看| 精品人妻熟女毛片av久久网站| 3wmmmm亚洲av在线观看| 夜夜骑夜夜射夜夜干| 精品久久久噜噜| 久久99一区二区三区| 国产精品三级大全| 中文字幕制服av| h日本视频在线播放| 中国三级夫妇交换| 少妇人妻精品综合一区二区| 国产成人精品福利久久| 最近最新中文字幕免费大全7| 中文字幕精品免费在线观看视频 | 国产日韩欧美亚洲二区| a级毛片在线看网站| 亚洲av电影在线观看一区二区三区| 免费看日本二区| 五月开心婷婷网| 精品一区二区免费观看| 亚洲国产精品成人久久小说| 欧美人与善性xxx| 人人妻人人澡人人看| 一级毛片我不卡| 欧美97在线视频| 亚洲精品自拍成人| 一级毛片久久久久久久久女| 丁香六月天网| 亚洲精品久久久久久婷婷小说| 久久久久久久精品精品| 免费少妇av软件| 最新的欧美精品一区二区| 精品少妇久久久久久888优播| 丝瓜视频免费看黄片| 高清午夜精品一区二区三区| 久久久精品94久久精品| 国产在线一区二区三区精| 亚洲国产精品一区二区三区在线| 激情五月婷婷亚洲| 黑人巨大精品欧美一区二区蜜桃 | 亚洲成人av在线免费| 精品午夜福利在线看| 亚洲欧美中文字幕日韩二区| 亚洲中文av在线| 免费不卡的大黄色大毛片视频在线观看| 精品一区二区三区视频在线| 高清视频免费观看一区二区| 国产精品国产三级国产av玫瑰| 日韩 亚洲 欧美在线| 国产伦在线观看视频一区| 国产伦精品一区二区三区视频9| 美女cb高潮喷水在线观看| 亚洲人成网站在线观看播放| 少妇人妻久久综合中文| 99热6这里只有精品| 国产在线免费精品| 激情五月婷婷亚洲| 亚洲不卡免费看| 亚洲成人av在线免费| 色哟哟·www| 亚洲成人一二三区av| 99精国产麻豆久久婷婷| 亚洲国产精品一区三区| 免费黄频网站在线观看国产| 国产真实伦视频高清在线观看| 精品少妇内射三级| 久久久欧美国产精品| 国产极品粉嫩免费观看在线 | 久久99热6这里只有精品| www.av在线官网国产| 亚洲一级一片aⅴ在线观看| 亚洲人与动物交配视频| 黑丝袜美女国产一区| 2018国产大陆天天弄谢| 中文字幕亚洲精品专区| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 亚洲国产毛片av蜜桃av| 国产一区二区三区av在线| videossex国产| 国产在线男女| 国产精品秋霞免费鲁丝片| 国产精品99久久99久久久不卡 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品日本国产第一区| 亚洲久久久国产精品| 一本—道久久a久久精品蜜桃钙片| 最近2019中文字幕mv第一页| 亚洲av二区三区四区| √禁漫天堂资源中文www| 丝袜在线中文字幕| 日韩人妻高清精品专区| 丰满少妇做爰视频| 亚洲久久久国产精品| 五月伊人婷婷丁香| 欧美日韩视频高清一区二区三区二| 大香蕉97超碰在线| 亚洲国产精品国产精品| 大香蕉久久网| 亚洲欧美精品专区久久| 免费av中文字幕在线| 秋霞伦理黄片| 天天操日日干夜夜撸| 在线观看三级黄色| 大香蕉97超碰在线| 欧美97在线视频| 交换朋友夫妻互换小说| 高清av免费在线| 丝袜脚勾引网站| 久久亚洲国产成人精品v| 天天躁夜夜躁狠狠久久av| 国产精品三级大全| 在线亚洲精品国产二区图片欧美 | 免费久久久久久久精品成人欧美视频 | 精品一区二区免费观看| 国产一级毛片在线| 久久久午夜欧美精品| 亚洲高清免费不卡视频| 十八禁网站网址无遮挡 | 国产精品不卡视频一区二区| 大码成人一级视频| 午夜精品国产一区二区电影| 欧美亚洲 丝袜 人妻 在线| 王馨瑶露胸无遮挡在线观看| 日韩在线高清观看一区二区三区| 国产成人免费观看mmmm| 久久99热6这里只有精品| 欧美日韩一区二区视频在线观看视频在线| 天美传媒精品一区二区| 18+在线观看网站| 亚洲综合色惰| 成年av动漫网址| 日韩av在线免费看完整版不卡| av国产久精品久网站免费入址| 午夜免费男女啪啪视频观看| 欧美精品国产亚洲| 伊人久久国产一区二区| 婷婷色综合大香蕉| 99热全是精品| 国产女主播在线喷水免费视频网站| 卡戴珊不雅视频在线播放| 人人妻人人看人人澡| 亚洲一级一片aⅴ在线观看| 男人狂女人下面高潮的视频| 亚洲精品亚洲一区二区| 色婷婷久久久亚洲欧美| 国产有黄有色有爽视频| 男男h啪啪无遮挡| 国产av一区二区精品久久| 亚洲av成人精品一二三区| 麻豆精品久久久久久蜜桃| 永久网站在线| 亚洲国产最新在线播放| h日本视频在线播放| 一级二级三级毛片免费看| 国产高清三级在线| 一区二区av电影网| 亚洲精品国产av蜜桃| 最新的欧美精品一区二区| 高清欧美精品videossex| 国产精品99久久久久久久久| 久久影院123| 免费少妇av软件| 人人妻人人澡人人看| 国产视频内射| 男女啪啪激烈高潮av片| 午夜老司机福利剧场| 熟女av电影| 在线观看三级黄色| 丝袜喷水一区| 亚洲av成人精品一二三区| 性色av一级| 成年美女黄网站色视频大全免费 | 亚洲av在线观看美女高潮| 国产黄色免费在线视频| 日韩 亚洲 欧美在线| 欧美三级亚洲精品| 午夜91福利影院| 欧美人与善性xxx| 夫妻午夜视频| 久久久久久久亚洲中文字幕| 最近最新中文字幕免费大全7| 在线亚洲精品国产二区图片欧美 | 综合色丁香网| 亚洲性久久影院| 国国产精品蜜臀av免费| 少妇精品久久久久久久| 日韩不卡一区二区三区视频在线| 中文字幕亚洲精品专区| 精华霜和精华液先用哪个| 久久精品国产亚洲网站| 最新的欧美精品一区二区| 男人爽女人下面视频在线观看| 大陆偷拍与自拍| 亚洲精品久久午夜乱码| 在线精品无人区一区二区三| 久久午夜综合久久蜜桃| 亚洲精品一二三| 一本色道久久久久久精品综合| 国产91av在线免费观看| 69精品国产乱码久久久| 如何舔出高潮| 中文欧美无线码| 亚洲欧美日韩另类电影网站| 国产黄色免费在线视频| 亚洲精品成人av观看孕妇| 纵有疾风起免费观看全集完整版| 日本欧美国产在线视频| 一级毛片黄色毛片免费观看视频| 一区二区三区免费毛片| 女性生殖器流出的白浆| 91精品伊人久久大香线蕉| 少妇的逼水好多| 大码成人一级视频| 黄色欧美视频在线观看| 啦啦啦在线观看免费高清www| videos熟女内射| 又大又黄又爽视频免费| 亚洲av不卡在线观看| 最黄视频免费看| 久久韩国三级中文字幕| 日韩熟女老妇一区二区性免费视频| 天美传媒精品一区二区| 国产在线视频一区二区| 亚洲精品久久午夜乱码| 色5月婷婷丁香| 欧美另类一区| 欧美变态另类bdsm刘玥| 国国产精品蜜臀av免费| 日韩欧美一区视频在线观看 | a级片在线免费高清观看视频| 中文精品一卡2卡3卡4更新| av线在线观看网站| 九九久久精品国产亚洲av麻豆| 午夜精品国产一区二区电影| 日韩不卡一区二区三区视频在线| 亚洲av在线观看美女高潮| 又大又黄又爽视频免费| 老熟女久久久| 精品一区二区免费观看| 久久久久久久久久久丰满| 3wmmmm亚洲av在线观看| 熟妇人妻不卡中文字幕| 99久久精品热视频| 国产在线一区二区三区精| 日韩不卡一区二区三区视频在线| 乱人伦中国视频| 日韩伦理黄色片| 欧美 日韩 精品 国产| 国产 一区精品| 国产视频内射| 亚洲熟女精品中文字幕| 欧美 日韩 精品 国产| 一区在线观看完整版| 免费在线观看成人毛片| 国产极品天堂在线| 色视频www国产| 免费人妻精品一区二区三区视频| 韩国高清视频一区二区三区| 黄片无遮挡物在线观看| 美女视频免费永久观看网站| av播播在线观看一区| 女性被躁到高潮视频| 国产成人精品久久久久久| 国产精品久久久久成人av| 欧美三级亚洲精品| 曰老女人黄片| 22中文网久久字幕| 青春草国产在线视频| 蜜臀久久99精品久久宅男| 国产乱来视频区| 亚洲av中文av极速乱| 亚洲精品中文字幕在线视频 | 国产精品一区www在线观看| av网站免费在线观看视频| 欧美日韩av久久| 天天操日日干夜夜撸| 亚洲精品日本国产第一区| 国产精品一区二区性色av| 亚洲四区av| 一级毛片黄色毛片免费观看视频| 中文字幕精品免费在线观看视频 | 女人久久www免费人成看片| 久久97久久精品| 五月开心婷婷网| 亚洲婷婷狠狠爱综合网| 精品熟女少妇av免费看| 乱系列少妇在线播放| 亚洲精品亚洲一区二区| 欧美3d第一页| 少妇人妻 视频| av黄色大香蕉| 在线观看三级黄色| 亚洲人成网站在线观看播放| 天堂8中文在线网| 一个人看视频在线观看www免费| 久久99热这里只频精品6学生| 99re6热这里在线精品视频| 久久97久久精品| 国产白丝娇喘喷水9色精品| 一区二区三区免费毛片| 寂寞人妻少妇视频99o| 亚洲国产精品成人久久小说| 汤姆久久久久久久影院中文字幕| 久热久热在线精品观看| 色婷婷av一区二区三区视频| 欧美97在线视频| 久久久精品94久久精品| 一个人看视频在线观看www免费| 一区二区三区四区激情视频| 人人妻人人澡人人爽人人夜夜| 亚洲高清免费不卡视频| 亚洲国产精品一区二区三区在线| 99热6这里只有精品| av播播在线观看一区| 18禁动态无遮挡网站| 精品国产一区二区久久| 久久久a久久爽久久v久久| 国产国拍精品亚洲av在线观看| 赤兔流量卡办理| 久久人人爽av亚洲精品天堂| 久久久国产精品麻豆| 午夜福利影视在线免费观看| 22中文网久久字幕| 成人二区视频| 亚洲综合色惰| av专区在线播放| 久久久久人妻精品一区果冻| 免费大片黄手机在线观看| 嫩草影院入口| 永久免费av网站大全| 秋霞在线观看毛片| 国产极品粉嫩免费观看在线 | 国产精品不卡视频一区二区| 一二三四中文在线观看免费高清| 肉色欧美久久久久久久蜜桃| 欧美xxxx性猛交bbbb| 一区二区三区精品91| 中文字幕制服av| av免费在线看不卡| 两个人免费观看高清视频 | 插逼视频在线观看| 亚洲国产精品成人久久小说| 国产精品无大码| 亚洲图色成人| 国产精品免费大片| 99国产精品免费福利视频| 99久久综合免费| 亚洲国产精品一区二区三区在线| 国产精品国产三级国产专区5o| 国产高清有码在线观看视频| 少妇高潮的动态图| av不卡在线播放| 亚洲精品视频女| 欧美 日韩 精品 国产| 国产成人91sexporn| 国产淫语在线视频| 偷拍熟女少妇极品色| 成人亚洲精品一区在线观看| 久久久久网色| 国产一区二区三区av在线| 亚洲av在线观看美女高潮| 啦啦啦在线观看免费高清www| 日韩在线高清观看一区二区三区| 91精品国产国语对白视频| 久久人妻熟女aⅴ| 国产成人aa在线观看| 一区在线观看完整版| 极品少妇高潮喷水抽搐| √禁漫天堂资源中文www| 亚洲在久久综合| 性色avwww在线观看| 国产精品一区二区性色av| 久久99热6这里只有精品| 亚洲av成人精品一二三区| 日韩视频在线欧美| 欧美日韩视频高清一区二区三区二| 国产精品一区二区在线观看99| 中文乱码字字幕精品一区二区三区| av在线观看视频网站免费| 欧美成人午夜免费资源| 精品久久久久久久久av| 久久久久国产网址| 五月天丁香电影| 亚洲欧美精品专区久久| 中文字幕人妻熟人妻熟丝袜美| √禁漫天堂资源中文www| 精品少妇黑人巨大在线播放| 亚洲国产最新在线播放| 亚洲av国产av综合av卡| av线在线观看网站| 成人特级av手机在线观看| 大香蕉久久网| 曰老女人黄片| 看非洲黑人一级黄片| 老女人水多毛片| 国产精品一区二区在线不卡| 亚洲成色77777| 久久久久网色| 日韩中字成人| 亚洲第一av免费看| 校园人妻丝袜中文字幕| 免费观看在线日韩| 麻豆精品久久久久久蜜桃| 精品久久久久久电影网| 国产精品国产av在线观看| 午夜福利视频精品| 日韩一本色道免费dvd| 狠狠精品人妻久久久久久综合| 黑人高潮一二区| 色5月婷婷丁香| 最近手机中文字幕大全| 只有这里有精品99| 老女人水多毛片| 成人漫画全彩无遮挡| 欧美日韩综合久久久久久| 国产中年淑女户外野战色| 草草在线视频免费看| 久久精品熟女亚洲av麻豆精品| 久久综合国产亚洲精品| 亚洲第一区二区三区不卡| 妹子高潮喷水视频| 中国三级夫妇交换| 伦精品一区二区三区| 一本大道久久a久久精品| 国产精品偷伦视频观看了| 丁香六月天网| 一级毛片久久久久久久久女| 久久人妻熟女aⅴ| 成人二区视频| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲四区av| 欧美日韩视频精品一区| 成人综合一区亚洲| 日日爽夜夜爽网站| 在线观看av片永久免费下载| av线在线观看网站| 中国国产av一级| 高清午夜精品一区二区三区| 国产69精品久久久久777片| 国产成人免费观看mmmm| 午夜福利,免费看| 午夜福利影视在线免费观看| 日本av免费视频播放| 久久久a久久爽久久v久久| 少妇裸体淫交视频免费看高清| 成年美女黄网站色视频大全免费 | 少妇高潮的动态图| 最近2019中文字幕mv第一页| 亚洲伊人久久精品综合| 亚洲精品一二三| 日韩免费高清中文字幕av| 国产午夜精品久久久久久一区二区三区| 又黄又爽又刺激的免费视频.| 亚洲情色 制服丝袜| 肉色欧美久久久久久久蜜桃| 丝袜在线中文字幕| 有码 亚洲区| 色94色欧美一区二区| 免费在线观看成人毛片| 黑人高潮一二区| 久久鲁丝午夜福利片| 亚洲欧美中文字幕日韩二区| 亚洲国产av新网站| 国产精品人妻久久久久久| 国产在视频线精品| 国产午夜精品久久久久久一区二区三区| 99热这里只有是精品50| 高清不卡的av网站| 亚洲精品日韩在线中文字幕| 免费黄频网站在线观看国产| 最近中文字幕高清免费大全6| 在线看a的网站| 亚洲av不卡在线观看| 久久久精品94久久精品| 全区人妻精品视频| 国产一级毛片在线| 夫妻午夜视频| 美女脱内裤让男人舔精品视频| 少妇人妻精品综合一区二区| 一本—道久久a久久精品蜜桃钙片| 不卡视频在线观看欧美| av女优亚洲男人天堂| 亚洲国产av新网站| 人人妻人人爽人人添夜夜欢视频 | 精品国产一区二区三区久久久樱花| 国产乱来视频区| 中文精品一卡2卡3卡4更新| h日本视频在线播放| 国产亚洲一区二区精品| 国产亚洲欧美精品永久| 女人久久www免费人成看片| 日韩一区二区三区影片| 国产男人的电影天堂91| 亚洲经典国产精华液单| 啦啦啦啦在线视频资源| 午夜激情久久久久久久| 成人免费观看视频高清| 日本黄色片子视频| 午夜久久久在线观看| 美女脱内裤让男人舔精品视频| 婷婷色av中文字幕| 国产亚洲午夜精品一区二区久久| 成人亚洲精品一区在线观看| 性色av一级| 亚洲欧洲日产国产| 一级毛片黄色毛片免费观看视频| 久久久久久久大尺度免费视频| 噜噜噜噜噜久久久久久91| 久久99一区二区三区| 精品亚洲乱码少妇综合久久|