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

    基于變慣性系數(shù)active demons算法的DTI多通道配準(zhǔn)研究

    2021-08-19 06:13:22鄭偉周楊李文華劉帥奇張曉丹馬澤鵬
    關(guān)鍵詞:驅(qū)動(dòng)力浮動(dòng)慣性

    鄭偉,周楊,李文華,劉帥奇,張曉丹,馬澤鵬

    (1.河北大學(xué) 電子信息工程學(xué)院,河北 保定 071002;2.河北省數(shù)字醫(yī)療工程重點(diǎn)實(shí)驗(yàn)室,河北 保定 071002;3.河北省機(jī)器視覺(jué)中心,河北 保定 071002;4.河北大學(xué)附屬醫(yī)院 CT-MRI診斷科,河北 保定 071000)

    彌散張量成像(diffision tensor imaging,DTI)是20世紀(jì)90年代新興的一種大腦結(jié)構(gòu)成像方法,它利用分子在大腦中的運(yùn)動(dòng)情況來(lái)繪成圖像,可以有效揭示阿爾茨海默病、精神分裂[1-2]、閱讀障礙等引起的大腦組織結(jié)構(gòu)反常變化,引導(dǎo)醫(yī)療人員進(jìn)行大腦手術(shù).在臨床上,對(duì)同一個(gè)病人不同時(shí)期DTI圖像研究分析和診斷具有十分重要的臨床價(jià)值.在做這些類型的比較之前,必須要進(jìn)行DTI圖像配準(zhǔn).DTI圖像配準(zhǔn)是指將2組或者多組DTI圖像進(jìn)行空間幾何變換,使各個(gè)像素或體素在相同解剖結(jié)構(gòu)上能夠相對(duì)應(yīng),而且DTI圖像配準(zhǔn)的結(jié)果將直接影響纖維束追蹤的效果.基于光流場(chǎng)理論[3]的demons算法由于其算法的高效性和理論基礎(chǔ)的完整性而被廣泛應(yīng)用于非剛性圖像配準(zhǔn)中[4],原始的demons算法由Thifion等在1998年提出[4],其原理是利用參考圖像的梯度信息以及浮動(dòng)圖像與參考圖像的差作為形變向量驅(qū)動(dòng)浮動(dòng)圖像朝著參考圖像的方向發(fā)生形變.原始的demons算法收斂速度慢,對(duì)于大形變的配準(zhǔn)取得的效果也不理想.針對(duì)以上問(wèn)題,一些學(xué)者往驅(qū)動(dòng)力公式中加入新的信息項(xiàng)或參數(shù),例如幾何形狀約束項(xiàng)[5],李果霖等[6]將系數(shù)f(θ)引入驅(qū)動(dòng)力公式解決了參考圖像和浮動(dòng)圖像梯度差異過(guò)大引起的過(guò)配準(zhǔn)問(wèn)題.蔣晨嬌等[7]為了解決傳統(tǒng)DTI圖像配準(zhǔn)方法容易忽略空間信息及高維特性的問(wèn)題,提出基于改進(jìn)雞群優(yōu)化算法的AD患者DTI圖像配準(zhǔn),提高了配準(zhǔn)精度,取得較好的配準(zhǔn)結(jié)果.為了保持配準(zhǔn)前后解剖結(jié)構(gòu)的一致性,Guo等[8]提出DTI配準(zhǔn)的高維空間標(biāo)準(zhǔn)化算法,先對(duì)待配準(zhǔn)圖像進(jìn)行空間標(biāo)準(zhǔn)化,然后對(duì)DTI圖像進(jìn)行高維標(biāo)準(zhǔn)化.通過(guò)對(duì)比實(shí)驗(yàn)驗(yàn)證了該算法的有效性.張桂梅等[9]提出自適應(yīng)分?jǐn)?shù)階demons算法解決了圖像容易陷入局部最優(yōu)的問(wèn)題.許磊等[10]針對(duì)原始demons算法不適合配準(zhǔn)大形變圖像的缺點(diǎn)對(duì)demons驅(qū)動(dòng)力公式進(jìn)行了重新定義.針對(duì)多模態(tài)醫(yī)學(xué)圖像配準(zhǔn)效果差精度低的問(wèn)題,Simona等[11]提出結(jié)合稀疏和稠密特征改進(jìn)腦DTI圖像多模態(tài)配準(zhǔn)方法,使用直方圖梯度的互信息作為配準(zhǔn)標(biāo)準(zhǔn),提高了多模態(tài)圖像配準(zhǔn)的精度.

    Wang等[12]提出active demons算法,該算法對(duì)于形變較大的圖像有良好的配準(zhǔn)效果,收斂速度和精度相對(duì)demons算法有很大提高,但圖像的拓?fù)浣Y(jié)構(gòu)容易發(fā)生改變.對(duì)于均化系數(shù)α的取值問(wèn)題,α不能兼顧圖像的大小形變,而在配準(zhǔn)過(guò)程中大小形變是同時(shí)存在的,針對(duì)這個(gè)問(wèn)題,薛鵬等[13]提出基于平衡系數(shù)的active demons配準(zhǔn)算法,引入新的參數(shù)平衡系數(shù)k,固定均化系數(shù),手動(dòng)選取最優(yōu)值,得到最優(yōu)平衡系數(shù)k.王昌等[14]提出基于局部聯(lián)合熵梯度的active demons算法,在驅(qū)動(dòng)力公式中加入新的信息項(xiàng)局部聯(lián)合熵梯度,并且引入了多分辨率的策略.薛文靜等[15]也將局部熵引入active demons算法,將改進(jìn)了active demons算法的驅(qū)動(dòng)力公式應(yīng)用于多模態(tài)的醫(yī)學(xué)圖像配準(zhǔn).徐曉瑩等[16]將平衡系數(shù)k當(dāng)做一個(gè)變量引入多通道DTI圖像配準(zhǔn),提出變參數(shù)active demons算法,令k取不同值,可以動(dòng)態(tài)改變形變向量u的大小,其中手動(dòng)調(diào)整均化系數(shù)α,系數(shù)的最優(yōu)值通過(guò)實(shí)驗(yàn)測(cè)得,但是收斂速度依然不是很高.

    本文在此基礎(chǔ)上提出了基于變慣性系數(shù)active demons算法的DTI多通道配準(zhǔn),將慣性系數(shù)b引入變參數(shù)active demons算法中,并且改變b的取值結(jié)合平衡系數(shù)一起循環(huán)迭代,實(shí)驗(yàn)結(jié)果顯示與變參數(shù)active demons算法相比進(jìn)一步提高了算法的收斂速度并且得到了最好的配準(zhǔn)精度,同時(shí)圖像的拓?fù)浣Y(jié)構(gòu)得到了很好的保持.

    1 變慣性系數(shù)active demons算法及在DTI多通道配準(zhǔn)中的應(yīng)用

    水分子在純凈液體中向各個(gè)方向彌散的程度相同,把這種性質(zhì)稱作各向同性.但是大腦中有各種各樣的組織會(huì)阻礙水分子的運(yùn)動(dòng),水分子向各方向的彌散存在差異,這種差異稱作各向異性.由于各向異性的存在,水分子在大腦中的移動(dòng)可以看作橢球形,這種性質(zhì)可以用數(shù)學(xué)上的對(duì)稱二階張量來(lái)表示,如式(1)所示:

    (1)

    對(duì)稱矩陣有9個(gè)元素,獨(dú)立的元素只有6個(gè),因此至少需要在6個(gè)不同的方向上測(cè)量彌散張量.

    連續(xù)運(yùn)動(dòng)的2幀圖像存在亮度的變化即光流現(xiàn)象,Thirion等[17]將光流的概念引入到圖像配準(zhǔn)中,把參考圖像和浮動(dòng)圖像看作連續(xù)運(yùn)動(dòng)的圖像,由于光流現(xiàn)象的存在,這2幅圖像存在灰度梯度的差異,利用這種差異構(gòu)造一個(gè)驅(qū)動(dòng)力來(lái)驅(qū)動(dòng)浮動(dòng)圖像向參考圖像發(fā)生形變來(lái)進(jìn)行圖像的配準(zhǔn).假定2幅圖像的亮度不變,則二維圖像的光流場(chǎng)方程如式(2)所示:

    (2)

    式中,u表示點(diǎn)(x,y)處的形變向量;M(x,y)和S(x,y)分別表示浮動(dòng)圖像和參考圖像的灰度值;S(x,y)表示參考圖像在(x,y)處的梯度.為防止S(x,y)趨于0或不存在時(shí)方程失去意義,將式(2)改為式(3):

    (3)

    式中,α為均化系數(shù).

    由demons算法的原理可知連續(xù)的2幀圖像形變向量較小,所以該算法不適用于大形變圖像配準(zhǔn),且收斂速度慢.

    針對(duì)以上demons算法出現(xiàn)的問(wèn)題,Rogelj等[18]提出active demons算法把浮動(dòng)圖像的梯度信息也引入形變方程見(jiàn)式4,參考圖像的梯度信息和浮動(dòng)圖像梯度信息的合力作為驅(qū)動(dòng)浮動(dòng)圖像發(fā)生形變的demons力,同時(shí)改變均化系數(shù)α的取值也可以改變驅(qū)動(dòng)力的大小,該算法有效解決了單一依靠參考圖像梯度信息的缺陷,收斂速度和精度得到提高.

    (4)

    均化系數(shù)α越大,形變向量u越小,收斂速度變慢,配準(zhǔn)時(shí)間較長(zhǎng),但精度更高,適合小形變區(qū)域的圖像配準(zhǔn);α越小,形變向量u越大,收斂速度更快,配準(zhǔn)時(shí)間較短,但精度低,適合大形變區(qū)域的圖像配準(zhǔn).

    對(duì)均化系數(shù)α分析可知,α不能同時(shí)兼顧大形變和小形變的圖像配準(zhǔn).薛鵬等[13]受參數(shù)α的啟發(fā),在分母的梯度模值平方項(xiàng)上加入一個(gè)新的參數(shù)平衡系數(shù)k,平衡系數(shù)k和均化系數(shù)α聯(lián)合調(diào)整驅(qū)動(dòng)力強(qiáng)度,如式(5)所示:

    (5)

    如果固定α,k值減小,u增大,形變向量增大,收斂速度加快;k值增大,u減小,形變向量減小,收斂速度變慢.

    將6個(gè)獨(dú)立的分量Dxx、Dyy、Dzz、Dxy、Dxz、Dyz分別作為算法的輸入疊加求和再取平均值作為驅(qū)動(dòng)參考圖像發(fā)生形變的demons力.對(duì)于DTI圖像中的任意體素其位移的迭代如式(6)所示:

    (6)

    針對(duì)active demons算法收斂速度慢的問(wèn)題,在un前加入慣性系數(shù)b,根據(jù)物理學(xué)中的慣性概念[19],慣性系數(shù)b為上次迭代對(duì)本次迭代的影響.當(dāng)上次結(jié)果對(duì)本次迭代影響較小時(shí),令b取較小的值,反之則增大b的取值.通過(guò)動(dòng)態(tài)改變慣性系數(shù)的取值來(lái)影響驅(qū)動(dòng)力強(qiáng)度是本文研究的核心.變慣性系數(shù)active demons算法的驅(qū)動(dòng)力迭代公式如式(7)所示.

    (7)

    利用式(7)計(jì)算浮動(dòng)圖像的形變向量,位移場(chǎng)的平滑使用高斯濾波器,插值采用三線性插值,利用主方向保留法(presreving principle diffusion,PPD)對(duì)張量的方向進(jìn)行重定向,具體流程如圖1所示:

    圖1 變慣性系數(shù)active demons算法配準(zhǔn)流程Fig.1 Active demons algorithm with variable inertia coefficient registration process

    2 實(shí)驗(yàn)及結(jié)果分析

    2.1 實(shí)驗(yàn)數(shù)據(jù)及預(yù)處理

    本文使用的DTI數(shù)據(jù)來(lái)自標(biāo)準(zhǔn)數(shù)據(jù)庫(kù)ADNI(https://ida.loni.usc.edu/login.jsp?Project=AD-NI),數(shù)據(jù)來(lái)自同一AD患者在2個(gè)不同時(shí)期采用相同設(shè)備和參數(shù)采集的DTI數(shù)據(jù),設(shè)備為3T的GE醫(yī)療系統(tǒng),重復(fù)時(shí)間(TR)為9 050 ms,回波時(shí)間(TE)為61.8 ms,翻轉(zhuǎn)角(FA)為90°,圖像矩陣大小為256×256×60,厚度為2.7 mm,反映附加梯度場(chǎng)快慢的參數(shù)b為1 000 s/mm2.

    采用FSL(http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Fsllnstallation)對(duì)DTI原始實(shí)驗(yàn)數(shù)據(jù)做預(yù)處理,原始數(shù)據(jù)格式轉(zhuǎn)換后先對(duì)DTI數(shù)據(jù)做渦流矯正輸出.nii.gz格式的DTI數(shù)據(jù),再去除非腦組織得到全腦mask文件,最后計(jì)算FA、MD等全張量數(shù)據(jù).原始的軸向圖見(jiàn)圖2.

    2.2 變慣性系數(shù) active demons算法下的多通道DTI配準(zhǔn)結(jié)果與分析

    引入慣性系數(shù)b,單一的b值在取值大于1的情況下可以暫時(shí)增加算法的收斂速度,但是隨著迭代次數(shù)的增加,配準(zhǔn)精度降低,不能得到好的實(shí)驗(yàn)結(jié)果,本文的改進(jìn)是動(dòng)態(tài)改變慣性系數(shù)的取值,使配準(zhǔn)前期的b取一個(gè)較大的值隨著迭代的增加b逐漸減小,配合平衡系數(shù)k一起調(diào)節(jié)驅(qū)動(dòng)力的大小.下面先討論單一的b值對(duì)實(shí)驗(yàn)結(jié)果影響.令均化系數(shù)α=0.007,k[16]的取值范圍為0.5

    a.軸向參考FA圖;b.軸向浮動(dòng)FA圖.圖2 原始軸向圖像Fig.2 Original axial image

    圖3 單一b值的MSE圖Fig.3 MSE diagram for a single b value

    縱坐標(biāo)MSE(mean square error,均方誤差)為均方誤差,橫坐標(biāo)為迭代次數(shù).當(dāng)b的取值小于1時(shí),圖像的收斂速度和精度明顯降低,這是因?yàn)榕錅?zhǔn)前期k的取值較小,這樣可以配準(zhǔn)大形變量,b的取值如果也較小會(huì)減小前期的形變量.因此b的取值小于1會(huì)降低收斂速度和精度.當(dāng)b等于1時(shí)就等價(jià)于變參數(shù)active demons 的迭代方程.當(dāng)b大于1時(shí),慣性系數(shù)b與平衡系數(shù)k一同加快收斂速度,隨著迭代的進(jìn)行,大形變得到配準(zhǔn),可是對(duì)于小形變出現(xiàn)過(guò)配準(zhǔn),圖像拓?fù)浣Y(jié)構(gòu)發(fā)生改變,主要原因是b值偏大,不能兼顧小形變.從圖3可以看出無(wú)論b值大于1還是小于1都不能達(dá)到理想的效果甚至出現(xiàn)過(guò)配準(zhǔn)現(xiàn)象.

    本文提出的基于變慣性系數(shù)的active demons算法可以解決以上出現(xiàn)的問(wèn)題,前期通過(guò)讓b取較大的值加快收斂速度,隨著迭代的進(jìn)行,大形變區(qū)域得到配準(zhǔn),再降低b值,配準(zhǔn)小形變區(qū)域.本實(shí)驗(yàn)中根據(jù)文獻(xiàn)[16],α=0.007,高斯平滑的標(biāo)準(zhǔn)差σ采用Park等[20]在demons算法下設(shè)定的11、9、5、3.標(biāo)準(zhǔn)差σ隨迭代次數(shù)的增加在迭代次數(shù)每1/4處減小,平衡系數(shù)k隨著σ按0.5、0.7、0.9、1的順序依次增大,b隨著標(biāo)準(zhǔn)差按1.15、1.05、1.01、1的順序依次減小.實(shí)驗(yàn)結(jié)果如圖4所示.

    圖4 3種算法的對(duì)比MSE曲線Fig.4 Comparison of MSE curves of three algorithms

    AD(active demons),IAD(improve active demons),IIAD(improve inertia coefficient active demons)分別為active demons算法、變參數(shù)active demons算法、變慣性系數(shù)active demons算法的MSE曲線,由圖4可知,本文算法迭代27次的結(jié)果與變參數(shù)active demons算法迭代50次的結(jié)果相同,本文算法的MSE值比變參數(shù)active demons算法提高13%,比active demons算法的MSE值提高19%,變慣性系數(shù)active demons算法相較于active demons算法和變參數(shù)active demons算法取得了較好的實(shí)驗(yàn)結(jié)果,見(jiàn)表1.

    表1 3種算法在相同迭代次數(shù)下的MSE值比較Tab.1 Comparison of MSE values of the three algorithms with the same number of iterations

    配準(zhǔn)后的FA圖及FA差值圖如圖5所示.

    a.原始demons算法;b.active demons算法;c.變參數(shù)active demons算法;d.變慣性系數(shù)active demons算法.圖5 4種算法的FA圖和FA差值圖Fig.5 FA images and FA difference graphs of the four algorithms

    可以看出變慣性系數(shù)active demons算法的配準(zhǔn)效果要遠(yuǎn)遠(yuǎn)好于demons算法,并且相對(duì)于其他2種算法也取得了最小的一個(gè)MSE值.

    3 總結(jié)

    在變參數(shù)active demons算法基礎(chǔ)上做改進(jìn),引入慣性系數(shù)b,討論b的取值對(duì)實(shí)驗(yàn)結(jié)果的影響,得出結(jié)論單一b值無(wú)法取得理想的結(jié)果,b在一定范圍內(nèi)取合適的變值可以提高算法的收斂速度和精度.比較配準(zhǔn)后圖像FA圖,發(fā)現(xiàn)配準(zhǔn)后圖像的拓?fù)浣Y(jié)構(gòu)得到很好的保持.本文算法也存在如下不足,由于改進(jìn)算法引入了新參數(shù),在進(jìn)行配準(zhǔn)實(shí)驗(yàn)時(shí),為了達(dá)到高配準(zhǔn)精度,針對(duì)不同的圖像類型,可能需要進(jìn)行一定的參數(shù)取值實(shí)驗(yàn)工作,以選擇最優(yōu)的配準(zhǔn)結(jié)果.

    猜你喜歡
    驅(qū)動(dòng)力浮動(dòng)慣性
    中國(guó)船級(jí)社(CCS)發(fā)布 《海上浮動(dòng)設(shè)施入級(jí)規(guī)范》(2023)
    你真的了解慣性嗎
    沖破『慣性』 看慣性
    油價(jià)上漲的供需驅(qū)動(dòng)力能否持續(xù)
    一種用于剪板機(jī)送料的液壓浮動(dòng)夾鉗
    溫暖厚實(shí),驅(qū)動(dòng)力強(qiáng)勁 秦朝 QM2018/QC2350前后級(jí)功放
    突出文化產(chǎn)業(yè)核心驅(qū)動(dòng)力
    以創(chuàng)新為驅(qū)動(dòng)力,兼具學(xué)院派的嚴(yán)謹(jǐn)態(tài)度 Q Acoustics
    帶有浮動(dòng)機(jī)構(gòu)的曲軸孔鏜刀應(yīng)用研究
    無(wú)處不在的慣性
    青春草亚洲视频在线观看| 国产亚洲午夜精品一区二区久久 | 久久久久久久国产电影| 一夜夜www| 亚洲精品国产av蜜桃| 80岁老熟妇乱子伦牲交| 国产 一区精品| 欧美不卡视频在线免费观看| 国产淫片久久久久久久久| 国内揄拍国产精品人妻在线| 天堂av国产一区二区熟女人妻| 简卡轻食公司| 免费观看a级毛片全部| 精品一区二区三卡| 国产精品女同一区二区软件| 亚洲精品aⅴ在线观看| 日日啪夜夜爽| 国产精品久久久久久精品电影| a级一级毛片免费在线观看| 美女cb高潮喷水在线观看| 精华霜和精华液先用哪个| 亚洲在线自拍视频| 亚洲成人久久爱视频| 久久综合国产亚洲精品| 日韩av在线免费看完整版不卡| 18+在线观看网站| 一个人观看的视频www高清免费观看| 2022亚洲国产成人精品| 国产综合懂色| 欧美xxxx性猛交bbbb| 亚洲av电影在线观看一区二区三区 | 深爱激情五月婷婷| 日本wwww免费看| 亚洲av中文av极速乱| 亚洲精品一二三| 亚洲精品国产av蜜桃| 日日摸夜夜添夜夜爱| 我的女老师完整版在线观看| 日韩av不卡免费在线播放| 少妇熟女欧美另类| 日本av手机在线免费观看| av在线蜜桃| 尾随美女入室| 日韩大片免费观看网站| 欧美变态另类bdsm刘玥| 久久6这里有精品| 日韩,欧美,国产一区二区三区| 日日干狠狠操夜夜爽| 亚洲国产精品专区欧美| 免费看光身美女| 国产亚洲av嫩草精品影院| 国产精品三级大全| 五月伊人婷婷丁香| 黄色配什么色好看| 在线 av 中文字幕| 啦啦啦中文免费视频观看日本| 亚洲av福利一区| 2021天堂中文幕一二区在线观| 97超视频在线观看视频| 国产乱人偷精品视频| 日日撸夜夜添| 三级毛片av免费| 少妇人妻一区二区三区视频| 91精品一卡2卡3卡4卡| 麻豆乱淫一区二区| 97人妻精品一区二区三区麻豆| 99re6热这里在线精品视频| 日韩一区二区三区影片| 国产精品久久久久久精品电影| 波野结衣二区三区在线| 简卡轻食公司| 麻豆国产97在线/欧美| 久久精品国产鲁丝片午夜精品| 亚洲成人中文字幕在线播放| 1000部很黄的大片| 欧美不卡视频在线免费观看| 久久久久久久久中文| 国产大屁股一区二区在线视频| 亚洲成人一二三区av| 亚洲精品久久午夜乱码| 国产av不卡久久| 亚洲真实伦在线观看| a级一级毛片免费在线观看| 午夜免费激情av| 国产 亚洲一区二区三区 | 国产在视频线在精品| 99热这里只有是精品在线观看| 青春草国产在线视频| 国产黄色小视频在线观看| 精品人妻视频免费看| 免费看日本二区| 午夜免费激情av| 六月丁香七月| 最近视频中文字幕2019在线8| 亚洲av二区三区四区| 亚洲国产精品成人综合色| 青青草视频在线视频观看| 男女边摸边吃奶| av在线播放精品| 少妇丰满av| 亚洲精品乱码久久久久久按摩| 国产毛片a区久久久久| 能在线免费看毛片的网站| 日本-黄色视频高清免费观看| 久久这里有精品视频免费| 免费观看性生交大片5| 国产精品国产三级国产专区5o| 久久这里只有精品中国| 五月天丁香电影| 久久精品久久久久久噜噜老黄| 久久久久国产网址| 国产亚洲av片在线观看秒播厂 | 国产高潮美女av| 久久久久久久久久久丰满| 七月丁香在线播放| 精品人妻一区二区三区麻豆| 噜噜噜噜噜久久久久久91| 国产黄片美女视频| 少妇猛男粗大的猛烈进出视频 | 一级黄片播放器| 欧美激情国产日韩精品一区| 久久6这里有精品| 国产精品1区2区在线观看.| 国模一区二区三区四区视频| 国产白丝娇喘喷水9色精品| 久久人人爽人人爽人人片va| a级一级毛片免费在线观看| 亚洲18禁久久av| 麻豆久久精品国产亚洲av| 国产淫片久久久久久久久| 亚洲aⅴ乱码一区二区在线播放| 国产 亚洲一区二区三区 | av黄色大香蕉| av国产久精品久网站免费入址| 久久久久久久久久黄片| 日本午夜av视频| 国产精品1区2区在线观看.| 天堂av国产一区二区熟女人妻| 人妻少妇偷人精品九色| 能在线免费看毛片的网站| 卡戴珊不雅视频在线播放| 亚洲综合色惰| 夫妻性生交免费视频一级片| 日韩国内少妇激情av| 国产探花极品一区二区| 午夜免费观看性视频| 国产国拍精品亚洲av在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产在视频线精品| 国产黄a三级三级三级人| 亚洲精品视频女| 亚洲精品自拍成人| av天堂中文字幕网| 人人妻人人澡人人爽人人夜夜 | 男人狂女人下面高潮的视频| 婷婷色综合www| 欧美日韩国产mv在线观看视频 | 天天躁夜夜躁狠狠久久av| 精品人妻熟女av久视频| 亚洲成人精品中文字幕电影| 精品久久久久久久久久久久久| 69av精品久久久久久| 久久久久久久久大av| 免费无遮挡裸体视频| 久久久久免费精品人妻一区二区| 男人舔女人下体高潮全视频| 国产乱人视频| 欧美日韩亚洲高清精品| 国产91av在线免费观看| 亚洲成人久久爱视频| 91久久精品国产一区二区三区| 3wmmmm亚洲av在线观看| 亚洲成人中文字幕在线播放| 中文欧美无线码| 亚洲av中文av极速乱| 成人毛片60女人毛片免费| 久久久久久久久久人人人人人人| 国产精品麻豆人妻色哟哟久久 | 熟妇人妻久久中文字幕3abv| 国产精品久久久久久久久免| 免费黄色在线免费观看| 卡戴珊不雅视频在线播放| 国产精品福利在线免费观看| 欧美变态另类bdsm刘玥| 成人av在线播放网站| 日本wwww免费看| 大话2 男鬼变身卡| 91狼人影院| 精品一区二区三卡| 欧美成人精品欧美一级黄| av国产久精品久网站免费入址| 少妇裸体淫交视频免费看高清| 日韩三级伦理在线观看| 两个人的视频大全免费| 欧美成人午夜免费资源| 久久草成人影院| 又爽又黄a免费视频| 国产高清三级在线| 久久久久久久久久久丰满| 久久久久久久久久久丰满| 老女人水多毛片| 十八禁国产超污无遮挡网站| 亚洲国产成人一精品久久久| 最近的中文字幕免费完整| 18禁在线无遮挡免费观看视频| or卡值多少钱| 美女主播在线视频| 免费观看av网站的网址| 少妇高潮的动态图| 高清毛片免费看| 女人被狂操c到高潮| 久久草成人影院| 亚洲美女搞黄在线观看| 直男gayav资源| 少妇熟女欧美另类| 美女高潮的动态| 禁无遮挡网站| 男人爽女人下面视频在线观看| 亚洲av免费高清在线观看| 在线 av 中文字幕| 国产黄色视频一区二区在线观看| 99久久精品国产国产毛片| 国产精品一区二区三区四区久久| 全区人妻精品视频| 极品教师在线视频| 免费av毛片视频| 国产 亚洲一区二区三区 | 网址你懂的国产日韩在线| 18禁动态无遮挡网站| 亚洲av.av天堂| 亚洲电影在线观看av| 国产成人精品一,二区| 老师上课跳d突然被开到最大视频| 夜夜爽夜夜爽视频| 日韩视频在线欧美| av黄色大香蕉| 伦精品一区二区三区| 亚洲国产精品国产精品| 久久久久久久国产电影| 麻豆乱淫一区二区| 久久韩国三级中文字幕| 欧美日韩国产mv在线观看视频 | 插阴视频在线观看视频| 欧美成人一区二区免费高清观看| 青青草视频在线视频观看| or卡值多少钱| 久久99热这里只有精品18| 97超视频在线观看视频| 五月伊人婷婷丁香| 国产免费一级a男人的天堂| 永久免费av网站大全| 亚洲欧美日韩东京热| 色5月婷婷丁香| 十八禁国产超污无遮挡网站| 女人被狂操c到高潮| 午夜激情欧美在线| 好男人视频免费观看在线| 老司机影院成人| 亚洲国产成人一精品久久久| 99热这里只有是精品在线观看| 国模一区二区三区四区视频| 日韩av在线大香蕉| 99久久中文字幕三级久久日本| 美女xxoo啪啪120秒动态图| 国产亚洲av嫩草精品影院| 免费黄色在线免费观看| 麻豆乱淫一区二区| 最近2019中文字幕mv第一页| 非洲黑人性xxxx精品又粗又长| 国产精品美女特级片免费视频播放器| 午夜激情欧美在线| 韩国高清视频一区二区三区| 在现免费观看毛片| 观看美女的网站| 26uuu在线亚洲综合色| 国产成人免费观看mmmm| 菩萨蛮人人尽说江南好唐韦庄| 日本午夜av视频| 亚洲av日韩在线播放| 人妻一区二区av| 一个人看视频在线观看www免费| 波野结衣二区三区在线| 国内精品美女久久久久久| 街头女战士在线观看网站| 日日啪夜夜撸| 欧美日韩精品成人综合77777| 91av网一区二区| 国产真实伦视频高清在线观看| 爱豆传媒免费全集在线观看| 国产伦精品一区二区三区视频9| 亚洲国产精品sss在线观看| 午夜福利网站1000一区二区三区| 男的添女的下面高潮视频| 日本熟妇午夜| 国产黄片美女视频| 国产黄片视频在线免费观看| 18禁在线无遮挡免费观看视频| 一本久久精品| 精品一区二区三区视频在线| 男人爽女人下面视频在线观看| 卡戴珊不雅视频在线播放| 一夜夜www| 美女黄网站色视频| 中文字幕免费在线视频6| 精品欧美国产一区二区三| 男人爽女人下面视频在线观看| 校园人妻丝袜中文字幕| 国产极品天堂在线| 免费观看av网站的网址| 一个人看视频在线观看www免费| 久久久久精品性色| 少妇猛男粗大的猛烈进出视频 | 伊人久久国产一区二区| 久久久久久久亚洲中文字幕| 99久久人妻综合| 久99久视频精品免费| 国产亚洲一区二区精品| 国产大屁股一区二区在线视频| 少妇的逼水好多| 日日啪夜夜爽| 中文字幕av在线有码专区| 国产成人一区二区在线| 久久草成人影院| 国产成年人精品一区二区| .国产精品久久| 亚洲久久久久久中文字幕| 午夜免费男女啪啪视频观看| 亚洲经典国产精华液单| 久久久精品免费免费高清| 国国产精品蜜臀av免费| 欧美日本视频| 超碰97精品在线观看| 免费不卡的大黄色大毛片视频在线观看 | 97精品久久久久久久久久精品| 久久久久精品性色| 一级毛片我不卡| 日韩成人av中文字幕在线观看| 国产精品99久久久久久久久| 亚洲欧洲日产国产| 永久免费av网站大全| 国产成人免费观看mmmm| 精品久久久精品久久久| 直男gayav资源| videos熟女内射| 午夜久久久久精精品| 国产综合懂色| 成人综合一区亚洲| 赤兔流量卡办理| 伦精品一区二区三区| 黄色配什么色好看| 免费无遮挡裸体视频| 九九在线视频观看精品| 欧美日韩一区二区视频在线观看视频在线 | 人妻少妇偷人精品九色| 在线免费十八禁| 人妻少妇偷人精品九色| 免费看av在线观看网站| 亚洲三级黄色毛片| 欧美激情国产日韩精品一区| 日本免费在线观看一区| 黄片无遮挡物在线观看| 久久久久久伊人网av| 亚州av有码| or卡值多少钱| 久久久色成人| 日韩欧美一区视频在线观看 | 中文字幕免费在线视频6| av在线亚洲专区| 不卡视频在线观看欧美| 亚洲国产精品sss在线观看| 热99在线观看视频| 天美传媒精品一区二区| 精品久久久精品久久久| 精品一区二区三卡| 亚洲av国产av综合av卡| 国产综合精华液| 欧美成人精品欧美一级黄| 精品不卡国产一区二区三区| av网站免费在线观看视频 | 国产极品天堂在线| or卡值多少钱| 毛片一级片免费看久久久久| 精品一区二区三区视频在线| 久热久热在线精品观看| 一级毛片黄色毛片免费观看视频| 日韩强制内射视频| 成年av动漫网址| 亚洲精品aⅴ在线观看| 中国国产av一级| 国产久久久一区二区三区| 亚洲精品aⅴ在线观看| 两个人视频免费观看高清| 日韩视频在线欧美| 日本午夜av视频| 国产精品嫩草影院av在线观看| 国产伦一二天堂av在线观看| 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| 精品不卡国产一区二区三区| 色哟哟·www| 国产精品蜜桃在线观看| 色综合色国产| 亚洲精品乱久久久久久| 亚洲精品国产av成人精品| 直男gayav资源| 国产精品一区二区三区四区免费观看| 久久久久精品久久久久真实原创| 精品一区在线观看国产| 婷婷色av中文字幕| 白带黄色成豆腐渣| 亚洲精品乱码久久久久久按摩| 久久人人爽人人爽人人片va| 久久久久精品性色| 小蜜桃在线观看免费完整版高清| 国产极品天堂在线| 午夜亚洲福利在线播放| 中文字幕制服av| 欧美日韩综合久久久久久| 高清日韩中文字幕在线| 久久久久精品性色| 联通29元200g的流量卡| 国产精品久久久久久久电影| 最新中文字幕久久久久| 国产精品熟女久久久久浪| 蜜桃久久精品国产亚洲av| kizo精华| 在线观看一区二区三区| 久久精品综合一区二区三区| 久久国内精品自在自线图片| 免费观看精品视频网站| 成人亚洲精品一区在线观看 | 日韩欧美 国产精品| 一级片'在线观看视频| 久久久久精品性色| 黄色一级大片看看| 亚洲内射少妇av| 久久人人爽人人片av| 欧美激情国产日韩精品一区| 最新中文字幕久久久久| 性色avwww在线观看| 国产精品麻豆人妻色哟哟久久 | 偷拍熟女少妇极品色| 国产精品国产三级国产专区5o| 男女那种视频在线观看| 国产欧美日韩精品一区二区| 国产精品1区2区在线观看.| 久久热精品热| 欧美日韩精品成人综合77777| 搡老乐熟女国产| 国产亚洲一区二区精品| 国产成人精品一,二区| 天天一区二区日本电影三级| 久久久久久久久久久免费av| 毛片女人毛片| 青春草国产在线视频| 国产伦一二天堂av在线观看| 看非洲黑人一级黄片| 男插女下体视频免费在线播放| 午夜免费观看性视频| 精品久久久久久久久亚洲| 国产美女午夜福利| 夫妻性生交免费视频一级片| 听说在线观看完整版免费高清| 成年av动漫网址| 纵有疾风起免费观看全集完整版 | 在线观看美女被高潮喷水网站| 熟妇人妻不卡中文字幕| 男的添女的下面高潮视频| 日韩av在线大香蕉| 国产伦精品一区二区三区四那| 精品久久久久久久久亚洲| 乱系列少妇在线播放| 国产精品.久久久| 日韩,欧美,国产一区二区三区| 深爱激情五月婷婷| 波多野结衣巨乳人妻| 美女黄网站色视频| 色尼玛亚洲综合影院| 日韩人妻高清精品专区| 欧美成人一区二区免费高清观看| 国产成人精品一,二区| 国产精品国产三级专区第一集| 性插视频无遮挡在线免费观看| 国产乱人偷精品视频| 国产爱豆传媒在线观看| 91久久精品国产一区二区成人| 麻豆成人av视频| 国产精品三级大全| 亚洲天堂国产精品一区在线| 91av网一区二区| 麻豆成人午夜福利视频| 久久久久久久久久人人人人人人| 青青草视频在线视频观看| a级一级毛片免费在线观看| 精品国产三级普通话版| 日韩一区二区三区影片| 美女脱内裤让男人舔精品视频| 国产黄片视频在线免费观看| 大又大粗又爽又黄少妇毛片口| 日韩不卡一区二区三区视频在线| 亚洲精品456在线播放app| 亚洲电影在线观看av| 日日干狠狠操夜夜爽| 欧美变态另类bdsm刘玥| xxx大片免费视频| 久久精品国产鲁丝片午夜精品| 亚洲精品视频女| 大片免费播放器 马上看| 久久精品熟女亚洲av麻豆精品 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲av电影不卡..在线观看| 精品人妻视频免费看| 亚洲精品成人av观看孕妇| 男的添女的下面高潮视频| 高清欧美精品videossex| 人体艺术视频欧美日本| 欧美zozozo另类| 色综合色国产| 91精品伊人久久大香线蕉| 亚洲在久久综合| 国产综合精华液| 午夜精品一区二区三区免费看| 国产精品一及| 亚洲人成网站在线播| 少妇人妻一区二区三区视频| 啦啦啦啦在线视频资源| 一区二区三区高清视频在线| 18禁在线无遮挡免费观看视频| 成年版毛片免费区| 最近中文字幕高清免费大全6| 熟妇人妻久久中文字幕3abv| 国产精品一区二区三区四区久久| 国产91av在线免费观看| 久久精品人妻少妇| 成人av在线播放网站| 国内少妇人妻偷人精品xxx网站| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| 少妇裸体淫交视频免费看高清| 久久久欧美国产精品| 黑人高潮一二区| 亚洲av在线观看美女高潮| 免费观看a级毛片全部| 51国产日韩欧美| 大香蕉97超碰在线| 大香蕉久久网| 超碰av人人做人人爽久久| kizo精华| 舔av片在线| 80岁老熟妇乱子伦牲交| 免费av毛片视频| 建设人人有责人人尽责人人享有的 | 日本欧美国产在线视频| 中文资源天堂在线| 国产成人福利小说| 天堂影院成人在线观看| 最新中文字幕久久久久| 一级av片app| 日本与韩国留学比较| 非洲黑人性xxxx精品又粗又长| 一区二区三区高清视频在线| 久久国内精品自在自线图片| 欧美最新免费一区二区三区| 免费观看无遮挡的男女| 国产探花极品一区二区| 视频中文字幕在线观看| 亚洲,欧美,日韩| 日韩强制内射视频| 在线播放无遮挡| 亚洲最大成人av| 国产精品人妻久久久久久| 久久久久久国产a免费观看| videos熟女内射| 日日撸夜夜添| 高清毛片免费看| 久久这里只有精品中国| 熟妇人妻不卡中文字幕| 男女下面进入的视频免费午夜| 成人鲁丝片一二三区免费| 97人妻精品一区二区三区麻豆| 国产精品久久久久久久久免| av播播在线观看一区| 国产乱来视频区| 久久久午夜欧美精品| 久99久视频精品免费| 久久久久久九九精品二区国产| 亚洲国产色片| 国产成人精品久久久久久| 亚洲四区av| 日本-黄色视频高清免费观看| 午夜免费激情av| 精品国产一区二区三区久久久樱花 | 精品久久国产蜜桃| 小蜜桃在线观看免费完整版高清| 七月丁香在线播放| 狠狠精品人妻久久久久久综合| 黄片无遮挡物在线观看| 午夜激情福利司机影院| 精品熟女少妇av免费看| 亚洲婷婷狠狠爱综合网| 色网站视频免费| 欧美zozozo另类| 中国美白少妇内射xxxbb| 一级二级三级毛片免费看| 青春草视频在线免费观看| 亚洲一级一片aⅴ在线观看| 直男gayav资源| 日本wwww免费看| 99热这里只有精品一区| 亚洲欧美精品专区久久| 久久99热6这里只有精品| 成人毛片a级毛片在线播放| 久久久久九九精品影院| 老司机影院毛片| 春色校园在线视频观看|