亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        雙重降維HOG結(jié)合SVM的快速手指靜脈識別

        2022-04-08 08:43:26褚洪佳陳光化汪凱旋
        紅外技術(shù) 2022年3期
        關(guān)鍵詞:降維手指區(qū)間

        褚洪佳,陳光化,2,汪凱旋

        雙重降維HOG結(jié)合SVM的快速手指靜脈識別

        褚洪佳1,陳光化1,2,汪凱旋1

        (1. 上海大學 微電子研究與開發(fā)中心,上海 200444;2. 上海大學機電工程與自動化學院,上海 200444)

        為減少手指靜脈識別時間,提出一種雙重降維方向梯度直方圖特征(Histogram of Oriented Gradient,HOG)結(jié)合支持向量機(Support Vector Machine,SVM)分類的手指靜脈識別方法。針對傳統(tǒng)HOG算法特征維數(shù)高的問題,首先通過Fisher準則衡量梯度方向區(qū)間HOG特征的分類能力,然后使用序列前向選擇法挑選出分類能力較優(yōu)異的梯度方向區(qū)間構(gòu)建部分方向區(qū)間HOG特征,最后使用主成分分析(Principal Component Analysis,PCA)降維。在公開的手指靜脈數(shù)據(jù)庫FV-USM和THU-FV上使用SVM多分類器進行分類識別,實驗結(jié)果表明:雙重降維HOG方法相較于HOG+PCA方法提取的特征維數(shù)降低了40%,識別時間減少了29.85%,識別準確率分別為99.17%和100%,等誤率分別為1.07%和0.01%。

        手指靜脈識別;方向梯度直方圖;特征選擇;主成分分析;支持向量機

        0 引言

        手指靜脈識別技術(shù)使用手指內(nèi)部的靜脈作為特征信息,使用近紅外光作為光源進行圖像采集,相比于傳統(tǒng)的人臉識別和指紋識別方式,具有活體性、高穩(wěn)定性和高安全性,在銀行、監(jiān)獄以及老年醫(yī)療等場合有著極大的優(yōu)勢。近年來手指靜脈識別技術(shù)成為生物特征識別領(lǐng)域的研究熱點。

        現(xiàn)有對手指靜脈識別的研究方法中,常采用基于紋路拓撲結(jié)構(gòu)的模板匹配方法[1-2],首先從手指靜脈圖像中分割出靜脈拓撲結(jié)構(gòu)作為特征信息,然后計算模板間的Hausdorff距離進行匹配識別。但是此方法魯棒性較差,識別準確率容易受光照和圖像背景噪聲等因素的影響。除了基于紋路拓撲結(jié)構(gòu)的模板匹配方法外,還有神經(jīng)網(wǎng)絡(luò)方法[3-4]和局部特征描述方法[5-6]?;谏窠?jīng)網(wǎng)絡(luò)的方法,無需對手指靜脈進行特征提取[7],可以根據(jù)已有的數(shù)據(jù)集學習到有價值的信息特征,然后對待識別手指靜脈圖像進行識別。但是神經(jīng)網(wǎng)絡(luò)方法對數(shù)據(jù)集的數(shù)據(jù)量要求比較高,當訓練數(shù)據(jù)集的數(shù)據(jù)量較少時,所訓練的網(wǎng)絡(luò)識別效果欠佳。局部特征描述方法主要包括局部二值模式[8](Local Binary Pattern,LBP)和方向梯度直方圖[9](histogram of oriented gradient,HOG),此類方法使用局部特征描述算子對手指靜脈特征進行編碼獲得圖像特征,然后使用支持向量機[10](support vector machines,SVM)、線性判別分析[11](linear discriminant analysis,LDA)、adaboost[12]等方法進行分類識別。文獻[8]使用改進LBP描述符提取手指靜脈特征,時間遠高于HOG方法。雖然HOG特征提取時間開銷少,但是采用部分區(qū)域重疊的方法逐區(qū)域進行特征提取[13]使得特征的維度非常高,導致后續(xù)特征識別過程的計算復雜度高,增加了識別過程的時間開銷。文獻[14]引入Harr積分圖簡化HOG計算過程,減少了特征提取時間,但沒有解決HOG特征維數(shù)高的問題。文獻[9]使用主成分分析(Principal Component Analysis,PCA)算法對手指靜脈HOG特征降維后進行分類,特征識別時間有明顯改善。但是PCA針對全局HOG特征進行降維,降維準則是保留對方差貢獻度大的特征,未考慮局部的HOG特性,存在降維不充分的問題。

        鑒于上述分析,本文提出一種雙重降維HOG特征結(jié)合SVM分類的手指靜脈識別模型,模型的總體結(jié)構(gòu)如圖1所示。首先,使用Sobel邊緣檢測算法提取手指靜脈圖像的感興趣區(qū)(region of interest,ROI)。其次,使用序列前向選擇法挑選出識別能力較優(yōu)異的梯度方向區(qū)間構(gòu)建部分方向區(qū)間HOG特征,并使用PCA降維輸出PCA子空間模型。然后,使用SVM對降維后的特征訓練輸出SVM分類模型。最后,使用訓練過程輸出的PCA子空間模型和SVM分類模型識別待識別圖像。

        1 雙重降維HOG手指靜脈特征提取

        1.1 ROI提取

        靜脈中的血紅蛋白對近紅外光的吸收能力較強,通常使用近紅外光作為光源來采集手指靜脈圖像。采集到的手指靜脈圖像進行預處理,定位手指靜脈圖像感興趣區(qū),去除背景區(qū)域,有利于后續(xù)進行特征提取。本文使用Sobel邊緣檢測算法定位手指邊緣,利用手指邊緣信息和手指指節(jié)信息提取感興趣區(qū)域,最后將感興趣區(qū)尺寸標準化,過程如圖2所示。

        圖1 手指靜脈識別總體結(jié)構(gòu)圖

        圖2 手指靜脈圖像預處理過程

        Sobel邊緣檢測算法使用區(qū)域模板與圖像進行卷積得到圖像在兩個方向上的梯度GG,然后對兩個方向上的梯度求模得到梯度幅值,通過梯度幅值與預設(shè)的閾值進行比較來確定圖像的邊緣。使用圖3中模板計算的過程如式(1)、(2)、(3)。

        G=(13+223+33)-(11+221+31) (1)

        G=(11+212+13)-(31+232+33) (2)

        1.2 部分方向區(qū)間HOG特征

        HOG特征是一種基于圖像局部梯度信息統(tǒng)計的特征描述符,其理論依據(jù)是一幅圖像中的局部形狀和邊緣特征能夠被梯度方向和梯度幅值所描述。

        使用圖4所示Sobel模板對圖像做卷積運算,計算像素點水平方向和垂直方向的梯度幅值G、G,根據(jù)兩個方向上的梯度計算像素點的梯度幅值和梯度方向

        將梯度方向范圍均勻分成9個方向區(qū)間,每個方向區(qū)間對應(yīng)cell直方圖的一個特征維度,根據(jù)像素點的梯度信息將其投影至直方圖的對應(yīng)特征維度中,投影權(quán)重為像素點的梯度幅值。按照cell、block和圖像三者間的構(gòu)成關(guān)系,分別將所有cell中相同的方向區(qū)間特征進行串聯(lián),得到9個單方向區(qū)間HOG特征向量。

        以Fisher準則作為衡量標準,將9個單方向區(qū)間HOG特征向量按照Fisher值從大到小排列。Fisher準則尋求一個最優(yōu)子空間使樣本類間散度盡量大,類內(nèi)散度盡量小,二者比值越大說明分類能力越強。Fisher計算公式如式(5)所示:

        式中:tr(·)是矩陣跡運算;是樣本進行線性變換的子空間;b是樣本的總類間散度;w是樣本的總類內(nèi)散度??傤愰g散度和總類內(nèi)散度的計算方法如式(6)和式(7)所示:

        使用序列前向選擇法選擇單方向區(qū)間HOG特征向量中的個構(gòu)建部分方向區(qū)間HOG特征向量。序列前向選擇法流程如圖5所示,依次增加排序后的單方向區(qū)間HOG特征的數(shù)目,直到分類準確率滿足要求。

        1.3 PCA降維

        PCA把原始的維特征映射到維空間上,映射后的維是在原有維特征的基礎(chǔ)上構(gòu)建出的全新的正交特征,在降低數(shù)據(jù)集維數(shù)的同時保留了對數(shù)據(jù)集的方差貢獻度最大的特征。

        根據(jù)公式(8)計算手指靜脈樣本庫的部分方向HOG特征向量的協(xié)方差矩陣。求取協(xié)方差矩陣的特征值=(1,2, …,)及其對應(yīng)的特征向量=(1,2, …,v),并根據(jù)特征值降序排列特征向量。

        式中:m為樣本數(shù)量;X為樣本特征向量;為樣本均值。

        根據(jù)公式(9)計算協(xié)方差矩陣特征值的累計貢獻率,貢獻率取值通常取85%~99%之間,根據(jù)貢獻率取值要求選擇前個較大的特征值所對應(yīng)的特征向量V=(1,2, …,v)作為主成分分析的特征子空間的投影矩陣pca=(1,2, …,v)。最后將輸入數(shù)據(jù)投影至子空間。

        2 基于SVM的手指靜脈識別

        支持向量機通過尋找一個超平面把兩類樣本完全分開并且兩類間的幾何間隔最大。設(shè)超平面的方程式為:·+=0,則可以把·+>0的數(shù)據(jù)歸為1類,把·+<0的數(shù)據(jù)歸為-1類。支持向量機求解超平面問題可以表示為式(10)所示的約束最優(yōu)化問題:

        式中:懲罰因子和松弛變量用于對目標函數(shù)和約束條件進行修正。引入拉格朗日乘子,超平面參數(shù)和可以使用對偶問題的解來確定,由此可以得到支持向量機的最優(yōu)分類函數(shù):

        式中:為核函數(shù),用于解決非線性分類問題,常用的核函數(shù)有線性核函數(shù)(Linear Kernel)、多項式核函數(shù)(polynomial kernel)、徑向基核函數(shù)(RBF kernel)和sigmoid核函數(shù)(sigmoid kernel)。在訓練集對4種核函數(shù)進行測試,相關(guān)指標如表1所示,線性核函數(shù)訓練的支持向量數(shù)量少于其他3種,識別準確率最高。此外線性核函數(shù)復雜度較低,故本文設(shè)計分類器使用線性核函數(shù)。

        表1 核函數(shù)分類訓練數(shù)據(jù)

        基于線性核函數(shù)的支持向量機還需要確定懲罰因子來獲取更好的識別性能。使用倍增的懲罰因子2-4, 2-3, 2-2, 2-1, 20, 21, 22, 23, 24進行測試,分類器識別準確率如圖6所示,在懲罰因子為4時識別準確率達到最高,因此選用=4作為支持向量機的懲罰因子。

        圖6 懲罰因子與準確率關(guān)系圖

        實際應(yīng)用中,手指靜脈識別過程為多分類問題,基本的二分類SVM模型無法滿足實際需求。本文對任意兩類手指靜脈都訓練一個二分類SVM模型,并根據(jù)決策函數(shù)確定二分類SVM的分類結(jié)果。最后,對所有的二分類SVM模型的分類結(jié)果進行投票表決,得票數(shù)最多的分類結(jié)果即為多分類SVM的分類結(jié)果。圖7所示為基于二分類SVM的四分類器,其中SVM[,]為使用第類和第類樣本訓練的二分類器。

        3 實驗與分析

        本文使用的計算機配置為8.0G內(nèi)存和1.6GHz主頻CPU,實驗平臺為Matlab2018a,使用的數(shù)據(jù)集為馬來西亞理工大學的FV-USM手指靜脈數(shù)據(jù)庫和清華大學的THU-FV手指靜脈數(shù)據(jù)庫。FV-USM數(shù)據(jù)庫中共有492個分類,每個分類在兩個階段共采集12張圖像。THU-FV數(shù)據(jù)庫中共有610個分類,每個分類采集8張圖像。本文實驗將每個分類中的前4張圖像作為訓練樣本,其余作為測試樣本,所有圖像尺寸標準化為128×64像素。

        圖7 基于SVM的四分類器

        3.1 方向區(qū)間數(shù)目的確定

        為確定部分方向區(qū)間HOG特征中包含的方向區(qū)間數(shù)目,本文使用FV-USM數(shù)據(jù)集中的手指靜脈圖像進行實驗,得到方向區(qū)間數(shù)目與手指靜脈識別準確率的關(guān)系如圖8所示。從圖中可知,選擇前3個方向區(qū)間時的識別準確率為98.83%,十分接近包含9個方向區(qū)間的傳統(tǒng)HOG的識別準確率99.02%,選擇前4個方向區(qū)間的識別準確率首次取得最大值99.17%。繼續(xù)增加方向區(qū)間數(shù)目,識別準確率非但沒有增加,反而在某些組合情況下出現(xiàn)輕微下降。原因在于不同的方向區(qū)間HOG特征向量之間存在相關(guān)性,后續(xù)增加的特征沒有提供新的更具區(qū)分能力的信息,反而可能引入噪聲信息,導致識別準確率降低。綜上所述,后續(xù)實驗中部分方向區(qū)間HOG特征包含的方向區(qū)間數(shù)目確定為4。

        圖8 方向區(qū)間數(shù)目與識別準確率關(guān)系圖

        3.2 時間開銷分析

        為評估本文提出的雙重降維HOG特征提取方法,將其與經(jīng)典的局部特征算法LBP和HOG以及HOG+PCA方法在FV-USM數(shù)據(jù)庫上進行實驗,結(jié)果如表2所示。LBP算法的特征提取時間遠高于HOG算法的提取時間,原因在于LBP算法對于每個像素點要與周圍的8個像素點進行比較后進行編碼和直方圖統(tǒng)計,計算量大。本文提出的雙重降維HOG方法的特征提取時間開銷略高于傳統(tǒng)HOG算法,與HOG+PCA方法基本持平,原因在于本文提出的雙重降維HOG特征沒有改變HOG算法的編碼原理,增加的僅是特征降維的過程。在識別時間方面,可以發(fā)現(xiàn)支持向量機識別時間與特征維數(shù)呈正相關(guān)關(guān)系,本文提出的雙重降維HOG方法提取的特征維數(shù)最少,相比于傳統(tǒng)HOG+PCA算法特征維數(shù)降低了40%,特征識別時間減少了29.85%。

        表2 時間開銷比較

        3.3 識別性能分析

        準確率(Accuracy)和等誤率(equal error rate,EER)是手指靜脈識別研究中常用的性能指標。準確率是測試時正確識別次數(shù)所占的比率,表示模型對當前數(shù)據(jù)集的適用程度。等誤率是拒識率(false rejection rate,F(xiàn)RR)和誤識率(false accept rate,F(xiàn)AR)近似相等時的拒識率或誤識率,等誤率的數(shù)值越小,說明識別效果越好。

        為評估本文提出的手指靜脈識別模型,將提出的算法模型與其他方法在FV-USM和THU-FV數(shù)據(jù)庫上進行測試,數(shù)據(jù)如表3所示,圖9為根據(jù)FV-USM數(shù)據(jù)庫的測試結(jié)果繪制的使用接受者操作特性(receive operating characteristic,ROC)曲線圖。從實驗數(shù)據(jù)中可以看出,本文提出的方法在兩個數(shù)據(jù)庫上都取得了最優(yōu)的準確率和等誤率,原因在于:在特征構(gòu)建過程中,依據(jù)Fisher準則選擇特征,使樣本類內(nèi)散度盡量小,類間散度盡量大;使用PCA保留對樣本方差貢獻度大的特征。表3中最后4組數(shù)據(jù)將歐式距離計算相似度的識別方法和SVM分類的識別方法進行對比,結(jié)果顯示使用SVM分類可以獲取更好的識別性能,原因在于SVM通過核函數(shù)將線性不可分的數(shù)據(jù)映射到高維空間進行求解,并通過改變懲罰因子的值來尋求更優(yōu)的分類能力。

        表3 不同方法實驗數(shù)據(jù)比較

        圖9 ROC曲線圖(FV-USM數(shù)據(jù)庫)

        實驗中各種方法在THU-FV數(shù)據(jù)庫的實驗數(shù)據(jù)優(yōu)于FV-USM數(shù)據(jù)庫的數(shù)據(jù),是因為THU-FV數(shù)據(jù)庫中的手指靜脈圖像質(zhì)量高并且來自同一手指的不同圖像間的差異極小,而FV-USM數(shù)據(jù)庫的圖像較為模糊,靜脈紋路不明顯,具體圖像如圖10所示。此外最大曲率+MHD方法和寬線檢測+LTS-HD方法都是基于靜脈紋路分割的手指靜脈識別方法,由于FV-USM數(shù)據(jù)庫中靜脈紋路比較模糊不易分割出完整的靜脈紋路拓撲,導致算法的識別準確率和等誤率數(shù)據(jù)較差。

        圖10 不同數(shù)據(jù)庫手指靜脈圖像

        4 結(jié)語

        本文提出一種雙重降維HOG特征結(jié)合SVM分類的手指靜脈識別模型。依據(jù)Fisher準則衡量各梯度方向區(qū)間的分類能力,并使用序列前向選擇法挑選出分類能力強的梯度方向區(qū)間的組合構(gòu)建部分方向區(qū)間HOG特征,實現(xiàn)HOG特征的第一次降維。對部分方向HOG特征使用PCA降維,實現(xiàn)HOG特征的第二次降維?;诙诸怱VM構(gòu)建多分類器用于所提取特征的分類識別。在公開的手指靜脈數(shù)據(jù)庫FV-USM和THU-FV上進行實驗的結(jié)果表明:本文提出的雙重降維HOG特征方法能夠有效降低特征維數(shù),使用SVM分類識別的時間有明顯改善,并且獲得了較高的識別準確率和較低的等誤率,能夠?qū)崿F(xiàn)手指靜脈的精準識別。

        [1] YANG L, YANG G P, YIN Y L, et al. Finger vein recognition with anatomy structure analysis[J]., 2018, 28(8): 1892-1905.

        [2] CHEN G H, DAI QH, TANG X, et al. An improved least trimmed square hausdorff distance finger vein recognition[C]//, 2018: 939-943.

        [3] LI S Y, ZHANG H G, YANG JF. Finger vein recognition based on local graph structural coding and CNN[C]//, 2019, 11069: 110693I-110693I-8.

        [4] ZHANG Y K, LI W J, ZHANG L P, et al. Adaptive Gabor convolutional neural networks for finger-vein recognition[C]//(Hpbd&Is), 2019: 219-222.

        [5] LIU H Y, YANG L, YANG G P, et al. Discriminative binary descriptor for finger vein recognition[J]., 2018, 6: 5795-5804.

        [6] WANG X, WANG H B, HE Y, et al. Novel Algorithm for finger vein recognition based on inception-Resnet module[J]., 2019, 11179: 111791D-111791D-9.

        [7] 陶志勇, 胡亞磊, 林森. 基于改進AlexNet的手指靜脈識別[J]. 激光與光電子學進展, 2020, 57(8): 58-66.

        TAO Zhiyong, HU Yalei, LIN Sen. Finger vein recognition based on improved AlexNet[J]., 2020, 57(8): 58-66.

        [8] 劉超, 王容川, 許曉偉, 等. 基于改進LBP的手指靜脈識別算法[J]. 計算機仿真, 2019, 36(1): 381-386.

        LIU Chao, WANG Rongchuan, XU Xiaowei, et al. Finger vein recognition algorithm based on improved LBP[J]., 2019, 36(1): 381-386.

        [9] 李菲, 李小霞, 周穎玥. 基于改進HOG特征和稀疏表示的手指靜脈識別[J]. 傳感器與微系統(tǒng), 2018, 37(11): 38-41.

        LI Fei, LI Xiaoxia, ZHOU Yingyue. Finger vein recognition based on improved HOG features and sparse representation[J]., 2018, 37(11): 38-41, 44.

        [10] Veluchamy S, Karlmarx L R. System for multimodal biometric recognition based on finger knuckle and finger vein using feature-level fusion and k-support vector machine classifier[J]., 2017, 6(3): 232-242.

        [11] 徐子豪, 陳光化, 傅志威. 改進型LDA結(jié)合LBP的手指靜脈識別[J]. 現(xiàn)代電子技術(shù), 2020, 43(12): 1-4.

        XU Zihao, CHEN Guanghua, FU Zhiwei. Finger vein recognition of improved LDA combined with LBP[J]., 2020, 43(12): 1-4.

        [12] 徐鑄業(yè), 趙小強. 基于Agast-Adaboost的圖像匹配算法[J]. 蘭州理工大學學報, 2020, 46(4): 110-115.

        XU Zhuye, ZHAO Xiaoqiang. Image matching algorithm based on Agast-Adaboost[J]., 2020, 46(6): 1-4.

        [13] 賈楚. 基于改進HOG特征的行人檢測算法研究[D]. 秦皇島: 燕山大學, 2016.

        JIA Chu. Research of Pedestrian Detection Based on Improved HOG Features[D]. Qinhuangdao: Yanshan University, 2016.

        [14] 蔣政, 程春玲. 基于Haar特性的改進HOG的人臉特征提取算法[J]. 計算機科學, 2017, 44(1): 303-307.

        JIANG Zheng, CHENG Chunling. Improved HOG face feature extraction algorithm based on Haar characteristics[J]., 2017, 44(1): 303-307.

        Fast Finger Vein Recognition Based on a Dual Dimension Reduction Histogram of Oriented Gradient and Support Vector Machine

        CHU Hongjia1,CHEN Guanghua1,2,WANG Kaixuan1

        (1.,,200444,;2.,,200444,)

        An identification model using a dual-dimension reduction histogram of oriented gradients (HOG) combined with a support vector machine (SVM) is proposed to reduce the time required for finger vein recognition. To solve the problem of high feature dimensionality in the traditional HOG algorithm, the classification ability of the gradient direction interval is first measured using the Fisher criterion. Next, the sequence forward selection method is used to select the gradient direction interval with optimal classification ability to construct a partial direction interval HOG feature. Finally, principal component analysis (PCA) is used to reduce the number of dimensions. An SVM multi-classifier was used for the classification of the FV-USM and THU-FV datasets. The experimental results demonstrate that compared to the HOG+PCA method, the feature dimensions extracted by the dual-dimensional reduction HOG method are reduced by 40%, the recognition time is reduced by 29.85%, the recognition accuracy is 99.17% and 100%, respectively, and the equal error rate is 1.07% and 0.01%, respectively.

        finger vein recognition, HOG, feature selection, PCA, SVM

        國家自然科學基金項目(61671285)。

        TP391.4

        A

        1001-8891(2022)03-0262-06

        2021-01-23;

        2021-04-08.

        褚洪佳(1994-),男,山東省棗莊市人,碩士研究生,研究方向為圖像處理、模式識別。E-mail: chu_hongjia@163.com

        猜你喜歡
        降維手指區(qū)間
        Three-Body’s epic scale and fiercely guarded fanbase present challenges to adaptations
        解兩類含參數(shù)的復合不等式有解與恒成立問題
        你學會“區(qū)間測速”了嗎
        手指背側(cè)復合組織塊離斷再植
        手指操
        降維打擊
        海峽姐妹(2019年12期)2020-01-14 03:24:40
        區(qū)間對象族的可鎮(zhèn)定性分析
        筆穿手指
        少年科學(2015年8期)2015-08-13 22:13:01
        拋物化Navier-Stokes方程的降維仿真模型
        計算物理(2014年1期)2014-03-11 17:00:18
        基于特征聯(lián)合和偏最小二乘降維的手勢識別
        无码精品一区二区三区免费16 | 污污污污污污污网站污| 一本久久综合亚洲鲁鲁五月天| 免费在线不卡黄色大片| 中文乱码字幕高清在线观看| 亚洲乱码av中文一区二区| 狠狠色丁香婷婷久久综合| 人妻少妇偷人精品视频| 精品国免费一区二区三区| 国产女合集小岁9三部| 国产md视频一区二区三区| 午夜精品久久99蜜桃| 国产一区二区三区韩国| 国产午夜福利小视频合集| 无码人妻精品一区二区三区9厂| 久久精品亚洲94久久精品| 在线视频一区二区在线观看| 亚洲第一看片| 亚洲欧美日韩中文在线制服| 国产成人精品日本亚洲专区61| 人妻经典中文字幕av| 伊人影院在线观看不卡| 国产乱人伦AV在线麻豆A| 亚洲av无码精品色午夜蛋壳| 精品国产av一区二区三区| 日韩av在线播放人妻| 国产精品美女一区二区av| 国产精品亚洲婷婷99久久精品 | 中文字幕一区二区三区综合网| 人妻无码∧V一区二区| 日韩区在线| 亲子乱aⅴ一区二区三区下载| 免费视频爱爱太爽了| 国产成人精品优优av| 久久精品国产亚洲av成人文字| 亚洲另类国产精品中文字幕| 欧美h久免费女| 在线精品日韩一区二区三区| 国产剧情国产精品一区| 久久精品国产www456c0m| 久久综合亚洲色hezyo国产|