重症监护病人心电导联信号质量评估、SQI

重症监护病人心电导联信号质量评估、SQI

上一篇内容总结了第三方库ecg-qc中计算的SQI 类型,本篇是在搜集资料过程中发现的一篇有关SQI的文章,在此记录一下本篇文章的SQI计算方法
文章:重症监护病人心电导联信号质量评估

1、bSQI

bSQI是用两种QRS检测算法,对同一个ECG信号进行QRS检测。如果两个算法检测出来的QRS波相符合,则信号质量较好,即bSQI数值较高。

b S Q I = N m a t c h e d ( k , w ) N a l l ( k , w ) bSQI = \frac{N_{matched}\left ( k,w\right)}{N_{all}\left ( k,w\right)} bSQI=Nall(k,w)Nmatched(k,w)
其中k 为当前的QRS 波‚w 为滑动分析窗口宽度10s‚以当前 QRS 波为中心,左右各5s, N m a t c h e d ( k , w ) N_{matched}\left ( k,w\right) Nmatched(k,w)为在 w 中两种QRS检测算法(DF、LT算法)检测出的 QRS 波匹配的数目, N a l l ( k , w ) N_{all}\left ( k,w\right) Nall(k,w)为在 w 中两种算法各自检测出的 QRS 波的数目总和。

2、iSQI

iSQI是用一种QRS检测算法,对同一时间同一病人的不同导联的ECG进行QRS检测。如果在同一时间不同导联检测的QRS相对应,即至少两导联信号都比较高时,iSQI数值较大。
i S Q I = m a x ( N m a t c h e d ( k , w ) i , j N a l l ( k , w ) i , j ) iSQI =max \left (\frac{N_{matched}\left ( k,w\right)_{i,j}}{N_{all}\left ( k,w\right)_{i,j}} \right ) iSQI=max(Nall(k,w)i,jNmatched(k,w)i,j)
其中,i为当前导联,j为不同于 i 的其它同步导联。 N m a t c h e d ( k , w ) N_{matched}\left ( k,w\right) Nmatched(k,w) N a l l ( k , w ) N_{all}\left ( k,w\right) Nall(k,w)为应用 同一QRS检测算法对不同导联分别进行QRS 检测时的匹配计数和总计数。

3、kSQI

kSQI是衡量信号的分布是否平均的指标。信号 x 的四阶矩( ν 4 \nu _{4} ν4)或称峰度(kurtosis,K)是信号高斯性的量度。一般来说,高斯分布的 K值为3;正常窦性心律的 K 值大于5;肌噪干扰的 K 值为5左右;基线漂移和50Hz 工频干扰的 K 值均小于5。
定义式为: k S Q I ( k ) = { 1 i f   k u r t o s i s ( k , w ) > 5 0 i f   k u r t o s i s ( k , w ) ≤ 5 kSQI(k) = \left\{\begin{matrix} 1 & if \ kurtosis\left ( k,w \right ) >5\\ 0 & if \ kurtosis\left ( k,w \right ) \leq 5 \end{matrix}\right. kSQI(k)={10if kurtosis(k,w)>5if kurtosis(k,w)5
当前波的 k u r t o s i s ( k , w ) 是 当 前 滑 动 窗 口 的 四 阶 矩 : kurtosis\left ( k,w \right )是当前滑动窗口的四阶矩: kurtosis(k,w) k _ s q i = E { ( x − μ x ) 4 } σ 4 k\_sqi=\frac{E\{(x-\mu_x)^4\}}{\sigma^4} k_sqi=σ4E{(xμx)4}

4、sSQI

定义式为:
s S Q I = { 1 i f    S D R ≥ 0.5    a n d    S D R   ≤ 0.8 0 i f    S D R < 0.5    o r    S D R   > 0.8 sSQI = \left\{\begin{matrix} 1 & if \ \ SDR\geq0.5 \ \ and\ \ SDR\ \leq0.8 \\ 0 & if \ \ SDR< 0.5 \ \ or\ \ SDR\ >0.8 \end{matrix}\right. sSQI={10if  SDR0.5  and  SDR 0.8if  SDR<0.5  or  SDR >0.8
其中SDR:
S D R ( k ) = ∑ f = 5 14 P S D ( k , w ) ∑ f = 5 50 P S D ( k , w ) SDR\left ( k \right ) = \frac{\sum_{f=5}^{14}PSD\left ( k,w \right )}{\sum_{f=5}^{50}PSD\left ( k,w \right )} SDR(k)=f=550PSD(k,w)f=514PSD(k,w)

5、ECG信号综合质量评估

E C G S Q I = { m a x ( b S Q I ( k ) , i S Q I ( k ) ) i f   k S Q I ( k ) = 1 &   s S Q I ( k ) = 1 b S Q I ( k ) i f   k S Q I ( k ) = 1 &   s S Q I ( k ) = 0 m a x ( b S Q I ( k ) , i S Q I ( k ) ) ∗ η i f   k S Q I ( k ) = 0 &   s S Q I ( k ) = 1 b S Q I ( k ) ∗ η i f   k S Q I ( k ) = 0 &   s S Q I ( k ) = 0 ECGSQI = \left\{\begin{matrix} max\left ( bSQI\left ( k \right ),iSQI\left ( k \right ) \right ) & if \ kSQI\left ( k \right )=1 & \& \ sSQI\left ( k \right )=1\\ bSQI\left ( k \right ) & if \ kSQI\left ( k \right )=1 & \& \ sSQI\left ( k \right )=0\\ max\left ( bSQI\left ( k \right ),iSQI\left ( k \right ) \right )*\eta & if \ kSQI\left ( k \right )=0 & \& \ sSQI\left ( k \right )=1\\ bSQI\left ( k \right )*\eta & if \ kSQI\left ( k \right )=0 & \& \ sSQI\left ( k \right )=0 \end{matrix}\right. ECGSQI=max(bSQI(k),iSQI(k))bSQI(k)max(bSQI(k),iSQI(k))ηbSQI(k)ηif kSQI(k)=1if kSQI(k)=1if kSQI(k)=0if kSQI(k)=0& sSQI(k)=1& sSQI(k)=0& sSQI(k)=1& sSQI(k)=0
其中,0< η \eta η<1是通过实验调整的经验系数,文章中取η=0.7

  • 2
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值