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

    肺部三維EIT模型構(gòu)建與圖像重建研究

    2017-12-28 12:56:54陳曉艷褚猛麗常曉敏章曉潔
    關(guān)鍵詞:電導(dǎo)率胸腔間距

    陳曉艷 褚猛麗 常曉敏 章曉潔

    (天津科技大學(xué)電子信息與自動化學(xué)院,天津 300222)

    肺部三維EIT模型構(gòu)建與圖像重建研究

    陳曉艷?褚猛麗 常曉敏 章曉潔

    (天津科技大學(xué)電子信息與自動化學(xué)院,天津 300222)

    肺部三維電阻抗圖像重建是電阻抗成像技術(shù)的重要應(yīng)用之一。利用3D光學(xué)云點(diǎn)技術(shù)對人體胸腔區(qū)域掃描,并融合X圖像提供的肺部結(jié)構(gòu)構(gòu)建肺部三維EIT仿真模型。根據(jù)肺部膨脹及收縮時的電導(dǎo)率先驗(yàn)知識,由COMSOL軟件求解獲得三維靈敏度矩陣,并在Matlab平臺下由共軛梯度迭代算法重建肺部EIT斷層圖像,進(jìn)行三維立體重構(gòu),獲得3D EIT圖像。為研究電流敏感場均勻性,設(shè)置4組不同電極層間距進(jìn)行比較仿真實(shí)驗(yàn),結(jié)果表明,對于33 cm胸腔區(qū)域,電極層間距為8 cm時,成像效果最佳,其最大相關(guān)系數(shù)為0.810 3,敏感場均勻性為1.869 6×103,結(jié)構(gòu)相似度為0.482 5,靈敏場的均勻性明顯改善,有利于圖像重建質(zhì)量的提高。

    電阻抗成像;肺部三維;圖像重建;共軛梯度迭代;先驗(yàn)信息

    引言

    EIT是一種新穎的無損成像技術(shù),它根據(jù)物體內(nèi)部組織電特性,對其表面施加安全電壓/電流激勵,通過測量物體邊界電信號,重建物體內(nèi)部電特性參數(shù)的分布及其變化情況[1-2]。電阻抗成像技術(shù)在醫(yī)學(xué)領(lǐng)域、工業(yè)領(lǐng)域的應(yīng)用日趨廣泛,尤其在醫(yī)學(xué)領(lǐng)域[3]更為社會各界所關(guān)注,如在肺通氣監(jiān)測[4-5]、腦組織變化跟蹤[6]、胃排空與收縮狀態(tài)檢測[7]等方面。目前在國內(nèi)肺部三維EIT成像研究中,天津大學(xué)王化祥小組采用人體胸腔輪廓,但僅局限于二維圖像重建;河北工業(yè)大學(xué)徐桂芝小組利用4層電極在圓柱體邊界模型下進(jìn)行肺部成像,實(shí)現(xiàn)了動態(tài)肺呼吸過程圖像重建[8];國外的曼徹斯特大學(xué)EIT課題小組利用MR圖像構(gòu)建胸腔輪廓模型,并結(jié)合分割方法與信息相互計(jì)算進(jìn)行肺部三維圖像重建[9]。上述圖像重建模型各有優(yōu)點(diǎn),但也各自存在缺陷,不斷探究更優(yōu)的圖像重建模型是三維EIT圖像重建的重要研究方向。

    EIT 電極分布[10-11]影響三維圖像重建質(zhì)量[12]。電極層間距離越大,敏感場均勻性越差;但是層間距離越小,測量的范圍就受限制,不能真實(shí)地反映胸腔內(nèi)肺部的空間三維分布情況。本研究對不同電極層間距離進(jìn)行對比、仿真成像試驗(yàn),采用相關(guān)系數(shù)、均勻場靈敏度性、結(jié)構(gòu)相似度3種評價指標(biāo)進(jìn)行電極層間距離的優(yōu)化。

    1 方法

    采用光學(xué)云點(diǎn)掃描構(gòu)建人體胸腔三維輪廓,并在胸腔模型中按真實(shí)人體肺臟構(gòu)建模型,以融合邊界結(jié)構(gòu)與肺臟部分布信息的方法提高三維EIT模型的準(zhǔn)確性。采用共軛梯度算法進(jìn)行三維圖像重建,在不同電極層間距條件下對圖像重建效果進(jìn)行評價判定。

    1.1 模型構(gòu)建

    1.1.1 胸腔輪廓提取

    為構(gòu)建準(zhǔn)確的人體胸廓模型,采用3D CaMega光學(xué)云點(diǎn)數(shù)字化掃描系統(tǒng),通過多方位、多角度拍攝掃描獲取三維云點(diǎn)數(shù)據(jù),由反求軟件與三維云點(diǎn)處理軟件自動生成高精度的人體胸腔輪廓。

    圖1為本研究中健康男性志愿者胸腔輪廓掃描實(shí)景,右側(cè)為掃描實(shí)物,左側(cè)為上位機(jī)顯示界面。

    圖1 3D光學(xué)掃描人體胸部實(shí)景Fig.1 Optical scanning human chest system

    圖2為上述系統(tǒng)掃描所得的初步模型。然后,利用三維軟件SolidWorks對其進(jìn)行適當(dāng)修正,以便降低后期圖像重建時逆問題的求解難度。

    圖2 3D光學(xué)掃描系統(tǒng)所得胸部初始模型Fig.2 The initial model by Optical scanning system

    1.1.2 肺臟模型

    由圖3(a)可見,人體兩片肺葉自上而下呈擴(kuò)張趨勢,心臟的存在導(dǎo)致左側(cè)肺葉面積要小于右側(cè)肺葉面積,且左右肺間距自上而下逐漸變大。

    通過X片可以獲取直觀有效的肺部生理結(jié)構(gòu)信息,并參考相關(guān)文獻(xiàn)、3D肺部圖片等信息,按照胸腔的尺寸,以一定的比例,利用3D軟件SolidWorks,構(gòu)建出橢圓形肺臟模型,分別為正視圖、側(cè)視圖、俯視圖,如圖3(b)所示。

    在呼吸過程中,肺組織輪廓存在較大差異。在吸氣末時,肺器官組織處于擴(kuò)張充盈狀態(tài),更便于觀察其內(nèi)部生理結(jié)構(gòu)信息。本研究選取吸氣末狀態(tài),作為研究肺部三維EIT圖像重建的胸腔模型。

    圖3 肺X片與仿真模型三視圖。(a)肺X圖片;(b)肺模型三視圖Fig.3 Lung X-ray images and simulation model.(a)Lung X-ray images;(b)Three view of lung simulation model

    1.1.3 胸肺組合模型

    將兩片肺葉模型按實(shí)際人體肺器臟狀態(tài)填充至胸腔模型內(nèi)部,形成胸肺組合模型。然后,將此胸肺組合模型導(dǎo)入有限元仿真軟件COMSOL 3.5a,在模型輪廓外邊界,均勻設(shè)置上下兩層各16個半徑為2.5 cm的圓形電極,如圖4所示。模型尺寸比例為掃描系統(tǒng)自適應(yīng)設(shè)置的最佳模式,模型縱向高33 cm,上部左右最寬為38 cm,下部左右最寬為22 cm。

    圖4 胸腔組合仿真模型Fig.4 Chest combination simulation model

    1.2 重建算法

    在邊界不規(guī)則、內(nèi)部結(jié)構(gòu)復(fù)雜的情況下,三維逆問題求解時系統(tǒng)矩陣維數(shù)增大,求解難度顯著增加。由于共軛梯度迭代(conjugate gradient,CG)算法[13]是通過每次迭代產(chǎn)生一個新的解,趨向逼近于真實(shí)解,無需進(jìn)行系統(tǒng)矩陣的分解,因此迭代時間較短,所以本研究采用CG算法進(jìn)行圖像重建。

    1.3 評價參數(shù)

    1)相關(guān)系數(shù)R[14]。該值靈敏地反映重建電導(dǎo)率和真實(shí)電導(dǎo)率之間相關(guān)關(guān)系的密切程度,它的取值范圍為[0,1],其值越大,逆問題求解質(zhì)量越好。

    2)敏感場靈敏度均勻性參數(shù)P[15]。在電學(xué)成像系統(tǒng)中,敏感場分布通常是不均勻的。P為評價敏感場分布參數(shù),其值越小,敏感場的分布越理想。

    3)結(jié)構(gòu)相似度SSIM[16]。該值表示重建圖像的非線性改變程度,主要考慮了圖像的亮度和對比度。SSIM的取值范圍為,值越大,圖像重建質(zhì)量越好。

    1.4 仿真實(shí)驗(yàn)

    仿真平臺參數(shù)如下:COMSOL3.5a,Matlab R2010b,Intel(R)Core i7 3.50GHz,8.00GB RAM。

    理論上的分析認(rèn)為,電極層數(shù)多,能更精準(zhǔn)反映三維的分布情況,圖像重建質(zhì)量會更好。但是,增加電極層數(shù)意味著圖像重建時間加長,計(jì)算的成本會呈指數(shù)型增長,對于仿真平臺的要求也會提高,在實(shí)現(xiàn)上存在一定的局限性。

    在本研究中,取電極層距離分別為6、7、8、9 cm,進(jìn)行4組仿真試驗(yàn),采用評價參數(shù)R、P和SSIM對成像結(jié)果對比分析。在有限元軟件COMSOL中,對構(gòu)建的空場模型(肺呼氣末狀態(tài))與滿場模型(肺吸氣末狀態(tài))分別施加大小為5 mA的直流電流激勵,考慮肺部電導(dǎo)率先驗(yàn)信息[17]的情況,空場下胸腔體與肺器官電導(dǎo)率均為0.037 s/m,滿場下胸腔體電導(dǎo)率與空場一致,肺部電導(dǎo)率設(shè)為0.138 s/m,采用對角激勵-相鄰測量模式采集成像數(shù)據(jù)[18-19]。

    2 結(jié)果

    2.1 流線圖對比

    在EIT測量中,電場的分布狀態(tài)對場域內(nèi)電流密度分布、場強(qiáng)大小及均勻性有著直接的影響。電流流線圖可以反映電流在場域內(nèi)方向、電流密度以及分布狀態(tài)。在本研究中,以不同電極層距離的(a)~(d)等4組試驗(yàn)的電流流線分布情況,探究其對成像效果的影響。為便于觀察,取每組仿真的一對激勵所得電流流線分布作為研究對象,4組流線分布狀態(tài)的仿真結(jié)果如圖5所示。由圖可見,每對激勵的電流流線均由電流輸入端到達(dá)電流的流出端。但是,在(a)、(b)和(d)中,有一部分電流線并不能到達(dá)流出端,表明電場分布不均勻;(c)的電流流線在整個場域內(nèi)各個部位的分布相對均勻,提示其電場分布較其他組均勻,與其相應(yīng)的靈敏度矩陣準(zhǔn)確度較高。

    2.2 圖像重建

    為了充分體現(xiàn)三維斷層成像效果,更好展示三維胸腔內(nèi)部不同斷層的生理結(jié)構(gòu)信息差異,模型自上而下取距最底部為25、22、17、15、11、8 cm 的 6 處位置斷層成像,對應(yīng)圖5中從(a)到(d)的順序,分別進(jìn)行圖像重建,成像結(jié)果如圖6(b)~(e)所示,(a)為仿真模型中6處位置的斷層圖像結(jié)果。由于仿真結(jié)果在迭代次數(shù)達(dá)150左右時均取得較好成像質(zhì)量,再增加迭代次數(shù)差別很小,因此在圖6中圖像統(tǒng)一取迭代次數(shù)為150。

    圖5 不同電極層間距時電場分布狀態(tài)仿真。(a)6 cm;(b)7 cm;(c)8 cm;(d)9 cmFig.5 The distribution of electric field under different electrode layer intervals.(a)6 cm;(b)7 cm;(c)8 cm;(d)9 cm

    圖6中每組成像結(jié)果均可體現(xiàn)胸腔內(nèi)肺葉自上而下的形狀變化、大小差異、位置距離改變等信息,說明成像結(jié)果在一定程度上可以反映實(shí)際人體模型肺部構(gòu)造。比較不同電極間隔距離的4組仿真成像效果,可見(d)組成像偽影少,目標(biāo)物位置準(zhǔn)確,邊界較清晰,成像質(zhì)量相對更好。

    表1為(b)~(e)4組仿真圖像的R、P和 SSIM參數(shù)評價結(jié)果,可見結(jié)果相差很小。但相對而言,(d)組的敏感場均勻性及結(jié)構(gòu)相似度較其他組更優(yōu),這定量驗(yàn)證了對圖6結(jié)果的觀察結(jié)論。

    表1 不同電極層間距的參數(shù)評價結(jié)果Tab.1 The evaluation indices of different electrode layer interval

    2.3 三維立體重建

    將上述6幅EIT斷層EIT圖像在Matlab軟件中進(jìn)行三維圖像構(gòu)建。首先,對EIT圖像進(jìn)行灰度化處理,提取邊界信息及感興趣區(qū)的輪廓曲線;其次,將預(yù)處理好的6幅EIT圖像導(dǎo)入Matlab,使用cat函數(shù)構(gòu)建三維矩陣D,經(jīng)smooth等函數(shù)對圖像進(jìn)行平滑、邊緣銳化及偽彩色增強(qiáng)等一系列可視化處理,最后由isosurface、isonormals、patch等函數(shù)實(shí)現(xiàn)三維立體重建,結(jié)果如圖7所示。

    圖6 實(shí)際模型與4種不同電極層間距下成像(每行從左到右分別為模型自上而下6處位置的斷層成像)。(a)實(shí)際模型;(b)~(e)電極層間距6~9 cmFig.6 The real model and imaging results with four different electrode layers spacings(From left to right in each row,tomography results of top-down with 6 different positions).(a)The real model;(b)~(e)Electrode layers spacing from 6~9 cm

    通過三維立體圖像可看出,由EIT二維斷層圖像所構(gòu)建的胸腔三維立體圖像,增強(qiáng)了邊界輪廓,豐富了結(jié)構(gòu)信息,減少了噪聲干擾,再現(xiàn)了圖3模型的三維可視化效果。在圖7中,邊界上一層層的輪廓相對還是比較粗糙的,如果導(dǎo)入的EIT圖像數(shù)量增加,將會得到邊界更加光滑的三維圖像。

    圖7 三維重建立體圖像Fig.7 The 3D reconstruction stereogram

    3 討論

    本研究結(jié)果表明,三維EIT電極分布并不是電極層間隔越大或越小成像效果就越好,這主要是由于電阻抗成像中注入激勵電流是一個三維電流場,場強(qiáng)分布不會限于同一個平面,而會向?qū)щ娪虻母鱾€方向分布。在三維肺部圖像重建過程中,電流場效應(yīng)不可忽視,電極層間距改變使得激勵電流形成的三維電流場隨之變化:如間距變大,場域內(nèi)的電流密度變小,場域分布均勻性越差,從而影響成像效果;如間距變小,兩電極層的電流場效應(yīng)就會相互產(chǎn)生影響,使得敏感場的非均勻性加重,最終成像效果也會受到影響。

    三維電流場效應(yīng)對被測場域的大小、形狀都存在著較強(qiáng)的依賴關(guān)系。本研究在建立的模型大小與形狀基礎(chǔ)上,通過仿真試驗(yàn)及參數(shù)評價方法,找出最佳的電流場域分布狀態(tài),形成相對較優(yōu)的電極層分布方案,為后續(xù)的三維胸腔電阻抗成像應(yīng)用研究奠定了理論基礎(chǔ)。本研究的不足之處在于:首先,只采用了4種不同電極層間距進(jìn)行研究,若增加層間距對比數(shù)量,細(xì)化模型求解域,可取得更為理想的成像效果;其次,電阻抗技術(shù)的廣泛應(yīng)用并不局限于肺部的成像檢測,在腦部、乳腺、腸胃以及其他應(yīng)用中也較為普遍。因而,對于不同生理結(jié)構(gòu)和生理狀態(tài)的組織器官,需要對其電極分布采取相應(yīng)合理有效的電極布置方案,這也是EIT技術(shù)要繼續(xù)研究的一個重要內(nèi)容。

    4 結(jié)論

    在本研究中,利用3D光學(xué)云點(diǎn)技術(shù)對自愿者的胸腔區(qū)域進(jìn)行掃描,并融合X圖片提供的肺部結(jié)構(gòu)先驗(yàn)信息和肺部電導(dǎo)率先驗(yàn)信息,構(gòu)建三維肺部EIT仿真模型,為三維胸腔模型構(gòu)建提供了又一新的重構(gòu)方式。通過圖像重建算法,以電極層間距分布為討論研究點(diǎn),重構(gòu)三維EIT的斷層圖像,經(jīng)評價參數(shù)的定量驗(yàn)證,重構(gòu)質(zhì)量較好,得出最佳電極層間距分布模式,且進(jìn)一步驗(yàn)證該模型的創(chuàng)新合理性、有效性。最后將胸腔區(qū)域的6層斷層圖像重建出三維立體圖像,提供了更為豐富的空間電導(dǎo)率分布信息。

    人體胸腔具有非常復(fù)雜的生理結(jié)構(gòu),其電導(dǎo)率分布情況也因個體差異而有所不同。因此,在下一步的研究中,將胸腔內(nèi)部的器官組織(如心臟、脊柱、皮下組織等)進(jìn)一步考慮到三維建模中,會獲取更為豐富的生物組織電導(dǎo)率分布狀態(tài);此外,增加模型樣本容量,提取更多的模型輪廓,提高模型擬合度,進(jìn)而獲取更為精準(zhǔn)的胸腔模型。

    [1] Brown BH.Electrical impedance tomography(EIT):A review[J].Med Eng Technol,2003,27(3):97-108.

    [2] Bayford RH,Bioimpedance tomography(electrical impedance tomography)[J].Annu Rev Biomed Eng,2006,8:63-91.

    [3] 董秀珍.生物電阻抗成像研究的現(xiàn)狀與挑戰(zhàn)[J].中國生物醫(yī)學(xué)工程學(xué)報,2008,27(5):641-649.

    [4] Chen Xiaoyan,Wang Huaxiang.Lung ventilation functional monitoring based on electrical impedance tomography [J].Trans Tianjin Univ,2009,15:007-012.

    [5] Heizmann S,Baumg?rtner M,Krüger-Ziolek S,et al.3-D lung visualization using electrical impedance tomography combined with body plethysmography [J].IFMBE Proc,2014,43(1):172-175.

    [6] Gilad O,Holder DS.Impedance changes recorded with scalp electrodes duringvisualevoked responses:implications for electrical impedance tomography of fast neural activity [J].Neuro Image,2009,47(2):514-522.

    [7] Gaggero PO,Adler A,Brunner J,et al.Electrical impedance tomography system based on active electrodes [J].Physiological Measurement,2012,33(5):831-847.

    [8] Zhang Shuai,Xu Guizhi,Zhang Xueying,et al.Computation of a 3-D model for lung imaging with electrical impedance tomography[C].IEEE Trans Magn,2012,48(2):651-654.

    [9] Crabb,MG,Davidson,JL,Little,R,et al.Mutual information as a measure of reconstruction quality in 3D dynamic lung EIT[C].J Phys Conf Ser,2013,434(1):012082

    [10] 王妍,任超世.3 D-EIT圖像重建的研究進(jìn)展 [J].國外醫(yī)學(xué)生物醫(yī)學(xué)分冊,2003,26(6):265-268.

    [11] 陳曉艷,李宏英,趙秋紅,等.3D-EIT系統(tǒng)電極優(yōu)化設(shè)計(jì)[J].中國生物醫(yī)學(xué)工程學(xué)報,2014,33(2):9-14.

    [12] Moreno MV,Chaset L,Bittner PA,et al.New reusable elastomer electrodes for assessing body composition [J].J Phys Conf Ser,2013,434(1):012026.

    [13] 侯雪,王超.基于共軛梯度算法EIT仿真[J].長春理工大學(xué)學(xué)報(自然科學(xué)版),2012,35(4):167-170.

    [14] 范文茹,郝魁紅,許燕斌.開放式電阻抗成像激勵策略研究[J]. 儀器儀表學(xué)報,2014,35(6):1269-1275.

    [15] 鄧娟,魏慶國,沙洪,等.三維電阻抗斷層成像仿真研究及圖像重建質(zhì)量評價[J].生物醫(yī)學(xué)工程學(xué)雜志,2013,30(6):1195-1199.

    [16] Xu Canhua,Dai Meng,You Fusheng,et al.An optimized strategy for real-time hemorrhage monitoring with electrical impedance tomography [J].Physiol Meas,2011,(32):585-598.

    [17] 王化祥,范文茹,馬雪翠.基于先驗(yàn)信息的肺部電阻抗成像算法[J]. 中國生物醫(yī)學(xué)工程學(xué)報,2009,28(5):680-685.

    [18] 陳曉艷,楊永政.一種肺呼吸電阻抗實(shí)時成像系統(tǒng)[J].電子技術(shù)應(yīng)用,2015,41(5):73-76.

    [19] 陳曉艷,高娜娜,黃華芳.數(shù)字式生物電阻抗與電阻抗斷層成像測量系統(tǒng)[J].生物醫(yī)學(xué)工程學(xué)雜志,2015,32(3):575-580.

    3D EIT Model Construction Based on Lung and Image Reconstruction Research

    Chen Xiaoyan Chu Mengli Chang Xiaomin Zhang Xiaojie
    (College of Information and Automation,Tianjin University of Science and Technology,Tianjin300222,China)

    electrical impedance imaging;three-dimensional lung;image reconstruction;CG algorithm;prior information

    R318.08 文獻(xiàn)標(biāo)志碼:D 文章編號:0258-8021(2017)05-0622-05

    10.3969 /j.issn.0258-8021.2017.05.016

    2016-05-24,錄用日期:2017-06-22

    國家自然科學(xué)基金(61301246);天津市應(yīng)用基礎(chǔ)及前沿技術(shù)研究計(jì)劃(12JCYBJC19300)

    ?通信作者(Corresponding author),E-mail:cxywxr@tust.edu.cn

    猜你喜歡
    電導(dǎo)率胸腔間距
    宋代君臣殿上間距考論
    胸腔巨大孤立性纖維瘤伴多發(fā)轉(zhuǎn)移1例
    胸外科術(shù)后胸腔引流管管理的研究進(jìn)展
    高速公路指揮中心小間距LED應(yīng)用探討
    基于比較測量法的冷卻循環(huán)水系統(tǒng)電導(dǎo)率檢測儀研究
    低溫脅迫葡萄新梢電導(dǎo)率和LT50值的研究
    算距離
    尿激酶聯(lián)合抗結(jié)核藥胸腔內(nèi)注入對結(jié)核性包裹性胸腔積液治療及其預(yù)后的影響
    高電導(dǎo)率改性聚苯胺的合成新工藝
    乳腺癌肺轉(zhuǎn)移胸腔積液的治療(附36例報告)
    一级a爱视频在线免费观看| 叶爱在线成人免费视频播放| 午夜影院日韩av| 熟妇人妻久久中文字幕3abv| 欧美一级毛片孕妇| 禁无遮挡网站| 国产亚洲精品久久久久久毛片| 女人被狂操c到高潮| 日本五十路高清| 99久久精品国产亚洲精品| 免费在线观看完整版高清| 两性午夜刺激爽爽歪歪视频在线观看 | 无人区码免费观看不卡| 欧美亚洲日本最大视频资源| 69av精品久久久久久| 国产乱人伦免费视频| 校园春色视频在线观看| 久久天堂一区二区三区四区| 韩国精品一区二区三区| 757午夜福利合集在线观看| 亚洲在线自拍视频| 极品人妻少妇av视频| 色精品久久人妻99蜜桃| 一夜夜www| 亚洲av第一区精品v没综合| 亚洲国产欧美一区二区综合| 亚洲av成人av| 婷婷精品国产亚洲av在线| 精品久久久久久久人妻蜜臀av | 欧洲精品卡2卡3卡4卡5卡区| 久久久久久久精品吃奶| 国产欧美日韩精品亚洲av| 日韩成人在线观看一区二区三区| 国产男靠女视频免费网站| 欧美日本亚洲视频在线播放| 精品一品国产午夜福利视频| 超碰成人久久| 亚洲在线自拍视频| www.精华液| 在线观看日韩欧美| 一个人观看的视频www高清免费观看 | av网站免费在线观看视频| 午夜a级毛片| 纯流量卡能插随身wifi吗| 亚洲人成网站在线播放欧美日韩| 中文字幕高清在线视频| 国产精品久久电影中文字幕| 夜夜夜夜夜久久久久| 精品欧美一区二区三区在线| 亚洲欧洲精品一区二区精品久久久| 99久久精品国产亚洲精品| 日韩有码中文字幕| 欧美成人一区二区免费高清观看 | 免费看美女性在线毛片视频| 午夜福利欧美成人| 亚洲精品国产区一区二| 欧美乱色亚洲激情| 看免费av毛片| 伦理电影免费视频| 精品午夜福利视频在线观看一区| 免费在线观看黄色视频的| 国产97色在线日韩免费| 精品午夜福利视频在线观看一区| 黄色a级毛片大全视频| 国产aⅴ精品一区二区三区波| 极品人妻少妇av视频| 国产精品永久免费网站| 中文字幕色久视频| 制服诱惑二区| 欧美在线一区亚洲| 国产成人精品久久二区二区免费| 亚洲一区二区三区不卡视频| 777久久人妻少妇嫩草av网站| 乱人伦中国视频| 精品国产超薄肉色丝袜足j| 一本综合久久免费| 精品国产乱子伦一区二区三区| 亚洲国产中文字幕在线视频| 久久天躁狠狠躁夜夜2o2o| 露出奶头的视频| 亚洲av美国av| 成人三级黄色视频| 在线观看午夜福利视频| 欧美激情 高清一区二区三区| 色在线成人网| 99久久久亚洲精品蜜臀av| 日本在线视频免费播放| 校园春色视频在线观看| 一进一出抽搐动态| 日韩大码丰满熟妇| 91精品三级在线观看| 午夜福利一区二区在线看| 亚洲国产精品sss在线观看| 亚洲全国av大片| 老汉色∧v一级毛片| 97人妻天天添夜夜摸| 亚洲一区中文字幕在线| 精品国内亚洲2022精品成人| 成人18禁高潮啪啪吃奶动态图| 自拍欧美九色日韩亚洲蝌蚪91| 操出白浆在线播放| 国产成人av教育| 久久久国产成人免费| 一级片免费观看大全| 成在线人永久免费视频| 99精品欧美一区二区三区四区| 久久亚洲精品不卡| 欧美日本亚洲视频在线播放| 免费少妇av软件| 人人妻人人澡欧美一区二区 | 精品免费久久久久久久清纯| 免费观看人在逋| 免费高清在线观看日韩| or卡值多少钱| 欧美激情高清一区二区三区| 免费高清在线观看日韩| 天天添夜夜摸| 露出奶头的视频| 日韩欧美免费精品| 制服丝袜大香蕉在线| 国产一区二区三区在线臀色熟女| 19禁男女啪啪无遮挡网站| 亚洲av电影在线进入| 日韩欧美三级三区| 老汉色av国产亚洲站长工具| 亚洲,欧美精品.| 91精品三级在线观看| 亚洲av美国av| АⅤ资源中文在线天堂| 岛国视频午夜一区免费看| 一区二区日韩欧美中文字幕| 日韩欧美一区视频在线观看| 精品人妻在线不人妻| 一级a爱片免费观看的视频| 亚洲av成人不卡在线观看播放网| 欧美+亚洲+日韩+国产| 午夜视频精品福利| 18禁观看日本| 天天躁夜夜躁狠狠躁躁| 乱人伦中国视频| 久久欧美精品欧美久久欧美| www.自偷自拍.com| 男人舔女人的私密视频| 伊人久久大香线蕉亚洲五| 人人妻,人人澡人人爽秒播| 人人妻,人人澡人人爽秒播| 国产av在哪里看| 91成人精品电影| 日韩精品青青久久久久久| 午夜视频精品福利| 女人高潮潮喷娇喘18禁视频| 伊人久久大香线蕉亚洲五| 久久欧美精品欧美久久欧美| 亚洲av成人av| 香蕉丝袜av| 99精品久久久久人妻精品| 亚洲精品久久国产高清桃花| 久久精品成人免费网站| 日韩大码丰满熟妇| 在线观看免费日韩欧美大片| 精品欧美国产一区二区三| 一级,二级,三级黄色视频| 欧美成人一区二区免费高清观看 | av在线天堂中文字幕| 国内毛片毛片毛片毛片毛片| а√天堂www在线а√下载| 亚洲五月天丁香| 亚洲色图 男人天堂 中文字幕| 人人妻人人澡人人看| 桃色一区二区三区在线观看| 久久伊人香网站| 欧美最黄视频在线播放免费| 99国产精品一区二区三区| 亚洲五月色婷婷综合| 日韩免费av在线播放| 男女午夜视频在线观看| 久久亚洲精品不卡| 久久青草综合色| 久久伊人香网站| 脱女人内裤的视频| 欧美一级毛片孕妇| 法律面前人人平等表现在哪些方面| 亚洲精品国产一区二区精华液| 久久这里只有精品19| 69精品国产乱码久久久| 国产精品九九99| 久久久精品欧美日韩精品| 日本免费a在线| 国产成人影院久久av| 精品熟女少妇八av免费久了| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲av一区麻豆| 国产av一区在线观看免费| 一级毛片女人18水好多| 老司机深夜福利视频在线观看| 成人国产综合亚洲| 国产精品免费视频内射| 午夜精品国产一区二区电影| 久久人人精品亚洲av| 日韩欧美一区视频在线观看| 99国产精品免费福利视频| 免费看美女性在线毛片视频| 日韩高清综合在线| 99香蕉大伊视频| 久久久久国产一级毛片高清牌| 日本三级黄在线观看| 欧美最黄视频在线播放免费| 女人爽到高潮嗷嗷叫在线视频| 欧美国产精品va在线观看不卡| 午夜激情av网站| 精品电影一区二区在线| 一级a爱片免费观看的视频| 美女高潮喷水抽搐中文字幕| 制服丝袜大香蕉在线| 桃色一区二区三区在线观看| 香蕉久久夜色| 亚洲第一av免费看| 国产精品二区激情视频| 香蕉丝袜av| 免费看十八禁软件| 免费在线观看亚洲国产| 婷婷丁香在线五月| 欧美久久黑人一区二区| 亚洲 欧美一区二区三区| 在线国产一区二区在线| 真人做人爱边吃奶动态| 午夜日韩欧美国产| 窝窝影院91人妻| 一区二区三区国产精品乱码| 真人做人爱边吃奶动态| 美国免费a级毛片| 777久久人妻少妇嫩草av网站| 999精品在线视频| 亚洲自拍偷在线| 日韩欧美一区视频在线观看| 男女床上黄色一级片免费看| 色综合站精品国产| 亚洲一码二码三码区别大吗| 青草久久国产| 非洲黑人性xxxx精品又粗又长| 久久久国产成人精品二区| 午夜福利成人在线免费观看| 日日摸夜夜添夜夜添小说| 一本大道久久a久久精品| 欧美精品亚洲一区二区| 久久久久久免费高清国产稀缺| 亚洲伊人色综图| 两个人免费观看高清视频| 亚洲av成人av| 国产精品av久久久久免费| 久久精品91无色码中文字幕| 人人妻人人澡人人看| 中出人妻视频一区二区| 国产日韩一区二区三区精品不卡| 国产成人精品久久二区二区免费| 激情在线观看视频在线高清| 黄色 视频免费看| 亚洲第一av免费看| 国产成年人精品一区二区| 美女高潮到喷水免费观看| 如日韩欧美国产精品一区二区三区| 桃色一区二区三区在线观看| 成人亚洲精品av一区二区| 丝袜美足系列| 欧美一区二区精品小视频在线| 国产91精品成人一区二区三区| 在线天堂中文资源库| 大型av网站在线播放| 这个男人来自地球电影免费观看| 国产极品粉嫩免费观看在线| 色综合欧美亚洲国产小说| 亚洲国产精品久久男人天堂| 又大又爽又粗| 欧美日本中文国产一区发布| 亚洲五月天丁香| 夜夜夜夜夜久久久久| a级毛片在线看网站| 69av精品久久久久久| 女警被强在线播放| 久久伊人香网站| 国产成+人综合+亚洲专区| 99久久久亚洲精品蜜臀av| 午夜福利一区二区在线看| 麻豆国产av国片精品| 女人被狂操c到高潮| 成人欧美大片| 亚洲成人免费电影在线观看| 久久久国产欧美日韩av| 国产主播在线观看一区二区| 大陆偷拍与自拍| 国产欧美日韩一区二区精品| 精品午夜福利视频在线观看一区| 波多野结衣巨乳人妻| 动漫黄色视频在线观看| 一a级毛片在线观看| 亚洲av日韩精品久久久久久密| 亚洲精华国产精华精| 亚洲欧美激情综合另类| 久久青草综合色| 激情在线观看视频在线高清| 久久久久精品国产欧美久久久| 国产激情欧美一区二区| 欧美精品亚洲一区二区| 亚洲成人精品中文字幕电影| 99国产精品免费福利视频| 亚洲av第一区精品v没综合| 亚洲色图 男人天堂 中文字幕| 久久天堂一区二区三区四区| 搡老妇女老女人老熟妇| 国产麻豆成人av免费视频| 国产aⅴ精品一区二区三区波| 国产免费男女视频| 国产精品久久久久久亚洲av鲁大| 亚洲精品美女久久久久99蜜臀| 久久天躁狠狠躁夜夜2o2o| 色av中文字幕| 一本大道久久a久久精品| 亚洲av成人不卡在线观看播放网| 亚洲男人天堂网一区| 日韩高清综合在线| 午夜久久久在线观看| 少妇熟女aⅴ在线视频| 亚洲精品在线观看二区| 亚洲精品国产精品久久久不卡| 午夜免费激情av| 美女扒开内裤让男人捅视频| 91在线观看av| 大型黄色视频在线免费观看| 18禁国产床啪视频网站| 电影成人av| 多毛熟女@视频| 激情在线观看视频在线高清| 天天躁狠狠躁夜夜躁狠狠躁| 99国产综合亚洲精品| 国内精品久久久久久久电影| 亚洲在线自拍视频| 中文亚洲av片在线观看爽| 精品人妻在线不人妻| 久久这里只有精品19| 日本撒尿小便嘘嘘汇集6| 一区二区三区国产精品乱码| 女人被躁到高潮嗷嗷叫费观| 亚洲精品美女久久久久99蜜臀| 久久性视频一级片| 成人三级做爰电影| 无限看片的www在线观看| 午夜福利18| 国产亚洲精品久久久久5区| 操美女的视频在线观看| 热re99久久国产66热| 老熟妇仑乱视频hdxx| 国产精品九九99| 国产欧美日韩综合在线一区二区| 亚洲 欧美一区二区三区| 一区二区三区高清视频在线| 欧美激情 高清一区二区三区| www.www免费av| 午夜久久久久精精品| 亚洲中文字幕日韩| 久久亚洲精品不卡| 亚洲五月天丁香| 在线播放国产精品三级| 在线观看免费视频日本深夜| 欧美激情 高清一区二区三区| 嫁个100分男人电影在线观看| svipshipincom国产片| 国产精品免费视频内射| 国产一卡二卡三卡精品| 一区二区三区精品91| 啦啦啦观看免费观看视频高清 | 老鸭窝网址在线观看| 成人三级做爰电影| 午夜激情av网站| 黄色女人牲交| 嫁个100分男人电影在线观看| 九色国产91popny在线| 成人av一区二区三区在线看| av在线天堂中文字幕| 亚洲成国产人片在线观看| 一边摸一边抽搐一进一出视频| 午夜a级毛片| 给我免费播放毛片高清在线观看| 窝窝影院91人妻| 女人高潮潮喷娇喘18禁视频| 亚洲欧美一区二区三区黑人| 欧美中文综合在线视频| 成人特级黄色片久久久久久久| 在线天堂中文资源库| 亚洲欧美激情综合另类| 久久久久国产精品人妻aⅴ院| 两个人免费观看高清视频| 一区二区三区精品91| 日日摸夜夜添夜夜添小说| 人人妻人人澡欧美一区二区 | 国产99久久九九免费精品| 亚洲精品在线观看二区| 大型av网站在线播放| 久久影院123| 18禁观看日本| 免费看美女性在线毛片视频| www.精华液| 亚洲一区二区三区色噜噜| 精品久久久久久久毛片微露脸| 国产精品一区二区免费欧美| 久久婷婷成人综合色麻豆| 亚洲人成77777在线视频| 亚洲第一欧美日韩一区二区三区| 777久久人妻少妇嫩草av网站| 深夜精品福利| 久99久视频精品免费| 国产精品亚洲美女久久久| 啦啦啦免费观看视频1| 丝袜人妻中文字幕| 精品高清国产在线一区| 久久国产乱子伦精品免费另类| 免费女性裸体啪啪无遮挡网站| 自拍欧美九色日韩亚洲蝌蚪91| 黄色成人免费大全| 亚洲,欧美精品.| 国产欧美日韩综合在线一区二区| 国产成人精品久久二区二区免费| 村上凉子中文字幕在线| 99国产精品一区二区蜜桃av| 国产人伦9x9x在线观看| 国产99白浆流出| 国产一区二区三区在线臀色熟女| 亚洲第一欧美日韩一区二区三区| 亚洲色图av天堂| 99国产精品免费福利视频| 国产麻豆69| 免费在线观看视频国产中文字幕亚洲| 亚洲国产毛片av蜜桃av| 免费看十八禁软件| 一级毛片高清免费大全| 亚洲欧美日韩无卡精品| www.www免费av| 91精品国产国语对白视频| 天天躁狠狠躁夜夜躁狠狠躁| 老鸭窝网址在线观看| 一级毛片女人18水好多| 欧美日本亚洲视频在线播放| 国产精品久久久久久人妻精品电影| 9色porny在线观看| 亚洲 国产 在线| 高清黄色对白视频在线免费看| 国产人伦9x9x在线观看| 老熟妇仑乱视频hdxx| 中文字幕另类日韩欧美亚洲嫩草| 亚洲最大成人中文| 亚洲成人久久性| 18禁美女被吸乳视频| 国产精品一区二区在线不卡| 国产一级毛片七仙女欲春2 | 亚洲av美国av| 欧美人与性动交α欧美精品济南到| 久久国产精品人妻蜜桃| 老熟妇仑乱视频hdxx| 国产亚洲精品久久久久5区| 国产av一区在线观看免费| 大型黄色视频在线免费观看| 欧美成人性av电影在线观看| 精品一区二区三区av网在线观看| 免费看a级黄色片| 日韩高清综合在线| 成人精品一区二区免费| 在线视频色国产色| 国产精华一区二区三区| 老司机深夜福利视频在线观看| 黄片小视频在线播放| 日本撒尿小便嘘嘘汇集6| 亚洲精品一卡2卡三卡4卡5卡| 亚洲熟妇熟女久久| 国产成人欧美在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久久国产精品麻豆| 中文字幕最新亚洲高清| 亚洲欧美日韩无卡精品| 两性夫妻黄色片| 搡老熟女国产l中国老女人| 日韩欧美三级三区| 欧美黄色淫秽网站| 国产午夜福利久久久久久| 香蕉丝袜av| 我的亚洲天堂| 美女大奶头视频| 好男人在线观看高清免费视频 | 精品久久久久久久人妻蜜臀av | 免费观看人在逋| 国产99久久九九免费精品| 久久狼人影院| 精品一区二区三区四区五区乱码| av免费在线观看网站| 一个人免费在线观看的高清视频| av天堂久久9| 村上凉子中文字幕在线| 久久性视频一级片| 国产精品国产高清国产av| 亚洲av第一区精品v没综合| 亚洲avbb在线观看| 女人被躁到高潮嗷嗷叫费观| 精品免费久久久久久久清纯| 日韩大尺度精品在线看网址 | 久久香蕉激情| 黄色成人免费大全| 日韩视频一区二区在线观看| 国产aⅴ精品一区二区三区波| 97碰自拍视频| 男人的好看免费观看在线视频 | 国产精品av久久久久免费| 久久国产精品人妻蜜桃| 国产xxxxx性猛交| 18禁黄网站禁片午夜丰满| 岛国在线观看网站| 久久香蕉国产精品| 日韩大尺度精品在线看网址 | 女人被狂操c到高潮| 国产精品综合久久久久久久免费 | 91麻豆av在线| 久久欧美精品欧美久久欧美| 一区在线观看完整版| 国语自产精品视频在线第100页| 欧美日韩乱码在线| 一级毛片高清免费大全| 欧美黄色片欧美黄色片| 法律面前人人平等表现在哪些方面| 色在线成人网| 国产麻豆成人av免费视频| 国产一级毛片七仙女欲春2 | 免费搜索国产男女视频| 久久精品亚洲精品国产色婷小说| 久久午夜综合久久蜜桃| av欧美777| 亚洲第一欧美日韩一区二区三区| 日本免费一区二区三区高清不卡 | 老熟妇乱子伦视频在线观看| 极品教师在线免费播放| 免费一级毛片在线播放高清视频 | 自拍欧美九色日韩亚洲蝌蚪91| 久久久久亚洲av毛片大全| 一区在线观看完整版| 一区二区三区激情视频| 最近最新免费中文字幕在线| 欧美日韩瑟瑟在线播放| 国内毛片毛片毛片毛片毛片| 一区二区三区国产精品乱码| 免费高清视频大片| av天堂久久9| 国产精品久久久人人做人人爽| 一卡2卡三卡四卡精品乱码亚洲| 久久国产乱子伦精品免费另类| 国产精品美女特级片免费视频播放器 | 男人的好看免费观看在线视频 | 色播亚洲综合网| 一本综合久久免费| 日韩视频一区二区在线观看| 一区二区三区国产精品乱码| 日韩三级视频一区二区三区| 怎么达到女性高潮| 亚洲av成人av| 制服人妻中文乱码| 99re在线观看精品视频| 可以在线观看的亚洲视频| 老熟妇乱子伦视频在线观看| 一二三四社区在线视频社区8| 夜夜夜夜夜久久久久| 黄网站色视频无遮挡免费观看| 美女高潮喷水抽搐中文字幕| 免费不卡黄色视频| 精品国产亚洲在线| 搡老妇女老女人老熟妇| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| aaaaa片日本免费| 妹子高潮喷水视频| 这个男人来自地球电影免费观看| 侵犯人妻中文字幕一二三四区| 亚洲精品久久国产高清桃花| 少妇 在线观看| 亚洲一码二码三码区别大吗| 欧美日韩瑟瑟在线播放| а√天堂www在线а√下载| 久久精品亚洲熟妇少妇任你| 99国产综合亚洲精品| 搡老熟女国产l中国老女人| 久久香蕉国产精品| 日日爽夜夜爽网站| 视频区欧美日本亚洲| 国产亚洲精品综合一区在线观看 | 天堂影院成人在线观看| 在线永久观看黄色视频| 天堂影院成人在线观看| 欧美激情 高清一区二区三区| 国产精品一区二区在线不卡| 日韩欧美三级三区| 久久国产精品人妻蜜桃| 1024香蕉在线观看| 制服人妻中文乱码| 美女高潮喷水抽搐中文字幕| 日韩大尺度精品在线看网址 | www.熟女人妻精品国产| 精品国产美女av久久久久小说| 天堂影院成人在线观看| 深夜精品福利| 啦啦啦韩国在线观看视频| 国产三级黄色录像| 色综合婷婷激情| 国产av精品麻豆| 老司机深夜福利视频在线观看| 国产高清有码在线观看视频 | avwww免费| 午夜久久久在线观看| 免费av毛片视频| 亚洲情色 制服丝袜| 久久精品影院6| 热re99久久国产66热|