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

    基于分組注意力殘差網(wǎng)絡(luò)的錐束XLCT成像

    2022-07-04 07:18:40周仲興郭司琪賈夢宇
    關(guān)鍵詞:方法

    周仲興,郭司琪,賈夢宇,張 林,吳 越

    基于分組注意力殘差網(wǎng)絡(luò)的錐束XLCT成像

    周仲興1, 2,郭司琪1, 2,賈夢宇1, 2,張 林2,吳 越2

    (1. 天津大學(xué)精密儀器與光電子工程學(xué)院,天津 300072;2. 天津市生物醫(yī)學(xué)檢測技術(shù)與儀器重點實驗室,天津300072)

    錐束X射線磷光斷層成像(cone-beam X-ray luminescence computed tomography,CB-XLCT)作為一種新興的混合成像技術(shù)可同時獲得解剖結(jié)構(gòu)信息和功能代謝信息,在疾病的早期診斷和靶向治療上具有廣泛發(fā)展前景.然而,由于錐束XLCT成像過程中的高度不適定性,重建結(jié)果的空間分辨率較低.針對該問題,提出了一種基于分組注意力殘差網(wǎng)絡(luò)的錐束XLCT定量磷光分布重建方法.該方法采用注意力機(jī)制優(yōu)化光子密度特征通道的權(quán)重系數(shù),并結(jié)合近似域約束損失函數(shù),從而改善殘差網(wǎng)絡(luò)的特征表達(dá)過程,提高圖像重建精度和模型穩(wěn)健性.結(jié)果表明:所提出的方法在仿真實驗中可清晰分辨并重建出靶心距離為7mm的雙目標(biāo)體,位置誤差(LE)約為0.37mm,戴斯相似性指數(shù)(Dice)達(dá)到了84%.在仿體實驗中,LE約為0.48mm,Dice達(dá)到了79%.因此,分組注意力殘差網(wǎng)絡(luò)方法能夠有效在定位精度、空間分辨率、穩(wěn)健性等多方面提高圖像重建質(zhì)量.

    錐束XLCT;深度學(xué)習(xí);分組注意力;圖像重建

    X射線激發(fā)傳感和光學(xué)分子傳感是利用兩種不同的相互作用方式和傳感方法激發(fā)生物組織中的納米分子探針,并對信號發(fā)射位置進(jìn)行檢測和定位的成像技術(shù).基于X射線的激發(fā)傳感技術(shù)在近十年時間內(nèi)發(fā)展迅速,為深入開展組織生理學(xué)及病理生理學(xué)研究奠定了基礎(chǔ)[1].X射線作為探針激發(fā)源的優(yōu)勢在于穿透性強(qiáng)并且具有廣泛的生物醫(yī)學(xué)成像場景.但X射線成像通?;诠怆娦?yīng),峰值衰減在keV能量范圍內(nèi),分子敏感性較低[2].相比而言,光學(xué)分子傳感成像提供了優(yōu)越的分子靈敏度.目前,隨著探測器和傳感器在探測范圍和光子靈敏度等方面的不斷發(fā)展,提高了光學(xué)傳感器件在信號檢測、采集方面的能力,為X射線在分子成像領(lǐng)域的應(yīng)用奠定了基礎(chǔ)[3].

    為進(jìn)一步提高分子敏感性,X射線磷光斷層成像(X-ray luminescence computed tomography,XLCT)作為一種混合成像模態(tài),通過使用傳統(tǒng)X射線的組織感知力和光學(xué)分子成像的對比能力,在深度外顯率和空間分辨率方面呈現(xiàn)出潛在優(yōu)勢[4].該技術(shù)通過X射線激發(fā)生物組織體內(nèi)的磷光探針產(chǎn)生可見光或近紅外,而后利用光敏探測器接收光學(xué)信號,最后通過圖像重建方法恢復(fù)出測量組織中磷光探針空間分布.目前,XLCT成像系統(tǒng)根據(jù)激發(fā)模式的不同主要分為窄束XLCT和錐束XLCT,兩類系統(tǒng)具有不同的硬件設(shè)計和成像性能[5].窄束XLCT通?;谟邢抟暯浅上窦夹g(shù),需要按照特定的順序激發(fā)X射線并獲取測量結(jié)果,具有一定的采樣復(fù)雜度.錐束XLCT基于全域激發(fā)可實現(xiàn)快速測量,因此在X射線利用率與數(shù)據(jù)采集模式上具有較大優(yōu)勢.然而,由于XLCT成像過程中的高度不適定性,重建圖像質(zhì)量相對較低.此外,錐束XLCT缺乏窄束XLCT中的X射線位置信息,進(jìn)一步加劇了成像目標(biāo)重建難度.目前,錐束XLCT成像仍存在空間分辨率不佳、重建精度受限等問題,從而限制了其面向高質(zhì)量快速組織體成像的應(yīng)用[6].因此,進(jìn)一步研究錐束XLCT圖像質(zhì)量的改善方法對于活體生物和生化過程的無創(chuàng)成像至關(guān)重要.

    目前,用于改善錐束XLCT成像不適定性的重建方法可分為兩類:一類為基于光學(xué)傳輸方程構(gòu)建代數(shù)迭代的傳統(tǒng)重建方法,以得到測量值與估計值之間的最佳匹配,如L1正則化、L2正則化方法等[7].該類算法通常具有嚴(yán)謹(jǐn)?shù)臄?shù)學(xué)理論支撐,并能得到穩(wěn)定且有效的解.然而,這類方法過度依賴于逆向問題求解模型的設(shè)計,且不同成像樣本也面臨反復(fù)調(diào)節(jié)參數(shù)以優(yōu)化重建結(jié)果的問題,重建效果受主觀因素影響較大,圖像重建質(zhì)量受限[8].另一類方法為深度學(xué)習(xí)重建方法,其核心思想為通過大型數(shù)據(jù)集學(xué)習(xí)逆問題的未知解,建立更符合實際環(huán)境的非線性模型.與傳統(tǒng)方法相比,深度學(xué)習(xí)方可有效避免研究學(xué)者將精力放在特征提取上,解決了人工提取特征不完整的問題[9].此外,隨著近年來大數(shù)據(jù)時代的到來及人工智能技術(shù)的發(fā)展,深度學(xué)習(xí)技術(shù)在光學(xué)分子成像領(lǐng)域發(fā)展如火如荼.如Gao等[10]提出了基于多層感知機(jī)的生物發(fā)光斷層成像(bioluminescence tomography,BLT)逆問題模擬方法,采用神經(jīng)網(wǎng)絡(luò)擬合組織表面光子密度與組織體內(nèi)發(fā)光源之間的非線性關(guān)系,實驗結(jié)果表明該方法顯著提高了BLT在活體腫瘤中的重建質(zhì)量.Lun等[11]使用殘差網(wǎng)絡(luò)(ResNet)實現(xiàn)了多個投影角度下的聚焦X射線磷光斷層成像,有效驗證了深度學(xué)習(xí)方法應(yīng)用于窄束XLCT重建的可行性.綜上所述,傳統(tǒng)方法通常存在模型精度誤差,且由于逆問題的嚴(yán)重病態(tài)性導(dǎo)致求解不穩(wěn)定,圖像重建的空間分辨率和量化度較差.深度學(xué)習(xí)方法發(fā)展迅速,具有廣泛的應(yīng)用潛力,但對于不同類型的光學(xué)成像技術(shù)需設(shè)計出符合實際需求的個性化網(wǎng)絡(luò).

    因此,為進(jìn)一步改善錐束XLCT圖像重建質(zhì)量,本文基于深度學(xué)習(xí)框架的自適應(yīng)參數(shù)學(xué)習(xí)特性,構(gòu)建了一種用于錐束XLCT定量磷光分布的分組注意力殘差網(wǎng)絡(luò)重建方法,該方法僅需通過簡單的訓(xùn)練、優(yōu)化,便可得到成像物體表面光子密度與內(nèi)部磷光探針濃度分布間的非線性映射.

    具體來說,該方法通過在殘差網(wǎng)絡(luò)中融合注意力機(jī)制,實現(xiàn)從錐束XLCT光學(xué)信號中高效學(xué)習(xí)到有助于目標(biāo)光源重建的有效特征表示,快速逼近成像過程中的近紅外光傳播函數(shù).此外,本文將磷光探針的濃度分布范圍作為錐束XLCT的先驗知識集成到網(wǎng)絡(luò)的損失函數(shù)中,使得從訓(xùn)練數(shù)據(jù)中學(xué)習(xí)的模型符合可接受的領(lǐng)域規(guī)則,增強(qiáng)網(wǎng)絡(luò)模型的魯棒性.最后通過仿真實驗及仿體實驗驗證了本文方法的可行性和有效性.

    1 理論方法

    1.1 成像模型

    一般來說,錐束X射線激發(fā)發(fā)光在生物組織 中的傳播可以描述為3個過程:X射線在生物組織 體內(nèi)的傳播模型、X射線激發(fā)磷光納米粒子產(chǎn)生近 紅外光子模型以及近紅外光子在生物組織體內(nèi)的傳播模型.

    首先,高能X射線源發(fā)射錐形X射線束并穿透生物組織.根據(jù)朗伯比爾定律,X射線在組織體內(nèi)的傳播模型[12]為

    其次,深埋在成像物體內(nèi)的磷光探針與X射線相互作用,受激產(chǎn)生的近紅外光強(qiáng)度為

    最后,近紅外光與組織體發(fā)生吸收、散射等相互作用后的能量變化可以用輻射傳輸方程(radiative transfer equation,RTE)描述[13].由于RTE方程求解較為復(fù)雜,光在高散射弱吸收組織中的傳輸過程通常近似為擴(kuò)散方程(diffusion equation,DE)來表示,即

    根據(jù)微分方程理論,求解擴(kuò)散方程需具備合適的邊界條件.本文采用羅賓邊界條件進(jìn)行補(bǔ)充[14],即

    1.2 基于分組注意力殘差網(wǎng)絡(luò)的錐束XLCT重建 方法

    1.2.1 分組注意力殘差網(wǎng)絡(luò)

    鑒于本文需對多個不同投影特征通道間的依賴關(guān)系進(jìn)行重要性解析.因此本文使用簡潔而有效的通道注意力機(jī)制(channel attention,CA)作為信息增強(qiáng)模塊[20],并根據(jù)錐束XLCT圖像重建特點,對注意力機(jī)制的作用域按照錐束XLCT測量模式進(jìn)行設(shè)計和優(yōu)化,進(jìn)一步處理RB單元中卷積模塊的輸出.如圖1(c)所示,CA作用于特征圖時,首先輸入維度為××的特征圖(為特征圖的高度,為特征圖的寬度,為特征圖的通道數(shù)),之后以兩種不同的池化形式充分提取通道信息,將各通道表示為最大池化特征和平均池化特征,得到兩個大小為1×1×的一維通道描述矢量.其次,將兩個一維矢量傳入卷積核數(shù)分別為/(為通道縮減率)和、卷積核大小為 1×1的兩層卷積層,學(xué)習(xí)各通道的依賴程度,分別生成兩個大小為1×1×的通道注意力特征向量.最后將向量對應(yīng)元素求和,使用Sigmoid激活函數(shù)生成數(shù)值介于0~1之間的通道注意力權(quán)重.將通道注意力權(quán)重值逐一加權(quán)到原特征通道上,根據(jù)依賴程度的不同更新各通道的特征權(quán)重,實現(xiàn)自適應(yīng)特征細(xì)化,增強(qiáng)特征的可分辨性,加權(quán)圖大小與原輸入特征圖大小保持一致.此外,本文采用基于多頭輸入的分組卷積[21]提取注意力信息,有效降低計算復(fù)雜度并擴(kuò)大通道注意力范圍,學(xué)習(xí)更為豐富的信息表示.如圖1(b)所示,將通道注意力的單輸入按通道維度方向進(jìn)行平均拆分,形成多組并行的4個輸入后分別進(jìn)行CA操作.最后,將處理后的分組特征按照通道維度連接(concatenate).此外,由于分組設(shè)計會帶來各組間信息孤立問題[22],本文在通道維度連接后引入卷積尺寸為1×1的卷積操作,避免引入過多參數(shù)的同時促進(jìn)組間信息流通與組內(nèi)信息融合.最后,將分組CA模塊嵌入RB中構(gòu)成RB+CA單元.與RB單元共同作為分組注意力殘差網(wǎng)絡(luò)中的基本單元(圖1(b)).因此,RB+CA單元可在網(wǎng)絡(luò)訓(xùn)練過程中實現(xiàn)對隱含在特征空間不同感知野范圍下的錐束XLCT光子密度信息進(jìn)行綜合挖掘,獲得各層次下對重建結(jié)果具有關(guān)鍵性影響的特征通道并進(jìn)行優(yōu)選,從而實現(xiàn)對于錐束XLCT四角度投影測量結(jié)果中有效部分的自動注意.

    圖1 網(wǎng)絡(luò)結(jié)構(gòu)圖及通道注意力模塊

    如圖1(a)所示,分組注意力殘差網(wǎng)絡(luò)主要由3個部分組成.首先將每個成像樣本在4個投影角度下所測量的光子密度信息作為網(wǎng)絡(luò)輸入,依次通過殘差基本單元RB及RB+CA組合基本單元,由淺至深進(jìn)行特征提取并捕獲關(guān)鍵信息通道.其次,從第2組RB單元開始,由于降維卷積的使用,特征圖每經(jīng)過一組RB單元,圖像大小減半,通道數(shù)為之前的2倍.最后,在一系列的特征提取后,連接1層用于緩解過擬合的全局池化層和1層用于生成重建仿體內(nèi)各個節(jié)點處磷光分布的全連接層.因本文所用仿體剖分后的節(jié)點數(shù)為9446個,因此全連接層所連接神經(jīng)元數(shù)量為9446.基于上述對于分組注意力殘差網(wǎng)絡(luò)的描述,本文通過此架構(gòu)以網(wǎng)絡(luò)的形式將給定數(shù)據(jù)集轉(zhuǎn)換為具有一系列網(wǎng)絡(luò)參數(shù)的非線性映射并動態(tài)學(xué)習(xí)映射內(nèi)各項參數(shù),完成從錐束XLCT測量光子信息中識別及提取有效特征表示并最終實現(xiàn)磷光探針分布重建的過程.

    1.2.2 領(lǐng)域知識損失函數(shù)

    為使得模型在弱信號、高噪聲的情況下也能產(chǎn)生良好的重建性能,提高模型的泛化能力及魯棒性,本文采用領(lǐng)域知識損失函數(shù)作為網(wǎng)絡(luò)訓(xùn)練過程中的損失函數(shù)[23],其主要原理為在TensorFlow庫中的平滑平均絕對誤差(Huber loss)基礎(chǔ)上引入近似域約束損失函數(shù),即將結(jié)果變量的定量范圍(磷光探針的濃度值分布范圍)作為自適應(yīng)領(lǐng)域知識,將其結(jié)合到深度前饋神經(jīng)網(wǎng)絡(luò)的訓(xùn)練損失函數(shù)中,與Huber loss共同作用于網(wǎng)絡(luò)的權(quán)值參數(shù)更新過程.領(lǐng)域知識損失函數(shù)計算公式為

    1.2.3 數(shù)據(jù)準(zhǔn)備與網(wǎng)絡(luò)訓(xùn)練

    通過計算機(jī)仿真模擬實際實驗獲取訓(xùn)練及測試數(shù)據(jù)集.本文采用由南加州大學(xué)與UCLA醫(yī)學(xué)院聯(lián)合開發(fā)的標(biāo)準(zhǔn)數(shù)字鼠模型[24]制作實驗所用仿體,該標(biāo)準(zhǔn)數(shù)字鼠橫斷面的最大直徑約為33mm,總長度大約為87.2mm.由于此標(biāo)準(zhǔn)數(shù)字鼠的結(jié)構(gòu)較為復(fù)雜,因此本文在此基礎(chǔ)上對數(shù)字鼠模型進(jìn)行一定的簡化處理,即僅研究軀干部分的幾個主要器官區(qū)域,例如心臟、肺部、肝臟、腎臟等.

    如圖2所示,本文設(shè)計了一個高為30mm、半徑為15mm的空心圓柱,并將數(shù)字鼠軀干簡化部分嵌入圓柱體成像腔內(nèi),構(gòu)建數(shù)字鼠軀干仿體.另外,在數(shù)字鼠軀干仿體與空心圓柱的縫隙區(qū)域,其光學(xué)參數(shù)被設(shè)定為與肌肉組織相同[25-26].采用有限元剖分軟件COMSOLTM將模型離散為9446個節(jié)點與51594個四面體.隨著血液循環(huán),磷光納米探針通常會在肝臟腫瘤中積聚.因此,在仿真實驗中,本文以肝臟區(qū)域的200個隨機(jī)點作為光源中心,構(gòu)造半徑為1.5~2.5mm,高為1mm的圓柱形光源.接下來,將單光源進(jìn)行兩兩組合,使得雙光源位置中心距離在7~9mm內(nèi)隨機(jī)變化.共構(gòu)建18850組雙目標(biāo)數(shù)據(jù)集.仿體內(nèi)光源區(qū)域磷光探針濃度隨機(jī)設(shè)置為1~3mg/mL,仿體內(nèi)非光源區(qū)域磷光探針濃度設(shè)置為0,根據(jù)上述雙光源特征進(jìn)行錐束XLCT前向仿真.其中,仿體內(nèi)不同器官光學(xué)特性參數(shù)如表1所示.

    圖2 數(shù)字鼠仿體

    表1 不同器官光學(xué)特性參數(shù)

    實驗硬件平臺采用Intel(R)Xeon(R)CPU E5-2650 v4 @2.20GHz處理器,NVIDIA GTX 1070 Ti顯卡(可使用CUDA加速模塊),深度學(xué)習(xí)框架采用基于Python的TensorFLow開源框架,Python版本為3.6.2,Tensorflow版本為2.3.0.

    2 實驗結(jié)果與分析

    2.1 評價指標(biāo)

    為了對圖像重建結(jié)果做出客觀公正的評價,本文實驗采用位置誤差(location error,LE)[27]、均方誤差(mean square error,MSE)[28]、戴斯相似性指數(shù)Dice[29]進(jìn)行定量分析,計算公式如下.

    LE用于評估圖像內(nèi)光源位置準(zhǔn)確度.LE值越小,表明重建目標(biāo)體的幾何中心與真實目標(biāo)體的幾何中心差異越小,其表達(dá)式為

    MSE用于評估圖像重建精度.MSE越小,表明重建的圖像更接近真實圖像.MSE的表達(dá)式為

    Dice用于衡量重建圖像與真實圖像的相似度.當(dāng)Dice為100%時,表示相同的圖像.其表達(dá)式為

    2.2 仿真實驗

    2.2.1 不同靶心距下重建性能比較

    實驗將設(shè)置7mm和9mm兩種靶心距測試不同方法重建精度,兩種情況對應(yīng)的雙光源的位置中心坐標(biāo)分別為(-5mm,0,14.5mm)、(2mm,0,14.5mm)和(-4mm,0,14.5mm)、(5mm,0,14.5mm).雙目標(biāo)體半徑均為2.5mm,磷光探針的真實濃度為1mg/mL.圖3 ART方法、ISPR-ART方法、Tikhonov方法、ResNet18方法、分組注意力殘差方法在=14.5mm切片處重建結(jié)果的二維截面圖.其中,黑色圓圈代表真實的納米發(fā)光目標(biāo)所在區(qū)域,紅色區(qū)域代表重建的納米發(fā)光目標(biāo),目標(biāo)體顏色越紅,表示重建得到的節(jié)點處磷光探針濃度值越高.從圖3中可以看出,傳統(tǒng)方法中的ART、Tikhonov幾乎無法清晰分辨雙目標(biāo)體,重建效果較差.ISPR-ART通過將可行性區(qū)域進(jìn)行迭代收縮,基本上可以重建出目標(biāo)體的濃度及位置,明顯改善了圖像重建質(zhì)量.但在靶心距為7mm時,ISPR-ART出現(xiàn)了一定的目標(biāo)黏連,即使在靶心距為9mm仍無法有效分離雙目標(biāo)體.

    相比3種傳統(tǒng)方法,ResNet18及分組注意力殘差網(wǎng)絡(luò)方法對于相鄰靶目標(biāo)的參數(shù)恢復(fù)準(zhǔn)確度有了較大提升,能夠明顯區(qū)分雙目標(biāo)體并重建出其位置及濃度,同時在視覺上保持了一個較好的目標(biāo)形狀.然而,ResNet18模型在靶心距為7mm時出現(xiàn)了一側(cè)目標(biāo)重建區(qū)域明顯偏離真實目標(biāo)區(qū)域的情況,很大程度上是由于ResNet18模型僅使用殘差單元無法充分聚焦局部關(guān)鍵特征,因此在重建質(zhì)量方面的性能受限,對于復(fù)雜雙目標(biāo)的模型收斂結(jié)果不如分組注意力殘差網(wǎng)絡(luò).

    圖4為不同靶心距下各方法在=0mm處的濃度曲線.可以看出,無論在何種靶心距時,傳統(tǒng)方法中僅ISPR-ART可以基本恢復(fù)出雙磷光靶體濃度,但在非目標(biāo)靶體區(qū)域的重建濃度與真實濃度誤差較大.其次,分組注意力殘差網(wǎng)絡(luò)的預(yù)測濃度值明顯優(yōu)于ResNet18網(wǎng)絡(luò)重建結(jié)果,與真實濃度間差異較?。硗?,本文提出的方法重建結(jié)果始終保持在合理的濃度約束范圍之內(nèi),充分驗證了本文在損失函數(shù)中加入近似域約束將重建結(jié)果限制在定量范圍的有效性,具有更為穩(wěn)定的重建性能.

    由表2可得,傳統(tǒng)方法中的ART算法、Tikhonov方法各項指標(biāo)相對較低.兩種深度學(xué)習(xí)方法與傳統(tǒng)方法相比,可以得到較為理想的定量指標(biāo)結(jié)果.在靶心距為7mm和9mm兩種情況下,采用分組注意力殘差網(wǎng)絡(luò)方法重建位置誤差最小,LE分別為0.37mm和0.30mm,均小于0.38mm,明顯優(yōu)于傳統(tǒng)重建方法和ResNet18.這說明分組注意力殘差方法可以在保證獲得復(fù)雜分類特征的同時,提高了細(xì)節(jié)特征提取能力,同時這也是其他算法無法比擬的優(yōu)勢.

    圖3 不同靶心距雙目標(biāo)重建結(jié)果剖面圖(z=14.5mm)

    圖4 不同靶心距各方法的重建濃度曲線比較

    表2 不同靶心距各方法量化重建結(jié)果

    2.2.2 不同噪聲水平下重建性能比較

    為評估本文方法在不同信噪比環(huán)境下的魯棒性,本組實驗分別在信噪比(SNR)為30dB、40dB、50dB 3種噪聲水平下測試,靶心距為9mm,目標(biāo)體半徑從2.5mm縮小至2.0mm,加大重建難度,光源的真實濃度為1mg/mL.雙目標(biāo)的中心坐標(biāo)分別為(4mm,0,14.5mm)和(5mm,0,14.5mm).圖5為不同信噪比下ResNet18及分組注意力殘差網(wǎng)絡(luò)在=14.5mm切片處重建得到的二維截面圖,圖6為信噪比為30dB時兩種深度學(xué)習(xí)方法的重建濃度曲線比較.結(jié)合圖5的定量重建結(jié)果及圖6的各方法重建濃度曲線可以得出:在靶目標(biāo)尺寸縮小的情況下,ResNet18網(wǎng)絡(luò)重建結(jié)果出現(xiàn)濃度預(yù)測值與真實濃度值偏差過大、不能均勻預(yù)測左側(cè)和右側(cè)目標(biāo)體的濃度等結(jié)果,磷光濃度分布重建能力較弱;而集成了注意力模塊的殘差網(wǎng)絡(luò),重建目標(biāo)的大小、濃度分布也更加均勻,對半徑為2.0mm的目標(biāo)依然可以得到清晰的磷光分布圖像.

    圖5 不同信噪比各方法重建結(jié)果剖面圖(z=14.5mm)

    圖6 信噪比為30dB時重建濃度曲線比較

    表3展示了不同信噪比環(huán)境下兩種深度學(xué)習(xí)方法的量化重建結(jié)果,可以得出,盡管隨著噪聲的增加會使得重建問題的病態(tài)性加重,但兩種深度學(xué)習(xí)算法的定量重建結(jié)果如LE、MSE等仍相對穩(wěn)定,進(jìn)一步突出了模型的魯棒性.其中,分組注意力殘差網(wǎng)絡(luò)的位置誤差相比于ResNet18網(wǎng)絡(luò)下降了約0.04mm,Dice提升了4%左右.由此可見,單一的堆疊殘差基本單元并不能較為有效地提升模型性能,在ResNet18主干網(wǎng)絡(luò)的基礎(chǔ)上加入分組注意力模塊增強(qiáng)了模型對于目標(biāo)的感知能力,尤其是小目標(biāo)體的磷光探針濃度分布恢復(fù)上.此外,將此組實驗結(jié)果與圖3中半徑為2.5mm的目標(biāo)重建實驗對比也可以得出,目標(biāo)體尺寸越大,位置定位越準(zhǔn)確,濃度恢復(fù)程度越高.

    仿真實驗表明,分組注意力殘差方法能夠有效分辨距離較近的靶目標(biāo).與ResNet18網(wǎng)絡(luò)重建結(jié)果相比,該方法增強(qiáng)了對于小目標(biāo)的感知能力,可以更為精確地恢復(fù)出磷光探針的空間分布,具有更好的指標(biāo)量化結(jié)果.此外,與傳統(tǒng)方法相比,該方法在定位精度、空間分辨率、穩(wěn)健性等多方面顯著提高了雙目標(biāo)光源成像質(zhì)量.

    表3 不同噪聲水平各方法量化重建結(jié)果

    2.3 仿體實驗

    為進(jìn)一步驗證本文方法的有效性和穩(wěn)健性,本文設(shè)計了一種圓柱形聚甲醛仿體(如圖7所示),該仿體的直徑為20mm,高為40mm,仿體對于近紅外光子的吸收系數(shù)為0.0072mm-1,散射系數(shù)約為0.72mm-1.使用注射器將由1%的脂肪乳、2%的瓊脂以及Gd2O2S:Eu3+配置而成的2mg/mL目標(biāo)溶液置入外徑為5mm的試管中,溶液深度為5mm.將試管嵌入仿體內(nèi),以仿體底面圓心為坐標(biāo)原點,雙目標(biāo)仿體中心位置坐標(biāo)分別為(-4.5mm,0,27.5mm)和(5.5mm,0,27.5mm).成像系統(tǒng)采用錐束XLCT成像系統(tǒng).X射線源的工作電壓為50kV,工作電流為1mA.CCD相機(jī)作為光學(xué)檢測儀器用于測量仿體在0°(初始位置設(shè)為0°)、90°、180°、270° 4個投影角度下的表面光子密度信息,用于后續(xù)的圖像重建.每次掃描CCD相機(jī)的曝光時間設(shè)置為2s,增益為260.

    圖7 仿體結(jié)構(gòu)示意

    圖8為分別采用ART、ISPR-ART、Tikhonov、ResNet18、分組注意力殘差網(wǎng)絡(luò)在=27.5mm處的重建剖面圖.圖9為仿體實驗各方法重建濃度曲線比較,結(jié)合二者可以看出,3種傳統(tǒng)重建算法在4個投影角度下的重建圖像對實驗環(huán)境中的測量噪聲敏感,重建目標(biāo)邊界模糊,無法準(zhǔn)確定位目標(biāo)體位置及恢復(fù)目標(biāo)體濃度,且濃度預(yù)測值相較于正確濃度存在較大差距.此外,因仿體實驗所用目標(biāo)體尺寸較大,所以在濃度恢復(fù)上兩類殘差模型都表現(xiàn)出了相對良好的重建精度,但ResNet18出現(xiàn)了雙目標(biāo)重建區(qū)域明顯超出真實光源區(qū)域的情況.

    圖8 仿體實驗各方法重建結(jié)果剖面圖(z=27.5mm)

    圖9 仿體實驗各方法重建濃度曲線比較

    結(jié)合表4可知,相較于ResNet18,分組注意力殘差方法重建光源中心更接近光源中心,其LE為0.48mm,明顯小于ResNet18的0.52mm.此外,混合了分組注意力機(jī)制的殘差網(wǎng)絡(luò)重建信號對比度明顯,MSE、Dice等量化指標(biāo)相比ResNet18均有一定的改善,更好地反映了目標(biāo)準(zhǔn)確收斂的特性.

    仿體實驗結(jié)果表明,相較于其他重建方法,分組注意力殘差網(wǎng)絡(luò)在恢復(fù)目標(biāo)體的中心位置、相似度上具有非常大的優(yōu)勢.面對實際實驗中不確定的噪聲,LE低至0.48mm且Dice仍能保持在79%的水平以上,驗證了本文所提出方法的泛化能力以及穩(wěn)健性.

    表4 仿體實驗量化重建結(jié)果

    3 結(jié) 語

    本文針對錐束XLCT圖像空間分辨率不高、傳統(tǒng)模型過度依賴于輻射傳輸方程的問題,提出了一種基于分組注意力殘差網(wǎng)絡(luò)的重建策略.相較于傳統(tǒng)殘差網(wǎng)絡(luò)ResNet18而言,該網(wǎng)絡(luò)可實現(xiàn)跨通道維度的有效特征識別及關(guān)鍵信息捕捉,作為一種高效的信息增強(qiáng)機(jī)制改善學(xué)習(xí)過程并提升網(wǎng)絡(luò)表達(dá)能力,為解決逆問題提供了一種更為精確的重建方法.結(jié)果表明,本文提出的方法在雙目標(biāo)數(shù)字鼠仿真實驗中,可清晰分辨半徑為2.0mm的小目標(biāo)體,LE均小于0.37mm,Dice可達(dá)到82%左右,且在不同噪聲水平下重建結(jié)果較為穩(wěn)定.綜上所述,本文方法相較于其他方法而言,可以精確、穩(wěn)健地重建出目標(biāo)體的形狀、大小、位置及濃度,具有優(yōu)異的圖像重建性能.

    雖然本文提出的重建方法在一定程度上提高了錐束XLCT的圖像質(zhì)量.但是,新的方法也面臨著新的挑戰(zhàn).在未來將關(guān)注以下挑戰(zhàn):首先,在網(wǎng)絡(luò)的架構(gòu)設(shè)計方面.現(xiàn)階段網(wǎng)絡(luò)可實現(xiàn)雙目標(biāo)光源重建任務(wù),考慮到未來臨床研究中腫瘤等病變部位數(shù)量及位置的不確定性、測量誤差及環(huán)境干擾等問題.因此從實際應(yīng)用角度出發(fā),未來將基于分組注意力殘差方法設(shè)計錐束XLCT多目標(biāo)成像方案,進(jìn)一步證明本文方法的有效性.同時,對于復(fù)雜目標(biāo)的重建研究也將有助于探索精細(xì)化的模型體系及加速收斂的優(yōu)化技巧,提高解決重建任務(wù)的能力.另外,在數(shù)據(jù)集的構(gòu)建方法上.現(xiàn)階段數(shù)據(jù)集無法窮盡所有目標(biāo)狀態(tài),即不同的成像對象需要根據(jù)相應(yīng)的輪廓建立不同的訓(xùn)練數(shù)據(jù)集,但真實人體內(nèi)病變目標(biāo)具有不同的大小和形狀,較難獲得足夠數(shù)量的病灶數(shù)據(jù),且仿真數(shù)據(jù)與實驗數(shù)據(jù)之間差異不可忽略.未來將采用從具有已知光學(xué)特性的病灶幾何形狀創(chuàng)建訓(xùn)練數(shù)據(jù)對,真實數(shù)據(jù)用于評估的數(shù)據(jù)集構(gòu)建策略,縮小仿真數(shù)據(jù)和實驗數(shù)據(jù)之間的差異性.最后,基于錐束XLCT成像在提高腫瘤檢測精度方面的臨床應(yīng)用價值.未來將基于本文方法開展大量生物分布成像研究如小鼠原位膠質(zhì)瘤實驗.通過在體實驗發(fā)現(xiàn)問題、解決問題,推動錐束XLCT在腫瘤早期檢測中的發(fā)展.

    [1] Bradley W G. History of medical imaging[J]. Proceed-ings of the American Philosophical Society,2008,152(3):349-361.

    [2] Ginat D T,Gupta R. Advances in computed tomography imaging technology[J]. Annual Review of Biomedical Engineering,2014,16(1):431-453.

    [3] Jaffer F A. Molecular imaging in the clinical arena[J]. JAMA,2005,293(7):855-862.

    [4] Pratx G,Carpenter C M,Sun C,et al. Tomographic molecular imaging of X-ray-excitable nanoparticles[J]. Optics Letters,2010,35(20):3345-3347.

    [5] 束月霞,趙麗麗,蔣皆恢,等. X射線發(fā)光光學(xué)斷層成像的研究進(jìn)展[J]. 科學(xué)通報,2017,62(33):3838-3850.

    Shu Yuexia,Zhao Lili,Jiang Jiehui,et al. Research progress of X-ray luminescence optical tomography[J]. Science Bulletin,2017,62(33):3838-3850(in Chi-nese).

    [6] Chen D,Zhu S,Yi H,et al. Cone beam X-ray lumi-nescence computed tomography:A feasibility study[J]. Medical Physics,2013,40(3):031111.

    [7] Zhao J,Guo H,Yu J,et al. A robust elastic net-l112 reconstruction method for X-ray luminescence computed tomography[J]. Physics in Medicine and Biology,2021,66(19):195005.

    [8] Arridge S R,Schotland J C. Optical tomography:Forward and inverse problems[J]. Inverse Problems,2009,25(12):1541-1548.

    [9] Lucas A,Iliadis M,Molina R,et al. Using deep neural networks for inverse problems in imaging:Beyond ana-lytical methods[J]. IEEE Signal Processing Magazine,2018,35(1):20-36.

    [10] Gao Y,Wang K,An Y,et al. Nonmodel-based bioluminescence tomography using a machine-learning reconstruction strategy[J]. Optica,2018,5(11):1451-1454.

    [11] Lun M C,Cong W,Arifuzzaman M,et al. Focused X-ray luminescence imaging system for small animals based on a rotary gantry[J]. Journal of Biomedical Op-tics,2021,26(3):036004.

    [12] 王正坤,王 鵬,徐昊昊. CT系統(tǒng)成像原理[J]. 信息記錄材料,2018,19(2):86-88.

    Wang Zhengkun,Wang Peng,Xu Haohao. Imaging principle of CT system[J]. Information Recording Mate-rials,2018,19(2):86-88(in Chinese).

    [13] 徐可欣,高 峰,趙會娟. 生物醫(yī)學(xué)光子學(xué)[M]. 北京:科學(xué)出版社,2011.

    Xu Kexin,Gao Feng,Zhao Huijuan. Biomedical Photonics[M]. Beijing:Science Press,2011(in Chi-nese).

    [14] Klose A D,Larsen E W. Light transport in biological tissue based on the simplified spherical harmonics equa-tions[J]. Journal of Computational Physics,2006,220(1):441-470.

    [15] Zhang Y,Lun M C,Li C,et al. Method for improving the spatial resolution of narrow X-ray beam based X-ray luminescence computed tomography imaging[J]. Journal of Biomedical Optics,2019,24(8):1-11.

    [16] 董 芳,侯榆青,余景景,等. 結(jié)合區(qū)域收縮和貪婪策略的熒光分子斷層成像[J]. 激光與光電子學(xué)進(jìn)展,2016,53(1):175-181.

    Dong Fang,Hou Yuqing,Yu Jingjing,et al. Fluores-cence molecular tomography combined with region con-traction and greedy strategy[J]. Progress in Laser and Optoelectronics,2016,53(1):175-181(in Chinese).

    [17] Zhang G,Liu F,Liu J,et al. Cone beam X-ray lumi-nescence computed tomography based on Bayesian method[J]. IEEE Transactions on Medical Imaging,2017,36(1):225-235.

    [18] He K,Zhang X,Ren S,et al. Deep residual learning for image recognition[C]// 2016 IEEE Conference on Computer Vision and Pattern Recognition(CVPR). Las Vegas,USA,2016:770-778.

    [19] He K,Zhang X,Ren S,et al. Identity mappings in deep residual networks[C]// European Conference on Computer Vision. Amsterdam,the Netherlands,2016:630-645.

    [20] Woo S,Park J,Lee J Y,et al. CBAM:Convolutional block attention module[C]// European Conference on Computer Vision,Glasgow. Munich,Germany,2018:3-19.

    [21] Krizhevsky A,Sutskever I,Hinton G. ImageNet classi-fication with deep convolutional neural networks[J]. Communications of the ACM,2017,60(6):84-90.

    [22] Zhang T,Qi G,Xiao B,et al. Interleaved group con-volutions for deep neural networks[C]// International Conference on Computer Vision. Venice,Italy,2017:4373-4382.

    [23] Muralidhar N,Islam M R,Marwah M,et al. Incorporating prior domain knowledge into deep neural networks[C]// 2018 IEEE International Conference on Big Data. Seattle,USA,2018:36-45.

    [24] Dogdas B,Stout D,Chatziioannou A F,et al. Digimouse:A 3D whole body mouse atlas from CT and cryosection data[J]. Physics in Medicine and Biology,2007,52(3):577-587.

    [25] Liu X,Liao Q,Wang H. Fast X-ray luminescence com-puted tomography imaging[J]. IEEE Transactions on Bio-medical Engineering,2014,61(6):1621-1627.

    [26] 侯榆青,賈 濤,易黃建,等. 基于改進(jìn)譜投影梯度算法的X射線發(fā)光斷層成像[J]. 光學(xué)精密工程,2017,25(1):42-49.

    Hou Yuqing,Jia Tao,Yi Huangjian,et al. X-ray lumi-nescence tomography based on improved spectral projection gradient algorithm[J]. Optical Precision Engineer-ing,2017,25(1):42-49(in Chinese).

    [27] Gao P,Rong J,Pu H,et al. Sparse view cone beam X-ray luminescence tomography based on truncated singular value decomposition[J]. Optics Express,2018,26(18):23233-23250.

    [28] Feng J,Sun Q,Li Z,et al. Back-propagation neural network based reconstruction algorithm for diffuse optical tomography[J]. Journal of Biomedical Optics,2018,24(5):1-12.

    [29] Zhang W,Zhu D,Lun M,et al. Multiple pinhole collimator based X-ray luminescence computed tomography[J]. Biomedical Optics Express,2016,7(7):2506.

    A Group Attention Residual Network for Image Reconstruction in Cone-Beam XLCT

    Zhou Zhongxing1, 2,Guo Siqi1, 2,Jia Mengyu1, 2,Zhang Lin2,Wu Yue2

    (1. School of Precision Instruments and Optoelectronics Engineering,Tianjin 300072,China;2. Tianjin Key Laboratory of Biomedical Detection Technology and Instrument,Tianjin 300072,China)

    Cone-beam X-ray luminescence computed tomography(CB-XLCT)is an emerging hybrid imaging modality that can simultaneously obtain both anatomical structure information and functional metabolism information,providing good prospects for the early diagnosis and targeted treatment of diseases. However,CB-XLCT generally suffers from high ill-posedness during image reconstruction,severely deteriorating the spatial resolution. To solve this problem,a quantitative phosphorescence distribution reconstruction method for CB-XLCT based on a group attention residual network is proposed. In this method,the attention module is used to optimize the weight of each photon density feature channel,which is combined with the approximate domain constraint loss function to improve the feature expression process and enhance the reconstruction accuracy and model robustness. Results show that the proposed method can clearly distinguish and reconstruct two targets with the target center distance of 7mm,yielding approximately 0.37mm location error(LE)and 84% Dice similarity index in the simulation experiment. In the phantom experiment,the LE is approximately 0.48mm,and the Dice reaches 79%. Therefore,the group attention residual network method can effectively improve the quality of image reconstruction in LE,spatial resolution,and robustness.

    cone-beam X-ray luminescence computed tomography(CB-XLCT);deep learning;group attention;image reconstruction

    TK448.21

    A

    0493-2137(2022)10-1082-11

    10.11784/tdxbz202108051

    2021-08-21;

    2021-12-23.

    周仲興(1979— ),男,博士,副教授.

    周仲興,zhouzhongxing@tju.edu.cn.

    國家自然科學(xué)基金資助項目(81971656,62175183);天津市自然科學(xué)基金資助項目(19JCYBJC28600).

    the National Natural Science Foundation of China(No. 81971656,No. 62175183),the Natural Science Foundation of Tianjin,China(No. 19JCYBJC28600).

    (責(zé)任編輯:孫立華)

    猜你喜歡
    方法
    中醫(yī)特有的急救方法
    中老年保健(2021年9期)2021-08-24 03:52:04
    高中數(shù)學(xué)教學(xué)改革的方法
    河北畫報(2021年2期)2021-05-25 02:07:46
    化學(xué)反應(yīng)多變幻 “虛擬”方法幫大忙
    變快的方法
    兒童繪本(2020年5期)2020-04-07 17:46:30
    學(xué)習(xí)方法
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    最有效的簡單方法
    山東青年(2016年1期)2016-02-28 14:25:23
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    香蕉丝袜av| 久99久视频精品免费| 亚洲人成电影免费在线| 久久久久国产一级毛片高清牌| 99热国产这里只有精品6| 很黄的视频免费| 久久久久视频综合| 最新美女视频免费是黄的| 免费在线观看日本一区| 中文字幕精品免费在线观看视频| av不卡在线播放| 老熟妇乱子伦视频在线观看| 一级a爱片免费观看的视频| 在线天堂中文资源库| 91国产中文字幕| 中文字幕人妻熟女乱码| 久久ye,这里只有精品| 久久 成人 亚洲| 免费在线观看视频国产中文字幕亚洲| 国产一区在线观看成人免费| 最近最新免费中文字幕在线| 久久影院123| 国产精品1区2区在线观看. | 免费人成视频x8x8入口观看| 欧美中文综合在线视频| 变态另类成人亚洲欧美熟女 | 99国产极品粉嫩在线观看| 久久香蕉国产精品| 精品一区二区三区四区五区乱码| 在线观看舔阴道视频| 韩国精品一区二区三区| 午夜福利在线免费观看网站| 亚洲精品乱久久久久久| 亚洲色图av天堂| 精品一区二区三卡| 视频在线观看一区二区三区| 淫妇啪啪啪对白视频| 久久久精品区二区三区| 女人久久www免费人成看片| 欧美人与性动交α欧美精品济南到| 两个人看的免费小视频| 飞空精品影院首页| 黄色女人牲交| 国产精品 国内视频| 国产男女内射视频| 久久香蕉精品热| 在线永久观看黄色视频| 青草久久国产| av网站在线播放免费| 久久久久国产一级毛片高清牌| 宅男免费午夜| 他把我摸到了高潮在线观看| 俄罗斯特黄特色一大片| 日韩有码中文字幕| 亚洲欧美日韩高清在线视频| 黄网站色视频无遮挡免费观看| 久久亚洲精品不卡| 满18在线观看网站| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 纯流量卡能插随身wifi吗| 日本精品一区二区三区蜜桃| 亚洲人成电影观看| svipshipincom国产片| 曰老女人黄片| 日韩 欧美 亚洲 中文字幕| 在线天堂中文资源库| 美女高潮喷水抽搐中文字幕| 亚洲精品美女久久久久99蜜臀| 精品国产乱码久久久久久男人| 18禁裸乳无遮挡免费网站照片 | 亚洲男人天堂网一区| 久久国产精品大桥未久av| 交换朋友夫妻互换小说| 欧美人与性动交α欧美精品济南到| 成人三级做爰电影| 国产精品电影一区二区三区 | 久久精品国产综合久久久| 免费日韩欧美在线观看| ponron亚洲| avwww免费| 国产av精品麻豆| 大香蕉久久成人网| 久久久久精品国产欧美久久久| 在线十欧美十亚洲十日本专区| 涩涩av久久男人的天堂| 久久国产精品大桥未久av| 欧美日韩福利视频一区二区| 免费av中文字幕在线| 91在线观看av| 国产精品电影一区二区三区 | 欧美激情高清一区二区三区| 亚洲一区二区三区欧美精品| 少妇的丰满在线观看| 亚洲情色 制服丝袜| 日韩免费高清中文字幕av| 成年人午夜在线观看视频| 999久久久国产精品视频| 国产主播在线观看一区二区| 波多野结衣一区麻豆| 丁香六月欧美| 国产一区有黄有色的免费视频| 狠狠狠狠99中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人免费电影在线观看| 午夜福利免费观看在线| 亚洲第一欧美日韩一区二区三区| 国产免费男女视频| 久久久久久人人人人人| 亚洲,欧美精品.| 亚洲aⅴ乱码一区二区在线播放 | 99在线人妻在线中文字幕 | 国产视频一区二区在线看| 99久久人妻综合| 亚洲av美国av| 少妇裸体淫交视频免费看高清 | 999久久久国产精品视频| 日韩 欧美 亚洲 中文字幕| 日韩三级视频一区二区三区| 亚洲一码二码三码区别大吗| 久久久国产成人免费| 中出人妻视频一区二区| 亚洲色图综合在线观看| 9热在线视频观看99| 国产亚洲精品第一综合不卡| 欧美日本中文国产一区发布| 日韩大码丰满熟妇| 精品久久蜜臀av无| 日日爽夜夜爽网站| 老汉色∧v一级毛片| 国产成人免费无遮挡视频| 久久 成人 亚洲| 久久久久国产精品人妻aⅴ院 | 亚洲一码二码三码区别大吗| 免费在线观看亚洲国产| 正在播放国产对白刺激| 成人精品一区二区免费| 999久久久精品免费观看国产| 黑人巨大精品欧美一区二区mp4| 丰满饥渴人妻一区二区三| 99re在线观看精品视频| 91av网站免费观看| 国产熟女午夜一区二区三区| 国产深夜福利视频在线观看| 国产在线观看jvid| 中文字幕人妻丝袜一区二区| 精品福利永久在线观看| 精品一区二区三区四区五区乱码| 精品第一国产精品| 亚洲一区中文字幕在线| av天堂在线播放| 黄片小视频在线播放| 免费高清在线观看日韩| 最新的欧美精品一区二区| 精品久久蜜臀av无| 黑人操中国人逼视频| 每晚都被弄得嗷嗷叫到高潮| 99国产精品免费福利视频| 久久性视频一级片| 亚洲在线自拍视频| 一级a爱视频在线免费观看| 日韩 欧美 亚洲 中文字幕| 嫁个100分男人电影在线观看| 久久 成人 亚洲| 免费少妇av软件| 日韩大码丰满熟妇| 丝袜美足系列| www.自偷自拍.com| 岛国毛片在线播放| 欧美在线一区亚洲| 高清在线国产一区| 国产一区二区三区视频了| 精品久久久久久久毛片微露脸| 久久久久国产一级毛片高清牌| 精品久久久精品久久久| 亚洲国产看品久久| 桃红色精品国产亚洲av| av一本久久久久| 黄网站色视频无遮挡免费观看| 亚洲久久久国产精品| 国产欧美日韩综合在线一区二区| 国产精品一区二区在线不卡| 12—13女人毛片做爰片一| 午夜福利在线免费观看网站| 女人精品久久久久毛片| 下体分泌物呈黄色| 一本一本久久a久久精品综合妖精| 亚洲熟女精品中文字幕| 国精品久久久久久国模美| 午夜免费成人在线视频| 中文字幕另类日韩欧美亚洲嫩草| 性少妇av在线| 9191精品国产免费久久| 久久国产乱子伦精品免费另类| 久久中文字幕一级| 亚洲第一欧美日韩一区二区三区| 亚洲国产精品合色在线| 80岁老熟妇乱子伦牲交| 十八禁人妻一区二区| 久久ye,这里只有精品| 又紧又爽又黄一区二区| 99热只有精品国产| 69精品国产乱码久久久| 欧美乱码精品一区二区三区| 精品人妻在线不人妻| 久久午夜亚洲精品久久| 王馨瑶露胸无遮挡在线观看| 欧美在线黄色| 免费看十八禁软件| 久久国产亚洲av麻豆专区| 日韩欧美国产一区二区入口| 99精国产麻豆久久婷婷| 久久久久视频综合| 精品少妇一区二区三区视频日本电影| 中文字幕人妻丝袜一区二区| 国产区一区二久久| 手机成人av网站| 亚洲成人国产一区在线观看| 亚洲五月婷婷丁香| 女人精品久久久久毛片| 女人爽到高潮嗷嗷叫在线视频| 好男人电影高清在线观看| 久久人妻福利社区极品人妻图片| x7x7x7水蜜桃| a级毛片黄视频| 99精品欧美一区二区三区四区| 日本vs欧美在线观看视频| 免费日韩欧美在线观看| 欧美黑人精品巨大| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品国产综合久久久| 日韩 欧美 亚洲 中文字幕| 免费看十八禁软件| 久久久久久久久免费视频了| 极品少妇高潮喷水抽搐| 欧美在线一区亚洲| 黄色怎么调成土黄色| 免费看a级黄色片| 人成视频在线观看免费观看| 9色porny在线观看| 大型黄色视频在线免费观看| 亚洲av第一区精品v没综合| 9色porny在线观看| 丝瓜视频免费看黄片| 免费人成视频x8x8入口观看| 精品一区二区三区av网在线观看| av线在线观看网站| 男女下面插进去视频免费观看| 精品一区二区三区av网在线观看| 国产野战对白在线观看| 国产精品 国内视频| 国产1区2区3区精品| 国产不卡一卡二| 久99久视频精品免费| 手机成人av网站| 欧美日本中文国产一区发布| 老司机福利观看| 在线观看66精品国产| 女人久久www免费人成看片| 自拍欧美九色日韩亚洲蝌蚪91| 久久香蕉精品热| 亚洲精品美女久久av网站| 国产欧美日韩综合在线一区二区| 91在线观看av| 国产欧美日韩一区二区三| 亚洲精品自拍成人| 国产熟女午夜一区二区三区| 日日爽夜夜爽网站| 欧美亚洲日本最大视频资源| 美女高潮到喷水免费观看| 美女扒开内裤让男人捅视频| 免费在线观看影片大全网站| 精品国产一区二区三区久久久樱花| www.熟女人妻精品国产| 嫩草影视91久久| 高清黄色对白视频在线免费看| 精品一区二区三区av网在线观看| 波多野结衣av一区二区av| 视频在线观看一区二区三区| 怎么达到女性高潮| 亚洲av成人不卡在线观看播放网| 国产精品久久久久久人妻精品电影| 国产精品电影一区二区三区 | 女同久久另类99精品国产91| 在线观看www视频免费| 99精国产麻豆久久婷婷| 久久热在线av| 久久久国产欧美日韩av| 99国产精品一区二区蜜桃av | av有码第一页| 久久草成人影院| 麻豆乱淫一区二区| 亚洲精品国产精品久久久不卡| 精品国产一区二区久久| 韩国av一区二区三区四区| 亚洲精品久久成人aⅴ小说| 99香蕉大伊视频| 精品久久久精品久久久| 在线观看免费视频网站a站| 欧美精品一区二区免费开放| 一级毛片女人18水好多| 涩涩av久久男人的天堂| 1024视频免费在线观看| 国产精品电影一区二区三区 | 欧美丝袜亚洲另类 | 亚洲国产毛片av蜜桃av| а√天堂www在线а√下载 | 午夜福利,免费看| 国产极品粉嫩免费观看在线| 高潮久久久久久久久久久不卡| 国产蜜桃级精品一区二区三区 | 日韩欧美一区二区三区在线观看 | 亚洲少妇的诱惑av| 免费在线观看亚洲国产| 人人妻,人人澡人人爽秒播| 欧美日韩福利视频一区二区| 久久久久久久国产电影| 午夜精品久久久久久毛片777| 国产亚洲一区二区精品| 一区二区三区精品91| av免费在线观看网站| 日本a在线网址| 香蕉久久夜色| 欧美激情高清一区二区三区| 男男h啪啪无遮挡| 91精品三级在线观看| 中文字幕最新亚洲高清| 亚洲精品一二三| 国产精品永久免费网站| 午夜亚洲福利在线播放| 香蕉久久夜色| 性色av乱码一区二区三区2| 国产成人免费无遮挡视频| 在线观看免费视频网站a站| 99在线人妻在线中文字幕 | 日韩一卡2卡3卡4卡2021年| 午夜福利影视在线免费观看| 欧美在线黄色| 99香蕉大伊视频| 亚洲午夜精品一区,二区,三区| 亚洲精品国产精品久久久不卡| 欧美日韩亚洲国产一区二区在线观看 | 精品少妇久久久久久888优播| 亚洲国产毛片av蜜桃av| 成年人黄色毛片网站| 成人永久免费在线观看视频| 久久中文看片网| 亚洲免费av在线视频| 亚洲中文字幕日韩| 国内毛片毛片毛片毛片毛片| 黄色丝袜av网址大全| 亚洲美女黄片视频| 黄片大片在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 一本一本久久a久久精品综合妖精| 91麻豆av在线| www.自偷自拍.com| 久久性视频一级片| 9热在线视频观看99| 亚洲性夜色夜夜综合| 国产高清videossex| 亚洲国产精品sss在线观看 | 黄片播放在线免费| 老汉色av国产亚洲站长工具| 一级毛片精品| 91成年电影在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲中文av在线| 看片在线看免费视频| 久久久久久久午夜电影 | 免费日韩欧美在线观看| 亚洲片人在线观看| 国产主播在线观看一区二区| 十八禁人妻一区二区| 国产在视频线精品| 他把我摸到了高潮在线观看| 午夜老司机福利片| 侵犯人妻中文字幕一二三四区| 日韩熟女老妇一区二区性免费视频| 无人区码免费观看不卡| 国产视频一区二区在线看| 欧美日韩亚洲综合一区二区三区_| 丰满迷人的少妇在线观看| 久久国产精品男人的天堂亚洲| 一级作爱视频免费观看| 免费观看精品视频网站| 三上悠亚av全集在线观看| 狠狠狠狠99中文字幕| 精品欧美一区二区三区在线| 国产视频一区二区在线看| 国产亚洲精品久久久久5区| 国产99久久九九免费精品| 久久亚洲精品不卡| 国产97色在线日韩免费| 黄色 视频免费看| 午夜福利在线观看吧| 好看av亚洲va欧美ⅴa在| 在线天堂中文资源库| 男女高潮啪啪啪动态图| 亚洲欧美色中文字幕在线| 99热网站在线观看| 午夜福利视频在线观看免费| 一级片'在线观看视频| 一a级毛片在线观看| 亚洲精品国产一区二区精华液| 国产精品香港三级国产av潘金莲| 午夜成年电影在线免费观看| 1024香蕉在线观看| 777米奇影视久久| 精品电影一区二区在线| 成人亚洲精品一区在线观看| 最新美女视频免费是黄的| 777久久人妻少妇嫩草av网站| 一二三四社区在线视频社区8| 亚洲五月婷婷丁香| 真人做人爱边吃奶动态| 久久久久久免费高清国产稀缺| 日韩熟女老妇一区二区性免费视频| 国产单亲对白刺激| 91九色精品人成在线观看| 国产在线观看jvid| 又黄又粗又硬又大视频| 一边摸一边做爽爽视频免费| 成年动漫av网址| 99riav亚洲国产免费| 啦啦啦视频在线资源免费观看| av中文乱码字幕在线| 啪啪无遮挡十八禁网站| 精品国内亚洲2022精品成人 | 亚洲熟女毛片儿| 99国产综合亚洲精品| 国产高清国产精品国产三级| 亚洲国产毛片av蜜桃av| 欧美日韩成人在线一区二区| 在线国产一区二区在线| 巨乳人妻的诱惑在线观看| 久久久久国内视频| 黄色女人牲交| 国产一卡二卡三卡精品| 久久久精品国产亚洲av高清涩受| 中文亚洲av片在线观看爽 | 窝窝影院91人妻| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 乱人伦中国视频| 久久香蕉激情| 午夜两性在线视频| 亚洲av欧美aⅴ国产| 老司机深夜福利视频在线观看| 成人免费观看视频高清| 人妻久久中文字幕网| 夜夜躁狠狠躁天天躁| 日韩精品免费视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区视频在线观看免费 | 黄色怎么调成土黄色| 丰满的人妻完整版| 久久久久国产一级毛片高清牌| 丝瓜视频免费看黄片| 99国产综合亚洲精品| 脱女人内裤的视频| a级片在线免费高清观看视频| 午夜久久久在线观看| 日本五十路高清| 在线观看日韩欧美| 在线观看www视频免费| 亚洲五月天丁香| 757午夜福利合集在线观看| 久久中文字幕人妻熟女| 99国产精品99久久久久| 黄色视频不卡| 久久国产精品人妻蜜桃| 久久精品aⅴ一区二区三区四区| 少妇的丰满在线观看| 国产在线一区二区三区精| 在线观看日韩欧美| 丰满饥渴人妻一区二区三| 久久天躁狠狠躁夜夜2o2o| 欧美性长视频在线观看| 丝袜美腿诱惑在线| 丝袜美足系列| 美女高潮到喷水免费观看| 欧美最黄视频在线播放免费 | 欧美精品一区二区免费开放| 美国免费a级毛片| 国产一区二区三区综合在线观看| 国产人伦9x9x在线观看| 精品视频人人做人人爽| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三区在线| 亚洲熟女精品中文字幕| 国产无遮挡羞羞视频在线观看| 国产黄色免费在线视频| 一级a爱视频在线免费观看| 久久久久久久精品吃奶| 超碰成人久久| 性色av乱码一区二区三区2| 国产成+人综合+亚洲专区| 视频区欧美日本亚洲| 91精品国产国语对白视频| www.999成人在线观看| 国产人伦9x9x在线观看| 国产精品国产av在线观看| 免费久久久久久久精品成人欧美视频| 美女午夜性视频免费| 久久久久视频综合| av中文乱码字幕在线| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 老司机午夜十八禁免费视频| 久久精品国产a三级三级三级| 极品人妻少妇av视频| 国产日韩欧美亚洲二区| 国产97色在线日韩免费| 色在线成人网| 国产成人欧美在线观看 | 成年动漫av网址| 久久久久国内视频| 黄色成人免费大全| 母亲3免费完整高清在线观看| 国产精品乱码一区二三区的特点 | 成人av一区二区三区在线看| 国产成人影院久久av| 动漫黄色视频在线观看| 精品熟女少妇八av免费久了| 韩国精品一区二区三区| 日本wwww免费看| 久久婷婷成人综合色麻豆| 18禁黄网站禁片午夜丰满| 女人高潮潮喷娇喘18禁视频| 可以免费在线观看a视频的电影网站| 午夜影院日韩av| 免费在线观看亚洲国产| 亚洲美女黄片视频| 18禁裸乳无遮挡动漫免费视频| 亚洲国产毛片av蜜桃av| 国产精品永久免费网站| 男女床上黄色一级片免费看| 91精品国产国语对白视频| 免费观看人在逋| 又黄又爽又免费观看的视频| 18在线观看网站| 悠悠久久av| 久99久视频精品免费| 涩涩av久久男人的天堂| 国产精品久久电影中文字幕 | 中文欧美无线码| 999久久久国产精品视频| 后天国语完整版免费观看| 欧美日韩亚洲高清精品| 国产成人av教育| 国产亚洲精品第一综合不卡| 亚洲综合色网址| 丰满的人妻完整版| 50天的宝宝边吃奶边哭怎么回事| 国产精品.久久久| 午夜福利影视在线免费观看| 69av精品久久久久久| 亚洲熟女毛片儿| 欧美日韩一级在线毛片| 亚洲美女黄片视频| 法律面前人人平等表现在哪些方面| 1024香蕉在线观看| 精品少妇久久久久久888优播| 99精国产麻豆久久婷婷| 国产免费男女视频| 国产成人精品久久二区二区91| 国产精品亚洲一级av第二区| 久久香蕉激情| 麻豆国产av国片精品| 男人舔女人的私密视频| 久久这里只有精品19| 精品久久久久久,| 日韩欧美在线二视频 | 99久久综合精品五月天人人| 国产欧美日韩一区二区三| 丰满迷人的少妇在线观看| 国产午夜精品久久久久久| 成人永久免费在线观看视频| 国产成人精品无人区| 亚洲欧美激情在线| 99香蕉大伊视频| av中文乱码字幕在线| 久久草成人影院| a级片在线免费高清观看视频| 男男h啪啪无遮挡| 99国产精品一区二区蜜桃av | 亚洲人成电影免费在线| 亚洲美女黄片视频| 777米奇影视久久| 日本欧美视频一区| 亚洲国产中文字幕在线视频| 99热只有精品国产| 久久 成人 亚洲| 俄罗斯特黄特色一大片| 亚洲色图av天堂| 波多野结衣一区麻豆| 男女午夜视频在线观看| 美女扒开内裤让男人捅视频| 美女高潮喷水抽搐中文字幕| 成人精品一区二区免费| а√天堂www在线а√下载 | 黄片大片在线免费观看| 91字幕亚洲| 一级片'在线观看视频| 人妻一区二区av| 99久久99久久久精品蜜桃| 在线看a的网站| 久久久精品区二区三区| 超色免费av| 最新美女视频免费是黄的| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| 亚洲五月婷婷丁香| 99久久人妻综合|