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

    基于梯度場的工業(yè)X射線圖像增強算法

    2019-11-15 04:49:03周沖劉歡趙愛玲張鵬程桂志國
    計算機應(yīng)用 2019年10期
    關(guān)鍵詞:圖像增強直方圖X射線

    周沖 劉歡 趙愛玲 張鵬程 劉 祎 桂志國

    摘 要:在X射線成像檢測厚薄不均構(gòu)件時,經(jīng)常會出現(xiàn)對比度低或?qū)Ρ榷炔痪约罢斩鹊偷膯栴},這會導(dǎo)致圖像顯示時構(gòu)件的一些細(xì)節(jié)難以被觀察與分析。針對這一問題,提出一種基于梯度場的X射線圖像增強算法。該算法以梯度場增強為核心,分為兩步:首先,提出一種基于對數(shù)變換的算法,壓縮圖像的灰度范圍、去除圖像冗余灰度信息、提升圖像對比度;然后,提出一種基于梯度場的算法,增強圖像細(xì)節(jié)、提升圖像局部對比度、提高圖像質(zhì)量,使構(gòu)件細(xì)節(jié)清晰顯示在檢測屏上。選擇一組厚薄不均構(gòu)件的X射線圖像進行了實驗,并與對比度受限自適應(yīng)直方圖均衡化(CLAHE)、同態(tài)濾波等算法進行了比較。實驗結(jié)果表明所提算法具有更明顯的增強效果,能更好地顯示構(gòu)件的細(xì)節(jié)信息,并且通過計算平均梯度和無參考結(jié)構(gòu)清晰度(NRSS)紋理分析的定量評價標(biāo)準(zhǔn)進一步表明了該算法的有效性。

    關(guān)鍵詞:圖像增強;X射線;梯度場;對數(shù)變換;直方圖

    中圖分類號:TP391.41

    文獻標(biāo)志碼:A

    Abstract:In the detection of components with uneven thickness by X-ray, the problems of low contrast or uneven contrast and low illumination often occur, which make it difficult to observe and analyze some details of components in the images obtained. To solve this problem, an X-ray image enhancement algorithm based on gradient field was proposed. The algorithm takes gradient field enhancement as the core and is divided into two steps. Firstly, an algorithm based on logarithmic transformation was proposed to compress the gray range of an image, remove redundant gray information of the image and improve image contrast. Then, an algorithm based on gradient field was proposed to enhance image details, improve local image contrast and image quality, so that the details of components were able to be clearly displayed on the detection screen. A group of X-ray images of components with uneven thickness were selected for experiments, and the comparisons with algorithms such as Contrast Limited Adaptive Histogram Equalization (CLAHE) and homomorphic filtering were carried out. Experimental? results show that the proposed algorithm has more obvious enhancement effect and can better display the detailed information of the components. The quantitative evaluation criteria of calculating average gradient and No-Reference Structural Sharpness (NRSS) texture analysis further demonstrate the effectiveness of this algorithm.Key words:? image enhancement; X-ray; gradient field; logarithmic transformation; histogram

    0 引言

    X射線成像已經(jīng)成為現(xiàn)代工業(yè)無損檢測不可或缺的一部分,其圖像質(zhì)量直接影響器件檢測效果。工業(yè)X射線圖像的低對比度和低照度特點,直接影響無損檢測人員的檢測速度與精準(zhǔn)度。隨著X射線探測器制造技術(shù)的不斷提高,采集圖像的動態(tài)范圍已經(jīng)由最初的8位發(fā)展到16位,然而顯示器至今依然停留在8位,少數(shù)專用顯示器能顯示12位圖像,但此類顯示器不僅會帶來視覺負(fù)擔(dān),而且其發(fā)展速度也無法追趕上X射線探測器制造技術(shù)的發(fā)展[1]。因此,增強工業(yè)X射線圖像以獲取更清晰和精確的缺陷信息很有必要[2]。

    為了解決X射線探測器制造技術(shù)與顯示器發(fā)展技術(shù)的不同步問題,傳統(tǒng)采用“窗位/窗寬”圖像增強技術(shù)顯示圖像。該技術(shù)是一種基于空域直接灰度變換的技術(shù),窗寬可以看成是一種線性拉伸圖像對比度的技術(shù),窗位結(jié)合窗寬控制圖像顯示在檢測屏的灰度范圍。這是一種全局對比度調(diào)節(jié)的技術(shù),它的缺點是:1)沒有考慮圖像局部信息,全局拉伸對比度會造成圖像細(xì)節(jié)信息丟失;2)調(diào)節(jié)工作繁瑣,無法兼顧整體與局部信息。因此引入圖像增強技術(shù)是一個重要的研究方向[2]。先進的圖像增強技術(shù)大致可以分為三類:1)圖像分解技術(shù)[3],即先將圖像分解為高頻、低頻信號,然后對不同頻域信號進行不同的濾波處理;2)基于變換的技術(shù)[4],即通過利用特定曲線函數(shù)變化性質(zhì)來改變圖像灰度,從而達到提高圖像的對比度的目的;3)直方圖修正技術(shù)[5],即通過對圖像灰度級的非線性修正達到提高圖像增強的目的。線性多尺度方法通過在空間頻率尺度上分解圖像來增強X射線圖像[6],但是由于強邊緣的平滑和噪聲與缺陷檢查相關(guān)的細(xì)節(jié)混合而引入光暈偽影;基于Retinex照射反射模型的算法[7-8]雖然可以提高圖像的視覺效果,但容易產(chǎn)生曝光過度的問題,常用于低照度彩色圖像增強,不適合直接用于低對比度、低照度的工業(yè)X射線圖像增強。基于變換的技術(shù)可以有效地提高圖像的全局對比度,但容易造成局部細(xì)節(jié)對比度低的問題。直方圖均衡化是一種十分有效的對比度增強技術(shù),但這是一種全局對比度修正的技術(shù),容易造成因為顯著改變圖像的亮度而引入光暈的問題。針對這一問題,近年來,很多改進的直方圖均衡化修正技術(shù)被提出來,這些修正技術(shù)也在特定圖像上表現(xiàn)出了更好的增強效果,但是不適合直接用于低對比度、低照度的工業(yè)X射線圖像增強。

    近年來,Li等[9]結(jié)合直方圖均衡化與小波變換提出了一種工業(yè)X射線圖像增強算法,該算法分為兩步:第一步使用對比度受限自適應(yīng)直方圖均衡化(Contrast Limited Adaptive Histogram Equalization, CLAHE)對圖像進行整體對比度調(diào)整,然后使用小波分解對圖像進行局部對比度調(diào)整。該算法存在一定的缺陷,CLAHE會顯著地改變圖像的亮度差異,后續(xù)局部對比度調(diào)整會引入人造光暈。Tan等[10]結(jié)合直方圖,提出了一種自適應(yīng)窗口的技術(shù),該算法首先計算圖像中短時幀直方圖的平均能量值;然后使用雙閾值來檢測包含有用信息的幀;最后,將檢測到的幀的端點灰度值視為窗口端點,通過比較和搜索圖像對比度的最大值,遍歷并確定最合適的幀長度和幀偏移。該算法不能顯示圖像的整體細(xì)節(jié),但可以突出局部感性區(qū)域的細(xì)節(jié)。劉萬程等[11-12]提出了圖像融合的方式增強工業(yè)X射線圖像,該算法通過改變射線能量獲取工件信息,然后通過融合不同能量的射線圖像以達到增強的效果;但是該算法存在操作復(fù)雜、耗時的缺點。為了兼顧圖像整體細(xì)節(jié)與局部細(xì)節(jié),同時避免復(fù)雜、耗時的操作以及光暈的產(chǎn)生,本文提出一種基于梯度場的增強算法,首先通過改進對數(shù)變換的映射算法,去除圖像冗余信息,提高圖像對比度;然后使用改進梯度場的增強算法,增強圖像細(xì)節(jié),提高圖像質(zhì)量;最后將構(gòu)件射線圖像清晰的顯示在檢測屏上,提高工業(yè)X射線檢測速率與精準(zhǔn)度。

    1 梯度場圖像增強模型

    通過操作圖像梯度來執(zhí)行動態(tài)范圍壓縮的梯度域算子幾乎沒有光暈和梯度反轉(zhuǎn)偽影,F(xiàn)attal等[13]通過對每個像素處乘以梯度衰減因子,減小高范圍圖像梯度;隨著梯度減小,通過最小能量泛函來得到最終的圖像,可以表示為:

    其中: I(x, y)表示輸入圖像(x, y)點的像素值。參數(shù)α通過使衰減因子為1來控制某些梯度不變,即梯度為α的地方;梯度小于α的會被輕微放大,梯度大于α的會被抑制。參數(shù)β可以控制衰減速率,β越小衰減得越快。衰減因子的變化隨參數(shù)α、β的變化可以表示如圖1所示。

    2 基于梯度場的工業(yè)X射線增強X探測器平板采集的圖像是16位的高范圍圖像,要把16位的圖像數(shù)據(jù)顯示在8位的顯示器上,就避免不了要對圖像像素范圍進行壓縮。如何在低動態(tài)范圍(Low Dynamic Range,LDR)顯示、提高高范圍圖像的所有細(xì)節(jié)的可見性是可視化技術(shù)的目標(biāo)[14]。傳統(tǒng)的可視化方式是窗位、窗寬線性映射,其公式可以表達如下:

    上述窗位、窗寬技術(shù)是一種全局的線性映射方式。由于沒有考慮局部特征信息,難免會出現(xiàn)局部細(xì)節(jié)信息丟失的問題。本章結(jié)合對數(shù)變換提出一種像素壓縮算法,對數(shù)映射已經(jīng)被應(yīng)用在多種高動態(tài)圖像顯示技術(shù)中[15]。像素壓縮算法可以抽象地表示為:

    其中:L(x, y)表示壓縮后的像素;P表示映射關(guān)系;S(x, y)表示原圖像像素。

    2.1 去除圖像的冗余灰度級

    工業(yè)X射線圖為16位圖像,其像素值灰度分布范圍廣,即0~65535,但是,在分布范圍內(nèi)有很多沒有被占用的灰度級。本節(jié)的目的是去除原圖像中沒有占用的灰度級,實現(xiàn)過程中使用了直方圖信息。算法過程如下:

    其中: p(i)表示輸入圖像中第i個灰度級的像素數(shù)量。經(jīng)過壓縮后,可能會帶來局部對比度降低的問題,因為去除沒有占用的灰度級,可能會導(dǎo)致某些局部細(xì)節(jié)的對比度變小即局部對比度降低。

    2.2 對數(shù)變換

    為了增加圖像整體亮度的同時彌補2.1節(jié)帶來的局部對比度降低問題,可以應(yīng)用對數(shù)函數(shù)。對數(shù)函數(shù)可以增加圖像對比度和亮度,尤其是在低亮度區(qū)域,這種技術(shù)已經(jīng)被多次用在高范圍圖像顯示技術(shù)中[22]。其數(shù)學(xué)公式表示為:

    其中:L(I)表示變換后的輸出像素;I表示輸入圖像;Imin表示輸入圖像的最小像素值;Imax表示輸入圖像的最大像素值;參數(shù)τ可以控制圖像的整體亮度,其值越小圖像越亮,其值越大圖像越暗。

    對2.1節(jié)像素壓縮后的圖片像素在對數(shù)域做歸一化處理。結(jié)合式(9)和式(11),得到本文提出的圖像壓縮映射函數(shù),其數(shù)學(xué)表達式如下:

    其中:I(x, y)表示輸入圖像(x, y)點的像素值;HL為2.1節(jié)生成的灰度映射表;LT(x, y)表示對數(shù)變換后(x, y)點的像素值,同時作為后續(xù)梯度場處理的輸入像素值。

    2.3 改進梯度場增強細(xì)節(jié)

    Fattal等[13]的衰減因子是全局的,沒有考慮局部信息,可以實現(xiàn)適當(dāng)?shù)膭討B(tài)范圍壓縮,但圖像中的細(xì)節(jié)沒有得到有效揭示。受文獻[16]算法的啟發(fā),在此基礎(chǔ)上結(jié)合低對比度、低照度工業(yè)X射線圖像的特征與人類視覺系統(tǒng)的對比感知特征[17],提出了一個自適應(yīng)衰減傳遞函數(shù),其表達式如下:

    前向和后向差分的這種組合確保了div G的近似與用于拉普拉斯算子的中心差分格式一致。有限差分格式產(chǎn)生了一個大的線性方程組——圖像中的每個像素對應(yīng)一個線性方程組,但是對應(yīng)的矩陣在每一行中只有5個非零元素,因為每個像素只與它的4個鄰居相連接。使用全多網(wǎng)格算法來解決這個系統(tǒng),采用高斯賽德爾平滑迭代法[18]。這可以以O(shè)(n)計算復(fù)雜度為代價到達近似解,其中n是圖像中的像素個數(shù)[13]。

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

    實驗的硬件平臺為Intel Core i7-7700 CPU @3.60GHz,軟件仿真環(huán)境是Matlab2016a,VS2015,其中Matlab2016a用于函數(shù)分析,算法仿真在VS2015上實現(xiàn)。采用由實驗室采集的射線圖作為輸入圖像,大小為1024×1024,即長、高都為1024。

    為了檢驗本文算法的有效性,本文選擇兩幅射線圖作為實驗對象。此外選擇文獻[19]的同態(tài)濾波算法、CLAHE、文獻[16]改進的梯度場作對比。對比分析分為主觀視覺對比與客觀評價指標(biāo)兩部分。

    3.1 主觀效果分析

    本節(jié)給出兩張實驗圖像的原圖,以及經(jīng)過上述3種對比算法和本文算法增強后的圖像。最后給出原圖以及處理后圖像的主觀分析。

    圖2~3主觀效果分析:對原圖(a),由于構(gòu)件的厚薄不均,射線能量太大會導(dǎo)致構(gòu)件薄的區(qū)域曝光過度,射線能量太低會導(dǎo)致構(gòu)件厚的區(qū)域無法透射,為了兼顧厚、薄區(qū)域,射線能量不能太高也不能太低,與此同時就會帶來對比度低、照度低的問題,導(dǎo)致圖像細(xì)節(jié)模糊不清;對圖(b),經(jīng)過文獻[19]的同態(tài)濾波算法處理后,圖像整體亮度提升的同時部分細(xì)節(jié)明顯突出,但出現(xiàn)的問題是某些局部細(xì)節(jié)突出得不夠明顯,同時某些局部背景也被改變了;在經(jīng)過CLAHE直方圖均衡化后的圖(c)中,圖像的部分細(xì)節(jié)被顯現(xiàn)出來,同時會出現(xiàn)局部細(xì)節(jié)突出不夠、在某些局部地區(qū)出現(xiàn)光暈現(xiàn)象;經(jīng)過文獻[16]改進梯度域增強算法后的圖(d)中,圖像大部分區(qū)域的細(xì)節(jié)增強比文獻[19]同態(tài)濾波算法處理后的效果好,但是其整體亮度太暗、視覺效果不好、某些局部細(xì)節(jié)也凸顯得不夠好;經(jīng)過本文算法處理后的圖(e)中,不僅能充分顯示各個部分的細(xì)節(jié)、同時整體亮度也符合人眼觀察、圖像細(xì)節(jié)層次感清晰。

    無參考結(jié)構(gòu)清晰度的計算方法簡述如下:1)為待評價圖像構(gòu)造參考圖像;2)提取待評價圖像和構(gòu)造參考圖像的梯度信息;3)找出待評價圖像的梯度圖像中梯度信息最豐富的N個圖像塊;4)計算無參考結(jié)構(gòu)清晰度。

    從表1可以看出,相較于所對比的算法,經(jīng)過本文算法處理后,圖像的平均梯度最大、NRSS系數(shù)最大,增強效果最好。綜合主觀視覺感受與客觀平均梯度指標(biāo),本文算法對工業(yè)低對比度X射線圖像的增強效果優(yōu)于所比較的算法。

    4 結(jié)語

    高范圍工業(yè)X射線圖像增強是一個比較新的研究課題。本文結(jié)合對數(shù)變換與梯度場,充分考慮圖像的整體與局部信息,組合成一個能夠應(yīng)用于低對比、低照度X工業(yè)射線增強的算法。實驗結(jié)果表明,本文算法不僅能突出圖像的細(xì)節(jié)信息,而且能讓高范圍工業(yè)X射線圖像符合人眼視覺系統(tǒng)的顯示在檢測屏幕設(shè)備上。但是本文算法需要根據(jù)不同情景進行調(diào)參,同時不能夠達到實時處理,即每秒30幀,無法直接用于X射線檢測設(shè)備的實時顯示上。后續(xù)的工作主要是自適應(yīng)參數(shù)選擇與提速,結(jié)合算法優(yōu)化、FPGA、GPU、SSE或CPU多線程并發(fā),達到實時處理顯示的目的。

    參考文獻(References)

    [1] 孫忠誠.射線數(shù)字成像技術(shù)[M]. 北京: 機械工業(yè)出版社, 2018: 122-124. (SUN Z C. X-ray Digital Imaging Technology[M]. Beijing: China Machine Press, 2018: 122-124.)

    [2] BURCH S F. Digital enhancement of video images for NDT[J]. NDT International, 1987, 20(1): 51-56.

    [3] PARIS S, HASINOFF S, KAUTZ J. Local Laplacian filters: edge-aware image processing with a Laplacian pyramid[J]. ACM Transactions on Graphics, 2011, 30(4): No.68.

    [4] GU K, ZHAI G, YANG X, et al. Automatic contrast enhancement technology with saliency preservation[J]. IEEE Transactions on Circuits & Systems for Video Technology, 2015, 25(9): 1480-1494.

    [5] DANIEL E, ANITHA J. Optimum green plane masking for the contrast enhancement of retinal images using enhanced genetic algorithm[J]. Optik — International Journal for Light and Electron Optics, 2015, 126(18): 1726-1730.

    [6] OSTOJIC V, STARCEVIC D, PETROVIC V. Artifact reduction in multiscale contrast enhancement for digital radiography[C]// Proceedings of the 22nd Telecommunications Forum Telfor. Piscataway: IEEE, 2014: 513-516.

    [7] LI M, LIU J, YANG W, et al. Structure-revealing low-light image enhancement via robust Retinex model[J]. IEEE Transactions on Image Processing, 2018, 27(6): 2828-2841.

    [8] 陳琛, 張建州. 改進的多尺度Retinex醫(yī)學(xué)X射線圖像增強算法[J]. 計算機工程與應(yīng)用, 2015, 51(9): 191-195. (CHEN C, ZHANG J C. Improved multi-scale Retinex enhancement algorithm of medical X ray images[J]. Computer Engineering & Applications, 2015, 51(9): 191-195.)

    [9] LI J, WANG Y, WANG L, et al. Industrial X-ray image enhancement algorithm based on adaptive histogram and wavelet[C]// Proceedings of the 6th International Forum on Strategic Technology. Piscataway: IEEE, 2011: 836-839.

    [10] TAN F, LI Y, GUAN M. Automatic windowing for highly dynamic industrial X-ray image based on short-term energy of gray histogram[J]. Journal of Nondestructive Evaluation, 2017, 36(3): 54-64.

    [11] 劉萬程. X射線圖像中目標(biāo)增強與檢測方法的研究[D]. 大連: 大連理工大學(xué), 2016: 16-28. (LIU W C. Research on target enhancement and detection based on X-ray image[D]. Dalian: Dalian University of Technology, 2016: 16-28.)

    [12] 韓永康. 高動態(tài)范圍X射線圖像融合方法研究[D]. 大連: 大連理工大學(xué), 2018: 20-35. (HAN Y K. Research on high dynamic range X-ray image fusion[D]. Dalian: Dalian University of Technology, 2018: 20-35.)

    [13] FATTAL R, LISCHINSKI D, WERMAN M. Gradient domain high dynamic range compression[C]// Proceedings of the 29th Conference on Computer Graphics & Interactive Techniques. New York: ACM, 2002: 249-256.

    [14] BRANCHITTA, FRANCESCO. New technique for the visualization of high dynamic range infrared images[J]. Optical Engineering, 2009, 48(9): No.096401.

    [15] DUAN J, DONG W, MU R, et al. Local contrast stretch based tone mapping for high dynamic range images[C]// Proceedings of the 2011 IEEE Symposium on Computational Intelligence for Multimedia, Signal and Vision Processing. Piscataway: IEEE, 2011: 26-32.

    [16] ZHANG F, XIE W, MA G, et al. High dynamic range compression and detail enhancement of infrared images in the gradient domain[J]. Infrared Physics & Technology, 2014, 67: 441-454.

    [17] KUNDU M K, PAL S K. Thresholding for edge detection using human psychovisual phenomena[J]. Pattern Recognition Letters, 1986, 4(6): 433-441.

    [18] PRESS W H, TEUKOLSKY S A, VETTERLING W T, et al. Numerical Recipes in C++: The Art of Scientific Computing[M]. Cambridge: Cambridge University Press, 2002: 39-44.

    [19] 梁傳暉, 王玉丹, 杜國浩, 等. 同步輻射X射線圖像對比度增強算法研究[J]. 光學(xué)學(xué)報, 2015, 35(3): 89-98. (LIANG C H, WANG Y D, DU G H, et al. Research on the contrast enhancement algorithm of synchrotron radiation X-ray image[J]. Acta Optica Sinica, 2015, 35(3): 89-98.)

    [20] 謝小甫, 周進, 吳欽章. 一種針對圖像模糊的無參考質(zhì)量評價指標(biāo)[J]. 計算機應(yīng)用, 2010, 30(4): 921-924. (XIE X F, ZHOU J, WU Q Z. No-reference quality index for image blur[J]. Journal of Computer Applications, 2010, 30(4): 921-924.)

    猜你喜歡
    圖像增強直方圖X射線
    統(tǒng)計頻率分布直方圖的備考全攻略
    “X射線”的那些事兒
    實驗室X射線管安全改造
    機電安全(2022年5期)2022-12-13 09:22:26
    符合差分隱私的流數(shù)據(jù)統(tǒng)計直方圖發(fā)布
    基于卷積神經(jīng)網(wǎng)絡(luò)的船舶圖像增強研究
    圖像增強技術(shù)在超跨聲葉柵紋影試驗中的應(yīng)用
    水下視覺SLAM圖像增強研究
    虛擬古生物學(xué):當(dāng)化石遇到X射線成像
    科學(xué)(2020年1期)2020-01-06 12:21:34
    用直方圖控制畫面影調(diào)
    基于圖像增強的無人機偵察圖像去霧方法
    精品日产1卡2卡| 999久久久精品免费观看国产| 99国产极品粉嫩在线观看| 99国产极品粉嫩在线观看| 宅男免费午夜| 免费人成视频x8x8入口观看| 精品午夜福利视频在线观看一区| 一区二区三区国产精品乱码| 婷婷亚洲欧美| 国产免费av片在线观看野外av| 国产黄色小视频在线观看| 精品无人区乱码1区二区| 我的亚洲天堂| 国产成人欧美| 婷婷六月久久综合丁香| 一边摸一边抽搐一进一小说| 国产片内射在线| 欧美绝顶高潮抽搐喷水| 国产精品久久久av美女十八| 精品国产乱码久久久久久男人| 黄色片一级片一级黄色片| 国产高清videossex| 757午夜福利合集在线观看| 亚洲真实伦在线观看| 亚洲av成人av| 欧美乱妇无乱码| 国产av一区在线观看免费| 女性被躁到高潮视频| 国产亚洲精品第一综合不卡| 久久久国产精品麻豆| 中文字幕人妻熟女乱码| 99国产精品99久久久久| 久久久久亚洲av毛片大全| 国产三级黄色录像| e午夜精品久久久久久久| 一级毛片精品| 久久精品人妻少妇| 91成年电影在线观看| 99在线视频只有这里精品首页| 正在播放国产对白刺激| 精品一区二区三区视频在线观看免费| 亚洲精品美女久久av网站| 十八禁人妻一区二区| 午夜成年电影在线免费观看| 日本成人三级电影网站| 亚洲男人天堂网一区| 99riav亚洲国产免费| 日本在线视频免费播放| 精品国产超薄肉色丝袜足j| 草草在线视频免费看| 久久久久久久精品吃奶| 人妻久久中文字幕网| 中文字幕精品免费在线观看视频| 亚洲国产毛片av蜜桃av| 一级片免费观看大全| 欧美另类亚洲清纯唯美| 一级黄色大片毛片| 国产精品一区二区三区四区久久 | 国产亚洲欧美精品永久| 少妇熟女aⅴ在线视频| 精品久久久久久,| 香蕉丝袜av| 国产亚洲欧美98| netflix在线观看网站| 神马国产精品三级电影在线观看 | 麻豆成人午夜福利视频| 人人妻人人看人人澡| 成人国产一区最新在线观看| 亚洲性夜色夜夜综合| 精品久久久久久久末码| 久久久久久亚洲精品国产蜜桃av| 欧美不卡视频在线免费观看 | 亚洲av成人av| 99久久无色码亚洲精品果冻| 精品欧美一区二区三区在线| 亚洲欧美日韩高清在线视频| 日本熟妇午夜| 妹子高潮喷水视频| 日韩免费av在线播放| 日韩精品青青久久久久久| 久久婷婷成人综合色麻豆| 国产野战对白在线观看| 在线观看一区二区三区| 欧美黑人欧美精品刺激| 免费看美女性在线毛片视频| 日韩大码丰满熟妇| 亚洲成a人片在线一区二区| 色av中文字幕| 午夜成年电影在线免费观看| 中文字幕高清在线视频| 免费无遮挡裸体视频| 欧美亚洲日本最大视频资源| www日本黄色视频网| 19禁男女啪啪无遮挡网站| 在线观看午夜福利视频| 久久久久国产一级毛片高清牌| 真人做人爱边吃奶动态| 国产成人系列免费观看| 国产成人精品久久二区二区91| 亚洲第一欧美日韩一区二区三区| 麻豆成人av在线观看| 国产欧美日韩一区二区三| 一本久久中文字幕| 成人精品一区二区免费| 精华霜和精华液先用哪个| 满18在线观看网站| 最新美女视频免费是黄的| 亚洲中文字幕一区二区三区有码在线看 | 一夜夜www| 国产成人系列免费观看| 日韩欧美三级三区| 最近最新中文字幕大全免费视频| a在线观看视频网站| 特大巨黑吊av在线直播 | 午夜影院日韩av| avwww免费| 久久久久久久午夜电影| 中国美女看黄片| 久久久久久亚洲精品国产蜜桃av| 国产精品一区二区免费欧美| 少妇熟女aⅴ在线视频| 午夜福利高清视频| 麻豆av在线久日| 日韩大码丰满熟妇| 午夜久久久在线观看| 国产精品电影一区二区三区| 亚洲国产精品合色在线| 熟女电影av网| 村上凉子中文字幕在线| 国产精品久久久久久精品电影 | 亚洲av片天天在线观看| 欧美日韩黄片免| 久久天堂一区二区三区四区| 制服丝袜大香蕉在线| 成人手机av| 91麻豆av在线| 国产99久久九九免费精品| 国产精品98久久久久久宅男小说| 久久久久亚洲av毛片大全| 精品乱码久久久久久99久播| 亚洲一码二码三码区别大吗| 国产真实乱freesex| 久久精品国产综合久久久| 国产又爽黄色视频| 后天国语完整版免费观看| 久久国产亚洲av麻豆专区| 日韩成人在线观看一区二区三区| 18禁裸乳无遮挡免费网站照片 | 91麻豆精品激情在线观看国产| 色精品久久人妻99蜜桃| 黄色毛片三级朝国网站| 啦啦啦免费观看视频1| 日韩大码丰满熟妇| 好看av亚洲va欧美ⅴa在| 少妇粗大呻吟视频| 深夜精品福利| 国产一区二区三区在线臀色熟女| 岛国在线观看网站| 午夜福利在线在线| 精品久久久久久久毛片微露脸| 一级作爱视频免费观看| 一级毛片精品| 男人舔奶头视频| 久久香蕉国产精品| 久久性视频一级片| 精品久久久久久久久久久久久 | 国产爱豆传媒在线观看 | 亚洲人成77777在线视频| 亚洲成av人片免费观看| 亚洲欧美日韩高清在线视频| 久久婷婷成人综合色麻豆| 黄片小视频在线播放| 日韩av在线大香蕉| 99精品久久久久人妻精品| 精品久久久久久久久久久久久 | 日本 av在线| 国产精品av久久久久免费| 麻豆成人午夜福利视频| 日韩欧美国产在线观看| 欧美日韩精品网址| 亚洲真实伦在线观看| 国产精品香港三级国产av潘金莲| 99精品欧美一区二区三区四区| 免费看十八禁软件| 国产精品久久电影中文字幕| 亚洲av成人不卡在线观看播放网| 国产精品久久久久久人妻精品电影| xxx96com| 精品国产一区二区三区四区第35| 三级毛片av免费| 岛国在线观看网站| 长腿黑丝高跟| 成年人黄色毛片网站| e午夜精品久久久久久久| 亚洲精品中文字幕一二三四区| 欧美人与性动交α欧美精品济南到| 午夜成年电影在线免费观看| 精品久久久久久久人妻蜜臀av| 日韩中文字幕欧美一区二区| 日本 av在线| 久久久久久久精品吃奶| 中出人妻视频一区二区| 91在线观看av| 男人舔奶头视频| av有码第一页| 99在线视频只有这里精品首页| 黄色丝袜av网址大全| 色老头精品视频在线观看| 久久人人精品亚洲av| 女性生殖器流出的白浆| 热99re8久久精品国产| 窝窝影院91人妻| 久久久久久九九精品二区国产 | 久久久久久国产a免费观看| 亚洲av电影不卡..在线观看| 97碰自拍视频| 欧美成人一区二区免费高清观看 | 一级毛片女人18水好多| 久久99热这里只有精品18| 国产精品美女特级片免费视频播放器 | 女警被强在线播放| 久久精品91蜜桃| 91av网站免费观看| 久久 成人 亚洲| 国产又色又爽无遮挡免费看| 99久久精品国产亚洲精品| 亚洲五月色婷婷综合| 亚洲无线在线观看| 精品久久久久久久久久免费视频| 一a级毛片在线观看| 日本一区二区免费在线视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲国产精品999在线| 成年人黄色毛片网站| www日本在线高清视频| АⅤ资源中文在线天堂| 亚洲第一电影网av| 一a级毛片在线观看| bbb黄色大片| 国产乱人伦免费视频| 国产精品爽爽va在线观看网站 | 久久国产乱子伦精品免费另类| 午夜免费激情av| 日本熟妇午夜| 老熟妇乱子伦视频在线观看| or卡值多少钱| 欧美大码av| 美女大奶头视频| 精品熟女少妇八av免费久了| 久久精品91蜜桃| 成人午夜高清在线视频 | 色综合站精品国产| 高潮久久久久久久久久久不卡| 久久国产精品男人的天堂亚洲| 99精品久久久久人妻精品| 一级a爱视频在线免费观看| 叶爱在线成人免费视频播放| 日日摸夜夜添夜夜添小说| 久久亚洲精品不卡| 国产精品亚洲av一区麻豆| 欧美亚洲日本最大视频资源| 亚洲五月婷婷丁香| 白带黄色成豆腐渣| 18禁美女被吸乳视频| 丁香六月欧美| 男女床上黄色一级片免费看| 制服人妻中文乱码| 国产真人三级小视频在线观看| 国产极品粉嫩免费观看在线| 妹子高潮喷水视频| 在线观看免费日韩欧美大片| 色综合婷婷激情| 欧美激情 高清一区二区三区| 色综合欧美亚洲国产小说| 国产爱豆传媒在线观看 | 国内精品久久久久精免费| 国产主播在线观看一区二区| 欧美日韩亚洲国产一区二区在线观看| 搡老岳熟女国产| 日韩有码中文字幕| 国产精品亚洲一级av第二区| 99久久无色码亚洲精品果冻| 成人国产综合亚洲| 亚洲av熟女| 热re99久久国产66热| 在线免费观看的www视频| 色av中文字幕| 欧美在线黄色| 91麻豆av在线| 欧美av亚洲av综合av国产av| av片东京热男人的天堂| 国产av在哪里看| 一本一本综合久久| av免费在线观看网站| 一夜夜www| 亚洲av片天天在线观看| a级毛片a级免费在线| 日韩大尺度精品在线看网址| 深夜精品福利| 夜夜躁狠狠躁天天躁| 他把我摸到了高潮在线观看| 法律面前人人平等表现在哪些方面| 午夜福利18| 久99久视频精品免费| 国产午夜福利久久久久久| 欧美成人一区二区免费高清观看 | 亚洲国产日韩欧美精品在线观看 | 亚洲av五月六月丁香网| 亚洲一码二码三码区别大吗| 欧美激情久久久久久爽电影| 男女之事视频高清在线观看| 日韩国内少妇激情av| 黄色视频,在线免费观看| 亚洲精品av麻豆狂野| 深夜精品福利| 高潮久久久久久久久久久不卡| 午夜免费观看网址| 午夜成年电影在线免费观看| 久久精品国产亚洲av香蕉五月| 日本撒尿小便嘘嘘汇集6| 一进一出抽搐gif免费好疼| 国产色视频综合| 精品不卡国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 国产免费av片在线观看野外av| 精品午夜福利视频在线观看一区| 日韩 欧美 亚洲 中文字幕| 性欧美人与动物交配| 亚洲真实伦在线观看| 午夜老司机福利片| 97碰自拍视频| 欧美色欧美亚洲另类二区| 老熟妇仑乱视频hdxx| 青草久久国产| 999精品在线视频| 18禁裸乳无遮挡免费网站照片 | 亚洲七黄色美女视频| 好看av亚洲va欧美ⅴa在| 露出奶头的视频| 99精品在免费线老司机午夜| 国产av不卡久久| 夜夜躁狠狠躁天天躁| 亚洲av美国av| 91老司机精品| 免费高清视频大片| 窝窝影院91人妻| 国产av不卡久久| 日韩大尺度精品在线看网址| 亚洲久久久国产精品| 91大片在线观看| 国产一卡二卡三卡精品| 制服人妻中文乱码| 波多野结衣高清作品| 国产aⅴ精品一区二区三区波| 国产精品久久久人人做人人爽| 欧美日韩亚洲国产一区二区在线观看| 一区二区日韩欧美中文字幕| 午夜成年电影在线免费观看| 国产91精品成人一区二区三区| 伦理电影免费视频| 88av欧美| 看免费av毛片| 亚洲av日韩精品久久久久久密| 中文字幕av电影在线播放| 亚洲精品国产精品久久久不卡| 中文字幕久久专区| 欧美日韩中文字幕国产精品一区二区三区| 波多野结衣巨乳人妻| 19禁男女啪啪无遮挡网站| 亚洲 欧美一区二区三区| 国产精品二区激情视频| 丁香六月欧美| 亚洲性夜色夜夜综合| 搡老熟女国产l中国老女人| 叶爱在线成人免费视频播放| 757午夜福利合集在线观看| 国产成人精品久久二区二区91| 亚洲一区中文字幕在线| 母亲3免费完整高清在线观看| 美女国产高潮福利片在线看| 中国美女看黄片| 一边摸一边做爽爽视频免费| tocl精华| 999久久久国产精品视频| 欧美人与性动交α欧美精品济南到| 成人av一区二区三区在线看| 国产麻豆成人av免费视频| 精品电影一区二区在线| 桃色一区二区三区在线观看| 美女午夜性视频免费| 国产成人欧美在线观看| av欧美777| www.精华液| 久久精品夜夜夜夜夜久久蜜豆 | 丁香六月欧美| 叶爱在线成人免费视频播放| 国产精品国产高清国产av| 欧美zozozo另类| 成人国产综合亚洲| 国产成人影院久久av| 亚洲三区欧美一区| 久99久视频精品免费| 久久欧美精品欧美久久欧美| 国产单亲对白刺激| 婷婷六月久久综合丁香| 精品久久久久久久久久免费视频| 久久精品影院6| 国产精品电影一区二区三区| 俄罗斯特黄特色一大片| 欧美精品啪啪一区二区三区| a级毛片在线看网站| 国产亚洲av高清不卡| 十八禁人妻一区二区| 亚洲精品中文字幕在线视频| 在线观看免费午夜福利视频| 色在线成人网| 国产精品 欧美亚洲| 18禁国产床啪视频网站| 日韩大尺度精品在线看网址| 日韩视频一区二区在线观看| 欧美性长视频在线观看| 在线观看午夜福利视频| 少妇粗大呻吟视频| 久久久久久亚洲精品国产蜜桃av| av欧美777| 91大片在线观看| 天堂√8在线中文| 色综合婷婷激情| 日本免费a在线| 久久久久国产精品人妻aⅴ院| 久久精品91无色码中文字幕| 男女午夜视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产高清激情床上av| 九色国产91popny在线| 在线天堂中文资源库| 窝窝影院91人妻| 国产精品精品国产色婷婷| 亚洲专区字幕在线| 久久精品91无色码中文字幕| 亚洲国产欧美日韩在线播放| 亚洲人成伊人成综合网2020| 国产成人系列免费观看| 日韩欧美国产在线观看| 国产精品九九99| 国产伦一二天堂av在线观看| 久久天躁狠狠躁夜夜2o2o| 手机成人av网站| 国产成人影院久久av| av电影中文网址| 嫩草影院精品99| 欧美精品啪啪一区二区三区| 一区二区三区精品91| 亚洲成人久久性| 非洲黑人性xxxx精品又粗又长| 日本五十路高清| 中出人妻视频一区二区| 国产精品,欧美在线| 精品国内亚洲2022精品成人| 欧美激情高清一区二区三区| 欧美精品亚洲一区二区| 国产真实乱freesex| 国产三级在线视频| 午夜久久久在线观看| 男女床上黄色一级片免费看| 91九色精品人成在线观看| 免费人成视频x8x8入口观看| 亚洲精品中文字幕在线视频| 最新在线观看一区二区三区| 国产免费av片在线观看野外av| 一级毛片女人18水好多| 国产亚洲精品久久久久5区| 黄色成人免费大全| 最近在线观看免费完整版| 波多野结衣高清无吗| 国产亚洲欧美精品永久| 亚洲av成人av| 亚洲自偷自拍图片 自拍| 免费看十八禁软件| 97碰自拍视频| 日韩欧美在线二视频| 亚洲国产欧洲综合997久久, | 夜夜爽天天搞| 日韩欧美 国产精品| 女人被狂操c到高潮| 黄色女人牲交| 脱女人内裤的视频| 免费观看人在逋| 婷婷精品国产亚洲av| 久久香蕉激情| 美女午夜性视频免费| 国产精品一区二区免费欧美| 丝袜美腿诱惑在线| 国产精品免费视频内射| 无人区码免费观看不卡| 美女免费视频网站| xxx96com| 欧美激情高清一区二区三区| 日韩中文字幕欧美一区二区| 亚洲成a人片在线一区二区| 国产av一区二区精品久久| 精品国产乱子伦一区二区三区| 亚洲九九香蕉| 免费无遮挡裸体视频| 麻豆成人午夜福利视频| 欧美在线黄色| 亚洲 欧美一区二区三区| 亚洲av第一区精品v没综合| 精品国产亚洲在线| 性欧美人与动物交配| 久久中文看片网| 欧美日本视频| 一本久久中文字幕| 婷婷精品国产亚洲av在线| 成人18禁在线播放| 91成年电影在线观看| 777久久人妻少妇嫩草av网站| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产毛片av蜜桃av| 欧美一级a爱片免费观看看 | 狂野欧美激情性xxxx| 1024香蕉在线观看| 女人爽到高潮嗷嗷叫在线视频| 久久天躁狠狠躁夜夜2o2o| 伊人久久大香线蕉亚洲五| 欧美中文综合在线视频| 久久久水蜜桃国产精品网| 99国产精品一区二区三区| 久久午夜亚洲精品久久| 精华霜和精华液先用哪个| 人妻丰满熟妇av一区二区三区| 国产黄色小视频在线观看| 精品一区二区三区视频在线观看免费| 女人被狂操c到高潮| 88av欧美| 在线十欧美十亚洲十日本专区| 国产色视频综合| 亚洲国产精品999在线| 午夜免费鲁丝| 久久久久久久精品吃奶| 大型黄色视频在线免费观看| 国产激情久久老熟女| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| 别揉我奶头~嗯~啊~动态视频| 亚洲国产看品久久| 一级毛片精品| 久久精品91无色码中文字幕| 午夜免费成人在线视频| 最近最新中文字幕大全电影3 | 波多野结衣巨乳人妻| 最近最新免费中文字幕在线| 97碰自拍视频| 国产又黄又爽又无遮挡在线| 免费观看人在逋| 黄色 视频免费看| 亚洲av电影不卡..在线观看| 国产亚洲欧美在线一区二区| 久久国产精品影院| 桃色一区二区三区在线观看| 男女床上黄色一级片免费看| 又黄又爽又免费观看的视频| 黄色片一级片一级黄色片| 国产精品久久久久久精品电影 | 又大又爽又粗| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全电影3 | 亚洲av片天天在线观看| a在线观看视频网站| 亚洲欧美日韩无卡精品| 啦啦啦韩国在线观看视频| 亚洲七黄色美女视频| 日本免费一区二区三区高清不卡| 欧美中文综合在线视频| 国产精品自产拍在线观看55亚洲| 精品国产超薄肉色丝袜足j| 色综合欧美亚洲国产小说| 成人免费观看视频高清| 久久香蕉精品热| 老司机午夜福利在线观看视频| 国产午夜精品久久久久久| 亚洲欧美精品综合一区二区三区| 久久这里只有精品19| 欧美午夜高清在线| 成人永久免费在线观看视频| 久久九九热精品免费| 一a级毛片在线观看| 精品卡一卡二卡四卡免费| 91国产中文字幕| 一a级毛片在线观看| 国产亚洲精品综合一区在线观看 | 我的亚洲天堂| 一个人观看的视频www高清免费观看 | 免费看日本二区| 欧美激情极品国产一区二区三区| 一级黄色大片毛片| 欧美中文日本在线观看视频| 久久香蕉国产精品| 熟女电影av网| 校园春色视频在线观看| 天天添夜夜摸| 丝袜在线中文字幕| 国产av一区在线观看免费| 91av网站免费观看| 日韩免费av在线播放| 91大片在线观看| 国产成人精品无人区| ponron亚洲| 最近最新中文字幕大全电影3 | 国产伦一二天堂av在线观看| 一级作爱视频免费观看| 91麻豆精品激情在线观看国产| 亚洲最大成人中文| 国产国语露脸激情在线看| 亚洲激情在线av|