一種多通道心電圖的qrs自動劃定方法
【專利摘要】本發(fā)明公開了一種多通道心電圖的QRS自動劃定方法，其采用小波變換生成每個通道的實(shí)時能量曲線，對每個通道的能量曲線計(jì)算其峰度和其他信號質(zhì)量的度量，對能量曲線取其絕對值，根據(jù)信號質(zhì)量的度量對多個通道取絕對值后的能量曲線計(jì)算其加權(quán)平均，該平均后的能量曲線作為綜合的心電能量曲線，并對其進(jìn)行劃定。由此本發(fā)明避免了由于個別通道信號劃定不同步而產(chǎn)生R峰檢測出錯的問題；同時在信號擬合的過程中參考了信號的峰度值，避免了由于某通道信號波形質(zhì)量太差而造成整體綜合能量曲線質(zhì)量下降的問題。因此本發(fā)明方法檢測得到的心電R峰效果穩(wěn)定準(zhǔn)確，特別是當(dāng)心電信號中部分通道信號質(zhì)量不過關(guān)時，也能準(zhǔn)確的檢測到R峰的位置。
【專利說明】
-種多通道心電圖的QRS自動劃定方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于屯、電檢測技術(shù)領(lǐng)域，具體設(shè)及一種多通道屯、電圖的QRS自動劃定方法。
【背景技術(shù)】
[0002] 屯、電圖化lectrocardiogram，ECG)是在屯、臟有規(guī)律的收縮和舒張的過程中，各部 分屯、肌細(xì)胞產(chǎn)生的動作電位綜合而成的電信號由電極從體表或胸腔測得，經(jīng)放大后顯示或 描記下來的波形。屯、電圖通常由一系列相同的波群構(gòu)成，主要包括W下：
[0003] P波：代表左右屯、房去極化過程的電位變化，歷時0.08~0111s，波幅不超過 0.25mv〇
[0004] QRS波群：代表左右屯、室去極化過程的電位變化，歷時0.06~0.10s;典型的QRS波 群包含Ξ個緊密相聯(lián)的波，第一個向下的波為Q波，其后向上的高而減的波為R波，繼R波之 后的一個向下的波為S波。
[000引PR間期：由P波起點(diǎn)至化RS波群起點(diǎn)之間的間期，代表自屯、房除極開始至屯、室除極 的時間，正常成人為0.12~0.20s。
[0006] T波:代表屯喧快速復(fù)極過程中的電位變化，波形圓純，歷時0.05~0.25s，波形的 前肢較長而后肢較短。
[0007] ST段：自QRS波終點(diǎn)至T波起點(diǎn)之間的線段，代表屯、室緩慢復(fù)極過程，正常的ST段多 為一等電位線，有時也可有輕微的偏移。在任一導(dǎo)聯(lián)，ST段下移一般不超過0.05mv;ST段上 抬在V1、V2導(dǎo)聯(lián)不超過0.3mv，在V3不超過0.5mv，在V4~V6導(dǎo)聯(lián)與肢體導(dǎo)聯(lián)不應(yīng)超過0.1 mv。 [000引在屯噸信號自動分析中，QRS波的檢測尤為重要。QRS波是ECG中最顯著的部分，它 反映了屯、室收縮時屯、臟的電行為，其能量占了屯、電信號能量的很大百分比。QRS波分部于屯、 電信號中、高頻區(qū)域，峰值落在10~20監(jiān)之間，并且幅度特征非常明顯，與其他波形區(qū)別顯 著，因而在屯、電特征波形的檢測中，常常先定位QRS波群，當(dāng)QRS波被確定后，才有可能計(jì)算 屯、率并進(jìn)行屯、率變異分析，才能檢測ST段的參數(shù)和分析屯、電的其他細(xì)節(jié)信息，由此QRS波的 檢測是ECG檢測中的首要問題，是診斷屯、率失常的重要依據(jù)。
[0009] 學(xué)者化η于1985年提出了一種雙闊值檢測方法，其首先讓屯、電信號一次經(jīng)過一個 帶通濾波和低通濾波的差分器，W增強(qiáng)5~15化的信號分量;接著，對信號進(jìn)行逐點(diǎn)平方和 移動窗口（150ms)積分運(yùn)算，得到變換后的信號。若差分后的信號和積分哈偶的信號同事滿 足闊值條件，則認(rèn)為檢測到一個QRS波，更新闊值。學(xué)者化ekelmans也提出過一種QRS檢測方 法，其使屯、電信號經(jīng)過500HZ的頻率采樣后，首先通過一個10點(diǎn)平均濾波器，W濾除50化工 頻干擾;將信號與一個經(jīng)過8Hz低通濾波的信號相減，形成帶通濾波。將各通道濾波后的各 對應(yīng)數(shù)據(jù)點(diǎn)的絕對值相加，得到信號A;將各通道信號的一階差分的絕對值相加得到信號D; 將信號A與信號D相加，W去除信號A中可能出現(xiàn)的凹陷；在離線處理部分，首先對數(shù)據(jù)進(jìn)行 統(tǒng)計(jì)，取得幅度闊值和斜率闊值;若幅度和前后斜率均超過闊值，則認(rèn)為檢測到一個QRS波。
[0010] 上述化η提出的算法適用于單通道屯、電圖，但多通道屯、電圖常常具有兩個問題： (1)每個通道的信號質(zhì)量不同；（2)各個通道的劃定不同步。如果隨機(jī)采取一個通道作為屯、 電圖的劃定的依據(jù)，則結(jié)果可能不準(zhǔn)確。Brekelmans雖然綜合各通道信號進(jìn)行判斷，但僅僅 是將各通道的信號簡單的疊加，若第一通道的導(dǎo)聯(lián)安放位置不好，導(dǎo)致信號低而噪聲大，而 第二個通道的信號質(zhì)量優(yōu)良，則對后一個通道來說，疊加后的數(shù)據(jù)信噪比反而降低。

【發(fā)明內(nèi)容】

[0011]針對現(xiàn)有技術(shù)所存在的上述技術(shù)缺陷，本發(fā)明提供了一種多通道屯、電圖的QRS自 動劃定方法，能夠避免通道間信號質(zhì)量差距所引起的檢測結(jié)果差異，提高R峰檢測的成功 率。
[0012] -種多通道屯噸圖的QRS自動劃定方法，包括如下步驟：
[0013] (1)對各通道的屯、電信號進(jìn)行小波變換得到各通道的能量曲線；
[0014] (2)對各通道的能量曲線進(jìn)行一階差分變換，得到各通道的能量差分曲線；
[0015] (3)計(jì)算各通道能量差分曲線的峰度，進(jìn)而確定各通道的權(quán)值；
[0016] (4)根據(jù)權(quán)值將各通道的能量差分曲線進(jìn)行累加得到綜合能量曲線，進(jìn)而采用窗 函數(shù)對綜合能量曲線進(jìn)行濾波；
[0017] (5)對濾波后的綜合能量曲線進(jìn)行R峰檢測，即完成對多通道屯噸圖的QRS劃定。
[0018] 所述的步驟(1)中根據(jù)W下公式對各通道的屯、電信號進(jìn)行小波變換：
[0019]
[0020] 其中：x(t)為任一通道的屯、電信號，z(t)為x(t)經(jīng)小波變換后得到的能量曲線，a 為縮放因子，b為平移參數(shù)，t為時間，餐^ )為墨西哥帽小波函數(shù)φ()的復(fù)數(shù)共輛。
[0021] 所述的步驟(2)中根據(jù)W下公式對各通道的能量曲線進(jìn)行一階差分變換：
[0022] y[n] =ζ[η]-ζ[η-1]
[0023] 其中：ζ[η]和ζ[η-1]分別為任一通道能量曲線中第η個采樣點(diǎn)和第η-1個采樣點(diǎn)的 信號值，y[n]為對應(yīng)通道能量差分曲線中第η個采樣點(diǎn)的信號值，η為大于0的自然數(shù)。
[0024] 所述的步驟(3)中根據(jù)W下公式計(jì)算各通道能量差分曲線的峰度：
[0025]
[0026] 其中：Κ為任一通道能量差分曲線的峰度，y[n]為對應(yīng)通道能量差分曲線中第η個 采樣點(diǎn)的信號值，歹為對應(yīng)通道能量差分曲線中所有采樣點(diǎn)的信號平均值，Ν為對應(yīng)通道能 量差分曲線中所有采樣點(diǎn)的個數(shù)。
[0027] 所述的步驟(3)中根據(jù)W下公式確定各通道的權(quán)值：
[002引
[0029] 其中:W功第i通道的權(quán)值，Κ功第i通道能量差分曲線的峰度，m為通道總數(shù)。
[0030] 所述的步驟(4)中根據(jù)W下公式將各通道的能量差分曲線進(jìn)行累加得到綜合能量 曲線：
[0031]
[0032] 其中：X為綜合能量曲線，wi為第i通道的權(quán)值，yi為第i通道的能量差分曲線，m為通 道總數(shù)。
[0033] 所述的步驟(4)中根據(jù)W下公式對綜合能量曲線進(jìn)行濾波：
[0034]
[0035] 其中：X[j-p]為綜合能量曲線中第j-p個采樣點(diǎn)的信號值，Xlj]為濾波后的綜合 能量曲線中第j個采樣點(diǎn)的信號值，Μ為窗口大小，j為大于0的自然數(shù)。
[0036] 所述的步驟(5)中對濾波后的綜合能量曲線進(jìn)行R峰檢測，具體過程如下：
[0037] 首先，提取濾波后的綜合能量曲線中前k秒每秒內(nèi)所有采樣點(diǎn)最大的信號值，共得 至化個信號值組成峰值序列；同時構(gòu)建同樣由k個信號值組成的噪聲序列W及由k個RR間期 組成的間期序列，初始化噪聲序列中的k個信號值均為0,初始化間期序列中的k個RR間期均 為1000ms，k為大于1的自然數(shù);進(jìn)而根據(jù)W下公式計(jì)算闊值t虹esh:
[0038] t 虹 esh= (medianI?p-medianNoise)*0
[0039] 其中:median化和medianNoise分別為峰值序列和噪聲序列的中值，β為預(yù)設(shè)系數(shù)；
[0040] 然后，逐個判斷濾波后的綜合能量曲線中每個采樣點(diǎn)：若當(dāng)前采樣點(diǎn)的信號值大 于其之后200ms內(nèi)所有采樣點(diǎn)的信號值，則令當(dāng)前采樣點(diǎn)為候選R峰；
[0041 ] 若候選R峰的信號值> thresh，則判定候選R峰為真正R峰，進(jìn)而移除峰值序列中對 應(yīng)時間最早的一個信號值且將該真正R峰的信號值填補(bǔ)至峰值序列中，同時移除間期序列 中對應(yīng)時間最早的一個RR間期且將該真正R峰與前一真正R峰的時間間隔作為新的RR間期 填補(bǔ)至間期序列中，并更新thresh;
[0042] 若候選R峰的信號值<thresh，則判定候選R峰為噪聲點(diǎn)，進(jìn)而移除噪聲序列中對 應(yīng)時間最早的一個信號值且將該噪聲點(diǎn)的信號值填補(bǔ)至噪聲序列中，并更新t虹esh;
[0043] 若超過1.5medianRR的時間段仍未找到下一個真正R峰，medianRR為間期序列的中 值，則判斷該時間段內(nèi)信號值最大的噪聲點(diǎn)，其信號值是否大于〇.5thresh:若是，則判定該 噪聲點(diǎn)為真正R峰，進(jìn)而移除峰值序列中對應(yīng)時間最早的一個信號值且將該真正R峰的信號 值填補(bǔ)至峰值序列中，同時移除間期序列中對應(yīng)時間最早的一個RR間期且將該真正R峰與 前一真正R峰的時間間隔作為新的RR間期填補(bǔ)至間期序列中，并更新t虹esh。
[0044] 本發(fā)明通過多通道屯、電信號擬合，避免了由于個別通道信號劃定不同步而產(chǎn)生R 峰檢測出錯的問題；同時在信號擬合的過程中參考了信號的峰度值，避免了由于某通道信 號波形質(zhì)量太差而造成整體綜合能量曲線質(zhì)量下降的問題。因此本發(fā)明方法檢測得到的屯、 電R峰效果穩(wěn)定準(zhǔn)確，特別是當(dāng)屯、電信號中部分通道信號質(zhì)量不過關(guān)時，也能準(zhǔn)確的檢測到 R峰的位置。
【附圖說明】
[0045] 圖1為本發(fā)明方法的步驟流程示意圖。
[0046] 圖2為Ξ通道屯、電信號及通過本發(fā)明加權(quán)擬合后檢測標(biāo)注有R峰屯、電信號的波形 示意圖。
【具體實(shí)施方式】
[0047] 為了更為具體地描述本發(fā)明，下面結(jié)合附圖及【具體實(shí)施方式】對本發(fā)明的技術(shù)方案 進(jìn)行詳細(xì)說明。
[0048] 如圖1所示，本發(fā)明多通道屯、電圖的QRS自動劃定方法，包括如下步驟：
[0049] (1)小波變換。
[0050] 屯、電圖的一個通道可W表示為一個時間序列函數(shù)x(t)。對該函數(shù)的小波變換的定 義如下：
[0化1 ]
[(K)對上式中：母小波柳t)是時頻域上的連續(xù)函數(shù)，餐*代表函數(shù)巧的復(fù)數(shù)共輛，a是縮放 因子，b是平移參數(shù)。在本實(shí)施方式中采用墨西哥帽小波來作為母小波用作屯、電的分析。小 波分析后，獲取20.細(xì)Z- 2甜Z頻段之間的能量曲線Z (t)。
[0053] (2) -階差分變換。
[0054] 對各通道小波能量曲線進(jìn)行一階差分變換：
[0055] y[n] =z[n]-z[n-l]
[0056] 其中：z[n]和z[n-l]分別為任一通道能量曲線中第η個采樣點(diǎn)和第n-1個采樣點(diǎn)的 信號值，y[n]為對應(yīng)通道能量差分曲線中第η個采樣點(diǎn)的信號值，η為大于0的自然數(shù)。
[0057] (3)計(jì)算峰度并確定權(quán)值。
[005引本實(shí)施方式用信號的峰度來代表每個通道的信號質(zhì)量，峰度又稱為峰態(tài)系數(shù)，在 統(tǒng)計(jì)學(xué)中亦成為四階標(biāo)準(zhǔn)矩。峰度用來表征分布曲線在平均值處峰值高低。用y[n]表示任 一通道能量差分曲線中第η個采樣點(diǎn)的信號值，用歹來表示能量差分曲線的信號均值，則峰 度Κ定義如下：
[0化9]
[0060] 根據(jù)上述得到每個通道的峰度值獲取各通道的權(quán)值：
[0061]
[0062] 其中:W功第i通道的權(quán)值，Κ功第i通道能量差分曲線的峰度，m為通道總數(shù)。
[0063] (4)絕對值加權(quán)累加。
[0064] 取各通道能量差分曲線的絕對值，根據(jù)步驟(3)中的權(quán)值進(jìn)行累加，得到綜合能量 曲線X:
[00 化]
[0066] 其中：yi為第i通道的能量差分曲線。
[0067] (5)窗口 平均。
[0068] 對取得的綜合能量曲線進(jìn)行80ms窗口平均，消除曲線上的小毛刺，便于檢測波峰；
[0069]
[0070] 其中：X[j-p]為綜合能量曲線中第j-p個采樣點(diǎn)的信號值，Xlj]為濾波后的綜合 能量曲線中第j個采樣點(diǎn)的信號值，j為大于0的自然數(shù)。
[0071] (6)初始化闊值。
[0072] 提取濾波后的綜合能量曲線中前8秒每秒中所有采樣點(diǎn)最大的信號值，共得到8個 信號值組成峰值序列；同時構(gòu)建同樣由8個信號值組成的噪聲序列W及由8個RR間期組成的 間期序列，初始化噪聲序列中的8個信號值均為0,初始化間期序列中的8個RR間期均為 1000ms;進(jìn)而根據(jù)W下公式計(jì)算闊值t虹esh:
[0073] t 虹 esh= (medianI?p-medianNoise)*0.45
[0074] 其中:median化和medianNoise分別為峰值序列和噪聲序列的中值。
[00巧](7)R峰檢測。
[0076] 逐個判斷濾波后的綜合能量曲線中每個采樣點(diǎn):若當(dāng)前采樣點(diǎn)的信號值大于其之 后200ms內(nèi)所有采樣點(diǎn)的信號值，則令當(dāng)前采樣點(diǎn)為候選R峰；
[0077] 若候選R峰的信號值 > 闊值t虹esh，則判定候選R峰為真正R峰，進(jìn)而移除峰值序列 中對應(yīng)時間最早的一個信號值且將該真正R峰的信號值填補(bǔ)至峰值序列中，同時移除間期 序列中對應(yīng)時間最早的一個RR間期且將該真正R峰與前一真正R峰的時間間隔作為新的RR 間期填補(bǔ)至間期序列中，并更新闊值t虹esh;
[0078] 若候選R峰的信號值<闊值thresh,則判定候選R峰為噪聲點(diǎn)，進(jìn)而移除噪聲序列 中對應(yīng)時間最早的一個信號值且將該噪聲點(diǎn)的信號值填補(bǔ)至噪聲序列中，并更新闊值 thresh；
[0079] 若超過1.5medianRR的時間段仍未找到下一個真正R峰，medianRR為間期序列的中 值，則判斷該時間段內(nèi)信號值最大的噪聲點(diǎn)，其信號值是否大于〇.5thresh:若是，則判定該 噪聲點(diǎn)為真正R峰，進(jìn)而移除峰值序列中對應(yīng)時間最早的一個信號值且將該真正R峰的信號 值填補(bǔ)至峰值序列中，同時移除間期序列中對應(yīng)時間最早的一個RR間期且將該真正R峰與 前一真正R峰的時間間隔作為新的RR間期填補(bǔ)至間期序列中，并更新t虹esh。
[0080] 圖2所示了一組Ξ導(dǎo)聯(lián)屯、電圖，如果隨意選取其中一個通道來做屯、率的劃定，則可 能得到不準(zhǔn)確的結(jié)果。比如，通道二和通道Ξ很容易給出錯誤的屯、率。在運(yùn)里，本實(shí)施方式 采用Ξ個導(dǎo)聯(lián)的能量基于峰度的加權(quán)，得到一個擬合后的加權(quán)屯、電能量信號，其中小圓圈 為采用加權(quán)屯、電能量信號標(biāo)注的R峰位置，可W看到加權(quán)的能量可W給出正確的屯、率。
[0081] 使傳統(tǒng)的自適應(yīng)差分闊值法對Mit-bih屯、律失常數(shù)據(jù)庫的檢測結(jié)果與本發(fā)明多通 道自動劃定算法結(jié)果進(jìn)行對比（共8個數(shù)據(jù)庫），從表1觀察檢測準(zhǔn)確率統(tǒng)計(jì)數(shù)據(jù)可W發(fā)現(xiàn)， 由于通道一和通道二數(shù)據(jù)的信號質(zhì)量存在差距，導(dǎo)致檢測結(jié)果差距十分明顯，而經(jīng)過多通 道加權(quán)平均后的劃定算法對每個庫的數(shù)據(jù)檢測的準(zhǔn)確率都很高比較穩(wěn)定。
[0082] 表 1
[0083]

[0084]~上述的對實(shí)施例的描述是為便于本技術(shù)領(lǐng)域的普通技術(shù)人員能理解和應(yīng)用本發(fā) 明。熟悉本領(lǐng)域技術(shù)的人員顯然可W容易地對上述實(shí)施例做出各種修改，并把在此說明的 一般原理應(yīng)用到其他實(shí)施例中而不必經(jīng)過創(chuàng)造性的勞動。因此，本發(fā)明不限于上述實(shí)施例， 本領(lǐng)域技術(shù)人員根據(jù)本發(fā)明的掲示，對于本發(fā)明做出的改進(jìn)和修改都應(yīng)該在本發(fā)明的保護(hù) 范圍之內(nèi)。
【主權(quán)項(xiàng)】
1. 一種多通道心電圖的QRS自動劃定方法，包括如下步驟： (1) 對各通道的心電信號進(jìn)行小波變換得到各通道的能量曲線； (2) 對各通道的能量曲線進(jìn)行一階差分變換，得到各通道的能量差分曲線； (3) 計(jì)算各通道能量差分曲線的峰度，進(jìn)而確定各通道的權(quán)值； (4) 根據(jù)權(quán)值將各通道的能量差分曲線進(jìn)行累加得到綜合能量曲線，進(jìn)而采用窗函數(shù) 對綜合能量曲線進(jìn)行濾波； (5) 對濾波后的綜合能量曲線進(jìn)行R峰檢測，即完成對多通道心電圖的QRS劃定。2. 根據(jù)權(quán)利要求1所述的QRS自動劃定方法，其特征在于:所述的步驟（1)中根據(jù)以下公 式對各通道的心電信號進(jìn)行小波變換：其中：x(t)為任一通道的心電信號，z(t)為x(t)經(jīng)小波變換后得到的能量曲線，a為縮 放因子，b為平移參數(shù)，t為時間，識1 )為墨西哥帽小波函數(shù)妒（）的復(fù)數(shù)共輒。3. 根據(jù)權(quán)利要求1所述的QRS自動劃定方法，其特征在于:所述的步驟(2)中根據(jù)以下公 式對各通道的能量曲線進(jìn)行一階差分變換： y[n] = z[n]-z[n-l] 其中：z[n]和z[n-l]分別為任一通道能量曲線中第η個采樣點(diǎn)和第n-1個采樣點(diǎn)的信號 值，y[n]為對應(yīng)通道能量差分曲線中第η個采樣點(diǎn)的信號值，η為大于〇的自然數(shù)。4. 根據(jù)權(quán)利要求1所述的QRS自動劃定方法，其特征在于:所述的步驟(3)中根據(jù)以下公 式計(jì)算各通道能量差分曲線的峰度：其中：Κ為任一通道能量差分曲線的峰度，y[n]為對應(yīng)通道能量差分曲線中第η個采樣 點(diǎn)的信號值，羅為對應(yīng)通道能量差分曲線中所有采樣點(diǎn)的信號平均值，Ν為對應(yīng)通道能量差 分曲線中所有采樣點(diǎn)的個數(shù)。5. 根據(jù)權(quán)利要求1所述的QRS自動劃定方法，其特征在于:所述的步驟(3)中根據(jù)以下公 式確定各通道的權(quán)值：其中:Wl為第i通道的權(quán)值，L為第i通道能量差分曲線的峰度，m為通道總數(shù)。6. 根據(jù)權(quán)利要求1所述的QRS自動劃定方法，其特征在于:所述的步驟(4)中根據(jù)以下公 式將各通道的能量差分曲線進(jìn)行累加得到綜合能量曲線：其中:X為綜合能量曲線，Wl為第i通道的權(quán)值，yi為第i通道的能量差分曲線，m為通道總 數(shù)。7. 根據(jù)權(quán)利要求1所述的QRS自動劃定方法，其特征在于:所述的步驟(4)中根據(jù)以下公 式對綜合能量曲線進(jìn)行濾波：其中:X[j_P]為綜合能量曲線中第j-ρ個采樣點(diǎn)的信號值，j]為濾波后的綜合能量曲 線中第j個采樣點(diǎn)的信號值，Μ為窗口大小，j為大于0的自然數(shù)。8.根據(jù)權(quán)利要求1所述的QRS自動劃定方法，其特征在于:所述的步驟(5)中對濾波后的 綜合能量曲線進(jìn)行R峰檢測，具體過程如下： 首先，提取濾波后的綜合能量曲線中前k秒每秒內(nèi)所有采樣點(diǎn)最大的信號值，共得到k 個信號值組成峰值序列；同時構(gòu)建同樣由k個信號值組成的噪聲序列以及由k個RR間期組成 的間期序列，初始化噪聲序列中的k個信號值均為0，初始化間期序列中的k個RR間期均為 1000ms，k為大于1的自然數(shù);進(jìn)而根據(jù)以下公式計(jì)算閾值thresh: thresh=(medianRp_medianNoise)*0 其中:medianRp和medianNoise分別為峰值序列和噪聲序列的中值，β為預(yù)設(shè)系數(shù)； 然后，逐個判斷濾波后的綜合能量曲線中每個采樣點(diǎn)：若當(dāng)前采樣點(diǎn)的信號值大于其 之后200ms內(nèi)所有采樣點(diǎn)的信號值，則令當(dāng)前采樣點(diǎn)為候選R峰； 若候選R峰的信號值多thresh，貝lj判定候選R峰為真正R峰，進(jìn)而移除峰值序列中對應(yīng)時 間最早的一個信號值且將該真正R峰的信號值填補(bǔ)至峰值序列中，同時移除間期序列中對 應(yīng)時間最早的一個RR間期且將該真正R峰與前一真正R峰的時間間隔作為新的RR間期填補(bǔ) 至間期序列中，并更新thresh; 若候選R峰的信號值< thre sh，則判定候選R峰為噪聲點(diǎn)，進(jìn)而移除噪聲序列中對應(yīng)時 間最早的一個信號值且將該噪聲點(diǎn)的信號值填補(bǔ)至噪聲序列中，并更新thresh; 若超過1.5medianRR的時間段仍未找到下一個真正R峰，medianRR為間期序列的中值， 則判斷該時間段內(nèi)信號值最大的噪聲點(diǎn)，其信號值是否大于〇.5threSh:若是，則判定該噪 聲點(diǎn)為真正R峰，進(jìn)而移除峰值序列中對應(yīng)時間最早的一個信號值且將該真正R峰的信號值 填補(bǔ)至峰值序列中，同時移除間期序列中對應(yīng)時間最早的一個RR間期且將該真正R峰與前 一真正R峰的時間間隔作為新的RR間期填補(bǔ)至間期序列中，并更新thresh。
【文檔編號】A61B5/0456GK105919590SQ201610388877
【公開日】2016年9月7日
【申請日】2016年6月2日
【發(fā)明人】趙曉鵬, 姚劍, 何挺挺, 姚志邦
【申請人】浙江銘眾科技有限公司, 浙江銘眾醫(yī)療器械有限公司, 浙江銘眾生物醫(yī)用材料與器械研究院