• <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ù)處理
    极品教师在线视频| 波多野结衣高清无吗| 亚洲综合色惰| 男女视频在线观看网站免费| 成年女人毛片免费观看观看9| 免费在线观看日本一区| 午夜福利在线观看免费完整高清在 | 国产精品自产拍在线观看55亚洲| 成人av在线播放网站| 一本精品99久久精品77| 日本黄色视频三级网站网址| 91麻豆精品激情在线观看国产| 91在线精品国自产拍蜜月| 久久99热6这里只有精品| 免费一级毛片在线播放高清视频| 中文字幕av成人在线电影| 天堂动漫精品| 麻豆成人午夜福利视频| 国产探花在线观看一区二区| ponron亚洲| 亚州av有码| 日韩av在线大香蕉| 精品日产1卡2卡| 伦精品一区二区三区| 成人国产一区最新在线观看| 亚洲最大成人av| 乱系列少妇在线播放| 免费观看精品视频网站| 91午夜精品亚洲一区二区三区 | 中国美白少妇内射xxxbb| 欧美一区二区精品小视频在线| 天堂动漫精品| 午夜激情欧美在线| 亚洲国产日韩欧美精品在线观看| 日韩欧美精品v在线| 最近视频中文字幕2019在线8| 可以在线观看的亚洲视频| 中文字幕久久专区| 国产三级中文精品| 变态另类丝袜制服| 久99久视频精品免费| 亚洲三级黄色毛片| 我的老师免费观看完整版| 亚洲av免费在线观看| 美女cb高潮喷水在线观看| 精品久久久久久久久久久久久| 国语自产精品视频在线第100页| 国产一区二区在线观看日韩| 99久久无色码亚洲精品果冻| 男女那种视频在线观看| 国产高清视频在线播放一区| 九九热线精品视视频播放| 国产精品美女特级片免费视频播放器| 在线国产一区二区在线| 午夜福利18| 极品教师在线免费播放| 无人区码免费观看不卡| 在线播放无遮挡| 亚洲精品在线观看二区| 亚洲熟妇中文字幕五十中出| 日韩 亚洲 欧美在线| 国产激情偷乱视频一区二区| 日日啪夜夜撸| 亚洲av熟女| 成熟少妇高潮喷水视频| 熟女电影av网| 国产精品,欧美在线| 成人无遮挡网站| 老司机福利观看| 日韩,欧美,国产一区二区三区 | 老女人水多毛片| 免费观看在线日韩| 亚洲自偷自拍三级| 中文字幕精品亚洲无线码一区| 看片在线看免费视频| ponron亚洲| 九九热线精品视视频播放| 韩国av在线不卡| 久久久久免费精品人妻一区二区| 少妇丰满av| 免费人成视频x8x8入口观看| 精品久久久噜噜| 精品一区二区三区人妻视频| 国产精品日韩av在线免费观看| 国产人妻一区二区三区在| 啦啦啦观看免费观看视频高清| 国产精品国产高清国产av| 亚洲国产高清在线一区二区三| 波多野结衣高清无吗| 国产精品久久久久久精品电影| 久久久久久大精品| 精品午夜福利视频在线观看一区| 国产精品一区二区三区四区免费观看 | 草草在线视频免费看| 天天躁日日操中文字幕| 久久精品国产亚洲av天美| 小说图片视频综合网站| 亚洲av成人av| 91精品国产九色| 少妇高潮的动态图| 88av欧美| 可以在线观看的亚洲视频| 日韩在线高清观看一区二区三区 | 国产精品美女特级片免费视频播放器| 日韩,欧美,国产一区二区三区 | 国产在线男女| 亚洲乱码一区二区免费版| 岛国在线免费视频观看| 一本精品99久久精品77| 大型黄色视频在线免费观看| 国产精品精品国产色婷婷| 成年女人毛片免费观看观看9| 精品福利观看| 亚洲内射少妇av| 精品无人区乱码1区二区| 亚洲性久久影院| 午夜福利欧美成人| 久久久色成人| 国产精品乱码一区二三区的特点| 免费不卡的大黄色大毛片视频在线观看 | 很黄的视频免费| 欧美人与善性xxx| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 国产伦精品一区二区三区四那| 日日啪夜夜撸| 亚洲av熟女| 国产综合懂色| 成人国产一区最新在线观看| 日本 欧美在线| 日本一本二区三区精品| 国产精品1区2区在线观看.| 国产精品福利在线免费观看| 国产一区二区亚洲精品在线观看| 麻豆av噜噜一区二区三区| 免费无遮挡裸体视频| 国产精品,欧美在线| 久99久视频精品免费| 日日摸夜夜添夜夜添av毛片 | 国产色爽女视频免费观看| 22中文网久久字幕| 亚洲国产精品成人综合色| 级片在线观看| 亚洲av中文av极速乱 | 身体一侧抽搐| а√天堂www在线а√下载| 亚洲不卡免费看| 国产极品精品免费视频能看的| 国产一区二区亚洲精品在线观看| 免费黄网站久久成人精品| 成人特级黄色片久久久久久久| 哪里可以看免费的av片| 欧美极品一区二区三区四区| 国内少妇人妻偷人精品xxx网站| 极品教师在线免费播放| 十八禁国产超污无遮挡网站| 精品人妻1区二区| av视频在线观看入口| 午夜久久久久精精品| 1000部很黄的大片| 亚洲国产精品成人综合色| 亚洲欧美日韩高清在线视频| 男人舔奶头视频| 精品免费久久久久久久清纯| 亚洲在线观看片| av天堂在线播放| 99国产精品一区二区蜜桃av| 真实男女啪啪啪动态图| 亚洲专区中文字幕在线| 亚洲精品久久国产高清桃花| 91狼人影院| 国产91精品成人一区二区三区| 乱系列少妇在线播放| 亚洲欧美日韩东京热| 国产精品免费一区二区三区在线| 成人二区视频| 国产乱人视频| 高清毛片免费观看视频网站| 中国美白少妇内射xxxbb| 日韩欧美一区二区三区在线观看| 午夜精品久久久久久毛片777| 国产成人一区二区在线| 长腿黑丝高跟| 欧美国产日韩亚洲一区| av在线观看视频网站免费| 性插视频无遮挡在线免费观看| 免费观看精品视频网站| 日韩欧美一区二区三区在线观看| 国产成人av教育| 国产毛片a区久久久久| 人妻丰满熟妇av一区二区三区| 午夜福利在线在线| 国模一区二区三区四区视频| 国产探花极品一区二区| 亚洲av日韩精品久久久久久密| 欧美成人免费av一区二区三区| 不卡一级毛片| 午夜亚洲福利在线播放| 国产探花在线观看一区二区| 久久精品国产自在天天线| 春色校园在线视频观看| 黄色配什么色好看| 色哟哟·www| 99热这里只有精品一区| 亚洲精品一区av在线观看| 最后的刺客免费高清国语| 久久久久久久久中文| 精品一区二区免费观看| 美女 人体艺术 gogo| 亚洲专区国产一区二区| 国产伦人伦偷精品视频| 免费一级毛片在线播放高清视频| 国产精品久久久久久精品电影| 热99re8久久精品国产| 99久久九九国产精品国产免费| 有码 亚洲区| 内地一区二区视频在线| 精品人妻1区二区| 男插女下体视频免费在线播放| 国产伦一二天堂av在线观看| 99热只有精品国产| 色av中文字幕| 久久人人精品亚洲av| 亚洲最大成人手机在线| h日本视频在线播放| 国产精品日韩av在线免费观看| 日本三级黄在线观看| 日韩人妻高清精品专区| 性色avwww在线观看| 欧美日韩国产亚洲二区| 精品久久久久久久久av| 九色国产91popny在线| 国产一区二区三区av在线 | 国产高清有码在线观看视频| 国产精品不卡视频一区二区| 十八禁网站免费在线| 校园春色视频在线观看| 精品久久久久久,| 最后的刺客免费高清国语| 久久久久久久久久成人| 亚洲18禁久久av| 成人二区视频| 亚洲真实伦在线观看| 亚洲欧美清纯卡通| 免费av不卡在线播放| a级毛片a级免费在线| 国产aⅴ精品一区二区三区波| 久久国产乱子免费精品| 国内精品久久久久久久电影| ponron亚洲| 国模一区二区三区四区视频| 亚洲色图av天堂| 听说在线观看完整版免费高清| 久久久精品大字幕| 热99在线观看视频| 国产精品伦人一区二区| 日韩人妻高清精品专区| 别揉我奶头~嗯~啊~动态视频| 黄色日韩在线| 欧美在线一区亚洲| 亚洲欧美日韩东京热| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久久久久久丰满 | 午夜激情福利司机影院| 免费观看在线日韩| 高清毛片免费观看视频网站| 久久久久国产精品人妻aⅴ院| 欧美日韩中文字幕国产精品一区二区三区| 亚洲在线自拍视频| 久久这里只有精品中国| 校园春色视频在线观看| www日本黄色视频网| 亚洲人成伊人成综合网2020| 国内少妇人妻偷人精品xxx网站| 国产成人影院久久av| 黄色配什么色好看| 91在线精品国自产拍蜜月| 18+在线观看网站| 88av欧美| 99热这里只有是精品50| 日日撸夜夜添| 日本在线视频免费播放| 久久国产乱子免费精品| 99热这里只有是精品在线观看| 我的老师免费观看完整版| 欧美xxxx黑人xx丫x性爽| 伦精品一区二区三区| 成年女人永久免费观看视频| 啦啦啦观看免费观看视频高清| 男女做爰动态图高潮gif福利片| xxxwww97欧美| 波多野结衣高清无吗| 真人一进一出gif抽搐免费| 在线a可以看的网站| 久久6这里有精品| 亚洲精华国产精华液的使用体验 | 国产91精品成人一区二区三区| 婷婷丁香在线五月| 美女免费视频网站| 亚洲欧美激情综合另类| 我的老师免费观看完整版| 最后的刺客免费高清国语| 国产精品一区二区三区四区免费观看 | 亚洲黑人精品在线| 熟妇人妻久久中文字幕3abv| 日本免费a在线| 18+在线观看网站| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 国产精品免费一区二区三区在线| 日韩欧美国产在线观看| 久久精品国产清高在天天线| 国内久久婷婷六月综合欲色啪| 3wmmmm亚洲av在线观看| 国产精品久久久久久久久免| 少妇人妻一区二区三区视频| 国产精品久久久久久久久免| www日本黄色视频网| 国产69精品久久久久777片| 一区二区三区高清视频在线| 能在线免费观看的黄片| 日韩一本色道免费dvd| 真人一进一出gif抽搐免费| 在线观看免费视频日本深夜| 俺也久久电影网| 在线观看66精品国产| 欧美成人性av电影在线观看| 亚洲va日本ⅴa欧美va伊人久久| a级毛片a级免费在线| 日日摸夜夜添夜夜添小说| 国产精品,欧美在线| 久久久久久久久久久丰满 | 干丝袜人妻中文字幕| 午夜福利高清视频| 国产精品美女特级片免费视频播放器| 日韩一本色道免费dvd| 亚洲国产色片| 永久网站在线| 国产人妻一区二区三区在| 日本-黄色视频高清免费观看| www.www免费av| 国产精品1区2区在线观看.| 亚洲aⅴ乱码一区二区在线播放| 国产伦在线观看视频一区| 狂野欧美激情性xxxx在线观看| 真实男女啪啪啪动态图| 国产男靠女视频免费网站| 99国产精品一区二区蜜桃av| 久久久精品大字幕| 国产真实伦视频高清在线观看 | 九九在线视频观看精品| 俺也久久电影网| 国产精品亚洲一级av第二区| 禁无遮挡网站| 欧美日韩黄片免| 一本久久中文字幕| 亚洲一区二区三区色噜噜| 午夜精品在线福利| 久久精品91蜜桃| 给我免费播放毛片高清在线观看| 永久网站在线| 精品人妻熟女av久视频| 99久久久亚洲精品蜜臀av| 一本久久中文字幕| 我的老师免费观看完整版| 欧美色视频一区免费| 精品一区二区三区视频在线观看免费| 日本欧美国产在线视频| 国产高潮美女av| 能在线免费观看的黄片| 男人舔奶头视频| 很黄的视频免费| 又粗又爽又猛毛片免费看| 真人一进一出gif抽搐免费| 国产一区二区三区视频了| 亚洲最大成人手机在线| 国内少妇人妻偷人精品xxx网站| 国产淫片久久久久久久久| 成人毛片a级毛片在线播放| 1024手机看黄色片| 久久精品人妻少妇| 99热这里只有是精品在线观看| 午夜免费成人在线视频| 国产v大片淫在线免费观看| www日本黄色视频网| 欧美日韩亚洲国产一区二区在线观看| 国内揄拍国产精品人妻在线| 国产又黄又爽又无遮挡在线| 精品午夜福利在线看| 联通29元200g的流量卡| 国产一区二区亚洲精品在线观看| 在线国产一区二区在线| 校园人妻丝袜中文字幕| 一个人看的www免费观看视频| 一区二区三区免费毛片| 日韩中文字幕欧美一区二区| 一本一本综合久久| 男人舔奶头视频| 亚洲成a人片在线一区二区| 国产av在哪里看| 桃色一区二区三区在线观看| 国产不卡一卡二| 亚洲图色成人| 免费av毛片视频| 观看免费一级毛片| 少妇的逼好多水| 免费一级毛片在线播放高清视频| 欧美一区二区亚洲| 精品福利观看| 久久九九热精品免费| 成人av一区二区三区在线看| 97超级碰碰碰精品色视频在线观看| 日本与韩国留学比较| 免费看av在线观看网站| 精品国产三级普通话版| 日韩欧美国产在线观看| 变态另类丝袜制服| 亚洲性久久影院| 久久久久久久精品吃奶| 一区二区三区激情视频| 成人性生交大片免费视频hd| 国产av在哪里看| 91麻豆av在线| 丰满的人妻完整版| 亚洲欧美日韩高清在线视频| or卡值多少钱| 色综合色国产| 亚洲四区av| 精品人妻偷拍中文字幕| 韩国av在线不卡| 97人妻精品一区二区三区麻豆| 亚洲av不卡在线观看| 免费av毛片视频| 大又大粗又爽又黄少妇毛片口| 久久久久久大精品| 黄色一级大片看看| 18禁裸乳无遮挡免费网站照片| 国产免费一级a男人的天堂| 一级黄片播放器| 国产一区二区三区av在线 | 久久草成人影院| 国产不卡一卡二| 一本精品99久久精品77| 亚洲一级一片aⅴ在线观看| 免费黄网站久久成人精品| 国产一区二区在线观看日韩| 国产高清三级在线| 国产精品人妻久久久影院| 99久久成人亚洲精品观看| av专区在线播放| 国产精品日韩av在线免费观看| 欧美日韩国产亚洲二区| 久久午夜福利片| 免费在线观看影片大全网站| 成人鲁丝片一二三区免费| 久久亚洲精品不卡| 国产探花在线观看一区二区| 亚洲av一区综合| 国产精品98久久久久久宅男小说| 久9热在线精品视频| 人人妻人人澡欧美一区二区| 欧美国产日韩亚洲一区| 精品久久久久久久久av| 伊人久久精品亚洲午夜| 99久久精品国产国产毛片| 久久精品人妻少妇| 欧美另类亚洲清纯唯美| 免费高清视频大片| 老司机深夜福利视频在线观看| 精品久久久久久久末码| 久久热精品热| 亚洲av电影不卡..在线观看| 亚洲精品色激情综合| 麻豆一二三区av精品| 搡老妇女老女人老熟妇| 亚洲人成网站在线播放欧美日韩| 午夜激情欧美在线| 精品国产三级普通话版| 可以在线观看毛片的网站| 欧美日韩乱码在线| 亚洲中文字幕日韩| 97超级碰碰碰精品色视频在线观看| 精品国产三级普通话版| 美女被艹到高潮喷水动态| 欧美日韩乱码在线| 日本黄色视频三级网站网址| 久久99热6这里只有精品| 草草在线视频免费看| 精品国内亚洲2022精品成人| 国产黄色小视频在线观看| 三级男女做爰猛烈吃奶摸视频| 日韩亚洲欧美综合| 村上凉子中文字幕在线| 高清在线国产一区| 亚洲精品国产成人久久av| 欧美不卡视频在线免费观看| 在线免费十八禁| 精华霜和精华液先用哪个| 又黄又爽又刺激的免费视频.| 噜噜噜噜噜久久久久久91| 欧美日韩瑟瑟在线播放| 1024手机看黄色片| 精华霜和精华液先用哪个| 精品一区二区三区av网在线观看| 听说在线观看完整版免费高清| 波野结衣二区三区在线| 三级毛片av免费| 18禁黄网站禁片免费观看直播| av黄色大香蕉| 亚洲性夜色夜夜综合| 一级黄色大片毛片| 亚洲欧美激情综合另类| 日日撸夜夜添| 久久久久免费精品人妻一区二区| 俺也久久电影网| 人人妻,人人澡人人爽秒播| 精品99又大又爽又粗少妇毛片 | 国内精品一区二区在线观看| 精品一区二区三区视频在线| 黄片wwwwww| 啦啦啦韩国在线观看视频| 日日干狠狠操夜夜爽| 人妻夜夜爽99麻豆av| 如何舔出高潮| 精品久久久久久久久av| 精品福利观看| 午夜激情福利司机影院| 久久久久久伊人网av| 日韩欧美在线乱码| 99riav亚洲国产免费| 国产av在哪里看| 乱码一卡2卡4卡精品| 99九九线精品视频在线观看视频| 亚洲国产高清在线一区二区三| 国产69精品久久久久777片| 97超视频在线观看视频| 国产免费男女视频| 亚洲欧美日韩高清专用| 欧美bdsm另类| 精品久久久久久久人妻蜜臀av| 欧美成人a在线观看| 日韩精品青青久久久久久| 九九爱精品视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 黄色日韩在线| 麻豆精品久久久久久蜜桃| 国产私拍福利视频在线观看| 国产精品不卡视频一区二区| 春色校园在线视频观看| 久久国产乱子免费精品| 免费人成在线观看视频色| 免费搜索国产男女视频| 国产精品久久久久久亚洲av鲁大| 别揉我奶头 嗯啊视频| 国产91精品成人一区二区三区| 久久久精品大字幕| 色综合站精品国产| 日韩国内少妇激情av| 久久久久九九精品影院| 国产一区二区三区av在线 | 久久九九热精品免费| 亚洲av中文av极速乱 | 狂野欧美激情性xxxx在线观看| 精品不卡国产一区二区三区| 亚洲一区高清亚洲精品| 亚洲精品色激情综合| 亚洲国产精品合色在线| 狂野欧美白嫩少妇大欣赏| 香蕉av资源在线| 在线观看舔阴道视频| 亚洲最大成人手机在线| 国产不卡一卡二| 人人妻人人看人人澡| 男人的好看免费观看在线视频| 国产精品福利在线免费观看| 成人综合一区亚洲| 天天躁日日操中文字幕| 免费看a级黄色片| 淫妇啪啪啪对白视频| 观看美女的网站| 老熟妇仑乱视频hdxx| 别揉我奶头 嗯啊视频| 深爱激情五月婷婷| 12—13女人毛片做爰片一| 别揉我奶头 嗯啊视频| 美女xxoo啪啪120秒动态图| 三级国产精品欧美在线观看| 亚洲av电影不卡..在线观看| 国产精品日韩av在线免费观看| 99视频精品全部免费 在线| 91午夜精品亚洲一区二区三区 | or卡值多少钱| 熟女电影av网| 亚洲av不卡在线观看| 日韩强制内射视频| 999久久久精品免费观看国产| 一卡2卡三卡四卡精品乱码亚洲| 日韩强制内射视频| 国产一区二区三区视频了| 欧美高清性xxxxhd video| av在线天堂中文字幕| 两个人的视频大全免费| 欧洲精品卡2卡3卡4卡5卡区| 美女xxoo啪啪120秒动态图| 国产高潮美女av| 一本精品99久久精品77| 岛国在线免费视频观看| 人人妻,人人澡人人爽秒播| 99热这里只有精品一区| 久久精品影院6| 伦精品一区二区三区| 精品人妻偷拍中文字幕| 国产真实乱freesex| 国产av在哪里看| 国产v大片淫在线免费观看|