[发明专利]一种微伏级T波交替的定量检测方法

权利要求书

1.一种微伏级T波交替的定量检测方法，其特征在于，包括以下步骤：

（1）获取心电图数据并进行心电信号预处理；

（2）心电信号特征点标定：对心电图中的QRS波群及T波进行标定；

（3）心拍选择：采取多个步骤的复合预处理对系列心拍进行选择；

（4）选取T波采样点：采用T波窗口，并以T波波峰对齐，选取128个T波段，每段T波提取7个采样点，构成待测的T波采样点序列；

（5）作T波交替散点图：对步骤（4）中的T波采样点序列一次差分后作散点图，并利用MATLAB7.0软件仿真绘图；

（6）确定散点图“边缘”，采取“横向搜索法”对散点图进行边缘提取，确定散点图的“边缘”并得到散点图边缘点集；

（7）提取散点图的定量参数，分析散点图外部形态与T波交替的关系，根据散点图“边缘”的几何特征，定义短轴、长轴、短长轴之比三个定量参数；

（8）确定T波交替的散点图定量检测指标，以步骤（7）所述短长轴之比作为定量检测指标，采用MATLAB7.0软件仿真并与谱分析法比较，找出合适阈值判断T波交替的存在与否。

2.根据权利要求1所述的检测方法，其特征在于，步骤（4）中，所述复合预处理包括心拍剔除。

3.根据权利要求1所述的检测方法，其特征在于，步骤（6）中所述采用“横向搜索法”对散点图进行边缘提取，具体步骤如下：

61）散点图预处理：将散点图转换成M×N的二值化的数字图像矩阵G，根据其横纵坐标序列X={x₂-x₁,x₃-x₂,…,x_i+1-x_i,…}、Y={x₃-x₂,x₄-x₃,…,x_i+2-x_i+1,…}，其中X和Y的长度都为L，及散点坐标(x_i+1-x_i,x_i+2-x_i+1)，首先将散点图坐标序列值取整后上移和右移，使其刚好完全移至平面直角坐标系第一象限内，然后通过如下方式进行转换：

G=g11g12···g1Ng21g22···g2N·········gM1gM2···gMN,M=max[Y(n)],N=max[X(n)],n=1:L,]]>

其中Y(n)，X(n)分别表示Y，X中的第n个元素值；max[Y(n)]，max[X(n)]分别表示X，Y中所有元素的最大值；矩阵元素g_i,j表示图像G中第i行，第j列的像素灰度值，其定义为：

得散点图的二值化数字图像G；

62）寻找边缘候选点：对图像G进行块操作，初步去除散点图中心非边缘点及边界外较稀疏的非边缘点，定义5×5的整数格，令整数格以g_i，j，其中，i为1~M之间的整数，j为1~N之间的整数，为中心像素开始对散点图进行行扫描，计算该整数格覆盖区域

的灰度值和a，a是A中各元素之和，若2≤a≤3，则将该中心像素标记为边缘候选点；

63）确定边缘点：对由上个步骤确定的边缘候选点进行筛选，去除散点图中心区域外的一些异常小点集，得到最终边缘点；首先以边缘候选点为中心像素构建21×21的整数格，则该整数格覆盖区域

将矩阵B按

划分，得矩阵B1、B2、B3、B4，分别计算其灰度值和b1、b2、b3、b4，b1、b2、b3、b4是相应矩阵B1、B2、B3、B4中各元素之和，若b1、b2、b3、b4的值都小于阈值6，则删除此候选点，将其余点保存在边缘点数组中，得散点图的边缘点集。