TWI809487B - 以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法 - Google Patents
以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法 Download PDFInfo
- Publication number
- TWI809487B TWI809487B TW110131859A TW110131859A TWI809487B TW I809487 B TWI809487 B TW I809487B TW 110131859 A TW110131859 A TW 110131859A TW 110131859 A TW110131859 A TW 110131859A TW I809487 B TWI809487 B TW I809487B
- Authority
- TW
- Taiwan
- Prior art keywords
- harmonic
- frequency
- ratio
- module
- peak
- Prior art date
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims description 23
- 230000037361 pathway Effects 0.000 title description 2
- 238000004891 communication Methods 0.000 claims abstract description 34
- 238000011156 evaluation Methods 0.000 claims abstract description 34
- 238000006243 chemical reaction Methods 0.000 claims abstract description 33
- 230000005540 biological transmission Effects 0.000 claims abstract description 21
- 238000010801 machine learning Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 230000017531 blood circulation Effects 0.000 claims description 48
- 238000012549 training Methods 0.000 claims description 32
- 238000012706 support-vector machine Methods 0.000 claims description 27
- 230000000875 corresponding effect Effects 0.000 description 40
- 230000035945 sensitivity Effects 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000012360 testing method Methods 0.000 description 9
- 208000031481 Pathologic Constriction Diseases 0.000 description 7
- 238000000502 dialysis Methods 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 7
- 208000037804 stenosis Diseases 0.000 description 7
- 230000036262 stenosis Effects 0.000 description 7
- 230000036772 blood pressure Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000001631 haemodialysis Methods 0.000 description 5
- 230000000322 hemodialysis Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 206010003226 Arteriovenous fistula Diseases 0.000 description 4
- 238000002790 cross-validation Methods 0.000 description 4
- 230000002792 vascular Effects 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 206010016717 Fistula Diseases 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000003890 fistula Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000002966 stenotic effect Effects 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 208000005189 Embolism Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 208000020832 chronic kidney disease Diseases 0.000 description 1
- 208000022831 chronic renal failure syndrome Diseases 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000002316 cosmetic surgery Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 238000010832 independent-sample T-test Methods 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/026—Measuring blood flow
- A61B5/0265—Measuring blood flow using electromagnetic means, e.g. electromagnetic flowmeter
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
- A61B5/0015—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
- A61B5/002—Monitoring the patient using a local or closed circuit, e.g. in a room or building
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/026—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/0507—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves using microwaves or terahertz waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/02042—Determining blood loss or bleeding, e.g. during a surgical procedure
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Cardiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Hematology (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Vascular Medicine (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Radar Systems Or Details Thereof (AREA)
- Complex Calculations (AREA)
Abstract
一種以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統，包含一擷取裝置及一評估裝置。該擷取裝置朝該廔管發射一無線電載波、接收由該無線電載波反射而成的信號並輸出為一傳輸信號。該評估裝置包括一通訊模組、一時頻域轉換模組、一比值運算模組及一機器學習模組。該評估裝置將該傳輸信號進行時頻域轉換，並取其中強度最大之頻率作為一基頻，運算該基頻的至少一諧波頻率區間所對應的強度值為至少一諧波強度值，運算其中一該諧波強度值與該基頻對應之強度值的比值或運算其中二該諧波強度值之比值，根據運算之比值輸出一評估結果。
Description
本發明是有關於一種醫療評估系統與方法，特別是指一種以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法。
慢性腎衰竭的病人必須先接受動靜脈廔管(arteriovenous fistula)的手術，才能進行洗腎治療，而由於動靜脈廔管的血流順暢度與洗腎效果密切相關，一旦廔管栓塞導致洗腎效果不佳或是無法洗腎，將造成病人須立即接受緊急處理甚或住院治療，因此，如何監控廔管通暢度以確保洗腎效果即為醫界重要研究目標。
目前醫院偵測廔管通路狀況，會使用Transonic公司之HD03血液透析監視儀，其使用方式是將兩根針頭插入廔管適當處，以量測血流量並評估廔管通路功能健全度，其量測精確度高，但儀器價格昂貴，且每次使用之針頭、導管等耗材亦不便宜，因此，不適合作為常規之廔管通路監測方式。
因此，本發明之目的，即在提供一種適用於常規監測廔管通路狀態的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統。
於是，本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統，適用於量測評估一受試者的一廔管的通路狀態，包含一擷取裝置及一評估裝置。
該擷取裝置包括一發射天線、一接收天線、一發射模組，及一接收模組，該發射模組經該發射天線朝該廔管發射一無線電載波，該接收模組經該接收天線接收由該無線電載波反射而成的一回波信號，該接收模組並將該回波信號輸出為一傳輸信號。
該評估裝置包括一通訊模組、一時頻域轉換模組、一比值運算模組，及一機器學習模組，該通訊模組信號連接該接收模組，接收該傳輸信號並輸出為一通訊信號，該時頻域轉換模組信號連接該通訊模組，接收該通訊信號並進行時頻域轉換以輸出一頻域資訊，該比值運算模組接收該頻域資訊並取該頻域資訊中強度最大之頻率作為一基頻，並運算該基頻的至少一諧波頻率區間所對應的強度值為至少一諧波強度值，運算其中一該諧波強度值與該基頻對應之強度值的比值或運算其中二該諧波強度值之比值，該機器學習模組根據運算之比值輸出一評估結果。
因此，本發明之目的，即在提供一種適用於常規監測廔管通路狀態的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法。
於是，本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法，適用於量測評估一受試者的一廔管的通路狀態，包含以下步驟：
(A)朝該廔管發射一無線電載波，並接收由該無線電載波反射而成的一回波信號。
(B)將該回波信號輸出為一傳輸信號。
(C)將該傳輸信號輸出為一通訊信號，並將該通訊信號進行時頻域轉換以輸出一頻域資訊。
(D)取該頻域資訊中強度最大之頻率作為一基頻，並運算該基頻的至少一諧波頻率區間所對應的強度值為至少一諧波強度值，運算其中一該諧波強度值與該基頻對應之強度值的比值或運算其中二該諧波強度值之比值。
(E)根據運算之比值輸出一評估結果。
本發明之功效在於：藉由發射該無線電載波並接收該回波信號，再搭配進行時頻域轉換取得該頻域資訊，及使用運算之比值供該機器學習模組進行評估，可以得到關於該廔管之血流量的預測結果，且本發明實施架構之價格實惠、操作方便、為非侵入性檢查，且無耗材產生，因此，十分適合作為醫院常規之廔管通路監測方式或是病患的居家護理監控。
參閱圖1與圖2，本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統之一實施例，適用於量測評估一受試者的一廔管91(例如，動靜脈廔管(arteriovenous fistula))的通路狀態，包含一擷取裝置2及一評估裝置3。
該擷取裝置2適用於對應於該廔管91裝設於該受試者身上，並適用於靠近該受試者之皮膚92設置，例如，將該擷取裝置2貼靠於該受試者鄰近該動靜脈廔管91的靜脈端出口處之皮膚92上。該擷取裝置2包括一發射天線21、一接收天線22、一發射模組23，及一接收模組24，該發射模組23經該發射天線21朝該廔管91發射一無線電載波，該接收模組24經該接收天線22接收由該無線電載波反射而成的一回波信號，該接收模組24並將該回波信號輸出為一傳輸信號。
該發射模組23具有一頻率可調方波產生電路231、一電連接該頻率可調方波產生電路231與該發射天線21的發射脈衝產生電路232，及一電連接該發射脈衝產生電路232的延遲脈衝產生電路233。
該接收模組24具有一電連接該接收天線22的解調與濾波電路241、一電連接該解調與濾波電路241的類比數位轉換電路242，及一電連接該類比數位轉換電路242的傳輸電路243。該傳輸電路243可使用短距離無線傳輸技術實施，例如使用藍牙(Bluetooth)技術實施。
該擷取裝置2之細節設置可見中華民國發明專利公告號I702938號「瘻管評估系統與方法」及發明專利申請號第109128981號「瘻管評估系統與方法」，因此，於本案中不再贅述。
該評估裝置3包括一通訊模組31、一時頻域轉換模組32、一比值運算模組33，及一機器學習模組34。其中，可使用一智慧型手機或平板等可攜式電子裝置實施該通訊模組31、該時頻域轉換模組32與該比值運算模組33，並使用一伺服器(server)實施該機器學習模組34，或亦可使用一伺服器直接實施該評估裝置3整體，於本實施例中，是以一智慧型手機搭配一伺服器進行說明，該手機可以使用預先安裝的APP(行動應用程式，mobile application)執行該通訊模組31、該時頻域轉換模組32與該比值運算模組33之功能運算。
該通訊模組31與該傳輸電路243間可使用藍牙技術進行無線傳輸，該通訊模組31接收該傳輸信號並輸出為一通訊信號。
該時頻域轉換模組32信號連接該通訊模組31，接收該通訊信號並進行數位濾波及時頻域轉換以輸出一頻域資訊。該時頻域轉換模組32可將該通訊信號內含的該廔管91的位移資訊儲存於手機中，且可顯示於該手機之螢幕以供使用者觀看。其中，該時頻域轉換模組32之數位濾波功能可例如使用FIR(有限脈衝響應，Finite impulse response)濾波器實施，時頻域轉換功能則例如使用快速傅立葉轉換(Fast Fourier Transform，縮寫為FFT)實施。該時頻域轉換模組32可濾取一預定頻段以輸出為該頻域資訊，於本實施例中，是將該通訊信號濾取0.2Hz~10Hz之頻段進行時頻域轉換以輸出該頻域資訊。
該比值運算模組33接收該頻域資訊並取該頻域資訊中強度最大之頻率作為一基頻，並運算該基頻的至少一諧波頻率區間所對應的強度值為至少一諧波強度值，運算其中一該諧波強度值與該基頻對應之強度值的比值或運算其中二該諧波強度值之比值。例如，該比值運算模組33於該預定頻段(0.2Hz~10Hz)中取該基頻的n-1個諧波頻率區間所分別對應的強度值作為n-1個諧波強度值，並分別運算
作為n-1個諧波峰值比值，其中，n為一預定值，且為大於1的正整數，
為該基頻對應之強度值(即，第1個諧波頻率區間所對應之強度值)，
為第n-1個諧波頻率區間所對應的諧波強度值，
為第n個諧波頻率區間所對應的諧波強度值。
該機器學習模組34預存至少一支援向量機341(support vector machine，縮寫為SVM)，該支援向量機341根據該等n-1個諧波峰值比值輸出該評估結果。
其中，該比值運算模組33與該機器學習模組34間(例如，該手機與該伺服器間)可以無線通訊網路技術(例如：GSM、Wi-Fi、藍牙技術等)進行訊號傳輸。該支援向量機341需先使用例如K折交叉驗證(K-fold Cross-Validation)方式進行訓練，以得到對應於一血流參數的預測模型。
參閱圖1、圖2、圖6及圖7，實際應用時，可搭配本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法之一實施例進行測試，該以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法包含以下步驟：
步驟51：朝該廔管91發射一無線電載波，並接收由該無線電載波反射而成的一回波信號。
步驟52：將該回波信號輸出為一傳輸信號。
步驟53：將該傳輸信號輸出為一通訊信號，並將該通訊信號進行時頻域轉換以輸出一頻域資訊。
其中，本實施例中，以每秒64次的頻率取樣1分鐘的該通訊信號，並先將該通訊信號進行數位濾波以濾取該預定頻段再進行處理，例如，可濾取0.2Hz~10Hz之頻段進行時頻域轉換以輸出該頻域資訊。請參照圖3~圖6，該通訊信號為時域訊號，其中一波形約如圖3或圖4所示，該通訊信號經由數位濾波與時頻域轉換後，則會形成如圖5或圖6所示之頻譜圖。
參閱圖1、圖2、圖6及圖7，步驟54：取該頻域資訊中強度最大之頻率作為一基頻，並運算該基頻的至少一諧波頻率區間所對應的強度值為至少一諧波強度值，運算其中一該諧波強度值與該基頻對應之強度值的比值或運算其中二該諧波強度值之比值。
其中，取該基頻的n-1個諧波頻率區間所分別對應的強度值作為n-1個諧波強度值，並分別運算
作為n-1個諧波峰值比值，n為一預定值，且為大於1的正整數，
為該基頻對應之強度值，
為第n-1個諧波頻率區間所對應的諧波強度值，
為第n個諧波頻率區間所對應的諧波強度值。值得一提的是，亦可僅運算部分的諧波峰值比值，例如，僅運算
，或僅運算
，或僅運算
等，或亦可使用
作為n-1個諧波峰值比值，不以此為限。
例如，以n=2說明，取該基頻的一個諧波頻率區間(例如，該基頻之二次諧波與鄰近區間，或亦可取三次諧波與鄰近區間、或四次諧波與鄰近區間、或五次諧波與鄰近區間、或六次諧波與鄰近區間…等)所對應的強度值作為該諧波強度值，並運算
作為該諧波峰值比值，
為該基頻對應之強度值，
為該基頻的第2個諧波頻率區間(即二次諧波與鄰近區間)所對應的諧波強度值。
以n=6說明，取該基頻的5個諧波頻率區間(例如，取該基頻之二次諧波與鄰近區間、三次諧波與鄰近區間、四次諧波與鄰近區間、五次諧波與鄰近區間、六次諧波與鄰近區間)所分別對應的強度值作為5個諧波強度值，並分別運算
、
、
、
、
作為5個諧波峰值比值，
為該基頻對應之強度值，
、
、
、
、
分別為第2~6個諧波頻率區間(即二次諧波與鄰近區間、三次諧波與鄰近區間、四次諧波與鄰近區間、五次諧波與鄰近區間、六次諧波與鄰近區間)所對應的諧波強度值。
其中，可參照圖6，於本實施例中，由於該預定頻段為0.2Hz ~10Hz，而強度最大之頻率(該基頻)一般落在1Hz~2Hz之間，因此，n值一般定為小於10的正整數(或是定為小於8的正整數)，於本實施例中，該比值運算模組33是以n=6進行運算。
其中，該等諧波強度值
~
之取法為，先取該基頻的2~n倍值分別作為n-1個參考諧波頻率，並以每一參考諧波頻率與其上下一預定區間作為該諧波頻率區間，取該諧波頻率區間中的強度最大值作為對應的該諧波強度值。例如，以取第2個諧波強度值說明，假設該基頻之頻率為1.3Hz，該預定區間為0.05Hz，則先取該基頻2倍值2.6Hz作為參考諧波頻率，再以2.6Hz與其上下該預定區間2.55Hz~2.65Hz作為該諧波頻率區間，再以2.55Hz~2.65Hz區間中的強度最大值作為對應的該諧波強度值。
步驟55：根據運算之比值輸出一評估結果。
其中，該支援向量機341根據該等n-1個諧波峰值比值輸出該評估結果，該評估結果為該支援向量機341的預測類別，一般為1或0之輸出。
於實際應用時，該機器學習模組34中可以內建一或多個支援向量機341，每一支援向量機341對應不同的該血流參數，例如，該機器學習模組34中可以只內建1個對應於600毫升/分鐘(ml/min)之流量的支援向量機341，或是可內建有7個支援向量機341，該等支援向量機341分別對應於500、600、650、750、900、1200、3500 ml/min之流量，或是可內建上述7個支援向量機341，但於販售時根據客戶的選擇而僅開放其中一個或多個支援向量機341的使用權限供客戶使用。
當要進行評估時，先以步驟51~步驟54量測該受試者以獲得該等諧波峰值比值
，接著，將該等諧波峰值比值
輸入所選擇(或具有使用權限)的該支援向量機341。例如，選擇對應於600 ml/min流量的該支援向量機341時，若該支援向量機341輸出之該評估結果為1，則表示該受試者之廔管91的預測血流量大於600 ml/min，若輸出之該評估結果為0，則表示該受試者之廔管91的預測血流量低於等於600 ml/min。醫護人員可根據需求而選擇對應合適該血流參數的該支援向量機341，並根據輸出的該評估結果判斷是否需要進行進一步的檢查。
該以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法，還包含下列訓練(training)步驟：
步驟61：收集至少一訓練資料組，每一訓練資料組對應不同之一血流參數且具有複數筆分別對應複數預試者的訓練資料，每一筆訓練資料具有一判定類別及n-1個諧波峰值比值。
其中，需先根據需求定義該至少一血流參數，並根據所定義的該血流參數收集該訓練資料組。例如，當定義單一個該血流參數(例如，600 ml/min)時，則須收集一組對應600 ml/min血流量的該訓練資料組，當定義多個該血流參數(例如，600 ml/min與750 ml/min)時，則須收集一組對應600 ml/min血流量的該訓練資料組、與一組對應750 ml/min血流量的該訓練資料組。以此類推，不再贅述。
以一組對應600 ml/min血流量的該訓練資料組作為說明，該訓練資料組具有複數筆分別對應複數預試者(於本實施例中，為45位)的訓練資料，每一筆訓練資料具有該判定類別及n-1個諧波峰值比值，該判定類別之資料為0或1，0表示低於等於該血流參數，1表示高於該血流參數。該等n-1個諧波峰值比值是使用步驟51~54對該等預試者進行量測所得。該判定類別是根據該血流參數與使用血流流量黃金標準儀器(例如，Transonic公司之HD03血液透析監視儀)量測該等預試者之血流量而得。例如，當其中一筆訓練資料為[0,
]時，表示該筆訓練資料所對應的該預試者所量測的血流量低於等於600 ml/min，而
為該預試者經步驟51~54量測及運算所得之諧波峰值比值。
步驟62：根據該至少一訓練資料組訓練至少一支援向量機341，每一支援向量機341對應不同的該血流參數。
其中，該支援向量機341可使用例如K折交叉驗證(K-fold Cross-Validation)方式進行訓練，以得到對應於該血流參數的預測模型。於本實施例中，設定K值為10，即，將該訓練資料組中的多筆訓練資料分為10組，每組可均分到相同數量之訓練資料，或分配不同數量之訓練資料，以其中任意9組作為訓練組，剩餘1組作為測試組，完成第1次訓練後，由訓練組中挑出一組未曾作為測試組的資料作為下一輪的測試組，其餘9組作為下一輪的訓練組，接著進行下一輪的訓練，如此直到每一組資料都作過測試組，即完成該支援向量機341之訓練。
本發明之原理說明如下：
該擷取裝置2適用於發射該無線電載波穿透皮膚92後從該廔管91表面反射，該廔管91因血壓脈動而產生週期位移，此位移所造成的都普勒效應(Doppler effect)會改變該無線電載波的頻率並形成該回波信號，該回波信號的特徵與廔管91通路的狹窄有高度相關性，其原因可基於流體力學理論進行說明，當血液通過廔管91管徑狹窄處，因廔管91狹窄處的血流速度會增加並產生亂流，因此，在廔管91靜脈端出口處偵測到的訊號會包含震動弦波成份，造成訊號非線性失真現象。參閱圖3~圖6，圖3及圖5為健康者之其中一血壓波形的圖片及將所量測的訊號進行時頻域轉換等處理後輸出的該頻域資訊，圖4及圖6則為被判斷為該廔管91狹窄之病患的其中一血壓波形的圖片，及將所量測的訊號進行時頻域轉換等處理後輸出的該頻域資訊。將兩組圖片進行比較，可以看出兩組圖片在波形與頻譜強度分布上皆具有差異(圖5及圖6中之虛線為波形43、44的端點近似曲線)。並且，經由將兩組該頻域資訊之諧波峰值比值進行獨立樣本T檢定(independent T test)，於血流參數為600 ml/min時，於
的T檢定得到之P值(P-value)為0.008，於血流參數為650 ml/min時，於
的T檢定得到之P值為0.035，於血流參數為750 ml/min時，於
的T檢定得到之P值為0.041，上述之P值皆小於0.05，足證兩者(正常與廔管91狹窄)之諧波峰值比值間具有顯著差異。
請參下表1~3，為使用本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法之實施例對45個受試者進行血流評估的結果，表1~3分別是以血流參數為600、650、750 ml/min進行評估，所採用之基準真相(Ground Truth)為，使用Transonic公司之HD03血液透析監視儀量測該受試者的血流量，低於等於對應之該血流參數，則列為0，高於對應之該血流參數，則列為1。
表1
表2
表3
| 血流參數： 600 ml/min | |||
| 基準真相 | 評估結果 | ||
| 0 ( 低 ) | 1 ( 高 ) | ||
| 0 ( 低 ) | 34 | 0 | 100% (特異度) |
| 1 ( 高 ) | 1 | 10 | 90.9% (靈敏度) |
| 血流參數： 650 ml/min | |||
| 基準真相 | 評估結果 | ||
| 0 ( 低 ) | 1 ( 高 ) | ||
| 0 ( 低 ) | 29 | 0 | 100% (特異度) |
| 1 ( 高 ) | 1 | 15 | 93.8% (靈敏度) |
| 血流參數： 750 ml/min | |||
| 基準真相 | 評估結果 | ||
| 0 ( 低 ) | 1 ( 高 ) | ||
| 0 ( 低 ) | 24 | 0 | 100% (特異度) |
| 1 ( 高 ) | 1 | 20 | 95.2% (靈敏度) |
由表1~3可見，本實施例在血流參數為600、650、750 ml/min時，特異度(Specificity)皆高達100%，而靈敏度(Sensitivity)也分別高達90.9%、93.8%、95.2%，足證本實施例的血流評估確實具有極佳的準確度。
請參閱圖9及下表4，曲線45為本實施例根據血流參數500、600、650、750、900、1200、1800、3500 ml/min之特異度與靈敏度所繪製的接收者操作特徵曲線（receiver operating characteristic curve，或ROC曲線），其ROC曲線下方的面積(Area under the Curve of ROC，或稱AUC)為0.997，其值趨近於1，足證本實施例之預測正確率極高。
表4
| 血流參數 | 500 | 600 | 650 | 750 | 900 | 1200 | 1800 | 3500 |
| 特異度 | 100% | 100% | 100% | 100% | 88.2% | 54.5% | 0% | 0% |
| 靈敏度 | 83.3% | 90.9% | 93.8% | 95.2% | 100% | 100% | 100% | 100% |
參閱圖1及圖2，經由以上的說明，本實施例的功效如下：
一、藉由發射該無線電載波並接收該回波信號，再搭配進行時頻域轉換取得該頻域資訊，及使用該諧波頻率區間與該基頻對應之強度值的比值供該機器學習模組34進行評估，可以得到關於該廔管91之血流量之預測結果。
並且，由於本實施例是使用無線電載波方式進行量測，其實施架構之價格遠較Transonic公司之HD03血液透析監視儀實惠，操作也更為方便，且並非侵入性檢查、無耗材產生，因此，十分適合作為醫院常規之廔管91通路評估方式或是病患的居家護理監控，再由醫護人員根據所預測的廔管91之血流量，初步評估病患之廔管91是否可能異常，或是否需要再使用Transonic公司之HD03血液透析監視儀(或超音波、血管攝影等方式)進行更精確的檢查。藉此，可方便醫護人員隨時量測病患之廔管91的血流量，提高病患在未發生栓塞前能早期發現廔管91是否狹窄化的機率，以及早安排精密診斷，評估是否需進行動靜脈廔管91整型術。故本實施例能供醫護人員評估掌控洗腎效果，有效降低病患的再入院率、減少醫療支出及提高生活品質。
再者，根據本實施例之量測驗證結果，本實施例在血流參數為600、650、750 ml/min時，特異度皆高達100%，而靈敏度也分別高達90.9%、93.8%、95.2%，在血流參數為900、1200 ml/min時，靈敏度皆高達100%，在血流參數為1800、3500 ml/min時，特異度與靈敏度皆高達100%，AUC之值則高達0.997，足證本實施例確實具有極佳的預測準確度。
綜上所述，本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法，確實能達成本發明的目的。
惟以上所述者，僅為本發明之實施例而已，當不能以此限定本發明實施之範圍，凡是依本發明申請專利範圍及專利說明書內容所作之簡單的等效變化與修飾，皆仍屬本發明專利涵蓋之範圍內。
2:擷取裝置
21:發射天線
22:接收天線
23:發射模組
231:頻率可調方波產生電路
232:發射脈衝產生電路
233:延遲脈衝產生電路
24:接收模組
241:解調與濾波電路
242:類比數位轉換電路
243:傳輸電路
3:評估裝置
31:通訊模組
32:時頻域轉換模組
33:比值運算模組
34:機器學習模組
341:支援向量機
41~44:波形
45:曲線
51~55:步驟
61~62:步驟
91:廔管
92:皮膚
~
:第1~6個諧波頻率區間的諧波強度值
本發明之其他的特徵及功效，將於參照圖式的實施方式中清楚地呈現，其中：
圖1是本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統的一實施例應用於量測一廔管的示意圖；
圖2是該實施例的一方塊示意圖；
圖3是該實施例所量測的一健康者之血壓波形圖；
圖4是該實施例所量測的一被判斷為廔管狹窄之病患的血壓波形圖；
圖5是該實施例所產生的一相關於健康者的頻域資訊；
圖6是該實施例所產生的一相關於廔管狹窄之病患的頻域資訊；
圖7、8是本發明以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法的一實施例的二流程圖；及
圖9為該實施例的一接收者操作特徵曲線圖。
2:擷取裝置
21:發射天線
22:接收天線
23:發射模組
231:頻率可調方波產生電路
232:發射脈衝產生電路
233:延遲脈衝產生電路
24:接收模組
241:解調與濾波電路
242:類比數位轉換電路
243:傳輸電路
3:評估裝置
31:通訊模組
32:時頻域轉換模組
33:比值運算模組
34:機器學習模組
341:支援向量機
Claims (8)
- 一種以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統，適用於量測評估一受試者的一廔管的通路狀態，包含： 一擷取裝置，包括一發射天線、一接收天線、一發射模組，及一接收模組，該發射模組經該發射天線朝該廔管發射一無線電載波，該接收模組經該接收天線接收由該無線電載波反射而成的一回波信號，該接收模組並將該回波信號輸出為一傳輸信號；及 一評估裝置，包括一通訊模組、一時頻域轉換模組、一比值運算模組，及一機器學習模組，該通訊模組信號連接該接收模組，接收該傳輸信號並輸出為一通訊信號，該時頻域轉換模組信號連接該通訊模組，接收該通訊信號並進行時頻域轉換以輸出一頻域資訊，該比值運算模組接收該頻域資訊並取該頻域資訊中強度最大之頻率作為一基頻，並運算該基頻的至少一諧波頻率區間所對應的強度值為至少一諧波強度值，運算其中一該諧波強度值與該基頻對應之強度值的比值或運算其中二該諧波強度值之比值，該機器學習模組根據運算之比值輸出一評估結果。
- 如請求項1所述的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統，其中，該比值運算模組取該基頻的n-1個諧波頻率區間所分別對應的強度值作為n-1個諧波強度值，並分別運算 作為n-1個諧波峰值比值，其中，n為一預定值，且為大於1的正整數， 為該基頻對應之強度值， 為第n-1個諧波頻率區間所對應的諧波強度值， 為第n個諧波頻率區間所對應的諧波強度值，該機器學習模組根據該等n-1個諧波峰值比值輸出該評估結果。
- 如請求項2所述的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統，其中，該時頻域轉換模組濾取一預定頻段輸出為該頻域資訊，且該比值運算模組於該預定頻段中取該等n-1個諧波頻率區間。
- 如請求項2所述的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統，其中，該機器學習模組預存至少一支援向量機，該至少一支援向量機根據該等n-1個諧波峰值比值輸出該評估結果。
- 一種以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法，適用於量測評估一受試者的一廔管的通路狀態，包含以下步驟： (A) 朝該廔管發射一無線電載波，並接收由該無線電載波反射而成的一回波信號； (B) 將該回波信號輸出為一傳輸信號； (C) 將該傳輸信號輸出為一通訊信號，並將該通訊信號進行時頻域轉換以輸出一頻域資訊； (D) 取該頻域資訊中強度最大之頻率作為一基頻，並運算該基頻的至少一諧波頻率區間所對應的強度值為至少一諧波強度值，運算其中一該諧波強度值與該基頻對應之強度值的比值或運算其中二該諧波強度值之比值；及 (E) 根據運算之比值輸出一評估結果。
- 如請求項5所述的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法，其中： 於步驟(D)中，取該基頻的n-1個諧波頻率區間所分別對應的強度值作為n-1個諧波強度值，並分別運算 作為n-1個諧波峰值比值，其中，n為一預定值，且為大於1的正整數， 為該基頻對應之強度值， 為第n-1個諧波頻率區間所對應的諧波強度值， 為第n個諧波頻率區間所對應的諧波強度值； 於步驟(E)中，根據該等n-1個諧波峰值比值輸出該評估結果。
- 如請求項6所述的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法，其中，於步驟(D)中，取該基頻的2~n倍值分別作為n-1個參考諧波頻率，並以每一參考諧波頻率與其上下一預定區間作為該諧波頻率區間，取該諧波頻率區間中的強度最大值作為對應的該諧波強度值。
- 如請求項6所述的以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之方法，還包含下列步驟； (F) 收集至少一訓練資料組，每一訓練資料組對應不同之一血流參數且具有複數筆分別對應複數預試者的訓練資料，每一筆訓練資料具有一判定類別及n-1個諧波峰值比值，該等n-1個諧波峰值比值是使用步驟(A)~(D)對該等預試者進行量測所得；及 (G) 根據該至少一訓練資料組訓練至少一支援向量機，每一支援向量機對應不同的該血流參數。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW110131859A TWI809487B (zh) | 2021-08-27 | 2021-08-27 | 以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法 |
| US17/890,015 US20230079269A1 (en) | 2021-08-27 | 2022-08-17 | Method and system for evaluating status of fistula |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW110131859A TWI809487B (zh) | 2021-08-27 | 2021-08-27 | 以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW202308555A TW202308555A (zh) | 2023-03-01 |
| TWI809487B true TWI809487B (zh) | 2023-07-21 |
Family
ID=85480117
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW110131859A TWI809487B (zh) | 2021-08-27 | 2021-08-27 | 以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US20230079269A1 (zh) |
| TW (1) | TWI809487B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20210212616A1 (en) * | 2017-03-29 | 2021-07-15 | Graftworx, Inc. | Wearable device with multimodal diagnostics |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI721395B (zh) * | 2018-04-21 | 2021-03-11 | 美商腎通生醫科技公司 | 動靜脈瘻管狹窄偵測系統及其方法及感測裝置 |
-
2021
- 2021-08-27 TW TW110131859A patent/TWI809487B/zh active
-
2022
- 2022-08-17 US US17/890,015 patent/US20230079269A1/en not_active Abandoned
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20210212616A1 (en) * | 2017-03-29 | 2021-07-15 | Graftworx, Inc. | Wearable device with multimodal diagnostics |
Also Published As
| Publication number | Publication date |
|---|---|
| US20230079269A1 (en) | 2023-03-16 |
| TW202308555A (zh) | 2023-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101577343B1 (ko) | 혈압 추정 방법 및 장치 | |
| CN104382571B (zh) | 一种基于桡动脉脉搏波传导时间的测量血压方法及装置 | |
| CN102988036B (zh) | 测量脉率的方法 | |
| US20190000326A1 (en) | Systems and methods for filtering noise and analyzing venous waveform signals | |
| CN112998690A (zh) | 一种基于脉搏波多特征融合的呼吸率提取方法 | |
| CN105105734B (zh) | 一种基于心音信号的无创连续血压测量装置 | |
| WO2012163738A1 (en) | Monitoring stenosis formation in an arteriovenous access | |
| Jaros et al. | Novel phonocardiography system for heartbeat detection from various locations | |
| CN107920768A (zh) | 用于测量由颈动脉中的涡流产生的声频的无创方法 | |
| TWI572327B (zh) | 利用音頻訊號判斷流體管路窄化的檢測裝置、電腦程式產品及電腦可讀取媒體 | |
| Panda et al. | Flexible, skin coupled microphone array for point of care vascular access monitoring | |
| CN117958770A (zh) | 一种动静脉内瘘血管健康监测方法及系统 | |
| JP2014530642A (ja) | 動脈の閉塞／再開を検出するための方法及び装置並びに収縮期血圧を測定するためのシステム | |
| CN105725983B (zh) | 一种周围动脉硬化的早期筛查方法及系统 | |
| TWI809487B (zh) | 以快速傅立葉轉換頻譜峰值比例評估廔管通路狀態之系統及方法 | |
| US20220000403A1 (en) | Neurotransmitter imbalance detection system and method of detecting a neurotransmitter imbalance | |
| Yılmaz et al. | Analysis of the Doppler signals using largest Lyapunov exponent and correlation dimension in healthy and stenosed internal carotid artery patients | |
| Pratiwi et al. | A review of equipment and signal processing of the automated auscultation for blood pressure measurement | |
| CN119488275B (zh) | 基于电子柯氏音法的智能血压测量方法和系统 | |
| CN106388788A (zh) | 一种兼具心功能信息采集的脉搏波诊断系统 | |
| CN111603151B (zh) | 一种基于时频联合分析的无创血液成分检测方法及系统 | |
| Ochi et al. | Non-contact Blood Pressure Estimation By Microwave Reflection Employing Machine Learning | |
| TWI755826B (zh) | 瘻管評估系統與方法 | |
| RU2656560C1 (ru) | Способ оценки риска осложнений заболеваний сердечно-сосудистой системы с сочетанной патологией | |
| KR20080043545A (ko) | 휴대용 맥진기 및 그 맥진방법 |