TWI667997B - 建構３ｄ偵測模型之心音處理方法及其系統 - Google Patents
建構３ｄ偵測模型之心音處理方法及其系統 Download PDFInfo
- Publication number
- TWI667997B TWI667997B TW106143404A TW106143404A TWI667997B TW I667997 B TWI667997 B TW I667997B TW 106143404 A TW106143404 A TW 106143404A TW 106143404 A TW106143404 A TW 106143404A TW I667997 B TWI667997 B TW I667997B
- Authority
- TW
- Taiwan
- Prior art keywords
- heart sound
- data
- model
- constructing
- stft
- Prior art date
Links
- 230000005236 sound signal Effects 0.000 claims abstract description 60
- 238000003672 processing method Methods 0.000 claims abstract description 23
- 238000001514 detection method Methods 0.000 claims description 63
- 230000002861 ventricular Effects 0.000 claims description 15
- 230000000747 cardiac effect Effects 0.000 claims description 14
- 206010010356 Congenital anomaly Diseases 0.000 claims description 12
- 230000001746 atrial effect Effects 0.000 claims description 8
- 238000000926 separation method Methods 0.000 claims description 8
- 238000013500 data storage Methods 0.000 claims description 7
- 230000001575 pathological effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000002604 ultrasonography Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 20
- 230000007170 pathology Effects 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000000034 method Methods 0.000 description 6
- 206010002915 Aortic valve incompetence Diseases 0.000 description 4
- 206010002906 aortic stenosis Diseases 0.000 description 4
- 201000002064 aortic valve insufficiency Diseases 0.000 description 4
- 208000019622 heart disease Diseases 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 208000035478 Interatrial communication Diseases 0.000 description 3
- 208000001910 Ventricular Heart Septal Defects Diseases 0.000 description 3
- 208000013914 atrial heart septal defect Diseases 0.000 description 3
- 206010003664 atrial septal defect Diseases 0.000 description 3
- 201000003130 ventricular septal defect Diseases 0.000 description 3
- 208000035211 Heart Murmurs Diseases 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 208000025339 heart septal defect Diseases 0.000 description 2
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000000709 aorta Anatomy 0.000 description 1
- 210000003157 atrial septum Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 210000003709 heart valve Anatomy 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
一種心音處理方法包含：將一心音訊號進行切割，以獲得數個心音訊號片段；將該心音訊號片段進行連續小波轉換，以獲得一CWT資料；將該心音訊號片段進行短時距傅立葉轉換，以獲得一STFT資料；利用該CWT資料及STFT資料計算產生至少一3D心音模型，且由該3D心音模型提供至少一3D心音模型資料。
Description
本發明係關於一種建構3D偵測模型之心音〔heart sound〕處理方法及其系統；特別是關於一種建構3D偵測模型之心臟病理〔cardiopathy〕偵測之心音處理方法及其系統，其用以提供一3D心音模型。
習用心臟病理偵測之心音處理之相關技術，例如：中華民國專利公開第201244691號之〝心音/心臟疾病風險之判別系統及其方法〞發明專利申請案，其揭示一種心音/心臟疾病風險之判別系統。該心音/心臟疾病風險之判別系統包含一訊號接收單元、一訊號處理單元、一儲存單元、一輸出單元及一顯示單元。
承上，該訊號接收單元用以接收一心音訊號A，且該心音訊號包含數個心音頻率。該訊號處理單元包含一第一運算單元、一濾波單元、一第二運算單元及一比對單元。該第一運算單元對心音訊號進行基於心音訊號的絕對值取自然對數之一數值後，將該數值與心音訊號之乘積之特殊函數運算〔例如，X=cAln|A'|，其中c可為任意數值或函數值〕，以產生一第一運算訊號。當A≠0時，A'=A；當A=0時，A'=R〔R≧1，且R為實數〕。該濾波單元對該第一運算訊號進行濾波，以產生一濾波訊號。該第二運算單元對該濾波訊號進行轉換運算，以產生一對應影像圖譜數據，且依該對應影像圖譜數據產生一影像圖譜資料。
承上，該儲存單元用以儲存有關心臟疾病之一預定心音圖譜資料，而該輸出單元為一無線傳送模組或有線傳輸模組。該比對單元進行比對該影像圖譜資料及預定心音圖譜資料，並產生一比對結果。接著，該輸出單元輸出該比對結果之訊號至該顯示單元，以顯示生理狀態及/或後續處理方法。
然而，前述第201244691號之心音處理技術必需先產生影像圖譜資料，再將該影像圖譜資料與預定心音圖譜資料進行比對。因此，前述第201244691號之心音處理技術仍存在進一步改良之需求。前述中華民國專利申請案僅為本發明技術背景之參考及說明目前技術發展狀態而已，其並非用以限制本發明之範圍。
有鑑於此，本發明為了滿足上述需求，其提供一種建構3D偵測模型之心音處理方法及其系統，將一心音訊號進行切割，以獲得數個心音訊號片段，且將該心音訊號片段進行連續小波轉換，以獲得一CWT資料，且將該心音訊號片段進行短時距傅立葉轉換，以獲得一STFT資料，利用該CWT資料或STFT資料計算產生至少一特徵值及至少一特徵向量或至少一3D心音模型，且由該特徵值產生一相關數值，或該特徵值與一樣本特徵值進行比對判別，以產生一比對結果，以便大幅提升心臟病理偵測之方便性。
本發明之主要目的係提供一種建構3D偵測模型之心音處理方法及其系統，其將一心音訊號進行切割，以獲得數個心音訊號片段，且將該心音訊號片段進行連續小波轉換，以獲得一CWT資料，且將該心音訊號片段進行短時距傅立葉轉換，以獲得一STFT資料，利用該CWT資料或STFT資料計算產生至少一特徵值及至少一特徵向
量或至少一3D心音模型，且由該特徵值產生一相關數值，或該特徵值與一樣本特徵值進行比對判別，以產生一比對結果，以達成提供由心音進行輔助心臟病理偵測之目的。
為了達成上述目的，本發明較佳實施例之建構3D偵測模型之心音處理方法包含：將一心音訊號進行切割，以獲得數個心音訊號片段；將該心音訊號片段進行連續小波轉換，以獲得一CWT資料；將該心音訊號片段進行短時距傅立葉轉換，以獲得一STFT資料；及利用該CWT資料及STFT資料計算產生至少一3D心音模型，且由該3D心音模型提供一3D心音模型資料。
本發明較佳實施例利用該CWT資料、STFT資料或3D心音模型資料計算產生至少一特徵值及至少一特徵向量，且由該特徵值產生一相關數值，或該特徵值與一樣本特徵值進行比對判別，以產生一比對結果。
本發明較佳實施例之該比對結果包含一心室分隔缺損、一心室分隔部之比對結果、一心房分隔缺損或一心房分隔部之比對結果。
本發明較佳實施例之該特徵值之相關數值包含一正相關性數值或一互相關性數值。
本發明較佳實施例之該心音訊號包含一第一心音S1訊號及一第二心音S2訊號，且該特徵值由該第二心音S2訊號產生。
本發明較佳實施例之該正相關性數值及互相關性數值用以建立一先天心室分隔模型。
本發明較佳實施例之該先天心室分隔模型為
一相關數值模型。
本發明較佳實施例利用該CWT資料、STFT資料與超音波資料進行預先比對，以尋找一最大頻率點聲音、一最高振幅位置聲音及至少一聲音間隔區間。
本發明較佳實施例採用皮爾遜積差相關係數計算該CWT資料及STFT資料。
本發明較佳實施例將數個該CWT資料及數個該STFT資料轉換成數個3D座標資料。
本發明較佳實施例之該3D心音模型提供一3D心音空間。
本發明較佳實施例之該3D心音模型顯示一3D轉換波形，以便立體顯示一心室分隔模型、一先天心室分隔模型、一心房分隔模型或一先天心房分隔模型。
為了達成上述目的，本發明另一較佳實施例之建構3D偵測模型之心音處理系統包含；一心音接收單元，其用以接收至少一心音訊號；一心音處理單元，其將該心音訊號進行切割，以獲得數個心音訊號片段，並將該心音訊號片段進行連續小波轉換，以獲得一CWT資料，且將該心音訊號片段進行短時距傅立葉轉換，以獲得一STFT資料；一資料儲存單元，其用以儲存至少一樣本特徵值；及一輸出單元，其用以輸出一相關數值；其中利用該CWT資料及STFT資料計算產生至少一3D心音模型，且由該3D心音模型提供一3D心音模型資料。
本發明較佳實施例利用該CWT資料、STFT資料或3D心音模型資料計算產生至少一特徵值及至少一
特徵向量，且由該特徵值產生一相關數值，或該特徵值與一樣本特徵值進行比對判別，以產生一比對結果。
本發明較佳實施例之該心音接收單元包含一第一接收單元及一第二接收單元。
本發明較佳實施例之該心音接收單元貼附於一人體皮膚表面上。
本發明較佳實施例之該特徵值之相關數值包含一正相關性數值或一互相關性數值。
本發明較佳實施例之該心音訊號包含一第一心音S1訊號及一第二心音S2訊號，且該特徵值由該第二心音S2訊號產生。
本發明較佳實施例之該正相關性數值及互相關性數值用以建立一先天心室分隔模型。
本發明較佳實施例之該先天心室分隔模型為一相關數值模型。
本發明較佳實施例之該心音訊號片段之頻率選取介於1Hz至100Hz之間。
本發明較佳實施例之該心音訊號片段包含一預定心音訊號數量，且每個該心音訊號片段之預定心音訊號數量相同。
本發明較佳實施例當該心音訊號片段包含一心音雜音訊號時，利用偵測該心音雜音訊號方式進行心臟病理偵測。
本發明較佳實施例將數個該CWT資料及數個該STFT資料轉換成數個3D座標資料。
本發明較佳實施例之該3D心音模型提供一3D心音空間。
本發明較佳實施例之該3D心音模型顯示一3D轉換波形，以便立體顯示一心室分隔模型、一先天心室分
隔模型、一心房分隔模型或一先天心房分隔模型。
S1‧‧‧步驟
S2‧‧‧步驟
S3‧‧‧步驟
S4A‧‧‧步驟
S4B‧‧‧步驟
10‧‧‧心音接收單元
20‧‧‧心音處理單元
21‧‧‧3D心音模型單元
30‧‧‧資料儲存單元
40‧‧‧輸出單元
第1圖：本發明較佳實施例之建構3D偵測模型之心音處理系統之方塊圖。
第2圖：本發明較佳實施例之建構3D偵測模型之心音處理方法之流程圖。
第2(A)圖：本發明較佳實施例之建構3D偵測模型之心音處理方法在轉換資料及建構3D心音模型之流程圖。
第3(a)及3(b)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統採用心音資料及心電圖資料之波形示意圖。
第4(A)及4(B)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第一心音資料及第二心音資料之波形示意圖。
第5(A)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第一病患之心音資料連續小波轉換後，顯示輕微大動脈回流及輕微大動脈狹窄狀況之波形示意圖。
第5(B)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第二病患之心音資料連續小波轉換後，顯示連續性心雜音及開放性動脈導管狀況之波形示意圖。
第6(A)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第一病患之心音資料短時距傅立葉轉換後，顯示輕微大動脈回流及大動脈狹窄輕微狀況之波形示意圖。
第6(B)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第二病患之心音資料短時距傅立葉轉
換後，顯示連續性心雜音及開放性動脈導管狀況之波形示意圖。
第7(A)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統將第一病患之原始心音資料、CWT資料及STFT資料進行比對之波形示意圖。
第7(B)圖：本發明較佳實施例之建構3D偵測模型之心音處理系統將第二病患之原始心音資料、CWT資料及STFT資料進行比對之波形示意圖。
第8圖：本發明較佳實施例之建構3D偵測模型之心音處理系統產生一3D心音模型之示意圖。
為了充分瞭解本發明，於下文將例舉較佳實施例並配合所附圖式作詳細說明，且其並非用以限定本發明。
本發明較佳實施例之建構3D偵測模型之心臟病理偵測之心音處理方法及其系統適用於各種心臟病理偵測裝置及其相關應用設備，例如：各類型居家照護〔home care〕系統、醫療器材自動控制系統〔例如：醫療檢查系統〕、遠距醫療照護系統或醫療教學系統，但其並非用以限定本發明之範圍。
第1圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統之方塊圖。請參照第1圖所示，舉例而言，本發明較佳實施例之用於心臟病理偵測之心音處理系統包含一心音接收單元10、一心音處理單元20、一3D心音模型單元21、一資料儲存單元30及一輸出單元40，且該心音接收單元10、資料儲存單元30及輸出單元40適當連接至該心音處理單元20及3D心音模型單元21。舉例而言，該心音接收單元10包含一第一接收單元及一第二接收單元，且該心音接收單元10貼附於一人體皮膚表面上。
請再參照第1圖所示，舉例而言，該資料儲存
單元30用以儲存數個超音波資料或心臟病理超音波資料樣本，以便比對連續小波轉換及短時距傅立葉轉換處理之心音資料，以建立時間與頻率的相關聯性樣本，如此日後據此可進行心臟病理偵測。
第2圖揭示本發明較佳實施例之建構3D偵測模型之心音處理方法之流程圖，其對應於第1圖之用於心臟病理偵測之心音處理系統。請參照第1及2圖所示，本發明較佳實施例之建構3D偵測模型之心音處理方法包含步驟S1：首先，將一心音訊號或一心音訊號片段及其它訊號〔例如：ECG，electrocardiogram〕以適當量測裝置〔例如：AUDICOR廠牌之量測裝置〕接收後進行切割，以獲得數個心音訊號片段〔fragment〕。
第3(a)及3(b)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統採用心音資料及心電圖資料之波形示意圖。請參照第3(a)圖所示，本發明較佳實施例之建構3D偵測模型之心音處理系統同時接收兩個接收單元之兩個心音訊號，因此其形成二個心音圖〔PCG，phonocardiogram〕。相對的，請參照第3(b)圖所示，本發明較佳實施例之用於心臟病理偵測之心音處理系統同時接收對應的心電圖〔ECG〕資料。
一般而言，心音包含第一心音S1、第二心音S2、第三心音S3及第四心音S4，且心音由心肌收縮、心臟瓣膜關閉、血液沖擊心室壁或大動脈壁所引起震動而產生的聲音。在心臟收縮及舒張時，產生第一心音S1及第二心音S2，因此其極易偵測〔例如：聽診器或儀器〕。至於第三心音S3，可在兒童時期或青少年時期可偵測，而第四心音S4則極不易偵測。
第4(A)及4(B)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第一心音資料及第二
心音資料之波形示意圖，其於心臟周圍區域不同部位以貼片〔patch〕偵測之心音資料。舉例而言，該心音訊號片段選擇包含一預定心音訊號數量〔例如：400筆〕，且每個該心音訊號片段之預定心音訊號數量選擇相同。另外，該心音訊號片段之頻率選取介於1Hz至100Hz之間或其它低頻區間。
請再參照第1及2圖所示，本發明較佳實施例之建構3D偵測模型之心音處理方法包含步驟S2：接著，將該心音訊號片段進行連續小波轉換〔continuous wavelet transformation，CWT〕處理，以獲得一CWT資料，其顯示一第一種時間與頻率的心臟病理波形。
第5(A)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第一病患之心音資料連續小波轉換後，顯示輕微大動脈回流〔AR，aortic regurgitation〕及輕微大動脈狹窄〔AS，aortic stenosis)〕狀況之波形示意圖。請參照第5(A)圖之所示，在第一病患之心音資料連續小波轉換後，獲得一第一連續小波轉換波形〔如第5(A)圖之上半部所示〕及一第二連續小波轉換波形〔如第5(A)圖之下半部所示〕。該第一連續小波轉換波形及第二連續小波轉換波形具有一第一最大振幅及一第二最大振幅，如兩個箭頭所示。
相對於第5(A)圖，第5(B)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第二病患之心音資料連續小波轉換後，顯示連續性心雜音〔continuous murmur〕及開放性動脈導管〔PDA，patent ductus arteriosum〕狀況之波形示意圖。請參照第5(B)圖之所示，在第二病患之心音資料連續小波轉換後，獲得一第一連續小波轉換波形〔如第5(B)圖之上半部所示〕及一第二連續小波轉換波形〔如第5(B)圖之下半部所示〕，且出現心音
雜音〔位於一主波之外〕，如較小波形所示。同樣的，該第一連續小波轉換波形及第二連續小波轉換波形具有一第一最大振幅及一第二最大振幅，如兩個箭頭所示，其位於該主波內。
請再參照第1及2圖所示，本發明較佳實施例之建構3D偵測模型之心音處理方法包含步驟S3：接著，將該心音訊號片段進行短時距傅立葉轉換〔short-time Fourier transformation，STFT〕處理，以獲得一STFT資料，其顯示一第二種時間與頻率的心臟病理波形。
第6(A)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第一病患之心音資料短時距傅立葉轉換後，顯示輕微大動脈回流及大動脈狹窄輕微狀況之波形示意圖。請參照第6(A)圖之所示，在第一病患之心音資料短時距傅立葉轉換後，獲得一第一短時距傅立葉轉換波形〔如第6(A)圖之上半部所示〕及一第二短時距傅立葉轉換波形〔如第6(A)圖之下半部所示〕。該第一短時距傅立葉轉換波形及第二短時距傅立葉轉換波形具有一第一最大振幅及一第二最大振幅，如兩個箭頭所示。
相對於第6(A)圖，第6(B)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統偵測第二病患之心音資料短時距傅立葉轉換後，顯示連續性心雜音及開放性動脈導管狀況之波形示意圖。請參照第6(B)圖之所示，
在第二病患之心音資料短時距傅立葉轉換後，獲得一第一短時距傅立葉轉換波形〔如第6(B)圖之上半部所示〕及一第二短時距傅立葉轉換波形〔如第6(B)圖之下半部所示〕，且出現心音雜音，如較小波形所示。同樣的，該第一短時距傅立葉轉換波形及第二短時距傅立葉轉換波形具有一第一最大振幅及一第二最大振幅，如兩個箭頭所示。
請參照第1及2圖所示，本發明較佳實施例之建構3D偵測模型之心音處理方法包含步驟S4A：接著，利用該CWT資料、STFT資料與一超音波資料進行比對，以尋找一時間與頻率的相關聯性數值狀態，以產生一相關判斷數值。利用該相關判斷數值進行判斷另一量測心音訊號之至少一CWT資料及至少一STFT資料，以產生一判斷結果。本發明較佳實施例之該判斷結果包含心室分隔缺損〔ventricular septal defect，VSD〕或心房分隔缺損〔atrial septal defect，ASD〕。當該心音訊號片段包含一心音雜音訊號時，利用偵測該心音雜音訊號方式進行心臟病理偵測。
第7(A)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統將第一病患之原始心音資料、CWT資料及STFT資料進行比對之波形示意圖。相對的，第7(B)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統將第二病患之原始心音資料、CWT資料及STFT資料進行比對之波形示意圖。請參照第7(A)及7(B)圖之所示，
將第一病患及及第二病患之原始心音資料、CWT資料及STFT資料進行適當對調整〔如第7(A)及7(B)圖之左側及右側所示〕，並選擇分別在該CWT資料之間及STFT資料之間進行波形樣式比對及計算相關性〔如表3所示〕，或選擇在該CWT資料及STFT資料之間進行波形樣式比對及計算相關性。
本發明較佳實施例採用皮爾遜〔Pearson〕積差相關係數〔product-moment coefficient〕計算該CWT資料及STFT資料。本發明較佳實施例採用皮爾遜積差相關係數之計算式如下：
其中z x 及z y 分別為x及y之標準化〔standardized〕z值，r為相關係數。
第2(A)圖揭示本發明較佳實施例之建構3D偵測模型之心音處理方法在轉換資料及建構3D心音模型之流程圖。請參照第1、2及2(A)圖所示，本發明較佳實施例之建構3D偵測模型之心音處理方法可選擇利用數個該CWT資料、數個該STFT資料或其它轉換資料以適當技術手段建構至少一3D心音模型。
請再參照第1、2及2(A)圖所示，本發明另一較佳實施例之建構3D偵測模型之心音處理方法可選擇將數個該CWT資料、數個該STFT資料或其它轉換資料先轉換成數個3D座標資料或其類似各種座標資料，並利用數個該3D座標資料建立該3D心音模型。
請再參照第1及2圖所示，本發明較佳實施例之建構3D偵測模型之心音處理方法包含步驟S4B：接著，利用該CWT資料或STFT資料計算產生至少一特徵值〔eigenvalue〕及至少一特徵向量〔eigenvector〕，且由該特徵值產生一相關數值，或該特徵值與一樣本〔sample〕特徵值〔例如：步驟S4A所示〕進行比對判別，以產生一比對結果，且該樣本特徵值、相關數值或比對結果可選擇儲存於該資料儲存單元30。
舉例而言，本發明另一較佳實施例之該特徵值之相關數值包含一正相關性〔auto-correlation〕數值或一互相關性〔cross-correlation〕數值，且該特徵值由該第二心音S2訊號產生。
本發明另一較佳實施例採用〝正相關性數值〞亦稱為〝序列相關〞，其為一個訊號於其自身在不同時間點的互相關。正相關性數值為兩次觀察之間的相似度對其本身之間的時間差的函數，其用以尋找重複模式〔例如：噪音或雜訊〔noise〕掩蓋的周期訊號〕，或辨識隱藏在訊號諧波頻率中消失的基頻的工具。
本發明另一較佳實施例採用〝互相關性數值〞亦稱為〝互協方差〞。在信號處理領域中互相關用以表示兩個信號之間相似性的度量工具，其通常利用與已知信號比較方式應用於尋找未知信號中的特性，其為兩個信號之間相對於時間的一個函數，有時亦稱為滑動點積。
請再參照第1及2圖所示，舉例而言，本發明較佳實施例在計算該CWT資料或STFT資料之特徵值之正相關性時，可選擇採用公式如下：
，但其並非用以限定本發明。
請再參照第1及2圖所示，舉例而言，本發明較佳實施例在計算該CWT資料或STFT資料之特徵值之互相關性時，可選擇採用公式如下：
，但其並非用以限定本發
明。
舉例而言，本發明另一較佳實施例之該正相關性數值及互相關性數值選擇用以建立一先天心室分隔模型〔congenital septal defect，CSD〕，且該先天心室分隔模型為一相關數值模型〔model〕或其它數學模型。
舉例而言，本發明較佳實施例之該心音訊號片段A與其特徵值之相關數值之正相關性數值及互相關性數值表示如下：
其中，對角線奇異值為正相關性數值，而其餘為互相關性數值；
且
由於多項式為，
因此該心音訊號片段A可表示為，並進一步可表示為|A-λ1|=0，如此，本發明較佳實施例之先天心室分隔模型
表示如下：
，但其並非用以限定本發
明。
請再參照第1、2及2(A)圖所示，本發明較佳實施例之建構3D偵測模型之心音處理方法可選擇利用該3D心音模型進一步計算至少一特徵值及至少一特徵向量，且由該特徵值產生一相關數值，或該特徵值與一樣本特徵值〔例如：步驟S4A所示〕進行比對判別，以產生一比對結果，且該樣本特徵值、相關數值或比對結果可選擇儲存於該資料儲存單元30。
第8圖揭示本發明較佳實施例之建構3D偵測模型之心音處理系統產生一3D心音模型之示意圖。請參照第8圖所示，本發明較佳實施例之建構3D偵測模型之心音處理系統產生至少一3D心音模型，而該3D心音模型包含X軸、Y軸、Z軸，且該3D心音模型提供一3D心音空間，且於該3D心音空間形成兩組3D轉換波形，如第8圖之上半部及下半部所示之兩個箭頭。
請再參照第8圖所示，本發明另一較佳實施例之建構3D偵測模型之心音處理系統可選擇在該3D心音空間之兩組3D轉換波形上進行波形樣式比對及計算相關性，或選擇在該3D心音空間之兩組3D轉換波形與樣本資料之間進行波形樣式比對及計算相關性。
請再參照第8圖所示，本發明另一較佳實施例之建構3D偵測模型之心音處理系統可選擇在該3D心音空間之兩組3D轉換波形上進行一特徵值分析作業或其它類似分析作業。
前述較佳實施例僅舉例說明本發明及其技術特徵，該實施例之技術仍可適當進行各種實質等效修飾及/或替換方式予以實施；因此，本發明之權利範圍須視後附申請專利範圍所界定之範圍為準。
Claims (10)
- 一種建構3D偵測模型之心音處理方法，其包含：將一心音訊號進行切割，以獲得數個心音訊號片段，且該心音訊號包含一第一心音S1、一第二心音S2、一第三心音S3及一第四心音S4；將該心音訊號片段進行連續小波轉換，以獲得一CWT資料；將該心音訊號片段進行短時距傅立葉轉換，以獲得一STFT資料；及利用該CWT資料及STFT資料計算產生至少一3D心音模型，且由該3D心音模型提供一3D心音模型資料。
- 依申請專利範圍第1項所述之建構3D偵測模型之心音處理方法，其中利用該CWT資料、STFT資料或3D心音模型資料計算產生至少一特徵值及至少一特徵向量，且由該特徵值產生一相關數值，或該特徵值與一樣本特徵值進行比對，以產生一比對結果。
- 依申請專利範圍第1項所述之建構3D偵測模型之心音處理方法，其中將數個該CWT資料及數個該STFT資料轉換成數個3D座標資料。
- 依申請專利範圍第1項所述之建構3D偵測模型之心音處理系統，其中該3D心音模型提供一3D心音空間。
- 依申請專利範圍第1項所述之建構3D偵測模型之心音處理系統，其中該3D心音模型顯示一3D轉換波形，以便立體顯示一心室分隔模型、一先天心室分隔模型、一心房分隔模型或一先天心房分隔模型。
- 一種建構3D偵測模型之心音處理系統，其包含：一心音接收單元，其用以接收至少一心音訊號，且該心音訊號包含一第一心音S1、一第二心音S2、一第三心音S3及一第四心音S4；一心音處理單元，其將該心音訊號進行切割，以獲得數 個心音訊號片段，並將該心音訊號片段進行連續小波轉換，以獲得一CWT資料，且將該心音訊號片段進行短時距傅立葉轉換，以獲得一STFT資料；一資料儲存單元，其用以儲存至少一有關心臟病理超音波資料樣本；及一輸出單元，其用以輸出一相關判斷數值；其中利用該CWT資料及STFT資料計算產生至少一3D心音模型，且由該3D心音模型提供一3D心音模型資料。
- 依申請專利範圍第6項所述之建構3D偵測模型之心音處理系統，其中利用該CWT資料、STFT資料或3D心音模型資料計算產生至少一特徵值及至少一特徵向量，且由該特徵值產生一相關數值，或該特徵值與一樣本特徵值進行比對，以產生一比對結果。
- 依申請專利範圍第6項所述之建構3D偵測模型之心音處理系統，其中將數個該CWT資料及數個該STFT資料轉換成數個3D座標資料。
- 依申請專利範圍第6項所述之建構3D偵測模型之心音處理系統，其中該3D心音模型提供一3D心音空間。
- 依申請專利範圍第6項所述之建構3D偵測模型之心音處理系統，其中該3D心音模型顯示一3D轉換波形，以便立體顯示一心室分隔模型、一先天心室分隔模型、一心房分隔模型或一先天心房分隔模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW106143404A TWI667997B (zh) | 2017-12-11 | 2017-12-11 | 建構３ｄ偵測模型之心音處理方法及其系統 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW106143404A TWI667997B (zh) | 2017-12-11 | 2017-12-11 | 建構３ｄ偵測模型之心音處理方法及其系統 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW201927243A TW201927243A (zh) | 2019-07-16 |
| TWI667997B true TWI667997B (zh) | 2019-08-11 |
Family
ID=68048984
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW106143404A TWI667997B (zh) | 2017-12-11 | 2017-12-11 | 建構３ｄ偵測模型之心音處理方法及其系統 |
Country Status (1)
| Country | Link |
|---|---|
| TW (1) | TWI667997B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8684960B2 (en) * | 2009-03-13 | 2014-04-01 | Stanley Batiste | Endothelial scaffold graft and method therefor |
-
2017
- 2017-12-11 TW TW106143404A patent/TWI667997B/zh not_active IP Right Cessation
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8684960B2 (en) * | 2009-03-13 | 2014-04-01 | Stanley Batiste | Endothelial scaffold graft and method therefor |
Also Published As
| Publication number | Publication date |
|---|---|
| TW201927243A (zh) | 2019-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10980429B2 (en) | Method and system for cuffless blood pressure estimation using photoplethysmogram features and pulse transit time | |
| US9968275B2 (en) | Non-invasive method and system for characterizing cardiovascular systems | |
| EP2459065B1 (en) | Method and apparatus for the analysis of a ballistocardiogram signal | |
| US20160045183A1 (en) | Electronic stethoscope apparatus, automatic diagnostic apparatus and method | |
| CN102458259B (zh) | 利用超声识别移动的解剖结构的方法和装置 | |
| Nigam et al. | Accessing heart dynamics to estimate durations of heart sounds | |
| JP6148620B2 (ja) | 経胸郭肺ドップラー超音波を使用した肺血圧測定 | |
| US11534133B2 (en) | Ultrasonic detection method and ultrasonic imaging system for fetal heart | |
| TW201244691A (en) | Heart sound signal/heart disease or cardiopathy distinguishing system and method | |
| JP2015198672A (ja) | 医用画像診断装置 | |
| Spitzer et al. | The role of automated 3D echocardiography for left ventricular ejection fraction assessment | |
| JPH06503483A (ja) | 自動デジタルイメージ定量化システム | |
| CN115804618A (zh) | 人体血管的应力测量装置、方法及存储介质 | |
| WO2022211972A1 (en) | Multisensor pulmonary artery and capillary pressure monitoring system | |
| JP6103591B2 (ja) | 聴診心音信号の処理方法、聴診心音信号の処理装置及び聴診心音信号を処理するためのプログラム | |
| Saeidi et al. | 3D heart sound source localization via combinational subspace methods for long-term heart monitoring | |
| TWI667997B (zh) | 建構３ｄ偵測模型之心音處理方法及其系統 | |
| WO2023086936A1 (en) | Systems and methods for estimating cardiac events | |
| US11730381B1 (en) | Systems, apparatuses, and methods for locating blood flow turbulence in the cardiovascular system | |
| TWI669101B (zh) | 採用特徵值偵測之心音處理方法及其系統 | |
| JP2021500182A (ja) | 冠動脈疾患のリスク判定 | |
| JP4870449B2 (ja) | 超音波診断装置及び超音波画像処理方法 | |
| JP5388520B2 (ja) | 超音波診断装置および超音波診断装置のデータ処理プログラム | |
| Thiyagaraja et al. | Smart phone monitoring of second heart sound split | |
| CN113226174B (zh) | 一种心脏舒张功能评估方法、设备和系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | Annulment or lapse of patent due to non-payment of fees |