• <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)遵循的倫理原則
    在线观看美女被高潮喷水网站| 人人妻人人添人人爽欧美一区卜| 国产精品.久久久| 秋霞在线观看毛片| 卡戴珊不雅视频在线播放| 国产成人免费观看mmmm| 22中文网久久字幕| 另类亚洲欧美激情| 在线观看人妻少妇| 99热这里只有是精品在线观看| 9色porny在线观看| 26uuu在线亚洲综合色| 两个人的视频大全免费| 欧美最新免费一区二区三区| 在线播放无遮挡| 久久精品国产亚洲av涩爱| 亚洲,欧美,日韩| 91成人精品电影| 中文字幕制服av| 精品少妇内射三级| 午夜激情久久久久久久| 精品国产乱码久久久久久小说| 99精国产麻豆久久婷婷| 高清视频免费观看一区二区| 欧美成人午夜免费资源| 最近2019中文字幕mv第一页| 国产成人av激情在线播放 | 特大巨黑吊av在线直播| 久久青草综合色| 久久久久视频综合| 高清黄色对白视频在线免费看| 伊人久久精品亚洲午夜| 国产成人午夜福利电影在线观看| 亚洲精品久久成人aⅴ小说 | h视频一区二区三区| 国产一级毛片在线| 简卡轻食公司| 久久99热这里只频精品6学生| 国产精品人妻久久久影院| 制服诱惑二区| 久久99一区二区三区| 女性生殖器流出的白浆| 久热这里只有精品99| 精品亚洲乱码少妇综合久久| 一边摸一边做爽爽视频免费| 免费黄色在线免费观看| 一本一本综合久久| 美女内射精品一级片tv| 久久人人爽人人爽人人片va| 久久久久精品久久久久真实原创| 国产有黄有色有爽视频| 日日摸夜夜添夜夜添av毛片| 在线观看人妻少妇| 亚洲精品456在线播放app| 久久久久久久精品精品| 国产男人的电影天堂91| 亚洲精品亚洲一区二区| 亚洲av国产av综合av卡| 亚洲av二区三区四区| 欧美日本中文国产一区发布| 赤兔流量卡办理| 黄色欧美视频在线观看| 日本欧美视频一区| 久久久久久久亚洲中文字幕| 少妇人妻精品综合一区二区| 亚洲国产最新在线播放| 观看美女的网站| 国产一区亚洲一区在线观看| 免费观看a级毛片全部| 国产精品 国内视频| 九九在线视频观看精品| 精品少妇黑人巨大在线播放| 国产男女超爽视频在线观看| 亚洲av综合色区一区| a级毛片黄视频| 一本久久精品| 99九九在线精品视频| videossex国产| 国产视频首页在线观看| 夫妻午夜视频| 亚州av有码| 一级片'在线观看视频| 日韩中文字幕视频在线看片| 日韩精品有码人妻一区| 亚洲熟女精品中文字幕| 日本wwww免费看| 视频在线观看一区二区三区| 午夜视频国产福利| 免费黄网站久久成人精品| 久久人人爽人人片av| av不卡在线播放| 久久国内精品自在自线图片| av不卡在线播放| 美女cb高潮喷水在线观看| 午夜老司机福利剧场| 夜夜骑夜夜射夜夜干| 国产乱来视频区| 亚洲欧美清纯卡通| 性色avwww在线观看| 久久 成人 亚洲| 久久ye,这里只有精品| 高清av免费在线| 精品午夜福利在线看| 9色porny在线观看| 免费不卡的大黄色大毛片视频在线观看| 2018国产大陆天天弄谢| 国产高清不卡午夜福利| 老司机影院毛片| 国产 一区精品| 一本—道久久a久久精品蜜桃钙片| 国产亚洲午夜精品一区二区久久| √禁漫天堂资源中文www| 国产亚洲最大av| 午夜久久久在线观看| 中文精品一卡2卡3卡4更新| 极品人妻少妇av视频| 国产国拍精品亚洲av在线观看| 国产亚洲午夜精品一区二区久久| 亚洲一级一片aⅴ在线观看| 一本色道久久久久久精品综合| 国产亚洲欧美精品永久| 有码 亚洲区| 色婷婷av一区二区三区视频| 亚州av有码| 国产综合精华液| 观看美女的网站| 观看美女的网站| videossex国产| 国产精品久久久久成人av| 国产在视频线精品| 在线免费观看不下载黄p国产| 国产黄色视频一区二区在线观看| 狂野欧美激情性bbbbbb| 久久久久国产精品人妻一区二区| 国产成人aa在线观看| 天堂俺去俺来也www色官网| 在线观看www视频免费| 九色亚洲精品在线播放| 黄色视频在线播放观看不卡| 午夜激情福利司机影院| 国产视频内射| 中文字幕亚洲精品专区| 在线观看人妻少妇| 亚洲精品久久成人aⅴ小说 | 99久久精品一区二区三区| 日韩视频在线欧美| 夜夜爽夜夜爽视频| 日韩在线高清观看一区二区三区| 日韩成人伦理影院| 日本欧美视频一区| 亚洲少妇的诱惑av| 中文字幕久久专区| 国产不卡av网站在线观看| 九草在线视频观看| 熟女电影av网| 国产 一区精品| 久久精品国产自在天天线| 免费播放大片免费观看视频在线观看| 成人国产麻豆网| 国产成人免费观看mmmm| 777米奇影视久久| 大码成人一级视频| 视频中文字幕在线观看| 国产精品女同一区二区软件| 青青草视频在线视频观看| 亚洲色图综合在线观看| 国产视频首页在线观看| 久久午夜综合久久蜜桃| 99久久精品国产国产毛片| 一本久久精品| 久热久热在线精品观看| 免费观看a级毛片全部| 全区人妻精品视频| 欧美精品一区二区大全| 国产一区二区在线观看日韩| 亚洲国产精品专区欧美| 亚洲欧洲精品一区二区精品久久久 | 国产精品一区二区在线观看99| 九草在线视频观看| 黑丝袜美女国产一区| 熟女电影av网| 中文欧美无线码| 看非洲黑人一级黄片| 欧美日韩视频精品一区| 日韩精品有码人妻一区| 久久久精品94久久精品| 激情五月婷婷亚洲| 亚洲欧美日韩卡通动漫| 免费日韩欧美在线观看| 丰满少妇做爰视频| 欧美 亚洲 国产 日韩一| 亚洲国产精品成人久久小说| 99热6这里只有精品| 国产日韩欧美视频二区| 精品亚洲乱码少妇综合久久| 精品视频人人做人人爽| 国产白丝娇喘喷水9色精品| 又大又黄又爽视频免费| 久久99热这里只频精品6学生| 国产69精品久久久久777片| 国产精品国产三级国产专区5o| 亚洲第一av免费看| www.av在线官网国产| 91午夜精品亚洲一区二区三区| 我的女老师完整版在线观看| 亚洲成人手机| 18禁观看日本| 日韩成人伦理影院| 亚洲国产精品一区二区三区在线| 国产精品秋霞免费鲁丝片| 免费不卡的大黄色大毛片视频在线观看| 伦理电影免费视频| 亚洲四区av| 久久精品国产亚洲网站| 久久精品国产鲁丝片午夜精品| 亚洲精品久久成人aⅴ小说 | 欧美最新免费一区二区三区| 我的老师免费观看完整版| 亚洲激情五月婷婷啪啪| 又大又黄又爽视频免费| 成人18禁高潮啪啪吃奶动态图 | 亚洲av免费高清在线观看| 久热久热在线精品观看| 成年女人在线观看亚洲视频| 你懂的网址亚洲精品在线观看| 蜜臀久久99精品久久宅男| 国产极品粉嫩免费观看在线 | 9色porny在线观看| 热re99久久国产66热| 亚洲国产精品成人久久小说| 女的被弄到高潮叫床怎么办| 中文字幕久久专区| 亚洲中文av在线| 2022亚洲国产成人精品| a级毛片免费高清观看在线播放| 大片免费播放器 马上看| 亚洲国产色片| 国产精品久久久久久久久免| 青春草视频在线免费观看| 女人久久www免费人成看片| 免费人妻精品一区二区三区视频| 欧美人与性动交α欧美精品济南到 | 人人妻人人澡人人看| av天堂久久9| 国产视频内射| 精品一品国产午夜福利视频| 成人毛片60女人毛片免费| 日韩不卡一区二区三区视频在线| 高清不卡的av网站| 精品卡一卡二卡四卡免费| 青春草国产在线视频| 亚洲无线观看免费| 午夜日本视频在线| 97超视频在线观看视频| 王馨瑶露胸无遮挡在线观看| 天天操日日干夜夜撸| 我的老师免费观看完整版| 少妇被粗大猛烈的视频| 国产男人的电影天堂91| 国产精品久久久久久av不卡| 香蕉精品网在线| 成人午夜精彩视频在线观看| 交换朋友夫妻互换小说| 777米奇影视久久| 99九九在线精品视频| 国产 精品1| 少妇 在线观看| 另类精品久久| 久久狼人影院| 日韩av不卡免费在线播放| 成人黄色视频免费在线看| 午夜av观看不卡| 亚洲av福利一区| 亚洲综合精品二区| 日韩视频在线欧美| 狠狠婷婷综合久久久久久88av| 国产精品国产三级专区第一集| 97超视频在线观看视频| 精品久久久精品久久久| 伊人久久国产一区二区| 国产不卡av网站在线观看| 精品熟女少妇av免费看| 亚洲精品日韩av片在线观看| 国产精品人妻久久久影院| 久久青草综合色| 在现免费观看毛片| 自线自在国产av| 麻豆成人av视频| 国产精品 国内视频| 2018国产大陆天天弄谢| 男人爽女人下面视频在线观看| 久久久a久久爽久久v久久| 一级毛片黄色毛片免费观看视频| 中国三级夫妇交换| 亚洲欧美清纯卡通| 国产深夜福利视频在线观看| 高清午夜精品一区二区三区| 少妇人妻 视频| 丰满少妇做爰视频| 精品亚洲成国产av| 国产精品.久久久| 搡女人真爽免费视频火全软件| 18禁在线播放成人免费| 超碰97精品在线观看| 狠狠精品人妻久久久久久综合| 国产高清不卡午夜福利| 黑人猛操日本美女一级片| 亚洲图色成人| 最近中文字幕2019免费版| 少妇丰满av| 亚洲高清免费不卡视频| 乱码一卡2卡4卡精品| 日本黄大片高清| 日韩强制内射视频| 丝瓜视频免费看黄片| 大话2 男鬼变身卡| 成人18禁高潮啪啪吃奶动态图 | 亚洲五月色婷婷综合| 日韩电影二区| 精品人妻偷拍中文字幕| 国产一区亚洲一区在线观看| av女优亚洲男人天堂| 免费看不卡的av| 久久精品久久久久久噜噜老黄| 男人爽女人下面视频在线观看| 国产又色又爽无遮挡免| 免费观看av网站的网址| 亚洲第一区二区三区不卡| 99热国产这里只有精品6| 国产精品一区二区在线不卡| 一级黄片播放器| 亚洲国产欧美日韩在线播放| 中国美白少妇内射xxxbb| 看十八女毛片水多多多| 欧美精品一区二区大全| 丰满迷人的少妇在线观看| 中文乱码字字幕精品一区二区三区| 极品少妇高潮喷水抽搐| 国产男女超爽视频在线观看| 日韩制服骚丝袜av| 国产成人免费无遮挡视频| 国产一级毛片在线| 成人免费观看视频高清| 一区二区三区精品91| 男女高潮啪啪啪动态图| 国产精品熟女久久久久浪| tube8黄色片| 精品熟女少妇av免费看| 亚洲一级一片aⅴ在线观看| 婷婷色综合www| 久久国产精品男人的天堂亚洲 | 免费高清在线观看日韩| 男女高潮啪啪啪动态图| 亚洲av不卡在线观看| 久久久久久久久久成人| 国产黄频视频在线观看| 狂野欧美激情性xxxx在线观看| 五月天丁香电影| 国产亚洲av片在线观看秒播厂| 在线观看免费日韩欧美大片 | 午夜福利视频在线观看免费| 欧美日韩综合久久久久久| 少妇被粗大猛烈的视频| 在线观看免费高清a一片| 日韩精品免费视频一区二区三区 | 亚洲少妇的诱惑av| 99国产综合亚洲精品| 一本色道久久久久久精品综合| 男人爽女人下面视频在线观看| 中文字幕久久专区| 国产免费一级a男人的天堂| 国产亚洲欧美精品永久| 中文乱码字字幕精品一区二区三区| h视频一区二区三区| 免费人妻精品一区二区三区视频| 午夜免费鲁丝| 新久久久久国产一级毛片| 交换朋友夫妻互换小说| 观看美女的网站| 久久精品国产a三级三级三级| 精品国产一区二区三区久久久樱花| 精品午夜福利在线看| 建设人人有责人人尽责人人享有的| av有码第一页| 亚洲精品色激情综合| 成人二区视频| 99视频精品全部免费 在线| 啦啦啦在线观看免费高清www| 久久精品久久久久久久性| 在线播放无遮挡| 街头女战士在线观看网站| 亚洲高清免费不卡视频| 80岁老熟妇乱子伦牲交| 男女边摸边吃奶| 久久 成人 亚洲| 青春草视频在线免费观看| 国产一区二区在线观看日韩| 一二三四中文在线观看免费高清| 999精品在线视频| 夜夜看夜夜爽夜夜摸| 狂野欧美激情性xxxx在线观看| 中文字幕免费在线视频6| 高清午夜精品一区二区三区| 伦理电影免费视频| 欧美成人午夜免费资源| 天天操日日干夜夜撸| 99九九在线精品视频| 国产亚洲av片在线观看秒播厂| 国产成人精品无人区| 青青草视频在线视频观看| 日本欧美视频一区| 成人手机av| 午夜福利影视在线免费观看| 亚洲精品一二三| 日本wwww免费看| 女性被躁到高潮视频| 欧美激情极品国产一区二区三区 | 亚洲精品亚洲一区二区| xxxhd国产人妻xxx| 97在线人人人人妻| 91午夜精品亚洲一区二区三区| 国产精品久久久久久av不卡| 亚洲精品久久成人aⅴ小说 | 精品一区在线观看国产| 国产爽快片一区二区三区| 亚洲av二区三区四区| 免费久久久久久久精品成人欧美视频 | 亚洲国产日韩一区二区| freevideosex欧美| 少妇猛男粗大的猛烈进出视频| 亚洲第一av免费看| 亚洲精品,欧美精品| av国产精品久久久久影院| 精品人妻一区二区三区麻豆| 亚洲av欧美aⅴ国产| 免费观看性生交大片5| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩在线中文字幕| 国产av国产精品国产| 久热久热在线精品观看| 国产成人精品无人区| 性色avwww在线观看| 亚洲欧美成人精品一区二区| 如日韩欧美国产精品一区二区三区 | 国产精品久久久久久av不卡| videos熟女内射| 少妇被粗大猛烈的视频| 国语对白做爰xxxⅹ性视频网站| 国产伦理片在线播放av一区| 一级毛片 在线播放| 久久久久久久大尺度免费视频| 日本免费在线观看一区| 丰满迷人的少妇在线观看| 九九在线视频观看精品| 精品久久久噜噜| 久久久久精品久久久久真实原创| 午夜福利视频精品| 精品少妇黑人巨大在线播放| 日本猛色少妇xxxxx猛交久久| 中文字幕制服av| 欧美精品一区二区大全| 能在线免费看毛片的网站| 成年人免费黄色播放视频| 男女边吃奶边做爰视频| 国产精品麻豆人妻色哟哟久久| 日韩制服骚丝袜av| 女的被弄到高潮叫床怎么办| av播播在线观看一区| 国产精品国产三级国产专区5o| 精品人妻偷拍中文字幕| 99九九在线精品视频| www.av在线官网国产| av国产精品久久久久影院| 久久婷婷青草| 九色成人免费人妻av| 久久精品国产a三级三级三级| 欧美bdsm另类| 婷婷成人精品国产| 亚洲美女黄色视频免费看| 久久99热6这里只有精品| 国产极品天堂在线| 色视频在线一区二区三区| 考比视频在线观看| 国产一区二区三区av在线| 人妻系列 视频| 最新的欧美精品一区二区| 国产片特级美女逼逼视频| 久久久国产精品麻豆| 制服人妻中文乱码| 久久这里有精品视频免费| 亚洲国产欧美在线一区| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 在现免费观看毛片| 亚洲人成77777在线视频| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 多毛熟女@视频| 欧美一级a爱片免费观看看| 国产成人91sexporn| 婷婷色av中文字幕| 日韩人妻高清精品专区| 丝袜美足系列| 精品人妻熟女毛片av久久网站| 狂野欧美激情性xxxx在线观看| 黄色欧美视频在线观看| 色婷婷av一区二区三区视频| 欧美丝袜亚洲另类| 久久久午夜欧美精品| 美女国产高潮福利片在线看| 人妻夜夜爽99麻豆av| 曰老女人黄片| 日本-黄色视频高清免费观看| 精品少妇内射三级| 99国产精品免费福利视频| 中文精品一卡2卡3卡4更新| 久久久久久人妻| 亚洲精品aⅴ在线观看| 成人无遮挡网站| 国产精品一区二区三区四区免费观看| 男女边吃奶边做爰视频| av在线app专区| 久久精品久久久久久噜噜老黄| 欧美xxxx性猛交bbbb| 日韩av不卡免费在线播放| 欧美日韩国产mv在线观看视频| 九九在线视频观看精品| 一本久久精品| 高清午夜精品一区二区三区| 在现免费观看毛片| 男的添女的下面高潮视频| 水蜜桃什么品种好| 亚洲精品乱码久久久久久按摩| 热re99久久国产66热| 国产欧美亚洲国产| 在线观看美女被高潮喷水网站| 蜜桃在线观看..| 日韩av免费高清视频| 97超碰精品成人国产| 激情五月婷婷亚洲| 男女无遮挡免费网站观看| 国产亚洲午夜精品一区二区久久| tube8黄色片| 成人国产av品久久久| 高清毛片免费看| av免费观看日本| a级毛色黄片| 亚洲精品久久久久久婷婷小说| 亚洲国产色片| 中文天堂在线官网| 51国产日韩欧美| 亚洲精品亚洲一区二区| 亚洲精品国产av成人精品| 插逼视频在线观看| 插阴视频在线观看视频| 中文字幕人妻丝袜制服| 亚洲成人手机| 精品少妇黑人巨大在线播放| 啦啦啦中文免费视频观看日本| 一个人看视频在线观看www免费| 搡老乐熟女国产| 日韩一区二区视频免费看| 免费久久久久久久精品成人欧美视频 | 我的老师免费观看完整版| 春色校园在线视频观看| 国产成人精品在线电影| 午夜福利视频精品| 大香蕉久久网| 久久99热这里只频精品6学生| 日韩电影二区| 我的女老师完整版在线观看| 日日啪夜夜爽| xxx大片免费视频| 免费人成在线观看视频色| 又粗又硬又长又爽又黄的视频| 少妇被粗大的猛进出69影院 | 美女国产视频在线观看| 国产一区有黄有色的免费视频| 啦啦啦在线观看免费高清www| 热99久久久久精品小说推荐| 日日摸夜夜添夜夜爱| 国产成人精品在线电影| 人妻人人澡人人爽人人| 99视频精品全部免费 在线| 老司机亚洲免费影院| 亚洲av中文av极速乱| 飞空精品影院首页| 菩萨蛮人人尽说江南好唐韦庄| 高清在线视频一区二区三区| 黄片无遮挡物在线观看| 在线观看www视频免费| 久久久久精品久久久久真实原创| 精品亚洲成国产av| 777米奇影视久久| 精品人妻在线不人妻| 亚洲内射少妇av| 黄色配什么色好看| 中文字幕精品免费在线观看视频 | 欧美精品人与动牲交sv欧美| 妹子高潮喷水视频| 亚洲伊人久久精品综合| 国产高清三级在线| 黑人高潮一二区| 国产免费又黄又爽又色| 在线观看三级黄色| 在线播放无遮挡| 国产成人a∨麻豆精品| 大话2 男鬼变身卡| 少妇熟女欧美另类| 又大又黄又爽视频免费| 亚洲av综合色区一区| 极品人妻少妇av视频| 久久这里有精品视频免费| 91在线精品国自产拍蜜月| 久久人妻熟女aⅴ|