心磁图生成方法及生成系统.pdf

《心磁图生成方法及生成系统.pdf》由会员分享，可在线阅读，更多相关《心磁图生成方法及生成系统.pdf（13页完整版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201710564362.7 (22)申请日 2017.07.12 (71)申请人 中国科学院上海微系统与信息技术 研究所 地址 200050 上海市长宁区长宁路865号 (72)发明人 孔祥燕 鲁丽 陈威 杨瑞虎 王佳磊 (74)专利代理机构 上海光华专利事务所(普通 合伙) 31219 代理人 徐秋平 (51)Int.Cl. A61B 5/04(2006.01) A61B 5/0452(2006.01) (54)发明名称 心磁图生成方法及生成系统 (57)摘要 本发明提供一种。

2、应用于心磁图仪的心磁图 生成方法及系统， 包括： 1） 首先截取各个通道一 个周期的数据， 并从截取的数据中选取某一时 刻， 取出各个通道该时刻的数据， 然后计算各通 道在胸部上方对应的位置坐标； 2） 让各通道某个 时刻的数据与坐标对应； 3） 对按位置排列好的二 维心磁数据进行二维插值； 4） 对插值后的二维数 据查询颜色表； 5） 画心磁图。 本发明生成的心磁 图， 图像清晰， 颜色分布均匀， 在T波段可清晰的 看到+-两个偶极子。 此外， 若从QRS波的S段到T段 按设定步长画不同时刻的心磁图， 则可以看到心 脏活动电流的时空演化过程， 医生可根据心磁图 的偶极子数目和心脏活动电流的演。

3、化方向诊断 该病人是否患心脏疾病。 权利要求书3页 说明书6页 附图3页 CN 107374610 A 2017.11.24 CN 107374610 A 1.一种心磁图的生成方法， 其特征在于， 包括： 从心磁图仪的各通道中分别截取一个周期的心磁数据， 从截取的心磁数据中选取某一 时刻， 取出各通道在该时刻的心磁数据； 计算各通道在胸腔表面对应的位置坐标； 令各通道在该时刻的心磁数据与计算得到的位置坐标对应， 并将对应后的各坐标按位 置排列以形成二维矩阵； 对所述二维矩阵的心磁数据进行二维插值处理； 在颜色表中查到与各插值后的心磁数据值相对应的颜色， 并在心磁图的对应位置显示 出来以绘制出彩。

4、色的心磁图。 2.根据权利要求1所述的方法， 其特征在于， 所述截取各通道一个周期的心磁数据， 包 括： 计算所有通道心磁数据的平均周期其中， Rti为第i个R 峰所在的时刻， N为心磁数据中R波的个数， K为通道个数； 以各通道心磁数据中的R峰所在时刻为参考点， 分别向前取分数倍个平均周期的数据， 向后取1-分数倍个平均周期的数据， 从而组成一个周期的心磁数据。 3.根据权利要求1所述的方法， 其特征在于， 所述计算各通道在胸腔表面对应的位置坐 标， 包括： 令： 所有通道在一个XY平面上， 第一个通道的坐标为(0,0)， 各通道间的距离相等且为 S； 则： 从原点(0,0)开始， 第一列通。

5、道的坐标为(0,0+iS)， 第二列通道的坐标为(S,0+iS)， 第三列通道的坐标为(2S,0+iS)， ， 第n列通道的坐标为(n-1)S,0+iS)； 其中， i为沿Y轴移 动的步数。 4.根据权利要求3所述的方法， 其特征在于， 所述将对应后的各坐标按位置排列以形成 二维矩阵， 包括： 令对应完成后的心磁数据沿Y轴正向及X轴正向排列， 将其坐标分布表示为C(x1,y1)、 (x2,y2)、 (x3,y3).(xK,yK)； 在坐标分布C中搜索与坐标(Xi， Yi)相等的坐标所对应的通道， 其中， (Xi， Yi)为通道i在 XY平面上的排列坐标； 当找到相应的坐标时跳出搜索循环， 并将。

6、通道i的心磁数据MCGi放到 一维数组的第k个位置， 其中， k为搜索循环的执行次数； 重复执行， 直到将所有的心磁数据 放到一维数组对应的位置上； 从获取的一维数组中取出前n个数据放在第一行， 再取出第n+1到第2n个数据放在第二 行， 以此类推， 直到将一维数组按Y轴负向及X轴正向排列成二维矩阵。 5.根据权利要求1所述的方法， 其特征在于， 所述对二维矩阵的心磁数据进行二维插值 处理， 包括： 令： 在XY平面上一系列坐标点(xi， yj)(i0,1,2,.,n； j0,1,2,.,m)上的心磁数 据值为aijf(xi,yj)； 构造二元函数zz(x,y)， 使z(xi,yj)的值逼近a。

7、ij， 此时， z(xi,yj)即是所求的坐标点 (xi， yj)处的插值。 6.根据权利要求5所述的方法， 其特征在于， 构造的二元函数z为： 权 利 要 求 书 1/3 页 2 CN 107374610 A 2 其中，为X轴坐标区间x0x1x2.xn上的非零三次基样条函数，j(y)为Y 轴坐标区间y0y1y2.ym上的非零三次基样条函数， uij为待定系数， 通过对z(x, y)求导得到。 7.根据权利要求1所述的方法， 其特征在于， 所述绘制彩色的心磁图， 还包括： 找出图像 中的最大值和最小值， 分别在最大值处画+号， 最小值处画-， 并把最大值和最小值用带箭头 的直线连接起来， 其中。

8、， 箭头指向最小值。 8.一种心磁图的生成系统， 其特征在于， 包括： 数据获取模块， 用于从心磁图仪的各通道中分别截取一个周期的心磁数据， 从截取的 心磁数据中选取某一时刻， 取出各通道在该时刻的心磁数据； 坐标计算模块， 用于计算各通道在胸腔表面对应的位置坐标； 对应处理模块， 用于令各通道在该时刻的心磁数据与计算得到的位置坐标对应， 并将 对应后的各坐标按位置排列以形成二维矩阵； 插值处理模块， 用于对所述二维矩阵的心磁数据进行二维插值处理； 图像生成模块， 用于在颜色表中查到与各插值后的心磁数据值相对应的颜色， 并在心 磁图的对应位置显示出来以绘制出彩色的心磁图。 9.根据权利要求8所。

9、述的系统， 其特征在于， 所述截取各通道一个周期的心磁数据， 包 括： 计算所有通道心磁数据的平均周期其中， Rti为第i个R 峰所在的时刻， N为心磁数据中R波的个数， K为通道个数； 以各通道心磁数据中的R峰所在时刻为参考点， 分别向前取分数倍个平均周期的数据， 向后取1-分数倍个平均周期的数据， 从而组成一个周期的心磁数据。 10.根据权利要求8所述的系统， 其特征在于， 所述计算各通道在胸腔表面对应的位置 坐标， 包括： 令： 所有通道在一个XY平面上， 第一个通道的坐标为(0,0)， 各通道间的距离相等且为 S； 则： 从原点(0,0)开始， 第一列通道的坐标为(0,0+iS)， 第。

10、二列通道的坐标为(S,0+iS)， 第三列通道的坐标为(2S,0+iS)， ， 第n列通道的坐标为(n-1)S,0+iS)； 其中， i为沿Y轴移 动的步数。 11.根据权利要求10所述的系统， 其特征在于， 所述将对应后的各坐标按位置排列以形 成二维矩阵， 包括： 令对应完成后的心磁数据沿Y轴正向及X轴正向排列， 将其坐标分布表示为C(x1,y1)、 (x2,y2)、 (x3,y3).(xK,yK)； 在坐标分布C中搜索与坐标(Xi， Yi)相等的坐标所对应的通道， 其中， (Xi， Yi)为通道i在 XY平面上的排列坐标； 当找到相应的坐标时跳出搜索循环， 并将通道i的心磁数据MCGi放到。

11、 一维数组的第k个位置， 其中， k为搜索循环的执行次数； 重复执行， 直到将所有的心磁数据 放到一维数组对应的位置上； 权 利 要 求 书 2/3 页 3 CN 107374610 A 3 从获取的一维数组中取出前n个数据放在第一行， 再取出第n+1到第2n个数据放在第二 行， 以此类推， 直到将一维数组按Y轴负向及X轴正向排列成二维矩阵。 12.根据权利要求8所述的系统， 其特征在于， 所述对二维矩阵的心磁数据进行二维插 值处理， 包括： 令： 在XY平面上一系列坐标点(xi， yj)(i0,1,2,.,n； j0,1,2,.,m)上的心磁数 据值为aijf(xi,yj)； 构造二元函数z。

12、z(x,y)， 使z(xi,yj)的值逼近aij， 此时， z(xi,yj)即是所求的坐标点 (xi， yj)处的插值。 13.根据权利要求12所述的系统， 其特征在于， 构造的二元函数z(x,y)为： 其中，为X轴坐标区间x0x1x2.xn上的非零三次基样条函数，j(y)为Y 轴坐标区间y0y1y2.ym上的非零三次基样条函数， uij为待定系数， 通过对z(x, y)求导得到。 14.根据权利要求8所述的系统， 其特征在于， 所述图像生成模块还用于： 找出图像中的 最大值和最小值， 分别在最大值处画+号， 最小值处画-， 并把最大值和最小值用带箭头的直 线连接起来， 其中， 箭头指向最小值。

13、。 权 利 要 求 书 3/3 页 4 CN 107374610 A 4 心磁图生成方法及生成系统 技术领域 0001 本发明属于生物医学信号分析领域， 特别涉及一种应用于心磁图仪中的心磁图生 成方法及系统。 背景技术 0002 心肌细胞周期性的极化与去极化过程形成了心电图(Electrocardiograph， ECG)， 但由于心电信号经人体组织传出后信号损失较大， 而磁信号的传导不受肌肉等组织的影 响， 因此心磁图(magnetocardiogram， MCG)成为最近研究的一个热点。 0003 心磁图是利用多通道心磁图仪测得人体心脏上方不同区域的磁场强度， 并把这些 测得的磁场值按位置。

14、排列以不同的颜色显示出来。 研究表明， 心磁图在冠心病、 心律失常等 心脏疾病的诊断、 风险分级等方面存在明显的优势。 Brisinda等人发现心磁图可用于鉴别 左束支传导阻滞患者(他们心磁图的R峰位置磁场演化方向与正常人相反)1。 Park等人的 研究表明， 对于预测急性胸痛和无ST段抬高的冠心病患者， MCG优于ECG、 超声心动图和肌钙 蛋白测试； 此外， MCG可以准确检测伴有急性胸痛的束支传导阻滞患者， 而ECG则无法诊 断2。 Fainziberg等人利用4通道心磁图仪研究了123个冠心病患者和124个正常人， 并对受 试者的心磁图进行分类(0正常4完全异常)， 发现若以分类值1.。

15、75为阈值， 则MCG分类的敏 感性和一致性分别达到了76和813。 0004 上述研究表明， 心磁图在心脏疾病的诊断、 预测、 风险评估等方面有潜在优势， 因 此， 如何通过心磁图算法产生高质量的心磁图成为临床研究的一个热点。 中国专利文献号 CN1552285A公开(公告)日2004.12.08， 公开了一种心磁波组图及其分析方法， 但该方法仅 是将不同位置同一时刻测得的心磁信号画在同一个坐标轴上， 无法反映心脏上方不同位置 的磁场分布情况。 中国专利文献号CN102682425B授权(公告)日2014.10.15， 公开了一种心 磁图系统和创建心磁图图像的方法， 该方法使用一种基于模型学。

16、习的算法， 即使用基于毕 奥-萨伐尔定律而随机产生的大量高分辨率MCG图像来构造模型， 通过为模型配备稀疏测 量， 创建高分辨MCG图像。 但该方法要通过创建模型， 并且为模型配备稀疏测量才能产生心 磁图图像， 算法过于复杂， 处理时间长， 无法满足临床上实时分析心磁图图谱的要求。 因此， 如何以有限点位置的磁场强度值产生高分辨率、 实时性好的心磁图是我们需要解决的一个 难题。 0005 参考文献： 0006 1D.Brisinda and R.Fenici,“Magnetocardiographic study of patients with right and left bundle b。

17、ranch blocks,“in International Congress Series, 2007,pp.451-454. 0007 2J .W .Park ,P .M .Hill ,N .Chung ,P .G .Hugenholtz ,and F .Jung ,“ Magnetocardiography predicts coronary artery disease in patients with acute chest pain,“Annals of noninvasive electrocardiology,vol.10,pp.312-323,2005. 0008 3L.Fa。

18、inzilberg,I.Chaikovsky,S.Auth-Eisernitz,B.Awolin,D.Ivaschenko, 说 明 书 1/6 页 5 CN 107374610 A 5 and B.Hailer,“Sensitivity and specificity of magnetocardiography,using computerized classification of current density vectors maps,in ischemic patients with normal ECG and echocardiogram,“in International C。

19、ongress Series,2007,pp.468-471. 发明内容 0009 鉴于以上所述现有技术的缺点， 本发明的目的在于提供一种应用于心磁图仪的心 磁图生成方法及系统， 用于解决现有技术中的上述问题， 从而为临床诊断心脏疾病提供参 考价值。 0010 为实现上述目的及其他相关目的， 本发明提供一种心磁图的生成方法， 包括： 从心 磁图仪的各通道中分别截取一个周期的心磁数据， 从截取的心磁数据中选取某一时刻， 取 出各通道在该时刻的心磁数据； 计算各通道在胸腔表面对应的位置坐标； 令各通道在该时 刻的心磁数据与计算得到的位置坐标对应， 并将对应后的各坐标按位置排列以形成二维矩 阵； 对。

20、所述二维矩阵的心磁数据进行二维插值处理； 在颜色表中查到与各插值后的心磁数 据值相对应的颜色， 并在心磁图的对应位置显示出来以绘制出彩色的心磁图。 0011 于本发明一实施例中， 所述截取各通道一个周期的心磁数据， 包括： 计算所有通道 心磁数据的平均周期其中， Rti为第i个R峰所在的时刻， N为 心磁数据中R波的个数， K为通道个数； 以各通道心磁数据中的R峰所在时刻为参考点， 分别 向前取分数倍个平均周期的数据， 向后取1-分数倍个平均周期的数据， 从而组成一个周期 的心磁数据。 0012 于本发明一实施例中， 所述计算各通道在胸腔表面对应的位置坐标， 包括： 令： 所 有通道在一个XY。

21、平面上， 第一个通道的坐标为(0,0)， 各通道间的距离相等且为S； 则： 从原 点(0,0)开始， 第一列通道的坐标为(0,0+iS)， 第二列通道的坐标为(S,0+iS)， 第三列通道 的坐标为(2S,0+iS)， ， 第n列通道的坐标为(n-1)S,0+iS)； 其中， i为沿Y轴移动的步数。 0013 于本发明一实施例中， 所述将对应后的各坐标按位置排列以形成二维矩阵， 包括： 令对应完成后的心磁数据沿Y轴负向及X轴正向排列， 将其坐标分布表示为C(x1,y1)、 (x2, y2)、 (x3,y3).(xK,yK)； 在坐标分布C中搜索与坐标(Xi， Yi)相等的坐标所对应的通道， 其。

22、 中， (Xi， Yi)为通道i在XY平面上的排列坐标； 当找到相应的坐标时跳出搜索循环， 并将通道i 的心磁数据MCGi放到一维数组的第k个位置， 其中， k为搜索循环的执行次数； 重复执行， 直 到将所有的心磁数据放到一维数组对应的位置上； 从获取的一维数组中取出前n个数据放 在第一行， 再取出第n+1到第2n个数据放在第二行， 以此类推， 直到将一维数组按Y轴负向及 X轴正向排列成二维矩阵。 0014 于本发明一实施例中， 所述对二维矩阵的心磁数据进行二维插值处理， 包括： 令： 在XY平面上一系列坐标点(xi， yj)(i0,1,2,.,n； j0,1,2,.,m)上的心磁数据值为 a。

23、ijf(xi,yj)； 构造二元函数zz(x,y)， 使z(xi,yj)的值逼近aij， 此时， z(xi,yj)即是所求 的坐标点(xi， yj)处的插值。 0015于本发明一实施例中， 构造的二元函数z为：其中， 说 明 书 2/6 页 6 CN 107374610 A 6 为X轴坐标区间x0x1x2.xn上的非零三次基样条函数， j(y)为Y轴坐标区 间y0y1y2.ym上的非零三次基样条函数， uij为待定系数， 通过对z(x,y)求导得 到。 0016 于本发明一实施例中， 所述绘制彩色的心磁图， 还包括： 找出图像中的最大值和最 小值， 分别在最大值处画+号， 最小值处画-， 并把。

24、最大值和最小值用带箭头的直线连接起 来， 其中， 箭头指向最小值。 0017 为实现上述目的及其他相关目的， 本发明提供一种心磁图的生成系统， 包括： 数据 获取模块， 用于从心磁图仪的各通道中分别截取一个周期的心磁数据， 从截取的心磁数据 中选取某一时刻， 取出各通道在该时刻的心磁数据； 坐标计算模块， 用于计算各通道在胸腔 表面对应的位置坐标； 对应处理模块， 用于令各通道在该时刻的心磁数据与计算得到的位 置坐标对应， 并将对应后的各坐标按位置排列以形成二维矩阵； 插值处理模块， 用于对所述 二维矩阵的心磁数据进行二维插值处理； 图像生成模块， 用于在颜色表中查到与各插值后 的心磁数据值相。

25、对应的颜色， 并在心磁图的对应位置显示出来以绘制出彩色的心磁图。 0018 于本发明一实施例中， 所述截取各通道一个周期的心磁数据， 包括： 计算所有通道 心磁数据的平均周期其中， Rti为第i个R峰所在的时刻， N为 心磁数据中R波的个数， K为通道个数； 以各通道心磁数据中的R峰所在时刻为参考点， 分别 向前取分数倍个平均周期的数据， 向后取1-分数倍个平均周期的数据， 从而组成一个周期 的心磁数据。 0019 于本发明一实施例中， 所述计算各通道在胸腔表面对应的位置坐标， 包括： 令： 所 有通道在一个XY平面上， 第一个通道的坐标为(0,0)， 各通道间的距离相等且为S； 则： 从原 。

26、点(0,0)开始， 第一列通道的坐标为(0,0+iS)， 第二列通道的坐标为(S,0+iS)， 第三列通道 的坐标为(2S,0+iS)， ， 第n列通道的坐标为(n-1)S,0+iS)； 其中， i为沿Y轴移动的步数。 0020 于本发明一实施例中， 所述将对应后的各坐标按位置排列以形成二维矩阵， 包括： 令对应完成后的心磁数据沿Y轴负向及X轴正向排列， 将其坐标分布表示为C(x1,y1)、 (x2, y2)、 (x3,y3).(xK,yK)； 在坐标分布C中搜索与坐标(Xi， Yi)相等的坐标所对应的通道， 其 中， (Xi， Yi)为通道i在XY平面上的排列坐标； 当找到相应的坐标时跳出搜。

27、索循环， 并将通道i 的心磁数据MCGi放到一维数组的第k个位置， 其中， k为搜索循环的执行次数； 重复执行， 直 到将所有的心磁数据放到一维数组对应的位置上； 从获取的一维数组中取出前n个数据放 在第一行， 再取出第n+1到第2n个数据放在第二行， 以此类推， 直到将一维数组按Y轴负向及 X轴正向排列成二维矩阵。 0021 于本发明一实施例中， 所述对二维矩阵的心磁数据进行二维插值处理， 包括： 令： 在XY平面上一系列坐标点(xi， yj)(i0,1,2,.,n； j0,1,2,.,m)上的心磁数据值为 aijf(xi,yj)； 构造二元函数zz(x,y)， 使z(xi,yj)的值逼近a。

28、ij， 此时， z(xi,yj)即是所求 的坐标点(xi， yj)处的插值。 0022于本发明一实施例中， 构造的二元函数z(x,y)为：其 中，为X轴坐标区间x0x1x2.xn上的非零三次基样条函数，j(y)为Y轴坐标 区间y0y1y2.ym上的非零三次基样条函数， uij为待定系数， 通过对z(x,y)求导 说 明 书 3/6 页 7 CN 107374610 A 7 得到。 0023 于本发明一实施例中， 所述图像生成模块还用于： 找出图像中的最大值和最小值， 分别在最大值处画+号， 最小值处画-， 并把最大值和最小值用带箭头的直线连接起来， 其 中， 箭头指向最小值。 0024 如上所。

29、述， 本发明的心磁图的生成方法及生成系统， 生成的心磁图图像清晰， 颜色 分布均匀， 在T波段可清晰的看到+-两个偶极子。 若从QRS波的S段到T段按设定步长画不同 时刻的心磁图， 则可以看到心脏活动电流的时空演化， 医生可根据心磁图的偶极子数目和 心脏活动电流的演化方向诊断该病人是否患心脏疾病。 附图说明 0025 图1显示为本发明一实施例中的心磁图的生成方法的流程图。 0026 图2显示为本发明一实施例中步长i5时K个通道在XY平面的坐标分布情况。 0027 图3显示为本发明一实施例中的心磁数据的插值算法的示意图。 0028 图4显示为用本发明的心磁图算法生成的正常人的心磁图。 0029 。

30、图5显示为用本发明的心磁图算法生成的心脏病患者的心磁图。 0030 图6显示为本发明一实施例中的心磁图的生成系统的模块图。 具体实施方式 0031 以下通过特定的具体实例说明本发明的实施方式， 本领域技术人员可由本说明书 所揭露的内容轻易地了解本发明的其他优点与功效。 本发明还可以通过另外不同的具体实 施方式加以实施或应用， 本说明书中的各项细节也可以基于不同观点与应用， 在没有背离 本发明的精神下进行各种修饰或改变。 需说明的是， 在不冲突的情况下， 以下实施例及实施 例中的特征可以相互组合。 0032 需要说明的是， 以下实施例中所提供的图示仅以示意方式说明本发明的基本构 想， 遂图式中仅。

31、显示与本发明中有关的组件而非按照实际实施时的组件数目、 形状及尺寸 绘制， 其实际实施时各组件的型态、 数量及比例可为一种随意的改变， 且其组件布局型态也 可能更为复杂。 0033 如图1所示， 本发明的心磁图生成方法， 具体包括5个步骤： 0034 步骤S1： 截取各个通道某个时刻的心磁数据。 具体实现方法为： 首先， 截取各个通 道一个周期的心磁数据MCG1MCGM(即从每个通道中取出的一个周期的心磁数据的个数为 M个， 对于36通道心磁采集系统， 通道总数为36， 每个通道一个周期的数据大概有800多个)， 并从截取的数据中选取某一时刻， 取出各个通道该时刻的数据； 然后， 计算各通道在。

32、胸部表 面对应的测量点的位置坐标。 0035 所述的截取各个通道一个周期的数据， 包括： 首先， 计算所有通道心磁数据的平均 周期然后， 以各通道心磁数据中的R峰所在时刻为参考点， 分别向前取2/5个平均周期的 数据， 向后取3/5个平均周期的数据， 组成一个周期。 0036所述的平均周期计算公式为，其中， Rti为第i个R峰所 在的时刻， N为心磁数据中R波的个数， K为总的通道个数。 说 明 书 4/6 页 8 CN 107374610 A 8 0037 所述的计算各通道在胸部上方对应的位置坐标： 假设所有通道(K个)在一个XY平 面上， 第一个通道的坐标为(0,0)， 各通道间的距离相等。

33、且为S， 则从原点开始沿X轴方向的 第一列通道的坐标为(0,0+iS)， i为沿Y轴方向向上移动的步数， i的取值为1K范围内的整 数； 第二列通道的坐标为(S,0+iS)， 以此类推， 可计算其它所有通道对应的XY平面坐标。 图2 显示了步长i5时， K个通道的坐标分布情况。 0038 步骤S2： 让各通道某个时刻的心磁数据与坐标对应， 即把各通道某个时刻的心磁 数据放到对应的位置坐标上， 形成二维矩阵。 具体实现方法为： 假设坐标对应完成后的数据 沿Y轴正向、 X轴正向排列， 其坐标分布为C(x1,y1)、 (x2,y2)、 (x3,y3).(xK,yK)； 假设通道i 在XY二维平面上的。

34、排列坐标为(Xi， Yi)。 首先， 在坐标分布C中寻找与坐标(Xi， Yi)相等的通 道， 当找到相应的坐标时， 跳出搜索循环， 假设此时循环次数为k， 则把通道i的心磁数据 MCGi放到一维数组的第k个位置； 然后， 按上述方法把所有的心磁数据放到一维数组对应的 位置上； 例如： 36通道心磁系统， 同一时刻每个通道就各取出此时刻的数据， 即每个通道取 出一个数据， 共取出36个数据， 在C上寻找与各个数据坐标(Xi,Yi)完全一样的坐标， 寻找到 即退出循环， 循环次数即是36个位置中对应的位置号1.36； 最后， 把上述一维数组按Y轴 负向、 X轴正向排列成二维矩阵， 即假设二维矩阵列。

35、数为n， 则首先取出一维数组的前n个数 据放在第一行， 然后取出一维数组的第n+1到第2n个数据放在第二行， 以此类推， 把所有数 据放在二维矩阵的相应位置。 0039 步骤S3： 对按位置排列好的二维心磁数据进行二维插值。 由于心磁图仪的信号通 道数有限， 例如： 36通道心磁图仪某一时刻仅有36个心磁数据， 用这36个心磁数据进行心磁 图成像， 图像分辨率太低， 因此， 需要对按位置排列好的二维心磁数据进行二维插值。 0040 请参阅图3， 本发明采用双三次样条插值算法对二维心磁数据进行插值， 具体实施 方法为： 假设在XY平面上一系列坐标点(xi， yj)(i0,1,2,.,n； j0,。

36、1,2,.,m； n和m的 取值为1K范围内的整数)上的心磁数据值为aijf(xi,yj)， 构造二元函数zz(x,y)， 使z (xi,yj)的值逼近aij， 即z(xi,yj)即是所求的坐标点(xi， yj)处的插值。 图3显示了插值前的采 集数据zi,j(xi， yj)与插值数据zim,n(xim,n， yim,n)的分布情况。 0041本发明中构造其中，为X轴坐标区间x0x1x2 .xn上的非零三次基样条函数，j(y)为Y轴坐标区间y0y1y2.ym上的非 零三次基样条函数， uij为待定系数， 可通过对z(x,y)求导计算得到。 0042 步骤S4： 对插值后的二维数据查询颜色表。 。

37、根据插值后的心磁数据值在颜色表中 查到相应的颜色在心磁图对应位置显示出来。 为提高心磁图的图像分辨率， 颜色表采用for 循环生成256色线性均匀颜色表， 每种颜色的差值可自行设定， 以达到较好的均匀性和连续 性。 0043 步骤S5： 画心磁图。 具体实现方法为： 首先， 设定图像的尺寸； 然后， 在图像的对应 位置上按心磁数据的值查找颜色表， 将相应的颜色值画在图像上； 最后， 找出图像中的最大 值和最小值， 分别在最大值处画+号， 最小值处画-， 并把最大值和最小值用带箭头的直线连 接起来， 箭头指向最小值。 0044 为验证心磁图生成算法的有效性， 我们利用9通道心磁图仪在徐汇区中心医。

38、院分 别采集了正常人和心脏病患者的心磁数据进行算法验证。 首先， 利用9通道心磁图仪在胸腔 说 明 书 5/6 页 9 CN 107374610 A 9 上方设定的四个固定点位置采集了四次心磁数据(每次采集9个通道的数据)， 共组成36通 道心磁数据； 然后， 按照本发明的心磁图生成算法， 对这36通道的心磁数据依次进行了截 取各个通道一个周期的数据， 并从截取的数据中选取某一时刻， 取出各个通道该时刻的数 据； 然后计算各通道在胸部上方对应的位置坐标。 让各通道某个时刻的数据与坐标对应。 对按位置排列好的二维心磁数据进行二维插值。 对插值后的二维数据查询颜色表。 画心磁图等5个步骤， 生成的。

39、正常人心磁图如附图4所示， 心脏病患者的心磁图如附图5所 示。 0045 如图6所示， 与上述方法实施例原理相似的是， 本发明还提供一种心磁图的生成系 统， 包括： 数据获取模块601、 坐标计算模块602、 对应处理模块603、 插值处理模块604、 图像 生成模块605。 由于前述方法实施例中的技术特征可以应用于本系统实施例， 因而不再重复 赘述。 0046 数据获取模块604从心磁图仪的各通道中分别截取一个周期的心磁数据， 从截取 的心磁数据中选取某一时刻， 取出各通道在该时刻的心磁数据。 坐标计算模块602计算各通 道在胸腔表面对应的位置坐标。 对应处理模块603令各通道在该时刻的心磁。

40、数据与计算得 到的位置坐标对应， 并将对应后的各坐标按位置排列以形成二维矩阵。 插值处理模块604对 所述二维矩阵的心磁数据进行二维插值处理。 图像生成模块605在颜色表中查到与各插值 后的心磁数据值相对应的颜色， 并在心磁图的对应位置显示出来以绘制出彩色的心磁图。 0047 综上所述， 本发明的心磁图生成方法及生成系统， 有效克服了现有技术中的种种 缺点而具高度产业利用价值。 0048 上述实施例仅例示性说明本发明的原理及其功效， 而非用于限制本发明。 任何熟 悉此技术的人士皆可在不违背本发明的精神及范畴下， 对上述实施例进行修饰或改变。 因 此， 举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完 成的一切等效修饰或改变， 仍应由本发明的权利要求所涵盖。 说 明 书 6/6 页 10 CN 107374610 A 10 图1 图2 说 明 书 附 图 1/3 页 11 CN 107374610 A 11 图3 图4 说 明 书 附 图 2/3 页 12 CN 107374610 A 12 图5 图6 说 明 书 附 图 3/3 页 13 CN 107374610 A 13 。