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

    全極化雷達(dá)遙感影像的迭代優(yōu)化非局部均值去噪法

    2019-09-04 00:35:24馬曉雙吳鵬海
    測(cè)繪學(xué)報(bào) 2019年8期
    關(guān)鍵詞:協(xié)方差鄰域相似性

    馬曉雙,吳鵬海

    安徽大學(xué)資源與環(huán)境工程學(xué)院,安徽 合肥 230601

    全極化SAR(polarimetric SAR,PolSAR)系統(tǒng)作為傳統(tǒng)單極化SAR系統(tǒng)的一種更高級(jí)形式,能同時(shí)發(fā)射、接收水平及垂直的極化雷達(dá)波,能夠獲取地物更為豐富的后向散射信息,更利于地物場(chǎng)景的解譯,已在諸多行業(yè)領(lǐng)域得到廣泛應(yīng)用。特別是,我國(guó)于2016年發(fā)射了國(guó)內(nèi)首顆民用全極化SAR衛(wèi)星——高分3號(hào),進(jìn)一步推動(dòng)了PolSAR的應(yīng)用前景。然而,由于成像機(jī)理的限制,使得PolSAR影像不可避免地存在相干斑噪聲。相干斑嚴(yán)重降低了SAR影像的質(zhì)量,也降低了利用SAR數(shù)據(jù)進(jìn)行地物解譯的能力[1-2]。因此,對(duì)于PolSAR數(shù)據(jù)的后續(xù)應(yīng)用而言,相干斑濾波是必不可少的預(yù)處理程序。

    對(duì)于全極化SAR數(shù)據(jù)而言,噪聲不僅存在于極化協(xié)方差矩陣的對(duì)角線元素上,同時(shí)也存在于非對(duì)角線元素上。并且,對(duì)角線元素受到乘性噪聲的影響,而非對(duì)角線元素則受到乘性和加性混合的噪聲的影響[3]。這使得全極化SAR的去噪問題要比單極化SAR要復(fù)雜得多。早期的全極化SAR去噪的方法大都是基于文獻(xiàn)[4]提出的極化白化濾波器。之后,文獻(xiàn)[1]提出了一種基于線性最小均方誤差(linear minimum mean square error,LMMSE)原則的PolSAR去噪算法(增強(qiáng)Lee濾波)。這一方法極大地推動(dòng)了對(duì)極化SAR去噪的研究,也引領(lǐng)了之后一系列基于LMMSE濾波器的算法[5-6]。除此之外,基于二叉樹的算法[7]、各向異性擴(kuò)散算法[8-9]、變分正則化算法[10-11]、模擬退火算法[12]等也被應(yīng)用于極化SAR去噪。近年來,基于非局部均值(nonlocal means,NLM)思想的算法已成為了PolSAR影像去噪的研究熱點(diǎn)[13-15],并取得了較好的去噪效果。NLM算法最先應(yīng)用于數(shù)字圖像處理[16]。其基本思想是認(rèn)為影像中往往都會(huì)存在重復(fù)出現(xiàn)的或相似的影像塊,因此像素的估計(jì)值可由它的非局部搜索窗口中的所有像素加權(quán)平均后得到,而其中每一個(gè)像素的權(quán)重由它的鄰域窗口與待處理像素鄰域窗口的差異性所獲得。差異越大,權(quán)重越小。

    然而,傳統(tǒng)PolSAR非局部均值去噪算法在度量像素或影像塊極化相似度時(shí)基本利用到的都是含噪聲影像的某種統(tǒng)計(jì)特性,而沒有考慮到原始像素間的差異性。并且,許多PolSAR非局部均值算法僅僅考慮的是強(qiáng)度信息相似性,而未考慮包括極化相位差在內(nèi)的全極化信息,導(dǎo)致難以獲得穩(wěn)健準(zhǔn)確的相似性測(cè)度。有鑒于此,本文提出一種迭代優(yōu)化的PolSAR NLM算法。該算法的核心思想是:在每次迭代去噪過程中,同時(shí)考慮原始影像全極化噪聲統(tǒng)計(jì)特性和前一次迭代所得的相對(duì)低噪聲度影像的全極化信息,以完善影像塊間極化相似性的度量。試驗(yàn)部分的模擬PolSAR和真實(shí)PolSAR影像去噪結(jié)果都表明,本文提出的算法取得了較好的去噪效果,噪聲得到了有效的抑制,影像邊緣細(xì)節(jié)和極化散射特性都得到了較好的保持。

    1 非局部均值及PolSAR相干斑統(tǒng)計(jì)特性

    1.1 非局部均值去噪思想

    NLM去噪算法濾波過程的通式可表示為

    (1)

    式中,f為原始噪聲影像;w代表一個(gè)權(quán)重函數(shù)。NLM算法認(rèn)為目標(biāo)像素的估計(jì)值由它的非局部搜索窗口中的像素加權(quán)平均后得到,而其中每一個(gè)像素的權(quán)重由該像素的鄰域窗口與待處理像素鄰域窗口的整體差異性所計(jì)算獲得。圖1展示了NLM的基本思想:對(duì)于以像素p為中心的鄰域塊而言,由于以q1為中心的鄰域塊的內(nèi)部結(jié)構(gòu)與之非常相似,所以在對(duì)p進(jìn)行去噪時(shí),像素q1的權(quán)重w(p,q1)較大;而由于以q2和q3為中心的鄰域塊的內(nèi)部結(jié)構(gòu)與之相差較大,所以在對(duì)p進(jìn)行去噪時(shí),像素q2和q3的權(quán)重很小??梢钥吹剑捎贜LM算法不僅考慮了單一像素間影像值的接近性,還考慮到了兩個(gè)像素鄰域塊的結(jié)構(gòu)布局的相似性,這使其比起傳統(tǒng)的鄰域?yàn)V波方法更能取得穩(wěn)定有效的濾波效果。顯而易見,決定NLM算法濾波效果的一個(gè)重要因素就是影像塊間相似性的度量,本文的主要?jiǎng)?chuàng)新就是提出了一種迭代優(yōu)化的全極化影像塊相似性度量方式。

    圖1 NLM算法的思路Fig.1 The diagram of the idea of nonlocal means

    1.2 PolSAR相干斑統(tǒng)計(jì)特性

    眾所周知,對(duì)于SAR的強(qiáng)度數(shù)據(jù)而言,其噪聲是一種服從Gamma分布的乘性噪聲[17]。然而,全極化數(shù)據(jù)的噪聲模型要復(fù)雜得多。對(duì)于噪聲完全發(fā)育的視數(shù)為L(zhǎng)的PolSAR影像而言,其極化協(xié)方差矩陣CL服從如下的復(fù)Wishart分布[18]

    (2)

    2 迭代優(yōu)化的PolSAR NLM去噪算法

    要完善影像塊間極化相似性的度量,就必須要考慮到全極化信息。本文提出的PolSAR NLM算法的創(chuàng)新性是在每次迭代去噪過程中,同時(shí)考慮原始影像全極化噪聲統(tǒng)計(jì)特性和前一次迭代所得的相對(duì)低噪聲度影像的全極化信息。因此,如何度量原始影像塊全極化相似性及迭代所得影像塊全極化相似性是本文的主要關(guān)注點(diǎn)。

    2.1 迭代優(yōu)化的PolSAR影像塊相似性度量

    (3)

    式中,m=n為影像的視數(shù)。對(duì)以上公式取自然對(duì)數(shù)并簡(jiǎn)化常數(shù)項(xiàng),可得到

    lnQ=6ln2+ln|X|+ln|Y|-2ln|X+Y|

    (4)

    式中,lnQ的取值范圍是(-∞,0]。像素越相似,則lnQ越大。

    隨著影像迭代去噪過程的進(jìn)行,所得影像的噪聲程度會(huì)逐漸減弱。因此,在每次迭代的過程中,都可以利用上一次迭代所得影像的信息來完善像素的相似性測(cè)度。這一迭代優(yōu)化的相似性度量思想可以表示為

    (5)

    (6)

    接下來可利用貝葉斯定理對(duì)式(6)右邊第1項(xiàng)進(jìn)行變化。貝葉斯公式如下

    (7)

    式中,P(B|A)代表?xiàng)l件A發(fā)生的前提下B發(fā)生的概率。據(jù)此,有

    (8)

    (9)

    迭代相似性測(cè)度最終表示為

    (10)

    對(duì)于PolSAR數(shù)據(jù)而言,Wishart距離可表示為[18]

    (11)

    若考慮上述Wishart距離的對(duì)稱性,對(duì)其進(jìn)行調(diào)整即可得到如下的Kullback-Leibler散度距離度量公式[20]

    (12)

    最終可推得,式(12)所示的相似測(cè)度值仍為負(fù)值。值越大,代表像素間越相似。需要說明的是:本文提出的迭代去噪方法的主要思路是在每次迭代的過程中完善權(quán)重的計(jì)算,而迭代時(shí)仍然將像素的原始值同權(quán)重結(jié)合進(jìn)行加權(quán)平均,所以本質(zhì)上并不屬于“二次濾波”。

    值得注意的是,由于式(4)中涉及矩陣行列式的對(duì)數(shù),為了避免奇異值的出現(xiàn),就應(yīng)當(dāng)保證|C|≠0,也即保證滿秩性。然而,當(dāng)影像視數(shù)小于協(xié)方差矩陣的維度時(shí)(即L<3),影像中的一些像素的極化協(xié)方差矩陣將不再是滿秩的,此時(shí)就會(huì)導(dǎo)致數(shù)學(xué)計(jì)算的問題。為了解決這一問題,本文算法加入一種對(duì)待去噪影像進(jìn)行預(yù)處理的方法[22],保證影像所有像素的協(xié)方差矩陣達(dá)到滿秩。這一協(xié)方差矩陣的預(yù)變換步驟可以表示為

    (13)

    式中,i和j代表協(xié)方差矩陣中元素的位置;γ=min(L/3,1)。這一過程的實(shí)質(zhì)就是保持協(xié)方差對(duì)角元素(即各通道的強(qiáng)度值)不變,而根據(jù)影像的視數(shù)對(duì)非對(duì)角元素進(jìn)行線性變換。研究表明[22],以上步驟處理后的影像仍然滿足復(fù)Wishart分布,也即基于復(fù)Wishart分布的似然比檢驗(yàn)法依然適用于變換后的影像。這里需要說明的是,以上預(yù)處理過程只是用于計(jì)算每個(gè)像素的權(quán)重,而在加權(quán)平均過程中使用的仍然是原始的數(shù)據(jù)。因此,不會(huì)存在改變數(shù)據(jù)極化信息的弊端。

    2.2 權(quán)重及濾波參數(shù)的設(shè)置

    每次迭代濾波過程中,對(duì)于每個(gè)參考影像塊x而言,窗口內(nèi)其他影像塊y的權(quán)重為

    (14)

    式中,N為影像塊的像素?cái)?shù);T在每次迭代時(shí)為一個(gè)固定參數(shù),本文采用文獻(xiàn)[23]提出的“噪聲估算器”來確定T的取值:每一次迭代時(shí),計(jì)算影像上所有相鄰像素的相似度,構(gòu)建由這些相似度所組成的直方圖,然后計(jì)算直方圖的90%分位數(shù),即將此值當(dāng)作當(dāng)前時(shí)刻的參數(shù)T。最終,對(duì)于某個(gè)像素而言,其極化協(xié)方差矩陣的估計(jì)值為

    (15)

    式中,M為濾波窗口的大小。

    本文提出的PolSAR NLM算法有3個(gè)參數(shù)需要預(yù)先設(shè)定,分別是迭代次數(shù)t、影像塊大小N和濾波窗口大小M。理論而言,隨著迭代次數(shù)增加,影像的去噪效果會(huì)越加優(yōu)化。但是,迭代次數(shù)增加也會(huì)帶來計(jì)算時(shí)間的增加。試驗(yàn)發(fā)現(xiàn),對(duì)于大多數(shù)影像,當(dāng)?shù)螖?shù)達(dá)到3次時(shí),就能達(dá)到優(yōu)良的去噪效果且繼續(xù)迭代后的效果提升非常有限。因此,本文提出的PolSAR NLM算法的迭代次數(shù)固定為3次。對(duì)于NLM算法而言,影像塊大小通常設(shè)置為3×3到7×7之間,而濾波窗口大小則通常設(shè)置為13×13到21×21之間。研究表明[14],隨著影像塊和濾波窗口大小的增加,影像的噪聲會(huì)得到更明顯地抑制,但細(xì)節(jié)丟失會(huì)逐漸加重。因此,為了達(dá)到去噪與保持影像細(xì)節(jié)之間的折中,本文試驗(yàn)部分將算法的影像塊和濾波窗口大小分別設(shè)為5×5和17×17。

    3 試驗(yàn)結(jié)果和分析

    為了驗(yàn)證本文提出的算法的濾波效果,本文采用一幅模擬的PolSAR數(shù)據(jù)和兩幅真實(shí)的PolSAR影像進(jìn)行去噪的試驗(yàn)。作為對(duì)比,試驗(yàn)部分也展示了增強(qiáng)Lee濾波[1]、IDAN濾波[5]及Pretest NLM濾波算法[14]的去噪效果。為了達(dá)到去噪與影像細(xì)節(jié)保持的平衡,本文所有試驗(yàn)中增強(qiáng)Lee濾波的濾波窗口大小設(shè)為文獻(xiàn)[1]中所推薦的7×7;IDAN濾波的像素搜索窗口大小設(shè)置為15×15;Pretest NLM方法的影像塊和濾波窗口大小同本文方法設(shè)置一致。

    3.1 模擬試驗(yàn)

    模擬試驗(yàn)部分利用Monte Carlo方法[18]生成了一幅PolSAR影像。為了定量評(píng)價(jià)算法對(duì)噪聲的抑制及邊緣保持的效果,本文采用等效視數(shù)(equivalent number of looks,ENL)[24]和基于平均比的邊緣保持度(edge preservation degree based on ratio of the average,EPD-ROA)[25]這兩個(gè)指數(shù)。

    ENL指數(shù)的定義為

    ENL=1/σ2

    (16)

    式中,σ代表強(qiáng)度影像在同質(zhì)區(qū)的方差系數(shù)。ENL的估計(jì)方法一般都是人工選取影像內(nèi)幾塊大的同質(zhì)區(qū)域,利用以上公式求得每塊區(qū)域的等效視數(shù)后,取平均得到影像全局的等效視數(shù)。ENL越大,代表影像的噪聲水平越低。

    EPD-ROA的定義如下

    (17)

    式中,m代表影像的大小;uD1(i)和uD2(i)分別表示去噪所得影像的某個(gè)方向上兩相鄰像素的強(qiáng)度值;類似地,uO1(i)和uO2(i)分別表示原始影像的某個(gè)方向上兩相鄰像素的強(qiáng)度值。EPD-ROA指數(shù)越接近于1,意味著濾波算法有著越好的邊緣保持能力。

    (18)

    最后對(duì)所有地物所有像素的ARB指數(shù)求均值即為參量θ最終的ARB指數(shù)。試驗(yàn)部分所采用的極化參數(shù)為Cloude極化分解參數(shù)(極化散射熵H,平均alpha角α,極化反熵A)[27],原因是該分解模型能夠直觀地反映地物的物理散射特性。

    圖2展示了幾種去噪算法的處理結(jié)果,表1為其定量評(píng)價(jià)結(jié)果。增強(qiáng)Lee濾波在較大程度上抑制了影像的噪聲。然而,仔細(xì)觀察可以發(fā)現(xiàn),這一算法處理后影像呈現(xiàn)出一定的斑塊效應(yīng),直接導(dǎo)致點(diǎn)目標(biāo)及邊緣受到輕度的模糊。產(chǎn)生這一問題的原因在于,該算法在濾波的過程中采用了一種選擇濾波窗口方向的策略,而非通常采用的方形窗口。與增強(qiáng)Lee濾波相比,雖然表1所示的ENL顯示出IDAN方法對(duì)噪聲的抑制能力有所提升,但其對(duì)于影像細(xì)節(jié)信息的模糊和弱化現(xiàn)象較為明顯。很顯然,與前面兩種傳統(tǒng)的局部濾波器相比,Pretest NLM濾波方法在噪聲抑制方面有了質(zhì)的提升,并且實(shí)現(xiàn)了對(duì)邊緣信息的優(yōu)良保持。然而,可以較明顯觀察到的是,Pretest NLM處理后的影像的點(diǎn)目標(biāo)模糊嚴(yán)重。相比之下,本文提出的去噪方法雖然ENL指數(shù)略微低于Pretest NLM算法,但影像的邊緣和細(xì)節(jié)信息都得到了較好的保持。就影像極化散射信息的保持性能而言,Pretest NLM方法和本文方法要明顯優(yōu)于其他兩種方法,并且這兩種方法對(duì)不同極化參數(shù)的估計(jì)精度各有優(yōu)劣。

    圖2 模擬影像去噪結(jié)果Fig.2 Filtering results on the simulated image

    表1 模擬影像的濾波定量評(píng)價(jià)結(jié)果Tab.1 Quantitative assessment on the simulated image

    3.2 真實(shí)試驗(yàn)

    除了模擬的PolSAR數(shù)據(jù),兩幅真實(shí)的PolSAR影像也被用來驗(yàn)證去噪算法的優(yōu)劣。這兩幅影像分別為C波段和L波段的AirSAR傳感器所獲得的數(shù)據(jù),視數(shù)為4。

    圖3展示了幾種去噪算法對(duì)C波段AirSAR數(shù)據(jù)的去噪結(jié)果??梢钥吹?,這幅影像場(chǎng)景較為復(fù)雜,既包括以表面散射為主的水體,也包含以體散射為主的森林以及散射機(jī)制較為復(fù)雜的城區(qū)。對(duì)于這類的影像,一個(gè)優(yōu)良的去噪方法應(yīng)當(dāng)在削弱噪聲的同時(shí)保持影像大部分的細(xì)節(jié)信息。很顯然,與模擬試驗(yàn)結(jié)果一樣,增強(qiáng)Lee濾波處理后的影像在同質(zhì)區(qū)存在著較明顯的斑塊效應(yīng),部分邊緣和點(diǎn)目標(biāo)丟失;IDAN雖然更好地抑制了影像的噪聲水平,但是過平滑問題較為嚴(yán)重;Pretest NLM算法的噪聲抑制能力在4種算法當(dāng)中表現(xiàn)最佳,但也存在輕微的斑塊效應(yīng),導(dǎo)致處理后影像城區(qū)部分一些紋理(主要是細(xì)小的道路)以及部分強(qiáng)散射點(diǎn)目標(biāo)丟失。相反,本文算法雖然噪聲去除能力略微遜于Pretest NLM,但細(xì)節(jié)保持能力有了較大的提升。值得注意的是,同Pretest NLM處理的結(jié)果相比,影像中機(jī)場(chǎng)跑道處的線狀地物(圖3(f)箭頭所指)隨著迭代的進(jìn)行反而變得更為模糊。出現(xiàn)這一問題的原因是:目視來看,該線狀目標(biāo)極化散射特性同機(jī)場(chǎng)周邊植被較為相近,因此初次濾波中大量權(quán)重略高的植被像素會(huì)參與對(duì)線狀地物的去噪。然而,初次濾波后這一目標(biāo)和植被像素的可區(qū)分性不僅沒有加大,目視來看反而變得更加接近。在這種情況下,迭代的后果就是這些像素的權(quán)重反而會(huì)提高,模糊問題不斷加深。可以看出,對(duì)于影像中的一些零星的細(xì)節(jié)目標(biāo),如果本文提出的方法在初次去噪過程中不能較好地將它們與其他異質(zhì)地物“區(qū)分”開來,反而有可能會(huì)出現(xiàn)不斷過平滑的問題。因此,如何完善初始估計(jì)的精確性是本文方法未來值得改進(jìn)之處。

    圖3 C波段AirSAR影像去噪結(jié)果Fig.3 Filtering results on the C band AirSAR image

    圖4展示了幾種去噪算法對(duì)一景L波段PolSAR影像的去噪結(jié)果。與之前的試驗(yàn)結(jié)果相似,增強(qiáng)Lee濾波帶來的塊狀效應(yīng)制約了其濾波效果;IDAN能更好地抑制噪聲,但影像細(xì)節(jié)信息出現(xiàn)了嚴(yán)重的丟失;Pretest NLM算法不僅更好地抑制了噪聲,而且對(duì)于影像邊緣的保持有了一定的提升;相比之下,本文提出的算法雖然ENL指標(biāo)不如PolSAR NLM方法,但EPD-ROA指標(biāo)反映出本文算法在邊緣和細(xì)節(jié)保持方面相比于前者的優(yōu)勢(shì)性。顯然,本文提出的算法更好地實(shí)現(xiàn)了去噪和細(xì)節(jié)保持之間的平衡。表2為幾種算法對(duì)C波段和L波段影像濾波的評(píng)價(jià)結(jié)果。

    圖4 L波段AirSAR影像去噪結(jié)果Fig.4 Filtering results on the L band AirSAR image

    表2 真實(shí)影像的濾波定量評(píng)價(jià)結(jié)果Tab.2 Quantitative assessment on the real images

    4 結(jié) 論

    相干斑濾波是PolSAR數(shù)據(jù)預(yù)處理中重要的一環(huán)。傳統(tǒng)PolSAR非局部均值濾波算法在度量像素或影像塊相似度時(shí)基本利用到的都是含噪聲影像的統(tǒng)計(jì)特性,而沒有考慮到原始像素間的差異性,而且往往無法利用到全極化信息,因此難以獲得穩(wěn)健準(zhǔn)確的相似性測(cè)度。本文提出了一種迭代優(yōu)化的PolSAR非局部均值方法。該方法的創(chuàng)新之處在于:在每次迭代去噪過程中,同時(shí)考慮原始影像全極化噪聲統(tǒng)計(jì)特性和前一次迭代所得的相對(duì)低噪聲度影像的全極化信息,以完善影像塊間極化相似性的度量。模擬的PolSAR數(shù)據(jù)和不同波段的真實(shí)PolSAR影像的去噪試驗(yàn)都表明,本文提出的算法既能有效地抑制相干斑噪聲,又能較好地保持影像的邊緣、強(qiáng)散射點(diǎn)目標(biāo)、極化散射特性等細(xì)節(jié)信息。最后需要指出的是,相比較于傳統(tǒng)的局部濾波器,運(yùn)行效率相對(duì)較低是非局部濾波方法的一大弊端,特別是,本文方法還采取了迭代去噪的策略。因此,如何加速算法的運(yùn)行效率從而提升算法的實(shí)用性,是未來值得研究的方向之一。

    猜你喜歡
    協(xié)方差鄰域相似性
    一類上三角算子矩陣的相似性與酉相似性
    淺析當(dāng)代中西方繪畫的相似性
    稀疏圖平方圖的染色數(shù)上界
    基于鄰域競(jìng)賽的多目標(biāo)優(yōu)化算法
    關(guān)于-型鄰域空間
    低滲透黏土中氯離子彌散作用離心模擬相似性
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    基于時(shí)序擴(kuò)展的鄰域保持嵌入算法及其在故障檢測(cè)中的應(yīng)用
    縱向數(shù)據(jù)分析中使用滑動(dòng)平均Cholesky分解對(duì)回歸均值和協(xié)方差矩陣進(jìn)行同時(shí)半?yún)?shù)建模
    男人舔奶头视频| 老司机深夜福利视频在线观看| 99国产精品99久久久久| 国产极品精品免费视频能看的| 美女大奶头视频| 亚洲18禁久久av| 成人午夜高清在线视频| 亚洲成人久久爱视频| 高清在线国产一区| 亚洲欧美精品综合一区二区三区| 国产伦精品一区二区三区视频9 | 一区福利在线观看| 精品乱码久久久久久99久播| 亚洲精华国产精华精| 婷婷丁香在线五月| 久久久久国产精品人妻aⅴ院| 久久人人精品亚洲av| 琪琪午夜伦伦电影理论片6080| 香蕉国产在线看| 99久久成人亚洲精品观看| 少妇熟女aⅴ在线视频| 国产三级中文精品| 国内毛片毛片毛片毛片毛片| e午夜精品久久久久久久| 人妻丰满熟妇av一区二区三区| 网址你懂的国产日韩在线| 美女大奶头视频| tocl精华| 亚洲精品乱码久久久v下载方式 | 婷婷六月久久综合丁香| 欧美成人一区二区免费高清观看 | 亚洲精品美女久久久久99蜜臀| 精品久久久久久久末码| 久久伊人香网站| 欧美不卡视频在线免费观看| 一进一出好大好爽视频| 男人和女人高潮做爰伦理| 级片在线观看| 国产99白浆流出| 免费电影在线观看免费观看| 禁无遮挡网站| 国产欧美日韩精品一区二区| 欧美一区二区精品小视频在线| 五月玫瑰六月丁香| 成人性生交大片免费视频hd| 欧美高清成人免费视频www| 国产三级中文精品| 日本免费a在线| 桃色一区二区三区在线观看| 天天一区二区日本电影三级| 国产主播在线观看一区二区| 桃红色精品国产亚洲av| 怎么达到女性高潮| 午夜福利高清视频| 国产精品精品国产色婷婷| 国产麻豆成人av免费视频| 国产精品电影一区二区三区| 久久久久亚洲av毛片大全| 90打野战视频偷拍视频| 首页视频小说图片口味搜索| 欧美激情久久久久久爽电影| 精品国产美女av久久久久小说| 久久天躁狠狠躁夜夜2o2o| 美女扒开内裤让男人捅视频| 亚洲国产高清在线一区二区三| 黄片小视频在线播放| 国产成人影院久久av| 丰满人妻熟妇乱又伦精品不卡| 黄色成人免费大全| 搡老岳熟女国产| 亚洲熟妇中文字幕五十中出| 国产精品女同一区二区软件 | 国产极品精品免费视频能看的| 两性午夜刺激爽爽歪歪视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 成年版毛片免费区| 国产高潮美女av| 久久国产乱子伦精品免费另类| 亚洲片人在线观看| 变态另类丝袜制服| 俄罗斯特黄特色一大片| 国产精品 国内视频| 欧美色视频一区免费| 成人鲁丝片一二三区免费| av黄色大香蕉| 一进一出抽搐gif免费好疼| 99在线人妻在线中文字幕| 国产欧美日韩一区二区精品| 国产成+人综合+亚洲专区| 村上凉子中文字幕在线| 欧美黑人巨大hd| 日本在线视频免费播放| 夜夜看夜夜爽夜夜摸| 精品乱码久久久久久99久播| 非洲黑人性xxxx精品又粗又长| 色综合婷婷激情| 亚洲美女黄片视频| 国产精品久久视频播放| 精品国产乱子伦一区二区三区| 国产精品综合久久久久久久免费| 精品国产乱码久久久久久男人| 精品久久蜜臀av无| 午夜两性在线视频| 国产黄a三级三级三级人| 国产成人av教育| 久久精品国产综合久久久| 午夜影院日韩av| 一夜夜www| 99热这里只有精品一区 | 免费大片18禁| 亚洲精品乱码久久久v下载方式 | a级毛片a级免费在线| aaaaa片日本免费| 国产精品亚洲美女久久久| 97人妻精品一区二区三区麻豆| 最近最新中文字幕大全免费视频| 日韩大尺度精品在线看网址| 亚洲国产精品sss在线观看| 最近最新中文字幕大全免费视频| 日韩欧美国产在线观看| 免费在线观看影片大全网站| 日本免费a在线| 99热这里只有精品一区 | 夜夜爽天天搞| 日日夜夜操网爽| 人人妻人人澡欧美一区二区| 中国美女看黄片| 国内精品久久久久精免费| 身体一侧抽搐| 久久天堂一区二区三区四区| 在线观看日韩欧美| 老司机深夜福利视频在线观看| 久久精品aⅴ一区二区三区四区| 午夜亚洲福利在线播放| 老司机午夜十八禁免费视频| 99国产精品一区二区三区| 精品国产超薄肉色丝袜足j| 色吧在线观看| 无人区码免费观看不卡| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| 99久久精品国产亚洲精品| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 桃色一区二区三区在线观看| 夜夜看夜夜爽夜夜摸| 琪琪午夜伦伦电影理论片6080| 欧美丝袜亚洲另类 | 美女cb高潮喷水在线观看 | av视频在线观看入口| 丰满人妻熟妇乱又伦精品不卡| 身体一侧抽搐| 国产成人精品久久二区二区91| 在线播放国产精品三级| 又粗又爽又猛毛片免费看| 久久中文字幕一级| 观看免费一级毛片| 亚洲av成人不卡在线观看播放网| 欧美成人性av电影在线观看| 精品国产乱子伦一区二区三区| 午夜免费成人在线视频| 成人无遮挡网站| 亚洲 欧美一区二区三区| 人人妻,人人澡人人爽秒播| 日日干狠狠操夜夜爽| www.熟女人妻精品国产| 久久久久久九九精品二区国产| 欧美性猛交╳xxx乱大交人| 成熟少妇高潮喷水视频| 99久久久亚洲精品蜜臀av| av在线天堂中文字幕| 亚洲avbb在线观看| 日韩精品青青久久久久久| 69av精品久久久久久| 男女那种视频在线观看| 99国产极品粉嫩在线观看| 91久久精品国产一区二区成人 | 久久久久久久久中文| 国产不卡一卡二| 毛片女人毛片| 国产精品亚洲美女久久久| 国产精品久久久人人做人人爽| 在线观看免费午夜福利视频| 亚洲欧美精品综合一区二区三区| 久久中文字幕一级| 脱女人内裤的视频| 久久精品aⅴ一区二区三区四区| 别揉我奶头~嗯~啊~动态视频| 国产免费男女视频| АⅤ资源中文在线天堂| 美女扒开内裤让男人捅视频| 日本五十路高清| 男女视频在线观看网站免费| 亚洲av成人一区二区三| 亚洲精品色激情综合| 日本撒尿小便嘘嘘汇集6| 免费无遮挡裸体视频| 天堂√8在线中文| 亚洲欧美精品综合久久99| 国产一级毛片七仙女欲春2| 听说在线观看完整版免费高清| 国产精品永久免费网站| 日韩欧美国产一区二区入口| 欧美乱妇无乱码| 亚洲国产日韩欧美精品在线观看 | 丝袜人妻中文字幕| 欧美绝顶高潮抽搐喷水| 国产精品 国内视频| 日日夜夜操网爽| 国产又色又爽无遮挡免费看| 精品一区二区三区四区五区乱码| 国产亚洲av高清不卡| 19禁男女啪啪无遮挡网站| 精品久久久久久成人av| 欧美日本亚洲视频在线播放| 亚洲成av人片免费观看| 色老头精品视频在线观看| 美女高潮喷水抽搐中文字幕| 成人亚洲精品av一区二区| 欧美一区二区精品小视频在线| 亚洲av熟女| 国产探花在线观看一区二区| 国产精品久久久人人做人人爽| 亚洲人成伊人成综合网2020| 午夜福利在线观看免费完整高清在 | 国产成人精品久久二区二区免费| 久久精品国产亚洲av香蕉五月| 别揉我奶头~嗯~啊~动态视频| 国产精品国产高清国产av| 国产欧美日韩精品亚洲av| 成人三级做爰电影| av天堂在线播放| 一边摸一边抽搐一进一小说| 757午夜福利合集在线观看| 床上黄色一级片| 午夜成年电影在线免费观看| 国产精品久久久久久亚洲av鲁大| 精品一区二区三区视频在线 | 欧美在线黄色| 精品熟女少妇八av免费久了| 日韩精品中文字幕看吧| 黄色视频,在线免费观看| 国产人伦9x9x在线观看| 啦啦啦免费观看视频1| 99精品久久久久人妻精品| 曰老女人黄片| 欧美日韩乱码在线| 小蜜桃在线观看免费完整版高清| 色吧在线观看| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人澡欧美一区二区| 黄色片一级片一级黄色片| 一个人看的www免费观看视频| 久久精品综合一区二区三区| 看免费av毛片| 免费电影在线观看免费观看| 日本三级黄在线观看| 久久久久亚洲av毛片大全| a在线观看视频网站| 免费在线观看视频国产中文字幕亚洲| 亚洲国产欧洲综合997久久,| 最近最新中文字幕大全电影3| 日本与韩国留学比较| 后天国语完整版免费观看| 亚洲午夜精品一区,二区,三区| 综合色av麻豆| 国产探花在线观看一区二区| 亚洲精品色激情综合| 一卡2卡三卡四卡精品乱码亚洲| 亚洲专区中文字幕在线| 欧美色视频一区免费| h日本视频在线播放| 久久久精品欧美日韩精品| 91老司机精品| 久久中文看片网| xxx96com| 亚洲 国产 在线| 成在线人永久免费视频| 成人欧美大片| 男人舔女人下体高潮全视频| 久久性视频一级片| 三级男女做爰猛烈吃奶摸视频| 叶爱在线成人免费视频播放| 99国产精品99久久久久| 亚洲精品粉嫩美女一区| 国产探花在线观看一区二区| 国产精品一区二区三区四区久久| 久久午夜综合久久蜜桃| 少妇丰满av| 欧美色欧美亚洲另类二区| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 激情在线观看视频在线高清| 精品免费久久久久久久清纯| 又黄又粗又硬又大视频| 久久中文字幕人妻熟女| 看黄色毛片网站| 国产成年人精品一区二区| 美女扒开内裤让男人捅视频| 日韩 欧美 亚洲 中文字幕| 99热6这里只有精品| 麻豆久久精品国产亚洲av| 国产精品 国内视频| 不卡av一区二区三区| 国产一区二区三区在线臀色熟女| 免费观看精品视频网站| 日韩 欧美 亚洲 中文字幕| 桃红色精品国产亚洲av| 国产成人精品久久二区二区免费| 床上黄色一级片| xxx96com| 深夜精品福利| 欧美大码av| 午夜a级毛片| 国产成人欧美在线观看| 亚洲av成人不卡在线观看播放网| 亚洲人成电影免费在线| 99久久久亚洲精品蜜臀av| 男人的好看免费观看在线视频| 又黄又爽又免费观看的视频| 18禁美女被吸乳视频| 99久久99久久久精品蜜桃| 免费高清视频大片| 久久人妻av系列| 色综合亚洲欧美另类图片| 欧美另类亚洲清纯唯美| 亚洲自拍偷在线| 亚洲欧美激情综合另类| 午夜福利在线在线| 免费在线观看亚洲国产| 欧美又色又爽又黄视频| av女优亚洲男人天堂 | 亚洲欧美激情综合另类| av中文乱码字幕在线| 亚洲人与动物交配视频| 国产乱人视频| 午夜激情福利司机影院| 成人国产一区最新在线观看| 国产麻豆成人av免费视频| 国产高潮美女av| 亚洲欧美日韩高清专用| 99久久久亚洲精品蜜臀av| 久久九九热精品免费| 全区人妻精品视频| 毛片女人毛片| 两个人的视频大全免费| 精品熟女少妇八av免费久了| 嫩草影视91久久| 在线观看免费午夜福利视频| 嫩草影视91久久| 一进一出抽搐动态| av天堂在线播放| x7x7x7水蜜桃| av中文乱码字幕在线| 亚洲欧美激情综合另类| 亚洲自拍偷在线| 美女高潮的动态| 亚洲乱码一区二区免费版| 欧美色欧美亚洲另类二区| 神马国产精品三级电影在线观看| 亚洲成人久久爱视频| 欧美一区二区精品小视频在线| 国产淫片久久久久久久久 | 一级毛片精品| 国产成人福利小说| 色综合婷婷激情| 午夜福利免费观看在线| 国产亚洲精品综合一区在线观看| 偷拍熟女少妇极品色| netflix在线观看网站| 久久久久久大精品| 91老司机精品| 午夜福利欧美成人| 美女扒开内裤让男人捅视频| 国产精品久久久人人做人人爽| 色尼玛亚洲综合影院| 国产亚洲精品一区二区www| 成人特级av手机在线观看| 免费一级毛片在线播放高清视频| 综合色av麻豆| 国产免费男女视频| 亚洲国产精品成人综合色| 超碰成人久久| 欧美激情在线99| 白带黄色成豆腐渣| 十八禁人妻一区二区| 久久精品91蜜桃| 午夜精品在线福利| 999久久久精品免费观看国产| 欧美乱妇无乱码| 999精品在线视频| tocl精华| 男人舔奶头视频| xxxwww97欧美| 国产一区在线观看成人免费| 国产av一区在线观看免费| 国产91精品成人一区二区三区| 欧美大码av| 国产成年人精品一区二区| 99久久成人亚洲精品观看| 色吧在线观看| 成人一区二区视频在线观看| 亚洲av日韩精品久久久久久密| 成人欧美大片| 午夜福利免费观看在线| 熟妇人妻久久中文字幕3abv| 男女视频在线观看网站免费| 国产精品久久久久久久电影 | 精品国内亚洲2022精品成人| 亚洲av日韩精品久久久久久密| 国产精品电影一区二区三区| a级毛片a级免费在线| 国产午夜精品论理片| 91在线精品国自产拍蜜月 | 色综合亚洲欧美另类图片| 国产综合懂色| 亚洲中文av在线| 亚洲成av人片在线播放无| 叶爱在线成人免费视频播放| 老熟妇乱子伦视频在线观看| 青草久久国产| 男人舔奶头视频| 美女 人体艺术 gogo| 国内精品美女久久久久久| 操出白浆在线播放| 午夜影院日韩av| 午夜福利成人在线免费观看| 国产蜜桃级精品一区二区三区| 色av中文字幕| 黄片小视频在线播放| 啦啦啦观看免费观看视频高清| 国产一区二区在线av高清观看| 亚洲成av人片在线播放无| 精品免费久久久久久久清纯| 欧美日韩综合久久久久久 | 国产午夜精品论理片| 国产高清有码在线观看视频| 久久香蕉精品热| av女优亚洲男人天堂 | 日本在线视频免费播放| 亚洲国产欧洲综合997久久,| 久99久视频精品免费| 免费无遮挡裸体视频| 免费人成视频x8x8入口观看| 女同久久另类99精品国产91| 久9热在线精品视频| 欧美日韩黄片免| 午夜精品在线福利| 亚洲 国产 在线| 久久久国产成人精品二区| 男女下面进入的视频免费午夜| 免费看日本二区| av欧美777| 亚洲天堂国产精品一区在线| 国产精品九九99| 国产亚洲精品av在线| 国产高清有码在线观看视频| 欧美三级亚洲精品| 非洲黑人性xxxx精品又粗又长| 丁香欧美五月| 两个人看的免费小视频| 欧美一区二区国产精品久久精品| 成人高潮视频无遮挡免费网站| 亚洲精品久久国产高清桃花| 老熟妇仑乱视频hdxx| 在线看三级毛片| 啦啦啦免费观看视频1| 中文字幕最新亚洲高清| 久久久久久国产a免费观看| 国产免费男女视频| 午夜免费激情av| 午夜视频精品福利| 亚洲第一欧美日韩一区二区三区| 在线永久观看黄色视频| 成人欧美大片| 最近最新中文字幕大全免费视频| 亚洲精品在线美女| e午夜精品久久久久久久| 久久亚洲精品不卡| 日韩大尺度精品在线看网址| 欧美日韩乱码在线| 亚洲国产高清在线一区二区三| 午夜激情欧美在线| 久久亚洲精品不卡| 成人欧美大片| 亚洲欧美日韩卡通动漫| 人妻夜夜爽99麻豆av| 久久这里只有精品19| 少妇丰满av| 夜夜躁狠狠躁天天躁| 欧美性猛交╳xxx乱大交人| 国产精品99久久99久久久不卡| 好男人在线观看高清免费视频| 国产野战对白在线观看| 中文字幕最新亚洲高清| 中文字幕熟女人妻在线| 变态另类成人亚洲欧美熟女| 久9热在线精品视频| 少妇人妻一区二区三区视频| 一本精品99久久精品77| 舔av片在线| 一二三四在线观看免费中文在| 男人和女人高潮做爰伦理| 午夜福利18| 国产精品永久免费网站| 午夜福利免费观看在线| 婷婷精品国产亚洲av| 1024手机看黄色片| 国产三级在线视频| 国产久久久一区二区三区| 麻豆国产97在线/欧美| 亚洲最大成人中文| 国产欧美日韩一区二区精品| 97人妻精品一区二区三区麻豆| 久久国产精品人妻蜜桃| 嫩草影院精品99| 久久国产精品影院| 老司机午夜十八禁免费视频| 丰满的人妻完整版| 99久久99久久久精品蜜桃| 精品久久久久久久末码| 波多野结衣高清无吗| 亚洲国产中文字幕在线视频| 亚洲精品美女久久av网站| 又黄又粗又硬又大视频| 欧美成狂野欧美在线观看| 亚洲真实伦在线观看| 天堂影院成人在线观看| 国产亚洲欧美在线一区二区| 国产精品99久久99久久久不卡| 女警被强在线播放| 男人舔女人下体高潮全视频| 国产av麻豆久久久久久久| 日韩av在线大香蕉| 亚洲成人久久性| 91麻豆av在线| 久久精品影院6| 亚洲国产精品合色在线| 精品熟女少妇八av免费久了| av欧美777| 麻豆一二三区av精品| 99国产精品一区二区蜜桃av| 欧美一区二区精品小视频在线| 熟女少妇亚洲综合色aaa.| 长腿黑丝高跟| 手机成人av网站| 亚洲人与动物交配视频| 久久99热这里只有精品18| 真人做人爱边吃奶动态| 亚洲精华国产精华精| 欧美又色又爽又黄视频| 久久午夜综合久久蜜桃| 在线观看美女被高潮喷水网站 | 观看免费一级毛片| 免费观看人在逋| 日本黄色视频三级网站网址| 成年免费大片在线观看| 欧美日韩一级在线毛片| 男女做爰动态图高潮gif福利片| 看黄色毛片网站| 性色avwww在线观看| 亚洲精品一区av在线观看| 精品久久久久久久久久久久久| 久久久久久大精品| 九九在线视频观看精品| 国产单亲对白刺激| 国产高清视频在线观看网站| 岛国在线观看网站| 99久久国产精品久久久| 九九久久精品国产亚洲av麻豆 | 少妇人妻一区二区三区视频| 蜜桃久久精品国产亚洲av| 免费电影在线观看免费观看| 啦啦啦韩国在线观看视频| 噜噜噜噜噜久久久久久91| 国产精品亚洲美女久久久| 搡老妇女老女人老熟妇| 欧美激情久久久久久爽电影| 中国美女看黄片| 午夜福利在线观看免费完整高清在 | 亚洲 国产 在线| 高清毛片免费观看视频网站| 国产精品久久久久久久电影 | 男插女下体视频免费在线播放| 999精品在线视频| 在线免费观看不下载黄p国产 | 天天添夜夜摸| 岛国在线观看网站| 18禁裸乳无遮挡免费网站照片| 老汉色∧v一级毛片| aaaaa片日本免费| 国产精品香港三级国产av潘金莲| 在线播放国产精品三级| 久久精品影院6| 国产精品亚洲美女久久久| 亚洲国产欧美一区二区综合| 夜夜看夜夜爽夜夜摸| 性欧美人与动物交配| 国产av不卡久久| 老司机在亚洲福利影院| e午夜精品久久久久久久| 长腿黑丝高跟| 亚洲专区国产一区二区| 欧美中文日本在线观看视频| 女生性感内裤真人,穿戴方法视频| 女同久久另类99精品国产91| 精品久久蜜臀av无| 成人特级av手机在线观看| 日本黄色视频三级网站网址| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕熟女人妻在线| 亚洲一区二区三区色噜噜| 美女扒开内裤让男人捅视频| 极品教师在线免费播放| 国产高清视频在线播放一区|