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

    基于極坐標的扇束CT周期性優(yōu)化算法

    2016-11-08 01:44:33桂志國
    測試技術(shù)學報 2016年5期
    關(guān)鍵詞:極坐標周期性余弦

    劉 輝,張 權(quán),劉 祎,桂志國,2

    (1.中北大學 電子測試技術(shù)重點實驗室,山西 太原 030051;2.中北大學 儀器科學與動態(tài)測試教育部重點實驗室,山西 太原 030051)

    ?

    基于極坐標的扇束CT周期性優(yōu)化算法

    劉輝1,張權(quán)1,劉祎1,桂志國1,2

    (1.中北大學 電子測試技術(shù)重點實驗室,山西 太原 030051;2.中北大學 儀器科學與動態(tài)測試教育部重點實驗室,山西 太原 030051)

    為了提高傳統(tǒng)扇束計算機斷層掃描(Computed Tomography,CT)卷積反投影重建算法(Conventional Back-Projection,CBP )的重建速度,提出了一種快速的扇束CT重建算法.該算法結(jié)合極坐標特性,利用三角函數(shù)一定程度上的周期性,對傳統(tǒng)CBP重建算法進行了改進,減少了正余弦函數(shù)的計算量和重建所需要的外圍循環(huán)次數(shù).同樣利用該特性對極坐標重建后的圖像坐標轉(zhuǎn)換過程進行改善,從而對CBP算法進行了改進,降低了重建時間.

    計算機斷層掃描;卷積反投影;三角函數(shù)特性;極坐標

    隨著計算機斷層成像技術(shù)(Computed Tomography,CT)的不斷發(fā)展,CT重建技術(shù)越來越多地運用到醫(yī)學診斷[1]、工業(yè)無損檢測等行業(yè)中[2].CT重建主要有解析算法和迭代算法[3-5].與迭代算法相比,解析重建法以其重建速度快、重建圖像質(zhì)量高的特點,已經(jīng)被廣泛應用到工程上.CBP重建算法是解析算法中最具代表性的重建算法之一[6-8].雖然CBP重建算法重建所需數(shù)據(jù)量較少,但是運算數(shù)據(jù)量依然很大,還需要進一步提高其重建速度,才能更好地運用到工業(yè)上.

    加權(quán)、卷積濾波和反投影是CBP算法的3個主要部分,其中反投影過程就占據(jù)了98.36%的重建時間[9],因此反投影重建過程的改進是主要研究方向.文獻[10]運用極坐標特性,提出了極坐標反投影重建算法(Polar Coordinates Back-Projection,PCBP),但極坐標向笛卡爾坐標的轉(zhuǎn)換過程中要使用插值運算,重建速度提高有限;文獻[11]利用正余弦函數(shù)特性,將16幅投影同時進行反投影重建,提出了對極坐標重建的優(yōu)化算法(Optimize Polar Conventional Back-Projection OPCBP),該算法將重建速度提高了4倍;文獻[12]利用了三角函數(shù)的對稱性,提出了極坐標反投影對稱優(yōu)化算法(Symmetry Polar Conventional Back-Projection SPCBP),同時對64幅預處理后的數(shù)據(jù)進行重建,將重建速度提高8倍.

    本文利用三角函數(shù)一定程度上的周期性和極坐標反投影重建的特點,將144幅反投影數(shù)據(jù)同時進行反投影重建,使外圍循環(huán)減少至傳統(tǒng)算法的1/144,同時大大減少了三角函數(shù)的運算量,從而實現(xiàn)對傳統(tǒng)算法的優(yōu)化.本文稱之為極坐標反投影周期性優(yōu)化算法(Periodicity Polar-coordinate Conventional Back-Projection,PPCBP).

    1 扇束CBP重建原理

    扇束CT掃描示意圖如圖1所示,圖中a(r,θ)是待重建圖像切片的射線衰減系數(shù)分布函數(shù);β為射線源投射角度;函數(shù)p(β,s)是該角度下的投影數(shù)據(jù);D為射線源到物體重建中心的距離,重建公式為

    (1)

    (2)

    (3)

    (4)

    (5)

    由式(1)~式(5)可知計算機實現(xiàn)CBP反投影重建時,首先通過式(3)給投影數(shù)據(jù)進行加權(quán);其次通過式(2)對加權(quán)過的數(shù)據(jù)進行卷積濾波;最后根據(jù)式(4)和式(5)確定的射線的反投影地址s和加權(quán)因子U對濾波后的數(shù)據(jù)進行累加反投影重建.假設重建的圖像的大小為N×N,投影角度為[0°,360°),投影視圖數(shù)為M,該算法的反投影過程流程如圖2所示,從圖2中可以看出CBP反投影重建中運最復雜的就是對反投影地址s和加權(quán)因子U求解.由式(4)和式(5)可知,每次求s和U需要9次乘法(含除法)和6次正余弦函數(shù)運算,反投影過程總共需要次乘法和次正余弦運算,因此主要考慮對s和U求解過程的簡化.

    圖1 扇束掃描示意圖Fig.1 The scanning schematic of fan-beam

    圖2 CBP反投影重建過程流程圖Fig.2 The flowchart of CBP back-projection reconstruction

    2 極坐標下進行重建

    (6)

    (7)

    圖3 PCBP反投影重建流程圖Fig.3 The flowchart of PCBP algorithm

    3 三角函數(shù)周期性優(yōu)化

    3.1反投影過程的優(yōu)化

    三角函數(shù)中30°,45°,60°和90°是幾個特殊的角度,三角函數(shù)值很容易計算得出.通過簡單的三角函數(shù)的誘導公式,可以得出120°,135°,150°,180°,210°,225°,240°,270°,300°,315°和330°的三角函數(shù)值,這些角度中除了45°及其誘導產(chǎn)生的角度之外,其余角度都是呈現(xiàn)出的等差關(guān)系,這些角度的正余弦值一定程度上表現(xiàn)出了周期性,給PCBP算法優(yōu)化提供了很好的條件.分別讓β-θ=30m(0≤m≤11,m∈Z),討論其正余弦值并研究各角度下的余弦值之間的關(guān)系

    (8)

    (9)

    (10)

    (11)

    圖4 PPCBP反投影重建流程圖Fig.4 The flowchart of PPCBP algorithm

    式(8)~式(10)之間關(guān)系只是差了一個正負號,而且A,B和A1,B1以及A2,B2具有式(12)所示的關(guān)系,因此只要求解出[0°,30°)范圍內(nèi)的正余弦值就可以求得所有掃描角度下β-θ的正余弦值,極坐標三角函數(shù)周期性優(yōu)化算法(PPCBP)算法的反投影過程流程如圖4所示.

    (12)

    通過以上分析,在一個投影角度β下通過三角函數(shù)的變換,只需要計算兩次正余弦值就可以計算出θ+30m(0≤m≤11,m∈z)下的s和U的值,因此一個循環(huán)可以完成12個投影視角下的144個投影數(shù)據(jù)參與重建,外層循環(huán)減少為極坐標下的1/144,同時將三角函數(shù)計算量降低為極坐標下的1/72,大大提高了重建速度.所以PPCBP算法無論在循環(huán),乘法或者是三角函數(shù)運算上都得到了很大的提高,在理論上證實了PPCBP的優(yōu)越性.

    3.2重建圖像坐標轉(zhuǎn)換的優(yōu)化

    本文在圖像的轉(zhuǎn)換過程中采用雙線性插值方法,如圖5所示,雙線性插值通過臨近目標點G的A,B,C和D 4點經(jīng)過3次線性插值從而求得G點像素值的過程,其原理如式(13)所示,其中P(x)代表了x點的像素值,a=BE/AB,b=DF/CD,c=FG/EF.

    (13)

    針對采用雙線性插值增加的計算量,本文使用正切函數(shù)對稱性對其簡化.圖6中所標記的點在極坐標中具有相同的值,且其正切值之間具有正負和互為倒數(shù)的關(guān)系.因此,如果知道極坐標下任意一個點的坐標,就能同時知道極坐標下其他7個點的坐標轉(zhuǎn)換,使外圍循環(huán)減少至改變前的1/8,同時使正余切函數(shù)值的計算量也減少至改變前的1/8,提高了坐標轉(zhuǎn)換的速度.

    圖5 雙線性插值示意圖Fig.5 The schematic of bilinear interpolation

    圖6 坐標轉(zhuǎn)換中像素點的選取示意圖Fig.6 Pixels chosen schematic during coordinate transformation

    4 實驗結(jié)果

    圖7 各種算法相對CBP重建時間加速比Fig.7 The speedup time compared with CBP

    為了驗證本文算法的可行性,將PTCBP算法和CBP算法、PCBP算法、OPCBP算法和SPCBP算法的重建結(jié)果進行對比.實現(xiàn)本文算法的軟件設備為Windows 7 32 b操作系統(tǒng)和VC++6.0編譯環(huán)境;硬件設備為英特爾 Pentium(奔騰)雙核 E5400@2.70 GHz處理器,DDR3 2G金士頓內(nèi)存.投影數(shù)據(jù)參數(shù)如下:射線源到旋轉(zhuǎn)中心的距離為800 mm,物體中心到探測器的距離為150 mm,射線為ISOVOLT 450,探測器為paxscan2520,探元大小0.127 mm,電壓110 kV,電流2.5 mA,投影數(shù)據(jù)大小為570×460,投影視圖數(shù)為360,重建圖像大小為512×512.

    4.1重建速度比較

    各算法的重建時間如表1所示,加速比如圖7所示:從重建的時間和加速比上可以看出,PTCBP算法的速度是CBP算法的將近12.5倍,是PCBP算法的近10倍,是SPCBP的近1.5倍,是OPCBP的2.1倍.因此,本文所提出的PPCBP算法大大提高了CBP算法的重建速度.由于使用的是三角函數(shù)一定程度上表現(xiàn)出來的周期性原理,使得重建過程中所占用的內(nèi)存也大大減少,在計算機編程中充分展示了PPCBP算法的優(yōu)越性.

    表1 各算法重建速度對比表

    4.2重建圖像質(zhì)量比較

    各種算法圖像如圖8所示,從圖中我們很難看出圖像之間的差別,圖9展示的是各算法重建出圖像的第145列像素大小之間的對比,從中我們發(fā)現(xiàn)各圖像之間差別很小,說明了利用PPCBP算法重建出的圖像與其他重建算法重建出的圖像在質(zhì)量上沒有明顯的差別.為了量化本文算法重建出的圖像與其他算法重建圖像之間的差別,本文采用歸一化均方距離 (Normalized Mean Square Distance,NMSD)和歸一化平均絕對距離(Normal Average Absolute Distance,NAAD)對圖像質(zhì)量進行評估,其公式分別為

    (14)

    (15)

    圖8 各種算法重建圖像質(zhì)量對比圖Fig.8 The quality comparison chart of the images constructed by different algorithms

    圖9 各算法第145列像素大小對比圖Fig.9 The comparison of pixels in the 145th row but from different images reconstructed by vary algorithms.

    算法CBPPCBPPTCBPOPTCBPSPCBPNMSD-0.3881030.3881080.3881040.405673NAAD-0.1984530.1984590.1984540.217913

    由表2可知,無論是從歸一化均方距離還是歸一化平均絕對距離上PTCBP算法跟PCBP,OPTCBP和SPCBP算法重建出的圖像在質(zhì)量上沒有明顯的差別,因此進一步的說明PTCBP算法在重建圖像時并沒有引入誤差,PTCBP算法在保證圖像質(zhì)量上有良好的適用性.

    5 結(jié) 論

    本文提出的PTCBP算法結(jié)合三角函數(shù)一定程度上的周期性以及對稱性分別對反投影過程和坐標轉(zhuǎn)換過程進行了改進.在保證了重建圖像質(zhì)量的前提下,很大程度上提高了傳統(tǒng)CBP算法的重建速度.

    [1]曾更生.醫(yī)學圖像重建[M].北京:高等教育出版社,2010.

    [2]Franco L,Tahoces P G,Martínez-Mera J A.Visualization software for CT:fan/cone beam and metrology applications [J].Procedia Engineering,2013,63(63):779-785.

    [3]莊天戈.CT原理與算法[M].上海:上海交通大學出版社,1992.

    [4]Rit S,Oliva M V,Brousmiche S,et al.The Reconstruction Toolkit (RTK),an open-source cone-beam CT reconstruction toolkit based on the Insight Toolkit (ITK)[C].Journal of Physics:Conference Series.2014,489:012079.

    [5]Wang A S,Stayman J W,Otake Y,et al.Soft-tissue imaging with C-arm cone-beam CT using statistical reconstruction[J].Physics in Medicine and Biology.2014,59(4):1005.

    [6]Zhang F,Du J,Jiang H,et al.Iterative geometric calibration in circular cone-beam computed tomography[J].Optik-International Journal for Light and Electron Optics,2014,125(11):2509-2514.

    [7]Zheng Fang,Xin Gu,Xing Liang Zhang,et al.Research on the relationship between projection number and image noise level in FBP reconstruction [J].Advanced Materials Research.2013,718-720:2324-2328.

    [8]Xin Y,Yang J,Xie Z Q,et al.An overlapping semantic community detection algorithm base on the ARTs multiple sampling models[J].Expert Systems with Applications,2014,42(7):3420-3432.

    [9]劉曉平.扇束卷積反投影法的程序優(yōu)化[J].CT 理論與應用研究.1996,5(1):35-37.

    Liu Xiaoping.Program optimization of swctor-beam convolution backprojection[J].Computerized Tomography Theory and Applications,1996,5(1):35-37.(in Chinese)

    [10]楊民,路宏年,黃朝志.用查找表和極坐標反投影法實現(xiàn)計算機斷層掃描快速重建[J].兵工學報,2004,35(4):476-479.

    Yang Min,Lu Hongnian,Huang Chaozhi.Fast computed tomography reconstruction with look-up table and polar coordonate back-projection[J].Acta Armamentarii,2004,35(4):476-479.(in Chinese)

    [11]張鑾,桂志國.扇束等距CT極坐標反投影重建算法的速度優(yōu)化[J].無損檢測,2010,32(2):95-98.

    Zhang Luan,Gui Zhiguo.Speed optimization of polar coordinate back-projection rectonstruction algorithm for fan beam collinear equispaced CT[J].Nondestructive Testing,2010,32(2):95-98.(in Chinese)

    [12]張晶,張權(quán),劉祎,等.扇束CT極坐標反投影重建算法的對稱優(yōu)化[J].計算應用.2014,34(6):1711-1714.

    Zhang Jing,Zhang Quan,Liu Yi,et al.Symmetry optimization of polar coordinate back-projection reconstruction algorithm for fan beam CT[J].Journal of Computer Applications,2014,34(6):1711-1714.(in Chinese)

    Periodic Optimization Algorithm Based on Fan-Beam CT Polar Coordinates

    LIU Hui1,ZHANG Quan1,LIU Yi1,GUI Zhiguo1,2

    (1.State Key Laboratory of Electronic Measurement Technology (North University),Taiyuan 030051,China;2.Instrument Science & Dynamic Measurement,Ministry of Education Key laboratory of (North University),Taiyuan 030051,China.)

    To improve the reconstruction speed of traditional fan-beam CT image reconstruction based on fan-beam conventional back projection (CBP),a fast polar coordinate back projection reconstruction algorithm was proposed.Combining with the features of polar coordinate,the proposed algorithm applied the trigonometric periodicity in some degree,and improved traditional CBP to reduce the amount of computation of trigonometric function and the times of external cycle in reconstruction process.Similarly,the mentioned properties were also used in the process of the conversion between polar coordinate and Cartesian coordinate to improve the CBP algorithm and the reconstruction time gets reduced.

    computed tomography; conventional back-projection; trigonometric function characteristics; polar coordinate

    1671-7449(2016)05-0394-06

    2016-03-01

    國家自然科學基金資助項目(61271357);國家重大科學儀器設備開發(fā)專項資助項目(2014YQ240445);山西省自然科學基金資助項目(2015011046)

    劉輝(1989-),男,碩士生,主要從事圖像重建的研究.

    TP391.4

    Adoi:10.3969/j.issn.1671-7449.2016.05.005

    猜你喜歡
    極坐標周期性余弦
    巧用極坐標解決圓錐曲線的一類定值問題
    數(shù)列中的周期性和模周期性
    極坐標視角下的圓錐曲線
    一類整數(shù)遞推數(shù)列的周期性
    兩個含余弦函數(shù)的三角母不等式及其推論
    不能忽視的極坐標
    基于擴頻碼周期性的單通道直擴通信半盲分離抗干擾算法
    分數(shù)階余弦變換的卷積定理
    圖像壓縮感知在分數(shù)階Fourier域、分數(shù)階余弦域的性能比較
    離散余弦小波包變換及語音信號壓縮感知
    一进一出抽搐gif免费好疼| 国产午夜精品久久久久久| 日本a在线网址| 中文字幕人妻熟女乱码| 国产aⅴ精品一区二区三区波| 日本熟妇午夜| 久久欧美精品欧美久久欧美| 一区福利在线观看| 成人国产一区最新在线观看| 99热6这里只有精品| 一区二区三区精品91| 日韩高清综合在线| 国产熟女午夜一区二区三区| 亚洲免费av在线视频| 黑人欧美特级aaaaaa片| 亚洲一区中文字幕在线| 一边摸一边做爽爽视频免费| 国产精品综合久久久久久久免费| 精品人妻1区二区| 男女午夜视频在线观看| 美女免费视频网站| 国产一卡二卡三卡精品| www.999成人在线观看| 国产午夜精品久久久久久| 欧美乱妇无乱码| 高清在线国产一区| 亚洲精品一卡2卡三卡4卡5卡| 国产又爽黄色视频| 久久久久久久精品吃奶| 美女扒开内裤让男人捅视频| 成人欧美大片| 日本熟妇午夜| 一区二区三区激情视频| 亚洲美女黄片视频| 亚洲无线在线观看| 欧美成人一区二区免费高清观看 | 欧美日韩瑟瑟在线播放| 少妇 在线观看| 成人特级黄色片久久久久久久| av有码第一页| 日本a在线网址| 大香蕉久久成人网| 99热6这里只有精品| 国产av一区在线观看免费| 啦啦啦观看免费观看视频高清| 久久精品亚洲精品国产色婷小说| 亚洲专区国产一区二区| 午夜视频精品福利| 午夜a级毛片| 国产精品二区激情视频| 日本 欧美在线| 国产激情久久老熟女| 搞女人的毛片| 99re在线观看精品视频| 免费搜索国产男女视频| 中文在线观看免费www的网站 | 草草在线视频免费看| 亚洲av熟女| 国产v大片淫在线免费观看| 亚洲专区字幕在线| 成人三级黄色视频| 啪啪无遮挡十八禁网站| 黄色丝袜av网址大全| 欧美精品亚洲一区二区| 欧美黑人精品巨大| 亚洲精品久久国产高清桃花| 91老司机精品| 亚洲 国产 在线| 久久久久久久午夜电影| 国产日本99.免费观看| 99久久99久久久精品蜜桃| 91字幕亚洲| 91成年电影在线观看| bbb黄色大片| 精品卡一卡二卡四卡免费| 国产精品久久久久久亚洲av鲁大| 香蕉丝袜av| 性欧美人与动物交配| 欧美乱色亚洲激情| 日韩精品免费视频一区二区三区| 亚洲国产欧美一区二区综合| 男女那种视频在线观看| 999久久久国产精品视频| 91成年电影在线观看| 成年免费大片在线观看| 日韩精品中文字幕看吧| 日韩欧美一区视频在线观看| 欧美成人一区二区免费高清观看 | 色av中文字幕| 少妇 在线观看| 亚洲,欧美精品.| 一边摸一边抽搐一进一小说| 人妻丰满熟妇av一区二区三区| 国产一区二区三区视频了| 免费高清视频大片| 免费人成视频x8x8入口观看| 最新在线观看一区二区三区| 亚洲国产欧美网| 国产乱人伦免费视频| 黄色女人牲交| 91国产中文字幕| 午夜免费鲁丝| 午夜免费成人在线视频| 亚洲av熟女| 国产av一区在线观看免费| 黄网站色视频无遮挡免费观看| 国产成年人精品一区二区| 真人一进一出gif抽搐免费| 日日摸夜夜添夜夜添小说| av天堂在线播放| 久久中文字幕一级| 亚洲男人的天堂狠狠| 久久国产精品影院| 亚洲av成人av| 两人在一起打扑克的视频| 亚洲专区字幕在线| 亚洲激情在线av| 国产欧美日韩一区二区精品| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品sss在线观看| 人成视频在线观看免费观看| 日本一本二区三区精品| 欧美又色又爽又黄视频| 欧美zozozo另类| 日本在线视频免费播放| 亚洲中文字幕一区二区三区有码在线看 | 亚洲片人在线观看| 哪里可以看免费的av片| 亚洲性夜色夜夜综合| 免费人成视频x8x8入口观看| 99在线视频只有这里精品首页| 色av中文字幕| 91在线观看av| 少妇粗大呻吟视频| 丰满的人妻完整版| 黄频高清免费视频| 免费在线观看黄色视频的| 日本三级黄在线观看| 久久精品亚洲精品国产色婷小说| 老熟妇仑乱视频hdxx| 国产一卡二卡三卡精品| 亚洲国产欧美一区二区综合| 美女扒开内裤让男人捅视频| 欧美国产精品va在线观看不卡| 日本五十路高清| 久久久久久久久久黄片| АⅤ资源中文在线天堂| 久久久久久九九精品二区国产 | 91字幕亚洲| 国产又色又爽无遮挡免费看| 亚洲狠狠婷婷综合久久图片| 在线观看免费日韩欧美大片| www.精华液| 久久青草综合色| 999久久久精品免费观看国产| 欧美午夜高清在线| 亚洲av电影在线进入| 好男人在线观看高清免费视频 | 亚洲国产日韩欧美精品在线观看 | 久久久精品欧美日韩精品| 色综合站精品国产| 精华霜和精华液先用哪个| 亚洲国产欧洲综合997久久, | 免费在线观看黄色视频的| 欧美黄色淫秽网站| 久久久精品国产亚洲av高清涩受| 久久精品国产99精品国产亚洲性色| 最新在线观看一区二区三区| 久久久久国产一级毛片高清牌| 三级毛片av免费| 母亲3免费完整高清在线观看| av有码第一页| 丝袜美腿诱惑在线| 琪琪午夜伦伦电影理论片6080| 午夜影院日韩av| 91麻豆精品激情在线观看国产| 久久人妻福利社区极品人妻图片| 三级毛片av免费| 50天的宝宝边吃奶边哭怎么回事| 黄网站色视频无遮挡免费观看| 精品国产一区二区三区四区第35| 国产一区二区三区在线臀色熟女| 香蕉国产在线看| 国产高清videossex| 欧美激情久久久久久爽电影| 亚洲美女黄片视频| 成在线人永久免费视频| 欧美日韩一级在线毛片| 亚洲在线自拍视频| 女同久久另类99精品国产91| 一区福利在线观看| cao死你这个sao货| 色av中文字幕| 国产精品98久久久久久宅男小说| 久久精品夜夜夜夜夜久久蜜豆 | videosex国产| www.精华液| 亚洲精品国产区一区二| 久久人妻福利社区极品人妻图片| 中文字幕久久专区| 国产av一区二区精品久久| 别揉我奶头~嗯~啊~动态视频| 欧美 亚洲 国产 日韩一| 美女高潮喷水抽搐中文字幕| 欧美另类亚洲清纯唯美| 波多野结衣高清无吗| 免费人成视频x8x8入口观看| 9191精品国产免费久久| 波多野结衣高清作品| 国产区一区二久久| 老汉色∧v一级毛片| 免费观看人在逋| 精品一区二区三区av网在线观看| 啦啦啦免费观看视频1| 成人手机av| 亚洲午夜理论影院| 1024香蕉在线观看| 精品午夜福利视频在线观看一区| 久久久久久九九精品二区国产 | 久久伊人香网站| 精品电影一区二区在线| 亚洲精华国产精华精| 久久国产亚洲av麻豆专区| 欧美av亚洲av综合av国产av| 一级作爱视频免费观看| 国产精品98久久久久久宅男小说| 丝袜美腿诱惑在线| 窝窝影院91人妻| av欧美777| 99精品久久久久人妻精品| 亚洲一区中文字幕在线| 听说在线观看完整版免费高清| 欧美国产日韩亚洲一区| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕人妻熟女乱码| 国产精品亚洲一级av第二区| 黄频高清免费视频| 人人妻人人看人人澡| 色综合亚洲欧美另类图片| aaaaa片日本免费| 久久香蕉精品热| 欧美激情高清一区二区三区| 真人做人爱边吃奶动态| 国产91精品成人一区二区三区| 桃色一区二区三区在线观看| 亚洲精品国产区一区二| 亚洲国产欧美一区二区综合| 搞女人的毛片| videosex国产| 成年免费大片在线观看| 国产精品乱码一区二三区的特点| 亚洲一区二区三区不卡视频| 国产在线精品亚洲第一网站| 亚洲一区二区三区不卡视频| 亚洲自拍偷在线| 男女之事视频高清在线观看| 精品久久久久久成人av| 国产精品98久久久久久宅男小说| 国产一卡二卡三卡精品| 99在线视频只有这里精品首页| 99国产精品99久久久久| 亚洲精品一区av在线观看| 最近最新中文字幕大全免费视频| 91麻豆av在线| 国产精品国产高清国产av| 午夜福利欧美成人| 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产在线观看| 欧美日韩黄片免| 国产亚洲欧美在线一区二区| 免费在线观看成人毛片| 国产一级毛片七仙女欲春2 | 日韩精品中文字幕看吧| 性色av乱码一区二区三区2| 久久 成人 亚洲| 国产精品久久久久久精品电影 | 亚洲精品色激情综合| 国产v大片淫在线免费观看| www.自偷自拍.com| 国产成人影院久久av| 搞女人的毛片| 亚洲无线在线观看| 真人做人爱边吃奶动态| 亚洲精品国产精品久久久不卡| 成人18禁在线播放| 久久精品aⅴ一区二区三区四区| 国产日本99.免费观看| 12—13女人毛片做爰片一| 欧美成人一区二区免费高清观看 | 色播在线永久视频| 亚洲av日韩精品久久久久久密| 欧美大码av| 在线观看一区二区三区| 免费在线观看亚洲国产| 老熟妇仑乱视频hdxx| 国产精华一区二区三区| 在线观看免费视频日本深夜| 色播亚洲综合网| 欧美乱色亚洲激情| 亚洲九九香蕉| 精品久久久久久久人妻蜜臀av| 免费在线观看完整版高清| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利在线在线| 免费在线观看成人毛片| 97碰自拍视频| 好男人在线观看高清免费视频 | 精品一区二区三区四区五区乱码| ponron亚洲| 久久性视频一级片| 亚洲熟女毛片儿| 日本 欧美在线| 757午夜福利合集在线观看| 亚洲九九香蕉| 又黄又粗又硬又大视频| 久久中文字幕一级| 熟妇人妻久久中文字幕3abv| 亚洲av电影在线进入| 丁香六月欧美| 欧美激情极品国产一区二区三区| 脱女人内裤的视频| 亚洲精品美女久久av网站| 两个人看的免费小视频| 久久久久久久久久黄片| 精品一区二区三区av网在线观看| 天堂动漫精品| 亚洲 欧美一区二区三区| 亚洲成人国产一区在线观看| 天堂影院成人在线观看| 国产免费av片在线观看野外av| 久久久久国产一级毛片高清牌| 亚洲精品一区av在线观看| 国产国语露脸激情在线看| 成人av一区二区三区在线看| 欧美成狂野欧美在线观看| 少妇 在线观看| 男女视频在线观看网站免费 | 男人舔女人下体高潮全视频| 亚洲欧美精品综合久久99| 亚洲国产精品999在线| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 午夜精品久久久久久毛片777| 妹子高潮喷水视频| 久久香蕉激情| ponron亚洲| 在线观看日韩欧美| 午夜福利高清视频| 99热6这里只有精品| 好男人电影高清在线观看| av福利片在线| 好男人电影高清在线观看| 久99久视频精品免费| 午夜视频精品福利| 无人区码免费观看不卡| www.自偷自拍.com| e午夜精品久久久久久久| www.自偷自拍.com| 最新在线观看一区二区三区| 一级a爱片免费观看的视频| 久久草成人影院| 国产亚洲精品第一综合不卡| 十分钟在线观看高清视频www| 国产野战对白在线观看| 特大巨黑吊av在线直播 | 最新美女视频免费是黄的| 91麻豆精品激情在线观看国产| 国产片内射在线| 婷婷六月久久综合丁香| 高清在线国产一区| 深夜精品福利| 黄色毛片三级朝国网站| 视频区欧美日本亚洲| 国产精品野战在线观看| 成在线人永久免费视频| 欧美国产精品va在线观看不卡| 久久精品国产亚洲av高清一级| 不卡av一区二区三区| 欧美日韩乱码在线| 老熟妇乱子伦视频在线观看| 又大又爽又粗| 精品国产亚洲在线| 99在线视频只有这里精品首页| 12—13女人毛片做爰片一| 免费看十八禁软件| 中文字幕精品免费在线观看视频| 18禁国产床啪视频网站| 麻豆国产av国片精品| 看片在线看免费视频| 亚洲精品一区av在线观看| 1024视频免费在线观看| 亚洲av中文字字幕乱码综合 | 色尼玛亚洲综合影院| 日韩精品青青久久久久久| 18禁裸乳无遮挡免费网站照片 | 1024香蕉在线观看| 狂野欧美激情性xxxx| 黄色 视频免费看| 最近最新中文字幕大全免费视频| 久久国产亚洲av麻豆专区| 国产91精品成人一区二区三区| 少妇的丰满在线观看| 国产一区二区三区视频了| 免费人成视频x8x8入口观看| 国产亚洲精品综合一区在线观看 | 国产精品爽爽va在线观看网站 | 一进一出抽搐gif免费好疼| 午夜福利在线观看吧| 一级黄色大片毛片| 国产色视频综合| 午夜视频精品福利| 成年免费大片在线观看| avwww免费| 天堂影院成人在线观看| 久久久水蜜桃国产精品网| 热re99久久国产66热| 叶爱在线成人免费视频播放| 一区二区三区国产精品乱码| 天天躁狠狠躁夜夜躁狠狠躁| 成人手机av| 亚洲午夜精品一区,二区,三区| 欧美黑人欧美精品刺激| 国产v大片淫在线免费观看| 男人舔女人的私密视频| 久久久久九九精品影院| 搞女人的毛片| 免费高清在线观看日韩| 50天的宝宝边吃奶边哭怎么回事| 两个人免费观看高清视频| 性欧美人与动物交配| 亚洲专区中文字幕在线| 男人的好看免费观看在线视频 | 日本成人三级电影网站| 欧美激情高清一区二区三区| 亚洲av中文字字幕乱码综合 | 欧美最黄视频在线播放免费| 黑丝袜美女国产一区| 久久精品亚洲精品国产色婷小说| 亚洲精华国产精华精| 国产99白浆流出| av福利片在线| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 欧美黄色片欧美黄色片| 亚洲av片天天在线观看| 久久精品国产亚洲av高清一级| 一级a爱视频在线免费观看| 高清毛片免费观看视频网站| 一区二区三区国产精品乱码| 国产精品 欧美亚洲| 老汉色∧v一级毛片| 色播亚洲综合网| 久久久久国产精品人妻aⅴ院| 免费女性裸体啪啪无遮挡网站| 搡老熟女国产l中国老女人| 91成年电影在线观看| 看免费av毛片| www.999成人在线观看| 国产精品亚洲美女久久久| АⅤ资源中文在线天堂| 亚洲欧美日韩高清在线视频| 国产精品,欧美在线| 91大片在线观看| 久久中文看片网| 国产精品一区二区免费欧美| 欧美激情 高清一区二区三区| 成人三级做爰电影| 自线自在国产av| 午夜免费激情av| 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲| 国产精品野战在线观看| 成人免费观看视频高清| 成人国产综合亚洲| 国产欧美日韩一区二区精品| 男女午夜视频在线观看| 精品欧美一区二区三区在线| 88av欧美| 一级毛片高清免费大全| 日本a在线网址| 亚洲成人免费电影在线观看| 亚洲欧美激情综合另类| 午夜a级毛片| 亚洲欧美日韩高清在线视频| 俄罗斯特黄特色一大片| 女同久久另类99精品国产91| 中文字幕精品亚洲无线码一区 | 国产精品久久视频播放| 男人舔女人的私密视频| 熟女少妇亚洲综合色aaa.| 波多野结衣高清作品| 搡老熟女国产l中国老女人| 中亚洲国语对白在线视频| 久久久久久人人人人人| 成熟少妇高潮喷水视频| 老汉色av国产亚洲站长工具| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美黑人欧美精品刺激| 久久精品国产综合久久久| 哪里可以看免费的av片| 日本熟妇午夜| 亚洲色图 男人天堂 中文字幕| 午夜影院日韩av| 亚洲欧洲精品一区二区精品久久久| 国产三级黄色录像| 岛国在线观看网站| 久9热在线精品视频| 国产精品久久久av美女十八| 欧美久久黑人一区二区| 欧美日韩亚洲综合一区二区三区_| 精品不卡国产一区二区三区| 国内精品久久久久久久电影| 国产高清激情床上av| 中文字幕av电影在线播放| 大型av网站在线播放| 国产成人欧美| 精品日产1卡2卡| 亚洲精品美女久久久久99蜜臀| 无限看片的www在线观看| 桃红色精品国产亚洲av| 久久久国产欧美日韩av| 中文在线观看免费www的网站 | 在线av久久热| 免费电影在线观看免费观看| 国产亚洲精品一区二区www| 国产成人影院久久av| 精品人妻1区二区| 国产人伦9x9x在线观看| 极品教师在线免费播放| 亚洲电影在线观看av| 女人被狂操c到高潮| 1024香蕉在线观看| 一本综合久久免费| 午夜精品在线福利| 啦啦啦免费观看视频1| 免费观看精品视频网站| а√天堂www在线а√下载| 精品国产超薄肉色丝袜足j| 亚洲成国产人片在线观看| 美女高潮喷水抽搐中文字幕| 免费在线观看视频国产中文字幕亚洲| 俺也久久电影网| 村上凉子中文字幕在线| 国产成人av激情在线播放| 亚洲欧美一区二区三区黑人| 亚洲精品国产区一区二| 成人亚洲精品av一区二区| 久久香蕉国产精品| 亚洲欧美精品综合一区二区三区| 老司机午夜十八禁免费视频| 18禁国产床啪视频网站| 欧美 亚洲 国产 日韩一| 麻豆国产av国片精品| 国产精品香港三级国产av潘金莲| 日日夜夜操网爽| 亚洲国产欧美一区二区综合| 婷婷精品国产亚洲av在线| 一本精品99久久精品77| 黄色视频,在线免费观看| 国产视频一区二区在线看| 欧美丝袜亚洲另类 | av中文乱码字幕在线| 丰满人妻熟妇乱又伦精品不卡| a级毛片a级免费在线| 久久香蕉精品热| 久久久国产成人精品二区| 久久精品影院6| 成年版毛片免费区| 亚洲人成伊人成综合网2020| 精品乱码久久久久久99久播| 熟妇人妻久久中文字幕3abv| 别揉我奶头~嗯~啊~动态视频| 国产高清激情床上av| 91在线观看av| 久久午夜亚洲精品久久| 一区福利在线观看| 91在线观看av| tocl精华| 久久精品夜夜夜夜夜久久蜜豆 | 一级a爱视频在线免费观看| xxxwww97欧美| 欧美黑人欧美精品刺激| 曰老女人黄片| 国产成人一区二区三区免费视频网站| 18禁国产床啪视频网站| 国产精品乱码一区二三区的特点| 欧美人与性动交α欧美精品济南到| 国产高清激情床上av| 韩国av一区二区三区四区| 久久午夜亚洲精品久久| 国产主播在线观看一区二区| 久久久久免费精品人妻一区二区 | 在线观看免费视频日本深夜| xxxwww97欧美| 国产av不卡久久| 久久国产精品影院| 少妇的丰满在线观看| av有码第一页| 成人特级黄色片久久久久久久| 成人三级做爰电影| 亚洲狠狠婷婷综合久久图片| 国产一区二区三区视频了| 午夜成年电影在线免费观看| 丝袜美腿诱惑在线| 欧美一区二区精品小视频在线| 中文字幕最新亚洲高清| av视频在线观看入口| 午夜福利在线在线| 精品国产乱子伦一区二区三区| 最近最新免费中文字幕在线| 亚洲国产精品sss在线观看| 波多野结衣巨乳人妻| 国产色视频综合| 两个人免费观看高清视频| 免费在线观看视频国产中文字幕亚洲|