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

    肺部聽診音數(shù)據(jù)庫建庫技術及方法研究

    2024-12-31 00:00:00張冬瑩葉培韜李洽勝簡文華梁振宇鄭勁平
    中國全科醫(yī)學 2024年36期
    關鍵詞:肺疾病支持向量機數(shù)據(jù)分析

    【摘要】 當前無論是物理聽診器亦或是電子聽診器的肺音聽診結果仍然主要依靠醫(yī)生專業(yè)的聽診鑒別能力,尚未能夠實現(xiàn)智能診斷判讀。當患者在家受到肺部疾病影響時,無法自行發(fā)現(xiàn)肺部異常而耽誤治療;當處于呼吸道傳染病救治過程中,入耳式的聽診器容易受到污染而造成院內感染。盡管聽診音包含了豐富的健康狀態(tài)信息,由于缺乏標準化的采集方法、分類標準和分析工具,使得聽診音的客觀分析和應用在實踐中受到了限制。本研究通過采用統(tǒng)一的聽診音采集設備和流程進行肺部聽診音數(shù)據(jù)采集、整理、數(shù)據(jù)庫設計,使用軟件MatlabR2017a進行數(shù)據(jù)管理和分析,建立了健康群體和肺部疾病患者群體的肺部聽診音數(shù)據(jù)庫,制訂一套標準的聽診音分類、標注規(guī)范、音頻特征信號參數(shù),構建一個用于存儲、管理和分析肺部聽診音數(shù)據(jù)的系統(tǒng),為肺部疾病的篩查、監(jiān)測以及醫(yī)學人工智能應用轉化等相關研究提供重要的數(shù)據(jù)支持。本研究為肺部聽診音音頻數(shù)據(jù)庫建庫積累了經(jīng)驗,為音頻類數(shù)據(jù)庫管理和分析提供有益的參考和借鑒,為支持后續(xù)醫(yī)學人工智能輔助聽診應用于肺部疾病篩查與監(jiān)測奠定基礎,具有重要的醫(yī)學價值和實際應用意義。

    【關鍵詞】 肺疾??;肺部聽診音;音頻數(shù)據(jù)庫;支持向量機;特征識別;數(shù)據(jù)分析

    【中圖分類號】 R 563 【文獻標識碼】 A DOI:10.12114/j.issn.1007-9572.2023.0863

    Study of Techniques and Methods for Building a Database of Lung Auscultation Sounds

    ZHANG Dongying1,2,YE Peitao3,LI Qiasheng2,JIAN Wenhua2,LIANG Zhenyu2,ZHENG Jinping2*

    1.Faculty of Medicine,Macau University of Science and Technology,Macau 999078,China

    2.The First Affiliated Hospital of Guangzhou Medical University/Guangzhou Institute of Respiratory Health/National Clinical Research Center for Respiratory Diseases,Guangzhou 510120,China

    3.Guangdong Second Provincial General Hospital,Guangzhou 510310,China

    *Corresponding author:ZHENG Jinping,Chief physician/Professor/Doctoral supervisor;E-mail:18928868238@163.com

    【Abstract】 Currently,the results of lung sound auscultation with either physical or electronic stethoscopes still rely mainly on the doctor's professional auscultation identification ability,which has not yet been able to realise intelligent diagnosis and interpretation. When patients are affected by lung diseases at home,they are unable to detect lung abnormalities on their own and delay treatment;when they are in the process of rescue and treatment of respiratory infectious diseases,in-ear stethoscopes are easily contaminated and cause nosocomial infections. Although stethoscopic sounds contain a wealth of information about health status,the lack of standardised collection methods,classification criteria and analysis tools has limited the objective analysis and application of stethoscopic sounds in practice. In this study,the data collection,arrangement and database design of the lung auscultation sound were carried out by using the unified auscultation sound collection equipment and process. The study used the software MetlabR2017a for data management and analysis to create a database of lung auscultation sounds in a healthy group and a group of patients with lung disease. A database of lung auscultation sounds was established for healthy groups and groups of patients with lung diseases. A standard set of classification of auscultatory tones,labelling specifications,audio characteristic signal parameters were developed. Building a system for storing,managing and analysing lung auscultation sound data to provide important data support for research related to the screening and monitoring of lung diseases and the translation of medical artificial intelligence applications. The study accumulated the experience of building an audio database of lung auscultation sounds,provided a useful reference for the management and analysis of the audio database,and laied the foundation for supporting the subsequent application of medical artificial intelligence-assisted auscultation in the screening and monitoring of lung diseases,which was of great medical value and practical application.

    【Key words】 Lung diseases;Lung auscultation sound;Audio database;Support vector machine;Feature recognition;Data analysis

    隨著醫(yī)學科學技術的發(fā)展,肺部疾病的診斷越來越依賴于醫(yī)學數(shù)據(jù),肺部聽診音作為肺部疾病診斷的重要依據(jù)之一,其標準音頻數(shù)據(jù)庫建立和數(shù)據(jù)分析顯得尤為重要。然而,建立一個高質量的肺部聽診音數(shù)據(jù)庫并非易事。醫(yī)務工作者需要考慮如何確保數(shù)據(jù)的準確性和可靠性,如何有效地組織和存儲數(shù)據(jù),如何利用人工智能(AI)和機器學習技術對數(shù)據(jù)進行處理和分析等問題[1]。近年來,AI和機器學習技術為解決這些問題提供了新的可能。這些技術能夠提高數(shù)據(jù)采集、處理和分析的效率和準確性,使人們從海量的醫(yī)療數(shù)據(jù)中提取出更有價值的信息[2]。因此,利用這些技術構建肺部聽診音數(shù)據(jù)庫系統(tǒng),成為當前肺部聽診音數(shù)據(jù)庫建庫技術及方法研究的重要方向。

    肺部聽診是醫(yī)生在診斷和監(jiān)測呼吸系統(tǒng)疾病時常用的重要手段。肺部音頻信號可以提供許多有價值的信息,如呼吸音、雜音等,可以輔助醫(yī)生診斷和判斷疾病的類型和嚴重程度[3]。然而,肺部聽診音的解讀和分析常需要醫(yī)生具備豐富的經(jīng)驗和專業(yè)知識。傳統(tǒng)的肺部聽診音主要依靠臨床實踐和經(jīng)驗積累,但這種方式的局限性在于其主觀性較強,存在個體差異,對醫(yī)生培訓和病情判斷的一致性有一定影響[4-5]。為了解決這些問題,國內外研究者探索利用機器學習和信號處理技術來對肺部聽診音數(shù)據(jù)進行分析和建模。通過機器學習和神經(jīng)網(wǎng)絡模型,AI可以對大量經(jīng)過聲音標注的呼吸音數(shù)據(jù)進行快速學習,從而實現(xiàn)對呼吸音的自動分析[6]。呼吸音的自動識別和分類可以分為兩步:特征提取和識別分類。然而,在呼吸音領域的AI發(fā)展過程中,面臨一個主要問題,即缺乏可用于研究算法和比較不同算法的公共高質量大型呼吸音數(shù)據(jù)庫,這也是AI在呼吸音領域發(fā)展的主要障礙[7]。建立肺部聽診音數(shù)據(jù)庫可以為自動化診斷系統(tǒng)、疾病監(jiān)測和健康評估等提供基礎。隨著智能醫(yī)療和遠程醫(yī)療的發(fā)展,建立一個公共的、大規(guī)模的肺部聽診音數(shù)據(jù)庫具有重要意義。這樣的數(shù)據(jù)庫可以促進醫(yī)療資源共享和協(xié)作研究,提高肺部疾病的診斷和治療水平。因此,對于肺部聽診音數(shù)據(jù)建庫技術及方法的研究可以幫助改善傳統(tǒng)肺部聽診方法的不足且推動肺部疾病的早期診斷和治療。

    1 肺部聽診音的采集

    1.1 音頻數(shù)據(jù)采集設備

    本研究采用美國3M公司的3MTM Littmann?3200聽診器型電子聽診器獲取肺部聽診音信息數(shù)據(jù),該聽診器有以下特點:通過數(shù)字信號傳導而非聲學振動來采集呼吸音數(shù)據(jù);可以收集更廣的音頻范圍,并進行數(shù)字信號轉換;可同步生成視覺聲音圖譜,即在聽到呼吸音的同時可看到呼吸音的波形圖譜。上述特點有助于使用者調整聽診技術,獲得更準確的聽診結果[8]。該類型聽診器與普通聽診探頭相比,可以聽到更低頻的聲音,音域更廣、響度更大,且具有良好的環(huán)境降噪效果。譬如其采集頻率范圍為20~2 000 Hz,采樣頻率為4 kHz,可將聲音放大24倍,并能降低聽診探頭端的雜音,提高聲學傳導質量。此外,其還具備錄音、存儲和數(shù)據(jù)傳輸功能,可以連接電腦,并可使用匹配的音頻軟件進行分析,是建立標準肺部聽診音的理想工具。

    1.2 數(shù)據(jù)獲取

    1.2.1 研究對象準備:本研究對象均來自廣州醫(yī)科大學附屬第一醫(yī)院(患者來源于呼吸科住院部,健康人群來源于醫(yī)院員工),入組前使其充分知情肺部聽診音建庫研究,告知采集流程和需其配合的環(huán)節(jié),并請其簽署知情同意書面材料。依據(jù)被采集人的年齡、性別、基礎病、被采集時的身體狀態(tài)等記錄詳細信息并做好分組。本項目是一項開放式的研究探索,參考既往相關研究,招募總計350名參與者作為研究樣本,這些參與者被劃分為7組,每組各50例,以便于進行詳盡的比較和分析。該設計標準可確保研究的準確性和有效性,同時也保護參與者的權益。本研究經(jīng)廣州醫(yī)科大學附屬第一醫(yī)院倫理審查通過(醫(yī)科倫審2017第82號)。

    1.2.1.1 納入標準:(1)志愿參與。(2)年齡≥18歲。(3)健康人群:①既往無肺部慢性疾病、無粉塵長期暴露史且無吸煙習慣,不受器質性心血管疾病的困擾者;②近3周內無急性上呼吸道感染及咳嗽、咳痰等呼吸道癥狀;③所有研究對象在過去6個月內進行的健康體檢中,胸部X線片未發(fā)現(xiàn)任何肺部病變,并且心電圖表現(xiàn)正常。COPD患者:①有呼吸困難、慢性咳嗽或咳痰,有危險因素暴露史;②肺功能檢查,吸入支氣管擴張劑后第1秒用力呼氣量(FEV1)/用力肺活量(FVC)lt;0.7;③無大量胸腔積液或氣胸并發(fā)癥;④臨床上確診為COPD;符合第1~3項或符合第3~4項。支氣管哮喘患者:①經(jīng)常出現(xiàn)反復喘息、氣急、胸悶或咳嗽,多與接觸變應原、冷空氣、物理、化學刺激、病毒性上呼吸道感染、運動等有關;②發(fā)作時在雙肺可聞及散在或彌漫性、以呼氣相為主的哮鳴音,呼氣相延長;③上述癥狀通常在使用平喘藥物治療后緩解或自行緩解;④除外其他疾病引起的喘息、氣急、胸悶或咳嗽;⑤對癥狀不典型者(如無明顯喘息或體征),應最少具備以下3項中至少1項陽性:(a)支氣管激發(fā)試驗(或運動激發(fā)試驗)陽性,(b)支氣管舒張試驗陽性,(c)晝夜呼氣峰值流量(PEF)變異率≥20%;⑥臨床上確診為支氣管哮喘;符合第1~4項或符合第4~5項或符合第6項者。支氣管擴張患者:①幼年時期有誘發(fā)支氣管擴張的呼吸道感染史,如麻疹、百日咳或流感后肺炎病史,或肺結核病史等;②出現(xiàn)長期慢性咳嗽、咳膿痰或反復咯血癥狀;③高分辨率計算機斷層掃描(HRCT)顯示支氣管擴張的異常影像學改變;④臨床上確診為支氣管擴張;符合第1~3項或符合第4項者。間質性肺疾病患者:①在雙側下肺區(qū)會有干啰音或濕啰音,尤其是在吸氣末期更加明顯;②胸部X線片顯示彌漫性的陰影或不常見的結節(jié)樣間質性陰影,以及支氣管壁增厚和血管周圍間隙的突出,小的規(guī)則或不規(guī)則陰影以及小的外周性環(huán)狀陰影;③肺功能檢測揭示了限制性通氣功能障礙和氣體交換功能減退;④通過對支氣管肺泡灌洗液(BALF)的細胞學、病原學、生化和炎癥遞質等檢測,可以依據(jù)BALF中炎癥免疫效應細胞的比例將間質性肺疾病分為淋巴細胞增多型和中性粒細胞增多型;⑤通過支氣管肺活檢或外科肺活檢獲取肺組織進行病理學檢查確診為間質性肺疾病;⑥臨床上確診為間質性肺疾病。肺炎患者:①咳嗽、咳痰,或原有呼吸道疾病癥狀加重,并出現(xiàn)膿性痰,伴或不伴胸痛;②發(fā)熱,體溫≥38 ℃;③在肺部實變病變或濕啰音中可以觀察到感染的跡象;④白細胞計數(shù)(WBC)gt;10×109/L或lt;4×109/L,伴隨或不伴隨核左移;⑤胸部X線檢查可能會顯示斑片狀浸潤陰影或間質性變化,時伴隨有胸腔積液;⑥臨床上確診為肺炎;符合第1~4項中任意一項和第5項。肺部腫瘤患者:①胸部CT提示肺部腫瘤,或提示肺部占位且考慮腫瘤可能性大;②肺部腫瘤直徑gt;2 cm;③臨床上考慮肺部腫瘤可能;④未合并大量胸腔積液和氣胸;⑤未行手術切除;⑥臨床上已確診為肺部腫瘤;符合第1~5項或符合第2、4、5、6項者。肺部腫瘤患者:①胸部CT提示肺部腫瘤,或提示肺部占位且考慮腫瘤可能性大;②肺部腫瘤直徑gt;2 cm;③臨床上考慮肺部腫瘤可能;④未合并大量胸腔積液和氣胸;⑤未行手術切除;⑥臨床上已確診為肺部腫瘤;符合第1~5項或符合第2、4、5、6項者。

    1.2.1.2 排除標準:(1)確診患有癲癇并且正在接受藥物治療以控制病情;(2)妊娠;(3)肺部切除手術的患者;(4)處于昏迷狀態(tài)或有意識障礙的患者;(5)無法配合研究流程的個體;(6)在重癥監(jiān)護病房(ICU)接受治療或病情處于危急狀態(tài)的患者;(7)未簽署知情同意書;(8)同時滿足以上2個或2個以上的排除標準的患者。

    1.2.2 錄音環(huán)境準備:要求在獨立安靜、溫度適宜的空間進行肺部聽診錄音,錄音時環(huán)境噪聲不高于40~50 dB。

    1.2.3 資料收集:收集研究對象的病史詳細信息及其相關檢查資料(包括年齡、身高、體質量、主要癥狀、吸煙史、呼吸系統(tǒng)相關主要診斷、胸部CT結果、肺功能檢查報告數(shù)據(jù)等)。

    1.2.4 錄音方法:囑研究對象放松心情,平靜呼吸數(shù)分鐘。選用3MTM Littmann?3200型電子聽診器,按人民衛(wèi)生出版社第9版《診斷學》肺部聽診要求對喉部及兩側上、中、下肺部、左右肩胛骨下共9個部位進行聽診(圖1),錄取所有入組研究對象的肺部聽診音,同時填寫資料表格。每個點至少采集15 s的呼吸音,研究對象保持安靜。聽診器探頭緊貼研究對象皮膚,開始采集后,囑研究對象深呼吸,確保每次采集的呼吸音包含2~3個完整的呼吸周期[9]。

    1.2.5 質量控制:由2名呼吸內科主治醫(yī)師分別聽取所采取的肺音進行錄音鑒別與診斷。如存在分歧,再組織第3名呼吸內科主治醫(yī)師聽取鑒別,消除歧義后錄用數(shù)據(jù)。

    2 數(shù)據(jù)傳輸與存儲

    2.1 數(shù)據(jù)傳輸

    3MTM Littmann?3200型電子聽診器內存有限,可存儲12段長度為30 s的音頻,因此錄制完每個研究對象需將肺音音頻通過藍牙和聽診器適配的軟件StatAssist傳輸?shù)诫娔X端,需注明每段音頻來源于具體研究對象所對應的肺部部位。

    2.2 數(shù)據(jù)存儲

    傳輸?shù)诫娔X端的肺音音頻以zsa格式保存,文件名會自動產(chǎn)生一串編碼,為了與相應研究對象對應,可在編碼后加上研究對象姓名以對應到具體的人,該段zsa格式的音頻即為原始數(shù)據(jù)。

    2.3 音頻導出

    將zsa格式音頻導出為wav格式保存,以便于后續(xù)音頻標注和分析。將同一研究對象的肺音音頻導入同一文件夾,文件夾以原文件名命名,按每段音頻所對應部位重命名,以便于每條肺音可以準確對應到某個研究對象的具體身體部位。

    2.4 數(shù)據(jù)脫敏

    包含研究對象姓名的肺音音頻數(shù)據(jù)屬于敏感信息,可能因設備遺失或電腦病毒而導致信息泄露,故需對數(shù)據(jù)進行脫敏。將文件名全部隨機重命名,并將重命名后的數(shù)據(jù)與研究對象姓名、基本信息等對應起來。

    2.5 云端數(shù)據(jù)庫的建立

    存儲在電腦端的數(shù)據(jù)有可能因電腦故障而丟失,故將脫敏數(shù)據(jù)傳入云端是很有必要。設立一個只供本地網(wǎng)絡訪問的云端數(shù)據(jù)庫,對數(shù)據(jù)的存儲及安全均有良好的保障。將未標注、標注后及審核后各階段數(shù)據(jù)分別保存在云端數(shù)據(jù)庫,方便統(tǒng)計數(shù)據(jù)納排情況及數(shù)據(jù)溯源。

    3 音頻數(shù)據(jù)的預處理

    3.1 肺音降噪

    肺音信號的預處理首先是去除肺音中的噪聲,目的是提高呼吸音信號的清晰度和質量,準確地捕捉和識別肺部聲音信號。電子聽診器雖然有一定的降噪能力,但在肺音的采集過程中會不可避免地引入一些噪聲,包括來源于外部環(huán)境雜聲、聽診器電路內部產(chǎn)生的電流干擾以及探頭與皮膚接觸面的摩擦聲和心音等人體內部噪聲等[10]。原始肺音信號的噪聲來源很多,如果不能有效地濾除噪聲,將對肺音的識別過程產(chǎn)生較大影響。

    3.1.1 肺音信號的高通濾波去噪:肺音信號的高通濾波去噪的目的是去除低頻噪聲,保留高頻部分,從而提取出肺音信號的有用信息。

    肺音信號中的呼吸音和心臟音一般集中在較高的頻率范圍內,而噪聲通常是低頻信號。通過應用高通濾波器,可以抑制低頻噪聲的干擾,使得肺音信號更加清晰。高通濾波的關鍵是選擇合適的截止頻率和濾波器參數(shù),以充分去除低頻噪聲,同時保留肺音信號中的有用信息。

    3.1.2 肺音信號的低通濾波去噪:根據(jù)肺音信號的特點,設計一個低通濾波器,可以將肺音信號中高頻部分去除,保留低頻部分。進行濾波操作,將設計好的低通濾波器應用到肺音信號上,通過濾波操作可以將高頻噪聲去除,同時保留肺音信號中的低頻成分。進行信號重構,經(jīng)過濾波操作后,得到的信號可能會有一定的變形,需要進行信號重構處理,使其恢復到原始形態(tài)。

    3.2 人工標注

    3.2.1 音頻標注軟件:Audacity是一款公開、免費的音頻處理軟件,用于錄音和編輯音頻,是自由、開放源代碼的軟件分析。該軟件具有的可視化工具可用于分析音頻或其他信號數(shù)據(jù),支持Python、Perl 或任何支持命名管道的語言編寫腳本,兼容性強,可在Mac OS X、Microsoft Windows、GNU/Linux和其他操作系統(tǒng)上運作。

    3.2.2 肺部聽診音標注者的選擇:肺部聽診音標注的準確性在很大程度上影響模型的識別效率,呼吸內科醫(yī)生對肺部聽診音的識別比較有經(jīng)驗,也有利于提高音頻數(shù)據(jù)標注的準確度,故建議呼吸內科醫(yī)生作為標注者。由于采錄設備錄制的聲音與臨床實際聽到的聲音有些許差別,所以每位標注者需經(jīng)過專業(yè)培訓后方能正式參與肺音標注。

    3.2.3 標注流程:將待標注音頻導入Audacity軟件,將默認的波形圖轉換為頻譜圖模式,頻譜圖可直觀顯示一個呼吸周期及該呼吸周期的吸氣相和呼氣相,有利于快速選區(qū)并標注,選定區(qū)域后輸入標簽。為了使標注更加精確,將整段音頻中每個呼吸周期先框選出來,再對吸氣相和呼氣相進行區(qū)分后作定性(標注為正常、濕啰音、哮鳴音、喘鳴音或鼾音),有助于呼吸音的特征提取。不完整的呼吸片段可能會降低特征提取的效果,故只標注較為完整的呼吸周期。

    標注的內容包括:呼吸周期的呼氣段和吸氣段,正常、濕啰音和干啰音發(fā)生的位置(精確到ms),干啰音細分的類型。還包括對呼吸音音頻質量等級標注,音頻質量等級按照音頻質量的好壞,依次分級為Ⅰ、Ⅱ、Ⅲ級,其中Ⅰ級最優(yōu)。評價標注包括音頻的信噪比、噪聲時長占比及人工主觀評價。標注完成后,帶標簽的音頻將以Audacity Project File(aup3)格式保存。標注時所有標注者使用同種品牌同種型號的收聽設備(鐵三角ATH-M20X頭戴式耳機),使所有標注者盡可能獲得較為一致的聽覺體驗。

    3.2.4 標簽審核:對于標注過程中有爭議的部分,需進行公開討論,最終獲得一致意見的,予以采納;若最終無法達成一致意見,則剔除該條數(shù)據(jù)。全部音頻標注完成后,由呼吸與危重癥學科高年資醫(yī)生對所有音頻進行審核,修改錯誤或不規(guī)范標簽,并剔除低質量數(shù)據(jù)(噪聲太多影響判斷或音量太小難以分辨的音頻)。經(jīng)過審核的標簽即可用于訓練AI模型。

    3.3 肺音識別機器學習算法及其模型建立過程

    3.3.1 濕啰音和哮鳴音識別及其持續(xù)時間:濕啰音和哮鳴音是2種常見的呼吸音,可以通過聽診器或其他聽診設備進行識別。濕啰音是由于呼吸道黏膜分泌物增多導致的濕潤聲音,其通常類似于呼吸道中的液體在被吸入或呼出時產(chǎn)生的咕嚕聲。濕啰音可以是細小而清澈的,也可以是較為粗糙和響亮的,其可能由于感染、充血、分泌物過多或阻塞引起。濕啰音的持續(xù)時間可以因患者情況而異,從持續(xù)數(shù)秒到幾分鐘不等[11]。哮鳴音是由于氣流通過狹窄的呼吸道導致的高音頻咳嗽聲或呼哧聲。這種聲音通常類似于從狹窄的空間經(jīng)過的高速氣流所產(chǎn)生的噪聲。哮鳴音可以是連續(xù)的或間斷的,并且強度和音調可能會變化;其可能由于肺部狹窄、支氣管痙攣、梗阻或病變引起。哮鳴音的持續(xù)時間也可以因患者情況而異,從數(shù)秒鐘到數(shù)分鐘不等[12]。

    濕啰音是偶發(fā)性的非連續(xù)信號,基于音頻數(shù)據(jù)分析,筆者發(fā)現(xiàn)其持續(xù)時間一般lt;20 ms。本研究團隊使用軟件Matlab,呈現(xiàn)出16個濕啰音頻譜圖,其中紅色的部分為以濕啰音峰值為中心的20 ms信號,詳見圖2。

    哮鳴音一般布滿整個呼氣相或吸氣相,其持續(xù)時間一般gt;250 ms。本研究采用紅色標出以哮鳴音片段峰值為中心的500 ms信號,詳見圖3。經(jīng)過大量數(shù)據(jù)的測試,本研究團隊發(fā)現(xiàn)500 ms在絕大多數(shù)情況下能夠覆蓋呼吸相或吸氣相的80%以上。

    3.3.2 濕啰音的波峰和哮鳴音的波峰定量分析:通過對超過1 500個濕啰音進行統(tǒng)計,發(fā)現(xiàn)濕啰音的波峰幅值是濕啰音平均幅值的2倍以上。使用軟件Matlab,呈現(xiàn)出了16個濕啰音頻譜圖,其中紅色的線是濕啰音平均幅值的2倍,詳見圖4。

    通過對超過100個哮鳴音測試,發(fā)現(xiàn)用200 ms長的漢明窗對哮鳴音能量值進行低通濾波,絕大部分低通濾波后的哮鳴音片段有160 ms的峰值部分會大于濾波后信號的平均值(圖5的紅線)。如圖5所示,本研究用該方法找到了該哮鳴音的所有哮鳴音片段峰值,并用綠色線標出。

    3.3.3 利用支持向量機(support vector machine,SVM)法進行濕啰音和哮鳴音的識別:SVM是一種常用的監(jiān)督學習算法,主要用于分類和回歸分析。其是一種二分類模型,通過最大化分類邊界與數(shù)據(jù)之間的間隔來對新的數(shù)據(jù)進行預測。SVM是基于統(tǒng)計學習理論,通過在高維空間中找到最優(yōu)分類超平面來對樣本進行分類,還可以用于回歸分析、異常檢測和時間序列預測等任務,是一種廣泛應用于機器學習和數(shù)據(jù)挖掘領域的算法,具有較高的準確性和魯棒性[13]。

    利用SVM進行肺部聽診音識別是一種基于機器學習和信號處理的技術,其基本原理是通過對肺部聽診音信號進行特征提取和分類,從而實現(xiàn)肺部聽診音的識別和分類。通過訓練SVM模型,可以根據(jù)已有的肺部聽診音數(shù)據(jù)和對應的標簽(如正常肺音、異常肺音)來建立分類模型[14]。這個模型可用于對肺部聽診音數(shù)據(jù)進行預測和分類,從而輔助臨床診斷。SVM在肺部聽診音識別中的優(yōu)勢在于,其可以處理高維數(shù)據(jù)和復雜的非線性關系,通過核函數(shù)將原始數(shù)據(jù)映射到高維特征空間中,從而避免了手動設計特征和選擇特征空間的復雜性。同時,SVM具有較高的準確性和魯棒性,能夠處理各種不同類型和質量的肺部聽診音數(shù)據(jù)。

    在本研究中,筆者團隊嘗試利用SVM法對肺部聽診音進行濕啰音和哮鳴音的識別,主要包括音頻數(shù)據(jù)庫的建立、信號特征提取、向量機訓練及向量機預測4個部分。與絕大多數(shù)研究不同的是,本研究的訓練數(shù)據(jù)集包含的并不是整個肺音采錄文件,而是經(jīng)過預處理的音頻數(shù)據(jù)片段,包括正常音頻片段和異常音頻片段。其中,濕啰音片段以濕啰音片段峰值為中心的20 ms信號組成,哮鳴音片段以哮鳴音片段峰值為中心的500 ms信號組成。這樣做的好處主要有兩個,其一,以峰值為中心的信號片段信號強,信噪比高,從而可以比較好地降低背景噪聲造成的干擾,提高算法的準確度;其二,以峰值為中心的信號片段數(shù)據(jù)量較小,從而明顯減少信號處理時間,提高算法的實時性、實用性。

    在采集肺部聽診音的過程中,本研究發(fā)現(xiàn)聽診器探頭薄膜隨著健康受試者吸氣胸腔擴張被動錯位或探頭在轉移過程中均會產(chǎn)生偽濕啰音,且現(xiàn)階段無法避免[15]。為了能夠提高識別濕啰音的準確度,筆者也對偽濕啰音進行人工標注和識別且建立偽濕啰音的音頻數(shù)據(jù)庫,為后續(xù)能夠實現(xiàn)對偽濕啰音的自動識別作準備。本研究提取的每個濕啰音片段時長為20 ms,通過算法及半人工參與的方式建立真性和偽濕啰音片段庫,詳見圖6,顯示的是真、偽濕啰音片段建庫過程。

    為了進一步處理和分析肺音信號,本研究運用小波變換方法,對濕啰音與偽濕啰音進行處理,并提取其頻域特征。這些特征值是從信號的基本波形和頻譜中提取出來的。將計算出的時域和頻域特征值作為輸入變量,利用經(jīng)過訓練的SVM模型進行分類。通過信號的小波變換分解、特征值的計算,以及SVM分類模型的應用能夠以較高的準確率來識別和區(qū)分正常與異常的呼吸音。具體的算法識別過程與結果如下。

    (1)本研究結合自動算法和人工審核技術,深入分析肺音數(shù)據(jù)。研究建立了兩個數(shù)據(jù)庫,一個收集真實濕啰音,另一個收集偽濕啰音。由于真實濕啰音通常不超過20 ms,研究確保所有聲音樣本長度為20 ms,以精確捕獲這些短暫信號。通過分析這些樣本,提取了關鍵的時域和頻域特征,這對分類濕啰音很重要,詳見圖7。本研究還使用SVM算法對特征進行訓練和優(yōu)化,成功區(qū)分真、假濕啰音。

    (2)研究首先對肺音信號進行預處理,然后通過帶通濾波器處理,強調100~2 000 Hz頻率范圍內的關鍵信號,該范圍被認為含有診斷相關的重要信息。使用20 ms的時間窗口逐步掃描整個信號,以連續(xù)搜索并捕捉瞬時濕啰音。每發(fā)現(xiàn)一個潛在的濕啰音,就從該段信號提取時域和頻域特征,并將這些特征用于SVM模型,以區(qū)分真實和偽濕啰音。

    通過算法及半人工參與,從肺音中提取兩個聲音庫,其中一個聲音庫由濕啰音組成,另外一個由偽濕啰音組成。由于濕啰音基本lt;20 ms,因此本研究提取的每個聲音檔為20 ms。

    提取真濕啰音和偽濕啰音的時域特征值及頻域特診值,然后利用向量機進行訓練,得出向量機參數(shù)。利用該向量機參數(shù)可以很好地分辨真濕啰音和偽濕啰音。先進行信號預處理,重采樣和識別弱信號,然后對呼吸音信號進行帶通濾波(肺音信號主要在100~2 000 Hz),然后計算智能閾值(每段呼吸音信號有自己的閾值,根據(jù)該閾值找出潛在濕啰音),然后用20 ms的窗沿著肺音信號移動,找出所有潛在濕啰音。提取潛在濕啰音的時域特征和頻域特征,并用向量機進一步判斷該潛在濕啰音是真濕啰音還是偽濕啰音,最后根據(jù)該呼吸音平均1 s所含濕啰音的數(shù)量來判斷該呼吸音是否為濕啰音或哮鳴音,圖中的藍色圓圈標示的是偽濕啰音,紅色圓圈標示的是真濕啰音,詳見圖8。

    由于哮鳴音會布滿整個呼氣相或吸氣相,所以本研究提取的每個哮鳴音片段為500 ms,絕大部分情況下,500 ms能夠覆蓋整個呼氣相或吸氣相80%以上。通過算法及半人工參與的方式建立哮鳴音片段庫,如圖9所示。應用訓練向量機識別濕啰音和哮鳴音的過程,詳見圖10~11。應用向量機識別肺音的完整流程,詳見圖12。

    3.3.4 利用小波分解法對肺部聽診音進行特征提?。盒〔ǚ纸馐且环N信號分析方法,用于將信號分解為不同尺度和頻率的成分,以便更好地分析信號的特點和特征。小波分解通過使用小波函數(shù)進行伸縮和平移操作,將信號分解為多個頻率成分,譬如濕啰音分解后的信號頻率分布,詳見表1。每個成分具有不同的尺度和形狀。這種分解方法可以在一定程度上消除信號中的噪聲,并且可以對信號進行壓縮,減小數(shù)據(jù)量。小波分解在信號處理、圖像處理、醫(yī)學影像等領域有著廣泛的應用[16]。

    肺部聽診音是一種復雜的生理信號,包含了多種頻率和特征的成分。小波分解應用于肺部聽診音特征提取的原理主要是基于信號的頻域特性。小波分解通過使用小波基函數(shù)對肺部聽診音進行分解,將信號分解為不同尺度和頻率的成分,從而更好地分析信號的特點和特征。這些成分可以反映肺部聽診音在不同時間和空間上的變化特征,從而為特征提取和分析提供更豐富的信息[17]。通過小波分解,可以提取肺部聽診音中的不同頻率成分,例如高頻成分可能代表呼吸急促、咳嗽等癥狀,低頻成分可能代表呼吸節(jié)律等特征。這些特征可以用于構建分類模型,用于識別不同的肺部疾病或癥狀。

    本研究將輸入聲音片段分解成和D1、D2、D3、D4、D5、D6、D7、A8共8段不同頻率的信號,小波分解肺音的過程,詳見圖13。濕啰音音頻的分解結果,詳見圖14。

    對于分解出的8段信號,本研究提取如下特征(共提取32個特征值):(1)提取每段信號絕對值的平均值,共8個平均值;(2)提取每段信號絕對值的方差,共8個方差;(3)提取每段信號絕對值的最大值,共8個最大值;(4)提取8個平均值相鄰值得比值。分別為mean(abs(s[x]))/mean(abs(D1),mean(abs(D1))/mean(abs(D2)),mean(abs(D2))/mean(abs(D3)),mean(abs(D3))/mean(abs(D4)),mean(abs(D4))/mean(abs(D5)),mean(abs(D5))/mean(abs(D6)),mean(abs(D6))/mean(abs(D7)),mean(abs(D7))/mean(abs(A7))。

    特征1~8(8個平均值)和特征17~24(8個最大值)代表了濕啰音在不同頻率的能量分布,特征9~16(8個方差)和25~32(相鄰值得比值)代表了濕啰音在不同頻率的能量變化情況。

    4 應用與結論

    應用上述采集肺部聽診音的操作流程與建立肺部聽診音音頻數(shù)據(jù)庫的技術方法,本研究采集了392例呼吸慢病患者,合格入庫3 362條音頻數(shù)據(jù),另有56例健康人的肺部聽診音500條音頻數(shù)據(jù),以及采用生物醫(yī)學電子與生物信息學研究生院公布的HF_Lung_V1、“國際”生物醫(yī)學與健康信息學會議公布的呼吸音數(shù)據(jù)庫及上海交通大學公布的SPRSound數(shù)據(jù)庫共9 469條音頻數(shù)據(jù)。應用軟件MetLab R2017a進行數(shù)據(jù)管理和分析,對音頻數(shù)據(jù)進行預處理、數(shù)據(jù)標注、向量機識別、特征提取等,基于呼吸音分析算法軟件,可以自動導入呼吸音并頻譜化顯示,并通過SVM算法建立的正常和異常呼吸音模型庫,自動識別出正常肺音,干、濕啰音,并在呼吸音頻譜圖上精確標注所在位置。應用程序界面識別肺部聽診音示例,詳見圖15。

    呼吸疾病的早發(fā)現(xiàn)、早診斷、早干預,有利于獲得呼吸慢病患者有效管理的關鍵節(jié)點。肺部聽診音數(shù)據(jù)庫的建立和相關方法研究為日后實現(xiàn)對肺部聽診音的智能判讀和居家應用健康管理產(chǎn)品的研發(fā)奠定基礎,甚至可以協(xié)助呼吸道傳染臨床診治,也為后期研制國產(chǎn)質優(yōu)價廉的智能判讀聽診器提供技術借鑒。

    本研究局限性:目前該數(shù)據(jù)庫中的病理性肺音片段相對較少,這可能是由以下幾個因素造成的,(1)本研究并沒有涵蓋那些病情嚴重的患者,這意味著那些可能會呈現(xiàn)更加明顯肺部異常聲音的患者群體并未被包括在內;(2)那些同時受多種呼吸系統(tǒng)疾病影響的患者也未能被納入本研究,復雜的臨床情況可能會造成多樣化的肺音變化,而這部分信息對于肺音數(shù)據(jù)庫而言是極其寶貴的;(3)醫(yī)生在更新患者病歷時可能會有所延遲,這可能導致在收集肺音的時候,患者的實際健康狀況已經(jīng)有所改善,比如,如果患者的病情在穩(wěn)定期,那么其肺部的異常聲音可能會相對較弱,甚至在檢測時不再明顯。這種情況下,收集到的肺音可能就不能準確反映患者在更加活躍病程階段的肺音特征。

    致謝:上海拓蕭智能科技有限公司協(xié)助算法開發(fā)。

    作者貢獻:張冬瑩提出肺音研究思路,設計研究方案,研究命題的提出、設計,包括肺音聽診對象分組,偽濕啰音的智能識別和目標導向的研究思路(智能判讀聽診器的推廣應用導向),以及項目研究過程中專利申報的思路等負責論文起草,負責最終版本修訂,對論文負責;張冬瑩、李洽勝、簡文華、梁振宇負責研究對象甄別入組,組織研究過程的實施,負責質量控制;張冬瑩、葉培韜負責數(shù)據(jù)收集、采集、清洗和統(tǒng)計學分析、繪制圖表等;鄭勁平為研究項目顧問,指導研究方案及論文修改。

    本文無利益沖突。

    參考文獻

    SRIVASTAVA A,JAIN S,MIRANDA R,et al. Deep learning based respiratory sound analysis for detection of chronic obstructive pulmonary disease[J]. PeerJ Comput Sci,2021,7:e369. DOI:10.7717/peerj-cs.369.

    BARDOU D,ZHANG K,AHMAD S M. Lung sounds classification using convolutional neural networks[J]. Artif Intell Med,2018,88:58-69. DOI:10.1016/j.artmed.2018.04.008.

    GOTTLIEB E R,ALIOTTA J M,TAMMARO D. Comparison of analogue and electronic stethoscopes for pulmonary auscultation by internal medicine residents[J]. Postgrad Med J,2018,

    94(1118):700-703. DOI:10.1136/postgradmedj-2018-136052.

    KWON A M,KANG K. A temporal dependency feature in lower dimension for lung sound signal classification[J]. Sci Rep,2022,12:7889. DOI:10.1038/s41598-022-11726-3.

    CHAMBERLAIN D,KODGULE R,GANELIN D,et al. Application of semi-supervised deep learning to lung sound analysis[J]. Annu Int Conf IEEE Eng Med Biol Soc,2016,2016:804-807. DOI:10.1109/EMBC.2016.7590823.

    葉培韜,張冬瑩,陳文雅,等. 人工智能在呼吸音方面的研究進展[J]. 國際呼吸雜志,2022,42(13):961-966. DOI:10.3760/cma.j.cn131368-20220120-00045.

    AVILES-SOLIS J C,VANBELLE S,HALVORSEN P A,et al. International perception of lung sounds:a comparison of classification across some European borders[J]. BMJ Open Respir Res,2017,4(1):e000250. DOI:10.1136/bmjresp-2017-000250.

    HOFFMANN C,F(xiàn)ALZONE E,VERRET C,et al. Brief report:pulmonary auscultation in the operating room:a prospective randomized blinded trial comparing electronic and conventional stethoscopes[J]. Anesth Analg,2013,117(3):646-648. DOI:10.1213/ANE.0b013e31829ec0a5.

    上海市醫(yī)學會兒科分會呼吸學組,福建省醫(yī)學會兒科分會呼吸學組,江蘇省醫(yī)學會兒科分會呼吸學組,等. 高質量呼吸音數(shù)據(jù)庫建設及管理規(guī)范專家共識[J]. 中華實用兒科臨床雜志,2023,38(7):481-485. DOI:10.3760/cma.j.cn101070-20230120-00054.

    姚小靜,王洪,李燕,等. 肺音信號分析及其識別方法的研究進展[J]. 重慶理工大學學報(自然科學),2013,27(12):95-100.

    ZHANG K X,WANG X F,HAN F F,et al. The detection of crackles based on mathematical morphology in spectrogram analysis[J]. Technol Health Care,2015,23(Suppl 2):S489-494. DOI:10.3233/THC-150986.

    REICHERT S,GASS R,BRANDT C,et al. Analysis of respiratory sounds:state of the art[J]. Clin Med Circ Respirat Pulm Med,2008,2:45-58. DOI:10.4137/ccrpm.s530.

    宋倩,林和平,周琪棟. 支持向量機理論研究[J]. 信息技術,2013,37(9):152-154,159. DOI:10.3969/j.issn.1009-2552.2013.09.042.

    劉國棟. 支持向量機在肺部病理音識別中的應用研究[D]. 天津:南開大學,2015.

    YE P T,LI Q S,JIAN W H,et al. Regularity and mechanism of fake crackle noise in an electronic stethoscope[J]. Front Physiol,2022,13:1079468. DOI:10.3389/fphys.2022.1079468.

    石陸魁,劉文浩,李站茹. 基于LDA和小波分解的肺音特征提取方法[J]. 計算機工程與應用,2017,53(22):116-120,149. DOI:10.3778/j.issn.1002-8331.1605-0330.

    劉文青. 結合EMD和小波變換的肺音信號重構及呼吸周期分割研究[D]. 武漢:江漢大學,2023.

    (收稿日期:2024-01-22;修回日期:2024-04-10)

    (本文編輯:毛亞敏)

    猜你喜歡
    肺疾病支持向量機數(shù)據(jù)分析
    呼吸機輔助呼吸對慢性阻塞性肺疾病急性加重期并發(fā)呼吸衰竭治療效果分析
    慢性阻塞性肺疾病患者緩解期護理過程中實施護理干預對康復的影響
    動態(tài)場景中的視覺目標識別方法分析
    論提高裝備故障預測準確度的方法途徑
    價值工程(2016年32期)2016-12-20 20:36:43
    基于熵技術的公共事業(yè)費最優(yōu)組合預測
    價值工程(2016年29期)2016-11-14 00:13:35
    Excel電子表格在財務日常工作中的應用
    淺析大數(shù)據(jù)時代背景下的市場營銷策略
    新常態(tài)下集團公司內部審計工作研究
    中國市場(2016年36期)2016-10-19 04:31:23
    淺析大數(shù)據(jù)時代對企業(yè)營銷模式的影響
    基于支持向量機的金融數(shù)據(jù)分析研究
    国产免费现黄频在线看| 18禁在线播放成人免费| 精品久久蜜臀av无| 一边亲一边摸免费视频| 亚洲av国产av综合av卡| 亚洲经典国产精华液单| 日本欧美视频一区| 少妇精品久久久久久久| 一边摸一边做爽爽视频免费| 伦理电影大哥的女人| 亚洲无线观看免费| 国产精品99久久久久久久久| 考比视频在线观看| 成人午夜精彩视频在线观看| 91午夜精品亚洲一区二区三区| 国产又色又爽无遮挡免| 国产乱人偷精品视频| 91成人精品电影| 51国产日韩欧美| 国产女主播在线喷水免费视频网站| 美女脱内裤让男人舔精品视频| 国产高清三级在线| 婷婷成人精品国产| 欧美成人午夜免费资源| 午夜影院在线不卡| 午夜视频国产福利| 精品人妻熟女毛片av久久网站| 精品国产一区二区三区久久久樱花| .国产精品久久| 你懂的网址亚洲精品在线观看| 岛国毛片在线播放| 蜜桃久久精品国产亚洲av| 22中文网久久字幕| 大香蕉久久网| 欧美变态另类bdsm刘玥| 国产成人av激情在线播放 | 又黄又爽又刺激的免费视频.| 亚洲精品国产av蜜桃| 国产精品熟女久久久久浪| av网站免费在线观看视频| 五月开心婷婷网| 日日啪夜夜爽| 18禁观看日本| 国产亚洲精品第一综合不卡 | 青春草国产在线视频| 啦啦啦中文免费视频观看日本| 婷婷色综合www| 国产片特级美女逼逼视频| 日日爽夜夜爽网站| 亚洲色图综合在线观看| 亚洲国产精品一区三区| av天堂久久9| 中文字幕av电影在线播放| 插阴视频在线观看视频| 国产精品欧美亚洲77777| 亚洲不卡免费看| 高清毛片免费看| av线在线观看网站| videos熟女内射| 国产亚洲精品第一综合不卡 | 午夜激情福利司机影院| av电影中文网址| 黄片播放在线免费| 综合色丁香网| 免费黄色在线免费观看| 91在线精品国自产拍蜜月| 91在线精品国自产拍蜜月| 国产男人的电影天堂91| 青春草视频在线免费观看| 纯流量卡能插随身wifi吗| 晚上一个人看的免费电影| 建设人人有责人人尽责人人享有的| 人成视频在线观看免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线观看免费高清a一片| 我要看黄色一级片免费的| 啦啦啦啦在线视频资源| 国产日韩欧美视频二区| 日本午夜av视频| freevideosex欧美| 亚洲人成网站在线观看播放| 女人久久www免费人成看片| 我的老师免费观看完整版| av在线观看视频网站免费| 亚洲一区二区三区欧美精品| 国产白丝娇喘喷水9色精品| av线在线观看网站| 成人无遮挡网站| 精品少妇内射三级| a 毛片基地| 欧美xxxx性猛交bbbb| 99久久精品国产国产毛片| 国产成人精品久久久久久| 日韩熟女老妇一区二区性免费视频| 亚洲欧美精品自产自拍| 亚洲色图 男人天堂 中文字幕 | 精品国产一区二区三区久久久樱花| 亚洲精品美女久久av网站| 免费不卡的大黄色大毛片视频在线观看| 少妇人妻 视频| 欧美精品国产亚洲| 精品少妇久久久久久888优播| 亚洲四区av| 久久久久久伊人网av| 中国三级夫妇交换| 亚洲成人手机| 狂野欧美激情性bbbbbb| 久久久国产欧美日韩av| 欧美人与善性xxx| 成年美女黄网站色视频大全免费 | 一级毛片电影观看| 久久97久久精品| 亚洲欧美日韩另类电影网站| 伊人久久国产一区二区| 国产伦理片在线播放av一区| 天堂中文最新版在线下载| av福利片在线| 久久午夜综合久久蜜桃| av国产久精品久网站免费入址| 韩国高清视频一区二区三区| 69精品国产乱码久久久| 亚洲熟女精品中文字幕| 大话2 男鬼变身卡| 久久精品久久久久久久性| 卡戴珊不雅视频在线播放| 久久久午夜欧美精品| 大香蕉久久成人网| 午夜激情久久久久久久| 亚洲精品av麻豆狂野| 久久久久国产精品人妻一区二区| 久久精品国产亚洲av天美| 亚洲久久久国产精品| 亚洲精品aⅴ在线观看| 国产精品秋霞免费鲁丝片| 青青草视频在线视频观看| 午夜激情av网站| 丝袜喷水一区| 亚洲国产日韩一区二区| 欧美亚洲 丝袜 人妻 在线| 精品久久国产蜜桃| 国模一区二区三区四区视频| 久热久热在线精品观看| 久久久精品免费免费高清| 在线天堂最新版资源| 国产一区二区在线观看日韩| 国产高清国产精品国产三级| 久久午夜福利片| 建设人人有责人人尽责人人享有的| 成年女人在线观看亚洲视频| 人体艺术视频欧美日本| 亚洲精品日韩av片在线观看| 9色porny在线观看| 女性被躁到高潮视频| videos熟女内射| 夫妻性生交免费视频一级片| 亚洲精品日韩av片在线观看| 国国产精品蜜臀av免费| 国内精品宾馆在线| 如何舔出高潮| 色5月婷婷丁香| 91午夜精品亚洲一区二区三区| 一本—道久久a久久精品蜜桃钙片| 午夜免费鲁丝| 日韩视频在线欧美| 国产精品欧美亚洲77777| 中文字幕av电影在线播放| 麻豆乱淫一区二区| 亚洲第一av免费看| 人人妻人人澡人人爽人人夜夜| 精品国产乱码久久久久久小说| 18禁在线播放成人免费| av女优亚洲男人天堂| 亚洲av在线观看美女高潮| 精品人妻熟女毛片av久久网站| 日产精品乱码卡一卡2卡三| 少妇的逼好多水| 99久久综合免费| 国产精品蜜桃在线观看| 高清视频免费观看一区二区| av网站免费在线观看视频| 777米奇影视久久| 99久久精品国产国产毛片| 伊人久久国产一区二区| www.av在线官网国产| 久久99蜜桃精品久久| 欧美国产精品一级二级三级| 五月天丁香电影| 精品亚洲成国产av| 看免费成人av毛片| 日本av免费视频播放| 多毛熟女@视频| 插阴视频在线观看视频| 国产欧美日韩一区二区三区在线 | 七月丁香在线播放| 黄色毛片三级朝国网站| 国产乱来视频区| 国产成人午夜福利电影在线观看| 成人毛片60女人毛片免费| 极品人妻少妇av视频| 丰满少妇做爰视频| 国产在线一区二区三区精| 国产成人精品在线电影| 观看av在线不卡| 美女xxoo啪啪120秒动态图| 国产 精品1| 久久久久久久久久久丰满| 色94色欧美一区二区| 亚洲怡红院男人天堂| 18禁在线播放成人免费| 韩国高清视频一区二区三区| 国产av一区二区精品久久| 伊人亚洲综合成人网| 大片电影免费在线观看免费| 午夜老司机福利剧场| 日本黄色片子视频| 免费看不卡的av| 男女啪啪激烈高潮av片| 日本色播在线视频| 国产成人一区二区在线| 国产av码专区亚洲av| 日日啪夜夜爽| 另类亚洲欧美激情| 纵有疾风起免费观看全集完整版| 精品国产一区二区久久| 国产精品一国产av| 99久久精品一区二区三区| 高清午夜精品一区二区三区| 国产一区二区在线观看日韩| 亚洲精品久久午夜乱码| 一区二区三区免费毛片| 欧美日韩亚洲高清精品| 丝袜在线中文字幕| 夜夜骑夜夜射夜夜干| 精品亚洲乱码少妇综合久久| 色94色欧美一区二区| 国产一区亚洲一区在线观看| 秋霞在线观看毛片| 免费大片黄手机在线观看| 久久久久久久久大av| 在现免费观看毛片| 欧美97在线视频| 日本黄大片高清| 一级a做视频免费观看| 水蜜桃什么品种好| 另类亚洲欧美激情| 成年美女黄网站色视频大全免费 | 人人妻人人添人人爽欧美一区卜| 国模一区二区三区四区视频| 亚洲av不卡在线观看| 国产免费现黄频在线看| 我的女老师完整版在线观看| 国产亚洲精品久久久com| 男人添女人高潮全过程视频| 亚洲av欧美aⅴ国产| 欧美少妇被猛烈插入视频| 午夜激情福利司机影院| 麻豆成人av视频| 99热国产这里只有精品6| 亚洲高清免费不卡视频| 精品人妻一区二区三区麻豆| 国产av码专区亚洲av| 午夜激情久久久久久久| 亚洲欧美日韩另类电影网站| 在线精品无人区一区二区三| 一区在线观看完整版| 国产免费视频播放在线视频| 精品亚洲成a人片在线观看| 最近中文字幕2019免费版| 大片电影免费在线观看免费| 一本大道久久a久久精品| 日本av手机在线免费观看| 999精品在线视频| 国产一区二区三区综合在线观看 | 国产又色又爽无遮挡免| 22中文网久久字幕| 天堂8中文在线网| 国产精品蜜桃在线观看| 大香蕉久久网| 久久久久久久久大av| 老司机影院成人| 精品一品国产午夜福利视频| 国产黄色免费在线视频| 国产日韩一区二区三区精品不卡 | 免费看不卡的av| 国产极品天堂在线| 黄片播放在线免费| 亚洲人成网站在线播| 免费观看av网站的网址| 精品久久蜜臀av无| 一级黄片播放器| 亚洲欧美成人精品一区二区| 免费高清在线观看视频在线观看| 啦啦啦中文免费视频观看日本| 国产亚洲精品第一综合不卡 | 日韩中字成人| 一区二区三区乱码不卡18| 国产淫语在线视频| 91久久精品电影网| 欧美日韩视频高清一区二区三区二| 飞空精品影院首页| 99国产综合亚洲精品| 黑人高潮一二区| 久久精品夜色国产| 国产精品99久久久久久久久| 国产一区二区在线观看日韩| 亚洲国产最新在线播放| 国产成人精品一,二区| 中文天堂在线官网| 久久人人爽av亚洲精品天堂| 纯流量卡能插随身wifi吗| 日本免费在线观看一区| 亚洲内射少妇av| 另类精品久久| 亚洲成色77777| 一级a做视频免费观看| 婷婷色av中文字幕| 女人久久www免费人成看片| 日韩中文字幕视频在线看片| 午夜激情福利司机影院| 九草在线视频观看| 99热网站在线观看| 黄色视频在线播放观看不卡| 两个人免费观看高清视频| 午夜av观看不卡| 亚洲高清免费不卡视频| 国产成人精品一,二区| 婷婷色麻豆天堂久久| 亚洲欧美日韩卡通动漫| 丁香六月天网| 美女大奶头黄色视频| 18禁在线播放成人免费| 久久久久久久国产电影| 欧美日韩一区二区视频在线观看视频在线| 男女边吃奶边做爰视频| 嘟嘟电影网在线观看| 亚洲成人av在线免费| 午夜日本视频在线| 天堂8中文在线网| 美女视频免费永久观看网站| 亚洲av二区三区四区| 777米奇影视久久| 在线观看www视频免费| 国产精品一区二区在线观看99| xxx大片免费视频| 欧美丝袜亚洲另类| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品国产色婷婷电影| 久久久精品区二区三区| 性高湖久久久久久久久免费观看| 中国三级夫妇交换| 日韩 亚洲 欧美在线| 夫妻性生交免费视频一级片| 97在线视频观看| 亚洲精品国产av蜜桃| 中文天堂在线官网| 国产成人av激情在线播放 | 免费观看av网站的网址| 欧美性感艳星| 新久久久久国产一级毛片| 国产极品天堂在线| 亚洲国产最新在线播放| 国产成人午夜福利电影在线观看| 草草在线视频免费看| 亚洲欧美中文字幕日韩二区| 黄色配什么色好看| 国产精品三级大全| 精品99又大又爽又粗少妇毛片| 午夜福利网站1000一区二区三区| 人人澡人人妻人| 又黄又爽又刺激的免费视频.| 免费人妻精品一区二区三区视频| 一级毛片我不卡| 精品人妻在线不人妻| 亚洲av成人精品一二三区| 国产精品嫩草影院av在线观看| 伦精品一区二区三区| 成人国产av品久久久| 女人精品久久久久毛片| 亚洲欧美日韩另类电影网站| 男人添女人高潮全过程视频| 男女免费视频国产| 丰满迷人的少妇在线观看| 亚洲av综合色区一区| 美女脱内裤让男人舔精品视频| 亚洲av成人精品一二三区| 丰满少妇做爰视频| 久久毛片免费看一区二区三区| 22中文网久久字幕| 久久人人爽av亚洲精品天堂| 久久久a久久爽久久v久久| 2022亚洲国产成人精品| 日日爽夜夜爽网站| 亚洲丝袜综合中文字幕| 美女主播在线视频| 18+在线观看网站| 日韩大片免费观看网站| 久久精品国产亚洲网站| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 九九爱精品视频在线观看| 色5月婷婷丁香| 亚洲精品乱久久久久久| 欧美3d第一页| 日本免费在线观看一区| 人人澡人人妻人| 免费少妇av软件| 免费观看a级毛片全部| 美女国产高潮福利片在线看| 女性生殖器流出的白浆| 丰满乱子伦码专区| 国产日韩欧美亚洲二区| av在线观看视频网站免费| 国产在视频线精品| 在线播放无遮挡| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品456在线播放app| 国产女主播在线喷水免费视频网站| 啦啦啦视频在线资源免费观看| 日韩欧美一区视频在线观看| 亚洲国产av影院在线观看| 人妻系列 视频| 免费观看av网站的网址| 国产av国产精品国产| 蜜桃在线观看..| 大片电影免费在线观看免费| 亚洲国产精品国产精品| 成人黄色视频免费在线看| 秋霞在线观看毛片| 高清黄色对白视频在线免费看| 亚洲高清免费不卡视频| 麻豆成人av视频| 男女边吃奶边做爰视频| 在线观看一区二区三区激情| 搡女人真爽免费视频火全软件| 一边摸一边做爽爽视频免费| 美女大奶头黄色视频| 午夜福利影视在线免费观看| 欧美日本中文国产一区发布| 亚洲欧洲日产国产| 国产亚洲av片在线观看秒播厂| 欧美亚洲 丝袜 人妻 在线| 日韩免费高清中文字幕av| 欧美+日韩+精品| 夫妻午夜视频| 亚洲国产毛片av蜜桃av| 国产一区二区三区av在线| 精品人妻熟女av久视频| 成人亚洲精品一区在线观看| 久久ye,这里只有精品| av.在线天堂| 午夜激情福利司机影院| 女人久久www免费人成看片| 夫妻午夜视频| 久久久精品免费免费高清| 日本猛色少妇xxxxx猛交久久| 丝袜喷水一区| 最新中文字幕久久久久| 韩国高清视频一区二区三区| 欧美97在线视频| 一级片'在线观看视频| 久久久久久久亚洲中文字幕| 大码成人一级视频| 日韩视频在线欧美| 黑人高潮一二区| 老女人水多毛片| 亚洲欧洲国产日韩| 狠狠精品人妻久久久久久综合| 91国产中文字幕| 妹子高潮喷水视频| av免费观看日本| 狠狠精品人妻久久久久久综合| 少妇丰满av| 欧美 日韩 精品 国产| 男人操女人黄网站| 国产精品久久久久久精品古装| 欧美亚洲 丝袜 人妻 在线| 黄色视频在线播放观看不卡| 成人毛片a级毛片在线播放| 男人爽女人下面视频在线观看| 精品少妇黑人巨大在线播放| 91久久精品国产一区二区三区| 国产免费又黄又爽又色| 一区二区av电影网| 免费久久久久久久精品成人欧美视频 | 91精品国产九色| 最新的欧美精品一区二区| 久久久精品94久久精品| 中文字幕av电影在线播放| 久久精品国产a三级三级三级| 婷婷色麻豆天堂久久| 在现免费观看毛片| 大陆偷拍与自拍| 午夜91福利影院| 高清视频免费观看一区二区| 亚洲av电影在线观看一区二区三区| 久久这里有精品视频免费| 欧美精品国产亚洲| 午夜影院在线不卡| 啦啦啦在线观看免费高清www| 丝袜喷水一区| 爱豆传媒免费全集在线观看| 久久久a久久爽久久v久久| 免费不卡的大黄色大毛片视频在线观看| 天天影视国产精品| 成人国语在线视频| 秋霞在线观看毛片| 22中文网久久字幕| 激情五月婷婷亚洲| 在线观看国产h片| 秋霞伦理黄片| 欧美日韩视频精品一区| videossex国产| 伦理电影大哥的女人| 日日摸夜夜添夜夜爱| 日本黄色日本黄色录像| 人人澡人人妻人| 久久久欧美国产精品| 免费大片18禁| 国产精品.久久久| 丝袜美足系列| 色94色欧美一区二区| 毛片一级片免费看久久久久| 国产精品不卡视频一区二区| av福利片在线| 日日摸夜夜添夜夜爱| 九草在线视频观看| av卡一久久| 欧美日韩综合久久久久久| 亚洲欧美日韩卡通动漫| 亚洲av二区三区四区| 午夜91福利影院| 老司机影院毛片| 亚洲精品国产av成人精品| 男女免费视频国产| 春色校园在线视频观看| 天堂中文最新版在线下载| 人妻少妇偷人精品九色| 99国产综合亚洲精品| 一本久久精品| 全区人妻精品视频| 亚洲综合精品二区| 日本黄大片高清| 黄片播放在线免费| 性色av一级| 国产黄频视频在线观看| a级毛色黄片| 久久久国产精品麻豆| 视频中文字幕在线观看| 免费看光身美女| 啦啦啦视频在线资源免费观看| 国产精品.久久久| 亚洲综合精品二区| 成年美女黄网站色视频大全免费 | 国产又色又爽无遮挡免| 自拍欧美九色日韩亚洲蝌蚪91| 考比视频在线观看| 天堂8中文在线网| 欧美日韩综合久久久久久| 国产片特级美女逼逼视频| 午夜日本视频在线| av不卡在线播放| 这个男人来自地球电影免费观看 | 久久精品熟女亚洲av麻豆精品| 激情五月婷婷亚洲| 在线精品无人区一区二区三| 午夜福利视频精品| 免费观看在线日韩| 欧美xxxx性猛交bbbb| 一区在线观看完整版| 国产日韩欧美亚洲二区| 日本猛色少妇xxxxx猛交久久| 麻豆成人av视频| av卡一久久| 亚洲av日韩在线播放| 在线免费观看不下载黄p国产| 国产一区二区三区av在线| 99久久人妻综合| 高清av免费在线| 麻豆成人av视频| 91久久精品国产一区二区成人| 成人亚洲欧美一区二区av| 国产精品久久久久成人av| av网站免费在线观看视频| 天堂俺去俺来也www色官网| 欧美bdsm另类| 久久久a久久爽久久v久久| 伦理电影免费视频| 日韩,欧美,国产一区二区三区| 久久午夜福利片| 在线天堂最新版资源| a 毛片基地| 一区二区三区精品91| 欧美成人午夜免费资源| 桃花免费在线播放| 国产不卡av网站在线观看| 在线观看免费高清a一片| 久久久a久久爽久久v久久| 狂野欧美激情性xxxx在线观看| av有码第一页| 女性被躁到高潮视频| 成人手机av| 极品少妇高潮喷水抽搐| 国产国拍精品亚洲av在线观看| 91精品一卡2卡3卡4卡| 啦啦啦中文免费视频观看日本| 人体艺术视频欧美日本| 午夜激情av网站| av免费观看日本| av在线老鸭窝| 久久ye,这里只有精品| 国产成人精品福利久久| 2021少妇久久久久久久久久久| 国产精品国产三级专区第一集| 亚洲精品乱码久久久v下载方式| 高清av免费在线|