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

    基于注意力區(qū)域不同組分特征的磁共振成像前列腺癌包膜侵犯診斷研究

    2021-03-06 05:53:12張一鴻侯瑩包婕王成龍宋陽張玉東楊光
    磁共振成像 2021年12期
    關鍵詞:前列腺癌特征模型

    張一鴻,侯瑩,包婕,王成龍,宋陽*,張玉東,楊光

    作者單位:1.華東師范大學上海市磁共振重點實驗室,上海200062;2.南京醫(yī)科大學第一附屬醫(yī)院,南京210029;3.蘇州大學第一附屬醫(yī)院,蘇州215006

    前列腺癌是美國男性發(fā)病率第一,病死率第二的癌癥,其準確診斷是制訂最佳治療計劃的保障[1,2]。前列腺癌發(fā)生包膜侵犯與根治性前列腺切除術后總生存率下降相關[3]。傳統(tǒng)經(jīng)直腸超聲引導的活檢往往低估了腫瘤的分期,導致預測不準確[4-6]。研究表明,多參數(shù)磁共振成像對包膜侵犯的檢測是有效的[7-9]。然而,傳統(tǒng)基于磁共振成像的研究,也不能準確地判斷侵犯發(fā)生,特別是對處于頂部的癌灶,特異度只有30%左右[10,11]。

    影像組學[12]是一種通過分析圖像ROI,提取定量特征,并對特征進行統(tǒng)計分析的方法,近幾年被廣泛運用于醫(yī)學影像研究之中[13,14]。在前列腺癌包膜侵犯方面,主要是針對癌灶內(nèi)的特征進行分析[15,16]。但包膜侵犯的發(fā)生,與癌灶和腺體包膜的位置有關,僅研究一個區(qū)域不能反映兩者關系。所以本文通過提取可能侵犯區(qū)域的組學特征,建立組學模型對前列腺癌是否發(fā)生包膜侵犯進行預測。相較于單一組織特征的模型,本文所建立的模型能獲得更好的診斷效果。

    1 材料和方法

    1.1 一般資料

    本研究回顧性分析了2015 年1 月至2019 年6 月在南京醫(yī)科大學第一附屬醫(yī)院進行磁共振掃描的病例,并通過醫(yī)院的倫理委員會批準(批準文號:2016-SRFA-093),免除受試者知情同意。數(shù)據(jù)納入標準:(1)前列腺癌癥患者進行了根治性前列腺切除術;(2)術前4周內(nèi)進行標準前列腺MRI檢查;(3)磁共振檢查圖像包括T2加權成像、彌散加權成像和對應的ADC圖像。排除標準:(1)既往有前列腺癌手術史;(2)進行過前列腺癌輔助治療(良性前列腺增生或膀胱外梗阻的干預被認為是可接受的)。最后納入718 例患者,其中未發(fā)生侵犯的527 例,發(fā)生侵犯的191 例。按照訓練集和測試集比例4∶1隨機對數(shù)據(jù)拆分,得到訓練集數(shù)據(jù)574例,測試集144例。

    1.2 圖像采集

    所有數(shù)據(jù)均使用德國西門子公司的3.0 T 磁共振掃描儀Skyra。掃描參數(shù):T2WI:TE=105 ms,TR=6000 ms,層厚=3.5 mm,掃描矩陣=384×384;DWI:TE=82 ms,TR=6000 ms,層厚=3.5 mm,掃描矩陣=128×128,含有 0、100、250、500、1000、1500、2000 mm2/s等多個b值。ADC由DWI單指數(shù)擬合模型計算得到。

    1.3 圖像分割

    兩位經(jīng)驗豐富的影像科醫(yī)師使用3D Slicer,對前列腺腺體和癌灶部分進行手動勾畫,得到腺體和癌灶的ROI。由于包膜侵犯是發(fā)生在癌灶與腺體包膜接觸的區(qū)域,本研究基于腺體ROI和癌灶ROI,根據(jù)以下規(guī)則生成注意力ROI,反映可能發(fā)生包膜侵犯的區(qū)域:(1)對于在癌灶內(nèi)部,且在前列腺腺體以外的體素,賦予最大注意力值1.0;(2)對于在癌灶以內(nèi),腺體以內(nèi)的體素,根據(jù)體素到腺體包膜的距離r,賦予注意力α/r 的注意力值,其中α 是調(diào)節(jié)參數(shù),本實驗中選取0.5,即圖像分辨率;(3)對上述兩步生成的區(qū)域往外做衰減直到像素值為0。

    1.4 組學特征提取及模型建立

    使用荷蘭鹿特丹伊拉斯謨醫(yī)療中心研發(fā)的Elastix 工具包[17]將ADC 圖像配準至T2 圖像上,進行重采樣操作,使得層內(nèi)分辨率為0.5 mm×0.5 mm。對生成的注意力圖,使用閾值從0.1 至0.5 以0.1 為間隔進行二值化操作,當閾值大于0.5 之后,注意力ROI 的范圍太小,不利于特征提取。由于包膜侵犯發(fā)生在癌灶與包膜接觸的位置,可以重點分析腺體包膜內(nèi)外差異,受影像組學中微環(huán)境(Habitats)[18]啟發(fā),本研究將二值化后的注意力ROI劃分為前列腺腺體區(qū)域(prostate,Pro)、前列腺癌區(qū)域(prostate cancer,PCa)和周圍組織背景區(qū)域(background,BG)(圖1)。分別對上述ROI 提取形狀特征14 個、T2 圖像和ADC圖像的灰度特征各18個。

    圖1 實驗流程圖Fig.1 The workflow of out study

    由于數(shù)據(jù)正負樣本比例約為1∶3,本研究采用合成少數(shù)過采集術(synthetic minority oversampling technique,SMOTE)算法平衡樣本。使用Z-Score 和Mean標準化的方法進行特征歸一化,皮爾森相關系數(shù)(Pearson correlation coefficients, PCC)對數(shù)據(jù)進行降維,特征遞歸消除(recursive feature elimination,RFE)和KW (Kruskal-Wallis test)的方法對降維后的特征進行進一步選擇,除去對結果貢獻小的特征。使用支持向量機(support vector machine,SVM)、線性判別分析(linear discriminant analysis,LDA)和線性回歸(linear regression,LR)作為分類器。運用5 折交叉驗證選取最優(yōu)模型參數(shù)。使用單位標準差(one-standard error,1-SE)方法進行模型選擇,模型用獨立測試集評估。

    以上特征提取和建模操作,使用開源軟件FeAture Explorer (FAE) v 0.3.6進行操作[19]。

    1.5 模型評價

    本文使用ROC曲線和AUC來衡量模型的預測結果,并用Bootstrap 方法計算95%置信區(qū)間(confidence interval,CI)。同時還計算了敏感度、特異度、陽性率和陰性率,并使用決策曲線分析(decision curve analysis,DCA)比較模型之間收益差異。

    1.6 統(tǒng)計學分析

    統(tǒng)計學分析使用Python 3.6和Scipy 1.3,年齡和前列腺特異性抗原(prostate specific antigen,PSA)使用曼-惠特尼U檢驗,病理Gleason 評分(biology Gleason Score,bGS)和病理侵犯統(tǒng)計數(shù)量及占比,使用χ2檢驗,來計算訓練集和測試集之間的統(tǒng)計分布。模型預測結果之間的差異使用威爾科克森符號秩檢驗。P<0.05表示差異具有統(tǒng)計學意義。

    2 結果

    2.1 數(shù)據(jù)描述

    718 例患者的訓練集和測試集數(shù)據(jù)的詳細臨床信息見表1,訓練集數(shù)據(jù)和測試集數(shù)據(jù)臨床信息差異無統(tǒng)計學意義。

    表1 訓練集和測試集患者臨床信息Tab.1 Clinical information of patients in the training and test cohorts

    2.2 模型結果

    2.2.1 閾值選擇

    使用不同閾值對注意力ROI進行二值化,用交叉驗證的方式進行建模比較,橫向?qū)Ρ乳撝祵τ柧毤Y果的影響,發(fā)現(xiàn)訓練集的AUC表現(xiàn)類似,見圖2。根據(jù)奧卡姆剃刀原則[20],選取特征數(shù)最少對應的模型,即閾值選取0.2。

    圖2 不同閾值注意力ROⅠ的建模結果。A:不同閾值模型的ROC曲線;B:不同模型AUC的箱型圖和建模使用的特征數(shù)量,其中箱線圖表示模型的AUC統(tǒng)計(左軸),綠色折線表示特征數(shù)(右軸) 圖3 測試集結果。圖A展示了測試集在腺體ROⅠ、癌灶ROⅠ、注意力ROⅠ和注意力ROⅠ子區(qū)域四個模型的ROC曲線;圖B為這四個模型的DCA曲線Fig. 2 Model results built by attention ROⅠwith different thresholds. A:The ROC curve of different threshold models;B:The box chart of AUC of different models(left axis),and the green plot of the number of selected features(right axis).Fig.3 Model performance on the test cohort.A denotes the ROC curves of the ModelPro,ModelPCa,ModelAtten,ModelRegion on the test cohort.B denotes the DCA curves of the four models.

    2.2.2 模型結果

    本研究同時使用癌灶ROI、腺體ROI、閾值0.2 的注意力ROI 和注意力ROI 子區(qū)域進行建模,結果見表2。

    表2 測試集在四個模型上的結果Tab.2 The results of the test cohorts on four models

    當使用腺體ROI建模時,使用交叉驗證和1-SE的方法進行模型確立,結果發(fā)現(xiàn)使用Mean標準化方法、KW 特征選擇方法和SVM 分類器時,選出16 個特征,模型在驗證集上的效果較好(記為Mean-KW-16-SVM),此時訓練集的AUC為0.740 (95%CI:0.690~0.788)、敏感度為0.603,特異度為0.787;測試集AUC 達到0.746 (95%CI:0.651~0.835)、敏感度為0.825,特異度為0.577。當使用癌灶ROI 時,選出的模型為Zscore-KW-12-LR,此時訓練集AUC為0.742 (95%CI:0.693~0.790),敏感度為0.702、特異度為0.709;測試集的AUC為0.755 (95%CI:0.670~0.838)、敏感度和特異度分別為0.725和0.692。

    使用注意力ROI時,Mean-RFE-2-SVM的模型得到訓練集AUC 為0.732 (95%CI:0.689~0.774),敏感度和特異度分別為0.722 和0.664,陽性率和陰性率分別為0.434和0.870;測試集AUC為0.766 (95%CI:0.684~0.846),敏感度可以達到0.900,特異度0.529,陽性率和陰性率分別為0.424和0.932。將注意力ROI中不同組分的特征分開提取并組合建模,模型Mean-RFE-14-SVM 在訓練集上的AUC 為0.794(95%CI:0.753~0.835),測試集上AUC 為0.792(95%CI:0.716~0.862)。訓練集上敏感度特異度分別為0.695 和0.785,陽性率和陰性率為0.536 和0.878。測試集上敏感度、特異度、陽性率和陰性率分別為0.800、0.721、0.525和0.904。

    四個模型在測試集上的ROC 曲線和DCA 曲線展示在圖3 中。子區(qū)域模型的AUC 最高,在概率閾值取0.2 至0.4 范圍內(nèi),子區(qū)域模型獲得的凈收益高于其他三個模型。

    2.3 結果分析

    常規(guī)使用腺體ROI進行組學分析時,模型挑出的特征包含了T2 和ADC 的一階特征,但并沒有形狀特征。對癌灶ROI進行特征提取建模,主要是形狀特征進行了貢獻,即通過分析癌灶形狀,可以輔助包膜侵犯診斷,特別是2D 最大直徑,在建模中權重占據(jù)最大。在使用注意力ROI進行建模時,只使用形狀特征的體積表面積比和ADC 的灰度最大值建模AUC 就可以達到0.766。注意力圖子區(qū)域建模,測試集AUC 達到了0.792,背景、腺體、癌灶的特征均有所貢獻。具體特征貢獻圖見圖4。

    圖4 四個模型建模使用的特征及權重,其中A~D分別代表腺體模型、癌灶模型、注意力模型、子區(qū)域模型Fig. 4 The weights of features contributed in the model built by prostate gland,prostate cancer,generated attention,and the sub-regions of the attention were shown in figure A to D,respectively.

    3 討論

    3.1 總結與發(fā)現(xiàn)

    本研究通過前列腺腺體ROI 和癌灶ROI 生成注意力ROI,并受影像組學微環(huán)境的啟發(fā),將二值化后的注意力ROI分成背景、腺體和癌灶三個子ROI,分開提取子ROI 的組學特征,經(jīng)過SMOTE 樣本均衡、Mean特征歸一化、PCC 降維、RFE 特征選擇后,使用SVM 建立模型,在測試集上得到AUC為0.792,說明此方法可以較好地預測前列腺癌患者是否發(fā)生包膜侵犯。

    3.2 相關研究比較

    相較于單純使用臨床信息,磁共振成像具有可以顯示癌灶位置和解剖學信息的優(yōu)勢[21]。除了跟臨床特征結合預測外,磁共振影像可以進行特征提取,并建立組學模型對前列腺癌包膜侵犯進行預測。Cuocolo等[22]對三個中心的前列腺癌患者的T2和ADC圖像癌灶內(nèi)部特征進行提取,使用單一中心數(shù)據(jù)訓練,另外兩個中心數(shù)據(jù)進行測試,測試集AUC 分別為0.80 和0.73,模型表現(xiàn)與放射科醫(yī)生差異無統(tǒng)計學意義。Xu 等[23]在結構像、彌散像和動態(tài)增強圖像上分析癌灶特征,在115 個數(shù)據(jù)上建立LASSO (least absolute shrinkage and selection operator)模型,AUC達到0.865,但使用需要注入對比劑的動態(tài)對比增強序列,增加了掃描成本。本研究僅使用了T2 和ADC 的序列,橫向比較了多種建模方法,最高可以在獨立測試集上達到0.792 的AUC,說明了基于磁共振成像的組學模型可以對包膜侵犯進行預測。

    3.3 注意力圖ROI的分析比較

    前列腺癌包膜侵犯問題,對于ROI的勾畫至關重要。Ma 等[15]對200 多例前列腺癌患者的T2 圖像的前列腺輪廓及輪廓內(nèi)約1 mm處進行ROI的勾畫,并納入被懷疑是惡性或者有侵犯趨勢的癌灶部分,使得ROI內(nèi)部能更多地提取包膜的整體細微結構特征,建立的模型在測試集上用LASSO 模型達到0.883 的AUC。Bai 等[16]對284 例的T2 和ADC 圖像進行癌灶勾畫,將癌灶區(qū)域進行擴張并減去其中的直腸、尿道、膀胱得到癌灶周邊區(qū)域,對其進行特征提取。癌灶區(qū)域、周邊區(qū)域和臨床特征分別和組合建模,癌灶周邊區(qū)域和臨床信息的結果最好,測試集上AUC 達到0.718。但這個過程是需要臨床醫(yī)生進行手動勾畫感興趣區(qū)域,勾畫結果依賴于醫(yī)生的臨床經(jīng)驗。本研究使用的是自動生成注意力圖的方式,借鑒了微環(huán)境的實現(xiàn)思路,將注意力圖劃分為更具有組織結構信息的背景、癌灶、腺體等三個區(qū)域,分別進行特征提取建模。

    3.4 組學特征分析

    相比于其他學者影像組學的工作,本研究使用的特征主要是形狀特征和灰階特征,其原因是考慮到彌散序列的低分辨率和低信噪比,相比于紋理特征,形狀和灰階特征的魯棒性更好[24,25]。本研究的注意力圖模型驗證發(fā)現(xiàn),不管是基于整體注意力ROI還是基于其子區(qū)域,最終結果表明體積表面積比和ADC圖像各組分的灰度值對前列腺癌的包膜侵犯的預測有較大的作用。原因可能是在注意力ROI生成時,選擇的是癌灶與腺體邊界相近的部分,而腺體邊界通常類似于球形,如果癌灶沒有突破包膜,那么生成的注意力ROI很可能是月牙狀,體積表面積比較低。而各個組分在ADC圖像上的整體灰度是有差異的,比如膀胱的亮度最高,其次是腺體,然后是癌灶,背景部分的亮度最低,因此當ADC圖像組分越多,灰度差異越大,組分越復雜,越可能發(fā)生侵犯。特別是針對子區(qū)域模型,特征來源清晰,解釋性較強,在臨床上具有較好的接受程度。

    3.5 局限性

    雖然研究使用自動化的方法生成注意力ROI,但是需要醫(yī)生手工對前列腺腺體和癌灶進行勾畫,勾畫過程也是異常煩瑣。后續(xù)可以使用深度學習對腺體和癌灶進行自動的分割,使得過程自動化。其次,本研究僅僅使用了單一中心的數(shù)據(jù),不具有普適性,后續(xù)可以在更多不同設備不同中心的數(shù)據(jù)上進行驗證和優(yōu)化。

    綜上所述,本研究提出了一種自動生成注意力圖的方法,對注意力ROI進行子區(qū)域劃分并提取特征進行影像組學建模,以輔助診斷磁共振前列腺癌包膜侵犯的診斷問題。本研究所提出的模型使用了解釋性較好的特征,有助于輔助臨床醫(yī)生進行診斷,在臨床決策中具有一定的參考價值,注意力圖的思路也能夠為其他相關診斷問題提供研究思路。

    作者利益沖突聲明:全體作者均聲明無利益沖突。

    猜你喜歡
    前列腺癌特征模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    前列腺癌復發(fā)和轉(zhuǎn)移的治療
    關注前列腺癌
    認識前列腺癌
    如何表達“特征”
    不忠誠的四個特征
    當代陜西(2019年10期)2019-06-03 10:12:04
    前列腺癌,這些蛛絲馬跡要重視
    抓住特征巧觀察
    制服人妻中文乱码| 欧美日韩亚洲高清精品| 99国产精品免费福利视频| 午夜福利网站1000一区二区三区| 色吧在线观看| 精品久久久精品久久久| 飞空精品影院首页| 免费看不卡的av| 亚洲国产成人一精品久久久| 成人午夜精彩视频在线观看| 久久免费观看电影| 一级a做视频免费观看| 王馨瑶露胸无遮挡在线观看| av线在线观看网站| 日韩熟女老妇一区二区性免费视频| 国产精品偷伦视频观看了| 亚洲国产欧美在线一区| 美女cb高潮喷水在线观看| 精品99又大又爽又粗少妇毛片| 熟女av电影| 欧美精品人与动牲交sv欧美| 色哟哟·www| 国产精品久久久久成人av| 99久久人妻综合| 欧美日本中文国产一区发布| 国产一区二区在线观看日韩| 日本爱情动作片www.在线观看| 少妇人妻精品综合一区二区| 蜜桃在线观看..| 国产免费一级a男人的天堂| 国产成人91sexporn| av国产精品久久久久影院| 狂野欧美白嫩少妇大欣赏| 亚洲三级黄色毛片| 精品人妻偷拍中文字幕| 制服人妻中文乱码| 国产精品久久久久久精品古装| 精品人妻在线不人妻| 大话2 男鬼变身卡| √禁漫天堂资源中文www| 在现免费观看毛片| 黄色欧美视频在线观看| 久久免费观看电影| 丰满迷人的少妇在线观看| 水蜜桃什么品种好| 午夜免费男女啪啪视频观看| 日韩人妻高清精品专区| 91久久精品国产一区二区三区| 色婷婷久久久亚洲欧美| 80岁老熟妇乱子伦牲交| 亚洲精品视频女| 国产亚洲午夜精品一区二区久久| 男人爽女人下面视频在线观看| 久久久国产欧美日韩av| 国产女主播在线喷水免费视频网站| 成人国语在线视频| 丰满少妇做爰视频| 日韩一区二区三区影片| 免费黄网站久久成人精品| 国产午夜精品一二区理论片| 久久鲁丝午夜福利片| 色网站视频免费| 国产免费福利视频在线观看| 色网站视频免费| 国产精品一区二区在线观看99| 伦理电影大哥的女人| 日韩熟女老妇一区二区性免费视频| 免费大片黄手机在线观看| 看免费成人av毛片| 久久女婷五月综合色啪小说| 大陆偷拍与自拍| 最新的欧美精品一区二区| 超色免费av| 国产精品嫩草影院av在线观看| av不卡在线播放| 国产老妇伦熟女老妇高清| 97在线人人人人妻| 五月开心婷婷网| 精品国产露脸久久av麻豆| 久久久久久久久久人人人人人人| 国产精品国产三级专区第一集| 国产成人91sexporn| 亚洲综合色网址| 91精品伊人久久大香线蕉| 简卡轻食公司| 亚洲欧美色中文字幕在线| 两个人免费观看高清视频| 国产免费福利视频在线观看| 我的女老师完整版在线观看| 久久久精品94久久精品| 青青草视频在线视频观看| 最近的中文字幕免费完整| 十八禁高潮呻吟视频| 欧美一级a爱片免费观看看| 欧美老熟妇乱子伦牲交| 五月开心婷婷网| 夫妻性生交免费视频一级片| 美女大奶头黄色视频| 91午夜精品亚洲一区二区三区| 久久99蜜桃精品久久| 精品久久久久久电影网| 菩萨蛮人人尽说江南好唐韦庄| 99久久中文字幕三级久久日本| 菩萨蛮人人尽说江南好唐韦庄| 2018国产大陆天天弄谢| 日韩成人av中文字幕在线观看| a级毛片黄视频| 蜜桃久久精品国产亚洲av| 色94色欧美一区二区| 日韩制服骚丝袜av| 国产 精品1| 国产又色又爽无遮挡免| 男人操女人黄网站| 午夜91福利影院| 久久久久久久亚洲中文字幕| 欧美日韩成人在线一区二区| 亚洲高清免费不卡视频| 男人添女人高潮全过程视频| 如日韩欧美国产精品一区二区三区 | 只有这里有精品99| 国精品久久久久久国模美| 国产深夜福利视频在线观看| 天天躁夜夜躁狠狠久久av| 国产日韩一区二区三区精品不卡 | 人妻夜夜爽99麻豆av| 国产精品欧美亚洲77777| 久久久精品区二区三区| 色网站视频免费| xxx大片免费视频| 亚洲美女搞黄在线观看| 永久网站在线| 日韩强制内射视频| 18禁观看日本| 伊人亚洲综合成人网| 丰满迷人的少妇在线观看| 亚洲精品国产av蜜桃| 大陆偷拍与自拍| 亚洲综合色网址| 国产日韩欧美视频二区| 黑丝袜美女国产一区| 大码成人一级视频| 亚洲精品中文字幕在线视频| 一级毛片黄色毛片免费观看视频| 亚洲情色 制服丝袜| 免费黄网站久久成人精品| 狠狠婷婷综合久久久久久88av| 国产亚洲欧美精品永久| 日本黄色片子视频| .国产精品久久| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产亚洲网站| 精品国产一区二区三区久久久樱花| 国产成人精品婷婷| 纯流量卡能插随身wifi吗| 免费高清在线观看日韩| 欧美另类一区| 999精品在线视频| 亚洲成色77777| 国产精品人妻久久久影院| 久久久久国产精品人妻一区二区| 色94色欧美一区二区| 亚洲综合色网址| 久久久精品区二区三区| 天堂中文最新版在线下载| 国产精品免费大片| 久久久久久人妻| 亚洲第一区二区三区不卡| 国产片内射在线| 香蕉精品网在线| 99国产综合亚洲精品| 国产男女内射视频| 自线自在国产av| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| av在线老鸭窝| 一区在线观看完整版| 能在线免费看毛片的网站| 久久久精品免费免费高清| 亚洲精品亚洲一区二区| 国产视频首页在线观看| av国产久精品久网站免费入址| 日本色播在线视频| 久久人人爽av亚洲精品天堂| 国产乱来视频区| 国产在线一区二区三区精| 夫妻性生交免费视频一级片| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲精品久久久com| 新久久久久国产一级毛片| 夫妻性生交免费视频一级片| 各种免费的搞黄视频| 国产亚洲精品久久久com| 亚洲第一av免费看| 欧美精品人与动牲交sv欧美| 一级毛片电影观看| 水蜜桃什么品种好| 黑人巨大精品欧美一区二区蜜桃 | 国产av一区二区精品久久| 欧美日韩成人在线一区二区| 一级二级三级毛片免费看| 在线 av 中文字幕| 久久久a久久爽久久v久久| 99国产精品免费福利视频| 日本免费在线观看一区| 哪个播放器可以免费观看大片| 在线观看免费视频网站a站| 国产男女超爽视频在线观看| 哪个播放器可以免费观看大片| 丰满乱子伦码专区| 精品久久久精品久久久| 99热这里只有精品一区| 国产伦理片在线播放av一区| 日本黄大片高清| 久久精品久久久久久噜噜老黄| 免费大片18禁| 亚洲精品一区蜜桃| 99国产综合亚洲精品| 99久久中文字幕三级久久日本| 精品国产国语对白av| 中文字幕亚洲精品专区| 亚洲经典国产精华液单| 成人国产av品久久久| 国产av国产精品国产| 亚洲不卡免费看| 99国产综合亚洲精品| 国国产精品蜜臀av免费| 国产视频内射| 91精品三级在线观看| 日韩伦理黄色片| 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品古装| 久久久久久久久久久免费av| 女人精品久久久久毛片| av国产久精品久网站免费入址| 国产成人精品一,二区| 男女免费视频国产| 国国产精品蜜臀av免费| 久热久热在线精品观看| 看十八女毛片水多多多| 在线观看三级黄色| 美女cb高潮喷水在线观看| 夜夜爽夜夜爽视频| 久久久久久久久久久久大奶| 寂寞人妻少妇视频99o| 日韩,欧美,国产一区二区三区| 亚洲熟女精品中文字幕| 国产日韩一区二区三区精品不卡 | 中文字幕最新亚洲高清| 三级国产精品片| 国产不卡av网站在线观看| 亚洲色图 男人天堂 中文字幕 | 国产一区亚洲一区在线观看| 美女cb高潮喷水在线观看| 国产精品麻豆人妻色哟哟久久| 欧美日韩一区二区视频在线观看视频在线| 亚洲人成77777在线视频| 欧美bdsm另类| 国产色爽女视频免费观看| 99视频精品全部免费 在线| 2018国产大陆天天弄谢| 国产亚洲av片在线观看秒播厂| 国产精品一国产av| 三级国产精品片| 国产成人精品久久久久久| xxxhd国产人妻xxx| 日日摸夜夜添夜夜添av毛片| 一级a做视频免费观看| 熟妇人妻不卡中文字幕| 亚洲精品日韩在线中文字幕| 2021少妇久久久久久久久久久| 久久鲁丝午夜福利片| 国产视频首页在线观看| 久久这里有精品视频免费| 国产精品秋霞免费鲁丝片| 青春草国产在线视频| 日本-黄色视频高清免费观看| 久久97久久精品| 考比视频在线观看| 一本色道久久久久久精品综合| 18禁动态无遮挡网站| 亚洲av国产av综合av卡| 亚洲欧美中文字幕日韩二区| 天堂中文最新版在线下载| 卡戴珊不雅视频在线播放| 97超碰精品成人国产| 日韩中字成人| 日韩免费高清中文字幕av| 性色avwww在线观看| 大香蕉97超碰在线| 又黄又爽又刺激的免费视频.| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产熟女欧美一区二区| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| 亚洲不卡免费看| 中文天堂在线官网| 久久99热这里只频精品6学生| 一级毛片aaaaaa免费看小| 菩萨蛮人人尽说江南好唐韦庄| 性色av一级| 日韩av不卡免费在线播放| 午夜福利影视在线免费观看| 人体艺术视频欧美日本| 午夜91福利影院| 黄色一级大片看看| 日日爽夜夜爽网站| 人妻制服诱惑在线中文字幕| 天堂8中文在线网| h视频一区二区三区| 成人国语在线视频| 高清在线视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 特大巨黑吊av在线直播| 五月伊人婷婷丁香| 嘟嘟电影网在线观看| 2022亚洲国产成人精品| 简卡轻食公司| 国产精品 国内视频| 老司机亚洲免费影院| 夫妻性生交免费视频一级片| 好男人视频免费观看在线| 岛国毛片在线播放| 91久久精品国产一区二区成人| 在线精品无人区一区二区三| 久久人人爽av亚洲精品天堂| 欧美激情国产日韩精品一区| 午夜免费鲁丝| 亚洲国产日韩一区二区| 人人澡人人妻人| 蜜桃在线观看..| 精品久久国产蜜桃| 国产av一区二区精品久久| 啦啦啦在线观看免费高清www| 丁香六月天网| 亚洲精品色激情综合| 亚洲av在线观看美女高潮| 3wmmmm亚洲av在线观看| 三上悠亚av全集在线观看| 欧美日韩视频精品一区| 欧美成人精品欧美一级黄| 大香蕉久久成人网| 国产69精品久久久久777片| 18在线观看网站| 日韩 亚洲 欧美在线| 寂寞人妻少妇视频99o| 婷婷色综合www| 精品国产一区二区三区久久久樱花| 热re99久久国产66热| 少妇高潮的动态图| 美女内射精品一级片tv| 99视频精品全部免费 在线| 精品卡一卡二卡四卡免费| 成人免费观看视频高清| 美女内射精品一级片tv| 精品国产一区二区三区久久久樱花| 毛片一级片免费看久久久久| 波野结衣二区三区在线| 丰满少妇做爰视频| 狠狠婷婷综合久久久久久88av| 精品久久久精品久久久| 欧美三级亚洲精品| 国产片特级美女逼逼视频| 人妻夜夜爽99麻豆av| 五月天丁香电影| 国产免费福利视频在线观看| a级毛片黄视频| 午夜福利网站1000一区二区三区| 熟女人妻精品中文字幕| 2022亚洲国产成人精品| 成人二区视频| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 各种免费的搞黄视频| av在线老鸭窝| 2018国产大陆天天弄谢| 另类亚洲欧美激情| 久久人人爽av亚洲精品天堂| 日韩,欧美,国产一区二区三区| av在线播放精品| 黄色一级大片看看| 精品国产国语对白av| 久久国产精品大桥未久av| 日日爽夜夜爽网站| 老司机影院毛片| 精品久久久精品久久久| 国产精品久久久久久久电影| 免费日韩欧美在线观看| 亚洲欧美一区二区三区黑人 | 高清不卡的av网站| 亚洲成色77777| videossex国产| 国产精品99久久久久久久久| 少妇的逼好多水| 啦啦啦中文免费视频观看日本| 男人爽女人下面视频在线观看| 午夜免费男女啪啪视频观看| 最近最新中文字幕免费大全7| 国产精品久久久久久精品古装| 国产高清不卡午夜福利| 日本wwww免费看| 亚洲成人一二三区av| 亚洲精品乱码久久久久久按摩| 久久久久网色| 99视频精品全部免费 在线| 国产成人a∨麻豆精品| 久久鲁丝午夜福利片| 国产成人av激情在线播放 | 内地一区二区视频在线| 亚洲精品aⅴ在线观看| 男男h啪啪无遮挡| 国产精品人妻久久久久久| 日本wwww免费看| 亚洲av综合色区一区| 亚洲欧美精品自产自拍| 制服人妻中文乱码| 亚洲av成人精品一区久久| 大香蕉久久成人网| 亚洲国产色片| 久久久久久久久久人人人人人人| 亚洲欧美成人综合另类久久久| 欧美日韩国产mv在线观看视频| 一级爰片在线观看| av网站免费在线观看视频| av一本久久久久| 天堂俺去俺来也www色官网| 国产探花极品一区二区| 纯流量卡能插随身wifi吗| 成年美女黄网站色视频大全免费 | av专区在线播放| 免费人妻精品一区二区三区视频| 久久97久久精品| 欧美成人午夜免费资源| 日本-黄色视频高清免费观看| 免费久久久久久久精品成人欧美视频 | 老司机影院毛片| 各种免费的搞黄视频| 亚洲欧美一区二区三区国产| 中文字幕制服av| 国产成人a∨麻豆精品| 精品一品国产午夜福利视频| 一本色道久久久久久精品综合| 18禁动态无遮挡网站| 久久99蜜桃精品久久| a级毛片黄视频| 国产有黄有色有爽视频| 久久午夜福利片| 尾随美女入室| 少妇的逼好多水| 只有这里有精品99| 欧美日韩av久久| 久久久国产精品麻豆| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 伊人久久精品亚洲午夜| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美中文字幕日韩二区| 99国产精品免费福利视频| 女人精品久久久久毛片| 内地一区二区视频在线| 亚洲精品aⅴ在线观看| 精品少妇久久久久久888优播| videos熟女内射| 最新的欧美精品一区二区| 精品少妇内射三级| 亚洲精品久久成人aⅴ小说 | 久久人人爽人人爽人人片va| 国产精品蜜桃在线观看| 日本黄大片高清| 在线观看美女被高潮喷水网站| videos熟女内射| 亚洲精品亚洲一区二区| 国产av国产精品国产| 精品久久久噜噜| 国产 一区精品| 美女视频免费永久观看网站| 天天操日日干夜夜撸| 国精品久久久久久国模美| 精品一区二区三卡| 国产精品成人在线| 日韩免费高清中文字幕av| 尾随美女入室| 一本色道久久久久久精品综合| 性色av一级| 22中文网久久字幕| videossex国产| 亚洲欧美成人精品一区二区| www.色视频.com| 狂野欧美白嫩少妇大欣赏| 久久精品人人爽人人爽视色| 91午夜精品亚洲一区二区三区| 中文字幕av电影在线播放| av专区在线播放| 色哟哟·www| 乱人伦中国视频| 2022亚洲国产成人精品| 一级毛片我不卡| 色5月婷婷丁香| 国产一区亚洲一区在线观看| 国产无遮挡羞羞视频在线观看| 国产熟女欧美一区二区| 国产精品女同一区二区软件| 午夜福利,免费看| 午夜激情久久久久久久| 韩国高清视频一区二区三区| 国产免费现黄频在线看| 日韩av免费高清视频| 欧美成人午夜免费资源| 国产黄色免费在线视频| 一区二区三区乱码不卡18| 欧美激情 高清一区二区三区| 久久ye,这里只有精品| 日韩视频在线欧美| 日韩中字成人| 亚洲精品第二区| 免费观看在线日韩| 国产精品免费大片| 亚洲精品国产av蜜桃| 亚洲精品日韩在线中文字幕| 中文字幕精品免费在线观看视频 | 国产成人午夜福利电影在线观看| 最近最新中文字幕免费大全7| 嘟嘟电影网在线观看| 日韩人妻高清精品专区| 午夜91福利影院| 十八禁网站网址无遮挡| 黑人猛操日本美女一级片| 国产高清三级在线| 一级毛片 在线播放| 伊人久久国产一区二区| 国产精品无大码| 国产精品久久久久久精品电影小说| 欧美最新免费一区二区三区| 国产成人av激情在线播放 | 精品熟女少妇av免费看| av国产精品久久久久影院| 18+在线观看网站| 国产高清国产精品国产三级| 日本欧美视频一区| 少妇人妻 视频| 国产精品久久久久久精品古装| 亚洲精品aⅴ在线观看| 春色校园在线视频观看| 欧美日本中文国产一区发布| 中国国产av一级| 欧美97在线视频| 女人精品久久久久毛片| 国产淫语在线视频| 午夜av观看不卡| 青春草国产在线视频| a级毛片在线看网站| 免费观看的影片在线观看| 亚洲经典国产精华液单| 丰满迷人的少妇在线观看| 久久这里有精品视频免费| 中国三级夫妇交换| 精品亚洲乱码少妇综合久久| av免费在线看不卡| 高清在线视频一区二区三区| 国产国语露脸激情在线看| 香蕉精品网在线| 亚洲天堂av无毛| 韩国高清视频一区二区三区| 男女啪啪激烈高潮av片| 国产精品久久久久久精品电影小说| 日本免费在线观看一区| 18禁在线无遮挡免费观看视频| 欧美激情极品国产一区二区三区 | 午夜免费鲁丝| 另类精品久久| 午夜福利网站1000一区二区三区| 亚洲av成人精品一二三区| 蜜桃久久精品国产亚洲av| av免费在线看不卡| 欧美精品高潮呻吟av久久| 夫妻午夜视频| 国产一区二区三区综合在线观看 | 午夜久久久在线观看| 亚洲无线观看免费| 黄色毛片三级朝国网站| 欧美三级亚洲精品| 久久精品夜色国产| 国产精品一国产av| 国产深夜福利视频在线观看| 中国三级夫妇交换| 日韩精品免费视频一区二区三区 | 国产一区二区三区av在线| 久久精品国产a三级三级三级| 国产一区亚洲一区在线观看| 青春草视频在线免费观看| 男女国产视频网站| 国产毛片在线视频| 一二三四中文在线观看免费高清| 美女国产高潮福利片在线看| 亚洲精品美女久久av网站| 国产综合精华液| 国产探花极品一区二区| 日韩成人av中文字幕在线观看| 免费看光身美女| 国产av国产精品国产| av有码第一页| 亚洲av免费高清在线观看| 亚洲av欧美aⅴ国产| a级毛片黄视频| 男女啪啪激烈高潮av片| 少妇熟女欧美另类| 国产精品成人在线| 晚上一个人看的免费电影| 午夜av观看不卡| 午夜精品国产一区二区电影| 日韩制服骚丝袜av| 日本91视频免费播放| 啦啦啦在线观看免费高清www| 欧美激情 高清一区二区三区| 国产高清有码在线观看视频| 伦理电影大哥的女人|