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

    基于檢測與跟蹤相互迭代的極暗弱目標搜索算法

    2021-11-05 01:30:06尹增山
    計算機應用 2021年10期
    關鍵詞:波包檢測器時域

    効 琦,尹增山*,高 爽

    (1.中國科學院微小衛(wèi)星創(chuàng)新研究院,上海 201203;2.上??萍即髮W信息科學與技術學院,上海 201210;3.中國科學院大學,北京 100049)

    0 引言

    暗弱運動小目標的檢測在控制制導、軌道規(guī)劃等航空航天的應用[1-3]中是極為關鍵的技術。在目標檢測[4]算法中,由于運動目標的小型化、高速化以及對更大范圍檢測的進一步要求,使得成像的小目標像素數(shù)量極少,圖像信噪比(Signalto-Noise Ratio,SNR)極低。此時暗弱小目標被強噪聲干擾,信息量不足導致二者難以區(qū)分。

    在目標檢測中可以利用目標的兩種信息:一是運動學特征,包括速度、方向和距離等;二是視覺特征,主要有紋理、顏色、尺寸和外形。暗弱目標,即只有運動學特征的信息可以利用的目標,所以基于利用視覺信息的幀差法[5]和三幀差法不適用于暗弱目標檢測,并且背景差分法[6]也會由于SNR 極低的原因無法分辨目標與構建的背景之間的差異,效果不理想。另外,利用單幀空域信息的目標檢測算法也需要視覺特征的信息[7]。近年來基于時空聯(lián)合[8]和時頻分析[9-11]的各種檢測算法被學者們不斷提出。這些算法利用了連續(xù)數(shù)幀的時域信息,可以檢測大約SNR 為3 dB 左右的目標,滿足絕大多數(shù)弱小目標的檢測任務。文獻[12]提出了一種基于時域高階的運動點目標檢測方法,該方法將其定義為基于雙譜均值和標準差的多變量測試問題。對于極低信噪比(extremely Low Signal-to-Noise Ratio,LSNR)的暗弱目標,檢測前跟蹤(Track Before Detect,TBD)[13]是目前比較流行的一種方法。TBD 算法利用目標在時域序列中的連續(xù)性,挑選和累積疑似軌跡,最后選中軌跡真實性最高的目標。不過TBD 計算復雜度極高,尤其在視頻尺寸大或者幀數(shù)多時更為嚴重。為降低計算復雜度,基于動態(tài)規(guī)劃(Dynamic Programming,DP)的TBD 算法DP-TBD 算法[14-15]被提出。傳統(tǒng)的目標檢測算法由于需要目標的視覺信息,無法檢測暗弱小目標?;跁r空聯(lián)合、時頻分析的方法和TBD 沒有充分挖掘目標的時域信息,只用數(shù)幀或者低幀率的圖像做檢測存在一定的局限性。

    要解決對LSNR 目標的檢測問題,首先要從增加目標的先驗信息考慮,可從提升空間分辨率和時間分辨率兩方面考慮,增加在成像過程中暗弱小目標的信息量。由于星載或機載攝像機載荷的限制,相較于增加鏡頭直徑,選擇具有更高采集速率的互補金屬氧化物半導體(Complementary Metal Oxide Semiconductor,CMOS)攝像機更加經(jīng)濟可行。高頻成像過程可以提供高時相信息,并且利用更多的幀數(shù)還可以提供更為充分的軌跡信息。綜合考慮以上現(xiàn)有算法在暗弱目標檢測效果方面的不足和高頻成像的優(yōu)勢,本文借鑒TBD 算法利用軌跡連續(xù)性的思想,提出了一種時空域中檢測與跟蹤相互迭代(Detection and Tracking mutual Iteration,DTI)的算法來對空間中運動的暗弱小目標進行檢測。

    DTI 運用以時間域信息補充空間信息的思想:使用高頻相機捕獲運動小目標,可以增加目標在像元處的駐留時間,利用相對較多的幀數(shù)可以構建更為準確的時域背景估計[6],累積更多待測目標的能量。這從本質(zhì)上擴展了先驗信息,解決了現(xiàn)有算法無法解決的先驗信息不充分導致的目標與背景不可分問題。另外,DTI 在跟蹤過程中利用動態(tài)規(guī)劃(DP)算法極大地縮小了搜索空間的尺度,解決了TBD 算法計算量大、不易實現(xiàn)的問題;并且,檢測與跟蹤的過程迭代進行,可以更加充分地交互利用時空域信息,以達到對LSNR 目標的探測的目的。

    1 數(shù)據(jù)模型定義及預處理

    1.1 檢測數(shù)據(jù)模型

    包含暗弱小目標的可見光成像圖序列f(x,y,k)可以被式(1)描述:

    式中:x、y表示成像平面的水平軸和豎直軸;k是視頻的幀序號;f、fs、fb和fn分別是觀測信號、目標信號、背景信號和噪聲,背景fb是一個緩慢變化的二維隨機過程。凝視拍攝的模式中背景變化極小,通常占據(jù)了整個觀測信號f的低頻部分。噪聲分量fn是傳感器電路產(chǎn)生的各類噪聲的總和,fn和fb不相關且在空間分布是隨機的。

    小目標的亮度模式可由光學系統(tǒng)的點擴散效應[16]建立,它的亮度分布遵從點擴散函數(shù)(Point Spread Function,PSF)約束。由能量守恒可知光子在像素的分布比例總和為1:

    其中h(x,y)是光子在感光器上的分布比例。由光學理論可知,衍射現(xiàn)象中能量的大約84%集中在中央艾里斑,其余能量依次分布在衍射環(huán)上。艾里斑強度分布服從1 階貝塞爾函數(shù)[17],強度分布可以由PSF表示:

    其中B和b由式(4)計算:

    其中:F為焦距,f為光圈大小,x是感光器上成像位置與目標中心的相對距離,J1(bx)為一階貝塞爾函數(shù)。本文用簡化的PSF仿真,如式(5):

    其中:τ是目標的強度幅值;Δx、Δy分別為目標在x、y方向的擴散寬度。

    1.2 基于自適應小波包的去噪

    圖像預處理可以改善原圖的SNR,此步驟可以明顯地提高檢測算法的檢測效率。預處理過程的主要目標是減少雜波信號、抑制背景、增強目標信號,可以選擇合適的空域圖像預處理方法。本文檢測器運行在時序一維信號中,因此本節(jié)主要介紹考慮信號的時域特征。本文使用小波包重建的方式對檢測數(shù)據(jù)逐像素地在一維進行去噪,依據(jù)目標和背景噪聲的小波系數(shù)的不同,對噪聲的小波系數(shù)進行抑制使得SNR 得到提升。

    1.2.1 小波包變換的去噪原理

    傳統(tǒng)頻率分析,如傅里葉變換、Z 變換等存在時頻信息不完整的問題。利用小波變換[18]可以得到不同頻率信息的信號分量在時域的分布情況,它使用一組正交的、迅速衰減的基函數(shù)對原信號進行擬合,基函數(shù)可通過其尺度變量和平移變量以獲得不同的頻率和時間位置。

    如圖1 所示,小波包分解進一步將小波變換沒有分解的高頻部分D1和DA2子頻帶空間進行分解,而正是這種多層次的分解提高了時頻分辨率。圖中信號s可由第1 級分解的系數(shù)集A1和D1重建恢復,也可由第2級分解系數(shù)集AA2、DA2、AD2、DD2 重建,類似地也可由第3 級分解系數(shù)集重建。分解層次越深,某一系數(shù)子集與原信號s中某種分量的對應關系越精細。一般而言,噪聲信號對應于分解系數(shù)高頻子帶,選擇合適策略減小該高頻子帶系數(shù)可抑制原信號s中的噪聲。

    圖1 小波包三層分解結(jié)構Fig.1 Three-layer decomposition structure of wavelet packet

    1.2.2 自適應小波包去噪

    圖像序列信號的時域表示,目標經(jīng)過特定像元時造成瞬態(tài)沖擊?;谛〔ò撝捣指畹娜ピ胨惴ǖ暮诵狞c是選擇合適的分割閾值,使噪聲分量盡量得到抑制而保留沖擊目標信號。令式(6)為觀測得到的時域一維信號的小波包變換的系數(shù)集。

    其中:k表示分解層數(shù),j則是該層的第j個系數(shù)。根據(jù)D可以無失真地重建被分解信號,記Dr為D經(jīng)過噪聲系數(shù)抑制的結(jié)果,則dr為小波包系數(shù)集中的元素?,F(xiàn)有的硬閾值和軟閾值小波系數(shù)分割函數(shù)分別為:

    其中:thr是手動指定的分割閾值,關于thr的選取,有通用閾值規(guī)則、啟發(fā)式閾值規(guī)則和極大極小準則閾值規(guī)則等。以上閾值選擇策略無法解決thr太大則目標信號被削弱、thr太小則噪聲抑制效果不顯著的矛盾。

    考慮到待測目標的沖擊特性,本文根據(jù)該特性提出了自適應調(diào)整的閾值分割策略,通過計算沖擊特性的某個特征值來自適應調(diào)整閾值變化,從而達到抑制噪聲的小波系數(shù)、增強目標的小波系數(shù)的目的。鑒于峭度這個統(tǒng)計量對沖激信號十分敏感,所以選擇信號的峭度因子作為表征沖擊大小的特征值。峭度因子是歸一化的四階中心矩,表達式為:

    定義基于Kdj的調(diào)整因子:

    假設小的都是由于噪聲引起的?;谶@個假設,高頻子帶的小波系數(shù)應賦予更小的調(diào)整因子。給出以峭度因子作為統(tǒng)計特征的自適應閾值分割表達式為:

    式(11)后面的指數(shù)因子是為了減小如軟閾值和硬閾值分割中截斷效應的影響,使得在固定調(diào)整因子αj的情況下,分割曲線的映射更加平滑。圖2 給出了本文提出的自適應分割函數(shù)(ada thr)和硬閾值分割函數(shù)(hard thr)、軟閾值分割函數(shù)(soft thr)的比較(固定調(diào)整因子a=1),更加平滑的閾值分割函數(shù)對原信號的小波重建系數(shù)畸變更小,因此對原信號中目標信號的成分影響更小,dr指經(jīng)過閾值分割函數(shù)變換后的系數(shù)值。

    圖2 閾值分割函數(shù)的對比Fig.2 Comparison of threshold segmentation functions

    2 DTI算法

    2.1 DTI算法框架

    本文假設待檢測目標的數(shù)據(jù)由攝像機凝視拍攝獲得。鑒于使用高幀頻相機捕獲待測目標,可假定目標做勻速直線運動。檢測系統(tǒng)整體框架由預處理模塊、檢測模塊和軌跡關聯(lián)(跟蹤)模塊組成,后兩者在運行過程中互相迭代進行??梢猿橄蟮赜肈和T分別標記檢測與跟蹤模塊,用A標記檢測器參數(shù)。記第n階段的檢測結(jié)果和跟蹤結(jié)果分別為Dn、Tn,且檢測器參數(shù)集合為An,則這三者在整個算法運行過程的依賴關系可以由式(13)解釋:

    式(13)表明現(xiàn)階段的跟蹤結(jié)果由此階段的檢測器輸出決定。后一階段檢測結(jié)果由被檢測器參數(shù)An+1決定,而An+1是由前一階段的跟蹤結(jié)果Tn決定的。

    本文算法中檢測與跟蹤過程迭代,時空域信息交互利用,最后輸出迭代最優(yōu)。檢測輸出作為跟蹤模塊的輸入,跟蹤模塊進行軌跡關聯(lián)并自適應調(diào)整軌跡處的檢測器參數(shù)。通過設置一個衰減因子來根據(jù)第n輪軌跡關聯(lián)的結(jié)果調(diào)整第n+1 輪檢測器參數(shù)An+1中的判決門限。調(diào)整策略是在第n+1 輪檢測中令檢測器傾向于把疑似目標軌跡延伸處的信號片段判定為目標信號經(jīng)過的片段。而疑似目標軌跡的延長(軌跡的能量累積),又使得該軌跡自己在下一輪軌跡的篩選中被保留下來的可能性增大。

    根據(jù)軌跡關聯(lián)的結(jié)果判斷所有候選軌跡的真實性,若判斷該軌跡是由偽目標造成,則需重新初始化檢測窗口的參數(shù),直到該參數(shù)超出預設的范圍值,而該范圍的邊界則與對真實目標的運動速度有關。圖3給出了本文算法流程的整個框架。

    圖3 DTI算法框架Fig.3 DTI algorithm framework

    DTI算法運行流程:

    1)輸入:多幀圖像f(x,y,k);

    2)從i=1,2,…,k對每幀f(x,y,i)預處理;

    3)從i=1,2,…,x,從j=1,2,…,y對每個像素點f(i,j,k)進行小波包去噪;

    4)初始化或者重置檢測器的參數(shù)集A;

    5)從i=1,2,…,k根據(jù)A檢測圖像數(shù)據(jù)f(x,y,i)得到檢測標記數(shù)據(jù);

    6)對檢測標記數(shù)據(jù)進行霍夫線檢測得到待選軌跡,將這些軌跡保存在表格中;

    7)對表格中各條軌跡的能量進行累積計算,并根據(jù)能量計算值大小動態(tài)地調(diào)整檢測器分割閾值Th;

    8)轉(zhuǎn)至步驟5),直到表格中的條目為1;

    9)如果軌跡判斷為真,則輸出結(jié)果,如果A超出范圍,則DTI算法退出;如果不滿足上面兩條,則轉(zhuǎn)至步驟4)。

    2.2 檢測模塊

    檢測模塊中,檢測算子計算窗口內(nèi)當前信號片段與背景的差別時使用核函數(shù)構建。本節(jié)簡單介紹核函數(shù)的基本原理和基于核函數(shù)的檢測模塊。

    2.2.1 核函數(shù)原理

    數(shù)據(jù)在低維空間中線性不可分時,可以用增加特征維度的方法在高維空間中尋找超平面劃分這些數(shù)據(jù)。利用核函數(shù)[19-20]則直接在輸入空間中計算映射到高維特征空間H(Hilbert空間)中的兩個數(shù)據(jù)的內(nèi)積。

    其中,Φ:x∈χ→Φ(x) ∈H為從輸入空間到特征空間的映射函數(shù)。

    2.2.2 基于核函數(shù)的時域檢測架構

    圖4 DTI的檢測器Fig.4 Detector of DTI

    運用假設檢驗的思想構建各像素采集點一維時域信號的檢測器,然后將檢測器逐像素地由空域順序遍歷全部像素點以完成一個階段的檢測過程。

    如圖5 所示,對信號Fxy(t)設置了長度為l的檢測窗口,檢測窗口的特征統(tǒng)計記為fl(x),則該窗口的直接后繼窗口為fl(x+d)(需要注意的是,這里的x是指某一像素點時間維t上的某個位置,不是該像素空間域的橫軸坐標),其中,l表示窗口長度,x為窗口的起始坐標,d(0

    圖5 中讓信號延遲d和窗口長度l相等,兩個檢測窗口的特征統(tǒng)計關聯(lián)為式(5),表示相鄰兩個信號片段的統(tǒng)計變化量:

    圖5 時域一維信號的檢測Fig.5 Detection of one-dimensional signal in time domain

    為了方便描述,可將式(15)簡潔為式(16):

    由希爾伯特核空間理論可知,fl(x,x+d)可以構造為核函數(shù)形式。由于暗弱小目標的SNR 極低,目標信號和背景噪聲在幅度上線性不可分,假設檢驗中簡單的閾值分割無法準確判別信號和背景,所以利用核技巧以達到在高維空間中解決線性不可分的目的。令y=x+d,給出單像素的特征統(tǒng)計表達式由高斯核函數(shù)的形式所定義:

    式(17)也可以被認為是坐標分別為x,y處的一種廣義距離度量。則在窗口長度l上的歸一化統(tǒng)計量為:

    至此,知道了兩個檢測窗口可以獲得一個關聯(lián)的特征統(tǒng)計值ki(i=1,2,…,n):

    這里假設相機拍攝了n幀圖像,i表示第i個統(tǒng)計特征值。將兩個檢測窗口同時在時域向t:1 →end 的方向滑動一個時域坐標直至最后一幀可以獲得一個n維的特征向量K={k1,k2,…,kn},則該像素背景的統(tǒng)計特征模型定義為K的期望,如式(20)。

    此檢測算法運行在時域,逐像素地進行檢測。將所有視頻幀的同一像素點的時序信號看成一維信號f(t),把檢測目標是否經(jīng)過、何時經(jīng)過該感光像素的模型簡化為檢測一維信號中未知出現(xiàn)位置、未知長度的信號檢測問題。把這個一維信號檢測問題作為一個假設檢驗問題,記有像素區(qū)域內(nèi)不存在目標(H0)和存在目標(H1)這兩種假設:

    式(22)是對式(1)在時域一維的抽象表示,f(t)、x(t)、b(t)和n(t)分別表示時域中的觀測信號、目標信號、背景信號和噪聲信號。利用凝視拍攝模式中背景恒定和噪聲的分布在短時間內(nèi)變化不大的特性,檢測目標是否存在。當有目標經(jīng)過時則會造成沖擊,改變原有的信號分布,使得式(23)成立:

    其中:Kt表示一種廣義的信號特征提取機制,本文中可以用特征向量K實現(xiàn);t1和t2是不同的兩個相鄰時刻,不等號左側(cè)用Kt提取觀測信號在這兩個時刻用一階差分形式關聯(lián)的統(tǒng)計特征。在不存在目標的情況下,一般來說各個時刻提取的統(tǒng)計值變化不大。當出現(xiàn)目標信號時必然使得式(23)不等號成立。DTI 檢測器利用目標會改變Kt提取的特征值的特性進行判斷檢測。在一定置信度準則下有式(24)~(25)成立:

    其中:B表示對該像素一維時域信號的背景估計,可以用各檢測窗口核空間相關性的期望計算,即式(20)。換句話說,在式(23)中,用Kt提取所有時刻關聯(lián)的特征值并進行平均就是B的估計值。ΔB是特征向量K的方差估計:

    Th是可預設初始值的分割閾值,較大的Th提高檢測率,但也會增加虛警率。注意,Th值是基于跟蹤模塊軌跡關聯(lián)與篩選結(jié)果在不同像素上自適應調(diào)整的參數(shù)。

    檢測窗口參數(shù)設置:假設待測目標的直徑為R,與攝像機距離為D,運動速度為V,且攝像機焦距為f。根據(jù)相機的光學成像原理,可得目標像素直徑和像素速度分別為:

    令該攝像機的采集幀頻為pHz,則目標信號在時域的持續(xù)時長為:

    窗口尺寸嚴重依賴于目標在時域信號的持續(xù)時長,如果二者不匹配則會極大降低信號對背景的有效區(qū)分性,故該檢測系統(tǒng)需要目標大致的運動速度和尺寸的范圍作為先驗條件。

    檢測模塊算法流程 首先,初始化檢測參數(shù);然后,用核函數(shù)提取每個檢測窗口的特征,在時域上滑動窗口,得到背景特征的估計;最后,將每個檢測窗口提取的特征與背景估計做閾值分割,得到檢測結(jié)果。

    2.3 跟蹤模塊

    2.3.1 基于動態(tài)規(guī)劃的跟蹤原理

    基于動態(tài)規(guī)劃[12]的跟蹤策略依賴于:

    1)沿著目標軌跡能量累積最大。

    2)只需要令剩下軌跡能量最大的最優(yōu)決策可以使得整個目標軌跡的能量最大化,而不用考慮過去的狀態(tài)和決策。

    令Ti表示第i幀中目標的坐標[x,y],定義目標的軌跡能量為:

    其中E(Ti)可以用該坐標處像元亮度表征。軌跡的能量依賴于軌跡的長度和軌跡內(nèi)各坐標點的像元亮度。式(30)將一個多階聯(lián)合決策過程弱化為多個一階決策過程。DP 將聯(lián)合過程分段化,有助于解決窮舉搜索整個決策空間而造成的計算量大的問題。窮舉搜索M×M×N的待測數(shù)據(jù),計算量正比于M2N,而使用動態(tài)規(guī)劃解決這個問題計算量與M2N成正比。

    狀態(tài)轉(zhuǎn)移是動態(tài)規(guī)劃算法里的關鍵環(huán)節(jié),也決定了動態(tài)規(guī)劃算法的優(yōu)劣。DTI 狀態(tài)轉(zhuǎn)移指將第n輪跟蹤過程中軌跡關聯(lián)的結(jié)果作為第n+1 輪檢測的輸入,第n+1 輪檢測器的輸出即是第n輪檢測的下一狀態(tài)。所以從動態(tài)規(guī)劃的角度可以定義檢測過程D、跟蹤過程T和參數(shù)更新規(guī)則A為狀態(tài)轉(zhuǎn)移方程。綜合式(13),有:

    基于DP-TBD 算法的目標可能出現(xiàn)的區(qū)域(狀態(tài)轉(zhuǎn)移的可達狀態(tài))一般以目標的當前坐標為中心、由目標運動速度等其他因素影響形成的環(huán)形區(qū)域。圖6 展示了DTI 狀態(tài)轉(zhuǎn)移和DP-TBD 狀態(tài)轉(zhuǎn)移的差異,圖中DTI的狀態(tài)轉(zhuǎn)移由關聯(lián)的軌跡的更迭表示。

    圖6 動態(tài)規(guī)劃的狀態(tài)轉(zhuǎn)移區(qū)域Fig.6 State transition region in dynamic programming

    2.3.2 跟蹤過程的軌跡關聯(lián)準則

    利用霍夫直線檢測算法檢測所有疑似目標的軌跡。將這些軌跡以能量從高到低的順序排列,并將這些信息以索引表的形式儲存起來。

    由圖7 所示,以軌跡的能量高低進行排序。根據(jù)被保留軌跡的起止點得到該軌跡的三維空間表達式。下一輪檢測過程中只檢測被保留軌跡覆蓋的像素。所有待檢測像素的分割標記值和特征向量Ki在第1 次做檢測時已經(jīng)被計算過,所以后續(xù)階段的檢測過程不需要重新計算。故只需將檢測器參數(shù)A根據(jù)軌跡關聯(lián)結(jié)果更新,本文采用衰減檢測器的判決閾值Th的策略:

    圖7 軌跡關聯(lián)中的軌跡篩選Fig.7 Trajectory selection in track association

    其中:Th,n是指第n輪檢測過程中某個軌跡所覆蓋的所有坐標處的局域檢測分割閾值。衰減因子α調(diào)節(jié)每條軌跡、每輪跟蹤過程中軌跡覆蓋點坐標的通過率。

    對于檢測閾值Th,每個軌跡以及每個軌跡的每個檢測階段都有不同的取值。后續(xù)迭代只需再做一次簡單的閾值分割。顯然,相較于第1 次粗檢測,后續(xù)檢測過程的計算量被極大地減少了。

    實際拍攝數(shù)據(jù)中由于LSNR 的原因,即使霍夫線檢測器中,能量累積最高的坐標點(霍夫變換空間中點對應原像平面中一條線)對應的軌跡中也有很多檢測缺失點。本文的核心理念是利用檢測與跟蹤過程相互迭代,交互利用空間域和時間域的信息得到可信度較高的目標軌跡。重復檢測與跟蹤過程直到將表格里的軌跡數(shù)量消除到只剩一條。最后判斷這條軌跡的真實性:若得到的軌跡不可信,則在提供的目標速度和尺寸范圍內(nèi)更新檢測窗口尺度,直到以一定步長遍歷整個區(qū)間;若最后還是沒找到可信任軌跡,則該檢測算法宣布無目標。

    3 實驗仿真與結(jié)果分析

    3.1 仿真數(shù)據(jù)分析

    3.1.1 一維信號的檢測

    由于算法的檢測器工作在像素的時域信號中,需要驗證去噪算法和檢測算法在一維信號中的檢測性能。在噪聲模式為高斯白噪聲、目標信號為時域高斯鐘形信號的仿真條件下,圖8 展示了自適應小波包去噪算法的SNR 增益??梢钥闯?,對于SNR 為0 dB 及以下的含噪聲信號,經(jīng)過去噪后SNR 可以提升大約4~6 dB。

    圖8 自適應小波包去噪的SNR增益Fig.8 SNR gain of adaptive wavelet packet denoising

    圖9 是在Matlab 信號處理工具箱內(nèi)置樣例信號leleccum上的表現(xiàn),可以看出軟閾值和硬閾值分割算法雖然可以消除高頻噪聲,但對沖擊信號沒有增強作用。自適應去噪算法對沖擊分量有極大的增強效果,在1 000、1 200 處沖擊被增強。實際應用中需先對信號做平滑處理,再使用該去噪算法。

    圖9 leleccum信號上不同小波閾值去噪算法的效果對比Fig.9 Effect comparison of different wavelet threshold denoising algorithms on leleccum signal

    模擬生成SNR為0的一維信號st,經(jīng)過基于小波包重建的自適應去噪算法獲得SNR 為4.95 dB 的信號,可由檢測器定位持續(xù)時長極短的目標信號片段。如圖10所示,對于0 dB左右的一維時域信號,檢測器能以較高的檢測率正確定位目標的位置。目標SNR 極低,無法直接在原信號中發(fā)現(xiàn)位置,圖10 用一根豎直的定位虛線在仿真的目標下面作位置標記,用水平的粗線標記檢測結(jié)果,矩形框標注目標所在的位置??梢钥吹轿唇?jīng)去噪的信號沒有檢測到目標。

    圖10 一維信號的去噪與檢測Fig.10 Denoising and detecting for one-dimensional signal

    3.1.2 仿真數(shù)據(jù)的檢測

    仿真生成1 000 幀100×100 的圖像構成視頻,該視頻中包含了一個從像平面左上角向右下角做勻速直線運動的點目標。該目標占1 pixel,其SNR 為0 dB,像素速度為0.1 pixel/s。圖11是從上述方法中生成的圖像序列中截取的一幀圖像。

    圖11 0 dB SNR仿真視頻的X-Y投影圖Fig.11 X-Y projection of simulation video under SNR of 0 dB

    圖12 為在其他參數(shù)不變,SNR 分別是0 dB、1 dB 和3 dB的情況下對仿真目標做檢測與軌跡跟蹤的展示。對仿真視頻(三維數(shù)據(jù))做檢測得到同尺寸的判別數(shù)據(jù),將判別數(shù)據(jù)在時間維度進行投影獲得像平面的空間軌跡如圖12 所示??梢园l(fā)現(xiàn),在0 dB 的LSNR 下,DTI 算法依舊可以檢測到明顯的軌跡輪廓;隨著SNR 的提升,對暗弱目標軌跡的檢測率會進一步提升,人眼觀察到的軌跡也更加明顯了。

    圖12 三種SNR下對仿真數(shù)據(jù)的檢測結(jié)果Fig.12 Detection results of simulation data with three SNRs

    DTI 算法與其他算法的比較結(jié)果如表1 所示。表1 中,Top-Hat[7]是一種只對圖像單幀進行形態(tài)學變換檢測目標的算法。背景差分[6]是利用連續(xù)少數(shù)幾幀圖像做圖像檢測的算法。這兩種算法設計簡單,結(jié)果非常不理想。與基于高階譜分析的高階統(tǒng)計量算法[12]、基于動態(tài)規(guī)劃的DP-TBD[14]和三維匹配濾波[11]等算法相比,本文算法DTI 能在較低虛警率的前提下以高檢測率檢測到目標。

    表1 0 dB SNR下不同算法檢測結(jié)果對比 單位:%Tab.1 Detection result comparison of different algorithms under SNR of 0 dB unit:%

    3.2 真實數(shù)據(jù)驗證

    為了驗證DTI 算法對實際暗弱小目標的檢測與跟蹤能力,在實驗室的暗室中拍攝小球作為待檢目標進行實測數(shù)據(jù)的驗證實驗。實驗設置在完全控制光照亮度的暗室中,拍攝自由降落的小物體,物體直徑為15 mm,拍攝距離為14.50 m。通過以上措施模擬LSNR 的天基空間環(huán)境。本實驗使用75 mm的長焦鏡頭和德國Mikrotron公司生產(chǎn)的一款高幀頻工業(yè)攝像機EoSens?Cube7 以1 000 fps(frames/s)進行拍攝,用獲得的480×480×3 000 的視頻數(shù)據(jù)測試DTI 對真實數(shù)據(jù)的檢測性能。

    圖13(a)為實拍數(shù)據(jù)的某幀截圖,添加的矩形框標示目標軌跡的區(qū)域。圖13(b)為截取矩形框位置的5 幀圖像拼接而成,為了目標的顯示清晰,對目標做了顯化處理。圖14(c)為對目標檢測與跟蹤后投影到二維的軌跡??梢园l(fā)現(xiàn),對于此LSNR 的小球目標,DTI 算法可以檢測到目標的運動軌跡;但是LSNR 下,軌跡中會有大量的缺失點,并且這種情況會在SNR小于0的條件下逐漸加劇。

    圖13 實拍小球的檢測Fig.13 Detecting for real captured small ball

    4 結(jié)語

    本文針對LSNR 暗弱小目標檢測提出了一種新的算法,該算法改進現(xiàn)有的TBD 算法,將檢測與跟蹤過程聯(lián)合、迭代以提高檢測能力。經(jīng)過仿真和實際數(shù)據(jù)的驗證分析驗證了DTI 算法的有效性。與TBD 算法相比較,在檢測率和虛警率上DTI算法的表現(xiàn)更優(yōu)。

    不過DTI 算法檢測窗口的長度依賴于待檢測目標的速度、尺度和距離等運動參數(shù)。這些先驗信息是LSNR 目標所不能提供的,因此本文采用首先估計待檢測目標運動參數(shù)的大致范圍,然后再滑動搜索該范圍的算法。為了避免算法運行之前人為指定運動參數(shù)偏差,所以接下來的工作是進一步解決自適應估計此參數(shù)的問題。

    猜你喜歡
    波包檢測器時域
    基于時域信號的三電平逆變器復合故障診斷
    測控技術(2018年11期)2018-12-07 05:49:02
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    測控技術(2018年8期)2018-11-25 07:42:08
    車道微波車輛檢測器的應用
    基于極大似然準則與滾動時域估計的自適應UKF算法
    一種霧霾檢測器的研究與設計
    基于小波包變換的電力系統(tǒng)諧波分析
    基于時域逆濾波的寬帶脈沖聲生成技術
    小波包理論與圖像小波包分解
    基于時域波形特征的輸電線雷擊識別
    電測與儀表(2015年2期)2015-04-09 11:28:50
    基于小波包的全信息解調(diào)方法及其應用
    高清午夜精品一区二区三区 | 国产精品久久视频播放| 在线观看av片永久免费下载| 久久精品影院6| 国内精品美女久久久久久| 欧美最黄视频在线播放免费| 亚洲在线自拍视频| 免费看av在线观看网站| 成人一区二区视频在线观看| 久久久国产成人精品二区| 99热全是精品| 插阴视频在线观看视频| 久久天躁狠狠躁夜夜2o2o| 久久6这里有精品| 我要搜黄色片| www.色视频.com| 亚洲在线自拍视频| 在线国产一区二区在线| 高清午夜精品一区二区三区 | 精品人妻视频免费看| 国产精品av视频在线免费观看| 天堂动漫精品| 亚洲精品国产成人久久av| 夜夜爽天天搞| 99热全是精品| 国产成人a区在线观看| 日日摸夜夜添夜夜添av毛片| 黄色日韩在线| 亚洲国产精品sss在线观看| 亚洲18禁久久av| 搡女人真爽免费视频火全软件 | 听说在线观看完整版免费高清| 亚洲av成人av| 听说在线观看完整版免费高清| 黄色视频,在线免费观看| 激情 狠狠 欧美| 欧美zozozo另类| 一个人免费在线观看电影| 校园人妻丝袜中文字幕| 美女xxoo啪啪120秒动态图| 免费观看人在逋| 91久久精品国产一区二区成人| 亚洲五月天丁香| 99久久无色码亚洲精品果冻| 日本一二三区视频观看| 国产探花在线观看一区二区| 国内精品美女久久久久久| 国产精品人妻久久久影院| 一级毛片我不卡| 尤物成人国产欧美一区二区三区| 精品一区二区三区视频在线观看免费| 亚洲熟妇熟女久久| 99热全是精品| 又黄又爽又免费观看的视频| 久久99热这里只有精品18| 夜夜爽天天搞| 国产成人a区在线观看| 亚洲成av人片在线播放无| 色哟哟哟哟哟哟| 国产av不卡久久| 日韩亚洲欧美综合| 午夜视频国产福利| 亚洲乱码一区二区免费版| 中国国产av一级| 美女 人体艺术 gogo| a级一级毛片免费在线观看| 亚洲激情五月婷婷啪啪| 搡老妇女老女人老熟妇| 午夜精品国产一区二区电影 | av天堂中文字幕网| 色吧在线观看| 别揉我奶头~嗯~啊~动态视频| 国产欧美日韩精品一区二区| 最好的美女福利视频网| 成人午夜高清在线视频| 色尼玛亚洲综合影院| 网址你懂的国产日韩在线| 精品国产三级普通话版| 久久国内精品自在自线图片| a级一级毛片免费在线观看| 男人舔奶头视频| 美女xxoo啪啪120秒动态图| 女人被狂操c到高潮| 国产精品99久久久久久久久| 精品一区二区免费观看| 国产欧美日韩一区二区精品| 在线播放无遮挡| 国产精品福利在线免费观看| 久久韩国三级中文字幕| 人妻丰满熟妇av一区二区三区| 午夜a级毛片| 日本免费a在线| 如何舔出高潮| 久久精品91蜜桃| 国产av不卡久久| 最近的中文字幕免费完整| 久久久久久久久中文| 亚洲精品日韩在线中文字幕 | 久久久成人免费电影| 91久久精品国产一区二区成人| 直男gayav资源| 成人亚洲欧美一区二区av| 日韩三级伦理在线观看| 中文字幕熟女人妻在线| 乱人视频在线观看| 成人无遮挡网站| 男女那种视频在线观看| 国产精品久久久久久亚洲av鲁大| 69人妻影院| 欧美国产日韩亚洲一区| 国产高潮美女av| 亚洲无线在线观看| 精品久久久久久成人av| 精品熟女少妇av免费看| 最近视频中文字幕2019在线8| 亚洲成a人片在线一区二区| 国产成人freesex在线 | 午夜日韩欧美国产| 日韩中字成人| 春色校园在线视频观看| 看黄色毛片网站| 人人妻人人澡欧美一区二区| 亚洲精品影视一区二区三区av| 亚洲精品乱码久久久v下载方式| 亚洲成人av在线免费| 精品乱码久久久久久99久播| 国产黄色视频一区二区在线观看 | 联通29元200g的流量卡| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 国产精品一区二区免费欧美| 国产精品无大码| 日韩av在线大香蕉| 免费观看精品视频网站| 久久国产乱子免费精品| 国产精华一区二区三区| 国产亚洲精品av在线| 黄片wwwwww| 91久久精品国产一区二区成人| eeuss影院久久| 91在线精品国自产拍蜜月| 中文字幕免费在线视频6| 三级男女做爰猛烈吃奶摸视频| 亚洲精品色激情综合| 能在线免费观看的黄片| 伊人久久精品亚洲午夜| 97超视频在线观看视频| 免费看av在线观看网站| 国产一级毛片七仙女欲春2| 有码 亚洲区| 又爽又黄无遮挡网站| 欧美最新免费一区二区三区| 国产精品三级大全| 亚洲av电影不卡..在线观看| 99久久久亚洲精品蜜臀av| 中国国产av一级| 免费黄网站久久成人精品| 亚洲电影在线观看av| 欧美性猛交黑人性爽| 如何舔出高潮| 日韩一本色道免费dvd| 我要搜黄色片| 国产探花在线观看一区二区| 热99在线观看视频| 麻豆国产av国片精品| 搡老熟女国产l中国老女人| 日本 av在线| 午夜免费激情av| 午夜精品一区二区三区免费看| 久久精品国产鲁丝片午夜精品| 婷婷精品国产亚洲av在线| 国产成人91sexporn| 日韩三级伦理在线观看| 亚洲av中文av极速乱| 国产白丝娇喘喷水9色精品| 一级毛片aaaaaa免费看小| 国产久久久一区二区三区| 日本a在线网址| 国产一区二区亚洲精品在线观看| avwww免费| 亚洲国产欧洲综合997久久,| 国产成人aa在线观看| 亚洲国产日韩欧美精品在线观看| 在线观看一区二区三区| 久久久久九九精品影院| 国产成人freesex在线 | 国产乱人偷精品视频| 日本-黄色视频高清免费观看| 亚洲精品粉嫩美女一区| 桃色一区二区三区在线观看| 免费av不卡在线播放| 99热网站在线观看| 麻豆乱淫一区二区| 高清午夜精品一区二区三区 | 日韩 亚洲 欧美在线| 一个人看视频在线观看www免费| 一区福利在线观看| 久久精品夜色国产| 在线a可以看的网站| 精品久久久久久久久av| 成人亚洲精品av一区二区| 日本黄色视频三级网站网址| а√天堂www在线а√下载| 久久久a久久爽久久v久久| 欧美区成人在线视频| 波多野结衣高清无吗| 国产不卡一卡二| 国产精品无大码| av在线观看视频网站免费| 97碰自拍视频| 日日撸夜夜添| 久久综合国产亚洲精品| 麻豆成人午夜福利视频| 波多野结衣巨乳人妻| 国产不卡一卡二| 国产午夜精品论理片| 国产精品av视频在线免费观看| 1024手机看黄色片| 在线观看66精品国产| 亚洲精品粉嫩美女一区| 如何舔出高潮| 亚洲激情五月婷婷啪啪| 免费观看的影片在线观看| 亚洲欧美成人精品一区二区| 日本精品一区二区三区蜜桃| 人妻久久中文字幕网| 国产极品精品免费视频能看的| 国产精品精品国产色婷婷| 亚洲人成网站在线播| 我的老师免费观看完整版| 国产一区二区在线av高清观看| 美女内射精品一级片tv| 乱系列少妇在线播放| 亚洲不卡免费看| 国产高清三级在线| 又黄又爽又刺激的免费视频.| 看片在线看免费视频| 国产黄片美女视频| 国产av在哪里看| 亚洲av熟女| 国内精品久久久久精免费| 成人三级黄色视频| av女优亚洲男人天堂| 黄色配什么色好看| 淫秽高清视频在线观看| 亚洲在线自拍视频| 婷婷亚洲欧美| 婷婷精品国产亚洲av在线| 亚洲自拍偷在线| 久久亚洲精品不卡| a级毛片免费高清观看在线播放| 色播亚洲综合网| 麻豆精品久久久久久蜜桃| 国产精品永久免费网站| 亚洲无线观看免费| 日韩人妻高清精品专区| 欧美又色又爽又黄视频| 国产高清激情床上av| 亚洲无线在线观看| 一本一本综合久久| av在线老鸭窝| 麻豆久久精品国产亚洲av| eeuss影院久久| 99久久精品热视频| 国产一区二区在线av高清观看| 国产亚洲精品久久久com| 国产日本99.免费观看| 蜜桃亚洲精品一区二区三区| 啦啦啦观看免费观看视频高清| 少妇熟女aⅴ在线视频| 国产三级中文精品| 少妇的逼好多水| 男女之事视频高清在线观看| 亚洲熟妇中文字幕五十中出| 黄色一级大片看看| 麻豆一二三区av精品| 亚洲在线自拍视频| 深爱激情五月婷婷| 丰满乱子伦码专区| 九九热线精品视视频播放| 免费搜索国产男女视频| 国产麻豆成人av免费视频| 亚洲精品一卡2卡三卡4卡5卡| 综合色av麻豆| 天堂av国产一区二区熟女人妻| 久久99热这里只有精品18| 久久久国产成人精品二区| 精品人妻偷拍中文字幕| 国内精品宾馆在线| 99视频精品全部免费 在线| 日韩欧美免费精品| 看黄色毛片网站| 免费不卡的大黄色大毛片视频在线观看 | 国产久久久一区二区三区| 99国产极品粉嫩在线观看| 亚洲欧美清纯卡通| 国产精品1区2区在线观看.| 99九九线精品视频在线观看视频| 久久久欧美国产精品| 精品午夜福利在线看| 晚上一个人看的免费电影| 97碰自拍视频| 久久精品夜色国产| 亚洲激情五月婷婷啪啪| а√天堂www在线а√下载| 麻豆国产av国片精品| 精品99又大又爽又粗少妇毛片| 日韩在线高清观看一区二区三区| 在线观看av片永久免费下载| 淫秽高清视频在线观看| 秋霞在线观看毛片| 国产精品1区2区在线观看.| 精品免费久久久久久久清纯| 国产免费男女视频| 99九九线精品视频在线观看视频| 亚洲成av人片在线播放无| 亚洲熟妇熟女久久| 国产成人aa在线观看| 最近手机中文字幕大全| 国产真实伦视频高清在线观看| 日韩一本色道免费dvd| 亚洲av免费在线观看| 一级av片app| 神马国产精品三级电影在线观看| 国产精品美女特级片免费视频播放器| 精品一区二区三区人妻视频| 最近手机中文字幕大全| 欧美日韩一区二区视频在线观看视频在线 | 3wmmmm亚洲av在线观看| 亚洲精品国产成人久久av| 国产一区二区亚洲精品在线观看| 国产精品久久视频播放| 99精品在免费线老司机午夜| avwww免费| 成年女人毛片免费观看观看9| 看片在线看免费视频| 成年版毛片免费区| 日韩人妻高清精品专区| 99热只有精品国产| 99riav亚洲国产免费| 天堂网av新在线| 人人妻,人人澡人人爽秒播| videossex国产| 久久人人精品亚洲av| 婷婷色综合大香蕉| 国产av一区在线观看免费| 国产精品久久久久久亚洲av鲁大| 亚洲最大成人av| 亚洲av成人精品一区久久| 亚洲电影在线观看av| 综合色丁香网| 欧美一级a爱片免费观看看| 国产午夜精品久久久久久一区二区三区 | 99久久中文字幕三级久久日本| 国产精品伦人一区二区| 男人舔奶头视频| 久久久色成人| 久久久久久久久久成人| 久久久久久久久大av| 给我免费播放毛片高清在线观看| 日韩欧美国产在线观看| 99国产极品粉嫩在线观看| 精品久久久久久久末码| 赤兔流量卡办理| 国产不卡一卡二| 欧美xxxx性猛交bbbb| 欧美人与善性xxx| 亚洲18禁久久av| 久久人妻av系列| 香蕉av资源在线| 国产在视频线在精品| 国产蜜桃级精品一区二区三区| 免费黄网站久久成人精品| 最新中文字幕久久久久| 天堂影院成人在线观看| 日本精品一区二区三区蜜桃| 国产欧美日韩精品一区二区| 欧美一区二区国产精品久久精品| 日日摸夜夜添夜夜添av毛片| 少妇裸体淫交视频免费看高清| 熟妇人妻久久中文字幕3abv| 两个人视频免费观看高清| 黄色一级大片看看| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 国产人妻一区二区三区在| 99久久久亚洲精品蜜臀av| 99久久精品国产国产毛片| av.在线天堂| 亚洲四区av| 国内精品美女久久久久久| 好男人在线观看高清免费视频| 亚洲中文字幕一区二区三区有码在线看| 午夜精品在线福利| 国产精品综合久久久久久久免费| aaaaa片日本免费| 日产精品乱码卡一卡2卡三| 麻豆av噜噜一区二区三区| 日韩欧美精品v在线| 国产精品伦人一区二区| 亚洲美女视频黄频| 国产又黄又爽又无遮挡在线| 亚洲精品粉嫩美女一区| 尾随美女入室| 两个人视频免费观看高清| 亚洲av二区三区四区| 18禁在线无遮挡免费观看视频 | 亚洲欧美中文字幕日韩二区| 嫩草影视91久久| 麻豆精品久久久久久蜜桃| 深夜a级毛片| 精品日产1卡2卡| 一a级毛片在线观看| 国产精品久久电影中文字幕| 99热只有精品国产| 午夜视频国产福利| 男女之事视频高清在线观看| av天堂中文字幕网| 久久午夜亚洲精品久久| 久久精品人妻少妇| 国产单亲对白刺激| 日本 av在线| 久久精品国产鲁丝片午夜精品| 国产精品国产三级国产av玫瑰| 啦啦啦观看免费观看视频高清| 国产一区二区亚洲精品在线观看| 特级一级黄色大片| 国产成人a∨麻豆精品| 国产熟女欧美一区二区| 一进一出抽搐动态| 久久人人精品亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 搡女人真爽免费视频火全软件 | 久久精品国产自在天天线| 亚洲国产精品久久男人天堂| 国产一区亚洲一区在线观看| 国产精品野战在线观看| 一个人免费在线观看电影| 少妇高潮的动态图| 毛片女人毛片| 蜜臀久久99精品久久宅男| 欧美成人精品欧美一级黄| 热99在线观看视频| 精品一区二区三区av网在线观看| 综合色av麻豆| 老熟妇乱子伦视频在线观看| 男人和女人高潮做爰伦理| 天天躁夜夜躁狠狠久久av| 免费看av在线观看网站| 久久久久久伊人网av| 亚洲在线观看片| 久久久精品94久久精品| 欧美绝顶高潮抽搐喷水| 国产精品久久视频播放| ponron亚洲| 性欧美人与动物交配| 久久精品久久久久久噜噜老黄 | 亚洲不卡免费看| 国产精品人妻久久久影院| 国产伦精品一区二区三区四那| 亚洲av一区综合| 亚洲精品一卡2卡三卡4卡5卡| 久久午夜亚洲精品久久| 2021天堂中文幕一二区在线观| 国产69精品久久久久777片| 婷婷色综合大香蕉| 国产在线精品亚洲第一网站| a级毛片免费高清观看在线播放| 日本熟妇午夜| 毛片一级片免费看久久久久| 12—13女人毛片做爰片一| 美女高潮的动态| 99久久精品热视频| 亚洲国产日韩欧美精品在线观看| 在线观看美女被高潮喷水网站| 网址你懂的国产日韩在线| 久久久色成人| 日韩欧美精品v在线| 亚洲国产精品国产精品| 97在线视频观看| 内射极品少妇av片p| 久久久久国内视频| 欧美+亚洲+日韩+国产| 国产v大片淫在线免费观看| 国产精品免费一区二区三区在线| 麻豆一二三区av精品| 两个人视频免费观看高清| 亚洲国产高清在线一区二区三| 秋霞在线观看毛片| 日韩制服骚丝袜av| 美女被艹到高潮喷水动态| 久久综合国产亚洲精品| 国产熟女欧美一区二区| 少妇熟女欧美另类| 久久亚洲精品不卡| 精品国内亚洲2022精品成人| 综合色丁香网| 精品熟女少妇av免费看| 日本撒尿小便嘘嘘汇集6| 国产伦精品一区二区三区视频9| 精品久久久噜噜| 欧美一区二区亚洲| 在线a可以看的网站| 欧美xxxx性猛交bbbb| 男女下面进入的视频免费午夜| 人妻丰满熟妇av一区二区三区| 超碰av人人做人人爽久久| 国产精品亚洲一级av第二区| 99热全是精品| 日日撸夜夜添| 中文资源天堂在线| 中国美女看黄片| 大型黄色视频在线免费观看| 亚洲欧美精品综合久久99| 日韩欧美三级三区| 成人永久免费在线观看视频| 一级毛片我不卡| 日韩av不卡免费在线播放| 免费黄网站久久成人精品| 有码 亚洲区| 欧美激情久久久久久爽电影| 在线免费十八禁| 国产一区二区三区av在线 | 欧美日韩一区二区视频在线观看视频在线 | 又爽又黄a免费视频| АⅤ资源中文在线天堂| 两性午夜刺激爽爽歪歪视频在线观看| 俺也久久电影网| 国产男靠女视频免费网站| 中文字幕人妻熟人妻熟丝袜美| 久久久久久国产a免费观看| 熟妇人妻久久中文字幕3abv| 男人和女人高潮做爰伦理| 97人妻精品一区二区三区麻豆| 国产精品一区二区免费欧美| 国产一区二区在线av高清观看| 午夜视频国产福利| 久久久久久大精品| 久久久久久久亚洲中文字幕| 日韩欧美一区二区三区在线观看| 久久久久久久久久久丰满| 99久久精品热视频| 蜜桃亚洲精品一区二区三区| 国产真实伦视频高清在线观看| 99热全是精品| 国产三级中文精品| 最近2019中文字幕mv第一页| 男女啪啪激烈高潮av片| 久久久久性生活片| 亚洲真实伦在线观看| 国产亚洲欧美98| 日本色播在线视频| 特大巨黑吊av在线直播| 国产高清激情床上av| 一本精品99久久精品77| 国产精品国产三级国产av玫瑰| 久久久欧美国产精品| 日韩欧美精品免费久久| av天堂在线播放| 男女那种视频在线观看| 亚洲七黄色美女视频| 欧美一区二区精品小视频在线| 亚洲精品乱码久久久v下载方式| 综合色av麻豆| 国产精品久久久久久久电影| 啦啦啦啦在线视频资源| 亚洲av免费在线观看| 少妇熟女aⅴ在线视频| 精品人妻熟女av久视频| 丰满的人妻完整版| 日韩成人伦理影院| 亚洲欧美清纯卡通| 五月伊人婷婷丁香| 成人精品一区二区免费| 俺也久久电影网| 亚洲国产精品sss在线观看| 99热精品在线国产| 亚洲高清免费不卡视频| 国产精品福利在线免费观看| 亚洲人成网站高清观看| 欧美zozozo另类| 天天躁夜夜躁狠狠久久av| 女人被狂操c到高潮| 少妇被粗大猛烈的视频| 精品人妻熟女av久视频| 亚洲三级黄色毛片| 亚洲性夜色夜夜综合| 国产免费男女视频| 亚洲高清免费不卡视频| 国产午夜福利久久久久久| 午夜激情欧美在线| 18禁黄网站禁片免费观看直播| 内地一区二区视频在线| 日日撸夜夜添| 中文字幕久久专区| 成年女人永久免费观看视频| 亚洲性久久影院| 亚洲婷婷狠狠爱综合网| 日日摸夜夜添夜夜添小说| av在线观看视频网站免费| 日本黄色片子视频| 天天躁夜夜躁狠狠久久av| 麻豆av噜噜一区二区三区| 久久这里只有精品中国| 亚洲人成网站高清观看| 亚洲精品国产成人久久av| 少妇高潮的动态图| 插阴视频在线观看视频| 日本五十路高清| 免费黄网站久久成人精品| 黑人高潮一二区| 尾随美女入室| 亚洲精品久久国产高清桃花| 国内精品美女久久久久久| 午夜爱爱视频在线播放|