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

    基于噪聲可見(jiàn)性函數(shù)的SAR圖像增強(qiáng)快速算法

    2018-03-08 08:52:25朱逸飛
    關(guān)鍵詞:細(xì)節(jié)

    朱逸飛,楊 國(guó),王 強(qiáng),吳 文

    (南京理工大學(xué) 近程高速目標(biāo)探測(cè)技術(shù)國(guó)防重點(diǎn)學(xué)科實(shí)驗(yàn)室,南京 210094)

    0 引言

    合成孔徑雷達(dá)SAR(synthetic aperture radar)是一種高分辨率成像雷達(dá),具有全天時(shí)、全天候、穿透性強(qiáng)等特點(diǎn),廣泛應(yīng)用于軍事和民用領(lǐng)域。但由于戰(zhàn)場(chǎng)環(huán)境復(fù)雜、大氣衰減、硬件設(shè)備干擾等原因,SAR成像結(jié)果會(huì)存在整體灰度分布低且集中、圖像對(duì)比度低、含有噪聲等問(wèn)題,造成目標(biāo)不清晰、圖像模糊,影響圖像識(shí)別效果,因此需要對(duì)SAR圖像進(jìn)行增強(qiáng)[1]。目前國(guó)內(nèi)外利用數(shù)據(jù)處理的方法提高SAR圖像質(zhì)量主要分為兩類[2]:一是在成像階段利用系統(tǒng)回波數(shù)據(jù),通過(guò)超分辨成像算法獲得高分辨率的SAR圖像,即SAR的超分辨成像方法;二是在成像后的圖像數(shù)據(jù)基礎(chǔ)上進(jìn)行圖像質(zhì)量提高的處理。其中后者是本文的研究?jī)?nèi)容,典型算法如Kuan算法、Lee算法、雙邊濾波算法、偏微分算法等,這些算法能夠?yàn)V除部分噪聲并保留細(xì)節(jié)邊緣和目標(biāo)特征[3]。但是,典型的SAR圖像增強(qiáng)算法復(fù)雜度較高,面對(duì)高分辨率圖像時(shí),需要的計(jì)算時(shí)間較長(zhǎng),在目標(biāo)識(shí)別、跟蹤及災(zāi)難評(píng)估等對(duì)實(shí)時(shí)性要求較高的領(lǐng)域,無(wú)法滿足系統(tǒng)要求。

    與中央處理器CPU(central processing unit)不同,圖形處理器GPU(graphic processing unit)專為密集型、高度并行化的計(jì)算而設(shè)計(jì),其中用于數(shù)據(jù)處理的晶體管數(shù)量遠(yuǎn)遠(yuǎn)大于緩存和邏輯控制部分,這使得GPU更適合處理無(wú)邏輯關(guān)系數(shù)據(jù)的并行運(yùn)算[4]。2007年,NVIDIA公司推出通用并行計(jì)算架構(gòu)CUDA(compute unified device architecture),目的是將GPU作為并行計(jì)算設(shè)備,進(jìn)行通用并行計(jì)算[5]。CUDA為開(kāi)發(fā)者提供了硬件的直接訪問(wèn)接口,可以方便地寫(xiě)出在GPU上執(zhí)行的程序,而不用像過(guò)去的GPGPU架構(gòu),如OpenGL等,將計(jì)算映射到圖形API中,大大降低了開(kāi)發(fā)門(mén)檻和難度。自推出后,CUDA被廣泛應(yīng)用于計(jì)算機(jī)可視化、音視頻編解碼、流體力學(xué)模擬等領(lǐng)域,都能獲得較高的加速比,有著良好的應(yīng)用前景[6]。

    本文提出一種基于噪聲可見(jiàn)性函數(shù)的SAR圖像增強(qiáng)快速算法及其基于GPU的并行化實(shí)現(xiàn)方法。該算法運(yùn)用圖像分層理論,分離SAR圖像細(xì)節(jié)與噪聲,在對(duì)圖像細(xì)節(jié)層進(jìn)行處理時(shí),結(jié)合人眼視覺(jué)特性,引入噪聲可見(jiàn)性函數(shù),控制細(xì)節(jié)層增益系數(shù),實(shí)現(xiàn)圖像增強(qiáng)。在提高圖像質(zhì)量、增強(qiáng)圖像細(xì)節(jié)的同時(shí),充分運(yùn)用GPU的并行計(jì)算特點(diǎn),具有較高的實(shí)時(shí)性。

    1 基于圖像分層的增強(qiáng)算法

    為了抑制圖像噪聲、增強(qiáng)圖像細(xì)節(jié)、提高圖像識(shí)別效果,Branchitta F等人提出了圖像分層理論[7]。其基本思想為:首先,選擇一種分層濾波器對(duì)原始輸入圖像進(jìn)行處理,抹平圖像細(xì)節(jié)和微小波動(dòng),得到圖像基本層和細(xì)節(jié)層;其次,分別對(duì)圖像基本層和細(xì)節(jié)層進(jìn)行處理,因?yàn)榛緦又胁话瑘D像細(xì)節(jié)信息,所以處理過(guò)程中不必考慮細(xì)節(jié)損失問(wèn)題,重點(diǎn)在于圖像對(duì)比度的調(diào)整;由于細(xì)節(jié)層不僅包含原始圖像的細(xì)節(jié)信息,還有噪聲的存在,因此需要對(duì)細(xì)節(jié)層進(jìn)行噪聲抑制和細(xì)節(jié)增強(qiáng);最后,將處理后的基本層和細(xì)節(jié)層融合并調(diào)整顯示范圍,得到輸出圖像。

    研究基于圖像分層的增強(qiáng)算法,其重點(diǎn)在于分層濾波器、基本層圖像處理方法和細(xì)節(jié)層圖像增強(qiáng)方法的選擇。本文采用雙邊濾波作為分層濾波器,對(duì)SAR圖像進(jìn)行雙邊濾波,得到抹平圖像細(xì)節(jié)和微小波動(dòng)后的基本層,并與原圖像相減,得到包含細(xì)節(jié)信息和大量噪聲的細(xì)節(jié)層;對(duì)基本層采用γ<1的伽馬變換;通過(guò)噪聲可見(jiàn)性函數(shù)控制細(xì)節(jié)層增益系數(shù)。算法框架如圖1所示。

    圖1 基于圖像分層的增強(qiáng)算法框架

    雙邊濾波(Bilateralfilter)結(jié)合空域信息和灰度相似性,對(duì)圖像進(jìn)行基于空間分布的高斯濾波,是一種非線性的濾波方法。鄰域內(nèi)離中心像素較遠(yuǎn)的像素占有的權(quán)重較小,不會(huì)對(duì)中心像素的值造成太大影響,這樣有利于圖像細(xì)節(jié)邊緣像素值的保存,達(dá)到保邊去噪的目的。

    設(shè)f(x,y)為原圖像,g(x,y)為原圖像經(jīng)雙邊濾波后的結(jié)果,即基本層,其表達(dá)式為:

    (1)

    式中,Sx,y表示以像素(x,y)為中心,半寬為N的鄰域,ω(i,j)為兩部分系數(shù)的乘積:

    ω(i,j)=ωs(i,j)ωr(i,j)

    (2)

    式中,ωs(i,j)為空域?yàn)V波器系數(shù),ωr(i,j)為值域?yàn)V波器系數(shù):

    (3)

    (4)

    對(duì)基本層g(x,y)采用γ<1的伽馬變換,拉伸圖像灰度范圍,提高對(duì)比度。伽馬變換,也叫做冪律變換,是一種非線性的灰度拉伸變換。設(shè)變換后的基本層為gγ(x,y),變換形式為:

    gγ(x,y)=C*g(x,y)γ

    (5)

    式中,γ為變換系數(shù),C為控制像素變化范圍的常數(shù)。

    2 基于NVF的細(xì)節(jié)層增強(qiáng)

    在研究水印估計(jì)問(wèn)題時(shí),S.Voloshnovskiy等人提出了噪聲可見(jiàn)性函數(shù)NVF(noise visibility function)[8-9]:

    (6)

    歸一化處理后,式(6)變?yōu)椋?/p>

    (7)

    式中,θ為常數(shù),m(x,y)為噪聲掩膜函數(shù),與圖像的局部方差有關(guān),可用于衡量圖像細(xì)節(jié)。

    由式(7)可以看出,NVF(x,y)可以表示圖像中各像素點(diǎn)對(duì)噪聲的敏感程度,與圖像局部能量成反比。越是平坦的區(qū)域,m(x,y)越小,NVF(x,y)越接近1,圖像在該像素處較為敏感,噪聲對(duì)人眼刺激較大,視覺(jué)效果降低;越是變化劇烈的區(qū)域,m(x,y)越大,NVF(x,y)越接近0,圖像在該像素處允許較大的噪聲,其對(duì)圖像視覺(jué)效果的影響不明顯,噪聲表現(xiàn)出低可見(jiàn)性,不易被察覺(jué)。

    將原圖像f(x,y)與基本層g(x,y)相減得到細(xì)節(jié)層圖像h(x,y)。h(x,y)中包含了原圖像中的細(xì)節(jié)信息,必須進(jìn)行保留或增強(qiáng),但細(xì)節(jié)層中又不可避免地存在原圖像中的大部分噪聲。因此,在對(duì)細(xì)節(jié)層進(jìn)行處理時(shí),不能直接簡(jiǎn)單地放大,這會(huì)造成平坦區(qū)域噪聲的過(guò)度放大。

    增強(qiáng)細(xì)節(jié)層圖像時(shí),通過(guò)NVF區(qū)分圖像的平坦與變化劇烈的區(qū)域,對(duì)不同的區(qū)域采用不同的增益系數(shù)。設(shè)細(xì)節(jié)層像素點(diǎn)(x,y)的增益系數(shù)為G(x,y),對(duì)于平坦區(qū)域的像素點(diǎn),NVF趨近于1,G(x,y)需要取較小值,避免噪聲被增強(qiáng);對(duì)于變化劇烈的區(qū)域,NVF趨近于0,G(x,y)可以取較大值,在增強(qiáng)細(xì)節(jié)時(shí)不用考慮噪聲的問(wèn)題。設(shè)Gmax為G(x,y)的最大值,Gmin為G(x,y)的最小值,構(gòu)造增益系數(shù)公式如下:

    G(x,y)=Gmin+(1-NVF(x,y))(Gmax-Gmin)

    (8)

    (9)

    且有:

    m(x,y)=k(x,y)-1-1

    (10)

    因?yàn)閗(x,y)∈[0,1],無(wú)需調(diào)整范圍,所以取θ=1。由式(7)、式(8)和式(10)可得:

    G(x,y)=Gmin+(1-k(x,y))(Gmax-Gmin)

    (11)

    增強(qiáng)后的細(xì)節(jié)層hNVF(x,y)可表示為:

    hNVF(x,y)=h(x,y)·G(x,y)

    (12)

    3 基于CUDA的圖像增強(qiáng)算法實(shí)現(xiàn)

    3.1 CUDA編程模型

    GPU擁有強(qiáng)大的并行計(jì)算能力,但早期的GPU編程存在很多困難,如圖形API、交互接口復(fù)雜等,這些都限制了GPU在通用計(jì)算上的發(fā)展[10]。CUDA的出現(xiàn)解決了這一問(wèn)題,其包含一整套軟硬件體系,支持通用并行計(jì)算。由于CUDA采用類C語(yǔ)言進(jìn)行開(kāi)發(fā),具有高性能計(jì)算指令開(kāi)發(fā)能力,開(kāi)發(fā)者能夠在GPU強(qiáng)大的并行計(jì)算能力基礎(chǔ)上建立更高效率的密集數(shù)據(jù)計(jì)算方案。

    一個(gè)完整的CUDA程序包含在CPU(主機(jī)端,host)中執(zhí)行的部分和在GPU(設(shè)備端,device)中執(zhí)行的部分。如圖2所示,主機(jī)端負(fù)責(zé)數(shù)據(jù)的輸入輸出、內(nèi)存分配、主機(jī)與設(shè)備間的數(shù)據(jù)傳輸;而在設(shè)備端中運(yùn)行的核函數(shù)(kernel)承擔(dān)著并行計(jì)算的任務(wù)[11]。由于CPU對(duì)顯存中的數(shù)據(jù)進(jìn)行存取時(shí)只能通過(guò)PCIExpress接口,其理論帶寬為雙向各4GB/s,速度較慢,若頻繁地進(jìn)行存取顯存的操作會(huì)降低GPU計(jì)算效率[12]。

    圖2 CUDA編程模型

    在并行線程總數(shù)確定的情況下,為了有效利用GPU的計(jì)算資源,需要對(duì)核函數(shù)進(jìn)行合理的線程配置,根據(jù)GPU規(guī)格和需要計(jì)算的數(shù)據(jù)規(guī)模,合理地設(shè)定線程塊Block的數(shù)量和每個(gè)線程塊中線程Thread的數(shù)量。每個(gè)線程塊中線程數(shù)量應(yīng)為該GPU的warp大小的整數(shù)倍[13]。warp為最小單指令多線程執(zhí)行單元,以32個(gè)線程為一組,并行創(chuàng)建、調(diào)度和執(zhí)行[14],考慮存儲(chǔ)器大小的限制,本文設(shè)定每個(gè)線程塊中線程數(shù)量為256,并根據(jù)不同的圖像大小劃分線程塊數(shù)量。

    3.2 并行化圖像增強(qiáng)算法實(shí)現(xiàn)

    上文中圖像增強(qiáng)算法復(fù)雜度較高,對(duì)于高分辨率圖像而言,串行實(shí)現(xiàn)方法實(shí)時(shí)性較低。該算法中像素點(diǎn)間的運(yùn)算相對(duì)獨(dú)立,數(shù)據(jù)依賴程度較低,具有較高的內(nèi)在并行性,適合GPU并行實(shí)現(xiàn)。在并行化雙邊濾波算法的基礎(chǔ)上[15]進(jìn)行線程內(nèi)的分層處理、基本層處理和基于NVF的細(xì)節(jié)層增強(qiáng),并行完成各像素的全部計(jì)算過(guò)程,實(shí)現(xiàn)算法并行化設(shè)計(jì)。

    具體流程如下:

    1)獲取圖像數(shù)據(jù),初始化變量,申請(qǐng)GPU顯存空間:cudaMalloc((void**)&Rdata,(width*height*sizeof(double)));

    2)完成CPU與GPU間的數(shù)據(jù)通信;

    3)定義線程配置:

    Dim3grid(blockx,blocky,1);

    Dim3block(threadx,thready,1);

    4)啟動(dòng)kernel函數(shù),進(jìn)行并行化雙邊濾波和分層,根據(jù)式(11)計(jì)算細(xì)節(jié)層增益系數(shù),完成細(xì)節(jié)層和基本層處理以及圖像融合;

    5)實(shí)現(xiàn)圖像輸出,釋放內(nèi)存和顯存空間。

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

    本文將從算法效果和加速性能兩方面進(jìn)行實(shí)驗(yàn)。硬件平臺(tái):CPU為IntelCorei5 6400,4核,主頻2.70GHz;GPU為NVIDIAGeForceGTX960,1024核,運(yùn)行頻率1.22GHz;編譯環(huán)境為VisualStudio2012 和CUDA7.5。

    4.1 算法效果

    本次實(shí)驗(yàn)中,考慮視覺(jué)效果和計(jì)算復(fù)雜度,取雙邊濾波參數(shù)σs=40、σr=20,鄰域半寬N=5,基本層伽馬變換系數(shù)γ=0.5,細(xì)節(jié)層增益系數(shù)最大值Gmax=1.5,最小值Gmin=1。算法實(shí)驗(yàn)效果如圖3所示。

    圖3 算法效果

    圖3(a)為1024×1024的原始SAR圖像,對(duì)圖像添加噪聲以模擬SAR成像結(jié)果,結(jié)果如圖3(b)所示。圖3(c)為經(jīng)雙邊濾波分層后得到的未處理的基本層,圖像的細(xì)節(jié)和微小波動(dòng)被抹平;圖3(d)為經(jīng)過(guò)伽馬變換后的基本層,圖像對(duì)比度得到增強(qiáng);圖3(e)為未處理的細(xì)節(jié)層,圖像細(xì)節(jié)和噪聲同時(shí)存在;圖3(f)為經(jīng)過(guò)增強(qiáng)后的細(xì)節(jié)層,可以看出,在細(xì)節(jié)增強(qiáng)的同時(shí),部分噪聲得到了抑制;圖3(g)和圖3(h)分別為本文所提算法的串行實(shí)現(xiàn)計(jì)算結(jié)果和并行化實(shí)現(xiàn)方法的計(jì)算結(jié)果,可以看到,圖像對(duì)比度提高的同時(shí),噪聲得到了抑制,農(nóng)田部分的紋理變得更加清晰,房屋、河流等輪廓也變得十分明顯。

    本文通過(guò)峰值信噪比PSNR來(lái)客觀評(píng)價(jià)圖像質(zhì)量[16],如式(13)所示。其中m、n為圖像維度,f0(i,j)表示原圖像,f1(i,j)表示待評(píng)價(jià)圖像。PSNR值越大,說(shuō)明失真越少,圖像質(zhì)量越高。結(jié)果如表1所示。

    (13)

    表1 圖像PSNR值比較

    從表1可以看出,本文所提算法能夠有效增強(qiáng)圖像細(xì)節(jié),提高目標(biāo)識(shí)別能力,有較好的增強(qiáng)效果,并且對(duì)圖像質(zhì)量的改善能力優(yōu)于單獨(dú)的雙邊濾波;串行實(shí)現(xiàn)和并行實(shí)現(xiàn)不會(huì)影響算法處理效果,兩者得到的圖像質(zhì)量是一致的。

    4.2 加速性能

    分別對(duì)不同尺寸的SAR圖像進(jìn)行加速性能測(cè)試,記錄算法串行實(shí)現(xiàn)和并行實(shí)現(xiàn)的計(jì)算時(shí)間。為了減小誤差,提高準(zhǔn)確程度,取10次運(yùn)算結(jié)果的平均值,計(jì)算加速比。實(shí)驗(yàn)結(jié)果如表2所示。

    表2 計(jì)算時(shí)間對(duì)比

    從表2可以看出,本文所提算法在GPU上的執(zhí)行時(shí)間明顯少于在CPU上的執(zhí)行時(shí)間,基于CUDA的并行化實(shí)現(xiàn)方法具有顯著的加速效果,1024×1024的圖像的加速比達(dá)到75.756,這是因?yàn)楸疚乃崴惴ň哂休^好的并行特性,能夠充分利用GPU的架構(gòu)和存儲(chǔ)結(jié)構(gòu)特點(diǎn),并行計(jì)算各像素的處理過(guò)程;隨著圖像尺寸的增大,加速比逐漸提高,這是因?yàn)橄到y(tǒng)調(diào)度和數(shù)據(jù)傳輸需要一定的時(shí)間開(kāi)銷,數(shù)據(jù)規(guī)模較大時(shí),雖然調(diào)用的線程較多,數(shù)據(jù)傳輸時(shí)間較長(zhǎng),但其在整體運(yùn)行時(shí)間中所占比重較小,因此加速比較大。

    5 結(jié)論

    本文提出一種基于噪聲可見(jiàn)性函數(shù)的SAR圖像增強(qiáng)快速算法及其基于GPU的并行化實(shí)現(xiàn)方法。實(shí)驗(yàn)結(jié)果表明,該算法能夠充分利用圖像分層的內(nèi)在并行性,實(shí)現(xiàn)并行加速,在有效提高圖像質(zhì)量,增強(qiáng)圖像細(xì)節(jié)的同時(shí),減少計(jì)算時(shí)間,進(jìn)一步滿足SAR圖像增強(qiáng)系統(tǒng)實(shí)時(shí)性要求。

    [1] 陳建飛,李躍華,王劍橋. 毫米波被動(dòng)合成孔徑成像系統(tǒng)[J].微波學(xué)報(bào),2012,S1:253-256.

    [2] 尹奎英. SAR圖像處理及地面目標(biāo)識(shí)別技術(shù)研究[D]. 西安:西安電子科技大學(xué),2011.

    [3] 陳 沫. SAR圖像軍事目標(biāo)識(shí)別方法研究[D]. 合肥:安徽大學(xué),2013.

    [4] 張 舒. GPU高性能運(yùn)算之CUDA[M]. 北京:中國(guó)水利水電出版社,2009.

    [5] 姚 旺,胡 欣,劉 飛,等. 基于GPU的高性能并行計(jì)算技術(shù)[J]. 計(jì)算機(jī)測(cè)量與控制,2014,22(12):4160-4162.

    [6] 宋曉麗,王 慶. 基于GPGPU的數(shù)字圖像并行化預(yù)處理[J]. 計(jì)算機(jī)測(cè)量與控制,2009,17(6):1169-1171.

    [7] Branchitta F,Diani M,Romagnoli M. New technique for the visualization of high dynamic range infrared images[J]. Optical Engineering,2009,48(9):6401.

    [8] Voloshynovskiy S,Pereira S,Iquise V,et al. Attack modelling: towards a second generation watermarking benchmark[J]. Signal Processing,2001,81(6):1177-1214.

    [9] 湯代佳,尚東方. 基于噪聲可見(jiàn)性函數(shù)的交通流量異常預(yù)警算法[J]. 上海電機(jī)學(xué)院學(xué)報(bào),2012,15(6):386-389.

    [10] 班陽(yáng)陽(yáng),張勁東,陳家瑞,等. 后向投影成像算法的 GPU 優(yōu)化方法研究[J]. 雷達(dá)科學(xué)與技術(shù),2014(6):659-665.

    [11] 張紹良,閆鈞華,劉 成,等. 基于CUDA的紅外圖像快速增強(qiáng)算法研究[J]. 電子設(shè)計(jì)工程,2012,20(17):153-157.

    [12] Zhao C,You W,Wang Y,et al. GPU implementation for real-time hyperspectral anomaly detection[A]. IEEE International Conference on Digital Signal Processing[C]. IEEE,2015:940-943.

    [13] Yiu B Y S,Tsang I K H,Yu A C H. Real-time GPU-based software beamformer designed for advanced imaging methods research[C]. Ultrasonics Symposium. IEEE,2011:1920-1923.

    [14] 李大禹,胡立發(fā),穆全全,等. CUDA架構(gòu)下的液晶自適應(yīng)波面數(shù)值解析[J]. 光學(xué)精密工程,2010,18(4):848-854.

    [15] 曾炫杰,陳強(qiáng),譚海鵬,等. 基于CUDA的加速雙邊濾波算法[J]. 計(jì)算機(jī)科學(xué),2015,42(s1).

    [16] 段 群,吳粉俠,劉小豫. 基于CUDA的椒鹽噪聲自適應(yīng)中值濾波算法[J]. 航空計(jì)算技術(shù),2013,43(2):60-63.

    猜你喜歡
    細(xì)節(jié)
    細(xì)節(jié)
    不讓細(xì)節(jié)打敗自己
    以細(xì)節(jié)取勝 Cambridge Audio AXR100/ FOCAL ARIA 906
    怎樣進(jìn)行細(xì)節(jié)描寫(xiě)
    留心細(xì)節(jié)處處美——《收集東·收集西》
    奇妙的細(xì)節(jié)
    用Adobe Camera Raw還原隱藏的細(xì)節(jié)
    夏日養(yǎng)生重細(xì)節(jié)
    幸福(2017年18期)2018-01-03 06:34:47
    關(guān)注48個(gè)生活細(xì)節(jié)
    細(xì)節(jié)取勝
    Coco薇(2016年10期)2016-11-29 19:59:58
    av免费在线看不卡| 99久久成人亚洲精品观看| 级片在线观看| 成人性生交大片免费视频hd| 嘟嘟电影网在线观看| 别揉我奶头 嗯啊视频| 亚洲欧美日韩无卡精品| 精品一区二区免费观看| 婷婷亚洲欧美| h日本视频在线播放| 欧美激情国产日韩精品一区| 国产欧美日韩精品一区二区| 桃色一区二区三区在线观看| 久久6这里有精品| 不卡一级毛片| 69人妻影院| 国产精品99久久久久久久久| 麻豆国产av国片精品| 村上凉子中文字幕在线| 欧美又色又爽又黄视频| 一区福利在线观看| 久久精品影院6| 久久久精品欧美日韩精品| 亚洲国产欧美人成| av在线亚洲专区| 国产黄a三级三级三级人| 99热这里只有是精品在线观看| 熟妇人妻久久中文字幕3abv| 激情 狠狠 欧美| 男女那种视频在线观看| 欧美成人一区二区免费高清观看| 免费观看人在逋| 91在线精品国自产拍蜜月| a级毛片a级免费在线| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品一区二区在线观看99 | av免费在线看不卡| 国产亚洲av片在线观看秒播厂 | 99riav亚洲国产免费| 亚洲七黄色美女视频| 亚洲高清免费不卡视频| av女优亚洲男人天堂| 如何舔出高潮| 性色avwww在线观看| 一区二区三区高清视频在线| 午夜免费激情av| 人妻制服诱惑在线中文字幕| 美女内射精品一级片tv| 91狼人影院| 国产一级毛片在线| 国产精品久久电影中文字幕| 国产综合懂色| 国产综合懂色| 成人欧美大片| 大又大粗又爽又黄少妇毛片口| 男插女下体视频免费在线播放| 一个人看的www免费观看视频| 亚洲国产欧美人成| 一个人看的www免费观看视频| 蜜桃亚洲精品一区二区三区| 亚洲图色成人| 欧美变态另类bdsm刘玥| av天堂中文字幕网| 天堂中文最新版在线下载 | 国产伦精品一区二区三区四那| 国语自产精品视频在线第100页| 淫秽高清视频在线观看| 久久久久久久久久久丰满| 麻豆国产av国片精品| 免费在线观看成人毛片| 精品人妻熟女av久视频| 欧美丝袜亚洲另类| 中文字幕精品亚洲无线码一区| 国产亚洲精品久久久com| 欧美在线一区亚洲| 18禁黄网站禁片免费观看直播| 精品日产1卡2卡| 99在线人妻在线中文字幕| 赤兔流量卡办理| 国产精品野战在线观看| 国语自产精品视频在线第100页| 欧美+亚洲+日韩+国产| 男女做爰动态图高潮gif福利片| 午夜免费男女啪啪视频观看| 久久久色成人| 国产真实伦视频高清在线观看| 国产乱人偷精品视频| 看非洲黑人一级黄片| 性欧美人与动物交配| 国产精品久久久久久精品电影| 一级毛片我不卡| av卡一久久| 成人无遮挡网站| 一进一出抽搐动态| 国产成人影院久久av| 一进一出抽搐动态| 在线免费观看不下载黄p国产| 午夜免费激情av| 久久精品综合一区二区三区| 国产色婷婷99| 欧美日韩乱码在线| av.在线天堂| 蜜桃久久精品国产亚洲av| 国产亚洲av片在线观看秒播厂 | 久久精品国产清高在天天线| 精品不卡国产一区二区三区| 久久国内精品自在自线图片| 亚洲av中文av极速乱| 国产69精品久久久久777片| 国产午夜福利久久久久久| 国产av不卡久久| 精品欧美国产一区二区三| 哪里可以看免费的av片| 国产一区二区在线观看日韩| 成人毛片a级毛片在线播放| 亚洲最大成人中文| 一级黄色大片毛片| 一本精品99久久精品77| 亚洲国产精品久久男人天堂| 3wmmmm亚洲av在线观看| 直男gayav资源| 91久久精品国产一区二区三区| 久久久久国产网址| 男人和女人高潮做爰伦理| av女优亚洲男人天堂| 别揉我奶头 嗯啊视频| 国产精品99久久久久久久久| av在线亚洲专区| 久久午夜亚洲精品久久| 看十八女毛片水多多多| 超碰av人人做人人爽久久| 国产精品不卡视频一区二区| 卡戴珊不雅视频在线播放| 伦理电影大哥的女人| 成人午夜高清在线视频| 嫩草影院入口| 一级av片app| 国产日韩欧美在线精品| 成熟少妇高潮喷水视频| 欧美成人精品欧美一级黄| 亚洲美女搞黄在线观看| 卡戴珊不雅视频在线播放| 国产蜜桃级精品一区二区三区| 少妇高潮的动态图| 男女那种视频在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲高清免费不卡视频| 一级av片app| 成人永久免费在线观看视频| 天天躁夜夜躁狠狠久久av| 深夜精品福利| 成人午夜高清在线视频| 国产精品免费一区二区三区在线| 欧美日本视频| 国产一区二区激情短视频| 一级毛片aaaaaa免费看小| 久久精品国产清高在天天线| 亚洲精华国产精华液的使用体验 | 99精品在免费线老司机午夜| 精品久久久久久久久久久久久| 人人妻人人看人人澡| 天堂av国产一区二区熟女人妻| 日韩一区二区视频免费看| 黄色一级大片看看| 九草在线视频观看| 三级国产精品欧美在线观看| 精品国产三级普通话版| 偷拍熟女少妇极品色| 国产淫片久久久久久久久| 国产v大片淫在线免费观看| 赤兔流量卡办理| 成人毛片60女人毛片免费| 我的女老师完整版在线观看| 欧美成人a在线观看| 久久九九热精品免费| 女的被弄到高潮叫床怎么办| 久久久久久久久中文| 最近最新中文字幕大全电影3| 国产色爽女视频免费观看| 国产色爽女视频免费观看| 美女国产视频在线观看| 国产av在哪里看| 国产91av在线免费观看| av在线播放精品| 性色avwww在线观看| 直男gayav资源| 日本成人三级电影网站| 男女下面进入的视频免费午夜| 欧美在线一区亚洲| 性插视频无遮挡在线免费观看| 中文精品一卡2卡3卡4更新| 寂寞人妻少妇视频99o| 亚洲最大成人中文| 国产精品永久免费网站| 久久久久网色| 久久中文看片网| 黄色一级大片看看| 国产人妻一区二区三区在| av在线天堂中文字幕| 男人和女人高潮做爰伦理| 99精品在免费线老司机午夜| 国产精品综合久久久久久久免费| 久久亚洲精品不卡| 国产视频首页在线观看| 1024手机看黄色片| 国产91av在线免费观看| 色5月婷婷丁香| 精品一区二区三区人妻视频| 99久久久亚洲精品蜜臀av| 亚洲中文字幕日韩| 久久久久久久亚洲中文字幕| 国产美女午夜福利| 一个人观看的视频www高清免费观看| 欧美精品国产亚洲| 给我免费播放毛片高清在线观看| 久久久久久久亚洲中文字幕| 1000部很黄的大片| 久久久久久久午夜电影| 国产亚洲精品久久久com| 国产精品久久久久久久久免| 国产黄色视频一区二区在线观看 | 欧美又色又爽又黄视频| 中文字幕av在线有码专区| 久久99热6这里只有精品| 国产精品久久视频播放| 舔av片在线| 久久婷婷人人爽人人干人人爱| 3wmmmm亚洲av在线观看| 99九九线精品视频在线观看视频| a级毛色黄片| www.av在线官网国产| 中文字幕av成人在线电影| 精品一区二区免费观看| a级毛片a级免费在线| 99久久精品国产国产毛片| 亚洲欧美日韩卡通动漫| 一边摸一边抽搐一进一小说| 在线免费观看的www视频| 日韩欧美一区二区三区在线观看| 久久精品人妻少妇| 三级经典国产精品| 国产伦精品一区二区三区视频9| 最后的刺客免费高清国语| 国产高清三级在线| 日本免费一区二区三区高清不卡| 婷婷色综合大香蕉| 3wmmmm亚洲av在线观看| 伦理电影大哥的女人| www.av在线官网国产| 国产v大片淫在线免费观看| 一级毛片久久久久久久久女| 欧美日韩在线观看h| 中文资源天堂在线| 亚洲精品乱码久久久v下载方式| 亚洲三级黄色毛片| 成人特级黄色片久久久久久久| 国产高清视频在线观看网站| 亚洲欧美清纯卡通| 国产色爽女视频免费观看| 国产高清不卡午夜福利| 日本一本二区三区精品| 欧美丝袜亚洲另类| 欧美成人免费av一区二区三区| 欧美成人免费av一区二区三区| 日本撒尿小便嘘嘘汇集6| 韩国av在线不卡| 国产人妻一区二区三区在| 日韩 亚洲 欧美在线| a级一级毛片免费在线观看| 久久人人精品亚洲av| 免费搜索国产男女视频| 精品久久久久久久久久免费视频| 少妇高潮的动态图| 国产精品国产高清国产av| 高清在线视频一区二区三区 | 一边亲一边摸免费视频| 黄片无遮挡物在线观看| 成年版毛片免费区| 国产在线精品亚洲第一网站| 成人性生交大片免费视频hd| 美女 人体艺术 gogo| 日韩在线高清观看一区二区三区| 久久热精品热| 乱人视频在线观看| 麻豆国产av国片精品| 一卡2卡三卡四卡精品乱码亚洲| 人妻制服诱惑在线中文字幕| 菩萨蛮人人尽说江南好唐韦庄 | 26uuu在线亚洲综合色| 久久久久久久久久黄片| 好男人视频免费观看在线| 毛片女人毛片| 麻豆av噜噜一区二区三区| 国产真实乱freesex| 欧美色视频一区免费| 国产人妻一区二区三区在| 国产黄色视频一区二区在线观看 | 蜜桃亚洲精品一区二区三区| 在线观看免费视频日本深夜| 国产一区二区在线av高清观看| 一级av片app| 你懂的网址亚洲精品在线观看 | 婷婷色综合大香蕉| 亚洲av男天堂| 18禁黄网站禁片免费观看直播| 国产爱豆传媒在线观看| 国产成人freesex在线| 精品久久久久久久久久免费视频| 免费在线观看成人毛片| 一区二区三区免费毛片| 狂野欧美激情性xxxx在线观看| 九九久久精品国产亚洲av麻豆| 天天躁日日操中文字幕| 天堂av国产一区二区熟女人妻| av国产免费在线观看| 久久午夜福利片| 亚洲不卡免费看| 国产成人freesex在线| 欧美日韩综合久久久久久| 成年女人看的毛片在线观看| 成人毛片60女人毛片免费| 国产精品无大码| 久久精品夜夜夜夜夜久久蜜豆| 看非洲黑人一级黄片| 女人十人毛片免费观看3o分钟| 国内精品久久久久精免费| 久久国产乱子免费精品| 男插女下体视频免费在线播放| 国产精品不卡视频一区二区| 97人妻精品一区二区三区麻豆| 亚洲国产欧洲综合997久久,| АⅤ资源中文在线天堂| or卡值多少钱| 97超视频在线观看视频| 男人舔女人下体高潮全视频| 亚洲人成网站在线播放欧美日韩| 如何舔出高潮| 18禁在线无遮挡免费观看视频| 观看美女的网站| 免费在线观看成人毛片| 欧美日韩综合久久久久久| 99九九线精品视频在线观看视频| 老女人水多毛片| 2022亚洲国产成人精品| 听说在线观看完整版免费高清| 亚洲一区高清亚洲精品| 亚洲欧美成人精品一区二区| 欧美日韩一区二区视频在线观看视频在线 | 国产极品精品免费视频能看的| 亚洲国产精品sss在线观看| 一本久久中文字幕| 中国美白少妇内射xxxbb| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产亚洲av天美| 校园春色视频在线观看| 国产久久久一区二区三区| 中文字幕免费在线视频6| 18禁在线播放成人免费| 国产精品免费一区二区三区在线| 久久99精品国语久久久| 九九久久精品国产亚洲av麻豆| 男女做爰动态图高潮gif福利片| 乱人视频在线观看| 亚洲精品乱码久久久v下载方式| 亚洲精品自拍成人| 国产亚洲av嫩草精品影院| 18+在线观看网站| 男的添女的下面高潮视频| 成人特级av手机在线观看| 国产精品一区www在线观看| 日本在线视频免费播放| 十八禁国产超污无遮挡网站| 久久久久免费精品人妻一区二区| 日本黄色视频三级网站网址| 2022亚洲国产成人精品| 美女黄网站色视频| 1024手机看黄色片| 网址你懂的国产日韩在线| 欧美性猛交黑人性爽| 真实男女啪啪啪动态图| 看非洲黑人一级黄片| 麻豆国产av国片精品| 91精品一卡2卡3卡4卡| 亚洲中文字幕一区二区三区有码在线看| 日本与韩国留学比较| 看非洲黑人一级黄片| 黑人高潮一二区| 美女xxoo啪啪120秒动态图| 国产精品久久久久久av不卡| 一级毛片我不卡| 免费在线观看成人毛片| 国产精品久久久久久久久免| 亚洲第一电影网av| 国产爱豆传媒在线观看| 久久午夜福利片| 国产在线精品亚洲第一网站| 日日啪夜夜撸| 久久午夜亚洲精品久久| 免费av不卡在线播放| 一个人免费在线观看电影| 小蜜桃在线观看免费完整版高清| www.色视频.com| 99久久精品国产国产毛片| 日韩欧美在线乱码| 精品人妻视频免费看| 91久久精品国产一区二区成人| 国产亚洲欧美98| 国产成人a区在线观看| av天堂中文字幕网| 久久久成人免费电影| 能在线免费看毛片的网站| 人妻久久中文字幕网| 22中文网久久字幕| 99国产精品一区二区蜜桃av| 国产精品久久久久久久久免| 伊人久久精品亚洲午夜| 亚洲国产欧美在线一区| 三级毛片av免费| 日韩视频在线欧美| 久久精品国产自在天天线| 麻豆乱淫一区二区| 一区二区三区高清视频在线| 六月丁香七月| 麻豆成人av视频| 国产精品一及| 性欧美人与动物交配| 1000部很黄的大片| 91久久精品国产一区二区成人| 久久久久久久久中文| 国产私拍福利视频在线观看| 久久99热6这里只有精品| 国产熟女欧美一区二区| 免费观看人在逋| 欧美3d第一页| 美女脱内裤让男人舔精品视频 | 一进一出抽搐动态| 悠悠久久av| 国内精品宾馆在线| 成人欧美大片| av在线观看视频网站免费| 精品久久久久久久末码| 大香蕉久久网| 97在线视频观看| 在线国产一区二区在线| 深夜a级毛片| 18禁裸乳无遮挡免费网站照片| 国产精品乱码一区二三区的特点| 国产激情偷乱视频一区二区| 高清午夜精品一区二区三区 | 精品国内亚洲2022精品成人| 乱人视频在线观看| 国产白丝娇喘喷水9色精品| 一区二区三区免费毛片| 久久久久久九九精品二区国产| a级毛色黄片| 久久精品夜夜夜夜夜久久蜜豆| av在线老鸭窝| 久久久精品欧美日韩精品| 黄色配什么色好看| 亚洲av免费高清在线观看| 日韩人妻高清精品专区| 国产精品伦人一区二区| 精品一区二区免费观看| 免费人成在线观看视频色| or卡值多少钱| 婷婷精品国产亚洲av| 1024手机看黄色片| 欧美日韩在线观看h| 国产乱人偷精品视频| 床上黄色一级片| 久久久久久久久久成人| 国产成人影院久久av| 午夜福利高清视频| 成熟少妇高潮喷水视频| 亚洲熟妇中文字幕五十中出| 久久亚洲国产成人精品v| 免费搜索国产男女视频| 国产高潮美女av| 欧美3d第一页| 国产精品,欧美在线| 中文字幕av在线有码专区| 女人十人毛片免费观看3o分钟| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美日韩高清在线视频| 久久综合国产亚洲精品| 美女脱内裤让男人舔精品视频 | 不卡一级毛片| 亚洲一区高清亚洲精品| 国产精品一及| 亚洲av中文字字幕乱码综合| 丝袜美腿在线中文| 日韩欧美 国产精品| 桃色一区二区三区在线观看| 久久久久久久久大av| 小蜜桃在线观看免费完整版高清| 免费搜索国产男女视频| ponron亚洲| 干丝袜人妻中文字幕| 高清日韩中文字幕在线| 少妇的逼水好多| 99久久无色码亚洲精品果冻| 我的老师免费观看完整版| 欧美3d第一页| 日韩欧美三级三区| 亚洲性久久影院| 一级毛片我不卡| 男人舔奶头视频| 欧洲精品卡2卡3卡4卡5卡区| 成年女人永久免费观看视频| 久久鲁丝午夜福利片| 男人和女人高潮做爰伦理| 亚洲自拍偷在线| 又爽又黄a免费视频| 久久久久久大精品| 三级男女做爰猛烈吃奶摸视频| 精品免费久久久久久久清纯| 亚洲第一区二区三区不卡| 少妇人妻精品综合一区二区 | 欧美丝袜亚洲另类| 亚洲av.av天堂| 97超碰精品成人国产| 一区二区三区免费毛片| 午夜福利在线在线| 一本精品99久久精品77| 深爱激情五月婷婷| 亚洲电影在线观看av| 日本黄色片子视频| 亚洲av免费在线观看| 麻豆成人av视频| 春色校园在线视频观看| 欧美日韩在线观看h| 久久久久久久亚洲中文字幕| 淫秽高清视频在线观看| 国产老妇女一区| 麻豆乱淫一区二区| 国产精品,欧美在线| 尾随美女入室| 久久久久九九精品影院| 欧美日韩在线观看h| 成人美女网站在线观看视频| 国产黄片视频在线免费观看| 国产一级毛片七仙女欲春2| 桃色一区二区三区在线观看| 97超视频在线观看视频| 精品少妇黑人巨大在线播放 | 一级毛片aaaaaa免费看小| 日本爱情动作片www.在线观看| 精品一区二区三区人妻视频| 美女 人体艺术 gogo| 国产蜜桃级精品一区二区三区| 欧美高清成人免费视频www| 日韩高清综合在线| 欧美日韩在线观看h| 黄色欧美视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产成人精品久久久久久| 亚洲欧美精品综合久久99| 变态另类丝袜制服| 亚洲经典国产精华液单| 最近的中文字幕免费完整| 亚洲欧美成人综合另类久久久 | 日本与韩国留学比较| 99热这里只有是精品50| 久久6这里有精品| 国产乱人视频| 又爽又黄a免费视频| 美女xxoo啪啪120秒动态图| 午夜福利在线观看免费完整高清在 | ponron亚洲| 色尼玛亚洲综合影院| 黄色视频,在线免费观看| 日本-黄色视频高清免费观看| 女人被狂操c到高潮| 国产av不卡久久| 一个人看的www免费观看视频| 人妻系列 视频| 久久中文看片网| 亚洲欧美日韩东京热| 国产美女午夜福利| 高清日韩中文字幕在线| 国产成人freesex在线| 国产精品一区www在线观看| 久久婷婷人人爽人人干人人爱| 亚洲国产日韩欧美精品在线观看| 亚洲最大成人手机在线| 深夜a级毛片| 久久这里有精品视频免费| 在线观看午夜福利视频| 少妇人妻一区二区三区视频| 亚洲丝袜综合中文字幕| 在线观看午夜福利视频| 自拍偷自拍亚洲精品老妇| 免费观看精品视频网站| 午夜福利在线观看吧| 亚洲欧美日韩无卡精品| 99久久精品一区二区三区| 美女内射精品一级片tv| 欧美日韩在线观看h| 久久这里有精品视频免费| 插阴视频在线观看视频| 欧美日韩精品成人综合77777| 美女被艹到高潮喷水动态| 嫩草影院入口| 亚洲欧美日韩无卡精品| 天堂av国产一区二区熟女人妻| 性色avwww在线观看| 日韩中字成人| 国产亚洲精品av在线| 午夜福利视频1000在线观看| 男女下面进入的视频免费午夜| 日韩在线高清观看一区二区三区| 男人和女人高潮做爰伦理| 国产伦精品一区二区三区视频9| 麻豆成人av视频| 国产黄色视频一区二区在线观看 |