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

    基于人體典型部位特征組合的行人檢測方法*

    2014-04-17 02:46:18葛平淑趙一兵趙宗艷
    汽車工程 2014年4期
    關鍵詞:腿部頭部行人

    郭 烈,葛平淑,趙一兵,趙宗艷

    (1.大連理工大學汽車工程學院,大連 116024;2.大連民族學院機電信息工程學院,大連 116600)

    前言

    為有效保護行人安全,提高汽車主動安全性能,研發(fā)行人檢測系統(tǒng)提醒駕駛員可能與鄰近的行人發(fā)生碰撞事故已成為減少道路交通中行人傷亡數(shù)量的重要手段[1]。

    在行人檢測過程中,需要定義一些特征來區(qū)分行人和背景,最直接的方法就是尋找圖像中的行人形狀特征對其進行識別,如梯度方向直方圖HOG特征[2]、Edgelet特征[3]和 Shapelet特征[4]等。由于每種特征的針對性不同,因此使用單一特征的檢測方法很難適應多場景下的行人檢測,國內外學者提出了多種特征信息融合的方法。文獻[5]中綜合使用HOG和Edgelet兩種局部特征,并通過SVM與Adaboost分別訓練各自的分類器,對紅外圖像中行人獲取較好的檢測效果;文獻[6]中提出了一種基于部位檢測和子結構組合的行人檢測方法,實現(xiàn)多姿態(tài)和部分遮擋的行人檢測;文獻[7]中分別在頭部和四肢等重點區(qū)域計算梯度方向直方圖,有效地減少了向量維數(shù),提高了檢測速度;文獻[8]中將基于局部特征的Adaboost級聯(lián)分類器和模板匹配相結合,實現(xiàn)視頻監(jiān)控領域的行人檢測。

    上述研究表明,多特征融合的方法可提高系統(tǒng)的檢測精度,但如何針對人體和人體典型部位提取有效特征,合理設計分類器結構,充分發(fā)揮各分類器的優(yōu)勢仍然是一個難題[9]。如果對于某些人體部位提取的特征并不能很好地區(qū)分行人與背景,那么這些無效的特征無疑會增加算法的負擔,影響實時性。為解決上述問題,提出一種基于部位特征組合的行人檢測方法。該方法對不同人體典型部位采用了不同的特征,針對各個部位的特點選擇最有效的特征來檢測,最后再將各自的檢測結果進行融合。

    1 算法描述

    1.1 基于HOG特征優(yōu)化的腿部檢測

    由于行人下半身所在的區(qū)域背景通常是路面圖像,而腿部占有較大面積,其輪廓邊緣比較明顯,在某些梯度方向和梯度強度上會出現(xiàn)有一定的峰值,其HOG特征與背景差別較大[10]。因此,提出針對行人腿部采用HOG特征進行檢測的方法。定義腿部區(qū)域檢測窗口大小為64×64像素,子窗口的寬度取為12、24、36、48、64 像素,子窗口的長寬比取為(1∶1)、(1∶2)和(2∶1),依據(jù)所設定的子窗口的寬度和長寬比便可以確定其數(shù)量和大小。通過設置多尺度的特征子窗口,能夠更加充分顯示腿部輪廓特征,各子窗口之間重疊的尺度為4、6、8個像素。采用Sobel算子從樣本集中提取所有的HOG特征,使用積分直方圖進行加速,最終從每個樣本中提取了1 386個HOG特征。

    經分析發(fā)現(xiàn),在眾多的特征中只有少部分特征對目標的分類起主導作用,其他大部分特征只含有少量的有用信息。通常這些冗余的特征不僅會增加算法的復雜性,使分類器訓練和檢測過程都非常耗時,而且會對最終的分類結果產生相反的作用[11]。因此,采用加權Fisher線性判別方法從中得到一部分能顯著區(qū)分正負樣本的HOG特征。

    加權Fisher線性判別能將樣本多維特征投影到一維直線上,投影方向能將各類樣本很好地分開,因此該投影過程具有一定的弱分類功能,同時實現(xiàn)特征降維[12]。投影方向的計算式為

    式中:wi為第i個樣本的權重;f(xi)為第i個樣本的特征為樣本類內均值為樣本類內加權協(xié)方差矩陣;a為最佳投影方向。

    為增強弱分類器的表達能力,采用查表的Gentle Adaboost算法[13]。在 Gentle AdaBoost中,弱分類器的輸出函數(shù)為

    式中:Pw(y=1|x)、Pw(y=-1|x)分別表示特征值為x的樣本時,其為正樣本或負樣本的加權概率。

    為提高分類器的精度,采用基于查表型弱分類器類型。查表的方法是將整個樣本集的特征值劃分為對應的n個不相交的子集,對于某個特征值為x的樣本,其弱分類器的輸出就表示為它在所屬子集中為正樣本與負樣本的加權概率差[14]。因此,查表型弱分類器的輸出為一個表示置信度的實數(shù)值,與離散型輸出相比,它包含了更多的信息并且更加適用于實際樣本的分布情況。查表型弱分類器的構造過程如下。

    經過加權Fisher線性判別后,HOG特征被降為1維,首先將它歸一化到[0,1]區(qū)間,將該區(qū)間均勻分成n個子區(qū)間:

    這樣,弱分類器可以定義為

    式中:W+1j和W-1j分別表示正、負樣本的特征值落在區(qū)間binj內的條件概率,表示為

    具體地說,令D+1j和D-1j分別表示第j個區(qū)間binj正負樣本權值的總和,即

    則條件概率為

    定義函數(shù)為

    則基于查表型Gentle AdaBoost弱分類器公式可表示為

    1.2 基于模板匹配的頭部檢測

    人體目標作為一種非剛體目標可能呈現(xiàn)出各種不同姿態(tài),很難用統(tǒng)一的模型對其進行描述。經過分析發(fā)現(xiàn),不管是行人的正面圖像還是側面圖像,其肩膀以上的頭部輪廓變動性最小,呈現(xiàn)出一定的形狀,只須定義不同大小的模板,就能滿足匹配的需要。另一方面,行人的頭部不易被遮擋,在檢測過程中很容易被檢測出來,因此模板匹配的方法更適用于行人頭部的檢測。針對頭部輪廓可以構建如圖1所示的“Ω”形位圖模板。對于64×128像素的行人樣本,定義頭部輪廓模板的大小為32×28像素。

    模板匹配可以通過距離變換將待檢測窗口轉換為DT圖像[15],其目的是將這種距離測度表示為模板變換參數(shù)的平滑函數(shù),一是方便快速搜索,二是提高匹配精度,三是通過對模板進行簡單縮放或者旋轉就能實現(xiàn)不同尺度、不同方向的匹配。因此在模板匹配之前,首先對待匹配圖像進行邊緣檢測,得到目標邊緣的二值圖像,然后對邊緣圖像進行距離轉換得到其DT圖像,匹配過程如圖2所示。

    將頭部模板與DT圖像做卷積,利用頭部模板圖像F中的特征點f所在的坐標相對應到DT圖像中的距離值相加,再求其平均值,即Chamfer距離,其計算公式為

    式中:DT(f)表示f對應于DT圖像的歐幾里得距離,取平均值的目的是為了對模板的大小進行歸一化。

    圖像越相似求得的Chamfer距離就越小,因此在待匹配的圖像中不包含頭部的區(qū)域,其Chamfer距離通常會遠大于有頭部的區(qū)域,模板匹配的原則是選取Chamfer距離最小的區(qū)域作為頭部最優(yōu)框的大小和位置。此外,當邊緣圖像的Chamfer距離值為檢測區(qū)域R中最小值時,須再做一次特征點數(shù)量上的判斷:

    式(11)可統(tǒng)計邊緣圖像的特征點數(shù)量,Threshold為特征點數(shù)判斷閾值,主要目的是確保不會因為特征點數(shù)稀少而降低平均距離。同時滿足上述兩個條件,才認為圖像匹配成功。

    1.3 部位特征組合的行人檢測

    在現(xiàn)實道路場景中,分類器的性能往往會受到復雜背景邊緣梯度的影響,進而造成誤判使得誤檢率增加,降低了算法精度。多特征融合的方法由于其檢測率較高、多場景適應性較好,目前應用最多。

    行人檢測過程中,分別對腿部和頭部輸出檢測結果之后,只有當頭部與腿部同時存在該區(qū)域中,才能認為所識別的目標是行人,這就要求對各部位都要有一個較為精確的判定。經過模板匹配后得到的是一個形似頭部的區(qū)域,須對其進行確認。根據(jù)行人頭部與背景的差異性特征,可將頭部輪廓的形狀特征轉化為數(shù)值特征并與支持向量機結合,訓練得到一個識別頭部的分類器。另外,通過部位相關性約束限制匹配區(qū)域也可以排除一部分虛警。最后將兩種特征的識別結果進行融合,正確判斷出目標為行人或非行人。

    通過觀察,行人頭部輪廓是一種描述目標全局的特征,它可以粗略感知行人是否存在,而HOG特征則偏重于對局部特征進行詳細的描述。特征組合就是建立在局部特征的基礎上,并且通過是否具有頭部全局特征來影響局部特征的判斷,可通過調整分類器的分類閾值實現(xiàn)。

    1.3.1 部位約束

    雖然模板匹配方法在一定程度上能夠檢測出人體的頭部,但是即使?jié)M足匹配條件的區(qū)域也會因為受到復雜背景的影響而出現(xiàn)誤判。為提高模板匹配的精度,通過人體部位之間的約束對匹配區(qū)域進行了限定。假設對于一個正確檢出的行人目標區(qū)域,行人通常位于圖像的正中央,而依據(jù)人體肢體的相對位置,頭部不可能在下半部分或者位于圖像的左、右邊緣。另外,由于對行人腿部的檢測在前,這也在一定程度上為后續(xù)的頭部檢測提供了必要的依據(jù),除了要避免多余背景的干擾,最主要的則是符合頭部在圖像中的實際位置。為此,建立如圖3所示的人體各部位約束區(qū)域,h、w分別為樣本高度和寬度,定義頭肩部位的高度為0.3h、軀干部位為0.48h、腿部為0.5h,將模板匹配范圍限定在腿部檢測框上方高為0.3h、寬為0.75w的矩形區(qū)域內。在限定了匹配區(qū)域后,就可以應用模板匹配方法對頭部進行檢測。

    1.3.2 頭部特征轉化

    經過模板匹配后,能獲得一個與行人頭部輪廓極為相似的區(qū)域,從而在形狀上找到兩者的區(qū)分。但卻很難通過單一的閾值去區(qū)分行人與非行人Chamfer距離間的差異,因此,可以考慮將行人頭部的這種形狀特征轉化為數(shù)值來描述。通過觀察,有頭部的區(qū)域輪廓往往類似于圓形,通過提取輪廓的類圓形特征,將形狀特征轉化為數(shù)值特征,進而對頭部匹配的結果在數(shù)值上做出精確判斷。

    針對行人頭部的圓形輪廓所提取的有效特征主要包括圓存在性[16]和邊界矩特征[17],并可將其組合共同描述行人頭部,從而將形狀特征轉化為數(shù)值特征。

    1.3.3 特征融合

    在通過上述圓形特征來判定模板匹配所產生的候選區(qū)域是否為頭部之后,最理想的情況是其判定結果與腿部一致,將目標確認為行人。然而實際上會出現(xiàn)兩個判定結果不一致的情況,因此,可以通過調整分類器的決策閾值,進而產生或寬松或嚴謹?shù)姆诸愋Ч麃韺Ψ诸惤Y果進行融合。

    圖4為閾值調整示意圖,當閾值取值為0時,特征點被判定為“-1”的非行人特征,在放寬閾值之后,則該特征點被判定為“+1”的行人特征,如圖4(a)所示;而增加閾值后,則將特征點劃分得比較嚴謹,行人特征視為非行人特征,如圖4(b)所示。

    頭部檢測的目的是排除腿部檢測產生的虛警,因此,閾值調整策略主要以頭部輪廓特征是否存在為主要依據(jù)。當頭部分類器沒有檢測到頭部時,應當采用較為嚴謹?shù)拈撝嫡{整策略,提高腿部分類器的決策值,進而重新對檢測出的腿部區(qū)域進行判定。但是此方法的前提是頭部區(qū)域判定十分精確。由于受邊緣提取算法的影響,信息的缺失也會使頭部產生誤判,為此需要對腿部分類器的決策值進一步分析。腿部檢測器的決策值越高則其為行人的可能性越大,若腿部分類器的決策值大于所有弱分類器權值的一半,則不進行調整,直接判斷腿部分類器分類正確,目標為行人。若腿部分類器的決策值介于0與所有弱分類器權值的一半時,說明腿部特征的決策值在一個模棱兩可的區(qū)域內,因此就要調高分類平面。經過調整之后,再對檢測出的腿部區(qū)域按照新的分類平面進行重新判斷,若判斷結果為假,則確定目標為非行人。若腿部判定結果仍然為真,此時則反過來對頭部分類平面進行寬松的調整,重新對頭部做一次判定,如果頭部判定結果為真,則可以確定目標為行人。最終,仍然存在一類情況,即腿部判定結果為真而頭部判定結果為假,由于特征組合是建立在局部特征存在的基礎上,此時判斷標準以腿部結果為主,因此,直接判定檢測結果為行人。

    設腿部分類器閾值調整前后的決策值分別為V1和V'1,頭部分類器閾值調整前后的決策值分別為 V2和 V'2,初始的分類閾值 Threshold1=0,Threshold2=0,弱分類器的權值之和為α,根據(jù)經驗確定腿部閾值調整的區(qū)間為2,頭部閾值調整的區(qū)間為0.5。部位組合閾值調整的具體步驟如下。

    2 試驗與分析

    2.1 數(shù)據(jù)集

    為驗證所提部位特征組合行人檢測算法的有效性,在Matlab9.0環(huán)境下進行了性能對比分析,采用Intel Core(TM)2的CPU和2G內存的計算機。所使用的試驗數(shù)據(jù)集共有兩個。數(shù)據(jù)集1包含兩個行人樣本庫,一是MIT行人數(shù)據(jù)庫,二是課題組自行采集的樣本集合,多為城市交通環(huán)境下的圖像。由于MIT數(shù)據(jù)庫與自行采集的數(shù)據(jù)庫的行人樣本比較簡單,將其合并為一個,共得到2 324個行人樣本,2 300個非行人樣本。分別取出其中的1 500個用于訓練,其余均為測試樣本。數(shù)據(jù)集2由INRIA行人庫組成,挑選出INRIA數(shù)據(jù)庫中腿部不存在嚴重遮擋的行人樣本,最終訓練集挑選出1 970張,測試集共890張。對于所有樣本,截取下半部分圖像并歸一化為64×64像素大小。

    2.2 腿部分類器測試結果

    通過對兩個試驗數(shù)據(jù)集進行訓練和測試,評價算法的標準采用檢測效果和檢測時間,檢測效果用ROC曲線來描述。

    由于對Adaboost中的弱分類器進行了改進,一是采用加權Fisher線性判別來代替線性SVM,從而降低HOG特征維數(shù),二是使用查表型弱分類器替代傳統(tǒng)閾值型弱分類器。因此,在試驗研究中比較了基于以下3種弱分類器訓練得到的腿部分類器在使用特征的數(shù)量上和檢測精度上的差別:線性SVM、加權Fisher線性判別的閾值型弱分類器和加權Fisher線性判別的查表型弱分類器。

    圖5比較了由兩個訓練集訓練得到的強分類器性能,訓練集1中線性SVM所使用的弱分類器達到了49個,而加權Fisher線性判別所使用的弱分類器只有15個,迭代次數(shù)大大減少。訓練集2中由于循環(huán)條件的限制線性SVM檢測率沒有達到100%,同樣加權Fisher線性判別所使用的弱分類器減少到了38個。加權Fisher線性判別不僅解決了每次迭代隨機選取部分HOG特征的問題,還能提高其分類速度。從圖5中還可看出,基于查表型弱分類器的收斂速度相對較快,使用較少的特征就能達到目標檢測率。

    圖6 比較了不同類型弱分類器形成的強分類器在兩個測試集上的檢測性能,即ROC曲線??梢钥闯?,使用查表型弱分類器后的Adaboost檢測器的精度與閾值型弱分類器相比,精度有所提高,能更進一步優(yōu)化強分類器,提升檢測性能。而由于線性SVM特征數(shù)較多,故檢測率比其他兩種方法稍高。

    表1列出了3種不同分類器在不同數(shù)據(jù)集上的時間消耗和分類性能。由表可見,加權Fisher與查表型弱分類器的分類性能分別達到了98.64%和93.32%,要優(yōu)于線性SVM與加權Fisher的閾值型分類器。在時間消耗方面,與線性SVM相比,使用加權Fisher判別的兩種弱分類器對數(shù)據(jù)集1的訓練時間分別減少了86.37%和86.43%,驗證了其在運算速度上的巨大優(yōu)勢。單獨分析查表型弱分類器,其檢測時間比簡單閾值型弱分類器稍長,主要是由于查表的方法在實際應用過程中,不僅須記錄查表的每個區(qū)間系數(shù),還須記錄特征的最大值和最小值。這些都會占用存儲空間,并增加一些運算量。

    表1 不同弱分類器時間和性能比較

    2.3 部位特征組合檢測試驗結果

    為驗證部位組合檢測方法的有效性,對樣本集中的圖像進行了處理,采用混淆矩陣來評價算法的檢測結果,主要比較了本文中所提出的部位組合算法、傳統(tǒng) HOG算法[2]和用積分圖像計算特征的HOG簡化算法。在樣本集的檢測性能,由于采用了積分圖像提取梯度方向的直方圖特征,省略了復雜的三維線性插值步驟,大大縮短了特征提取的時間。圖7和圖8分別給出這3種不同方法在兩個測試集的檢測結果。

    從圖7看出,傳統(tǒng)HOG算法在734個樣本中有4個行人誤認為非行人,檢測率達到了99.4%,檢測結果近乎完美;特征組合算法有24個行人誤認為非行人,檢測率也達到了96.7%,與傳統(tǒng)HOG算法相比略有下降;而HOG簡化算法檢測率只有76.1%,可見省略三維線性插值步驟降低了算法的精度,通過對特征和分類器進行一系列改進后檢測精度得到明顯提高。

    從圖8看出,特征組合算法的檢測率比HOG簡化算法還低,這是因為測試集2中的行人頭部所處的環(huán)境比較復雜,存在很多遮擋情況,在一定程度上影響了系統(tǒng)的性能。另一方面,采用特征組合算法的檢測率主要取決于各部位檢測率的乘積,因此會低于整體的檢測率。單獨分析誤檢率可以看到,特征組合算法比HOG簡化算法要好,進一步驗證了采用特征組合算法能夠有效排除虛警的優(yōu)越性。

    圖9示出特征組合算法對不同道路環(huán)境圖像的檢測結果。上面兩張圖為對INRIA數(shù)據(jù)庫中圖像的測試結果,下面兩張圖為室外拍攝的道路圖像測試結果。從檢測結果可以看出,特征組合算法對不同道路場景下的行人都能取得較好的檢測效果,對于在單人和多人的情況同樣適用。但是,試驗中也發(fā)現(xiàn)一些問題,道路兩側的背景復雜使得算法出現(xiàn)誤判,將道路兩側的樹木誤認為是行人目標,復雜背景對邊緣提取算法產生影響,使行人頭部定位不準確。

    3 結論

    針對人體的不同部位使用不同的特征進行檢測,再通過一定的方法將檢測結果融合,實現(xiàn)基于人體典型部位特征組合的行人檢測。部位特征組合是建立在人體部位約束、頭部確認與分類閾值調整三者結合的基礎上,三者之間存在著逐級遞進關系。部位約束不但減少了模板匹配的范圍,提高了實時性,而且提高了匹配精度。通過頭部新特征的提取將頭部形狀特征轉化為數(shù)值特征,并結合支持向量機對頭部進行最終確認,從而有效排除虛警。根據(jù)不同人體部位融合的特點,通過簡單的閾值調整便可實現(xiàn),降低算法的復雜性。試驗表明,部位特征組合的方法能夠有效地檢測出現(xiàn)實場景中的行人,并能排除一定的虛警目標。

    試驗中也發(fā)現(xiàn)頭部定位不精確和復雜背景造成的行人誤識別等問題。頭部識別的關鍵在于獲取區(qū)分度較高的特征,為提高頭部檢測的準確率,須進一步研究更加可靠、穩(wěn)定的頭部目標特征。此外,還可以將級聯(lián)分類機制引入到分類器中以提高系統(tǒng)的實時性。

    [1] Zegeer CV,Bushell M.Pedestrian Crash Trends and Potential Countermeasures from Around the World[J].Accident Analysis and Prevention,2012,44(1):3 -11.

    [2] Dalal N,Triggs B.Histograms of Oriented Gradients for Human Detection[C].IEEE Conf.on Computer Vision and Pattern Recognition,San Diego,CA,USA:IEEE,2005,1:886 -893.

    [3] Wu B,Nevatia R.Detection and Tracking of Multiple,Partially Occluded Humans by Bayesian Combination of Edgelet Based Part Detectors[J].Int.Journal of Computer Vision,2007,75(2):247-266.

    [4] Sabzmeydani P,Mori G.Detecting Pedestrians by Learning Shapelet Features[C].IEEE Conf.on Computer Vision and Pattern Recognition,Minneapolis,Minnesota,USA:IEEE,2007.

    [5] Zhang L,Wu B,Nevatia R.Pedestrian Detection in Infrared Images Based on Local Shape Features[C].IEEE Conf.Computer Vision and Pattern Recognition,Minneapolis,Minnesota,USA:IEEE,2007.

    [6] 胡斌,王生進,丁曉青.基于部位檢測和子結構組合的行人檢測方法[J].計算機科學,2009,36(11):242 -246.

    [7] 曾春,李曉華,周激流.基于感興趣區(qū)梯度方向直方圖的行人檢測[J].計算機工程,2009,35(24):182 -184.

    [8] 周晨卉,王生進,丁曉青.基于局部特征級聯(lián)分類器和模板匹配的行人檢測[J].中國圖象圖形學報,2010,15(5):824 -829.

    [9] Dollár P,Wojek C,Schiele B,et al.Pedestrian Detection:An E-valuation of the State of the Art[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2012,34(4):743 -761.

    [10] Enzweiler M,Gavrila D M.Monocular Pedestrian Detection:Survey and Experiments[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2009,31(12):2179 -2195.

    [11] Paisitkriangkcrai S,Shen C H,Zhang J.Fast Pedestrian Detection Using a Cascade of Boosted Covariance Features[J].IEEE Transactions on Circuits and Systems for Video Technology,2008,18(8):1140-1151.

    [12] Laptev I.Improvements of Object Detection Using Boosted Histograms[J].Journal of Image and Vision Computing,2009,27:535-544.

    [13] Freund Y,Schapire E R.A Decision-Theoretic Generalization of on Line Learning and an Application to Boosting[C].Proceedings of the Second European Conference on Computational Learning Theory,Barcelona,Spain:Springer,1995:23 -37.

    [14] Wu B,Ai H Z,Huang C.LUT-based Adaboost for Gender Classification[C].International Conference on Audio and Video-Based Biometric Person Authentication,Guildford,UK:Springer,2003:104 -110.

    [15] Gavrila D M.A Bayesian,Exemplar-based Approach to Hierarchical Shape Matching[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2007,29(8):1408 -1421.

    [16] 趙敏.單目視覺多行人目標檢測與跟蹤技術研究[D].重慶:重慶大學,2010.

    [17] Chen C C.Improved Moment Invariant for Shape Discrimination[J].Pattern Recognition,1993,26(5):683 -686.

    猜你喜歡
    腿部頭部行人
    分體式腿部氣動按摩裝置設計
    玩具世界(2022年2期)2022-06-15 07:35:58
    頭部按摩治療老伴失憶
    毒舌出沒,行人避讓
    意林(2021年5期)2021-04-18 12:21:17
    火箭的頭部為什么是圓鈍形?
    軍事文摘(2020年22期)2021-01-04 02:16:38
    路不為尋找者而設
    揚子江(2019年1期)2019-03-08 02:52:34
    我是行人
    自適應統(tǒng)計迭代重建算法在頭部低劑量CT掃描中的應用
    腳下放靠墊睡覺更解乏
    腳下放靠墊睡覺更解乏
    WONCA研究論文摘要匯編:青霉素對腿部疏松結締組織炎反復發(fā)作的預防作用
    少妇的逼水好多| 不卡视频在线观看欧美| 内地一区二区视频在线| 特大巨黑吊av在线直播| 美女被艹到高潮喷水动态| 人妻制服诱惑在线中文字幕| 草草在线视频免费看| 91午夜精品亚洲一区二区三区 | 乱系列少妇在线播放| 干丝袜人妻中文字幕| 亚洲不卡免费看| 久久人妻av系列| 中文资源天堂在线| 日韩欧美 国产精品| 搡老岳熟女国产| 99久久精品热视频| 国产一区二区三区在线臀色熟女| 午夜老司机福利剧场| 麻豆av噜噜一区二区三区| 日韩欧美国产在线观看| 欧美精品啪啪一区二区三区| 亚洲 国产 在线| 国产免费av片在线观看野外av| 欧美一区二区国产精品久久精品| 成年人黄色毛片网站| 久久精品国产99精品国产亚洲性色| 国产单亲对白刺激| 国产精品亚洲美女久久久| 91狼人影院| 免费观看在线日韩| 91麻豆精品激情在线观看国产| 露出奶头的视频| 最后的刺客免费高清国语| 久久精品国产亚洲网站| 免费在线观看日本一区| 日韩一区二区视频免费看| 婷婷亚洲欧美| 精品久久久久久,| 精品乱码久久久久久99久播| 一个人免费在线观看电影| 嫩草影院入口| 免费看日本二区| 99热精品在线国产| 亚洲电影在线观看av| 无遮挡黄片免费观看| 99热只有精品国产| 偷拍熟女少妇极品色| 我要看日韩黄色一级片| 高清日韩中文字幕在线| 免费电影在线观看免费观看| 美女大奶头视频| 久久久久久伊人网av| 国产精品自产拍在线观看55亚洲| 午夜爱爱视频在线播放| 直男gayav资源| 国产亚洲精品综合一区在线观看| 国内少妇人妻偷人精品xxx网站| 热99在线观看视频| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩高清专用| 淫妇啪啪啪对白视频| 在线观看66精品国产| 五月伊人婷婷丁香| 欧美激情在线99| 成人鲁丝片一二三区免费| 99久久精品一区二区三区| 身体一侧抽搐| 十八禁网站免费在线| 欧美在线一区亚洲| 99久国产av精品| 尤物成人国产欧美一区二区三区| 免费无遮挡裸体视频| 欧美三级亚洲精品| 草草在线视频免费看| 深爱激情五月婷婷| 久久人人精品亚洲av| 亚洲av成人av| 桃红色精品国产亚洲av| 国产av一区在线观看免费| 韩国av一区二区三区四区| 身体一侧抽搐| eeuss影院久久| 久久久久久伊人网av| 在线观看美女被高潮喷水网站| 男人和女人高潮做爰伦理| 蜜桃亚洲精品一区二区三区| 亚洲性久久影院| 国产精品综合久久久久久久免费| 天堂动漫精品| 国产精华一区二区三区| 99久久成人亚洲精品观看| 不卡视频在线观看欧美| 欧美日韩精品成人综合77777| 亚洲狠狠婷婷综合久久图片| 久久精品国产亚洲网站| 天堂av国产一区二区熟女人妻| 亚洲成人免费电影在线观看| 午夜a级毛片| 一个人观看的视频www高清免费观看| 亚洲va在线va天堂va国产| 国产精品1区2区在线观看.| 三级男女做爰猛烈吃奶摸视频| 亚洲自拍偷在线| 国产精品爽爽va在线观看网站| 在现免费观看毛片| 亚洲在线观看片| 在线a可以看的网站| 亚洲精品成人久久久久久| 中文字幕熟女人妻在线| 国产人妻一区二区三区在| 亚洲天堂国产精品一区在线| 99在线视频只有这里精品首页| 亚洲男人的天堂狠狠| 婷婷亚洲欧美| 色吧在线观看| 在线免费观看不下载黄p国产 | 国产精品爽爽va在线观看网站| 日本撒尿小便嘘嘘汇集6| 国产午夜精品论理片| 午夜影院日韩av| 偷拍熟女少妇极品色| 免费看日本二区| 国内精品一区二区在线观看| 一区福利在线观看| 国产亚洲av嫩草精品影院| av在线老鸭窝| 波多野结衣高清无吗| 亚洲狠狠婷婷综合久久图片| 久久精品国产99精品国产亚洲性色| 国产精品乱码一区二三区的特点| 免费av毛片视频| 久久国内精品自在自线图片| 色综合婷婷激情| 国产午夜精品久久久久久一区二区三区 | 亚洲精华国产精华液的使用体验 | 国产女主播在线喷水免费视频网站 | 噜噜噜噜噜久久久久久91| 很黄的视频免费| 韩国av一区二区三区四区| 97人妻精品一区二区三区麻豆| 女人十人毛片免费观看3o分钟| 国产精品久久电影中文字幕| 日韩欧美 国产精品| 亚洲性夜色夜夜综合| 国产单亲对白刺激| 成年女人看的毛片在线观看| 国内毛片毛片毛片毛片毛片| 国产大屁股一区二区在线视频| 国产午夜福利久久久久久| 色哟哟·www| 一级黄片播放器| 搡老熟女国产l中国老女人| 美女大奶头视频| 国产精品久久视频播放| 国产欧美日韩精品一区二区| 久久国产乱子免费精品| 亚洲四区av| a级一级毛片免费在线观看| 国产av在哪里看| 日本a在线网址| 精品久久久久久成人av| 成人av一区二区三区在线看| 国产精品女同一区二区软件 | 99热只有精品国产| 又爽又黄无遮挡网站| 99久久九九国产精品国产免费| 国产蜜桃级精品一区二区三区| 国模一区二区三区四区视频| 美女被艹到高潮喷水动态| 欧美又色又爽又黄视频| 一个人观看的视频www高清免费观看| 真实男女啪啪啪动态图| 日韩精品青青久久久久久| 欧美一区二区亚洲| 老师上课跳d突然被开到最大视频| 国产女主播在线喷水免费视频网站 | 在线天堂最新版资源| 丰满的人妻完整版| 亚洲av熟女| 热99re8久久精品国产| 欧美潮喷喷水| 在线观看免费视频日本深夜| 国产精品一区二区三区四区久久| 欧美黑人巨大hd| 成人亚洲精品av一区二区| 婷婷亚洲欧美| 丰满的人妻完整版| 久久久久九九精品影院| 最近最新中文字幕大全电影3| 黄色日韩在线| 久久亚洲精品不卡| 丝袜美腿在线中文| 91麻豆av在线| 国产免费男女视频| 日韩欧美国产一区二区入口| 国产蜜桃级精品一区二区三区| 日日干狠狠操夜夜爽| 最近最新免费中文字幕在线| 久久国内精品自在自线图片| 高清在线国产一区| 亚洲欧美日韩东京热| 最好的美女福利视频网| 成年女人永久免费观看视频| 美女cb高潮喷水在线观看| 色尼玛亚洲综合影院| 欧美+日韩+精品| 日日啪夜夜撸| 欧美性猛交黑人性爽| 国产成人av教育| 欧美最黄视频在线播放免费| 国产真实伦视频高清在线观看 | 国产黄片美女视频| 91久久精品国产一区二区三区| 无人区码免费观看不卡| 国内精品一区二区在线观看| 久久久久久久午夜电影| 在线播放国产精品三级| 2021天堂中文幕一二区在线观| 性插视频无遮挡在线免费观看| 黄色欧美视频在线观看| 精品久久久久久久久亚洲 | 久久久午夜欧美精品| 69av精品久久久久久| 久久久精品大字幕| 欧美bdsm另类| 国产淫片久久久久久久久| 欧美日韩国产亚洲二区| 十八禁网站免费在线| 美女高潮的动态| 欧美日韩瑟瑟在线播放| 九九热线精品视视频播放| 国产人妻一区二区三区在| 九九爱精品视频在线观看| 在线免费观看的www视频| 琪琪午夜伦伦电影理论片6080| 国产中年淑女户外野战色| 欧美日本视频| 国产一区二区三区av在线 | 黄色配什么色好看| 精华霜和精华液先用哪个| 精品久久久久久久末码| 精品久久久久久久久久久久久| 免费无遮挡裸体视频| 久久精品国产清高在天天线| 久久国内精品自在自线图片| 亚洲国产精品久久男人天堂| 亚洲专区国产一区二区| 我的女老师完整版在线观看| 免费看av在线观看网站| 一进一出抽搐gif免费好疼| 亚洲av第一区精品v没综合| 99热这里只有是精品50| 搡女人真爽免费视频火全软件 | 女人被狂操c到高潮| 中文亚洲av片在线观看爽| 听说在线观看完整版免费高清| 亚洲avbb在线观看| 国产av在哪里看| 日韩亚洲欧美综合| 国产视频一区二区在线看| 亚洲狠狠婷婷综合久久图片| 99热精品在线国产| 国产高清视频在线观看网站| av在线天堂中文字幕| 51国产日韩欧美| 黄色日韩在线| 欧美性猛交╳xxx乱大交人| 国产在线精品亚洲第一网站| 婷婷六月久久综合丁香| 成熟少妇高潮喷水视频| 色在线成人网| 日本欧美国产在线视频| 国产精品1区2区在线观看.| 亚洲性夜色夜夜综合| 亚洲三级黄色毛片| 国产私拍福利视频在线观看| 亚洲一区高清亚洲精品| 人人妻,人人澡人人爽秒播| 成年人黄色毛片网站| 婷婷亚洲欧美| 国内毛片毛片毛片毛片毛片| 国产亚洲精品av在线| 波多野结衣高清作品| 人妻久久中文字幕网| 神马国产精品三级电影在线观看| 亚洲国产色片| 日本黄大片高清| 亚洲aⅴ乱码一区二区在线播放| a级一级毛片免费在线观看| 丰满的人妻完整版| 99热这里只有精品一区| 日本五十路高清| 禁无遮挡网站| 国产乱人伦免费视频| 久久久久久久午夜电影| 亚洲av成人精品一区久久| 搡老岳熟女国产| 欧美xxxx性猛交bbbb| 在线观看66精品国产| 我要看日韩黄色一级片| 日日夜夜操网爽| 少妇被粗大猛烈的视频| 成人二区视频| 午夜福利在线在线| 亚洲欧美日韩高清专用| 久久香蕉精品热| 校园人妻丝袜中文字幕| 久久久久免费精品人妻一区二区| 亚洲国产欧美人成| 一本精品99久久精品77| 久久精品国产亚洲av香蕉五月| 国产中年淑女户外野战色| 久久精品91蜜桃| 亚洲七黄色美女视频| 可以在线观看毛片的网站| 99久久精品热视频| 蜜桃久久精品国产亚洲av| 亚洲真实伦在线观看| 国产亚洲精品久久久com| 国内精品一区二区在线观看| 一个人免费在线观看电影| 小蜜桃在线观看免费完整版高清| 成年女人永久免费观看视频| 麻豆一二三区av精品| 97超视频在线观看视频| 老女人水多毛片| 一a级毛片在线观看| 欧美丝袜亚洲另类 | 久久久国产成人精品二区| 22中文网久久字幕| 最近在线观看免费完整版| 精品不卡国产一区二区三区| 欧美成人a在线观看| 日本熟妇午夜| 日韩欧美免费精品| 熟妇人妻久久中文字幕3abv| 一级毛片久久久久久久久女| 国产午夜精品论理片| 舔av片在线| 99热精品在线国产| 床上黄色一级片| 亚洲在线自拍视频| 99在线人妻在线中文字幕| 久久精品国产99精品国产亚洲性色| 国产黄片美女视频| 欧美激情国产日韩精品一区| 日韩欧美国产一区二区入口| 男人的好看免费观看在线视频| 高清毛片免费观看视频网站| 国产av一区在线观看免费| 亚洲中文字幕日韩| 国产精品久久久久久亚洲av鲁大| 成人午夜高清在线视频| av专区在线播放| 国产高清三级在线| 亚洲av成人精品一区久久| 少妇丰满av| 午夜爱爱视频在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 精品午夜福利在线看| 欧美极品一区二区三区四区| 日韩人妻高清精品专区| 精品久久久久久久末码| 国产欧美日韩精品一区二区| 日本免费一区二区三区高清不卡| 最近视频中文字幕2019在线8| 亚洲人成网站在线播| 国产极品精品免费视频能看的| or卡值多少钱| 99热这里只有精品一区| 国产精品一区二区性色av| 一个人免费在线观看电影| 亚洲天堂国产精品一区在线| 婷婷精品国产亚洲av在线| 草草在线视频免费看| 春色校园在线视频观看| 日韩在线高清观看一区二区三区 | 国产精品一区二区三区四区久久| 精品久久久久久久久久免费视频| 日本精品一区二区三区蜜桃| 天堂网av新在线| 在线免费观看的www视频| 精品久久久久久久久久久久久| 国产伦精品一区二区三区视频9| 极品教师在线视频| 久久99热这里只有精品18| 欧美国产日韩亚洲一区| 又粗又爽又猛毛片免费看| 日韩高清综合在线| 小蜜桃在线观看免费完整版高清| 成人三级黄色视频| 99九九线精品视频在线观看视频| 99精品在免费线老司机午夜| 久久精品国产亚洲av香蕉五月| 成年女人毛片免费观看观看9| 国产高清视频在线播放一区| 亚洲性夜色夜夜综合| av专区在线播放| 88av欧美| 最近最新中文字幕大全电影3| 欧美bdsm另类| 成人二区视频| 欧美日本亚洲视频在线播放| 国产人妻一区二区三区在| 久久香蕉精品热| 亚洲黑人精品在线| 日日摸夜夜添夜夜添av毛片 | 真实男女啪啪啪动态图| 搡老岳熟女国产| 国产精品美女特级片免费视频播放器| 亚洲专区国产一区二区| 国内毛片毛片毛片毛片毛片| 日本三级黄在线观看| 精品人妻1区二区| 99riav亚洲国产免费| 又黄又爽又免费观看的视频| 亚洲最大成人av| 高清在线国产一区| 亚洲精品亚洲一区二区| 超碰av人人做人人爽久久| 亚洲不卡免费看| 两个人的视频大全免费| 亚洲最大成人av| 国产久久久一区二区三区| 夜夜看夜夜爽夜夜摸| 国产亚洲91精品色在线| 韩国av在线不卡| 国产老妇女一区| 乱人视频在线观看| 蜜桃久久精品国产亚洲av| 12—13女人毛片做爰片一| 黄色女人牲交| 欧美黑人巨大hd| av在线老鸭窝| 亚洲精品一区av在线观看| 啪啪无遮挡十八禁网站| 成人美女网站在线观看视频| 精品人妻1区二区| 中国美女看黄片| 看片在线看免费视频| 美女黄网站色视频| 午夜爱爱视频在线播放| 伦理电影大哥的女人| 久久精品影院6| 亚洲男人的天堂狠狠| 国产av在哪里看| 欧美色视频一区免费| 三级男女做爰猛烈吃奶摸视频| 偷拍熟女少妇极品色| 可以在线观看毛片的网站| 我的女老师完整版在线观看| 亚洲熟妇中文字幕五十中出| 国产综合懂色| 国产一区二区在线av高清观看| 日韩欧美在线乱码| 日日干狠狠操夜夜爽| av天堂中文字幕网| 亚洲国产色片| 午夜福利18| 精品久久久久久成人av| 亚洲自偷自拍三级| 中文字幕高清在线视频| 丰满乱子伦码专区| 联通29元200g的流量卡| 黄片wwwwww| 久9热在线精品视频| 熟女电影av网| 亚洲五月天丁香| 成人国产综合亚洲| 在线观看美女被高潮喷水网站| 久久久精品大字幕| 午夜免费成人在线视频| bbb黄色大片| 22中文网久久字幕| 人妻少妇偷人精品九色| 欧美日韩综合久久久久久 | 男人的好看免费观看在线视频| 亚洲成av人片在线播放无| 国产精品免费一区二区三区在线| 国产v大片淫在线免费观看| 成人国产一区最新在线观看| 99riav亚洲国产免费| 久久久久国产精品人妻aⅴ院| 国产精品自产拍在线观看55亚洲| 国产色爽女视频免费观看| 少妇人妻一区二区三区视频| 欧美成人性av电影在线观看| 久久精品国产亚洲av天美| 欧美三级亚洲精品| 日本 欧美在线| 久久久久国产精品人妻aⅴ院| 国产精品女同一区二区软件 | 国产私拍福利视频在线观看| 窝窝影院91人妻| 我要看日韩黄色一级片| h日本视频在线播放| 在现免费观看毛片| 2021天堂中文幕一二区在线观| 99热这里只有精品一区| 长腿黑丝高跟| 久久精品国产亚洲av涩爱 | а√天堂www在线а√下载| av天堂中文字幕网| 能在线免费观看的黄片| 美女免费视频网站| 又黄又爽又刺激的免费视频.| 国产一区二区激情短视频| 搞女人的毛片| 午夜老司机福利剧场| 亚洲专区中文字幕在线| 午夜福利成人在线免费观看| 亚洲熟妇中文字幕五十中出| 亚洲,欧美,日韩| 亚洲精品久久国产高清桃花| 亚洲精品亚洲一区二区| 精品久久久久久成人av| 亚洲熟妇熟女久久| 久久久久久久久久成人| 亚洲欧美清纯卡通| 婷婷精品国产亚洲av在线| 99热这里只有精品一区| 少妇的逼好多水| 欧美日韩国产亚洲二区| 国产成人影院久久av| 此物有八面人人有两片| 九九爱精品视频在线观看| 长腿黑丝高跟| 日韩国内少妇激情av| 久久久久国内视频| 99国产极品粉嫩在线观看| 在线播放国产精品三级| 综合色av麻豆| 国语自产精品视频在线第100页| 国产精品嫩草影院av在线观看 | av女优亚洲男人天堂| 五月伊人婷婷丁香| 窝窝影院91人妻| 免费人成视频x8x8入口观看| av国产免费在线观看| 99热这里只有精品一区| 国产成人福利小说| 成人国产综合亚洲| 99精品久久久久人妻精品| 一级a爱片免费观看的视频| 又黄又爽又刺激的免费视频.| 听说在线观看完整版免费高清| 亚洲人成网站在线播放欧美日韩| 少妇人妻精品综合一区二区 | av福利片在线观看| 美女 人体艺术 gogo| 久久精品国产自在天天线| 午夜精品一区二区三区免费看| 观看美女的网站| 三级国产精品欧美在线观看| 99久国产av精品| 99热只有精品国产| 中文资源天堂在线| 亚洲av免费高清在线观看| 亚洲欧美日韩无卡精品| 淫妇啪啪啪对白视频| 免费观看人在逋| 国产免费一级a男人的天堂| 成人鲁丝片一二三区免费| 免费不卡的大黄色大毛片视频在线观看 | 在线播放无遮挡| 最后的刺客免费高清国语| 人人妻,人人澡人人爽秒播| 级片在线观看| 国内毛片毛片毛片毛片毛片| 动漫黄色视频在线观看| 美女高潮的动态| 深夜精品福利| 成人国产综合亚洲| 亚洲人与动物交配视频| 国产蜜桃级精品一区二区三区| 国产精品人妻久久久久久| 一区二区三区免费毛片| 日本五十路高清| 国产aⅴ精品一区二区三区波| 中文字幕熟女人妻在线| 在线天堂最新版资源| 美女高潮的动态| 亚州av有码| 国产探花极品一区二区| 真人一进一出gif抽搐免费| 美女免费视频网站| 99热精品在线国产| 91久久精品国产一区二区成人| 老熟妇仑乱视频hdxx| 久久精品国产亚洲av天美| 亚洲狠狠婷婷综合久久图片| 国产视频一区二区在线看| 97超级碰碰碰精品色视频在线观看| 国产精品一及| 亚洲,欧美,日韩| 国产精品嫩草影院av在线观看 | 亚洲人成伊人成综合网2020| 国产精品一及| 久久精品国产亚洲av涩爱 | 久久午夜福利片| av天堂在线播放| av在线老鸭窝| 搡老妇女老女人老熟妇| 悠悠久久av| 亚洲图色成人| 99九九线精品视频在线观看视频| 亚洲国产高清在线一区二区三| 欧美3d第一页| 国产欧美日韩精品亚洲av| 伦理电影大哥的女人| 非洲黑人性xxxx精品又粗又长| 免费电影在线观看免费观看| 日本熟妇午夜| 在线免费十八禁| 黄色女人牲交|