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

        ?

        基于小波變換與脈沖耦合神經(jīng)網(wǎng)絡(luò)的人臉識別①

        2017-10-20 03:09:04何世強(qiáng)劉金清蔡淑寬陳存弟周曉童鄧淑敏吳慶祥
        關(guān)鍵詞:小波人臉識別人臉

        何世強(qiáng),劉金清,劉 引,蔡淑寬,陳存弟,周曉童,鄧淑敏,吳慶祥

        (福建師范大學(xué) 光電與信息工程學(xué)院 醫(yī)學(xué)光電科學(xué)與技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,福州 350007)

        基于小波變換與脈沖耦合神經(jīng)網(wǎng)絡(luò)的人臉識別①

        何世強(qiáng),劉金清,劉 引,蔡淑寬,陳存弟,周曉童,鄧淑敏,吳慶祥

        (福建師范大學(xué) 光電與信息工程學(xué)院 醫(yī)學(xué)光電科學(xué)與技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,福州 350007)

        脈沖耦合神經(jīng)網(wǎng)絡(luò)(Pulse Coupled Neural Network,PCNN)是基于生物視覺特性而提出的新一代人工神經(jīng)網(wǎng)絡(luò),它在數(shù)字圖像處理及人工智能等領(lǐng)域具有廣泛應(yīng)用前景.本文通過研究PCNN理論模型及其工作特性的基礎(chǔ)上提出了一種提取人臉特征的方法.首先利用小波變換提取人臉圖像低頻特征,降低人臉圖像的維度,然后利用簡化的PCNN提取小波低頻系數(shù)重構(gòu)后的人臉圖像的相應(yīng)時間序列,并以此作為人臉識別的特征序列.最后利用時間序列和歐式距離完成人臉的識別過程.本文通過ORL人臉庫進(jìn)行實(shí)驗(yàn)證明了該方法的有效性.

        脈沖神經(jīng)網(wǎng)絡(luò); 人臉識別; 小波變換; 時間序列; 歐氏距離

        人臉識別,是一種基于人臉部信息進(jìn)行識別的生物特征識別技術(shù).人們在鑒定一個人的身份時,一般都是通過人臉來判斷的,人們從人臉獲得的判別信息達(dá)到90%以上.這給人臉識別技術(shù)的發(fā)展提供了理論依據(jù).人臉識別相對于其它生物特征識別(掌紋、指紋、虹膜等)具有方便、友好、不用接觸就能獲取信息等優(yōu)點(diǎn),人臉識別廣泛應(yīng)用于刑偵破案、門禁系統(tǒng)、網(wǎng)絡(luò)應(yīng)用、支付系統(tǒng)等等,是數(shù)字圖像處理和機(jī)器視覺領(lǐng)域的研究熱點(diǎn)[1].在國內(nèi)外有許多的研究人員在人臉識別上做了大量的研究并提出了許多有意義的算法.雖然人臉識別技術(shù)已經(jīng)取得可喜的成果,但是由于人臉的復(fù)雜性,在實(shí)際應(yīng)用中仍然面臨著許多嚴(yán)峻的問題,因此對人臉識別進(jìn)行研究具有重要意義.目前主要的人臉識別方法有:基于幾何特征的方法[2]、基于子空間分析的方法[3]和基于神經(jīng)網(wǎng)絡(luò)的方法[4]等.

        人臉主要由眼睛、鼻子、嘴巴、下巴等器官構(gòu)成.正是由于這些器官的形狀、大小以及結(jié)構(gòu)上的差異使得人臉千差萬別,也正是這些差異為人臉識別幾何特征提取奠定了理論基礎(chǔ),它主要是根據(jù)人臉的輪廓及曲線來確定若干顯著點(diǎn),并利用這些特征點(diǎn)作為識別特征,該方法在表情變化、姿態(tài)變化和光照變化條件下很難提取穩(wěn)定的特征點(diǎn).子空間方法作為一種流行的人臉識別方法,它利用線性或非線性的空間變換,使數(shù)據(jù)在映射在相應(yīng)的子空間內(nèi)降低數(shù)據(jù)的維度以便實(shí)現(xiàn)更好的模式分類,如基于線性變換的PCA[5]、ICA[6]和LDA[7]等.但這些方法對圖像旋轉(zhuǎn)和畸變的適應(yīng)性較差.神經(jīng)網(wǎng)絡(luò)的方法利用神經(jīng)網(wǎng)絡(luò)具有學(xué)習(xí)和非線性分類的能力來實(shí)現(xiàn)人臉識別,它相對于其它識別方法具有獨(dú)到的優(yōu)點(diǎn),它避免了復(fù)雜的特征提取工作,可以用神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過程來獲得人臉的規(guī)律特征,常用神經(jīng)網(wǎng)絡(luò)有BP神經(jīng)網(wǎng)絡(luò)和SOM神經(jīng)網(wǎng)絡(luò)等.但這些傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)模型在樣本學(xué)習(xí)過程中,容易出現(xiàn)過擬合的現(xiàn)象,而且需要較長的網(wǎng)絡(luò)訓(xùn)練時間.

        脈沖耦合神經(jīng)網(wǎng)絡(luò)[8]是基于觀測貓的大腦視覺皮層與視覺感知特征相關(guān)的神經(jīng)元同步行為而提出的一種新型人工神經(jīng)網(wǎng)絡(luò)模型,它是一種單層的,以迭代為主的自監(jiān)督自學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)模型.與傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)相比它不需要提前進(jìn)行訓(xùn)練并且具有很好的旋轉(zhuǎn)、平移、縮放不變特性[9].本文通過分析簡化的脈沖耦合神經(jīng)網(wǎng)絡(luò)模型提出了一種基于小波變化與脈沖耦合神經(jīng)網(wǎng)絡(luò)相結(jié)合的人臉特征提取方法.該方法利用脈沖耦合神經(jīng)網(wǎng)絡(luò)模擬生物的視覺感知過程,將人臉分解為若干個二值圖像,提取每幅圖像的時間序列[10]作為人臉識別的特征,利用歐式距離進(jìn)行分類識別,仿真實(shí)驗(yàn)證明了該方法的有效性.

        1 基于小波變換特征提取

        1.1 小波變換原理

        小波分析是數(shù)學(xué)領(lǐng)域的一個新的分支,它是泛函分析、Fourier分析、樣條分析、數(shù)值分析的完美結(jié)晶.小波變換的實(shí)質(zhì)是對信號用一組不同尺度的高通和低通濾波器組進(jìn)行濾波,將信號的高頻和低頻信息分解到不同的頻帶上進(jìn)行分析處理,然后重復(fù)以上過程,直到達(dá)到設(shè)定的閾值為止.與傅里葉變換相比,小波變換是時域和頻域的局部變換,因而能從信號中有效地提取信息,通過伸縮和平移等運(yùn)算對信號進(jìn)行多尺度分析,克服了傅里葉變換的窗口大小不能隨頻率變換的缺點(diǎn),被譽(yù)為“數(shù)學(xué)顯微鏡”.小波分析在信號處理、圖像分析、非線性科學(xué)等領(lǐng)域已有重大突破,成為了20世紀(jì)最輝煌的科學(xué)成就之一[11].

        小波變換是通過對基本小波進(jìn)行尺度伸縮和位移得到的.基本小波是一個具有特殊性質(zhì)的實(shí)值函數(shù),其振蕩快速衰減,且在數(shù)學(xué)上滿足積分為零的條件,即:

        其頻譜滿足條件:

        即基本小波在頻域也具有較好的衰減性質(zhì).

        小波基函數(shù)是通過尺度因子和位移因子由基本小波產(chǎn)生的,即:

        其中,a 為尺度因子,b 為平移量.由表達(dá)式可以得出,在連續(xù)變化的a,b值下具有很強(qiáng)的關(guān)聯(lián)性,因此,小波變換系數(shù)的信息量是冗余的,而冗余信息不利于對信號的分析和處理,所以對尺度和時移參數(shù)進(jìn)行離散化后,再進(jìn)行小波重構(gòu).

        1.2 小波的圖像壓縮技術(shù)

        所謂的圖像壓縮就是去掉冗余的信息而保留重要的信息.在給定的原始圖像中往往存在著許多冗余的信息(如空間冗余、時域冗余、頻域冗余、信息熵冗余等等).對人臉圖像進(jìn)行小波壓縮,不僅可以降低人臉圖像的數(shù)據(jù)維度,加快后續(xù)過程的處理速度,同時能夠提取圖像的主要特征.小波基的選擇、分解層數(shù)的確定和分解子圖的選取是小波變換應(yīng)用于人臉識別中的幾個關(guān)鍵問題,對識別效果有著非常重要的影響.針對以上幾個問題,當(dāng)前還沒有很好的理論指導(dǎo)和解決方案,本文在的實(shí)驗(yàn)的基礎(chǔ)上本文通過實(shí)驗(yàn)比較選取的小波濾波器組為db1小波.

        一幅二維圖像經(jīng)過一次小波變換之后,可以得到低頻區(qū)域LL,高頻區(qū)域HH、LH、HL四個區(qū)域,如圖1(a)所示,其中低頻區(qū)域表示近似分量,為原圖像的平滑圖像,高頻區(qū)域分別表示水平分量、垂直分量和對角分量.

        對于圖像信號,小波變換要經(jīng)過一次行變換和一次列變換來完成.對于圖像,小波系數(shù)的空間分布同原始圖像的空間分布的對應(yīng)關(guān)系為:LL頻帶是圖像內(nèi)容的縮略圖,它是圖像數(shù)據(jù)能量集中的頻帶.而LH、HL和HH頻帶存放的是圖像的細(xì)節(jié)信息.

        二維小波分解遞推公式如下所示:

        公式中 Cj,m,n是圖像的低頻分量,它包含了圖像的主要能量,表示的是圖像的主要信息.和是圖像的高頻分量,表示了圖像在水平方向、垂直方向和對角方向上的細(xì)節(jié)分量.人臉圖像經(jīng)過一次小波分解及重構(gòu)后的圖像如圖1(b-d)所示.

        圖1 小波分解及低頻重構(gòu)圖

        2 脈沖耦合神經(jīng)網(wǎng)絡(luò)

        1990年,Eckhorn在研究貓視覺皮層研究中發(fā)現(xiàn)了視覺皮層神經(jīng)元同步脈沖發(fā)放現(xiàn)象; 1993年Johnson在Eckhorn的基礎(chǔ)上提出了脈沖耦合神經(jīng)網(wǎng)絡(luò)模型,被廣泛應(yīng)用于圖像分割[12]、圖像融合[13]及目標(biāo)檢測[14]等領(lǐng)域.由于原始的脈沖耦合神經(jīng)網(wǎng)絡(luò)參數(shù)多,不易優(yōu)化,許多科學(xué)研究者在保留脈沖耦合神經(jīng)網(wǎng)絡(luò)的視覺信息處理的基礎(chǔ)上對脈沖耦合神經(jīng)網(wǎng)絡(luò)進(jìn)行了簡化,其中比較流行的一種模型如圖2所示.

        簡化的脈沖耦合神經(jīng)網(wǎng)絡(luò)模型的數(shù)學(xué)表達(dá)式如下所示:

        由圖2可知簡化脈沖神經(jīng)網(wǎng)絡(luò)主要由三部分組成:接收域、調(diào)制域和脈沖產(chǎn)生器三部分組成.相對于原始的脈沖神經(jīng)網(wǎng)絡(luò)模型,簡化的脈沖耦合神經(jīng)網(wǎng)絡(luò)模型的參數(shù)個數(shù)較少,保留了原始模型的雙通道調(diào)制、動態(tài)閾值和時空累加特性.可以更快更有效地對圖像進(jìn)行處理.在簡化的PCNN中,接受域把外部信號激勵I(lǐng)ij作為神經(jīng)元的反饋通道的輸入項,神經(jīng)元的脈沖發(fā)生器部分仍然由內(nèi)部活動項Uij和變化的閾值θij構(gòu)成,當(dāng)閾值大于神經(jīng)元內(nèi)部活動值時,脈沖發(fā)放器被關(guān)閉,閾值按照指數(shù)規(guī)律衰減.當(dāng)閾值衰減到小于內(nèi)部狀態(tài)值時,脈沖發(fā)生器就會被重新被激活,此時神經(jīng)元就會再次輸出脈沖.所以總的來說,簡化的脈沖耦合神經(jīng)網(wǎng)絡(luò)具有非線性的雙通道調(diào)制和變閾值特性,在這點(diǎn)上它和標(biāo)準(zhǔn)的脈沖耦合神經(jīng)網(wǎng)絡(luò)是一致的,也是和傳統(tǒng)神經(jīng)網(wǎng)絡(luò)的重要區(qū)別.

        圖2 簡化 PCNN 模型

        3 脈沖耦合神經(jīng)網(wǎng)絡(luò)人臉特征提取

        將人臉圖像的每個像素點(diǎn)都與一個神經(jīng)元相對應(yīng)建立一個與人臉圖像大小相同的神經(jīng)網(wǎng)絡(luò),一幅二維圖像輸入到PCNN后會輸出一幅代表二值圖像的脈沖序列,這些圖像序列像攜帶了原始圖像中大量的有效特征信息,但它們的維度較大,不適宜直接進(jìn)行模式分類.如果對這些二值圖像依照迭代順序進(jìn)行求和運(yùn)算,則可將這些二值圖像信息轉(zhuǎn)化為一個一維的時間序列G(n),稱為振蕩時間序列,其數(shù)學(xué)表達(dá)式如下:

        時間序列將二維圖像信息轉(zhuǎn)化為了長度為n的一維向量信息,在實(shí)現(xiàn)圖像特征提取的同時也實(shí)現(xiàn)了數(shù)據(jù)降維.Johnson在實(shí)驗(yàn)中證明了每一幅圖像具有唯一性,這為基于圖像的時間序列進(jìn)行人臉識別提供了可能.

        圖3中的人臉圖像S10_1和S10_2屬于同一個人,而人臉S09_1為另一個人,由圖中可知PCNN的時間序列不僅實(shí)現(xiàn)了數(shù)據(jù)的降維,而且同一個人的不同圖像的時間序列表現(xiàn)出較好的一致性,而不同人臉圖像的時間序列曲線有很大的差異.

        圖3 人臉圖像及時間序列圖

        將同一個人的人臉圖像經(jīng)過縮放、旋轉(zhuǎn)幾何變換所對應(yīng)的時間序列如圖4所示,由圖中可知同一幅圖像經(jīng)過縮放、旋轉(zhuǎn)幾何變換后它們的時間序列形狀幾乎沒發(fā)生任何變化.由于時間序列是對脈沖數(shù)的統(tǒng)計,它的幅值與圖像的大小有很大的關(guān)系,所以當(dāng)對圖像進(jìn)行縮放時,時間序列的幅值發(fā)生了變換,但是形狀是相對穩(wěn)定的.這個實(shí)驗(yàn)說明了圖像的時間序列在對圖像進(jìn)行放大、縮小和旋轉(zhuǎn)具有很好的適應(yīng)性.

        4 人臉識別過程

        由前面幾節(jié)的分析可知,人臉圖像經(jīng)脈沖耦合神經(jīng)網(wǎng)絡(luò)提取的時間序列是人臉有效的鑒別特征.而脈沖耦合神經(jīng)網(wǎng)絡(luò)需要對整幅圖像的像素點(diǎn)進(jìn)行迭代處理,而小波變換在保持圖像主要能量的基礎(chǔ)上,不僅能達(dá)到對圖像數(shù)據(jù)的壓縮,加快后續(xù)過程的處理速度,而且還能模糊人臉表情的變化提高識別率.因此本文提出了如圖5所示的基于小波變換及脈沖耦合神經(jīng)網(wǎng)絡(luò)的人臉識別系統(tǒng).

        圖4 人臉圖像幾何變換時間序列圖

        圖5 人臉識別系統(tǒng)框圖

        圖5 所示的系統(tǒng)中,首先將人臉圖像的訓(xùn)練集和測試集圖片先經(jīng)過小波的濾波器組提取人臉圖像的低頻成分,然后送入脈沖耦合神經(jīng)網(wǎng)絡(luò)中提取相應(yīng)的時間序列,然后計算測試圖像和訓(xùn)練圖像的時間序列之間的歐式距離最小值所屬的類,該類即為測試的人臉圖像所屬的類.

        與傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)相比,脈沖耦合神經(jīng)網(wǎng)絡(luò)的神經(jīng)元不需要通過人臉圖像來進(jìn)行訓(xùn)練,而且脈沖耦合神經(jīng)網(wǎng)絡(luò)達(dá)到了人臉特征的提取和數(shù)據(jù)的降維作用.

        5 實(shí)驗(yàn)仿真及結(jié)果分析

        仿真實(shí)驗(yàn)平臺為PC和MATLAB2014b軟件編程實(shí)現(xiàn),在ORL人臉庫(包括40個人共200幅的256級灰度圖像,這些圖像具有不同的角度、縮放和面部表情變化)中選擇每個人的前n幅圖像作為訓(xùn)練集,剩下的10-n幅作為測試圖像進(jìn)行實(shí)驗(yàn).脈沖耦合神經(jīng)網(wǎng)絡(luò)神經(jīng)元的參數(shù)設(shè)置為 W=[0.5 1 0.5; 1 0.5 1; 0.5 1 0.5],α=0.1,β=0.257,VL=1,Vθ=1.1847,θ(0)=1.1,迭代次數(shù)N=25.

        為了考察本文算法與不經(jīng)過小波變換處理直接利用脈沖耦合神經(jīng)網(wǎng)絡(luò)進(jìn)行特征提取時間復(fù)雜度和識別率的差異,分別對兩種方法在訓(xùn)練樣本數(shù)為n=5的情況下進(jìn)行了實(shí)驗(yàn)比較,具體結(jié)果如表1所示.

        另外,本文取單個人臉訓(xùn)練樣本數(shù)為2、3、4、5,四種情況下將歐式距離改為余弦距離進(jìn)行分析比較,結(jié)果如表2所示.

        表1 本文算法與 PCNN 算法比較

        表2 余弦距離與歐式距離的比較

        從表2中的結(jié)果顯示可知,只有當(dāng)n=2時,余弦距離測度下人臉識別率要高于歐式距離的識別率,而在其它情況下利用歐式距離測度函數(shù)識別率都明顯高于余弦距離測度函數(shù).

        最后,在訓(xùn)練樣本為5的條件下,本文算法與傳統(tǒng)算法再識別率上進(jìn)行比較,其結(jié)果如表3所示,顯然本文算法的識別率要優(yōu)于表中所列.

        表3 本文算法與其它算法的比較

        實(shí)驗(yàn)結(jié)果表明,本文提出的算法相對于單獨(dú)使用PCNN不僅降低了時間成本,而且在識別率上有所提高.而且相對于某些傳統(tǒng)的人臉識別算法(PCA、DCT和LBP)表現(xiàn)出了更好的識別效果,驗(yàn)證了本文的算法的有效性.

        1毛亮,李立琛.淺談人臉識別技術(shù).中國公共安全,2014,(10):121–123.

        2Lu GC,Lin CY.PCA based immune networks for human face recognition.Applied Soft Computing,2011,11(2):1743–1752.[doi:10.1016/j.asoc.2010.05.017]

        3Shao H,Chen S,Zhao JY,et al.Face recognition based on subset selection via metric learning on manifold.Frontiers of Information Technology & Electronic Engineering,2015,16(12):1046–1058.

        4祁彥慶,汪烈軍,吳生武.一種基于稀疏表達(dá)和神經(jīng)網(wǎng)絡(luò)的人臉識別算法.計算機(jī)應(yīng)用與軟件,2016,33(10):172–175.[doi:10.3969/j.issn.1000-386x.2016.10.038]

        5Shchegoleva NL,Kukharev GA.Application of twodimensional principal component analysis for recognition of face images.Pattern Recognition and Image Analysis,2010,20(4):513–527.[doi:10.1134/S1054661810040127]

        6郭雅楠,李鴻燕.基于ICA和GSPSO-SVM的人臉識別方法.計算機(jī)工程與設(shè)計,2014,35(12):4302–4305,4310.[doi:10.3969/j.issn.1000-7024.2014.12.047]

        7張燕平,竇蓉蓉,趙姝,等.基于集成學(xué)習(xí)的規(guī)范化 LDA 人臉識別.計算機(jī)工程,2010,36(14):144–146.[doi:10.3969/j.issn.1000-3428.2010.14.052]

        8馬義德,李廉,王亞馥,等.脈沖耦合神經(jīng)網(wǎng)絡(luò)原理及其應(yīng)用.北京:科學(xué)出版社,2006.

        9聶仁燦,姚紹文,周冬明,等.基于簡化脈沖耦合神經(jīng)網(wǎng)絡(luò)的人臉識別.計算機(jī)科學(xué),2014,41(2):297–301.

        10趙光蘭,周冬明,趙東風(fēng),等.脈沖耦合神經(jīng)網(wǎng)絡(luò)的人臉識別方法.云南大學(xué)學(xué)報(自然科學(xué)版),2011,33(2):141–146.

        11岑翼剛,陳曉方,岑麗輝,等.基于單層小波變換的壓縮感知圖像處理.通信學(xué)報,2010,31(8A):52–55.

        12魏偉一,李戰(zhàn)明.基于改進(jìn)PCNN和互信息熵的自動圖像分割.計算機(jī)工程,2010,36(13):199–200.[doi:10.3969/j.issn.1000-3428.2010.13.071]

        13陳粟宋.一種融合PCNN和Watershed變換的圖像分割方法.武漢大學(xué)學(xué)報(工學(xué)版),2011,44(2):269–272.

        14嚴(yán)春滿,郭寶龍,易盟.基于改進(jìn)LP變換及自適應(yīng)PCNN的多聚焦圖像融合方法.控制與決策,2012,27(5):703–707,712.

        Face Recognition Based on Wavelet Transform and Pulse Coupled Neural Network

        HE Shi-Qiang,LIU Jin-Qing,LIU Yin,CAI Shu-Kuan,CHEN Cun-Di,ZHOU Xiao-Tong,DENG Shu-Min,WU Qing-Xiang

        (Key Aboratory of OptoElectronic Science and Technology for Medicine of Ministry of Education,College of Photonic and Electronic Engineering,Fujian Normal University,Fuzhou 350007,China)

        Pulse Coupled Neural Network (PCNN)is a new generation artificial neural network (ANN)based on biological vision.It has wide application prospects in the field of digital image processing and artificial intelligence.In this paper,we propose a method to extract face features by studying PCNN theoretical model and its working characteristics.Firstly,the low frequency feature of face image is extracted by wavelet transform.Then,the simplified PCNN is used to extract the corresponding time series of face image reconstructed by wavelet low-frequency coefficient,which is used as the feature sequence of face recognition.Finally,the face recognition process is completed with time series and Euclidean distance.In this paper,we demonstrate the effectiveness of the method with ORL face database.

        pulsed neural network; face recognition; wavelet transform; time series; Euclidean distance

        劉金清,E-mail:jqliu8208@fjnu.edu.cn

        何世強(qiáng),劉金清,劉引,蔡淑寬,陳存弟,周曉童,鄧淑敏,吳慶祥.基于簡化脈沖耦合神經(jīng)網(wǎng)絡(luò)的人臉識別.計算機(jī)系統(tǒng)應(yīng)用,2017,26(9):214–218.http://www.c-s-a.org.cn/1003-3254/5944.html

        國家自然科學(xué)基金(61179011); 福建教育廳項目(JAS151254); 福建師大項目(I201502019)

        2016-12-25; 修改時間:2017-01-23; 采用時間:2017-08-15

        猜你喜歡
        小波人臉識別人臉
        人臉識別 等
        有特點(diǎn)的人臉
        構(gòu)造Daubechies小波的一些注記
        揭開人臉識別的神秘面紗
        基于MATLAB的小波降噪研究
        電子制作(2019年13期)2020-01-14 03:15:32
        三國漫——人臉解鎖
        動漫星空(2018年9期)2018-10-26 01:17:14
        基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
        基于類獨(dú)立核稀疏表示的魯棒人臉識別
        馬面部與人臉相似度驚人
        長得象人臉的十種動物
        奇聞怪事(2014年5期)2014-05-13 21:43:01
        96免费精品视频在线观看| 色婷婷一区二区三区四区成人网| 女人被狂躁c到高潮| 久久国产人妻一区二区| 国産精品久久久久久久| 日韩av在线毛片| av网站一区二区三区| 亚洲男人天堂黄色av| 午夜无码国产理论在线| 亚洲最大成av人网站| 亚洲国产精品二区三区| 99久久国产精品网站| 在线人成免费视频69国产| 久久精品—区二区三区无码伊人色 | 成l人在线观看线路1| h在线国产| 精品国产夫妻自拍av| 亚洲国产婷婷六月丁香| 97精品人妻一区二区三区香蕉| 亚洲av人妖一区二区三区| 国产一区二区亚洲一区| 久爱www人成免费网站| 国农村精品国产自线拍| 亚洲AV无码久久久久调教| 久草手机视频在线观看| 亚洲一卡2卡3卡4卡5卡精品| 小12箩利洗澡无码视频网站| 日本一区二区三区激视频| 最新国产毛2卡3卡4卡| 天天影视色香欲综合久久| 亚洲区精品久久一区二区三区女同| 亚洲一区二区三区99| 久久精品国产亚洲av麻豆| 亚洲黄色尤物视频| 国产中文字幕亚洲国产| 99久久精品国产一区二区| 国产一区二区三区在线观看精品| 久久精品av一区二区免费| 日韩人妻不卡一区二区三区| 国产成人亚洲日韩欧美| 麻豆AV无码久久精品蜜桃久久|