• <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例報告)
    韩国高清视频一区二区三区| 国产激情久久老熟女| 中文字幕另类日韩欧美亚洲嫩草| 69av精品久久久久久 | h视频一区二区三区| 美女午夜性视频免费| 日韩电影二区| 18禁国产床啪视频网站| 在线观看一区二区三区激情| 午夜免费观看性视频| 日韩精品免费视频一区二区三区| 男女床上黄色一级片免费看| 欧美黑人欧美精品刺激| 中文字幕人妻丝袜一区二区| 日韩熟女老妇一区二区性免费视频| 国产伦人伦偷精品视频| 99国产综合亚洲精品| 亚洲精品粉嫩美女一区| 成人免费观看视频高清| 大陆偷拍与自拍| 亚洲五月婷婷丁香| 免费观看a级毛片全部| 亚洲国产中文字幕在线视频| 亚洲综合色网址| 桃花免费在线播放| 桃花免费在线播放| 老司机午夜福利在线观看视频 | 亚洲av国产av综合av卡| 欧美人与性动交α欧美精品济南到| 99热国产这里只有精品6| 日本一区二区免费在线视频| 国产97色在线日韩免费| 十八禁人妻一区二区| a在线观看视频网站| 免费高清在线观看视频在线观看| 超色免费av| 国产熟女午夜一区二区三区| 汤姆久久久久久久影院中文字幕| 国产一区二区三区av在线| 又大又爽又粗| 永久免费av网站大全| 久久精品国产亚洲av高清一级| 汤姆久久久久久久影院中文字幕| 黄片大片在线免费观看| 精品一区二区三区四区五区乱码| 国产一区二区在线观看av| 国产麻豆69| 成年女人毛片免费观看观看9 | 久久精品国产亚洲av香蕉五月 | 亚洲自偷自拍图片 自拍| 免费在线观看日本一区| 精品久久久精品久久久| 亚洲三区欧美一区| 动漫黄色视频在线观看| 黄色 视频免费看| 国产日韩欧美亚洲二区| 少妇被粗大的猛进出69影院| av在线app专区| 午夜福利,免费看| 亚洲精品美女久久久久99蜜臀| 日本精品一区二区三区蜜桃| 欧美激情极品国产一区二区三区| 亚洲精品久久久久久婷婷小说| 两人在一起打扑克的视频| av国产精品久久久久影院| 在线 av 中文字幕| 九色亚洲精品在线播放| 热re99久久精品国产66热6| 国产亚洲av片在线观看秒播厂| 欧美激情 高清一区二区三区| 天堂中文最新版在线下载| 人成视频在线观看免费观看| 久久久国产成人免费| 精品一区二区三区四区五区乱码| 久久久国产一区二区| 日日摸夜夜添夜夜添小说| 亚洲va日本ⅴa欧美va伊人久久 | www.自偷自拍.com| 搡老岳熟女国产| 色94色欧美一区二区| 性色av乱码一区二区三区2| 丝袜美腿诱惑在线| 性色av乱码一区二区三区2| 亚洲九九香蕉| 狠狠婷婷综合久久久久久88av| 国产伦人伦偷精品视频| 性色av乱码一区二区三区2| 亚洲精品日韩在线中文字幕| 国产一区二区在线观看av| 窝窝影院91人妻| 久久国产亚洲av麻豆专区| 中文字幕高清在线视频| 日韩中文字幕欧美一区二区| 男女午夜视频在线观看| h视频一区二区三区| 欧美激情极品国产一区二区三区| www.av在线官网国产| 肉色欧美久久久久久久蜜桃| 国产精品 国内视频| 精品熟女少妇八av免费久了| 久久久国产一区二区| 国产亚洲精品一区二区www | 69av精品久久久久久 | 母亲3免费完整高清在线观看| 亚洲天堂av无毛| 视频区图区小说| 久久精品久久久久久噜噜老黄| 男人操女人黄网站| 国产男女内射视频| 午夜免费鲁丝| 亚洲精品成人av观看孕妇| 18禁观看日本| 国产欧美日韩一区二区三 | 满18在线观看网站| 欧美日韩成人在线一区二区| 一区二区av电影网| 91麻豆av在线| 超碰成人久久| 国产精品免费视频内射| 纯流量卡能插随身wifi吗| 好男人电影高清在线观看| 国产主播在线观看一区二区| 免费高清在线观看日韩| 亚洲精华国产精华精| 欧美日韩一级在线毛片| 日本撒尿小便嘘嘘汇集6| 亚洲国产欧美在线一区| 日韩免费高清中文字幕av| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黑人欧美精品刺激| 肉色欧美久久久久久久蜜桃| 国产91精品成人一区二区三区 | 在线观看免费视频网站a站| 亚洲欧美精品自产自拍| 女警被强在线播放| 中国国产av一级| 在线十欧美十亚洲十日本专区| 啦啦啦免费观看视频1| 国产免费现黄频在线看| 久久免费观看电影| 亚洲伊人久久精品综合| 91老司机精品| 曰老女人黄片| 欧美日韩亚洲高清精品| 91成人精品电影| 成年人免费黄色播放视频| av超薄肉色丝袜交足视频| e午夜精品久久久久久久| 啦啦啦在线免费观看视频4| 少妇 在线观看| 日韩视频在线欧美| 久久中文字幕一级| 国产熟女午夜一区二区三区| 国产精品一区二区在线不卡| 美女国产高潮福利片在线看| 丁香六月天网| 午夜福利,免费看| 热99re8久久精品国产| 欧美国产精品va在线观看不卡| 热99久久久久精品小说推荐| e午夜精品久久久久久久| 欧美精品亚洲一区二区| 欧美日本中文国产一区发布| 91麻豆av在线| 纯流量卡能插随身wifi吗| 亚洲欧美精品自产自拍| 日日夜夜操网爽| 国产精品久久久av美女十八| 在线观看免费日韩欧美大片| 两性午夜刺激爽爽歪歪视频在线观看 | 免费黄频网站在线观看国产| 中亚洲国语对白在线视频| 青青草视频在线视频观看| 国产深夜福利视频在线观看| 桃红色精品国产亚洲av| 一本综合久久免费| 少妇被粗大的猛进出69影院| 一区二区av电影网| 建设人人有责人人尽责人人享有的| videosex国产| 久久影院123| 精品亚洲乱码少妇综合久久| 久久精品亚洲av国产电影网| 国产老妇伦熟女老妇高清| 亚洲第一青青草原| 91麻豆av在线| 丝瓜视频免费看黄片| 纵有疾风起免费观看全集完整版| 色播在线永久视频| 久久毛片免费看一区二区三区| 99热国产这里只有精品6| 岛国在线观看网站| 精品久久久久久久毛片微露脸 | 好男人电影高清在线观看| 日本a在线网址| 国产精品一区二区在线不卡| 亚洲精品国产av蜜桃| 大陆偷拍与自拍| 久久久国产一区二区| 亚洲成人免费电影在线观看| 最近中文字幕2019免费版| av网站在线播放免费| a级毛片黄视频| 91麻豆精品激情在线观看国产 | 久久久久久久久免费视频了| 搡老岳熟女国产| 亚洲精华国产精华精| 久久久久久久精品精品| 日韩 欧美 亚洲 中文字幕| 亚洲人成电影观看| 天天躁日日躁夜夜躁夜夜| 一区福利在线观看| 一区二区三区精品91| 国产1区2区3区精品| 国产亚洲av片在线观看秒播厂| 亚洲国产欧美一区二区综合| 大片免费播放器 马上看| 在线观看免费视频网站a站| 桃红色精品国产亚洲av| 午夜福利视频精品| av视频免费观看在线观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美精品自产自拍| 黄片播放在线免费| 国产日韩欧美在线精品| 如日韩欧美国产精品一区二区三区| 亚洲av成人一区二区三| 宅男免费午夜| 亚洲伊人久久精品综合| 亚洲中文字幕日韩| 免费观看av网站的网址| 亚洲精品久久午夜乱码| 欧美av亚洲av综合av国产av| 国产欧美日韩一区二区精品| 国产老妇伦熟女老妇高清| 国产一区二区三区综合在线观看| 亚洲av男天堂| 丝袜在线中文字幕| cao死你这个sao货| 在线观看免费高清a一片| 午夜福利视频在线观看免费| 男女下面插进去视频免费观看| 黄片播放在线免费| 一区福利在线观看| 99久久精品国产亚洲精品| 少妇 在线观看| 黄色视频,在线免费观看| 高清在线国产一区| 十八禁高潮呻吟视频| 母亲3免费完整高清在线观看| 两个人看的免费小视频| 超碰97精品在线观看| 亚洲成人国产一区在线观看| 国产免费一区二区三区四区乱码| 亚洲国产欧美日韩在线播放| 久久天堂一区二区三区四区| 亚洲熟女精品中文字幕| 我的亚洲天堂| 午夜精品国产一区二区电影| 黄色视频,在线免费观看| 在线av久久热| 老司机亚洲免费影院| 国产又爽黄色视频| 国产福利在线免费观看视频| 欧美日韩福利视频一区二区| 亚洲国产欧美网| 国产日韩欧美亚洲二区| 国产无遮挡羞羞视频在线观看| 99精品欧美一区二区三区四区| 美女主播在线视频| 99久久综合免费| 国产亚洲精品第一综合不卡| 亚洲精品中文字幕一二三四区 | 99热网站在线观看| 在线十欧美十亚洲十日本专区| 国产亚洲av片在线观看秒播厂| 高清黄色对白视频在线免费看| 热re99久久国产66热| 国产伦理片在线播放av一区| 国产高清国产精品国产三级| 国产免费av片在线观看野外av| 99热网站在线观看| 精品亚洲乱码少妇综合久久| 三级毛片av免费| 天堂俺去俺来也www色官网| 国产福利在线免费观看视频| 欧美国产精品va在线观看不卡| 国产成人精品在线电影| 一级毛片精品| 亚洲精品美女久久久久99蜜臀| 精品一区在线观看国产| 午夜福利一区二区在线看| 首页视频小说图片口味搜索| 极品少妇高潮喷水抽搐| 国产一区二区在线观看av| 蜜桃国产av成人99| 青草久久国产| 岛国在线观看网站| 久久中文字幕一级| 午夜福利影视在线免费观看| 老司机影院毛片| 久久久久久久精品精品| 亚洲,欧美精品.| 多毛熟女@视频| 国产又爽黄色视频| 精品欧美一区二区三区在线| 午夜影院在线不卡| 婷婷丁香在线五月| 搡老岳熟女国产| 亚洲av男天堂| 乱人伦中国视频| 99精品久久久久人妻精品| 丰满饥渴人妻一区二区三| 久久天躁狠狠躁夜夜2o2o| 女性被躁到高潮视频| 国产精品自产拍在线观看55亚洲 | 亚洲av电影在线观看一区二区三区| 啦啦啦免费观看视频1| 国产区一区二久久| 精品亚洲成a人片在线观看| 精品熟女少妇八av免费久了| 999久久久精品免费观看国产| 18在线观看网站| 九色亚洲精品在线播放| 久久久久久久久久久久大奶| 国产精品一区二区精品视频观看| 交换朋友夫妻互换小说| 欧美精品亚洲一区二区| 国产精品 国内视频| 五月天丁香电影| 国产亚洲精品久久久久5区| 在线观看舔阴道视频| 91成年电影在线观看| 国产1区2区3区精品| 欧美日本中文国产一区发布| 久久 成人 亚洲| 99re6热这里在线精品视频| 久久影院123| 91麻豆精品激情在线观看国产 | 成人亚洲精品一区在线观看| 国产成人a∨麻豆精品| 精品人妻熟女毛片av久久网站| 亚洲国产精品一区三区| 中文字幕人妻丝袜一区二区| 精品国内亚洲2022精品成人 | 亚洲精品中文字幕在线视频| 亚洲美女黄色视频免费看| 久久人人爽人人片av| 欧美精品人与动牲交sv欧美| 国产精品免费视频内射| 欧美日韩成人在线一区二区| 亚洲精品国产区一区二| 欧美日韩福利视频一区二区| av线在线观看网站| 国产97色在线日韩免费| 免费少妇av软件| 可以免费在线观看a视频的电影网站| 亚洲国产中文字幕在线视频| 97在线人人人人妻| 中文字幕色久视频| 久久狼人影院| 欧美另类亚洲清纯唯美| 久热这里只有精品99| av网站在线播放免费| 欧美97在线视频| 性少妇av在线| 国产精品久久久av美女十八| 人妻久久中文字幕网| 免费少妇av软件| 十分钟在线观看高清视频www| 国产野战对白在线观看| 少妇的丰满在线观看| 在线观看人妻少妇| 成人三级做爰电影| 国产精品久久久人人做人人爽| 男女国产视频网站| 人人妻人人爽人人添夜夜欢视频| 人人妻人人添人人爽欧美一区卜| 99re6热这里在线精品视频| 午夜免费鲁丝| 亚洲av日韩精品久久久久久密| 午夜福利视频精品| 久久人妻福利社区极品人妻图片| 国产老妇伦熟女老妇高清| 人成视频在线观看免费观看| 亚洲五月色婷婷综合| 欧美亚洲 丝袜 人妻 在线| 久久国产精品影院| 免费人妻精品一区二区三区视频| 亚洲国产欧美日韩在线播放| 亚洲国产av影院在线观看| 国产深夜福利视频在线观看| 日韩大片免费观看网站| 天天影视国产精品| 久久久精品国产亚洲av高清涩受| 天堂俺去俺来也www色官网| 男女边摸边吃奶| 超碰成人久久| 无限看片的www在线观看| 少妇精品久久久久久久| 亚洲九九香蕉| 欧美精品av麻豆av| 午夜福利视频在线观看免费| a级毛片在线看网站| 欧美 亚洲 国产 日韩一| 国产有黄有色有爽视频| 精品国产国语对白av| 啦啦啦 在线观看视频| 色视频在线一区二区三区| 国产av一区二区精品久久| 久久久久视频综合| 午夜视频精品福利| 中国美女看黄片| 国产成+人综合+亚洲专区| 国产在线免费精品| www.熟女人妻精品国产| 少妇精品久久久久久久| av福利片在线| 亚洲欧美精品综合一区二区三区| 免费高清在线观看日韩| 精品国内亚洲2022精品成人 | 黄片大片在线免费观看| 在线 av 中文字幕| 91麻豆精品激情在线观看国产 | 女人爽到高潮嗷嗷叫在线视频| 午夜福利乱码中文字幕| 国产成人免费观看mmmm| 亚洲精品国产区一区二| 国产精品 国内视频| 热99国产精品久久久久久7| 成人免费观看视频高清| 青青草视频在线视频观看| 在线观看免费视频网站a站| 国产av一区二区精品久久| 首页视频小说图片口味搜索| 亚洲专区中文字幕在线| 亚洲五月婷婷丁香| 日韩 亚洲 欧美在线| 黄频高清免费视频| 一级毛片女人18水好多| 亚洲精品一二三| av一本久久久久| 久久天堂一区二区三区四区| 日日夜夜操网爽| 精品久久久久久久毛片微露脸 | 国产一区二区三区在线臀色熟女 | 免费人妻精品一区二区三区视频| 国产av精品麻豆| 搡老熟女国产l中国老女人| 狠狠狠狠99中文字幕| 又紧又爽又黄一区二区| 午夜精品久久久久久毛片777| 乱人伦中国视频| 中文字幕色久视频| 国产一卡二卡三卡精品| 国产成人a∨麻豆精品| 精品国产乱码久久久久久男人| 下体分泌物呈黄色| 成人手机av| 自线自在国产av| 成人国语在线视频| 热99久久久久精品小说推荐| 啦啦啦中文免费视频观看日本| 最近中文字幕2019免费版| 丝袜脚勾引网站| 十八禁人妻一区二区| 成年av动漫网址| 18禁裸乳无遮挡动漫免费视频| 久久精品国产亚洲av香蕉五月 | 丰满人妻熟妇乱又伦精品不卡| 亚洲,欧美精品.| 国产1区2区3区精品| 精品免费久久久久久久清纯 | 久久久久网色| 国产黄频视频在线观看| 天天躁夜夜躁狠狠躁躁| 大片电影免费在线观看免费| 十八禁网站免费在线| 精品卡一卡二卡四卡免费| 亚洲一码二码三码区别大吗| 国产亚洲精品一区二区www | av免费在线观看网站| 午夜福利免费观看在线| 黄色毛片三级朝国网站| 国产在线免费精品| 国产日韩一区二区三区精品不卡| 国产精品国产av在线观看| 别揉我奶头~嗯~啊~动态视频 | 啪啪无遮挡十八禁网站| 精品久久久精品久久久| 国产精品.久久久| 国产1区2区3区精品| 成年人免费黄色播放视频| 看免费av毛片| 99国产精品免费福利视频| 制服人妻中文乱码| h视频一区二区三区| av网站在线播放免费| 亚洲国产看品久久| 欧美一级毛片孕妇| 欧美av亚洲av综合av国产av| 两人在一起打扑克的视频| 国产在线一区二区三区精| 一进一出抽搐动态| 国产一级毛片在线| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 国产极品粉嫩免费观看在线| videos熟女内射| 欧美少妇被猛烈插入视频| 国产亚洲欧美在线一区二区| 老司机亚洲免费影院| 日本av免费视频播放| 亚洲成人免费av在线播放| 法律面前人人平等表现在哪些方面 | 狂野欧美激情性xxxx| 午夜福利一区二区在线看| 亚洲 欧美一区二区三区| 亚洲欧美精品自产自拍| 动漫黄色视频在线观看| 少妇的丰满在线观看| 美女高潮喷水抽搐中文字幕| 一区二区三区乱码不卡18| 又大又爽又粗| 人人妻,人人澡人人爽秒播| www.熟女人妻精品国产| 久久精品人人爽人人爽视色| 一级a爱视频在线免费观看| 精品国产乱码久久久久久小说| 日本欧美视频一区| 亚洲精品久久午夜乱码| 51午夜福利影视在线观看| 婷婷色av中文字幕| 成人国产av品久久久| 一个人免费在线观看的高清视频 | 狠狠婷婷综合久久久久久88av| 国产视频一区二区在线看| 操出白浆在线播放| 国产成人免费无遮挡视频| 天堂中文最新版在线下载| 亚洲精品中文字幕在线视频| 成人三级做爰电影| 国产三级黄色录像| 水蜜桃什么品种好| 视频区图区小说| 波多野结衣av一区二区av| 亚洲第一欧美日韩一区二区三区 | 亚洲人成77777在线视频| 国产日韩一区二区三区精品不卡| a级片在线免费高清观看视频| 欧美少妇被猛烈插入视频| 亚洲第一青青草原| 91老司机精品| 岛国毛片在线播放| 91精品国产国语对白视频| 俄罗斯特黄特色一大片| 一本大道久久a久久精品| 精品免费久久久久久久清纯 | 一区二区av电影网| 无遮挡黄片免费观看| 久久久久视频综合| 一区福利在线观看| 亚洲欧洲精品一区二区精品久久久| 高清视频免费观看一区二区| 一级毛片女人18水好多| 国产一区二区三区在线臀色熟女 | 最新的欧美精品一区二区| 亚洲av成人不卡在线观看播放网 | 91成年电影在线观看| 少妇 在线观看| 国产97色在线日韩免费| 一本一本久久a久久精品综合妖精| 久久精品亚洲熟妇少妇任你| 久久精品久久久久久噜噜老黄| 午夜免费观看性视频| netflix在线观看网站| 免费高清在线观看视频在线观看| 亚洲国产毛片av蜜桃av| 精品国内亚洲2022精品成人 | 他把我摸到了高潮在线观看 | 午夜影院在线不卡| 国产99久久九九免费精品| 日韩人妻精品一区2区三区| 女人高潮潮喷娇喘18禁视频| 9色porny在线观看| 69av精品久久久久久 | 国产深夜福利视频在线观看| 欧美黑人欧美精品刺激| 首页视频小说图片口味搜索| 国产淫语在线视频| 亚洲精品国产一区二区精华液| 国产一卡二卡三卡精品| 亚洲情色 制服丝袜| 90打野战视频偷拍视频| 女人爽到高潮嗷嗷叫在线视频| av在线老鸭窝| 国产精品免费大片| 亚洲成人国产一区在线观看| 国精品久久久久久国模美| 久久精品国产亚洲av高清一级| 国产精品久久久久久精品古装| 两人在一起打扑克的视频| 在线永久观看黄色视频| 窝窝影院91人妻| 黄频高清免费视频| 秋霞在线观看毛片| 99久久国产精品久久久| 欧美97在线视频| 天堂8中文在线网| 黑丝袜美女国产一区| 国产精品免费视频内射| 在线观看免费午夜福利视频| 亚洲国产毛片av蜜桃av| 制服人妻中文乱码| 一个人免费在线观看的高清视频 |