两个人的电影免费视频_国产精品久久久久久久久成人_97视频在线观看播放_久久这里只有精品777_亚洲熟女少妇二三区_4438x8成人网亚洲av_内谢国产内射夫妻免费视频_人妻精品久久久久中国字幕

三維大腦磁共振圖像大腦皮層表面最大主方向場彌散方法

文檔序號:1184709閱讀:287來源:國知局
專利名稱:三維大腦磁共振圖像大腦皮層表面最大主方向場彌散方法
技術(shù)領(lǐng)域
本發(fā)明涉及一種三維大腦磁共振圖像大腦皮層表面最大主方向場彌散方法,利用 alpha-expansion圖割(graph cuts)方法的三維大腦磁共振圖像的大腦皮層表面上最大 主方向場彌散方法,屬于醫(yī)學(xué)圖像處理領(lǐng)域。
背景技術(shù)
從三維大腦磁共振圖像重構(gòu)出來的大腦皮層表面是嵌入在三維空間的、極其卷曲 和褶皺的二維流型?;诖竽X皮層表面的分析非常適合于刻畫和分析大腦皮層的本質(zhì)結(jié) 構(gòu),例如曲率、測地距離、皮層厚度等。大腦皮層表面上的方向場包含重要的信息,例如,最 大主方向場和最小主方向場分別指示主曲率達到最大和最小的方向。但是,可靠計算大腦 皮層表面上的主方向場并不容易,例如,在平坦的皮層表面區(qū)域所計算的主方向場通常噪 聲較大和不可靠,因為在平坦區(qū)域曲率在各個方向上的變化不大,微小的噪聲就可能嚴(yán)重 影響所計算的主方向場。因此,為了利用主方向場對大腦皮層表面進行分析,必須實施主方 向場彌散,其目的是平滑主方向場中的噪聲,同時保持主方向場中主要的幾何結(jié)構(gòu)和不連 續(xù)性。彌散后的主方向場可以用于基于流場跟蹤的腦溝或腦回盆地分割和方向?qū)虻钠?褶皺形態(tài)分析等。目前已有的大腦皮層表面上的最大主方向場彌散方法將該問題看做一個能量最 小化問題,該能量函數(shù)包含一個數(shù)據(jù)項和一個平滑項。數(shù)據(jù)項用于使彌散后的最大主方向 場與原始計算的最大主方向場盡可能相似,平滑項用于使彌散后的最大主方向場平滑地變 化。該能量函數(shù)通過變分方法或梯度下降方式求解。已有的大腦皮層表面上的最大主方向 場彌散方法具有以下的主要缺陷其一、不同皮層表面區(qū)域的權(quán)重設(shè)為相同,導(dǎo)致彌散的最 大主方向場在腦溝底部和腦回頂部過于平滑。其二、目前的能量最小化方法容易陷入到低 劣的局部最優(yōu)。

發(fā)明內(nèi)容
要解決的技術(shù)問題為了避免現(xiàn)有方法的不足之處,本發(fā)明提出一種三維大腦磁共振圖像大腦皮層表 面最大主方向場彌散方法,適用于人類三維大腦磁共振圖像重構(gòu)出的、由三角形表示的大 腦皮層表面上最大主方向場彌散。技術(shù)方案本發(fā)明的基本思想是將大腦皮層表面上的最大主方向場彌散問題公式化為一個 能量最小化問題,進而轉(zhuǎn)化為一個離散標(biāo)記問題,即給大腦皮層表面上的每一個頂點賦予 一個標(biāo)記來表示彌散的最大主方向,該離散標(biāo)記問題可以利用alpha-expansion圖割方法 有效的求解,因為圖割方法可以保證得到某些類型的能量函數(shù)的強局部最優(yōu)。一種三維大腦磁共振圖像的大腦皮層表面上最大主方向場彌散方法,其特征在于 步驟如下
3
步驟1對三維大腦磁共振圖像進行預(yù)處理和大腦皮層表面重建利用變形模型方 法去除腦殼,利用配準(zhǔn)方法去除非大腦組織,利用高斯混合模型和馬爾科夫隨機場方法對 圖像進行組織分割,得到白質(zhì)、灰質(zhì)和腦脊髓液三種組織類型表示的圖像,利用基于圖的方 法對白質(zhì)圖像進行拓撲校正,利用Marching Cubes方法重建由三角形表示的大腦皮層表步驟2計算大腦皮層表面上頂點的主曲率、主方向以及主曲率的導(dǎo)數(shù)利用有限 差分方法計算大腦皮層表面上頂點的主曲率、主方向以及主曲率的導(dǎo)數(shù),當(dāng)每個頂點的最 大主曲率在最大主方向的導(dǎo)數(shù)為正值時,將該頂點的最大主方向改為相反方向;步驟3對最大主方向場進行彌散利用alpha-expansion圖割方法最小化能量函 數(shù)
腦皮層表面上的頂點坐標(biāo),p(x)為頂點X處原始的最大主方向;(X,y) G N表示X和y為 相鄰頂點的集合,S為大腦皮層表面上的所有頂點的集合;U | |為2范數(shù)為頂點x彌 散的最大主方向,lx為頂點x彌散的最大主方向在解空間的標(biāo)識且lx e L ;^為頂點y彌 散的最大主方向,ly為頂點y彌散的最大主方向在解空間的標(biāo)識且ly e L ;所述的L = {l1,I2,. . .,ln}為離散的彌散的最大主方向場的解空間 = {v1, v2,... ,vn}的標(biāo)識,其中n = n0 X (n,-2)+2 ;ne為x_y平面中離散角度的數(shù)目為12 36 ; n,為z軸方向中離散角度的數(shù)目為9 18 ;所述解空間 中7= (0,0,1), vn= (0,0,-1); 所述g(x) = exp("A c(x) |),所述 g(y) = exp (-A c (y) |),所述 h(x)= l-g(x),其中\(zhòng)為權(quán)重參數(shù)為5.0 10.0,c(0為最大主曲率;利用alpha-expansion圖割方法求解,得到彌散的最大主方向場;步驟4將彌散的最大主方向場投影到每個頂點的切平面內(nèi)。有益效果本發(fā)明提出的三維大腦磁共振圖像大腦皮層表面最大主方向場彌散方法,使得三 維大腦磁共振圖像中大腦皮層表面上最大主方向場彌散方法的可行性體現(xiàn)在首先,隨著 磁共振成像設(shè)備的精度不斷提高和三維大腦磁共振圖像的預(yù)處理方法進一步成熟,獲取幾 何結(jié)構(gòu)精確、拓撲結(jié)構(gòu)正確的大腦皮層表面相對容易;同時,將最大主方向場彌散看做一個 能量最小化問題,并轉(zhuǎn)換為一個離散標(biāo)記問題,利用圖割方法有效求解該離散標(biāo)記問題是 可行的。本發(fā)明相對于其它方法具有以下優(yōu)點1、該方法在大腦皮層表面上不同區(qū)域的平 滑項和數(shù)據(jù)項的權(quán)重設(shè)為不同,這樣可以更好的保持最大主方向場中固有的幾何結(jié)構(gòu)和不 連續(xù)性;2、將最大主方向場彌散看做一個能量最小化問題,并轉(zhuǎn)換為一個離散標(biāo)記問題,利 用圖割方法可以有效求解該能量函數(shù)的強局部最優(yōu)。


圖1 本發(fā)明實施例的三維大腦磁共振圖像大腦皮層表面上最大主方向場彌散,(a)顯示了實施例圖像中大腦皮層表面上的最大主曲率,(b)顯示了圖1(a)中矩形區(qū)域中大腦皮層表面上的最大主方向場,(c)顯示了圖1(b)中對應(yīng)的大腦皮層表面上的彌散的最大主方向場;圖2 仿真表面上的最大主方向場彌散的一個例子,(a)顯示了三維大腦磁共振圖像仿真表面上的理想最大主方向場,(b)顯示了三維大腦磁共振圖像仿真噪聲的最大主方向場,(c)顯示了三維大腦磁共振圖像仿真彌散的最大主方向;圖3 三維大腦磁共振圖像大腦皮層表面上的最大主方向場彌散前和彌散后的一 致性比較實例,(a)顯示實施例中三維大腦磁共振圖像大腦皮層表面上的最大主方向場的一致性(b)顯示實施例中三維大腦磁共振圖像大腦皮層表面上彌散的最大主方向場的一 致性;圖4 本實例大腦皮層表面上的最大主方向場彌散前和彌散后的一致性的分布直 方圖;圖5 :12個大腦皮層表面上的最大主方向場彌散前和彌散后的平均一致性的比 較;
具體實施例方式現(xiàn)結(jié)合實施例、附圖對本發(fā)明作進一步描述根據(jù)本發(fā)明提出的基于圖割方法的大腦皮層表面上的最大主方向場彌散方法,我 們用C++語言實現(xiàn)了一個該方法的原型系統(tǒng)。圖像數(shù)據(jù)的來源是實際中正常人的三維大 腦磁共振圖像。具體的實施步驟如下1.預(yù)處理和大腦皮層表面重建利用變形模型方法去除腦殼,利用配準(zhǔn)方法去除非大腦組織,利用高斯混合模型 和馬爾科夫隨機場方法對圖像進行組織分割,得到白質(zhì)、灰質(zhì)和腦脊髓液三種組織類型表 示的圖像,利用基于圖的方法對白質(zhì)圖像進行拓撲校正,利用Marching Cubes方法重建由 三角形表示的大腦皮層表面;2.計算大腦皮層表面上頂點的主曲率、主方向以及主曲率的導(dǎo)數(shù)利用有限差分方法計算大腦皮層表面上頂點的主曲率、主方向以及主曲率的導(dǎo) 數(shù),檢驗每個頂點的最大主曲率在最大主方向的導(dǎo)數(shù)是否為正值,如果為正值,則翻轉(zhuǎn)該頂 點的最大主方向為相反方向;圖1(a)顯示了一例大腦皮層表面上的最大主曲率;圖1(b) 顯示了圖1(a)中矩形區(qū)域中大腦皮層表面上的最大主方向場。3.對最大主方向場進行彌散大腦皮層表面上最大主方向場的彌散可以表示為如下的能量函數(shù)
5 其中,v(x) = (u(x), v(x), w(x))是頂點x處彌散的最大主方向,v(y)是頂點y 處彌散的最大主方向,P(X)是頂點x處原始的最大主方向。(x,y) GN表示X和y為相 鄰頂點的集合,S表示大腦皮層表面上的所有頂點的集合。| | | |表示2范數(shù)。g(x)= exp (- A c(x) ),g(y) = exp (- A c (y) ), h(x) = l_g(x),其中入為權(quán)重參數(shù)為 8.0,c(.)為最大主曲率。在大腦皮層表面上強烈彎曲的區(qū)域,對應(yīng)于腦溝底部和腦回冠 部區(qū)域,該處c( )的絕對值較大,最大主方向場可以相對可靠的計算;而在平坦的區(qū)域所 計算的最大主方向場的噪聲較大,比較不可靠。因此g( )設(shè)為在腦溝底部和腦回冠部的
值較小,而h(_)在腦溝底部和腦回冠部的值較大。!>(x)lv(x) —P(x)ll為數(shù)據(jù)項,用于使彌 散的最大主方向場和原始的最大主方向場盡可能相似。為平
(x.y)eA'
滑項,用于使彌散的最大主方向場平滑變化??梢钥闯?,在腦溝底部和腦回冠部,數(shù)據(jù)項的 權(quán)重較大,迫使彌散的最大主方向場與原始最大主方向場相似;在平坦的皮層區(qū)域,平滑項 的權(quán)重較大,迫使彌散的最大主方向場平滑的變化。為了將最大主方問場彌散轉(zhuǎn)化為一個 離散標(biāo)記問題,將 V(x)表示為 v(x) =
其中
表示 0 中離散角度的數(shù)目, 設(shè)為12。令小=j* Ji/n,, j G
,、表示小中離散角度的數(shù)目,設(shè)為9。令L = {l1,I2,. . .,ln}表示離散的彌散的最大主方向場的解空間 = {^,^,...,^^,其中??!
。以上的能量函數(shù)可以轉(zhuǎn)化為
alpha-expansion圖割方法有效的求解該能量函數(shù)。步驟4將彌散的最大主方向場投影到切平面內(nèi)。圖1(c)顯示了圖1(b)中對應(yīng)的大腦皮層表面上的彌散的最大主方向場。為了定量測試大腦皮層表面上最大主方向場彌散方法的效果,將該方法應(yīng)用到一 個仿真的表面上。在仿真的表面上可以精確解析計算每個頂點的理想最大主方向,然后對 每個頂點的位置添加隨機噪聲,并利用有限差分方法計算添加噪聲后每個頂點的最大主方 向。圖2(a)、(b)和(c)分別顯示了仿真表面上的理想最大主方向場、噪聲的最大主方向場 和彌散的最大主方向。通過以下公式計算彌散后的最大主方向場、添加噪聲的最大主方向 場與理想最大主方向場的平均角度差 其中,p(x)表示理想最大主方向,v(x)表示彌散的最大主方向或添加噪聲的最大 主方向。添加噪聲的最大主方向場與理想最大主方向場的平均角度差為67. 7度,而彌散后 的最大主方向場與理想最大主方向場的平均角度差為6. 3度,表明該方法達到了很好的彌 散結(jié)果。將該方法應(yīng)用到12個正常人的大腦皮層表面上。定義一致性來評價彌散結(jié)果 其中,v(x)表示頂點x的方向,(x,y) G N表示相鄰頂點集合。大的一致性值表 明方向場更加平滑和一致。因此,該值可以在某種程度上反映彌散結(jié)果。圖3(a)和(b)分 別顯示一例大腦皮層表面上的最大主方向場和彌散的最大主方向場的一致性。在理想情況 下,在彌散的最大主方向場中,一致性在平坦的皮層區(qū)域應(yīng)該接近于1.0,一致性在腦回冠 部和腦溝底部接近于0??梢悦黠@的看出,和原始的最大主方向場相比,彌散的最大主方向 場中,一致性在平坦的皮層區(qū)域增加了,而在腦溝底部和腦回冠部區(qū)域很好的保持。同時計 算了大腦皮層表面上一致性的分布直方圖,如圖4所示,可以看出相比于原始的最大主方 向場,彌散的最大主方向場中一致性明顯增大。圖5顯示了在12個大腦皮層表面上,平均 一致性在最大主方向場彌散前和彌散后的變化,可以看出平均一致性明顯的增大,表明了 本方法的有效性。
權(quán)利要求
一種三維大腦磁共振圖像的大腦皮層表面上最大主方向場彌散方法,其特征在于步驟如下步驟1對三維大腦磁共振圖像進行預(yù)處理和大腦皮層表面重建利用變形模型方法去除腦殼,利用配準(zhǔn)方法去除非大腦組織,利用高斯混合模型和馬爾科夫隨機場方法對圖像進行組織分割,得到白質(zhì)、灰質(zhì)和腦脊髓液三種組織類型表示的圖像,利用基于圖的方法對白質(zhì)圖像進行拓撲校正,利用Marching Cubes方法重建由三角形表示的大腦皮層表面;步驟2計算大腦皮層表面上頂點的主曲率、主方向以及主曲率的導(dǎo)數(shù)利用有限差分方法計算大腦皮層表面上頂點的主曲率、主方向以及主曲率的導(dǎo)數(shù),當(dāng)每個頂點的最大主曲率在最大主方向的導(dǎo)數(shù)為正值時,將該頂點的最大主方向改為相反方向;步驟3對最大主方向場進行彌散利用alpha-expansion圖割方法最小化能量函數(shù) <mrow><mi>E</mi><mo>=</mo><munder> <mi>&Sigma;</mi> <mrow><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow><mo>&Element;</mo><mi>N</mi> </mrow></munder><msub> <mi>V</mi> <mrow><mi>x</mi><mo>,</mo><mi>y</mi> </mrow></msub><mrow> <mo>(</mo> <msub><mi>l</mi><mi>x</mi> </msub> <mo>,</mo> <msub><mi>l</mi><mi>y</mi> </msub> <mo>)</mo></mrow><mo>+</mo><munder> <mi>&Sigma;</mi> <mrow><mi>x</mi><mo>&Element;</mo><mi>S</mi> </mrow></munder><msub> <mi>D</mi> <mi>x</mi></msub><mrow> <mo>(</mo> <msub><mi>l</mi><mi>x</mi> </msub> <mo>)</mo></mrow> </mrow>其中,x和y為大腦皮層表面上的頂點坐標(biāo),p(x)為頂點x處原始的最大主方向;(x,y)∈N表示x和y為相鄰頂點的集合,S為大腦皮層表面上的所有頂點的集合;||·||為2范數(shù);為頂點x彌散的最大主方向,lx為頂點x彌散的最大主方向在解空間的標(biāo)識且lx∈L;為頂點y彌散的最大主方向,ly為頂點y彌散的最大主方向在解空間的標(biāo)識且ly∈L;所述的L={l1,l2,...,ln}為離散的彌散的最大主方向場的解空間Θ={v1,v2,...,vn}的標(biāo)識,其中n=nθ×(nφ-2)+2;nθ為x-y平面中離散角度的數(shù)目為12~36;nφ為z軸方向中離散角度的數(shù)目為9~18;所述解空間Θ中v1=(0,0,1),vn=(0,0,-1);vi=(cos(((i-2)%nθ)·2π/nθ)sin((1+(i-2)/nθ)·π/nφ),sin(((i-2)%nθ)·2π/nθ)sin((1+(i-2)/nθ)·π/nφ),cos((1+(i-2)/nθ)·π/nφ)),1<i<n;所述g(x)=exp(-λ·|c(x)|),所述g(y)=exp(-λ·|c(y)|),所述h(x)=1-g(x),其中λ為權(quán)重參數(shù)為5.0~10.0,c(·)為最大主曲率;利用alpha-expansion圖割方法求解,得到彌散的最大主方向場;步驟4將彌散的最大主方向場投影到每個頂點的切平面內(nèi)。FSA00000156047900012.tif,FSA00000156047900013.tif,FSA00000156047900014.tif,FSA00000156047900015.tif
全文摘要
本發(fā)明涉及一種三維大腦磁共振圖像大腦皮層表面最大主方向場彌散方法,技術(shù)特征在于對三維大腦磁共振圖像進行預(yù)處理和大腦皮層表面重建,計算大腦皮層表面上頂點的主曲率、主方向以及主曲率的導(dǎo)數(shù);利用alpha-expanon圖割方法,最小化一個能量函數(shù)對最大主方向場進行彌散;將彌散的最大主方向場投影到切平面內(nèi)。本發(fā)明的優(yōu)點1、該方法在大腦皮層表面上不同區(qū)域的平滑項和數(shù)據(jù)項的權(quán)重設(shè)為不同,這樣可以更好的保持最大主方向場中固有的幾何結(jié)構(gòu)和不連續(xù)性;2、將最大主方向場彌散看做一個能量最小化問題,并轉(zhuǎn)換為一個離散標(biāo)記問題,利用alpha-expansion圖割方法可以有效求解該能量函數(shù)的強局部最優(yōu)。
文檔編號A61B5/055GK101866485SQ201010197300
公開日2010年10月20日 申請日期2010年6月10日 優(yōu)先權(quán)日2010年6月10日
發(fā)明者李剛, 聶晶鑫, 胡新韜, 趙天云, 郭雷, 韓軍偉 申請人:西北工業(yè)大學(xué)
網(wǎng)友詢問留言 已有0條留言
  • 還沒有人留言評論。精彩留言會獲得點贊!
1
久久精品久久久久久久性| 亚洲一区高清亚洲精品| 男女国产视频网站| 久久韩国三级中文字幕| 水蜜桃什么品种好| 国产91av在线免费观看| 丝袜喷水一区| 亚洲国产日韩欧美精品在线观看| 国产老妇女一区| 日韩av在线免费看完整版不卡| 国产精品福利在线免费观看| 五月伊人婷婷丁香| 99热网站在线观看| 中国国产av一级| 搡老妇女老女人老熟妇| 18禁裸乳无遮挡免费网站照片| 国产亚洲一区二区精品| 夜夜爽夜夜爽视频| 国产精品1区2区在线观看.| 日本免费在线观看一区| 国产不卡一卡二| 床上黄色一级片| 激情 狠狠 欧美| 久久午夜福利片| 三级毛片av免费| 亚洲人成网站在线播| 免费人成在线观看视频色| 伊人久久精品亚洲午夜| 精品一区二区三区人妻视频| 欧美激情久久久久久爽电影| 亚洲精品,欧美精品| 精品国产三级普通话版| 九草在线视频观看| 97超碰精品成人国产| 夜夜看夜夜爽夜夜摸| 欧美最新免费一区二区三区| 日韩欧美精品v在线| 欧美日韩精品成人综合77777| 国产黄a三级三级三级人| 一区二区三区四区激情视频| 亚洲精品成人av观看孕妇| 九九爱精品视频在线观看| 久久这里有精品视频免费| 少妇高潮的动态图| 亚洲欧美日韩无卡精品| 中文精品一卡2卡3卡4更新| 国产黄色小视频在线观看| 高清欧美精品videossex| 一级av片app| 亚洲三级黄色毛片| 人人妻人人看人人澡| 老司机影院毛片| 亚洲欧美日韩无卡精品| 听说在线观看完整版免费高清| 日韩欧美精品免费久久| 亚洲精品国产av蜜桃| 欧美高清性xxxxhd video| 青春草视频在线免费观看| 久久久久国产网址| 人体艺术视频欧美日本| 国产精品av视频在线免费观看| 亚洲欧洲日产国产| 国产精品嫩草影院av在线观看| 国产精品蜜桃在线观看| 人妻一区二区av| 国产精品一及| 亚洲自偷自拍三级| 亚洲va在线va天堂va国产| 亚洲精品中文字幕在线视频 | 91av网一区二区| 黑人高潮一二区| 亚洲精品久久午夜乱码| 国产在视频线精品| 91精品国产九色| 国产伦精品一区二区三区四那| 美女国产视频在线观看| 国产欧美另类精品又又久久亚洲欧美| av在线天堂中文字幕| 久久国内精品自在自线图片| 成人亚洲精品av一区二区| 小蜜桃在线观看免费完整版高清| 中文天堂在线官网| 国产精品嫩草影院av在线观看| 毛片女人毛片| av黄色大香蕉| 麻豆成人av视频| 国产老妇女一区| 91av网一区二区| 九色成人免费人妻av| 搡老乐熟女国产| 久久人人爽人人片av| 青春草视频在线免费观看| 欧美最新免费一区二区三区| 日韩不卡一区二区三区视频在线| 超碰av人人做人人爽久久| av女优亚洲男人天堂| 狠狠精品人妻久久久久久综合| 99热这里只有精品一区| 久久久精品免费免费高清| 国产毛片a区久久久久| 2022亚洲国产成人精品| 丰满乱子伦码专区| av女优亚洲男人天堂| 中文在线观看免费www的网站| 色5月婷婷丁香| av国产免费在线观看| 亚洲欧美日韩卡通动漫| 亚洲在线观看片| 中文天堂在线官网| 一区二区三区四区激情视频| 尤物成人国产欧美一区二区三区| 麻豆久久精品国产亚洲av| 亚洲国产欧美人成| 女人被狂操c到高潮| 国产高潮美女av| 久久99蜜桃精品久久| 国产黄a三级三级三级人| 性色avwww在线观看| 欧美极品一区二区三区四区| 中文资源天堂在线| 国产极品天堂在线| 欧美另类一区| 成年人午夜在线观看视频 | 国产又色又爽无遮挡免| 中文天堂在线官网| 精品午夜福利在线看| 久久久久久久午夜电影| 国产精品久久久久久精品电影小说 | 国产成人精品一,二区| 亚洲成色77777| 日韩人妻高清精品专区| 国产男人的电影天堂91| 国产人妻一区二区三区在| 欧美一级a爱片免费观看看| 国产老妇女一区| 久久久成人免费电影| 国产av国产精品国产| 又爽又黄无遮挡网站| 91精品一卡2卡3卡4卡| 亚洲激情五月婷婷啪啪| 亚洲欧洲国产日韩| 插阴视频在线观看视频| 久久久久久久亚洲中文字幕| 免费观看在线日韩| 国产精品人妻久久久影院| 一级毛片我不卡| 精品久久久久久成人av| 中文字幕av在线有码专区| 欧美+日韩+精品| 一级a做视频免费观看| 嫩草影院入口| 亚洲久久久久久中文字幕| 亚洲最大成人中文| 亚洲真实伦在线观看| 噜噜噜噜噜久久久久久91| 亚洲人与动物交配视频| 日韩成人伦理影院| 神马国产精品三级电影在线观看| 亚洲va在线va天堂va国产| 老女人水多毛片| 欧美日韩综合久久久久久| 亚洲欧美成人综合另类久久久| av在线播放精品| 欧美成人精品欧美一级黄| 不卡视频在线观看欧美| 国产一区二区在线观看日韩| 一本久久精品| 老女人水多毛片| 九九在线视频观看精品| 菩萨蛮人人尽说江南好唐韦庄| 亚洲人成网站高清观看| 国产精品一区二区三区四区免费观看| 久久精品夜夜夜夜夜久久蜜豆| 高清av免费在线| 国产美女午夜福利| 亚洲精品一二三| 日韩一本色道免费dvd| 精品久久久久久久久av| 最近2019中文字幕mv第一页| 国产精品人妻久久久久久| 22中文网久久字幕| 免费大片18禁| 色视频www国产| 精品久久久久久久久久久久久| 人人妻人人澡人人爽人人夜夜 | av在线播放精品| 国产精品久久久久久久久免| 成年版毛片免费区| 国产欧美另类精品又又久久亚洲欧美| 简卡轻食公司| 免费观看精品视频网站| 成人毛片a级毛片在线播放| 不卡视频在线观看欧美| 91av网一区二区| 精品一区二区免费观看| 91精品伊人久久大香线蕉| 97精品久久久久久久久久精品| 久久久久性生活片| 国产精品人妻久久久久久| 婷婷色麻豆天堂久久| 国国产精品蜜臀av免费| 蜜桃亚洲精品一区二区三区| 最近中文字幕2019免费版| 国产午夜精品久久久久久一区二区三区| 一个人免费在线观看电影| 国产亚洲一区二区精品| av天堂中文字幕网| 精品午夜福利在线看| 久99久视频精品免费| 亚洲人成网站在线观看播放| 日韩av在线免费看完整版不卡| 亚洲国产欧美在线一区| 天堂俺去俺来也www色官网 | 成人亚洲精品一区在线观看 | 色视频www国产| 日韩,欧美,国产一区二区三区| 亚洲自拍偷在线| 人人妻人人澡欧美一区二区| 国产视频首页在线观看| 国产成人精品一,二区| 亚洲成人中文字幕在线播放| 国产精品久久久久久精品电影小说 | 久久午夜福利片| 人人妻人人看人人澡| 久久久久久伊人网av| 国产中年淑女户外野战色| 亚洲成色77777| or卡值多少钱| 联通29元200g的流量卡| 免费av毛片视频| 青春草国产在线视频| 最近中文字幕2019免费版| 一级毛片久久久久久久久女| 可以在线观看毛片的网站| 少妇丰满av| 三级毛片av免费| 精品国产一区二区三区久久久樱花 | 少妇的逼好多水| 久久久色成人| 亚洲欧美日韩东京热| 国产亚洲91精品色在线| 伊人久久国产一区二区| 精品久久久久久久久亚洲| 亚洲国产精品成人久久小说| 欧美区成人在线视频| 我要看日韩黄色一级片| 国产淫语在线视频| 亚洲一级一片aⅴ在线观看| 韩国高清视频一区二区三区| 久久热精品热| 久久久久网色| 日本免费在线观看一区| 国精品久久久久久国模美| 99久久精品一区二区三区| 久久久久精品久久久久真实原创| 韩国av在线不卡| 在线天堂最新版资源| 日日干狠狠操夜夜爽| 亚洲丝袜综合中文字幕| 日韩三级伦理在线观看| 高清在线视频一区二区三区| 少妇丰满av| 日韩一区二区三区影片| 中文天堂在线官网| 好男人在线观看高清免费视频| 亚洲最大成人中文| 少妇裸体淫交视频免费看高清| 中文字幕久久专区| 99久久精品热视频| 久久久久久久大尺度免费视频| 两个人的视频大全免费| 99久久精品国产国产毛片| 欧美不卡视频在线免费观看| 亚洲精品国产成人久久av| 国产精品麻豆人妻色哟哟久久 | 国产精品麻豆人妻色哟哟久久 | 啦啦啦韩国在线观看视频| 草草在线视频免费看| 久久人人爽人人片av| 大话2 男鬼变身卡| 春色校园在线视频观看| 精品不卡国产一区二区三区| 好男人视频免费观看在线| av在线亚洲专区| 国产国拍精品亚洲av在线观看| 丝瓜视频免费看黄片| 人人妻人人澡欧美一区二区| 少妇熟女欧美另类| 亚洲精品亚洲一区二区| 五月天丁香电影| 色网站视频免费| 久久6这里有精品| 久久亚洲国产成人精品v| 久久久久久久久久黄片| videos熟女内射| 欧美日韩精品成人综合77777| 人妻系列 视频| 国产男女超爽视频在线观看| 一二三四中文在线观看免费高清| 亚洲性久久影院| 久久精品综合一区二区三区| 亚洲人成网站高清观看| 亚洲精品乱久久久久久| 免费黄频网站在线观看国产| 久久久久久伊人网av| 狠狠精品人妻久久久久久综合| 国产黄频视频在线观看| 韩国高清视频一区二区三区| 国产三级在线视频| 99热这里只有是精品在线观看| 大话2 男鬼变身卡| 五月天丁香电影| 你懂的网址亚洲精品在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲性久久影院| 狂野欧美激情性xxxx在线观看| 天堂影院成人在线观看| 精品一区二区免费观看| 日日啪夜夜爽| 高清日韩中文字幕在线| 免费黄频网站在线观看国产| 亚洲一级一片aⅴ在线观看| 99热这里只有精品一区| 亚洲av成人精品一二三区| 欧美成人精品欧美一级黄| 亚洲乱码一区二区免费版| 狂野欧美白嫩少妇大欣赏| 国产女主播在线喷水免费视频网站 | 特级一级黄色大片| 国产精品精品国产色婷婷| 赤兔流量卡办理| 午夜激情久久久久久久| 成人美女网站在线观看视频| 国产精品人妻久久久影院| 国产大屁股一区二区在线视频| 亚洲国产av新网站| 色尼玛亚洲综合影院| 久久久久精品久久久久真实原创| 九九爱精品视频在线观看| 国产一级毛片七仙女欲春2| 久久午夜福利片| 看十八女毛片水多多多| 午夜福利视频精品| 夫妻性生交免费视频一级片| 国模一区二区三区四区视频| 老司机影院毛片| 六月丁香七月| 插逼视频在线观看| 久久久久九九精品影院| 免费观看在线日韩| 精品一区二区三区视频在线| 纵有疾风起免费观看全集完整版 | 成年版毛片免费区| 嫩草影院精品99| 狠狠精品人妻久久久久久综合| 美女脱内裤让男人舔精品视频| 亚洲熟女精品中文字幕| 日本av手机在线免费观看| 超碰97精品在线观看| 久久久久久久久久人人人人人人| 丰满人妻一区二区三区视频av| 韩国av在线不卡| 国产精品.久久久| 日韩大片免费观看网站| 国产美女午夜福利| 欧美激情久久久久久爽电影| www.色视频.com| 精品欧美国产一区二区三| 国产精品不卡视频一区二区| av在线观看视频网站免费| 狂野欧美激情性xxxx在线观看| 久久精品国产亚洲av天美| 成人毛片60女人毛片免费| 国产成人a∨麻豆精品| 少妇的逼水好多| 亚洲熟女精品中文字幕| 免费观看性生交大片5| 精品亚洲乱码少妇综合久久| 2021天堂中文幕一二区在线观| 精品久久久精品久久久| 日韩欧美一区视频在线观看 | or卡值多少钱| 一级爰片在线观看| 欧美bdsm另类| 熟女人妻精品中文字幕| 久久精品国产亚洲av涩爱| 日韩三级伦理在线观看| 久久草成人影院| 国产熟女欧美一区二区| 美女xxoo啪啪120秒动态图| 乱码一卡2卡4卡精品| 亚洲国产精品sss在线观看| 一本一本综合久久| 人人妻人人澡人人爽人人夜夜 | 水蜜桃什么品种好| 黄色配什么色好看| 亚洲熟妇中文字幕五十中出| 国产乱来视频区| 99视频精品全部免费 在线| 我要看日韩黄色一级片| 乱系列少妇在线播放| 一个人看视频在线观看www免费| 国产伦精品一区二区三区视频9| 亚洲婷婷狠狠爱综合网| 中文字幕av在线有码专区| 国产 亚洲一区二区三区 | 亚洲欧美中文字幕日韩二区| 成年版毛片免费区| 18禁在线播放成人免费| 欧美97在线视频| 亚洲欧洲国产日韩| 国产视频内射| 欧美xxⅹ黑人| 亚洲aⅴ乱码一区二区在线播放| 午夜免费男女啪啪视频观看| 亚州av有码| 亚洲精品乱码久久久v下载方式| 久久99精品国语久久久| 精品欧美国产一区二区三| 免费av不卡在线播放| 久久精品久久精品一区二区三区| 18+在线观看网站| 九九久久精品国产亚洲av麻豆| 久久人人爽人人片av| 成人亚洲欧美一区二区av| or卡值多少钱| 成人综合一区亚洲| 亚洲怡红院男人天堂| 国产69精品久久久久777片| 麻豆久久精品国产亚洲av| 久久午夜福利片| 日本熟妇午夜| 久久精品综合一区二区三区| 国产色爽女视频免费观看| 日韩不卡一区二区三区视频在线| 国产69精品久久久久777片| 精品人妻偷拍中文字幕| 精品久久久精品久久久| 色综合色国产| 久久综合国产亚洲精品| 80岁老熟妇乱子伦牲交| 777米奇影视久久| 精品久久久久久久久久久久久| 2022亚洲国产成人精品| 女的被弄到高潮叫床怎么办| 国产精品久久久久久精品电影| 国产日韩欧美在线精品| 日韩强制内射视频| 国内精品一区二区在线观看| 最近手机中文字幕大全| 午夜福利网站1000一区二区三区| 国产乱来视频区| 欧美bdsm另类| 精品久久久噜噜| 丰满乱子伦码专区| 国产精品一区二区三区四区免费观看| 午夜福利视频1000在线观看| 搞女人的毛片| 一区二区三区四区激情视频| 黄色日韩在线| 麻豆成人午夜福利视频| 国产一区二区三区av在线| 成年av动漫网址| 免费电影在线观看免费观看| 少妇的逼水好多| 日日干狠狠操夜夜爽| 美女脱内裤让男人舔精品视频| 亚洲欧美中文字幕日韩二区| 国产真实伦视频高清在线观看| 亚洲久久久久久中文字幕| 黑人高潮一二区| kizo精华| 亚洲婷婷狠狠爱综合网| 男的添女的下面高潮视频| 亚洲精品第二区| 一区二区三区高清视频在线| 男插女下体视频免费在线播放| 亚洲经典国产精华液单| 六月丁香七月| 成人毛片a级毛片在线播放| 成人鲁丝片一二三区免费| 热99在线观看视频| 国产探花在线观看一区二区| 91久久精品国产一区二区三区| 99久久精品一区二区三区| 男女边吃奶边做爰视频| 午夜久久久久精精品| 日本wwww免费看| 丝瓜视频免费看黄片| 亚洲怡红院男人天堂| 免费观看精品视频网站| 一个人免费在线观看电影| 欧美区成人在线视频| 欧美xxxx性猛交bbbb| av国产免费在线观看| 十八禁国产超污无遮挡网站| av在线天堂中文字幕| 波野结衣二区三区在线| 日韩 亚洲 欧美在线| 国产一区二区在线观看日韩| 成人综合一区亚洲| 日本爱情动作片www.在线观看| 全区人妻精品视频| 99热这里只有是精品在线观看| 久久99热这里只有精品18| 精品一区二区三卡| 一夜夜www| 亚洲电影在线观看av| 99热全是精品| 日本wwww免费看| 99久久精品一区二区三区| 两个人视频免费观看高清| 街头女战士在线观看网站| 国产91av在线免费观看| 一级a做视频免费观看| 人妻制服诱惑在线中文字幕| 国产老妇伦熟女老妇高清| 久久99蜜桃精品久久| 中文在线观看免费www的网站| 综合色av麻豆| 大陆偷拍与自拍| 一级黄片播放器| 亚洲国产av新网站| 国产精品人妻久久久久久| 亚洲婷婷狠狠爱综合网| 日韩精品青青久久久久久| 亚洲图色成人| 日韩强制内射视频| 精品不卡国产一区二区三区| 免费黄色在线免费观看| 国产精品久久久久久精品电影| 看黄色毛片网站| 18禁在线无遮挡免费观看视频| 日韩一本色道免费dvd| 高清毛片免费看| 精品久久国产蜜桃| 欧美不卡视频在线免费观看| 亚洲成人久久爱视频| 国产一区二区在线观看日韩| 亚洲丝袜综合中文字幕| 一个人观看的视频www高清免费观看| 成人毛片60女人毛片免费| 女的被弄到高潮叫床怎么办| 国产成人福利小说| 三级经典国产精品| 一级黄片播放器| 国产精品一及| 久久久久久久久久成人| 久久久色成人| 日韩伦理黄色片| 成人午夜精彩视频在线观看| 久久午夜福利片| 国产av在哪里看| av在线老鸭窝| 少妇裸体淫交视频免费看高清| 国产午夜福利久久久久久| 免费看日本二区| 国产国拍精品亚洲av在线观看| 国产精品一区二区三区四区免费观看| 日韩 亚洲 欧美在线| 亚洲伊人久久精品综合| 国产精品人妻久久久久久| 国语对白做爰xxxⅹ性视频网站| 久久99蜜桃精品久久| 成人美女网站在线观看视频| 国产av国产精品国产| 国内少妇人妻偷人精品xxx网站| 国产91av在线免费观看| 国产精品爽爽va在线观看网站| 久久久久久久久久人人人人人人| 久久99热6这里只有精品| 亚洲精品乱久久久久久| 久久久久国产网址| 亚洲成人久久爱视频| 色综合亚洲欧美另类图片| 久久久久免费精品人妻一区二区| 禁无遮挡网站| 亚洲成人中文字幕在线播放| 国产欧美日韩精品一区二区| 亚洲精品影视一区二区三区av| 亚洲精品乱码久久久久久按摩| 精品久久久久久久久久久久久| 免费大片18禁| 欧美精品一区二区大全| 亚洲国产精品国产精品| 欧美性感艳星| 纵有疾风起免费观看全集完整版 | 亚洲精品自拍成人| 最近最新中文字幕大全电影3| 久久99热这里只频精品6学生| av网站免费在线观看视频 | 1000部很黄的大片| 国产视频内射| 九九久久精品国产亚洲av麻豆| 成人综合一区亚洲| 搡老乐熟女国产| 一个人看视频在线观看www免费| 美女内射精品一级片tv| 身体一侧抽搐| 久久精品久久久久久久性| 国产伦精品一区二区三区视频9| 国产免费又黄又爽又色| 最近视频中文字幕2019在线8| 不卡视频在线观看欧美| 丰满少妇做爰视频| 日韩人妻高清精品专区| 欧美成人精品欧美一级黄| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 51国产日韩欧美| 久久精品国产自在天天线| 十八禁国产超污无遮挡网站| a级毛片免费高清观看在线播放| 天天一区二区日本电影三级| 能在线免费看毛片的网站| 中文字幕av成人在线电影|