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

    GPU在眼科FD-OCT系統(tǒng)實時圖像顯示及數(shù)據(jù)處理中的應用

    2013-12-05 09:38:06劉巧艷李躍杰徐秋晶趙金城王立偉高用賀
    中國醫(yī)療器械雜志 2013年1期
    關鍵詞:掃描模式插值數(shù)據(jù)處理

    【作 者】劉巧艷,李躍杰,徐秋晶,趙金城,王立偉,高用賀

    1 北京協(xié)和醫(yī)學院,北京市,100730

    2 中國醫(yī)學科學院生物醫(yī)學工程研究所,天津市,300192

    3 天津市眼科醫(yī)學設備技術工程中心,天津市,300384

    0 引言

    光學相干層析成像技術(Optical Coherence Tomography),自上世紀90年代被成功用于眼科疾病診斷[1]后,得到了迅速發(fā)展:由時域OCT階段發(fā)展到頻域OCT階段[2];由組織結構成像向組織功能成像發(fā)展[3-4];由眼科診斷拓展到心血管、皮膚、口腔、組織工程等領域中的應用[5-6]。

    隨著超高速CMOS線陣掃描相機的發(fā)展,頻域OCT光譜譜線轉換及線采樣率已經可以達到300 k線/s[7],為臨床OCT系統(tǒng)實時成像提供了前提。目前影響眼科OCT系統(tǒng)實時成像的技術瓶頸是需要先將采樣數(shù)據(jù)進行頻譜域空間(λ空間)到波數(shù)空間(K空間)變換,進行插值變換和FFT變換,然后再將變換后的數(shù)據(jù)進行2D或3D成像。由于成像的數(shù)據(jù)量很大,特別是進行C模式掃描成像(如眼底視網膜en-face成像模式)時,需要先將獲得的3D圖像數(shù)據(jù)進行處理后,再將得到的數(shù)據(jù)成像。因此,如何提高數(shù)據(jù)處理的速度,成為眼科OCT系統(tǒng)實現(xiàn)實時成像的關鍵。

    為了克服眼科OCT系統(tǒng)實時成像的技術瓶頸,近年來研究人員進行了大量研究,并提出了一些相關的解決方法。有研究人員將多核CPU引入到OCT系統(tǒng)進行采樣數(shù)據(jù)的并行處理,對于1024個采樣點模式OCT系統(tǒng),非均勻K空間數(shù)據(jù)處理的速度可達到80 k線/s[8],均勻K空間數(shù)據(jù)的處理速度可達到207 k線/s[9]。也有研究人員通過在原有系統(tǒng)中增加DSPs[10]或者FPGAs[11]硬件模塊,加快實現(xiàn)數(shù)據(jù)處理過程。隨著圖形處理單元(GPU)技術的發(fā)展,其運算能力和可編程性得到大幅度提高,它除了在傳統(tǒng)的圖像處理領域應用繼續(xù)保持優(yōu)勢外,作為通用的并行計算處理器已被廣泛用于科學計算、工程、金融以及其他領域中。目前,有研究人員將GPU引入到OCT系統(tǒng)中,利用它強大的并行計算能力來加速采樣數(shù)據(jù)處理過程,解決OCT系統(tǒng)實時成像的技術瓶頸。Kang Zhang和Jin U.Kang將GPU引入到OCT系統(tǒng)中,在專業(yè)圖像工作站平臺上采用線性插值算法,在1024個采樣點模式下數(shù)據(jù)處理速度達到680 k線/s,2048個采樣點模式下數(shù)據(jù)處理速度可達到320 k線/s,達到了3D實時成像對于處理速度的要求[12]?,F(xiàn)階段將GPU應用于OCT系統(tǒng)還處于實驗室研究階段,為獲得更高處理速度,往往需要配置高性能的專用圖像工作站(配置多核處理器)和GPU,硬件成本較高。

    本課題在不改變實驗室現(xiàn)有硬件平臺的條件下,把低成本GPU引入到我們正在開發(fā)應用的眼科OCT儀器中,實現(xiàn)儀器性能的大幅度提高,解決眼科OCT系統(tǒng)實時成像的問題。

    1 系統(tǒng)組成和工作原理

    圖1 FD-OCT成像系統(tǒng)組成Fig.1 FD-OCT imaging system configuration

    1.1 系統(tǒng)組成

    本系統(tǒng)組成如圖1所示,主要包括眼科OCT信號采集系統(tǒng)和信號處理系統(tǒng)兩大部分。其中眼科OCT信號采集系統(tǒng)采用寬帶超亮發(fā)光二極管(SLD)(λ0=840 nm,△λ= 50 nm)作為系統(tǒng)光源,以2048像素CMOS線陣相機(采樣速率為70 k線/s)作為光譜儀的檢測器。信號處理系統(tǒng)由計算機和GPU組成。計算機CPU為 Intel Celeron Dual-Core E3400 @ 2.60 GHz,2 G內存;GPU采用 NVIDIA 公司的GeForce GTX 460顯卡(336個流處理器,1 GB顯存)。

    1.2 系統(tǒng)工作原理

    光源發(fā)出的光經過50:50的光纖分束器后,被均勻分成兩束光,分別進入OCT系統(tǒng)的參考臂和樣品臂。從樣品臂反射回來的信號光和從參考臂返回的參考光再次經過光纖分束器匯合后發(fā)生干涉。包含樣品不同深度信息的干涉信號光譜,經光譜儀的CMOS線陣掃描相機采集,并通過相機數(shù)據(jù)線傳輸?shù)接嬎銠C,由計算機內的圖像采集卡對干涉信號光譜進行A/D轉換,并將轉換結果作為采樣數(shù)據(jù)存入到計算機內存中。將采樣數(shù)據(jù)通過PCIE×16總線傳輸?shù)紾PU顯存,借助GPU強大的并行數(shù)據(jù)處理能力進行數(shù)據(jù)處理,并將處理好的結果數(shù)據(jù)送回計算機進行圖像顯示。顯示的圖像包含了檢測樣品不同深度的結構信息。

    2 基于GPU加速技術的數(shù)據(jù)處理流程

    2.1 CUDA架構

    CUDA(Compute Unified Device Architecture)是一種由NVIDIA公司推出的通用并行計算架構,該架構使GPU能夠解決復雜的計算問題。在CUDA架構下,開發(fā)人員可以通過CUDA C語言(CUDA C語言是對標準C語言的一種簡單擴展)對GPU編程[13]。

    CUDA架構中,CPU作為主機(Host),GPU作為協(xié)處理器或者設備(Device)。在一個系統(tǒng)中可以存在一個主機和多個設備。 CPU 主要負責進行邏輯性強的事物處理和串行計算,GPU 則專注于執(zhí)行高度線程化的并行處理任務。CPU 、GPU 各自擁有相互獨立的存儲器地址空間:主機端的內存和設備端的顯存;在CUDA程序中,將運行在GPU上一個可以被并行執(zhí)行的步驟稱為kernel(內核函數(shù))。

    圖2 CPU-GPU系統(tǒng)數(shù)據(jù)處理流程圖Fig.2 Signal processing flow chart of CPU-GPU hybrid system architecture

    2.2 數(shù)據(jù)處理過程

    本系統(tǒng)的數(shù)據(jù)處理過程如圖2所示,其中實箭頭所指方向代表數(shù)據(jù)在不同設備間的流動方向,空箭頭所指方向代表數(shù)據(jù)在GPU內部的流動方向。系統(tǒng)數(shù)據(jù)處理過程包括對采樣數(shù)據(jù)進行預處理、頻譜域空間(λ空間)到波數(shù)空間(K空間)變換、插值變換、FFT變換和POST-FFT變換。

    在頻域OCT系統(tǒng)中,采樣數(shù)據(jù)是通過對OCT的光路系統(tǒng)掃描由相機采集到的,掃描一次得到一列數(shù)據(jù)(一個A-SCAN)。處理時是一列一列數(shù)據(jù)進行處理的。針對每列數(shù)據(jù)彼此相互獨立、可以并行處理的特點,利用CUDA架構將OCT系統(tǒng)整個數(shù)據(jù)處理過程改寫成適合在GPU上執(zhí)行的kernel函數(shù),大大提高了數(shù)據(jù)處理速度,從而達到系統(tǒng)實時成像的要求。

    2.3 程序設計

    首先確定系統(tǒng)數(shù)據(jù)處理過程中的串行部分和并行部分,選擇合適的算法,并按照算法確定數(shù)據(jù)和任務的劃分方式,將每個需要并行實現(xiàn)的步驟映射為一個滿足CUDA兩層并行模型的內核函數(shù)。其中頻譜域空間(λ空間)到波數(shù)空間(K空間)變換模塊在主程序中執(zhí)行;數(shù)據(jù)預處理、插值運算、POST-FFT運算需要重新改寫成適合在GPU上執(zhí)行的kernel函數(shù);FFT運算部分采用CUDA自帶的CUFFT庫函數(shù)來進行計算,CUDA4.1版本目前已經支持FFT雙精度運算。整個程序設計如下:

    (1)給采樣數(shù)據(jù)和計算過程中的各個變量分配空間,包括內存空間和顯存空間,在Host端CPU程序中,準備好計算要用到的數(shù)據(jù),將數(shù)據(jù)從主機內存?zhèn)鬏數(shù)斤@存中。

    (2)數(shù)據(jù)預處理包括數(shù)據(jù)類型轉換和去噪運算。去噪運算是將每一列數(shù)據(jù)組成的數(shù)組(一個A-SCAN)都減去一列同大小的噪聲數(shù)組。

    (3)插值運算實現(xiàn)了在進行FFT變換之前,采樣數(shù)據(jù)K空間分布的均勻化。系統(tǒng)最常用的插值算法包括最鄰近插值算法、線性插值算法和三次樣條插值算法。

    (4)利用CUDA自帶的CUFFT庫函數(shù)對插值運算結果進行FFT變換。

    (5)POST-FFT運算包括對FFT運算結果取模、取對數(shù)并進行歸一化。POST-FFT運算結果存儲在體積數(shù)組中。

    (6)根據(jù)不同成像平面的需要,抽取體積數(shù)組數(shù)據(jù)作為GPU結果數(shù)據(jù)。

    (7)將結果數(shù)據(jù)從顯存?zhèn)鬏數(shù)街鳈C內存中。

    2.4 優(yōu)化分析

    在用CUDA對GPU上的kernel函數(shù)進行改寫時,采用了一些優(yōu)化方法:

    (1)為了使內存和顯存之間的數(shù)據(jù)傳輸速度更快,在分配主機端內存的時候采用了pinned memory。

    (2)使用流運算來隱藏主機端和設備端數(shù)據(jù)通信的時間。

    (3)合理利用讀取速度更快的shared memory來存放計算過程中的中間變量。

    (4)將已經計算好的在插值運算過程中值保持不變的采樣數(shù)據(jù)的等間隔和非等間隔橫坐標值存放在常數(shù)寄存器里面。

    (5)根據(jù)系統(tǒng)資源進行grid和block維度設計。

    (6)使用CUDA profiler對CUDA程序進行性能測試,對耗時較長的模塊進行優(yōu)化。

    3 實驗結果

    本系統(tǒng)采用Microsoft Visual Studio2010中集成CUDA Toolkit 32 bit 4.1、CUDA SDK 32 bit 4.1和Nvidia Driver for Windows732 bit為開發(fā)環(huán)境,對采樣數(shù)據(jù)進行B掃描模式和C掃描模式成像。

    3.1 B掃描模式成像

    B掃描模式成像圖像能提供視網膜斷層結構,能清晰地顯示視網膜各層細微結構及病理改變,并作出定性或定量分析,目前已成為視網膜疾病和青光眼強有力的診斷工具[14]。

    我們采用100幀共計195 MB數(shù)據(jù)(每幀數(shù)據(jù)大小為500線×2048 像素/線×2 字節(jié)/像素)進行B掃描模式成像。分別采用線性插值算法和三次樣條插值算法,利用CUDA提供的計時函數(shù)分別對CPU模式和CPU-GPU模式下系統(tǒng)單幀B掃描模式圖像成像時間進行計時(計算100幀圖像成像時間取平均),實驗結果如表1所示。 從表1可知,采用GPU+CPU模式執(zhí)行成像數(shù)據(jù)處理的速度較CPU模式執(zhí)行同樣數(shù)據(jù)處理的速度提高超過一個數(shù)量級,其中采用線性插值算法速度提高了60倍,采用三次樣條插值算法速度提高了35倍。實驗成像效果圖如圖3所示,其中(a)為采用線性插值算法在CPU模式下系統(tǒng)成像圖像;(b)為采用線性插值算法在CPU+GPU模式下系統(tǒng)成像圖像;(c)為采用三次樣條插值算法在CPU模式下系統(tǒng)成像圖像;(d)為采用三次樣條插值算法在CPU+GPU模式下系統(tǒng)成像圖像。由圖3可知,采用相同的插值算法,在CPU模式和CPU-GPU模式下系統(tǒng)成像圖像質量無差異,采用三次樣條插值算法成像圖像質量比采用線性插值算法成像圖像質量要好。

    圖3 視網膜B掃描成像圖像Fig.3 B-scan imaging images of the retina

    插值運算是系統(tǒng)數(shù)據(jù)處理過程中計算量最大、耗時最長的一個環(huán)節(jié)。插值算法的選取不僅影響系統(tǒng)加速比,也影響系統(tǒng)成像速度和成像圖像質量。采用線性插值算法,成像速度快,成像圖像質量差;采用三次樣條算法,成像速度慢、成像圖像質量好。

    表1 CPU和GPU 成像速度對比Tab.1 Comparison of CPU and GPU imaging speed

    3.2 C掃描模式成像

    C掃描模式成像圖像能直觀地顯示眼底視網膜血管和黃斑等眼底組織的結構信息,臨床上可用于診斷眼底病變,如黃斑病變[15]、眼底神經組織結構變化等,可用于對視網膜下新生血管的深度和層次進行準確定位[16]。

    進行C掃描模式成像時,我們先將采樣數(shù)據(jù)進行B掃描模式處理,并將處理后的結果數(shù)據(jù)組織成3D紋理數(shù)組存放在顯存之中。根據(jù)不同的需要,對3D數(shù)組進行切片提取顯示en-face單層圖像,或者通過多切片疊加取平均顯示眼底組織結構信息。

    我們采用的3D數(shù)據(jù)塊為480幀共計204 MB數(shù)據(jù)(每幀數(shù)據(jù)大小為249線×896 像素/線×2 字節(jié)/像素)進行C掃描模式成像。3D數(shù)據(jù)塊平均成像速度為1.8 s,能快速實現(xiàn)視網膜en-face單層切片成像或en-face多切片疊加平均成像。成像效果如圖4所示,其中圖(a)~(c)分別顯示了視網膜en-face單層切片圖像;沿縱軸方向看,(b)切片比圖4(a)切片深約30 μm,(c)切片比(b)切片深約30 μm。圖4(a)~(c)清晰地顯示了視網膜血管、微血管、黃斑等眼底組織細微結構。視網膜en-face單層切片成像的軸向分辨率可達到5 μm。與眼底相機只能對眼底復合結構成像相比,視網膜en-face單層切片成像可以對眼底視網膜下的細微組織的深度和層次進行準確定位。圖4(d)~(f)分別顯示了視網膜11個en-face 切片疊加后再取平均的成像圖像。沿縱軸方向看,圖4(d)~(f)、分別以圖4(a)~(c)的切片為中心,前后各取5切片疊加后取平均所成圖像。從圖4(d)~(f)能清晰地觀察到視網膜血管、微血管、黃斑等眼底組織細微結構,與圖4(a)~(c)單切片成像圖相比,圖4(d)~(f)包含了更多的復合結構信息。

    圖4 視網膜en-face成像圖像Fig.4 En-face imaging images of the retina

    實驗結果表明:在成像圖像質量不變的前提下,在眼科OCT系統(tǒng)中,采用GPU+CPU模式實現(xiàn)成像數(shù)據(jù)處理的速度,較傳統(tǒng)的基于CPU平臺的串行計算和成像模式執(zhí)行同樣數(shù)據(jù)處理的速度提高超過一個數(shù)量級。

    4 結語

    本文利用實驗室現(xiàn)有標準OCT系統(tǒng)平臺,在未增加硬件的基礎上,利用計算機通用顯卡GPU,并將基于GPU的統(tǒng)一計算設備架構(CUDA)引入到眼科OCT系統(tǒng)成像中的數(shù)據(jù)處理過程,借助GPU強大的并行數(shù)據(jù)處理能力和浮點計算能力,用CUDA對OCT系統(tǒng)數(shù)據(jù)處理過程進行改寫,使得眼科OCT系統(tǒng)的成像速度較基于CPU平臺處理成像速度提高了數(shù)十倍,達到了臨床2D實時成像的要求,為眼科3D實時成像打下了基礎。

    [1]D.Huang,E.A.Swanson,C.P.Lin,et al.Optical coherence tomography[J].Science,1991,254(5035): 1178-1181.

    [2]R.Leitgeb, C.K.Hitzenberger, A.F.Fercher.Performance of fourier domain vs.time domain optical coherence tomography[J].Opt Express,2003,11(8): 889–894.

    [3]Ruikang,K.Wang.In vivo structural and flow imaging: US,US20100027857A1[P],2010-2-4.

    [4]B.Cense,T.C.Chen,B.H.Park,et al.Thickness and birefringence of healthy retinal nerve fiber layer tissue measure measured with polarization-sensitive optical coherence tomography[J].Invest Ophth Vis Sci,2004,45(12): 2606-2612.

    [5]N.D.Gladkova,G.A.Petrova,N.K.Nikulin,et al.In vivo optical coherence tomography imaging of human skin: norm and pathology[J].Skin Res Technol,2000,6(1): 6-16.

    [6]L.Vabre,A.Dubois,A.C.Boccara.Thermal-light full-field optical coherence tomography[J].Opt Lett,2002,27(7): 530-532.

    [7]B.Potsaid,I.Gorczynska,V.J.Srinivasan,et al.Ultrahigh speed spectral/Fourier domain OCT ophthalmic imaging at 70,000 to 312,500 axial scans per second[J].Opt Lett,2008,16(19): 15149–15169.

    [8]G.Liu,J.Zhang,L.Yu,et al.Real-time polarization-sensitive optical coherence tomography data processing with parallel computing[J]. Appl Opt,2009,48(32): 6365–6370.

    [9]J.Probst,P.Koch,G.Huttmann.Real-time 3D rendering of optical coherence tomography volumetric data[C].Proc SPIE,2009,7372,73720Q.

    [10]A.W.Schaefer,J.J.Reynolds,D.L.Marks,et al.Real-time digitalsignal processing-based optical coherence tomography and Doppler optical coherence tomography[J].IEEE Trans Biomed Eng,2004,54: 186-190.

    [11]T.E.Ustun,N.V.Iftimia,R.D.Ferguson,et al.Real-time processing for Fourier domain optical coherence tomography using a field programmable gate array[J].Rev Sci Instrum,2008,79:114301-114310.

    [12]Kang Zhang,Jin U.Kang,Real-time 4D signal processing and visualization using graphics processing unit on a regular nonlinear-k Fourier-domain OCT system[J], Opt Express,2010,18(11): 11772-11784.

    [13]Jason Sanders,Edward Kandrot.GPU高性能編程CUDA實戰(zhàn)[M].北京: 機械工業(yè)出版社,2011

    [14]J.S.Schuman,M.R.Hee,C.A.Puliafito.et al.Quantification of nerve fiber layer thickness in normal and glaucomatous eyes using optical coherence tomography[J].Arch Ophth,1995,113(5): 586-596.

    [15]M.Altaweel,M.Ip.Macular hole: improved understanding of pathogenesis,staging,and management based on optical coherence tomography[J].Semin Ophth,2003,18(2): 58-66.

    [16]紀淑興,張軍軍,唐健,等.中心性滲出性脈絡膜視網膜病變的光學相干斷層掃描圖像[J].中華眼底病雜志,2002,18(2): 121-124.

    猜你喜歡
    掃描模式插值數(shù)據(jù)處理
    認知診斷缺失數(shù)據(jù)處理方法的比較:零替換、多重插補與極大似然估計法*
    心理學報(2022年4期)2022-04-12 07:38:02
    ILWT-EEMD數(shù)據(jù)處理的ELM滾動軸承故障診斷
    水泵技術(2021年3期)2021-08-14 02:09:20
    氣質聯(lián)用中性丟失掃描定性含氯化合物的方法
    基于Sinc插值與相關譜的縱橫波速度比掃描方法
    雙光能X射線骨密度儀測量腰椎骨密度不同掃描模式的對比研究*
    一種改進FFT多譜線插值諧波分析方法
    基于四項最低旁瓣Nuttall窗的插值FFT諧波分析
    基于希爾伯特- 黃變換的去噪法在外測數(shù)據(jù)處理中的應用
    Blackman-Harris窗的插值FFT諧波分析與應用
    基于POS AV610與PPP的車輛導航數(shù)據(jù)處理
    亚洲国产毛片av蜜桃av| 在线免费观看不下载黄p国产| 香蕉精品网在线| 国产男女超爽视频在线观看| 一级爰片在线观看| av国产精品久久久久影院| 亚洲精品aⅴ在线观看| 午夜免费男女啪啪视频观看| 国产亚洲av片在线观看秒播厂| 国产成人aa在线观看| 婷婷色综合www| 少妇人妻 视频| 日韩av在线免费看完整版不卡| 国产色爽女视频免费观看| 黑人欧美特级aaaaaa片| 黑人高潮一二区| 免费播放大片免费观看视频在线观看| 亚洲av国产av综合av卡| 亚洲精品一区蜜桃| 中文乱码字字幕精品一区二区三区| 亚洲色图综合在线观看| 大片电影免费在线观看免费| 国产精品一国产av| 久久久久久久久久久久大奶| 高清黄色对白视频在线免费看| 久久久国产精品麻豆| 捣出白浆h1v1| 爱豆传媒免费全集在线观看| 夫妻性生交免费视频一级片| 尾随美女入室| 国产av一区二区精品久久| 免费少妇av软件| 巨乳人妻的诱惑在线观看| 一级a做视频免费观看| 日本wwww免费看| 看非洲黑人一级黄片| 中文字幕人妻熟女乱码| 一级片免费观看大全| 你懂的网址亚洲精品在线观看| 亚洲人成网站在线观看播放| 自线自在国产av| 亚洲精品456在线播放app| 最近2019中文字幕mv第一页| 丝袜美足系列| 久久人人爽av亚洲精品天堂| 国产精品国产三级国产av玫瑰| 搡女人真爽免费视频火全软件| 久久精品国产自在天天线| 丝袜人妻中文字幕| 亚洲激情五月婷婷啪啪| 国产有黄有色有爽视频| 国产探花极品一区二区| 久久亚洲国产成人精品v| 丝袜喷水一区| 久久久a久久爽久久v久久| 建设人人有责人人尽责人人享有的| 亚洲av免费高清在线观看| 国产av国产精品国产| 777米奇影视久久| 亚洲精品,欧美精品| 夫妻性生交免费视频一级片| 亚洲,欧美精品.| 午夜视频国产福利| 日韩制服丝袜自拍偷拍| 中文字幕制服av| 999精品在线视频| 国产69精品久久久久777片| 亚洲欧美日韩另类电影网站| 午夜激情av网站| 在线 av 中文字幕| 日韩欧美精品免费久久| 久久韩国三级中文字幕| 久久影院123| 最近最新中文字幕免费大全7| 少妇人妻 视频| 又大又黄又爽视频免费| 黄片无遮挡物在线观看| 美女脱内裤让男人舔精品视频| 18+在线观看网站| 一级爰片在线观看| 纯流量卡能插随身wifi吗| 亚洲成人一二三区av| 狂野欧美激情性bbbbbb| 欧美人与性动交α欧美软件 | 国产精品欧美亚洲77777| 国产一区二区三区av在线| 日韩av不卡免费在线播放| 免费观看在线日韩| 看免费成人av毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 热99国产精品久久久久久7| 一级毛片我不卡| 两个人看的免费小视频| 久久久精品94久久精品| 高清av免费在线| 日本猛色少妇xxxxx猛交久久| 丰满迷人的少妇在线观看| 成人漫画全彩无遮挡| av国产久精品久网站免费入址| 欧美人与性动交α欧美软件 | 波野结衣二区三区在线| av.在线天堂| 久久人妻熟女aⅴ| 寂寞人妻少妇视频99o| 国产免费福利视频在线观看| 亚洲精品美女久久av网站| 捣出白浆h1v1| 色94色欧美一区二区| 大香蕉久久网| 精品人妻熟女毛片av久久网站| 精品少妇黑人巨大在线播放| 免费看光身美女| 少妇高潮的动态图| 日韩熟女老妇一区二区性免费视频| 久久国产精品大桥未久av| 国产精品久久久久久久电影| 在线天堂中文资源库| 久久久国产精品麻豆| 国产av码专区亚洲av| 五月开心婷婷网| 免费看光身美女| 又黄又爽又刺激的免费视频.| 自线自在国产av| 久久久a久久爽久久v久久| 国产色婷婷99| 赤兔流量卡办理| 国产又色又爽无遮挡免| 欧美97在线视频| 精品国产一区二区三区四区第35| 久久婷婷青草| 亚洲在久久综合| www.av在线官网国产| 免费av不卡在线播放| 免费观看av网站的网址| 一级,二级,三级黄色视频| 国产黄色免费在线视频| 亚洲激情五月婷婷啪啪| 日韩中字成人| 国产精品秋霞免费鲁丝片| 日韩中文字幕视频在线看片| 日韩三级伦理在线观看| 亚洲综合色网址| 久久青草综合色| 啦啦啦中文免费视频观看日本| 校园人妻丝袜中文字幕| 啦啦啦视频在线资源免费观看| 9色porny在线观看| 久久97久久精品| 天堂8中文在线网| 久久精品久久精品一区二区三区| 亚洲高清免费不卡视频| 我要看黄色一级片免费的| 老司机亚洲免费影院| 青春草国产在线视频| 黄色 视频免费看| 欧美亚洲日本最大视频资源| 丝瓜视频免费看黄片| 日韩电影二区| 国产白丝娇喘喷水9色精品| 国产福利在线免费观看视频| 丝袜人妻中文字幕| 99香蕉大伊视频| 亚洲伊人久久精品综合| 国产精品免费大片| 国产成人免费无遮挡视频| 亚洲四区av| 热99国产精品久久久久久7| 丰满饥渴人妻一区二区三| 三上悠亚av全集在线观看| 亚洲图色成人| 久久国产精品男人的天堂亚洲 | 亚洲国产av影院在线观看| 国产片内射在线| 午夜91福利影院| 三级国产精品片| 国产在线一区二区三区精| 亚洲国产精品999| 美女脱内裤让男人舔精品视频| 成年人免费黄色播放视频| 在线精品无人区一区二区三| 欧美日韩一区二区视频在线观看视频在线| av在线观看视频网站免费| 欧美日韩精品成人综合77777| 欧美少妇被猛烈插入视频| 毛片一级片免费看久久久久| 黄网站色视频无遮挡免费观看| 夜夜爽夜夜爽视频| 国产成人精品在线电影| 国产免费一级a男人的天堂| 成人黄色视频免费在线看| 亚洲欧洲精品一区二区精品久久久 | 成人综合一区亚洲| 伦精品一区二区三区| 亚洲欧美一区二区三区黑人 | 亚洲熟女精品中文字幕| 免费av中文字幕在线| 男男h啪啪无遮挡| 97人妻天天添夜夜摸| 日日啪夜夜爽| 国产成人午夜福利电影在线观看| 国产免费现黄频在线看| 18在线观看网站| 亚洲一级一片aⅴ在线观看| av国产久精品久网站免费入址| 天堂俺去俺来也www色官网| 欧美精品高潮呻吟av久久| 久久久久久久久久成人| 在线观看免费高清a一片| 久久韩国三级中文字幕| 母亲3免费完整高清在线观看 | 秋霞伦理黄片| 伦精品一区二区三区| 天天操日日干夜夜撸| 一区二区av电影网| 久久久国产一区二区| 日本爱情动作片www.在线观看| 亚洲av免费高清在线观看| 少妇人妻精品综合一区二区| 婷婷色av中文字幕| 国产男女超爽视频在线观看| 肉色欧美久久久久久久蜜桃| 亚洲av欧美aⅴ国产| 免费高清在线观看视频在线观看| 精品国产乱码久久久久久小说| 久久人人爽人人爽人人片va| 综合色丁香网| av电影中文网址| 丰满饥渴人妻一区二区三| av在线app专区| 国产成人一区二区在线| 寂寞人妻少妇视频99o| 日本av免费视频播放| 免费播放大片免费观看视频在线观看| 有码 亚洲区| 国产精品三级大全| 丝袜美足系列| 国产老妇伦熟女老妇高清| 欧美亚洲日本最大视频资源| 午夜免费观看性视频| 免费看光身美女| 成年美女黄网站色视频大全免费| 波多野结衣一区麻豆| 在线观看国产h片| 久久精品国产综合久久久 | 欧美精品高潮呻吟av久久| 91精品三级在线观看| 香蕉丝袜av| 26uuu在线亚洲综合色| 免费大片18禁| 三级国产精品片| 91成人精品电影| 欧美老熟妇乱子伦牲交| 尾随美女入室| av在线观看视频网站免费| 草草在线视频免费看| 国产 一区精品| 亚洲av.av天堂| 亚洲欧洲精品一区二区精品久久久 | 日韩熟女老妇一区二区性免费视频| 欧美精品一区二区大全| 久久99蜜桃精品久久| 亚洲综合色网址| 天堂8中文在线网| 五月天丁香电影| 国产亚洲精品第一综合不卡 | 爱豆传媒免费全集在线观看| 两个人看的免费小视频| 亚洲av免费高清在线观看| 99热这里只有是精品在线观看| 国产精品久久久久久久久免| av黄色大香蕉| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人国语在线视频| 一边亲一边摸免费视频| 美女xxoo啪啪120秒动态图| 日韩欧美精品免费久久| 人妻少妇偷人精品九色| 91精品国产国语对白视频| 高清欧美精品videossex| 建设人人有责人人尽责人人享有的| 国产精品一二三区在线看| 又黄又爽又刺激的免费视频.| av有码第一页| 插逼视频在线观看| 国产精品一国产av| 国产亚洲欧美精品永久| 丰满少妇做爰视频| 最近的中文字幕免费完整| 亚洲欧美色中文字幕在线| 欧美日本中文国产一区发布| 国产熟女午夜一区二区三区| 99热国产这里只有精品6| 两性夫妻黄色片 | 久久精品久久精品一区二区三区| 狠狠精品人妻久久久久久综合| 寂寞人妻少妇视频99o| 久久久久精品性色| 夜夜骑夜夜射夜夜干| 欧美日韩国产mv在线观看视频| 久久久久国产精品人妻一区二区| 麻豆精品久久久久久蜜桃| 久久精品国产a三级三级三级| 黄色配什么色好看| 亚洲欧美精品自产自拍| 亚洲av成人精品一二三区| 亚洲欧洲精品一区二区精品久久久 | 日韩av免费高清视频| 亚洲国产精品一区三区| 最后的刺客免费高清国语| 性色av一级| 精品国产一区二区三区四区第35| 亚洲经典国产精华液单| 国产精品国产三级专区第一集| 色婷婷av一区二区三区视频| 2018国产大陆天天弄谢| 精品少妇内射三级| 国产69精品久久久久777片| 青春草视频在线免费观看| 亚洲综合色惰| av电影中文网址| 国产69精品久久久久777片| 夜夜爽夜夜爽视频| av有码第一页| 1024视频免费在线观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品一区三区| 女性被躁到高潮视频| 日韩一区二区三区影片| 成人午夜精彩视频在线观看| 九九在线视频观看精品| 91aial.com中文字幕在线观看| av又黄又爽大尺度在线免费看| 久久人人爽av亚洲精品天堂| 久久精品夜色国产| 成年女人在线观看亚洲视频| 久久久国产精品麻豆| 日本av免费视频播放| 久久青草综合色| 久久久久久人人人人人| 中国美白少妇内射xxxbb| 亚洲av电影在线进入| 又大又黄又爽视频免费| www.熟女人妻精品国产 | 久久久久久久久久人人人人人人| 精品熟女少妇av免费看| 久久久久精品久久久久真实原创| xxxhd国产人妻xxx| 婷婷成人精品国产| 宅男免费午夜| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 一边摸一边做爽爽视频免费| 精品一区二区三卡| 18禁国产床啪视频网站| 国产福利在线免费观看视频| 免费观看a级毛片全部| 久久综合国产亚洲精品| 亚洲国产欧美日韩在线播放| 伦精品一区二区三区| 国产一区二区三区av在线| 免费观看无遮挡的男女| 免费看av在线观看网站| 熟女av电影| 在线观看免费日韩欧美大片| 美女xxoo啪啪120秒动态图| 丝袜脚勾引网站| 久久精品熟女亚洲av麻豆精品| 99视频精品全部免费 在线| www.色视频.com| 在现免费观看毛片| 欧美精品高潮呻吟av久久| 波野结衣二区三区在线| 日韩大片免费观看网站| 18禁在线无遮挡免费观看视频| 狠狠精品人妻久久久久久综合| 国产伦理片在线播放av一区| 91国产中文字幕| 大片电影免费在线观看免费| 一区二区三区乱码不卡18| 久久久久久久国产电影| 人人妻人人澡人人爽人人夜夜| 看非洲黑人一级黄片| 亚洲色图综合在线观看| 久久久精品区二区三区| 精品国产一区二区三区四区第35| 免费高清在线观看日韩| 亚洲欧美成人综合另类久久久| 日韩av免费高清视频| 亚洲人与动物交配视频| 午夜视频国产福利| 香蕉国产在线看| 在线观看免费高清a一片| 精品人妻一区二区三区麻豆| 国产不卡av网站在线观看| 永久网站在线| 亚洲精品色激情综合| 免费观看性生交大片5| 热re99久久精品国产66热6| 亚洲,欧美精品.| 香蕉国产在线看| 国产在视频线精品| 午夜福利视频在线观看免费| 国产成人精品无人区| 国产精品女同一区二区软件| 国产 精品1| 精品人妻在线不人妻| 免费高清在线观看日韩| 在线观看三级黄色| av电影中文网址| 日日啪夜夜爽| 免费女性裸体啪啪无遮挡网站| 中文精品一卡2卡3卡4更新| 久久久久精品性色| 在线观看免费视频网站a站| 全区人妻精品视频| 欧美另类一区| 97精品久久久久久久久久精品| 蜜桃在线观看..| 日本vs欧美在线观看视频| 少妇猛男粗大的猛烈进出视频| 日韩制服丝袜自拍偷拍| 成人免费观看视频高清| 亚洲av男天堂| 亚洲欧美日韩另类电影网站| 色5月婷婷丁香| 亚洲国产av新网站| 五月天丁香电影| 国产午夜精品一二区理论片| 人人妻人人澡人人看| 丰满少妇做爰视频| 国产成人午夜福利电影在线观看| 五月伊人婷婷丁香| 欧美国产精品va在线观看不卡| a级片在线免费高清观看视频| 99热国产这里只有精品6| 日韩欧美精品免费久久| 久久久a久久爽久久v久久| freevideosex欧美| 久久97久久精品| 精品亚洲成a人片在线观看| 久久99一区二区三区| 久久午夜综合久久蜜桃| 欧美 日韩 精品 国产| 亚洲成色77777| 国产男女内射视频| 亚洲精品中文字幕在线视频| 国产成人精品福利久久| 一级a做视频免费观看| 日韩制服丝袜自拍偷拍| 成人毛片a级毛片在线播放| 国产成人av激情在线播放| 亚洲综合精品二区| 人妻 亚洲 视频| 少妇被粗大猛烈的视频| 国语对白做爰xxxⅹ性视频网站| 国产成人午夜福利电影在线观看| 婷婷色综合大香蕉| 久久国内精品自在自线图片| 亚洲欧美色中文字幕在线| 免费少妇av软件| 亚洲精品乱码久久久久久按摩| 中文字幕av电影在线播放| 欧美变态另类bdsm刘玥| 制服诱惑二区| 国产精品嫩草影院av在线观看| 亚洲精品美女久久av网站| 五月玫瑰六月丁香| 国产亚洲一区二区精品| 中文字幕人妻丝袜制服| 免费av中文字幕在线| 如何舔出高潮| 人成视频在线观看免费观看| 久久影院123| 高清毛片免费看| 久久久久精品人妻al黑| 最近最新中文字幕免费大全7| 妹子高潮喷水视频| 成年美女黄网站色视频大全免费| 男女国产视频网站| 在线观看www视频免费| 亚洲人成77777在线视频| 在线观看国产h片| 伊人久久国产一区二区| 91午夜精品亚洲一区二区三区| 黄色怎么调成土黄色| 国产日韩欧美亚洲二区| 大话2 男鬼变身卡| 国产1区2区3区精品| 女人被躁到高潮嗷嗷叫费观| 国产精品久久久久久久电影| 亚洲情色 制服丝袜| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 丰满乱子伦码专区| 精品99又大又爽又粗少妇毛片| 国产熟女欧美一区二区| 黑人高潮一二区| 女的被弄到高潮叫床怎么办| 国产女主播在线喷水免费视频网站| 亚洲内射少妇av| 久久99热6这里只有精品| 国产精品.久久久| 免费久久久久久久精品成人欧美视频 | 精品午夜福利在线看| 亚洲av综合色区一区| 男女啪啪激烈高潮av片| 91精品伊人久久大香线蕉| 国产在线一区二区三区精| 国产精品偷伦视频观看了| 精品卡一卡二卡四卡免费| 国产午夜精品一二区理论片| 又粗又硬又长又爽又黄的视频| 天天躁夜夜躁狠狠躁躁| 超色免费av| 婷婷色综合www| 人妻一区二区av| 国产精品三级大全| 黑人巨大精品欧美一区二区蜜桃 | 欧美精品av麻豆av| 熟妇人妻不卡中文字幕| 亚洲av日韩在线播放| 亚洲国产最新在线播放| 久久久久久久久久久免费av| 搡老乐熟女国产| 日韩视频在线欧美| 亚洲美女视频黄频| a级毛色黄片| 日韩av在线免费看完整版不卡| 精品卡一卡二卡四卡免费| 两个人免费观看高清视频| 精品一区二区三区四区五区乱码 | 青春草国产在线视频| 最近最新中文字幕大全免费视频 | 夜夜爽夜夜爽视频| 边亲边吃奶的免费视频| 99热国产这里只有精品6| 久久久久国产精品人妻一区二区| 久久人妻熟女aⅴ| 精品少妇内射三级| 97超碰精品成人国产| 最新中文字幕久久久久| 侵犯人妻中文字幕一二三四区| 亚洲成人av在线免费| 日日摸夜夜添夜夜爱| 蜜桃在线观看..| 少妇人妻久久综合中文| 欧美xxxx性猛交bbbb| 男女免费视频国产| 免费观看性生交大片5| 久久婷婷青草| 一级毛片 在线播放| 女的被弄到高潮叫床怎么办| 人体艺术视频欧美日本| 亚洲情色 制服丝袜| 精品人妻在线不人妻| 亚洲av福利一区| 国产成人欧美| 寂寞人妻少妇视频99o| 久久精品国产综合久久久 | 欧美日韩成人在线一区二区| 精品一区二区三区视频在线| 欧美老熟妇乱子伦牲交| 国产国拍精品亚洲av在线观看| 一个人免费看片子| 国产一区二区在线观看日韩| 秋霞在线观看毛片| 国产极品天堂在线| 亚洲高清免费不卡视频| 99re6热这里在线精品视频| 国产片内射在线| 亚洲欧美色中文字幕在线| 国产熟女午夜一区二区三区| 久久久精品区二区三区| 久久99蜜桃精品久久| 久久久精品区二区三区| 满18在线观看网站| 美女xxoo啪啪120秒动态图| 免费在线观看完整版高清| 男女高潮啪啪啪动态图| 久久久欧美国产精品| 2022亚洲国产成人精品| 女性被躁到高潮视频| 成人亚洲欧美一区二区av| 久久精品国产综合久久久 | 国产精品无大码| 亚洲av综合色区一区| 1024视频免费在线观看| 麻豆精品久久久久久蜜桃| 美女大奶头黄色视频| 王馨瑶露胸无遮挡在线观看| 久久久久久久久久久免费av| 视频中文字幕在线观看| 日本黄色日本黄色录像| 成人影院久久| 久久精品aⅴ一区二区三区四区 | 国产成人欧美| 日本av免费视频播放| 欧美亚洲日本最大视频资源| 免费人妻精品一区二区三区视频| 色婷婷久久久亚洲欧美| 亚洲一码二码三码区别大吗| 春色校园在线视频观看| 欧美变态另类bdsm刘玥| 中文精品一卡2卡3卡4更新| 欧美精品高潮呻吟av久久| 中文精品一卡2卡3卡4更新| av在线老鸭窝| 午夜av观看不卡| 免费看av在线观看网站| 日韩精品有码人妻一区| 晚上一个人看的免费电影| 精品久久久久久电影网| 老熟女久久久| 国产欧美另类精品又又久久亚洲欧美| 精品人妻在线不人妻|