《電子技術(shù)應(yīng)用》
您所在的位置:首頁(yè) > 顯示光電 > 設(shè)計(jì)應(yīng)用 > 基于OpenCV和LSSVM的數(shù)字儀表讀數(shù)自動(dòng)識(shí)別
基于OpenCV和LSSVM的數(shù)字儀表讀數(shù)自動(dòng)識(shí)別
2017年微型機(jī)與應(yīng)用第2期
林劍萍1,廖一鵬2
1. 陽(yáng)光學(xué)院,福建 福州 350015;2. 福州大學(xué) 物理與信息工程學(xué)院,福建 福州 350000
摘要: 為了提高儀表讀數(shù)及記錄的工作效率,減少人工干預(yù)并降低出錯(cuò)概率,提出了一種基于OpenCV和LSSVM的數(shù)字儀表自動(dòng)識(shí)別方法。首先選取免驅(qū)的USB攝像頭采集儀表圖像,然后采用Canny檢測(cè)與Hough變換相結(jié)合的算法對(duì)儀表圖像的傾斜角度進(jìn)行校正;利用數(shù)學(xué)形態(tài)學(xué)、Otsu算法以及連通域提取對(duì)圖像進(jìn)行預(yù)處理和字符定位分割;最后,進(jìn)行字符特征提取并利用最小二乘支持向量機(jī)算法對(duì)儀表讀數(shù)進(jìn)行識(shí)別。經(jīng)過(guò)大量的實(shí)驗(yàn)驗(yàn)證,與現(xiàn)有的儀表讀數(shù)識(shí)別方法相比,該方法具有更好的準(zhǔn)確性,實(shí)現(xiàn)了數(shù)字儀表的自動(dòng)識(shí)別功能。
Abstract:
Key words :

  林劍萍1,廖一鵬2

  (1. 陽(yáng)光學(xué)院,福建 福州 350015;2. 福州大學(xué) 物理與信息工程學(xué)院,福建 福州 350000)

       摘要:為了提高儀表讀數(shù)及記錄的工作效率,減少人工干預(yù)并降低出錯(cuò)概率,提出了一種基于OpenCV和LSSVM的數(shù)字儀表自動(dòng)識(shí)別方法。首先選取免驅(qū)的USB攝像頭采集儀表圖像,然后采用Canny檢測(cè)與Hough變換相結(jié)合的算法對(duì)儀表圖像的傾斜角度進(jìn)行校正;利用數(shù)學(xué)形態(tài)學(xué)、Otsu算法以及連通域提取對(duì)圖像進(jìn)行預(yù)處理和字符定位分割;最后,進(jìn)行字符特征提取并利用最小二乘支持向量機(jī)算法對(duì)儀表讀數(shù)進(jìn)行識(shí)別。經(jīng)過(guò)大量的實(shí)驗(yàn)驗(yàn)證,與現(xiàn)有的儀表讀數(shù)識(shí)別方法相比,該方法具有更好的準(zhǔn)確性,實(shí)現(xiàn)了數(shù)字儀表的自動(dòng)識(shí)別功能。

  關(guān)鍵詞機(jī)器視覺(jué)圖像處理;OpenCV;最小二乘支持向量機(jī)

  中圖分類號(hào):TP394.1文獻(xiàn)標(biāo)識(shí)碼:ADOI: 10.19358/j.issn.1674-7720.2017.02.013

  引用格式:林劍萍,廖一鵬.基于OpenCV和LSSVM的數(shù)字儀表讀數(shù)自動(dòng)識(shí)別[J].微型機(jī)與應(yīng)用,2017,36(2):37-40.

0引言

  目前,在很多領(lǐng)域受工作條件限制及經(jīng)濟(jì)性、便捷性的要求,很多用于計(jì)量的儀器儀表都沒(méi)有專門的數(shù)據(jù)通信接口,導(dǎo)致無(wú)法自動(dòng)識(shí)別讀數(shù),對(duì)儀表數(shù)字的讀取還是要靠人工來(lái)實(shí)現(xiàn),但是在一些特定場(chǎng)合,例如高溫高壓、化工冶金、高山峭壁、有核輻射等人體不能適應(yīng)的地方,人工抄表存在很大的危險(xiǎn)性及不可行性,而機(jī)器視覺(jué)卻可以無(wú)視這些不利因素,廣泛地用于長(zhǎng)時(shí)間惡劣的工作環(huán)境,這就為儀表讀數(shù)的讀取提供了便利,同時(shí)可以提高工作效率,更好地保障人身安全。儀表讀數(shù)識(shí)別也可以說(shuō)是自動(dòng)抄表,自動(dòng)抄表的試點(diǎn)和運(yùn)用始于20世紀(jì)80年代,在20世紀(jì)90年代得到較快的發(fā)展[1]。目前,國(guó)內(nèi)外許多科技人員在基于機(jī)器視覺(jué)的數(shù)字識(shí)別上做了很多深入研究。從圖像預(yù)處理、傾斜校正、字符分割及數(shù)字識(shí)別等方面各有側(cè)重點(diǎn)地做了不同方法的研究比較。在儀表數(shù)字識(shí)別方面,工業(yè)和信息化部電子第五研究所的申中鴻、蔣春旭等人利用機(jī)器視覺(jué)技術(shù),結(jié)合計(jì)算幾何軌跡識(shí)別方法與D-S證據(jù)理論實(shí)現(xiàn)數(shù)字顯示儀表的自動(dòng)讀數(shù);中國(guó)計(jì)量科學(xué)研究院的劉科等人采用投影變換的方法對(duì)數(shù)字圖像進(jìn)行分割、編碼,進(jìn)而實(shí)現(xiàn)了對(duì)該類儀表示值的自動(dòng)讀取;山東管理學(xué)院的常曉瑋利用光學(xué)字符識(shí)別技術(shù)實(shí)現(xiàn)了遠(yuǎn)程識(shí)別儀表上顯示的數(shù)字信息[25]。綜合各方研究,實(shí)現(xiàn)基于機(jī)器視覺(jué)的數(shù)字儀表讀數(shù)自動(dòng)識(shí)別的方式多樣,不僅在圖像預(yù)處理上可根據(jù)不同環(huán)境不同儀表選擇合適的算法,在數(shù)字的定位識(shí)別上也有多種方法可供選擇。本文采用Canny檢測(cè)與Hough變換相結(jié)合的算法對(duì)攝像頭采集到的儀表圖像的傾斜角度進(jìn)行校正,利用數(shù)學(xué)形態(tài)學(xué)方法、Otsu(大津法或最大類間方差法)算法等技術(shù)對(duì)校正后的圖像進(jìn)行預(yù)處理[67],之后對(duì)字符進(jìn)行分割,采用最小二乘支持向量機(jī)(LSSVM)算法對(duì)儀表讀數(shù)進(jìn)行識(shí)別。經(jīng)過(guò)大量的實(shí)驗(yàn)驗(yàn)證,該方法具有良好的準(zhǔn)確性,實(shí)現(xiàn)了數(shù)字儀表的自動(dòng)識(shí)別功能。

1數(shù)字儀表圖像的預(yù)處理

  1.1圖像的采集

  圖像采集是基于機(jī)器視覺(jué)進(jìn)行圖像處理的第一步,也是很重要的一步。采集到質(zhì)量較高的圖像可以在很大程度上降低處理的難度。本文采用分辨率為480×680的CCD免驅(qū)USB攝像頭對(duì)儀表數(shù)據(jù)進(jìn)行采集,該攝像機(jī)可以直接捕捉影像傳輸處理,運(yùn)用更加便捷。

  1.2儀表圖像的傾斜校正

  進(jìn)行數(shù)字儀表圖像采集時(shí),采集到的數(shù)字儀表圖像經(jīng)常會(huì)發(fā)生傾斜。為使儀表端正以方便后期處理,需要對(duì)傾斜的儀表圖像進(jìn)行校正。本文采用的是Canny邊緣檢測(cè)與Hough變換相結(jié)合的傾斜校正方法。用Canny算子對(duì)儀表圖像進(jìn)行邊緣檢測(cè),然后利用Hough變換檢測(cè)數(shù)字儀表邊框平行線的傾斜角度,再利用旋轉(zhuǎn)函數(shù)進(jìn)行校正,得到端正水平的數(shù)字儀表圖像。

001.jpg

  攝像頭采集到的圖像如圖1所示,Canny邊緣檢測(cè)及Hough變換檢測(cè)得到最長(zhǎng)直線圖如圖2所示,傾斜校正后的圖像如圖3所示。

002.jpg

  1.3圖像的形態(tài)學(xué)處理

  為簡(jiǎn)化圖像數(shù)據(jù),并除去不相干的結(jié)構(gòu),本系統(tǒng)對(duì)圖像進(jìn)行腐蝕(Erosion)膨脹(Dilation)等數(shù)學(xué)形態(tài)學(xué)運(yùn)算,對(duì)圖像進(jìn)行濾波處理,平滑邊界。提高數(shù)字儀表圖像的質(zhì)量,并填充儀表讀數(shù)中的細(xì)小空洞,避免單個(gè)數(shù)字各段碼之間的中斷。

  膨脹的定義:

  D=X⊕B={(x,y)|Bxy∩X≠φ}(1)

  腐蝕的定義:

  E=XB={(x,y)|BxyX}(2)

  式(1)和(2)中,B為結(jié)構(gòu)元素,X為原始圖像。

  1.4圖像的二值化

  為了字符的準(zhǔn)確識(shí)別,在字符識(shí)別之前對(duì)圖像做二值化處理,將圖像分割成背景和目標(biāo)兩部分。本文采用最經(jīng)典的全局閾值的Otsu算法(最大類間方差法)。為了減少對(duì)比度和光照不均勻?qū)Χ祷挠绊懀诙祷幚碇跋葘?duì)圖像進(jìn)行頂帽處理,讓二值化達(dá)到更好的效果。Otsu算法是在最小二乘法原理基礎(chǔ)上推導(dǎo)出來(lái)的,它的基本原理是:令連通區(qū)域S,其中像素灰度范圍為[0,255],點(diǎn)(x,y)的灰度表示為I(x,y),灰度級(jí)x的點(diǎn)的概率為P(x),則Otsu就是求下式達(dá)到最大值的灰度級(jí)m:

  PS6~3(_HM}(F0[IKY~@93LD.png

  

003.jpg

  灰度級(jí)m即為最佳分割閾值。利用此方法二值化處理后的圖像如圖4所示。

2數(shù)字儀表圖像的特征提取

  經(jīng)過(guò)預(yù)處理之后,可以得到質(zhì)量較高的二值化圖像,接下來(lái)就需要對(duì)數(shù)字儀表的讀數(shù)進(jìn)行定位與分割,為數(shù)字識(shí)別做好準(zhǔn)備。定位分割的準(zhǔn)確性,直接影響數(shù)字識(shí)別的難易程序及識(shí)別率。

  2.1儀表數(shù)字區(qū)域的特征分析

  在數(shù)字儀表讀數(shù)的定位分割過(guò)程中,必須根據(jù)數(shù)字區(qū)域與背景儀表的不同特征進(jìn)行提取,才能達(dá)到定位分割的目的。而數(shù)字儀表的特征總體是有跡可循的,可以將數(shù)字區(qū)域與背景分別出來(lái)。從視覺(jué)角度出發(fā),儀表數(shù)字區(qū)域的特征可以從其形狀、顏色和亮度三個(gè)方面進(jìn)行分析。

  2.2儀表數(shù)字定位分割方法

  儀表數(shù)字定位分割一般有基于連通域和基于投影法的儀表數(shù)字定位與分割。但是不同的光照條件,儀表的影子對(duì)水平投影的影響過(guò)大,在有影子存在的情況下,僅根據(jù)水平投影不足以將兩行讀數(shù)分割開(kāi)來(lái),而水平投影結(jié)果不理想導(dǎo)致分割不準(zhǔn)確,也就使得垂直投影不能理想地分割出每行中的每一個(gè)字符。因此本文采用的是基于連通域的方法來(lái)進(jìn)行數(shù)字的定位與分割。

  在圖像當(dāng)中用f(x,y)表示位于圖像陣列中第x行、第y列的像素的值,一幅m×n的圖像具有m行n列,f(0,0)表示圖像左上角的像素值,f(m-1,n-1)表示圖像右下角的像素值。遍歷圖像,通過(guò)4鄰域連通標(biāo)準(zhǔn)查找出所有符合的連通域,這樣就可以把圖像中每一個(gè)分離的區(qū)域檢測(cè)出來(lái)。

  在查找出所有連通域之后,按照連通域的面積、讀數(shù)的形狀、長(zhǎng)寬比來(lái)進(jìn)行篩選,排除不符合條件的區(qū)域,這樣就可以直接定位并分割出數(shù)字儀表的每個(gè)讀數(shù)。算法實(shí)現(xiàn)步驟如下。

  (1)計(jì)算連通域面積:連通域A的面積就是像素值為255的點(diǎn)的數(shù)目,即區(qū)域的邊界內(nèi)包含的像素點(diǎn)數(shù)。面積的計(jì)算公式為:

  5YYA`{8EMSQUC$2]OEUP$`L.png

  (2)計(jì)算連通域的長(zhǎng)寬比:連通域的寬/連通域的長(zhǎng)。

  (3)遍歷所有連通域,如果連通域面積大于500小于5 000且連通域的長(zhǎng)寬比大于1.5小于2.5則提取,如果不滿足以上條件則舍棄。這樣便可篩選出滿足要求的連通域,即數(shù)字儀表每個(gè)讀數(shù)。圖5儀表讀數(shù)分割定位結(jié)果

  數(shù)字儀表讀數(shù)的定位分割結(jié)果如圖5所示。

004.jpg

  2.3儀表數(shù)字特征提取

  數(shù)字的特征主要可以從形狀、顏色、紋理、結(jié)構(gòu)等方面入手來(lái)提取,而大部分的數(shù)字儀表的讀數(shù)區(qū)域是基于七段數(shù)碼管來(lái)實(shí)現(xiàn)顯示,因此也可以根據(jù)數(shù)碼管的特征來(lái)提取。本文用的提取特征的方法就是基于七段數(shù)碼管,根據(jù)每段碼的像素值來(lái)判斷識(shí)別。將定位好的數(shù)碼管每位進(jìn)行4×4分割,統(tǒng)計(jì)黑點(diǎn)和白點(diǎn)的個(gè)數(shù),每塊的像素用x0~x15表示,提取結(jié)果如圖6所示。

005.jpg

 3儀表讀數(shù)的模式識(shí)別

  數(shù)字儀表的讀數(shù)識(shí)別是基于機(jī)器視覺(jué)的數(shù)字儀表讀數(shù)自動(dòng)識(shí)別的重點(diǎn),只有這一步做好了才能正確識(shí)別讀數(shù)并且保證一定的識(shí)別率。最常用的數(shù)字字符識(shí)別算法主要有模板匹配法、基于人工神經(jīng)網(wǎng)絡(luò)法、穿線法等[810]。本文采用的是最小二乘支持向量機(jī)(LSSVM)算法[1113]。

  設(shè)訓(xùn)練集為(xi,yi)(i=1,2,…,N),N表示訓(xùn)練樣本數(shù),xi∈Rm為樣本輸入,yi為輸出,LSSVM在高維特征空間的線性函數(shù)為:

  f(x)=wTφ(x)+b(5)

  式(5)中,w為權(quán)值向量,b為偏置量,φ(x)為非線性映射函數(shù)。根據(jù)結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則,綜合考慮擬合誤差和算法復(fù)雜度,轉(zhuǎn)換為一個(gè)等式約束的優(yōu)化問(wèn)題:

  I2~~CF]Y`10}WP)9PY7$IJK.png

  式(6)中,e2i為誤差變量,d為調(diào)整參數(shù)因子。

  為求解上述優(yōu)化問(wèn)題,引入拉格朗日等式,ai是拉格朗日乘子,把約束優(yōu)化問(wèn)題變?yōu)闊o(wú)約束優(yōu)化問(wèn)題,如式(7)所示:

  CTYKZLJ)1V4_G)UDA}[050X.png

  根據(jù)非線性最優(yōu)規(guī)劃(KarushKuhnTucker,KKT)條件求得參數(shù)ai,b代入公式(5),得到支持向量機(jī)的輸出:

  1IE`~){2H3SKKG0GYDBS40P.png

  最后,選擇RBF(Radial Basis Function)作為L(zhǎng)SSVM的核函數(shù),得到LSSVM的預(yù)測(cè)函數(shù):

  FZ9RLBXU$3_[6KDWK1_VF~Y.png

  式中,σ為核寬度,反映數(shù)據(jù)集的特性。

4系統(tǒng)軟件實(shí)現(xiàn)與測(cè)試

  4.1系統(tǒng)軟件實(shí)現(xiàn)

  本系統(tǒng)利用OpenCV開(kāi)源庫(kù)[1415]在VS2010環(huán)境中編寫(xiě)MFC程序,實(shí)現(xiàn)數(shù)字儀表讀數(shù)的實(shí)時(shí)采集、處理及識(shí)別輸出。通過(guò)點(diǎn)擊界面內(nèi)按鈕即可完成相應(yīng)的處理,最終將數(shù)字儀表原圖及讀數(shù)識(shí)別結(jié)果都直觀地顯示在窗口界面上。

  根據(jù)本文系統(tǒng)采取的方法對(duì)實(shí)時(shí)采集到的數(shù)字儀表圖像進(jìn)行處理、分割及識(shí)別,以采集到的其中一張測(cè)試圖為例,系統(tǒng)從實(shí)時(shí)采集、傾斜校正、預(yù)處理(濾波、二值化)、定位分割到識(shí)別的每一步效果圖如圖7所示。

 

006.jpg

  4.2LSSVM與其他算法的比較

  本文采用最小二乘支持向量機(jī)算法對(duì)數(shù)字儀表的讀數(shù)進(jìn)行識(shí)別,為了驗(yàn)證該算法的識(shí)別率,總共實(shí)時(shí)采集了200幅儀表圖像近800個(gè)數(shù)字字符,用模板匹配法、基于人工神經(jīng)網(wǎng)絡(luò)法、穿線法、最小二乘支持向量機(jī)分別進(jìn)行了測(cè)試,對(duì)這800個(gè)數(shù)字字符進(jìn)行分類統(tǒng)計(jì),測(cè)試結(jié)果如表1所示。

  統(tǒng)計(jì)分析發(fā)現(xiàn),這幾種算法對(duì)每個(gè)數(shù)字的識(shí)別率都高于90.0%,LSSVM對(duì)每個(gè)數(shù)字的識(shí)別率都高于95.0%,對(duì)10個(gè)字符的平均識(shí)別率達(dá)到98.2%,比其他幾種識(shí)別模式的平均識(shí)別率有明顯提高,識(shí)別效果良好,達(dá)到了預(yù)期的目標(biāo)與要求。

5結(jié)論

  根據(jù)數(shù)字儀表顯示的特征,本文首先采用攝像頭實(shí)時(shí)采集并顯示數(shù)字儀表圖像;然后將Canny邊緣檢測(cè)算子與Hough變換相結(jié)合,對(duì)儀表圖像進(jìn)行傾斜校正;利用數(shù)學(xué)形態(tài)學(xué)方法對(duì)圖像光照不均勻等進(jìn)行處理,通過(guò)濾波增強(qiáng)有效信息;根據(jù)數(shù)字儀表讀數(shù)特征,采取查找最小連通域后篩選的方法對(duì)儀表讀數(shù)進(jìn)行定位分割;采用和聲搜索法和最小二乘支持向量機(jī)的算法對(duì)儀表的數(shù)字讀數(shù)進(jìn)行模式識(shí)別;利用OpenCV開(kāi)源庫(kù)在VS2010環(huán)境中編寫(xiě)MFC程序,實(shí)現(xiàn)數(shù)字儀表讀數(shù)的實(shí)時(shí)采集、處理及識(shí)別輸出。本設(shè)計(jì)在實(shí)驗(yàn)室環(huán)境下具有良好的效果,下一步的研究方向是進(jìn)一步改進(jìn)完善,并將其推廣應(yīng)用到實(shí)際環(huán)境中。  

007.jpg

參考文獻(xiàn)

  [1] 譚志強(qiáng),黃懿.自動(dòng)抄表技術(shù)的發(fā)展[J].電測(cè)與儀表, 2009,46(1):1-5.

  [2] 申中鴻,蔣春旭,徐華偉,等.基于機(jī)器視覺(jué)的數(shù)字顯示儀表自動(dòng)讀數(shù)[J].電子產(chǎn)品可靠性與環(huán)境試驗(yàn),2013,31(A01):110-115.

  [3] 申中鴻,劉群興,蔣春旭,等.基于機(jī)器視覺(jué)技術(shù)和D—S證據(jù)理論的數(shù)字顯示儀表自動(dòng)讀數(shù)[J].信息工程期刊(中英文版),2013,3(6):122-129.

  [4] 劉科,侯立新,卞昕.基于機(jī)器視覺(jué)的儀表示值識(shí)別算法研究[J].計(jì)量學(xué)報(bào),2013,34(5):425429.

  [5] 常曉瑋.基于機(jī)器視覺(jué)的數(shù)字化儀表字符識(shí)別[J].計(jì)算機(jī)測(cè)量與控制,2014,22(9):30713073.[6] 丁文,李勃,陳啟美.Canny色彩邊緣檢測(cè)及分類新方法[J].北京郵電大學(xué)學(xué)報(bào),2012,35(1):115119.[7] 黃登山,王曉楠.基于Canny和HillbertHuang變換的改進(jìn)的SIFT算法研究[J].西北工業(yè)大學(xué)學(xué)報(bào),2015(6):962-965.

  [8] 朱顥東,李紅嬋.基于特征加權(quán)模糊模板匹配的字符識(shí)別[J].蘭州理工大學(xué)學(xué)報(bào),2013,39(1) :68-71.

  [9] 劉宗林,田瑩,李雪瀅.一種數(shù)顯儀表字符快速自動(dòng)識(shí)別方

  法[J].現(xiàn)代計(jì)算機(jī)(專業(yè)版),2015(18):37-40.

  [10] 王明勛.基于神經(jīng)網(wǎng)絡(luò)的液晶數(shù)字識(shí)別在船用設(shè)備管理中的應(yīng)用[J].艦船電子工程,2012,32(7):95-97.

  [11] 劉永平,郭小波.基于最小二乘支持向量機(jī)車牌字符特征識(shí)別[J].數(shù)字技術(shù)與應(yīng)用,2015(7):119-120.

  [12] 劉靜.基于最小二乘支持向量機(jī)的車牌字符特征分類研究[J].計(jì)算機(jī)與數(shù)字工程,2015(7):1315-1319.

  [13] 喬宗良,張蕾,周建新,等.一種改進(jìn)的CPSOLSSVM軟測(cè)量模型及其應(yīng)用[J].儀器儀表學(xué)報(bào),2014,35(1):234-240.

  [14] 王舒憬,杜皎潔.基于OPENCV的數(shù)字萬(wàn)用表數(shù)字識(shí)別方法[J].自動(dòng)化與儀器儀表, 2014(1):62-63.

  [15]  徐歡,李振璧,姜媛媛,等.基于OpenCV和改進(jìn)Canny算子的路面裂縫檢測(cè)[J].計(jì)算機(jī)工程與設(shè)計(jì),2014,35(12):4254-4258.


此內(nèi)容為AET網(wǎng)站原創(chuàng),未經(jīng)授權(quán)禁止轉(zhuǎn)載。
欧美激情办公室aⅴ_国产欧美综合一区二区三区_欧美午夜精品久久久久免费视_福利视频欧美一区二区三区

          国产在线精品二区| 国产精品国产三级欧美二区| 久久久久久久波多野高潮日日| 91久久久久| 国产欧美日韩一区二区三区| 国产精品一国产精品k频道56| 一本一本久久a久久精品综合妖精| 亚洲精品欧洲精品| 国产精品免费一区二区三区在线观看| 校园激情久久| 国产精品初高中精品久久| 亚洲午夜精品福利| 日韩视频在线一区二区三区| 国产精品一区二区三区四区五区| 欧美资源在线| 激情偷拍久久| 香蕉亚洲视频| 亚洲视频一二| 久久av一区二区三区| 亚洲视频高清| 麻豆久久久9性大片| 亚洲天堂黄色| 毛片一区二区| 日韩天天综合| 国产精品v亚洲精品v日韩精品 | 国产欧美日韩综合一区在线观看 | 国产人成精品一区二区三| 亚洲欧美日韩视频二区| 国产精品mm| 美女久久一区| 国产精品老牛| 亚洲欧洲精品一区二区三区波多野1战4 | 久久av在线| 精品白丝av| 欧美 日韩 国产一区二区在线视频| 在线观看一区| 国产精品videosex极品| 久久午夜精品一区二区| 在线综合亚洲| 亚洲经典三级| 亚洲国产日韩欧美一区二区三区| 亚洲欧美文学| 久久一区二区精品| 午夜在线观看免费一区| 日韩一区二区久久| 激情欧美日韩一区| 伊人久久综合| 狠狠色综合网| 伊人久久婷婷| 在线不卡欧美| 最新日韩在线| 一区二区三区国产盗摄| 在线午夜精品| 国产日韩精品久久| 国产亚洲福利| 亚洲欧美日韩精品综合在线观看| 亚洲视频二区| 亚洲永久免费精品| 西西人体一区二区| 久久国产高清| 欧美久久综合| 黄色av成人| 亚洲精品影视| 性欧美xxxx大乳国产app| 西西裸体人体做爰大胆久久久| 亚洲免费在线精品一区| 久久精品网址| 欧美婷婷在线| 91久久黄色| 中文在线一区| 欧美一区精品| 亚洲欧洲午夜| 美日韩精品免费| 欧美深夜福利| 国产一区二区三区的电影 | 欧美日韩国内| 一区二区视频欧美| 国产欧美日韩综合一区在线播放 | 午夜精品亚洲| 韩国一区二区三区美女美女秀| 亚洲日本精品国产第一区| 中国成人亚色综合网站| 亚洲一区二区三区精品在线观看| 久久久久看片| 欧美日韩99| 99综合精品| 欧美另类亚洲| 国产精品久久久一区二区三区 | 国产精品视频免费一区| 欧美在线一二三区| 在线日韩av| 狂野欧美一区| 一区二区三区精品视频在线观看| 久久视频一区| 亚洲精品123区| 欧美ab在线视频| 夜夜夜久久久| 国内外成人免费视频| 亚洲一区二区三区精品动漫| 黄色综合网站| 午夜精品亚洲| 久久国产66| 在线一区欧美| 亚洲国产激情| 欧美日韩一区二区三区四区在线观看 | 亚洲九九精品| 午夜日韩在线| 久久九九免费| 亚洲永久免费| 99亚洲视频| 亚洲天堂久久| 欧美视频官网| 欧美欧美全黄| 欧美激情第六页| 欧美一区三区二区在线观看| 午夜在线视频观看日韩17c| 亚洲巨乳在线| 亚洲久久一区| 一区二区三区四区五区精品| 亚洲成人在线视频网站| 狠狠88综合久久久久综合网| 午夜精品国产| 欧美精品三区| 欧美全黄视频| 欧美三区视频| 激情久久综合| 在线播放不卡| 影音先锋亚洲精品| 亚洲激情女人| 99成人免费视频| 亚洲少妇一区| 亚洲一区免费| 老司机精品视频网站| 美女日韩在线中文字幕| 久久久久久色| 欧美a级片网站| 国产一区自拍视频| 亚洲欧洲久久| 亚洲一区精彩视频| 久久久精品国产一区二区三区| 久久久综合网| 国产伊人精品| 亚洲视频播放| 老牛嫩草一区二区三区日本| 欧美日韩国产精品一区二区亚洲| 国内精品视频在线播放| 99精品热视频只有精品10| 亚洲欧美日韩在线观看a三区| 久久综合九色综合久99| 欧美视频福利| 99亚洲一区二区| 性xx色xx综合久久久xx| 欧美日韩在线高清| 中文日韩欧美| 欧美日韩 国产精品| 夜夜精品视频| 欧美 亚欧 日韩视频在线| 一区在线播放| 久久精品卡一| 亚洲另类黄色| 欧美日韩一区综合| 午夜在线一区| 亚洲黄色一区| 女同性一区二区三区人了人一 | 午夜精品剧场| 亚洲美女少妇无套啪啪呻吟| 久久蜜桃精品| 国产日本精品| 国产精品国码视频| 免费在线播放第一区高清av| 在线成人www免费观看视频| 噜噜爱69成人精品| 制服诱惑一区二区| 亚洲一级特黄| 欧美另类高清视频在线| 欧美亚洲免费在线| 亚洲久久一区二区| 韩日成人在线| 欧美精品亚洲精品| 久久成人精品| 亚洲一区二区三区高清| 亚洲激情另类| 激情综合亚洲| 国产主播一区二区三区四区| 久久一区免费| 久久精品午夜| 麻豆精品网站| 亚洲中午字幕| 一区二区精品| 亚洲伦伦在线| 亚洲精品视频一区二区三区| 亚洲天堂男人| 影院欧美亚洲| 亚洲国产精品一区二区第一页| 国产精品国产一区二区| 欧美特黄视频| 欧美日本三区| 欧美人成在线| 国内精品久久久久久久97牛牛| 欧美精品综合| 国产一区二区三区四区老人| 欧美日韩在线精品| 国产精品观看| 影音先锋亚洲一区| 亚洲欧洲在线一区| 在线视频免费在线观看一区二区| 亚洲日本精品国产第一区| 亚洲人体大胆视频| 国产精品亚洲产品| 久久久久久精| 午夜精品久久99蜜桃的功能介绍| 欧美日韩亚洲一区二区三区在线| 国产精品xxx在线观看www| 激情视频一区二区三区| 一道本一区二区| 六月丁香综合| 国内久久精品| 国产日韩精品久久| 久久国产直播| 国产综合欧美| 亚洲视频www| 噜噜噜噜噜久久久久久91| 欧美日韩国产高清| 一本一本久久a久久精品综合妖精| 国产婷婷精品| 久久综合久久久| 好看的日韩av电影| 国产精品一卡| 欧美激情一区| 亚洲精品偷拍| 久久午夜精品一区二区| 亚洲国产精品久久久久婷婷老年| 亚洲一区二区三区精品视频| 欧美精品啪啪| 亚洲美女色禁图| 午夜久久福利| 亚洲欧美成人| 在线看片一区| 久久精品女人| 一本久久知道综合久久| 欧美日韩午夜| 性欧美暴力猛交另类hd| 亚洲国产一区二区三区a毛片| 麻豆久久久9性大片| 亚洲国产精品一区二区第一页| 美日韩精品视频| 亚洲伦理精品| 亚洲一级高清| 欧美精品啪啪| 久久精品五月婷婷| 国产精品夜夜夜一区二区三区尤| 亚洲先锋成人| 欧美日韩一区二区高清| 久久精品电影| 亚洲一区三区电影在线观看| 亚洲日本无吗高清不卡| 欧美日韩岛国| 女女同性女同一区二区三区91| 国产视频一区在线观看一区免费| 国语对白精品一区二区| 欧美国产专区| 久久午夜精品一区二区| 亚洲在线国产日韩欧美| 国产一区成人| 国产精品久久九九| 国产视频精品网| 国产亚洲毛片在线| 中文在线不卡| 99精品国产99久久久久久福利| 亚洲图片欧洲图片日韩av| 国产精品成人观看视频免费| 欧美日韩在线一二三| 欧美伊人影院| 欧美日本精品| 欧美日韩福利| 国产精品红桃| 国内精品国语自产拍在线观看| 欧美精品一卡| 狠狠色丁香久久综合频道| 国产精品theporn| 狠狠色狠狠色综合日日tαg| 亚洲一级高清| 亚洲精品一区二区三区樱花| 夜夜嗨网站十八久久| 一区二区精品在线观看| 亚洲欧美日韩国产| 欧美aⅴ99久久黑人专区| 欧美日韩高清在线一区| 极品中文字幕一区| 在线亚洲国产精品网站| 亚洲综合精品四区| 欧美aa国产视频| 欧美日韩理论| 亚洲精品色图| 性8sex亚洲区入口| 欧美日韩在线一二三| 亚洲国产综合在线看不卡| 国产精品美女黄网| 欧美国产高清| 在线日韩av| 亚洲一区二区三区欧美| 午夜国产一区| 亚洲精品社区| 久久久精品五月天| 亚洲小说区图片区| 亚洲一区激情| 黑人一区二区三区四区五区| 国产亚洲欧美一区二区| 久久亚洲午夜电影| 在线成人亚洲| 久久精品亚洲| 亚洲日本久久| 欧美aa国产视频| 亚洲美女毛片| 欧美日韩在线不卡一区| 国产亚洲一区在线播放| 欧美 日韩 国产在线| 一本一道久久综合狠狠老精东影业| 久久久国产精品一区二区中文| 国产精品magnet| 亚洲欧美日韩另类精品一区二区三区| 欧美黄色aaaa| 先锋影音久久久| 亚洲黄色毛片| 牛夜精品久久久久久久99黑人| 99精品视频免费观看| 欧美日韩另类丝袜其他| 噜噜噜在线观看免费视频日韩| 亚洲欧洲日韩综合二区| 欧美日韩免费观看一区| 香蕉久久夜色精品国产| 亚洲精品在线视频观看| 欧美福利网址| 亚洲一区二区三区高清| 亚洲激情成人| 精品动漫3d一区二区三区免费版| 欧美一区激情| 午夜亚洲福利在线老司机| 欧美日韩免费观看一区=区三区| 欧美日韩综合网| 国产一区二区三区的电影| 欧美日韩在线精品一区二区三区| 欧美一级久久| 国产精品日韩欧美一区二区| 亚洲久久在线| 亚洲国产精品一区二区第一页| 欧美日韩在线精品一区二区三区| 亚洲一区精彩视频| 一区二区三区四区五区精品视频| 国产自产在线视频一区| 午夜精品一区二区在线观看| 久久国产精品一区二区三区四区| 一区二区福利| 在线亚洲免费| 国产精品一区在线观看| 国产精品久久久免费| 亚洲精品在线观看免费| 亚洲第一伊人| 亚洲精品专区| 99亚洲视频| 国产精品乱看| 欧美中文字幕| 久久久久久久久久码影片| 久久精品成人| 欧美日本亚洲| 狠狠爱www人成狠狠爱综合网| 狠狠88综合久久久久综合网| 尹人成人综合网| 亚洲经典自拍| 亚洲专区一区| 欧美一区网站| 在线观看欧美亚洲| 国产精品免费在线| 久久亚洲色图| 尤物在线精品| 国产九九精品| 欧美日韩午夜| 一本色道久久99精品综合| 午夜在线视频观看日韩17c| 女女同性女同一区二区三区91| 好看的亚洲午夜视频在线| 国产欧美另类| 欧美激情日韩| 国产情侣一区| 欧美精品首页| 亚洲经典在线| 久久综合一区| 日韩午夜在线| 欧美日韩亚洲一区二区三区在线| 日韩视频一区| 亚洲欧美综合| 国产精品美女xx| 欧美私人啪啪vps| 亚洲欧美不卡| 亚洲国产电影| 欧美精品在线一区| 午夜在线精品| 日韩一区二区久久| 国产精品激情电影| 久久久福利视频|