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

    三維點擴散函數(shù)層距與圖像復(fù)原的關(guān)系

    2016-07-28 06:51:57梁日柳莫春球
    關(guān)鍵詞:采樣圖像復(fù)原

    梁日柳,陳 華,莫春球,聶 雄

    (廣西大學(xué)計算機與電子信息學(xué)院, 廣西南寧530004)

    ?

    三維點擴散函數(shù)層距與圖像復(fù)原的關(guān)系

    梁日柳,陳華,莫春球,聶雄

    (廣西大學(xué)計算機與電子信息學(xué)院, 廣西南寧530004)

    摘要:為研究三維點擴散函數(shù)層距與圖像復(fù)原的關(guān)系,對數(shù)字共焦光學(xué)顯微成像系統(tǒng)3D-PSF的構(gòu)成、3D-PSF不同層距與圖像復(fù)原效果的關(guān)系進行理論分析。通過仿真實驗,研究了不同放大倍數(shù)的3D-PSF,在不同層距,相同層數(shù)和徑向大小的情況下與圖像復(fù)原效果的關(guān)系。實驗表明:層距小于0.6 μm時,放大倍數(shù)為40倍,圖像在層距為0.5 μm處取得最佳復(fù)原效果,放大倍數(shù)為20倍和10倍時,圖像在層距為0.025 μm處取得較好的復(fù)原效果。研究結(jié)果為在細胞切片的采集和三維顯微圖像復(fù)原實際應(yīng)用中對3D-PSF的層距選取提供了理論依據(jù)和選取方法。

    關(guān)鍵詞:三維點擴散函數(shù);圖像復(fù)原;層距;采樣;改善信噪比

    0引言

    由于可以在傳統(tǒng)的熒光顯微鏡基礎(chǔ)上,利用數(shù)學(xué)算法去除圖像中的離焦模糊成分,提高圖像分辨率,數(shù)字共焦顯微技術(shù)近年來迅速成為研究熱點,以該技術(shù)發(fā)展的數(shù)字共焦顯微儀在遺傳學(xué)、醫(yī)學(xué)、生物植物學(xué)等領(lǐng)域具有廣泛的應(yīng)用前景[1]。在該技術(shù)中,沿著z軸方向(即光軸或者軸向)移動生物細胞樣本通過生物光學(xué)顯微鏡物鏡焦平面,獲取序列二維切片圖像,然后使用三維顯微圖像復(fù)原算法進行處理[2],獲取清晰的二維圖像序列和三維圖像[3-6]。

    在三維顯微圖像處理中,三維點擴散函數(shù)(three dimensional point spread function, 3D-PSF)的選取與圖像復(fù)原效果密切相關(guān)。研究表明,3D-PSF層數(shù)越多,復(fù)原效果越好,但同時處理的時間也越長。當(dāng)3D-PSF層數(shù)一定時,3D-PSF層距的選取與圖像復(fù)原效果存在何種關(guān)系,一直是有待研究解決的問題。

    本文在三種不同的放大倍數(shù)下,通過仿真實驗,使用四種不同層數(shù)的細胞樣本切片,研究不同層距,相同層數(shù)和徑向大小的3D-PSF與圖像復(fù)原的關(guān)系,為更好地獲取高分辨率圖像提供理論依據(jù)。

    1光學(xué)顯微成像系統(tǒng)3D-PSF

    1.1三維點擴散函數(shù)結(jié)構(gòu)

    圓形孔徑小散焦光學(xué)傳遞函數(shù)OTF的數(shù)學(xué)公式[7-9]如下:

    (1)

    其中w為散焦誤差,q為頻率,fc為系統(tǒng)截止頻率。計算不同散焦量的光學(xué)傳遞函數(shù),可以獲得相應(yīng)散焦量的2D-PSF,進而構(gòu)成雙錐體型3D-PSF[10]。

    3D-PSF的空間大小包括軸向大小z和徑向大小(x-y),其在離散空域里可視為一個三維矩陣。沿著3D-PSF雙錐體中心軸(z軸)不同的橫截面(x-y面),對應(yīng)著一系列不同散焦量的徑向二維點擴散函數(shù)(two dimensional point spread function, 2D-PSF),其中z=0處的中間截面為焦面2D-PSF。

    研究表明,3D-PSF的能量分布與放大倍數(shù)密切相關(guān)。三維點擴散函數(shù)的能量主要集中在雙錐體中部的錐頂附近區(qū)域,放大倍數(shù)越大,在z軸方向,能量擴散程度越小,遠離焦面z=0時能量迅速衰減,在徑向(即x-y面)方向上,能量擴散程度越大,衰減速度緩慢,隨著放大倍數(shù)減小,三維點擴散函數(shù)的能量在z軸方向擴散程度增大,衰減速度變緩,在徑向方向上,三維點擴散函數(shù)的能量擴散程度小,主要分布在z軸周圍,遠離z軸時能量迅速衰減。

    1.2不同層距3D-PSF與圖像復(fù)原關(guān)系分析

    圖像復(fù)原處理采用的3D-PSF空間大小決定著復(fù)原效果和復(fù)原時間,3D-PSF空間大小選取越大,復(fù)原效果越好,同時處理的時間也越長。研究表明,3D-PSF的能量主要集中雙錐體中部區(qū)域,因此在進行圖像復(fù)原時,通常舍去3D-PSF周圍大部分能量稀少的空間區(qū)域,而選取3D-PSF能量較為集中的中部區(qū)域。

    將選取的3D-PSF以錐頂為中心,沿著z軸兩端等間隔離散采樣,獲得一系列的2D-PSF構(gòu)成新的3D-PSF,并對其進行適當(dāng)修正,以免破壞3D-PSF的整體性。其中,每個3D-PSF僅層距不同,層數(shù)和徑向大小相同。

    可以預(yù)測,3D-PSF的層距太小或者太大都無法得到最佳復(fù)原效果,復(fù)原效果與層距的關(guān)系曲線呈現(xiàn)倒U的形狀。即在相同層數(shù)和徑向大小的3D-PSF下,圖像復(fù)原效果隨著3D-PSF層距增大而提高,當(dāng)3D-PSF層距增大到臨界值時,圖像復(fù)原效果隨著3D-PSF層距增大而下降。這是因為,當(dāng)3D-PSF的層數(shù)固定時,層距太小使得3D-PSF的空間區(qū)域過小,只包含靠近焦面2D-PSF兩側(cè)的2D-PSF,無法有效地去除焦面的鄰近平面造成的散焦信息,取得最佳復(fù)原效果,層距太大使得3D-PSF的空間區(qū)域過大,3D-PSF兩端的部分2D-PSF遠離焦面2D-PSF,無法有效還原焦面的細節(jié),復(fù)原效果作用不大,甚至引起失真。

    根據(jù)上述的分析,本文在三種不同的放大倍數(shù)下,使用不同層數(shù)的細胞樣本切片,研究不同層距,相同層數(shù)和徑向大小的3D-PSF與圖像復(fù)原效果的關(guān)系。

    2實驗與分析

    2.1制作3D-PSF

    設(shè)顯微鏡物鏡參數(shù):放大倍數(shù)M=40,數(shù)值孔徑NA=0.65;顯微鏡機械鏡筒長度為160 mm;光源波長為550 nm;CCD參數(shù)為1/3英寸,像素值為640×480。

    2.1.1制作顯微鏡物鏡光學(xué)系統(tǒng)3D-PSF

    根據(jù)光學(xué)傳遞函數(shù)制作一個直徑為91(91×91)、層數(shù)為91、軸向采樣間隔Δz(以下簡稱層距)為0.025 μm的光學(xué)系統(tǒng)3D-PSF,命名為h10,其中h10的40指放大倍數(shù),大小為91×91×91。

    2.1.2制作用于復(fù)原的層距不同層數(shù)相同的3D-PSF

    對h10進行等間隔采樣,獲得徑向大小、層數(shù)均為15的3D-PSF(15×15×15),層距分別為Δz1=0.025 μm、Δz2=0.05 μm,…,Δz24=0.6 μm,分別命名為h10_0025、h10_005,…,h10_06,其中h10_0025中的40指放大倍數(shù),0025指層距0.025 μm,得到一系列不同層距的h15.15.15[11-12]。

    2.2制作三維仿真樣本

    圖1為清晰圖像、三種不同放大倍數(shù)下的模糊圖像和復(fù)原圖像。

    使用圖1(a)所示的清晰圖像(151×151,256灰度級)作為初始的二維樣本,通過旋轉(zhuǎn)方式制作出351幅相互關(guān)聯(lián)的二維圖像,并將這351幅圖像構(gòu)建為三維清晰圖像f(151×151×351),分別將f與h10卷積得到三維仿真模糊切片圖像g40,大小為151×151×351,層距為0.025 μm,圖1(b)是位于g40中間層的圖像。

    分別用層距為Δz1、Δz2,…,Δz24對g40進行等間隔軸向采樣,采樣層數(shù)均為13層,得到第一組不同層距、相同層數(shù)和徑向大小的三維模糊切片圖像,分別命名為g40_13_0025、g40_13_005,…,g40_13_06,其中g(shù)40_13_0025的40指放大倍數(shù),13指模糊切片層數(shù),0025指層距0.025 μm,大小均為151×151×13。

    同理,將采樣層數(shù)增加為19層、25層、31層時,得到得到另外三組模糊切片圖像,最大的層距分別為0.4 μm、0.35 μm和0.25 μm。

    2.3仿真實驗和結(jié)果分析

    2.3.1圖像復(fù)原

    分別將四組三維模糊切片圖像與相應(yīng)層距的h15.15.15進行去卷積復(fù)原處理。采用經(jīng)典的復(fù)原算法:最大似然法[13-14],迭代次數(shù)為500次,分別記錄復(fù)原時間,復(fù)原結(jié)果圖像分別命名為fi40_13_0025、fi40_13_005等,fi40_13_0025中的40指放大倍數(shù),13指模糊切片層數(shù)。0025指層距0.025 μm,圖1(c)為fi40_13_05的中間層復(fù)原圖像。

    同理,將顯微鏡物鏡參數(shù)依次修改為:①放大倍數(shù)M=20,數(shù)值孔徑NA=0.45。②放大倍數(shù)M=10,數(shù)值孔徑NA=0.25,其他條件不變再次進行仿真實驗。

    2.3.2不同層距的3D-PSF復(fù)原結(jié)果圖像的改善信噪比

    對得到的三維復(fù)原結(jié)果圖像使用公式(2)進行改善信噪比ISNR計算。

    (2)

    其中:f、g和fi分別是三維清晰圖像、三維模糊圖像和三維復(fù)原結(jié)果圖像(151×151)。ISNR值越大,表明圖像的復(fù)原效果越好。

    (a) f

    (b) g40

    (c) fi40_13_05

    (d) g20

    (e) fi20_13_0025

    (f) fi20_13_06

    (g) g10

    (h) fi10_13_0025

    (i) fi10_31_0175

    圖2為不同放大倍數(shù)、不同層數(shù)模糊切片圖像下,相同層數(shù)和徑向大小的3D-PSF層距與三維復(fù)原圖像的ISNR關(guān)系曲線。

    2.3.3結(jié)果分析

    ①對比圖2(a)~圖2(c)的曲線圖可以看出,使用相同層距的h15.15.15對圖像復(fù)原時,放大倍數(shù)越大,復(fù)原效果越好。同時,40倍時,ISNR曲線最先呈現(xiàn)出倒U的形狀。這是因為在3D-PSF能量分布中,相同空間大小的3D-PSF,放大倍數(shù)越大,它所包含的能量越多,復(fù)原效果越好。

    ②從圖2(a)可以看出,M=40倍時,ISNR曲線在Δz=0.1 μm處,得到最小值,隨著層距的增加,ISNR曲線上升,在層距為0.5 μm處取得最大值,復(fù)原圖像如圖1(c)所示,繼續(xù)增大層距,ISNR曲線下降,總體趨勢呈現(xiàn)一個倒U的形狀。這說明適當(dāng)增大3D-PSF的層距,可以提高圖像的改善信噪比。當(dāng)層距過小時,h15.15.15的空間區(qū)域過小,只包含靠近焦面2D-PSF兩側(cè)的2D-PSF,無法有效去除焦面附近平面造成的模糊效果,復(fù)原效果不理想。增大層距,h15.15.15的空間區(qū)域隨之增大,包含的2D-PSF攜帶更多的復(fù)原信息,圖像清晰度逐漸改善最終得到最佳復(fù)原效果。再增大層距時,構(gòu)成h15.15.15的部分2D-PSF遠離焦面2D-PSF,復(fù)原作用減小,復(fù)原效果開始下降。

    (a)M=40

    (b) M=20

    (c) M=10

    ③圖2(b)中的曲線可以看出,M=20倍時,使用相同的h15.15.15復(fù)原19層、25層和31層的模糊切片,ISNR曲線在Δz=0.025 μm處取得最大值,在Δz=0.3 μm附近取得最小值,隨著層距的增加,ISNR曲線呈上升趨勢。這是因為放大倍數(shù)越小,3D-PSF能量在z軸方向上擴散程度越大,h15.15.15需要更大的層距才能包含有用的2D-PSF用于復(fù)原焦面圖像,所以ISNR曲線逐漸上升,同時ISNR最大值將在較大的層距(Δz>0.5 μm)取得,層距繼續(xù)增大,h15.15.15包含的2D-PSF無法有效復(fù)原焦面圖像,ISNR曲線下降。

    ④在圖2(c)中,10倍的復(fù)原效果曲線總體呈下降趨勢,在Δz<0.6 μm時沒有出現(xiàn)明顯的上升趨勢,這是因為相比40倍和20倍,10倍的3D-PSF能量在z軸方向上擴散程度最大,能量更加分散。在層距較小時,采樣的2D-PSF過于相近,無法精確反映3D-PSF的能量分布。此時,將h15.15.15與相應(yīng)層距的三維切片模糊圖像進行處理時,復(fù)原效果惡化,ISNR曲線下降。依據(jù)M=10倍時3D-PSF的能量分布,隨著層距的增大,10倍的ISNR曲線將會上升,同時取得最大值的層距比20倍的層距更大。

    ⑤圖2(b)中13層模糊切片與其他不同切片層數(shù)的ISNR曲線不同,一直呈下降趨勢,這是因為,圖像的復(fù)原不僅與3D-PSF函數(shù)相關(guān),也與切片的采樣數(shù)量有關(guān)。如圖1(f)所示在Δz=0.6 μm處,由于切片采樣過少,3D-PSF無法有效還原鄰近切片對焦面引起的散焦信息,復(fù)原效果不好。所以合適增加切片的采樣,有助于模糊圖像的復(fù)原。同理,在圖2(c)中,由于切片采樣層數(shù)多,使用相同的3D-PSF復(fù)原時,25層、31層切片的中間層模糊圖像在Δz=0.2 μm附近取得較好的復(fù)原效果,如圖1(i)所示。

    研究結(jié)果表明,M=40時,3D-PSF層距為0.5 μm,圖像可以得到最佳的復(fù)原效果,M=20和M=10時,3D-PSF層距為0.025 μm,可以獲得較高的改善信噪比,如圖1(e)和圖1(h)。放大倍數(shù)為20倍、10倍時,ISNR曲線具體在多大的層距完整呈現(xiàn)出現(xiàn)倒U形,將在后續(xù)的實驗進一步研究。

    3結(jié)語

    為研究三維點擴散函數(shù)層距與圖像復(fù)原的關(guān)系,通過仿真實驗,在三種不同放大倍數(shù)下,使用不同層距,相同層數(shù)和徑向大小的3D-PSF與相應(yīng)層距的三維仿真模糊圖像去卷積,獲得3D-PSF層距與三維圖像復(fù)原效果的關(guān)系。實驗表明,層距小于0.6 μm時,放大倍數(shù)為40倍,圖像在層距為0.5 μm處取得最佳復(fù)原效果,放大倍數(shù)為20倍和10倍時,圖像在層距為0.025 μm處取得較好的復(fù)原效果。本文研究結(jié)果為生物細胞切片的采樣和三維顯微圖像復(fù)原實際應(yīng)用中對3D-PSF的層距選取提供了參考,進一步提高3D-PSF空間大小與采樣密度選取的科學(xué)性和確定性。

    參考文獻:

    [1]CONCHELLO J A, LICHTMAN J W.Optical sectioning microscopy[J]. Nature methods,2005,2(12):920-931.

    [2]WALLACEL W, SCHAEFER L H, SWEDLOW J R.A working person’s guide to deconvolution in light microscopy[J]. Bio Techniques,2001,31:1076-1097.

    [3]陳華,金偉其,蘇秉華,等.數(shù)字共焦顯微技術(shù)及其處理算法的進展[J]. 光電子技術(shù)與信息, 2005, 18(2):1-5.

    [4]袁景和,方曉紅.超分辨熒光顯微鏡-顯納鏡[J]. 化學(xué)通報, 2014, 77(11):1-4.

    [5]蘇秉華.超分辨率圖像復(fù)原及其進展[J]. 光學(xué)技術(shù), 2001, 26(5):1-3.

    [6]林廣升,陳華.基于LTC6090的橋式壓電物鏡驅(qū)動器驅(qū)動電源設(shè)計[J]. 廣西大學(xué)學(xué)報(自然科學(xué)版), 2015, 40(3):1-6.

    [7]唐玉科,何小海,陶青川,等.基于澤尼克多項式的顯微鏡點擴展函數(shù)研究[J]. 光學(xué)學(xué)報, 2009,29(1):169-175.

    [8]李蕊,陶青川,何小海,等.高斯型點擴散函數(shù)估計的最近鄰算法[J]. 光電工程, 2007,34(6):1-6.

    [9]MARKHAM J, CONEHELLE J A.Fast maximun likelihood image-restoration algorithms for three-dimensional fluorescence microscopy[J]. J opt soc am a,2001,18(5):1062-1072.

    [10]滕奇志,陶青川,趙佳,等.基于深度變化成像模型的圖像估計[J]. 計算機工程與應(yīng)用, 2005(4):32-34.

    [11]戴志華,王金剛,趙星,等.光學(xué)顯微集成成像三維圖像快速獲取技術(shù)[J]. 中國科學(xué):技術(shù)科學(xué), 2012,42(4):1-5.

    [12]聶雄,陳華.數(shù)字共焦顯微儀序列光學(xué)切片自動采集方法研究[J]. 儀器儀表學(xué)報, 2010,9(31):1-4.

    [13]PREZA C, CONCHELLO J A.Depth-variant maximun-likelihood restoration for three-dimensional fluorescence microscopy[J]. J opt soc am a,2004,21(9):1593-1601.

    [14]HOLMES T J.Maximum-likelihood image restoration adapted for non-coherent optical imaging[J]. Journal of the Optical Society of Anerica (A),1988,5(5):666-673.

    (責(zé)任編輯梁碧芬)

    收稿日期:2016-01-07;

    修訂日期:2016-03-01

    基金項目:國家自然科學(xué)基金資助項目(61164019);廣西科學(xué)基金資助項目(2012GXNSFAA053216)

    通訊作者:陳華(1958—),廣西桂平人,廣西大學(xué)教授;E-mail: chenhua_cn@163.com。

    doi:10.13624/j.cnki.issn.1001-7445.2016.0810

    中圖分類號:TP391

    文獻標(biāo)識碼:A

    文章編號:1001-7445(2016)03-0810-06

    On the relationship between the layer distance of three dimensional point spread function and image restoration

    LIANG Ri-liu, CHEN Hua, MO Chun-qiu, NIE Xiong

    (School of Computer, Electronics and Information, Guangxi University, Nanning 530004, China)

    Abstract:In order to study the relationship between the three dimensional point spread function and the image restoration, the structure of 3D-PSF and the relationship between the effect of different layer distance and image restoration of digital confocal optical microscope imaging system are analyzed in theory. The relationship between the 3D-PSF of different distance and same layer number and radial amount is studied in different magnifications. Simulation results show that the best image restoration effect is obtained at the layer distance of 0.5 μm when the distance is less than 0.6 μm and the magnification is 40 times. When the magnification is 20 times or 10 times, better image restoration is achieved at the layer distance of 0.025 μm. The research results provide a theoretical basis for the selection of 3D-PSF layer distance in the practical application of the collection of cell sections and the practical application of three dimensional microscopic image restoration.

    Key words:three dimensional point spread function(3D-PSF); image restoration; layer distance;sampling; improved signal to noise ratio (ISNR)

    引文格式: 梁日柳,陳 華,莫春球.三維點擴散函數(shù)層距與圖像復(fù)原的關(guān)系[J].廣西大學(xué)學(xué)報(自然科學(xué)版),2016,41(3):810-815.

    猜你喜歡
    采樣圖像復(fù)原
    基于MTF的實踐九號衛(wèi)星圖像復(fù)原方法研究
    談食品安全事故的調(diào)查取證
    祖國(2017年5期)2017-03-22 13:06:04
    一種基于顯著性邊緣的運動模糊圖像復(fù)原方法
    基層動物抗體監(jiān)測采樣注意事項
    農(nóng)業(yè)地質(zhì)調(diào)查中土壤樣品采集淺析
    一種煤炭機械化采制樣裝置的設(shè)計
    科技視界(2016年22期)2016-10-18 15:44:39
    鍋爐智能監(jiān)控儀表在參數(shù)自動檢測中的應(yīng)用
    一種基于采樣法的LCR測量儀的研制
    基于MTFC的遙感圖像復(fù)原方法
    模糊圖像復(fù)原的高階全變差正則化模型構(gòu)建
    天堂动漫精品| 看片在线看免费视频| 曰老女人黄片| 午夜免费成人在线视频| 18禁国产床啪视频网站| 在线观看免费日韩欧美大片| 国产在线精品亚洲第一网站| 国产精品乱码一区二三区的特点 | 久久久久久亚洲精品国产蜜桃av| 女人爽到高潮嗷嗷叫在线视频| 成人特级黄色片久久久久久久| 涩涩av久久男人的天堂| 国产亚洲欧美在线一区二区| 桃红色精品国产亚洲av| 最新美女视频免费是黄的| 精品福利观看| 色综合亚洲欧美另类图片| 高清黄色对白视频在线免费看| 黄色片一级片一级黄色片| 亚洲天堂国产精品一区在线| 欧美成人午夜精品| 丰满的人妻完整版| 久久天躁狠狠躁夜夜2o2o| 老司机午夜福利在线观看视频| 一级毛片精品| 国产精品综合久久久久久久免费 | 少妇的丰满在线观看| 国产成人免费无遮挡视频| 亚洲欧美激情在线| 久久亚洲精品不卡| 中文字幕另类日韩欧美亚洲嫩草| 国产麻豆成人av免费视频| 亚洲国产精品久久男人天堂| 手机成人av网站| 免费看a级黄色片| 啦啦啦 在线观看视频| 日韩大码丰满熟妇| av天堂久久9| 久久精品亚洲精品国产色婷小说| 亚洲精品一区av在线观看| 日本精品一区二区三区蜜桃| 欧美乱妇无乱码| 中文字幕色久视频| 老司机深夜福利视频在线观看| 国产精品影院久久| 天堂√8在线中文| 国产成人精品无人区| 日韩一卡2卡3卡4卡2021年| 国产激情欧美一区二区| 老熟妇仑乱视频hdxx| 国产精品一区二区三区四区久久 | 国内精品久久久久久久电影| 免费一级毛片在线播放高清视频 | 极品人妻少妇av视频| 欧美久久黑人一区二区| 色哟哟哟哟哟哟| 91麻豆精品激情在线观看国产| av网站免费在线观看视频| 亚洲免费av在线视频| 亚洲国产欧美网| 午夜精品国产一区二区电影| 日本精品一区二区三区蜜桃| 操出白浆在线播放| 欧美最黄视频在线播放免费| 成人av一区二区三区在线看| 女生性感内裤真人,穿戴方法视频| 好男人在线观看高清免费视频 | 99久久久亚洲精品蜜臀av| 宅男免费午夜| 麻豆av在线久日| 色尼玛亚洲综合影院| 日韩欧美一区二区三区在线观看| 久久久久国内视频| 91大片在线观看| 女同久久另类99精品国产91| 亚洲av成人av| 国产色视频综合| 午夜久久久久精精品| 国语自产精品视频在线第100页| 美女免费视频网站| 亚洲精品中文字幕一二三四区| 桃色一区二区三区在线观看| 黑人操中国人逼视频| 好男人电影高清在线观看| 久久久国产精品麻豆| 99香蕉大伊视频| 欧美一级a爱片免费观看看 | 久久香蕉精品热| 欧美午夜高清在线| 香蕉国产在线看| 国产在线精品亚洲第一网站| 天天躁夜夜躁狠狠躁躁| 极品教师在线免费播放| 欧美另类亚洲清纯唯美| 女同久久另类99精品国产91| 久久久水蜜桃国产精品网| 欧美成狂野欧美在线观看| 黄色a级毛片大全视频| 亚洲电影在线观看av| 国产成年人精品一区二区| 大香蕉久久成人网| 亚洲欧美日韩高清在线视频| 12—13女人毛片做爰片一| 国产精品免费一区二区三区在线| 一级作爱视频免费观看| 亚洲中文日韩欧美视频| 日本欧美视频一区| a在线观看视频网站| 国产aⅴ精品一区二区三区波| 国产亚洲精品久久久久久毛片| 亚洲激情在线av| 久久性视频一级片| 国产精品野战在线观看| 日本免费一区二区三区高清不卡 | 国产精品一区二区精品视频观看| 亚洲熟妇中文字幕五十中出| 无人区码免费观看不卡| 757午夜福利合集在线观看| 亚洲精品一区av在线观看| 久久精品人人爽人人爽视色| 狂野欧美激情性xxxx| 岛国在线观看网站| 亚洲男人天堂网一区| 99在线视频只有这里精品首页| 制服人妻中文乱码| 999久久久精品免费观看国产| www.999成人在线观看| 久久午夜亚洲精品久久| 午夜日韩欧美国产| 久久久国产欧美日韩av| 夜夜爽天天搞| 国产成人av激情在线播放| 少妇粗大呻吟视频| 国产一区二区激情短视频| 免费在线观看黄色视频的| 日日摸夜夜添夜夜添小说| 国产高清视频在线播放一区| 日韩欧美免费精品| 欧美中文综合在线视频| 亚洲三区欧美一区| 亚洲成av片中文字幕在线观看| 国产熟女午夜一区二区三区| 国产熟女xx| 国产又爽黄色视频| 欧美成人一区二区免费高清观看 | 国产欧美日韩精品亚洲av| 午夜精品久久久久久毛片777| 麻豆一二三区av精品| 日韩欧美国产一区二区入口| 母亲3免费完整高清在线观看| 高清在线国产一区| 久久欧美精品欧美久久欧美| 99国产精品一区二区蜜桃av| 欧美久久黑人一区二区| 精品电影一区二区在线| 满18在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 身体一侧抽搐| 免费看美女性在线毛片视频| 亚洲天堂国产精品一区在线| 自线自在国产av| 一夜夜www| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩黄片免| or卡值多少钱| 高清在线国产一区| 国产亚洲精品第一综合不卡| 桃红色精品国产亚洲av| 一本大道久久a久久精品| 精品国产美女av久久久久小说| 国产熟女午夜一区二区三区| 午夜激情av网站| 岛国视频午夜一区免费看| 99精品久久久久人妻精品| 精品卡一卡二卡四卡免费| 女人高潮潮喷娇喘18禁视频| 亚洲五月天丁香| 91大片在线观看| 悠悠久久av| 黑丝袜美女国产一区| 精品日产1卡2卡| 国产精品亚洲一级av第二区| 亚洲无线在线观看| 天堂影院成人在线观看| 日韩精品青青久久久久久| 欧美在线一区亚洲| 国产激情欧美一区二区| 他把我摸到了高潮在线观看| 国产亚洲欧美98| 少妇 在线观看| 女生性感内裤真人,穿戴方法视频| 精品久久蜜臀av无| 一进一出抽搐gif免费好疼| 久久人妻熟女aⅴ| 午夜久久久久精精品| 丝袜在线中文字幕| 啪啪无遮挡十八禁网站| 一区二区三区国产精品乱码| 国产免费av片在线观看野外av| 给我免费播放毛片高清在线观看| 亚洲国产精品合色在线| 亚洲av日韩精品久久久久久密| 午夜福利成人在线免费观看| 嫩草影视91久久| 母亲3免费完整高清在线观看| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 日本vs欧美在线观看视频| 九色亚洲精品在线播放| 精品国产乱码久久久久久男人| 久久国产精品人妻蜜桃| 国产精品久久久av美女十八| www.www免费av| 国产成人精品久久二区二区免费| 欧美av亚洲av综合av国产av| 欧美人与性动交α欧美精品济南到| 色婷婷久久久亚洲欧美| 99热只有精品国产| 精品午夜福利视频在线观看一区| 男人舔女人的私密视频| av超薄肉色丝袜交足视频| 精品久久久久久久久久免费视频| 一级a爱片免费观看的视频| 桃红色精品国产亚洲av| 美女国产高潮福利片在线看| or卡值多少钱| 日日夜夜操网爽| avwww免费| 一进一出抽搐动态| 精品少妇一区二区三区视频日本电影| 成在线人永久免费视频| 亚洲男人的天堂狠狠| 日本免费一区二区三区高清不卡 | 12—13女人毛片做爰片一| 大型av网站在线播放| 国产精品电影一区二区三区| 国产91精品成人一区二区三区| 在线av久久热| 老汉色∧v一级毛片| 国产成人精品久久二区二区91| 看免费av毛片| 午夜免费观看网址| 亚洲av电影不卡..在线观看| 女人高潮潮喷娇喘18禁视频| 制服诱惑二区| 女同久久另类99精品国产91| av网站免费在线观看视频| 9191精品国产免费久久| 黄色视频不卡| 久久久久精品国产欧美久久久| www.自偷自拍.com| 精品欧美国产一区二区三| 在线十欧美十亚洲十日本专区| 老汉色av国产亚洲站长工具| 中文字幕最新亚洲高清| 亚洲国产精品sss在线观看| 男人操女人黄网站| 激情在线观看视频在线高清| 免费看美女性在线毛片视频| 国产亚洲精品综合一区在线观看 | 日本 av在线| 久久久精品国产亚洲av高清涩受| 不卡一级毛片| 日韩大码丰满熟妇| 亚洲欧美精品综合一区二区三区| 亚洲一码二码三码区别大吗| 午夜日韩欧美国产| 女性生殖器流出的白浆| 看免费av毛片| www.999成人在线观看| 日韩欧美在线二视频| 日本五十路高清| 如日韩欧美国产精品一区二区三区| 国语自产精品视频在线第100页| 精品久久久久久久毛片微露脸| 久久人妻av系列| 美女免费视频网站| 成人国语在线视频| 国产av在哪里看| 男人操女人黄网站| 国产一级毛片七仙女欲春2 | 久久热在线av| 色综合亚洲欧美另类图片| www国产在线视频色| 伊人久久大香线蕉亚洲五| 嫁个100分男人电影在线观看| 又黄又爽又免费观看的视频| 搡老妇女老女人老熟妇| 香蕉久久夜色| av有码第一页| 免费一级毛片在线播放高清视频 | 啦啦啦观看免费观看视频高清 | 日日摸夜夜添夜夜添小说| 国产精品,欧美在线| 欧美一级毛片孕妇| 此物有八面人人有两片| 日韩有码中文字幕| 亚洲成人免费电影在线观看| 一级片免费观看大全| 亚洲专区字幕在线| 视频区欧美日本亚洲| 精品久久久久久久久久免费视频| 男女做爰动态图高潮gif福利片 | 欧美av亚洲av综合av国产av| 女人被狂操c到高潮| 日韩精品青青久久久久久| 欧美色视频一区免费| 一边摸一边做爽爽视频免费| 日韩大码丰满熟妇| 午夜两性在线视频| 禁无遮挡网站| 亚洲精品在线美女| 99riav亚洲国产免费| 色老头精品视频在线观看| 黑丝袜美女国产一区| 黄色毛片三级朝国网站| 久久久水蜜桃国产精品网| 久久人人爽av亚洲精品天堂| 在线观看免费视频网站a站| 在线观看免费日韩欧美大片| 国产区一区二久久| 国产97色在线日韩免费| 在线国产一区二区在线| 99久久久亚洲精品蜜臀av| 亚洲无线在线观看| 欧美丝袜亚洲另类 | 黄片播放在线免费| 国产精品久久久久久亚洲av鲁大| 午夜亚洲福利在线播放| 久久午夜综合久久蜜桃| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲精品第一综合不卡| 电影成人av| 国产成人免费无遮挡视频| 美女免费视频网站| 国产色视频综合| 精品国产国语对白av| 精品国产美女av久久久久小说| 久久草成人影院| 亚洲专区中文字幕在线| 亚洲狠狠婷婷综合久久图片| 一级毛片女人18水好多| 99在线人妻在线中文字幕| 夜夜躁狠狠躁天天躁| 午夜福利,免费看| 国产人伦9x9x在线观看| 在线永久观看黄色视频| 国产片内射在线| 久久精品aⅴ一区二区三区四区| 国产亚洲欧美精品永久| av在线播放免费不卡| 亚洲成国产人片在线观看| 国产野战对白在线观看| 国产区一区二久久| 亚洲第一av免费看| www.精华液| 久久亚洲精品不卡| 正在播放国产对白刺激| 人人妻人人爽人人添夜夜欢视频| 黄片播放在线免费| 黄色视频不卡| 亚洲av成人一区二区三| 日本免费a在线| 中文字幕精品免费在线观看视频| 一夜夜www| 亚洲va日本ⅴa欧美va伊人久久| 国产伦人伦偷精品视频| 手机成人av网站| av网站免费在线观看视频| 老熟妇乱子伦视频在线观看| 1024视频免费在线观看| 亚洲精品在线美女| 操出白浆在线播放| 亚洲精品一区av在线观看| 亚洲成a人片在线一区二区| av天堂在线播放| aaaaa片日本免费| av中文乱码字幕在线| 日本vs欧美在线观看视频| 亚洲精品一区av在线观看| 亚洲免费av在线视频| 午夜影院日韩av| 性欧美人与动物交配| 国产在线观看jvid| 国产精品 欧美亚洲| 免费av毛片视频| 亚洲精品国产色婷婷电影| 久久天堂一区二区三区四区| 国产精品一区二区免费欧美| 国产一区二区在线av高清观看| 中文字幕av电影在线播放| 久久久国产精品麻豆| 9191精品国产免费久久| 日本a在线网址| 亚洲一区中文字幕在线| 18美女黄网站色大片免费观看| av免费在线观看网站| 18美女黄网站色大片免费观看| 在线观看免费视频日本深夜| 国产色视频综合| 99久久久亚洲精品蜜臀av| 国产精品美女特级片免费视频播放器 | 电影成人av| 久热爱精品视频在线9| 高清毛片免费观看视频网站| 丝袜在线中文字幕| 很黄的视频免费| av中文乱码字幕在线| 久久精品国产99精品国产亚洲性色 | 欧美激情高清一区二区三区| 久久久久国产精品人妻aⅴ院| 99国产精品一区二区蜜桃av| 中文亚洲av片在线观看爽| 如日韩欧美国产精品一区二区三区| 少妇粗大呻吟视频| 操美女的视频在线观看| 在线观看一区二区三区| 两个人看的免费小视频| 国产一区二区在线av高清观看| 色播在线永久视频| 满18在线观看网站| 亚洲三区欧美一区| 国产麻豆成人av免费视频| 亚洲精品一区av在线观看| 欧美一级毛片孕妇| 国产一卡二卡三卡精品| 日日摸夜夜添夜夜添小说| 国产av精品麻豆| 一本大道久久a久久精品| cao死你这个sao货| 精品国产超薄肉色丝袜足j| 国产熟女午夜一区二区三区| 国产三级在线视频| 校园春色视频在线观看| 深夜精品福利| 男男h啪啪无遮挡| 国产精品av久久久久免费| 成人国产综合亚洲| 给我免费播放毛片高清在线观看| 神马国产精品三级电影在线观看 | 国产欧美日韩综合在线一区二区| 久久亚洲精品不卡| 91麻豆av在线| 欧美老熟妇乱子伦牲交| 我的亚洲天堂| 亚洲人成77777在线视频| 亚洲专区字幕在线| 中国美女看黄片| 麻豆久久精品国产亚洲av| 精品无人区乱码1区二区| 亚洲欧美日韩无卡精品| 国产高清videossex| 男女床上黄色一级片免费看| 国产亚洲精品av在线| 男女下面插进去视频免费观看| 亚洲人成网站在线播放欧美日韩| 午夜精品久久久久久毛片777| 午夜免费成人在线视频| 法律面前人人平等表现在哪些方面| 1024香蕉在线观看| 成人精品一区二区免费| 精品免费久久久久久久清纯| 欧美精品亚洲一区二区| 久9热在线精品视频| 人人妻人人澡欧美一区二区 | 亚洲久久久国产精品| 国产成人精品久久二区二区免费| 99国产精品一区二区蜜桃av| 身体一侧抽搐| 女人被狂操c到高潮| av视频免费观看在线观看| 91老司机精品| 精品久久久久久久毛片微露脸| 欧美激情久久久久久爽电影 | 亚洲少妇的诱惑av| 午夜福利成人在线免费观看| 国产精品一区二区三区四区久久 | 久久国产精品影院| 亚洲一卡2卡3卡4卡5卡精品中文| 成人三级黄色视频| 久久热在线av| 国产精品av久久久久免费| 欧美精品亚洲一区二区| 在线观看日韩欧美| 丁香欧美五月| 91老司机精品| 制服丝袜大香蕉在线| 女同久久另类99精品国产91| 精品国产一区二区久久| 久久狼人影院| 日韩大尺度精品在线看网址 | 亚洲熟女毛片儿| 在线观看免费日韩欧美大片| 亚洲三区欧美一区| 久久久久九九精品影院| 亚洲国产精品合色在线| 最近最新中文字幕大全免费视频| 高潮久久久久久久久久久不卡| 精品国产美女av久久久久小说| 校园春色视频在线观看| 一级a爱片免费观看的视频| 少妇的丰满在线观看| 亚洲自拍偷在线| 性少妇av在线| 精品人妻1区二区| 日本撒尿小便嘘嘘汇集6| 婷婷精品国产亚洲av在线| 无人区码免费观看不卡| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩一级在线毛片| 久久久久久人人人人人| 99在线视频只有这里精品首页| 亚洲精品中文字幕在线视频| 精品欧美国产一区二区三| 每晚都被弄得嗷嗷叫到高潮| 亚洲色图 男人天堂 中文字幕| 在线观看www视频免费| 久久九九热精品免费| 不卡一级毛片| 久99久视频精品免费| 天天躁夜夜躁狠狠躁躁| 精品一区二区三区四区五区乱码| 此物有八面人人有两片| 欧美丝袜亚洲另类 | 一区二区三区国产精品乱码| 最近最新免费中文字幕在线| 此物有八面人人有两片| 国产高清videossex| 精品久久久久久久毛片微露脸| av免费在线观看网站| 无人区码免费观看不卡| 丰满的人妻完整版| 国产亚洲精品av在线| 又大又爽又粗| av中文乱码字幕在线| 制服人妻中文乱码| 91av网站免费观看| 亚洲国产高清在线一区二区三 | 久久影院123| 精品久久久久久久毛片微露脸| 少妇 在线观看| 日本撒尿小便嘘嘘汇集6| 91麻豆精品激情在线观看国产| 韩国av一区二区三区四区| 日韩高清综合在线| 女人高潮潮喷娇喘18禁视频| 欧美久久黑人一区二区| 搡老妇女老女人老熟妇| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久成人av| 亚洲中文av在线| 精品福利观看| 日韩一卡2卡3卡4卡2021年| 亚洲成人国产一区在线观看| 中文字幕精品免费在线观看视频| 国产av又大| av在线播放免费不卡| 国产1区2区3区精品| 久久亚洲真实| 亚洲伊人色综图| 久久影院123| 免费女性裸体啪啪无遮挡网站| 麻豆久久精品国产亚洲av| 欧美另类亚洲清纯唯美| 淫秽高清视频在线观看| 夜夜看夜夜爽夜夜摸| 久久精品国产99精品国产亚洲性色 | 久久性视频一级片| 日本欧美视频一区| 黄色毛片三级朝国网站| 这个男人来自地球电影免费观看| 亚洲成人精品中文字幕电影| 久久青草综合色| 亚洲国产精品成人综合色| 18禁观看日本| 亚洲第一电影网av| 香蕉丝袜av| 少妇粗大呻吟视频| 国产精品永久免费网站| 免费看美女性在线毛片视频| 免费在线观看影片大全网站| 欧洲精品卡2卡3卡4卡5卡区| 天天添夜夜摸| 亚洲国产高清在线一区二区三 | 我的亚洲天堂| 免费搜索国产男女视频| 女人高潮潮喷娇喘18禁视频| 午夜免费鲁丝| av超薄肉色丝袜交足视频| 国产男靠女视频免费网站| 9色porny在线观看| 欧美在线一区亚洲| 精品人妻1区二区| 精品国产乱子伦一区二区三区| 夜夜爽天天搞| 国产三级在线视频| 久久青草综合色| 欧美成人一区二区免费高清观看 | 国产在线精品亚洲第一网站| 99国产精品一区二区蜜桃av| 亚洲伊人色综图| 午夜免费鲁丝| 一级,二级,三级黄色视频| 亚洲美女黄片视频| 日日爽夜夜爽网站| 一边摸一边抽搐一进一小说| 国产三级在线视频| 日日爽夜夜爽网站| 一本久久中文字幕| 男人舔女人下体高潮全视频| 亚洲午夜理论影院| 精品欧美国产一区二区三| 黄色女人牲交| 国产成+人综合+亚洲专区| 99香蕉大伊视频|