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

    基于連續(xù)成像的頸椎運(yùn)動特征

    2017-07-31 21:54:17裴欣欣羅嗣頻萬業(yè)達(dá)王學(xué)民天津大學(xué)精密儀器與光電子工程學(xué)院天津30007天津醫(yī)院放射科天津300天津市生物醫(yī)學(xué)檢測技術(shù)與儀器重點(diǎn)實(shí)驗(yàn)室天津30007
    關(guān)鍵詞:頸椎病頸椎受試者

    裴欣欣,羅嗣頻,萬業(yè)達(dá),王學(xué)民,3,周 鵬,3*(.天津大學(xué)精密儀器與光電子工程學(xué)院,天津 30007;.天津醫(yī)院放射科,天津 300; 3.天津市生物醫(yī)學(xué)檢測技術(shù)與儀器重點(diǎn)實(shí)驗(yàn)室,天津 30007)

    基于連續(xù)成像的頸椎運(yùn)動特征

    裴欣欣1,羅嗣頻2,萬業(yè)達(dá)2,王學(xué)民1,3,周 鵬1,3*
    (1.天津大學(xué)精密儀器與光電子工程學(xué)院,天津 300072;2.天津醫(yī)院放射科,天津 300211; 3.天津市生物醫(yī)學(xué)檢測技術(shù)與儀器重點(diǎn)實(shí)驗(yàn)室,天津 300072)

    目的 在獲取X線頸椎屈伸運(yùn)動連續(xù)圖像的基礎(chǔ)上探討頸椎的運(yùn)動特征。方法 對頸椎病患者及正常人頸椎過屈位至過伸位自發(fā)連續(xù)運(yùn)動的多幀矢狀位X線圖像進(jìn)行圖像預(yù)處理、配準(zhǔn)后,記錄單椎體的運(yùn)動軌跡,提取各椎體的特征點(diǎn),并計算出活躍椎體的相對幾何學(xué)參數(shù)、運(yùn)動速率。結(jié)果 患者組與正常組C4、5、6椎體動態(tài)軌跡不同。與正常組比較,過屈位時患者組C4、5椎體左側(cè)邊緣曲線與基準(zhǔn)線的夾角較小(P均<0.05)?;颊呓MC4椎體的運(yùn)動存在不穩(wěn)定性,角度變化的波動性、速率改變的幅度及頻率均較大。且患者組相鄰單椎體相對位置的改變更小。結(jié)論 通過本研究圖像預(yù)處理、配準(zhǔn)、參數(shù)提取及結(jié)果分析方法,可真實(shí)反映頸椎病變帶來的功能學(xué)特征改變,為基于X線圖像的頸椎動態(tài)特征分析提供新思路。

    頸椎;動態(tài)特征;運(yùn)動速率;單椎體運(yùn)動

    傳統(tǒng)的頸椎病檢查多采用X線進(jìn)行形態(tài)學(xué)特征分析[1-2],而頸椎病初期是功能性的,靜態(tài)觀察不能很好地反映其特征。以往通過有限元法[3]可追蹤到靜態(tài)圖像椎體邊緣,基于雙目立體視覺的頸椎運(yùn)動測量[4]、光學(xué)追蹤[5]和電磁追蹤[6-7]等均依賴于外部標(biāo)記運(yùn)動追蹤,所得數(shù)據(jù)誤差較大。Anderst等[8]將直徑1.6 mm的鈦珠植入研究對象的椎體和棘突內(nèi)以獲取椎體運(yùn)動特征,但方法具有侵入性。葉斌[9]通過攝取頸椎張口位、左右旋轉(zhuǎn)后張口位研究頸椎的動態(tài)運(yùn)動。但此類研究多采用不同位置的靜態(tài)矢狀位圖像,分別對圖像進(jìn)行測量分析,其動態(tài)特征不明顯。本研究通過訓(xùn)練受試者自主完成頸椎屈伸運(yùn)動,選取頸椎過屈位至過伸位自發(fā)連續(xù)運(yùn)動的多幀矢狀位X線圖像,經(jīng)預(yù)處理、配準(zhǔn)后,提取頸椎椎體特征點(diǎn),通過與相鄰圖像的連續(xù)處理,描繪每一椎體的整體運(yùn)動的軌跡,計算椎體運(yùn)動速率,分析感興趣椎體動態(tài)功能特征和相鄰兩椎體在運(yùn)動過程中相對位置的改變,以期更真實(shí)地反映頸椎病變引起的頸椎動態(tài)功能的改變,進(jìn)一步加深對生理載荷下頸椎運(yùn)動狀況的認(rèn)識。

    1 試驗(yàn)流程

    在獲取頸椎屈伸運(yùn)動連續(xù)圖像的基礎(chǔ)上,對圖像進(jìn)行預(yù)處理,以增強(qiáng)頸椎的結(jié)構(gòu),提取出每一椎骨的特征點(diǎn)。由于受試者進(jìn)行頸椎屈伸運(yùn)動時,C7椎體常被遮擋,X線圖像中難以區(qū)分,為清楚地了解各椎骨的相互關(guān)系,以C6椎體圖像進(jìn)行配準(zhǔn)處理。完成上述處理后,對相鄰的圖像進(jìn)行參數(shù)提取和計算,分析正常人群與頸椎病患者頸椎運(yùn)動速率的改變,并記錄單椎體的運(yùn)動軌跡和相對位置的變化,試驗(yàn)流程見圖1。

    圖1 試驗(yàn)流程圖

    2 圖像獲取

    于2015年10月—12月征集21名受試者納入研究,包括臨床診斷為頸椎病患者14例(患者組)及健康志愿者7名(正常組)?;颊呓M中,男4例,女10例,年齡22~54歲,中位年齡37.0歲;正常組中,男2例,女5例,年齡23~45歲,中位年齡34.0歲。所有受試者對試驗(yàn)風(fēng)險均知情同意,自愿參與試驗(yàn)。

    采用島津Sonial vision SafireⅡ平板數(shù)字化X線多功能透視攝影系統(tǒng)進(jìn)行頸椎動態(tài)檢測。設(shè)置管電流為2 mA,管電壓為60 kV,X球管與受試者頸部皮膚間的距離為1.2 m,調(diào)整成像中心于受試者頸椎C4椎體。為觀察頸椎的運(yùn)動狀態(tài),X線攝片時囑受試者側(cè)位站立,雙手叉腰,肘向后靠攏,首先進(jìn)行低頭至過屈位,而后在15 s內(nèi)從過屈位勻速運(yùn)動至過伸位。將X線機(jī)置于“準(zhǔn)備”態(tài),在受試者頸椎以勻速從過屈位運(yùn)動至過伸位的過程中進(jìn)行X線曝光,圖像采集速率為6幀/秒,共采集15 s[10]。所獲圖像為1179×782像素,將圖像按順序編號命名,取每一受試者8幀不同位置圖像,見圖2。

    3 圖像預(yù)處理及配準(zhǔn)

    3.1選取ROI 為方便后續(xù)的參數(shù)提取,對所獲取圖像進(jìn)行預(yù)處理,在Matlab中讀取已保存的圖像,選取涵蓋頸椎過屈位至過伸位的全部位置區(qū)域的最小圖像作為ROI,記為W,大小為u×v像素,其中u為橫向像素個數(shù),v為縱向像素個數(shù),對每一受試者選取的ROI圖像大小不同,但需使每個位置的頸椎均在ROI內(nèi),圖3所示即為從圖2的圖像中選取的ROI。

    3.2圖像增強(qiáng) 為提高圖像對比度,調(diào)整圖像的灰度直方圖分布,使頸椎部分更為突出且減少其周圍組織和噪聲的影響,根據(jù)圖像自身的特點(diǎn),選取基于形態(tài)學(xué)的頂帽、底帽變換圖像增強(qiáng)算法和基于Retinex色彩理論[11]的圖像增強(qiáng)算法進(jìn)行雙增強(qiáng)。

    頂帽、底帽變換的一個重要用途是校正不均勻光照[12],在多張圖像的處理中,該算法對背景光照進(jìn)行了校正,見圖3。Retinex色彩理論是基于人類視覺基礎(chǔ)的理論,該算法具有銳化圖像、顏色恒常性、動態(tài)壓縮范圍大、可對不同圖像進(jìn)行自適應(yīng)增強(qiáng)等優(yōu)點(diǎn)[11],但對噪聲較敏感[12]。聯(lián)合兩種方法可使圖像增強(qiáng)效果更好(圖4),以便頸椎邊緣的完整提取。

    3.3邊緣提取 為提取頸椎各椎體左側(cè)邊緣曲線得到頸椎運(yùn)動軌跡,從而觀察整體運(yùn)動趨勢,采用性能優(yōu)良Canny算子[13]初步提取邊緣,見圖5。

    圖2 受試者頸椎過屈位(A~D)運(yùn)動至過伸位(E~H)X線圖像

    圖3 頂帽、底帽變換圖像增強(qiáng)算法后圖像 圖4 Retinex色彩理論增強(qiáng)后圖像 圖5 采用Canny算子初步提取頸椎邊緣曲線 圖6 修正后提取的椎體左側(cè)邊緣曲線

    在第1幀圖像的頸椎椎體邊緣選取m個種子點(diǎn),并將各點(diǎn)坐標(biāo)保存于m×2的矩陣A(x,y)內(nèi)。在第2幀圖像內(nèi)遍歷像素點(diǎn),在m個種子點(diǎn)的增長區(qū)域內(nèi)尋找像素值為1的點(diǎn),并將此點(diǎn)的坐標(biāo)保存在與圖像大小相等的矩陣內(nèi)。種子點(diǎn)(xij,yij)的增長區(qū)域?yàn)椋?/p>

    橫坐標(biāo):(xij-a1,xij+a2),i=1,2,3,…,n;j=1,2,3,…,m

    縱坐標(biāo):(yij-b1,yij+b2),i=1,2,3,…,n;j=1,2,3,…,m

    其中,a1,a2,b1,b2由多次實(shí)驗(yàn)得出,范圍為3~7,不同受試者的參數(shù)選擇不同。n為圖像的幀數(shù),以包括過屈位和過伸位為邊界,一般為80~90幀。初次完成后,對提取的邊緣作出修正,此時可較光滑地提取頸椎椎體邊緣的左側(cè)曲線,見圖6。左側(cè)曲線的像素點(diǎn)即為新的種子點(diǎn),將其保存于矩陣A內(nèi),作為下一幀圖像的遍歷的種子點(diǎn),各幀圖像均以前一幀圖像頸椎左側(cè)邊緣曲線為種子點(diǎn)進(jìn)行遍歷,各幀圖像依次進(jìn)行,并將各幀所得矩陣A內(nèi)的種子點(diǎn)坐標(biāo)均保存于數(shù)組C內(nèi),數(shù)組C大小與W相同。

    3.4邊緣曲線配準(zhǔn) 為觀察頸椎整體運(yùn)動的范圍和趨勢變化,減少在拍攝時受試者肩部位置改變的影響,選取第1幀圖像頸椎的C6椎體左下側(cè)頂點(diǎn)為基點(diǎn)進(jìn)行配準(zhǔn),在提取完頸椎椎體左側(cè)曲線后,此點(diǎn)即為縱坐標(biāo)y值的最大點(diǎn),坐標(biāo)記為M1(X1,Y1),第2幀圖像相同位置的點(diǎn)坐標(biāo)記為M2(X2,Y2),以此類推,最后一幀相同位置點(diǎn)坐標(biāo)記為Mn(Xn,Yn)。

    配準(zhǔn)距離差為:Di(X,Y)=Mi(Xi,Yi)-M1(X1,Y1);i=2,3,4…,n

    配準(zhǔn)過程為:Bi(xij,yij)=Ai(xij,yij)-Di(X,Y);i=2,3,4,…,n;j=1,2,3,…,m

    表1 患者組與正常組C4、5椎體過屈位、過伸位的角度比較(°)

    圖7 C4、5椎體左邊緣曲線 A.正常人; B.頸椎病患者

    此時,頸椎椎體的左側(cè)曲線全部以C6椎體的左下側(cè)頂點(diǎn)為基準(zhǔn)點(diǎn)配準(zhǔn),觀察頸椎整體的運(yùn)動趨勢、頸椎的生理曲線狀況及頸椎后伸位的最大位置所在,見圖7。

    4 參數(shù)提取與結(jié)果分析

    將各受試者的圖像歸一化為86幀,以C4、5椎體為感興趣節(jié)段,分別以其第1幀左側(cè)邊緣曲線為起點(diǎn)選擇種子點(diǎn),用上述方法處理(圖7),紅色代表配準(zhǔn)前所在位置,其他不同顏色代表各幀圖像各椎體左側(cè)邊緣所在位置。

    患者組與正常組C4、5、6椎體動態(tài)軌跡不同,患者組頸椎生理彎曲消失,且C4、5椎體間隙較小,且在頸椎屈伸運(yùn)動過程中椎體間隙大小發(fā)生明顯變化。

    記錄單椎體運(yùn)動,選取數(shù)組C內(nèi)第5~10行所有數(shù)據(jù),按列求取平均值的相對固定點(diǎn)坐標(biāo)Ni(xi,yi),

    i=2,3,…,n。為觀察單椎體的運(yùn)動趨勢,以配準(zhǔn)點(diǎn)的橫坐標(biāo)作基準(zhǔn)線,計算椎體左側(cè)邊緣曲線與基準(zhǔn)線的夾角:

    i=2,3,4,…,n

    患者組及正常組的C4與C5椎體的夾角及其在運(yùn)動過程中的變化見圖8、9。

    以C4椎體分析頸椎的運(yùn)動速率,其角度變化見圖8B;對離散數(shù)據(jù)進(jìn)行差分運(yùn)算,其速率變化見圖10。由于頸椎病患者不能很好地控制頸椎的運(yùn)動從而產(chǎn)生被動的速率變化,表現(xiàn)為運(yùn)動的不穩(wěn)定性?;颊呓M角度變化的波動性大,其速率改變的幅度及改變頻率均較大(圖10)。患者組與正常組頸椎運(yùn)動速率方差對比見圖11,患者組方差為0.462,正常組為0.331。

    觀察單椎體運(yùn)動軌跡(圖8、9),分析患者組與正常組C4、5椎體左側(cè)邊緣曲線與基準(zhǔn)線的夾角,對患者組和正常組C4、5椎體過屈位、過伸位的角度采用SPSS軟件進(jìn)行兩獨(dú)立樣本t檢驗(yàn)結(jié)果見表1?;颊呓M與正常組C4、5椎體在過屈位時角度有差異有統(tǒng)計學(xué)意義(P均<0.05),而在過伸位時差異無統(tǒng)計學(xué)意義(P均>0.05)。

    分析相鄰單椎體相對位置的改變時,將C4、5椎體角度的變化相減,計算C4、5椎體相對位置變化的數(shù)據(jù)并擬合出模板?;颊呓M及正常組C4、5椎體相對位置變化數(shù)據(jù)比較見圖12,紅色為正常組模板數(shù)據(jù),藍(lán)色為患者組的數(shù)據(jù),可見患者組椎體相對位置的改變較小,曲線較正常組平緩且最低點(diǎn)絕對值較小。

    圖8 C4椎體運(yùn)動軌跡 A.C4椎體夾角示意圖; B.C4椎體角度變化對比圖 圖9 C5椎體運(yùn)動軌跡 A.C5椎體夾角示意圖; B.C5椎體角度變化對比圖

    圖10 正常人與患者椎體角度變化速率對比圖 圖11 正常人與患者椎體角度變化速率方差對比圖 圖12 正常人與患者C4與C5椎體相對位置變化

    5 結(jié)語

    本研究在對圖像進(jìn)行預(yù)處理的基礎(chǔ)上對頸椎椎體左側(cè)邊緣進(jìn)行提取分析,以配準(zhǔn)點(diǎn)橫坐標(biāo)為基準(zhǔn)線,觀察C4、5椎體在自發(fā)性頸椎屈伸運(yùn)動過程中椎體左側(cè)邊緣曲線與基準(zhǔn)線夾角的變化趨勢,分析頸椎運(yùn)動速率,發(fā)現(xiàn)頸椎病患者頸椎運(yùn)動過程中速率呈跳躍性變化。此外,本研究以C4、5椎體左側(cè)邊緣與基準(zhǔn)線的夾角的差值反映單椎體在運(yùn)動過程中相對位置的改變,基于連續(xù)圖像的動態(tài)特征分析,為單椎體運(yùn)動的研究提供更多數(shù)據(jù),也為頸椎病診斷提供更多參考信息。但本研究的樣本量較少,且因存在個體差異,不同患者C4、5椎體在運(yùn)動過程中相對位置變化曲線不同,該差異還有待進(jìn)一步深入研究。

    [1] Nagamoto Y, Ishii T, Sakaura H, et al. In vivo three-dimensional kinematics of the cervical spine during head rotation in patients with cervical spondylosis. Spine (Phila Pa 1976), 2011,36(10):778-783.

    [2] Nagamoto Y, Iwasaki M, Sugiura T, et al. In vivo 3D kinematic changes in the cervical spine after laminoplasty for cervical spondylotic myelopathy. J Neurosurg Spine, 2014,21(3):417-424.

    [3] Xu X, Hao HW, Yin XC, et al. Automatic segmentation of cervical vertebrae in X-ray images. International Joint Conference on Neural Networks, 2012,20:1-8.

    [4] 鄭偉龍,王志勇,王世斌,等.頸椎三維運(yùn)動測量及生物力學(xué)特性研究//北方七省市區(qū)力學(xué)學(xué)術(shù)會議,2014:70-71.

    [5] Lee JH, Kim JS, Lee JH, et al. Comparison of cervical kinematics between patients with cervical artificial disc replacement and anterior cervical discectomy and fusion for cervical disc herniation. Spine J, 2014,14(7):1199-1204.

    [6] Bell KM, Bechara BP, Hartman RA, et al. Influence of number of operated levels and postoperative time on active range of motion following anterior cervical decompression and fusion procedures. Spine (Phila Pa 1976), 2011,36(4):263-268.

    [7] Bechara BP, Bell KM, Hartman RA, et al. In vivo analysis of cervical range of motion after 4- and 5-level subaxial cervical spine fusion. Spine (Phila Pa 1976), 2012,37(1):E23-E29.

    [8] Anderst WJ, Vaidya R, Tashman S. A technique to measure three-dimensional in vivo rotation of fused and adjacent lumbar vertebrae. Spine J, 2008,8(6):991-997.

    [9] 葉斌.動態(tài)X線片在頸椎病早期診斷中的作用.現(xiàn)代診斷與治療,2013,24(5):1130-1131.

    [10] Bifulco P, Cesarelli M, Romano M, et al. Measurement of intervertebral cervical motion by means of dynamic X-ray image processing and data interpolation. Int J Biomed Imaging, 2013,2013:152920.

    [11] 陳琛,張建州.改進(jìn)的多尺度Retinex醫(yī)學(xué)X射線圖像增強(qiáng)算法.計算機(jī)工程與應(yīng)用,2015,51(9):191-195.

    [12] 馬超玉.光照不均勻條件下圖像增強(qiáng)算法研究.長春:長春理工大學(xué), 2014:4-13.

    [13] Wang XJ, Liu XM, Yong G. Image edge detection algorithm based on improved canny operator. Computer Engineering, 2012,38(14):168-172.

    Motion characteristics of cervical spine based on continuous images

    PEIXinxin1,LUOSipin2,WANYeda2,WANGXuemin1,3,ZHOUPeng1,3*
    (1.SchoolofPrecisionInstrumentandOpto-ElectronicsEngineeringofTianjinUniversity,Tianjin300072,China; 2.DepartmentofRadiology,TianjinHospital,Tianjin300211,China; 3.TianjinKeyLaboratoryofBiomedicalTestingTechnologyandInstruments,Tianjin300072,China)

    Objective To explore the motion features of cervical spine based on continuous X-ray images. Methods The cervical spontaneous continuous multi-frame sagittal images from flexion to extention positions were selected from cervical spondylosis patients (patients group) and healthy adult (normal group). After preprocessing and registration, the trajectory of single vertebral body were recorded, and the feature points of each vertebral body were extracted. Meanwhile, the relative geometry parameters and movement rate were calculated. Results The motion trajectory of patients' cervical spine C4—6 were different from healthy people. The angles between the left edge curve of the vertebral body (C4 and C5 vertebra) and the baseline of patients group were smaller than those of normal group in flexion position (allP<0.05). There were instability in the movement of C4 vertebral body in patient group, and the volatility of the angle change, the rate of change and the frequency were larger. The relative position change of the adjacent single vertebral body in the patient group are smaller. Conclusion Through the preprocessing, registration, parameter extraction and result analysis, the changes of functional features in cervical spondylosis patients are truly reflected. And it also provides a new idea for dynamic analysis of cervical vertebrae based on X-ray images.

    Cervical vertebral; Dynamic characteristics; Movement speed; Single vertebral body motion

    裴欣欣(1991—),女,河北邢臺人,在讀碩士。研究方向:醫(yī)學(xué)圖像處理。E-mail: peixx@tju.edu.cn

    周鵬,天津大學(xué)精密儀器與光電子工程學(xué)院,300072; 天津市生物醫(yī)學(xué)檢測技術(shù)與儀器重點(diǎn)實(shí)驗(yàn)室, 300072。E-mail: zpzp@tju.edu.cn

    2016-10-17

    2017-04-10

    10.13929/j.1003-3289.201610065

    TP391.5

    A

    1003-3289(2017)07-1090-05

    猜你喜歡
    頸椎病頸椎受試者
    頸椎不能承受之重
    大自然探索(2024年1期)2024-02-29 09:11:26
    涉及人的生物醫(yī)學(xué)研究應(yīng)遵循的倫理原則
    涉及人的生物醫(yī)學(xué)研究應(yīng)遵循的倫理原則
    頸椎為什么易發(fā)病
    養(yǎng)護(hù)頸椎,守住健康
    頸椎病與老年癡呆
    涉及人的生物醫(yī)學(xué)研究應(yīng)遵循的倫理原則
    游泳 趕走頸椎病
    拿什么拯救你,我的頸椎
    涉及人的生物醫(yī)學(xué)研究應(yīng)遵循的倫理原則
    午夜激情福利司机影院| 欧美日韩亚洲综合一区二区三区_| 久久久国产精品麻豆| 制服丝袜大香蕉在线| 一级片免费观看大全| 最近在线观看免费完整版| 国产亚洲欧美98| 在线播放国产精品三级| 可以在线观看毛片的网站| 韩国av一区二区三区四区| 免费在线观看视频国产中文字幕亚洲| 久热爱精品视频在线9| 一区二区三区激情视频| 人人妻,人人澡人人爽秒播| 男女下面进入的视频免费午夜| 久久中文看片网| 极品教师在线免费播放| 男女之事视频高清在线观看| 国产免费男女视频| 高清在线国产一区| 久久国产精品人妻蜜桃| 午夜日韩欧美国产| 国产乱人伦免费视频| 国产99白浆流出| 久久久久亚洲av毛片大全| av国产免费在线观看| 国产精品av视频在线免费观看| 黄色视频不卡| 岛国在线免费视频观看| 青草久久国产| 国产区一区二久久| 精品久久久久久久久久久久久| 国产精品电影一区二区三区| 久久久久国产精品人妻aⅴ院| 又紧又爽又黄一区二区| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲av香蕉五月| 男女床上黄色一级片免费看| 丝袜美腿诱惑在线| 又大又爽又粗| e午夜精品久久久久久久| 久久久久精品国产欧美久久久| 免费一级毛片在线播放高清视频| 久久久久国产精品人妻aⅴ院| 我的老师免费观看完整版| 久久久久亚洲av毛片大全| 亚洲国产高清在线一区二区三| 日本免费一区二区三区高清不卡| 伊人久久大香线蕉亚洲五| 在线观看一区二区三区| 黄色视频不卡| 老鸭窝网址在线观看| 一级作爱视频免费观看| 久久久久国产一级毛片高清牌| 亚洲av成人精品一区久久| 香蕉久久夜色| 精品国产超薄肉色丝袜足j| 男女午夜视频在线观看| 国产精品,欧美在线| 极品教师在线免费播放| 午夜福利视频1000在线观看| 欧美3d第一页| 亚洲在线自拍视频| www日本在线高清视频| 亚洲成人久久性| 国产又色又爽无遮挡免费看| 日韩三级视频一区二区三区| 免费人成视频x8x8入口观看| 精品熟女少妇八av免费久了| www.精华液| 亚洲午夜理论影院| 亚洲熟女毛片儿| 亚洲国产精品sss在线观看| 在线免费观看的www视频| 久久人妻福利社区极品人妻图片| 国产亚洲av嫩草精品影院| 又爽又黄无遮挡网站| 脱女人内裤的视频| av免费在线观看网站| 婷婷精品国产亚洲av在线| 国产精品久久久久久精品电影| 亚洲av电影在线进入| 俺也久久电影网| 老司机午夜福利在线观看视频| 中亚洲国语对白在线视频| 老司机午夜十八禁免费视频| 制服人妻中文乱码| 欧美日韩国产亚洲二区| 亚洲精品中文字幕在线视频| 国产不卡一卡二| 91成年电影在线观看| 十八禁人妻一区二区| 可以免费在线观看a视频的电影网站| 欧美一级a爱片免费观看看 | 久久精品aⅴ一区二区三区四区| or卡值多少钱| 91麻豆精品激情在线观看国产| 全区人妻精品视频| 精品国产亚洲在线| 国产精品一区二区三区四区免费观看 | 一区福利在线观看| 亚洲一区中文字幕在线| 午夜两性在线视频| 亚洲国产看品久久| 欧美日韩乱码在线| 国产日本99.免费观看| 国产精品电影一区二区三区| 一级黄色大片毛片| av天堂在线播放| 欧美av亚洲av综合av国产av| 国产精品一及| 久久久国产精品麻豆| a在线观看视频网站| 欧美国产日韩亚洲一区| 男男h啪啪无遮挡| 欧美精品啪啪一区二区三区| 午夜视频精品福利| 久久精品夜夜夜夜夜久久蜜豆 | 欧美zozozo另类| 国产欧美日韩一区二区精品| 狂野欧美激情性xxxx| 日日爽夜夜爽网站| 国产精品久久久久久亚洲av鲁大| 中文字幕人妻丝袜一区二区| 久久久水蜜桃国产精品网| 国产亚洲av高清不卡| 999久久久精品免费观看国产| 草草在线视频免费看| 亚洲欧美精品综合一区二区三区| 人妻久久中文字幕网| a级毛片a级免费在线| 免费在线观看完整版高清| 国产亚洲精品久久久久5区| 免费搜索国产男女视频| 桃色一区二区三区在线观看| 免费一级毛片在线播放高清视频| 欧美日韩精品网址| 国内少妇人妻偷人精品xxx网站 | 国产亚洲精品第一综合不卡| 亚洲av日韩精品久久久久久密| 麻豆国产av国片精品| xxx96com| 亚洲熟女毛片儿| 久久这里只有精品中国| 国产黄色小视频在线观看| 后天国语完整版免费观看| 午夜久久久久精精品| 免费无遮挡裸体视频| 美女大奶头视频| 黄色视频,在线免费观看| 给我免费播放毛片高清在线观看| 午夜成年电影在线免费观看| 熟妇人妻久久中文字幕3abv| 波多野结衣高清无吗| 亚洲熟妇熟女久久| 日韩大码丰满熟妇| svipshipincom国产片| www日本在线高清视频| 日韩欧美精品v在线| 国产精品爽爽va在线观看网站| 好男人在线观看高清免费视频| 国产又黄又爽又无遮挡在线| 精品福利观看| 欧美在线黄色| 精品国产乱码久久久久久男人| 国产成人一区二区三区免费视频网站| 欧美一级a爱片免费观看看 | 午夜福利视频1000在线观看| 亚洲欧洲精品一区二区精品久久久| 日本熟妇午夜| 此物有八面人人有两片| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清专用| 国产高清videossex| 曰老女人黄片| 亚洲国产看品久久| 黄频高清免费视频| 精品久久蜜臀av无| 女人爽到高潮嗷嗷叫在线视频| 我的老师免费观看完整版| 久久天堂一区二区三区四区| 在线观看美女被高潮喷水网站 | 久久性视频一级片| 亚洲人成伊人成综合网2020| www.熟女人妻精品国产| 亚洲av中文字字幕乱码综合| 韩国av一区二区三区四区| 久久久久久免费高清国产稀缺| 丁香欧美五月| 九九热线精品视视频播放| 制服诱惑二区| 国产69精品久久久久777片 | 动漫黄色视频在线观看| 日韩欧美三级三区| 精品日产1卡2卡| 亚洲一区二区三区色噜噜| 欧美日韩亚洲综合一区二区三区_| 国产人伦9x9x在线观看| 国产在线观看jvid| 18禁黄网站禁片午夜丰满| 欧美日韩瑟瑟在线播放| a级毛片a级免费在线| www.熟女人妻精品国产| 国产精品久久久av美女十八| 欧美日韩福利视频一区二区| 国产视频一区二区在线看| 淫秽高清视频在线观看| 午夜精品一区二区三区免费看| 亚洲美女黄片视频| 欧美日韩福利视频一区二区| 亚洲av日韩精品久久久久久密| www国产在线视频色| 亚洲一区二区三区色噜噜| 精品免费久久久久久久清纯| 日韩欧美在线乱码| 免费高清视频大片| 色综合欧美亚洲国产小说| 一区二区三区高清视频在线| 1024视频免费在线观看| 国产成+人综合+亚洲专区| 国产精品,欧美在线| 国产真实乱freesex| 男男h啪啪无遮挡| 国产免费av片在线观看野外av| 国产成人一区二区三区免费视频网站| 麻豆成人av在线观看| 国产精品自产拍在线观看55亚洲| 啦啦啦观看免费观看视频高清| 一级a爱片免费观看的视频| 欧美色视频一区免费| 又黄又爽又免费观看的视频| 欧美日韩亚洲综合一区二区三区_| 日本一本二区三区精品| 国内揄拍国产精品人妻在线| 免费看日本二区| 一进一出抽搐gif免费好疼| 大型av网站在线播放| 黄色毛片三级朝国网站| 人人妻,人人澡人人爽秒播| av在线天堂中文字幕| 久久这里只有精品19| 久久热在线av| 国产免费男女视频| 大型黄色视频在线免费观看| 亚洲18禁久久av| 好男人在线观看高清免费视频| 久久伊人香网站| 亚洲人成网站在线播放欧美日韩| 巨乳人妻的诱惑在线观看| 午夜亚洲福利在线播放| 亚洲成a人片在线一区二区| 天堂av国产一区二区熟女人妻 | 国产免费av片在线观看野外av| 一级片免费观看大全| 欧美大码av| 国产精品亚洲一级av第二区| 一a级毛片在线观看| 国产三级黄色录像| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲美女久久久| 久久亚洲精品不卡| 国产av又大| 国产熟女xx| 亚洲午夜精品一区,二区,三区| 妹子高潮喷水视频| 精品久久久久久久毛片微露脸| 国产激情久久老熟女| 老司机深夜福利视频在线观看| 草草在线视频免费看| 99riav亚洲国产免费| 九色国产91popny在线| 国产高清视频在线播放一区| 91av网站免费观看| 级片在线观看| 亚洲精品在线美女| 国产精品精品国产色婷婷| 成人亚洲精品av一区二区| 久久香蕉精品热| 久久精品国产亚洲av香蕉五月| 欧美黄色片欧美黄色片| 免费在线观看亚洲国产| 在线a可以看的网站| 国产三级黄色录像| 日本在线视频免费播放| 身体一侧抽搐| 久久久精品欧美日韩精品| 欧美丝袜亚洲另类 | 国产精品久久视频播放| 久久久国产成人精品二区| 国产成+人综合+亚洲专区| 精品久久久久久,| 国产精品精品国产色婷婷| 欧美丝袜亚洲另类 | 国产精品一区二区免费欧美| 亚洲一区二区三区色噜噜| 天堂动漫精品| 亚洲成人久久爱视频| 一区二区三区国产精品乱码| 此物有八面人人有两片| 狠狠狠狠99中文字幕| 淫妇啪啪啪对白视频| 亚洲精品久久成人aⅴ小说| 又爽又黄无遮挡网站| 好男人电影高清在线观看| 成年人黄色毛片网站| 在线观看66精品国产| 亚洲aⅴ乱码一区二区在线播放 | 亚洲欧美日韩高清专用| 国产99白浆流出| www.www免费av| 午夜激情av网站| 国产高清有码在线观看视频 | 哪里可以看免费的av片| 久久国产精品人妻蜜桃| 美女高潮喷水抽搐中文字幕| 国产精品一区二区免费欧美| 五月玫瑰六月丁香| 久久久精品大字幕| 免费高清视频大片| 国产激情欧美一区二区| 亚洲精品美女久久久久99蜜臀| 国产片内射在线| 男女下面进入的视频免费午夜| 狂野欧美激情性xxxx| 国产av在哪里看| 国产成人av教育| 午夜视频精品福利| 国产av又大| 亚洲 欧美 日韩 在线 免费| 香蕉国产在线看| 亚洲国产高清在线一区二区三| 免费看十八禁软件| 亚洲中文字幕日韩| 成人手机av| 亚洲中文字幕日韩| 香蕉国产在线看| 日本黄色视频三级网站网址| 欧美成狂野欧美在线观看| 女同久久另类99精品国产91| 男插女下体视频免费在线播放| 亚洲专区国产一区二区| 高清毛片免费观看视频网站| 亚洲电影在线观看av| 99riav亚洲国产免费| 在线a可以看的网站| av超薄肉色丝袜交足视频| 色老头精品视频在线观看| 亚洲七黄色美女视频| 特级一级黄色大片| 国产成人精品无人区| 在线十欧美十亚洲十日本专区| 丝袜美腿诱惑在线| 亚洲最大成人中文| 亚洲国产欧美一区二区综合| 国产亚洲av高清不卡| 又大又爽又粗| 亚洲中文字幕一区二区三区有码在线看 | 久久久久九九精品影院| 全区人妻精品视频| 夜夜躁狠狠躁天天躁| 99精品欧美一区二区三区四区| 观看免费一级毛片| 国产av不卡久久| 看免费av毛片| 国产精品九九99| 日日摸夜夜添夜夜添小说| 国产黄a三级三级三级人| 欧美最黄视频在线播放免费| 九色成人免费人妻av| 国产探花在线观看一区二区| 九九热线精品视视频播放| 亚洲国产欧美网| 亚洲性夜色夜夜综合| 中文资源天堂在线| 国产黄色小视频在线观看| 欧美 亚洲 国产 日韩一| 日韩欧美免费精品| 亚洲午夜精品一区,二区,三区| 色播亚洲综合网| 婷婷精品国产亚洲av| 最近最新免费中文字幕在线| 夜夜看夜夜爽夜夜摸| 免费看十八禁软件| 亚洲成人免费电影在线观看| 国产欧美日韩一区二区精品| www.自偷自拍.com| 国产成人欧美在线观看| 国产主播在线观看一区二区| 制服人妻中文乱码| 日韩三级视频一区二区三区| netflix在线观看网站| 俄罗斯特黄特色一大片| 国产精品国产高清国产av| 69av精品久久久久久| av免费在线观看网站| 日韩高清综合在线| 老司机靠b影院| 在线观看日韩欧美| 日本五十路高清| 狂野欧美激情性xxxx| 午夜福利欧美成人| 亚洲欧洲精品一区二区精品久久久| 久久中文看片网| 日韩大尺度精品在线看网址| 90打野战视频偷拍视频| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产亚洲av高清一级| 老熟妇乱子伦视频在线观看| 桃红色精品国产亚洲av| 日本免费一区二区三区高清不卡| 久久久久久亚洲精品国产蜜桃av| avwww免费| 亚洲欧洲精品一区二区精品久久久| 美女黄网站色视频| 真人做人爱边吃奶动态| 一个人免费在线观看电影 | 国产精品av视频在线免费观看| 正在播放国产对白刺激| www.999成人在线观看| 亚洲精品av麻豆狂野| a级毛片a级免费在线| 久久精品亚洲精品国产色婷小说| 亚洲国产欧美一区二区综合| 色哟哟哟哟哟哟| 亚洲成人久久爱视频| 亚洲精品一区av在线观看| 久久草成人影院| av欧美777| 国产aⅴ精品一区二区三区波| 一二三四社区在线视频社区8| 男男h啪啪无遮挡| 叶爱在线成人免费视频播放| 亚洲一区中文字幕在线| 久久性视频一级片| 97碰自拍视频| 成人亚洲精品av一区二区| 久久天堂一区二区三区四区| 日日摸夜夜添夜夜添小说| 国产精品av久久久久免费| 天天一区二区日本电影三级| 免费看十八禁软件| 日韩免费av在线播放| 久久久久性生活片| 日韩精品青青久久久久久| 久久久国产成人免费| 国产成人aa在线观看| 欧美精品啪啪一区二区三区| 久久精品综合一区二区三区| 国产精品 欧美亚洲| 精品人妻1区二区| 日本一区二区免费在线视频| 好看av亚洲va欧美ⅴa在| 亚洲av日韩精品久久久久久密| 中文在线观看免费www的网站 | 亚洲一卡2卡3卡4卡5卡精品中文| 正在播放国产对白刺激| 999久久久精品免费观看国产| 久久精品91无色码中文字幕| 欧美高清成人免费视频www| 亚洲18禁久久av| 午夜久久久久精精品| 国产主播在线观看一区二区| 可以免费在线观看a视频的电影网站| 此物有八面人人有两片| 99精品久久久久人妻精品| 精品国产超薄肉色丝袜足j| 日韩高清综合在线| 国产精品98久久久久久宅男小说| 久久精品国产亚洲av高清一级| 中文资源天堂在线| 国产精品爽爽va在线观看网站| 国产激情欧美一区二区| 久久精品影院6| 欧美一级a爱片免费观看看 | 亚洲人成网站高清观看| 久久天堂一区二区三区四区| 成人18禁在线播放| 久久婷婷人人爽人人干人人爱| 99久久综合精品五月天人人| 亚洲自拍偷在线| 桃色一区二区三区在线观看| 人人妻,人人澡人人爽秒播| 无人区码免费观看不卡| 久久久久久九九精品二区国产 | 午夜福利高清视频| 久9热在线精品视频| 人妻夜夜爽99麻豆av| 久久性视频一级片| 午夜两性在线视频| 国产成人影院久久av| 亚洲欧美精品综合一区二区三区| 欧美在线一区亚洲| 日韩国内少妇激情av| 国产成人精品无人区| 真人做人爱边吃奶动态| 国产一区二区在线av高清观看| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 久久热在线av| 久久久久久久久久黄片| 两个人看的免费小视频| 丰满的人妻完整版| 正在播放国产对白刺激| 午夜免费观看网址| 成人18禁高潮啪啪吃奶动态图| 一区福利在线观看| 国产久久久一区二区三区| 亚洲国产看品久久| 校园春色视频在线观看| 男女下面进入的视频免费午夜| 欧美另类亚洲清纯唯美| 在线观看一区二区三区| www.精华液| 国产av一区二区精品久久| 日本在线视频免费播放| 丰满的人妻完整版| 国产又黄又爽又无遮挡在线| 久久中文看片网| 91老司机精品| 长腿黑丝高跟| 国产激情欧美一区二区| 国产亚洲欧美在线一区二区| 国产av麻豆久久久久久久| 国产精品爽爽va在线观看网站| 欧美绝顶高潮抽搐喷水| 亚洲av熟女| 国产精品98久久久久久宅男小说| 久久精品成人免费网站| 手机成人av网站| 一边摸一边抽搐一进一小说| 欧美成人午夜精品| 午夜免费激情av| 999久久久精品免费观看国产| 毛片女人毛片| 两个人的视频大全免费| 日本一区二区免费在线视频| 午夜福利在线在线| 免费电影在线观看免费观看| 欧美久久黑人一区二区| 一二三四社区在线视频社区8| 丁香六月欧美| 国产精品精品国产色婷婷| 亚洲va日本ⅴa欧美va伊人久久| 69av精品久久久久久| 久热爱精品视频在线9| 两性夫妻黄色片| 18禁观看日本| 免费在线观看成人毛片| 亚洲乱码一区二区免费版| 国产69精品久久久久777片 | 不卡一级毛片| 1024手机看黄色片| 97碰自拍视频| 成人亚洲精品av一区二区| 久久亚洲精品不卡| 欧美丝袜亚洲另类 | 国产精品久久久av美女十八| 国产精品美女特级片免费视频播放器 | 黄色视频,在线免费观看| 国产aⅴ精品一区二区三区波| 天堂√8在线中文| 人人妻人人看人人澡| 亚洲av美国av| 可以在线观看毛片的网站| 国产精品九九99| 免费一级毛片在线播放高清视频| 视频区欧美日本亚洲| 亚洲专区字幕在线| cao死你这个sao货| 久久香蕉精品热| 亚洲人成电影免费在线| 国产精品久久久久久亚洲av鲁大| 夜夜夜夜夜久久久久| 国产真人三级小视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品色激情综合| 在线观看舔阴道视频| 欧美成人性av电影在线观看| 在线a可以看的网站| 香蕉国产在线看| 久久欧美精品欧美久久欧美| 国内精品久久久久久久电影| 搡老岳熟女国产| 伦理电影免费视频| 哪里可以看免费的av片| 韩国av一区二区三区四区| 色播亚洲综合网| 欧美日韩黄片免| 桃色一区二区三区在线观看| 久久精品夜夜夜夜夜久久蜜豆 | 欧美激情久久久久久爽电影| 精品欧美国产一区二区三| 国产不卡一卡二| 亚洲av片天天在线观看| 两性夫妻黄色片| 亚洲片人在线观看| 中文资源天堂在线| 欧美人与性动交α欧美精品济南到| av欧美777| 国产亚洲精品一区二区www| 久久香蕉激情| 午夜免费观看网址| 在线十欧美十亚洲十日本专区| 三级毛片av免费| 欧美黄色片欧美黄色片| 激情在线观看视频在线高清| 可以在线观看的亚洲视频| 亚洲成人国产一区在线观看| 日韩国内少妇激情av| 国产av又大| 亚洲精品久久成人aⅴ小说| 欧美日本视频| 午夜精品久久久久久毛片777| 草草在线视频免费看| 91九色精品人成在线观看|