CN101002674A - 结合lsqr和遗传算法的心外膜电位分布的检测方法 - Google Patents
结合lsqr和遗传算法的心外膜电位分布的检测方法 Download PDFInfo
- Publication number
- CN101002674A CN101002674A CNA2007100667388A CN200710066738A CN101002674A CN 101002674 A CN101002674 A CN 101002674A CN A2007100667388 A CNA2007100667388 A CN A2007100667388A CN 200710066738 A CN200710066738 A CN 200710066738A CN 101002674 A CN101002674 A CN 101002674A
- Authority
- CN
- China
- Prior art keywords
- lsqr
- potential
- genetic algorithm
- epicardial
- body surface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000002068 genetic effect Effects 0.000 claims abstract description 27
- 238000012546 transfer Methods 0.000 claims abstract description 6
- 238000001514 detection method Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 7
- 210000003516 pericardium Anatomy 0.000 claims description 5
- 239000004020 conductor Substances 0.000 claims description 3
- 238000002224 dissection Methods 0.000 claims description 2
- 238000011160 research Methods 0.000 description 8
- 230000000747 cardiac effect Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 3
- 210000000115 thoracic cavity Anatomy 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种结合LSQR和遗传算法的心外膜电位分布的检测方法。根据CT获取的解剖几何信息以及躯干容积导体的电磁特性，建立人体的心脏－躯干模型；采用边界元方法建立心外膜与体表电位之间的传递矩阵TBH，边界元法是有限元法在区域边界的一种典型应用；获取躯干的体表电位分布φB；采用LSQR迭代方法求解心外膜电位φH的解；将LSQR迭代方法得到的心外膜电位分布作为遗传算法的初始种群，采用遗传算法优化心外膜电位解的精度。本发明提出的基于LSQR和遗传算法相结合的心电逆问题求解方法，能够有效的解决心电逆问题中的病态特性，获取较精确的心外膜电位分布，正确的检测出心脏电活动信息，为心脏疾病的预防和治疗提供了临床信息。
Description
技术领域
本发明涉及用于诊断目的的探测方法，尤其是涉及一种结合LSQR和遗传算法的心外膜电位分布的检测方法。
背景技术
心电逆问题就是指从体表电位分布推断出其心脏内的电活动情况。临床上广为应用的心电图诊断其实就是一种心电逆问题的求解过程，只不过它是一种基于经验知识的定性诊断，而这里所指的心电逆问题，则是根据体表电位的分布，人体的几何形状以及躯干容积导体的电磁特性，通过数学物理方法求得心脏电活动的定量解。
由体表电位可无创伤地推出心脏状态，但是体表电位只是心脏电兴奋事件在人体表面的粗略投影，而且体表电位受到躯体的不均匀性、胸腔组织电导率不均匀性，胸腔模型的不规则性等各个方面的影响，因此很难由体表电位分布图推断心脏电兴奋的具体细节。相比较而言心外膜电位能够反映心电兴奋的细节，是心电源的精确反映，而且心外膜电位受躯体的不均匀性和胸腔各参数的影响较体表电位小地多，一些不能用体表电位解释的心脏电兴奋事件可用心外膜电位给出直接的解释。根据体表电位计算出心外膜电位，再根据心外膜电位推断心脏状态这种无创伤地研究方法对研究心脏正常和异常情况下的电活动具有重要意义，也是用来无创伤地得到心脏电活动信息的一个有效的临床工具，具有极大的临床和实验价值。基于心外膜电位的逆问题研究方法与等效心电源法相比有如下一些优点：(1)只要心外膜表面的几何形状是固定的，则用以心外膜电位表示的解就是唯一的；(2)心外膜电位直接与生理过程有关系；(3)不必对心电源增加许多限制；(4)在电位公式中隐含地考虑了胸腔内的血液不均匀性。因此基于心外膜电位的研究方法成为心电逆问题研究的主要手段。
用一个表示心外膜和躯干间区域的躯干模型，从体表电位计算出心外膜电位分布。如果，N点体表电位B和M点心外膜电位H，(N＞M)之间的理论关系可以通过线性矩阵TBH由方程给定：
B＝TBHH (1)
由B求出H就可得出逆问题的解，但(1)是病态的，即在体表微小的测量误差和几何噪声将对心外膜电位求解的干扰很大，因此对(1)式直接求逆是不合适的，因此需要采用正则化技术来求解心外膜电位的逆问题解。常见的正则化方法有：截断奇异值分解法(Truncated Singular-value Decomposition，TSVD)、Tikhonov正则化、共轭梯度最小二乘法(Conjugate Gradients Least Squares，CGLS)，LSQR迭代方法(Least Squares QR，LSQR)等。正则化技术的研究在心电逆问题的研究中有着非常重要的意义，国内外在这方面做了大量的研究。
发明内容
本发明的目的是提供一种结合LSQR和遗传算法的心外膜电位分布的检测方法，对由体表电位推断心外膜电位的分布的求解算法进行研究，为国内外的心电功能成像的研究提供一种稳定、精确的求解方法。
本发明解决其技术问题所采用的技术方案，该方法包括以下各步骤：
(1)根据CT获取的解剖几何信息以及躯干容积导体的电磁特性，建立人体的心脏-躯干模型；
(2)采用边界元方法建立心外膜与体表电位之间的传递矩阵TBH，边界元法是有限元法在区域边界的一种典型应用；
(3)获取躯干的体表电位分布B；
(4)采用LSQR迭代方法求解心外膜电位H的解；
(5)将LSQR迭代方法得到的心外膜电位分布作为遗传算法的初始种群，采用遗传算法进一步优化心外膜电位解的精度。
其中所述的将遗传算法与LSQR迭代方法相结合，以LSQR迭代求解得到的心电逆问题解作为遗传算法的初始种群，利用遗传算法进一步优化心电逆问题的解。
本发明具有的有益效果是：
本发明提出的基于LSQR和遗传算法相结合的心电逆问题求解方法，能够有效的解决心电逆问题中的病态特性，获取较精确的心外膜电位分布，正确的检测出心脏电活动信息，为心脏疾病的预防和治疗提供了很有意义的临床信息。
附图说明
图1是心脏-躯干模型。
图2是体表电位分布；(A)前胸电位分布，(B)背胸电位分布。
图3是采用LSQR获取的心外膜电位分布；(A)前面心外膜电位分布，(B)背面心外膜电位分布。
图4是经过遗传算法优化后的心外膜电位分布；(A)前面心外膜电位分布，(B)背面心外膜电位分布。
具体实施方式
本发明提出的一种结合LSQR和遗传算法的心外膜电位分布的检测方法，具体的实施步骤如下：
(1)由CT获取人体的心脏、肺、躯干的几何信息，进行三角形网格划分，建立心脏-躯干模型。体表划分成412个节点，820个网格；心脏表面由279个节点，531个的三角形网格组成；肺由297个节点，586个的三角形网格组成。如图1所示。
(2)采用边界元方法建立心外膜与体表电位之间的传递矩阵TBH。考虑到体表电位标测系统电极数量有限，将传递矩阵降维为：220×154。
(3)体表电位获取。在本发明中，设定心外电位分布H呈偶极子分布，通过传递矩阵B＝TBHH计算得到的体表电位，再叠加一定高斯随机噪声(信噪比(SNR)＝30dB)构成体表电位分布B。如图2所示。
(4)采用LSQR迭代方法(Least Squares QR，LSQR)求解心外膜电位H的解。LSQR方法是基于Lanczos过程的一种方法，Lanczos过程可将矩阵化为双对角阵，将这一双对角化思想应用于线性方程组和最小二乘问题，形成了LSQR(Least Squares QR)方法。LSQR方法的主要思想，即将一个复杂的最小二乘问题转化为一个简单的最小二乘问题，然后采用了标准的QR因子分解的方法求解最小二乘问题。采用L曲线方法获取LSQR方法所需要迭代的次数获取最优解。得到心电逆问题解的相对误差为(RE)：0.6691，相关系数为(CC)：0.8083。如图3所示。
(5)将LSQR迭代方法得到的心外膜电位分布作为遗传算法的初始种群，采用遗传算法进一步优化心外膜电位解的精度。由于遗传算法是一种随机最优化方法，通过繁殖、交叉能够得到一个更加精确的解。遗传算法的实现主要包括：参数的编码，初始种群的设定，适应度函数的设计，遗传操作，算法控制参数的设定和约束条件的处理等。根据遗传算法的特性将其利用到心电逆问题的求解中。采用遗传算法的参数设置如下：
经过遗传算法优化后得到心电逆问解的相对误差为(RE)：0.3601，相关系数为(CC)：0.9348。如图4所示。
Claims (2)
1、一种结合LSQR和遗传算法的心外膜电位分布的检测方法，其特征在于该方法包括以下各步骤：
(1)根据CT获取的解剖几何信息以及躯干容积导体的电磁特性，建立人体的心脏-躯干模型；
(2)采用边界元方法建立心外膜与体表电位之间的传递矩阵TBH，边界元法是有限元法在区域边界的一种典型应用；
(3)获取躯干的体表电位分布B；
(4)采用LSQR迭代方法求解心外膜电位H的解；
(5)将LSQR迭代方法得到的心外膜电位分布作为遗传算法的初始种群，采用遗传算法进一步优化心外膜电位解的精度。
2、根据权利要求1所述的一种结合LSQR和遗传算法的心外膜电位分布的检测方法，其特征在于：其中所述的将遗传算法与LSQR迭代方法相结合，以LSQR迭代求解得到的心电逆问题解作为遗传算法的初始种群，利用遗传算法进一步优化心电逆问题的解。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100667388A CN100464699C (zh) | 2007-01-16 | 2007-01-16 | 结合lsqr和遗传算法的心外膜电位分布的检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100667388A CN100464699C (zh) | 2007-01-16 | 2007-01-16 | 结合lsqr和遗传算法的心外膜电位分布的检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101002674A true CN101002674A (zh) | 2007-07-25 |
| CN100464699C CN100464699C (zh) | 2009-03-04 |
Family
ID=38702220
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2007100667388A Expired - Fee Related CN100464699C (zh) | 2007-01-16 | 2007-01-16 | 结合lsqr和遗传算法的心外膜电位分布的检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100464699C (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105025783A (zh) * | 2013-03-15 | 2015-11-04 | 圣犹达医疗用品电生理部门有限公司 | 用于采用医疗装置非接触标测的正则化方案 |
| CN105184089A (zh) * | 2015-09-22 | 2015-12-23 | 赵蒙海 | 异常病例在小样本中的多维量化检出方法 |
| CN105640543A (zh) * | 2016-03-14 | 2016-06-08 | 恩识医疗科技(上海)有限公司 | 一种心外膜心电成像方法及其成像系统 |
| CN105796094A (zh) * | 2016-05-13 | 2016-07-27 | 浙江大学 | 一种基于ecgi的室性早搏异常激动点定位方法 |
| CN106691438A (zh) * | 2016-12-07 | 2017-05-24 | 首都医科大学附属北京安贞医院 | 用于复杂心律失常的整体心脏三维标测系统 |
| CN106880353A (zh) * | 2015-12-16 | 2017-06-23 | 厦门大学 | 心电逆处理方法及装置 |
| CN108324263A (zh) * | 2018-01-11 | 2018-07-27 | 浙江大学 | 一种基于低秩稀疏约束的无创心脏电生理反演方法 |
-
2007
- 2007-01-16 CN CNB2007100667388A patent/CN100464699C/zh not_active Expired - Fee Related
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105025783A (zh) * | 2013-03-15 | 2015-11-04 | 圣犹达医疗用品电生理部门有限公司 | 用于采用医疗装置非接触标测的正则化方案 |
| CN105025783B (zh) * | 2013-03-15 | 2017-05-31 | 圣犹达医疗用品电生理部门有限公司 | 用于采用医疗装置非接触标测的正则化方案 |
| US10049771B2 (en) | 2013-03-15 | 2018-08-14 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Laplacian and Tikhonov regularization for voltage mapping with a medical device |
| CN105184089A (zh) * | 2015-09-22 | 2015-12-23 | 赵蒙海 | 异常病例在小样本中的多维量化检出方法 |
| CN106880353A (zh) * | 2015-12-16 | 2017-06-23 | 厦门大学 | 心电逆处理方法及装置 |
| CN105640543A (zh) * | 2016-03-14 | 2016-06-08 | 恩识医疗科技(上海)有限公司 | 一种心外膜心电成像方法及其成像系统 |
| CN105796094A (zh) * | 2016-05-13 | 2016-07-27 | 浙江大学 | 一种基于ecgi的室性早搏异常激动点定位方法 |
| CN106691438A (zh) * | 2016-12-07 | 2017-05-24 | 首都医科大学附属北京安贞医院 | 用于复杂心律失常的整体心脏三维标测系统 |
| CN106691438B (zh) * | 2016-12-07 | 2022-05-31 | 首都医科大学附属北京安贞医院 | 用于复杂心律失常的整体心脏三维标测系统 |
| CN108324263A (zh) * | 2018-01-11 | 2018-07-27 | 浙江大学 | 一种基于低秩稀疏约束的无创心脏电生理反演方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN100464699C (zh) | 2009-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100464699C (zh) | 结合lsqr和遗传算法的心外膜电位分布的检测方法 | |
| US10376171B2 (en) | Method and device for determining and presenting surface charge and dipole densities on cardiac walls | |
| Li et al. | Localization of the site of origin of cardiac activation by means of a heart-model-based electrocardiographic imaging approach | |
| TWI461180B (zh) | 增進電阻抗斷層攝影影像解析度的方法 | |
| US20070270703A1 (en) | Methods and apparatus of three dimensional cardiac electrophysiological imaging | |
| CN108324263B (zh) | 一种基于低秩稀疏约束的无创心脏电生理反演方法 | |
| JP7051707B2 (ja) | モバイル心臓モニタリング装置 | |
| RU2008146996A (ru) | Способ неинвазивного электрофизиологического исследования сердца | |
| JP2012179352A (ja) | 電流双極子を構築するシステムおよび電流双極子を構築する方法 | |
| Stenroos et al. | Boundary element computations in the forward and inverse problems of electrocardiography: comparison of collocation and Galerkin weightings | |
| CN110811596B (zh) | 基于低秩与稀疏约束和非局部全变分的无创心脏电位重建方法 | |
| KR101525485B1 (ko) | 3차원 심장 윤곽 재구성 방법 | |
| Beretta et al. | An asymptotic formula for boundary potential perturbations in a semilinear elliptic equation related to cardiac electrophysiology | |
| Semenov et al. | Microwave tomographic imaging of the heart in intact swine | |
| CN110393522B (zh) | 一种基于图总变分约束的无创心脏电生理反演方法 | |
| Horáček | The Forward and Inverse Problems of Electrocardiography | |
| CN101991412B (zh) | 一种心脏表面透壁电位分布的检测方法 | |
| Xie et al. | Non‐invasive reconstruction of dynamic myocardial transmembrane potential with graph‐based total variation constraints | |
| Sanjo et al. | Sensitivity of Electrocardiogram on Electrode-Pair Locations for Wearable Devices: Computational Analysis of Amplitude and Waveform Distortion | |
| Wong | Effects of unilateral and bilateral cochlear implantation on cortical activity measured by an EEG neuroimaging method in children | |
| Ferree et al. | Primer on electroencephalography for functional connectivity | |
| CN117523138A (zh) | 构造肿瘤拟合曲面的方法、电阻抗成像逆问题求解方法、检测系统 | |
| Ahmad | Inverse electrocardiography by simultaneous and iterative imposition of multiple constraints. | |
| Wang et al. | Noninvasive imaging of 3D cardiac electrophysiology | |
| Jiang et al. | Localization of the origin of ventricular premature beats by reconstruction of electrical sources using spatio-temporal map-based regularization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| ASS | Succession or assignment of patent right |
Owner name: ZHEJIANG HUAJI TIANFU INTELLIGENT TECHNOLOGY CO., Free format text: FORMER OWNER: ZHEJIANG UNIVERSITY Effective date: 20101221 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| COR | Change of bibliographic data |
Free format text: CORRECT: ADDRESS; FROM: 310027 NO. 38, ZHEDA ROAD, XIHU DISTRICT, HANGZHOU CITY, ZHEJIANG PROVINCE TO: 310013 8/F, TOWER B, LIANHE BUILDING, NO. 2, ZIJINGHUA ROAD, XIHU DISTRICT, HANGZHOU CITY, ZHEJIANG PROVINCE |
|
| TR01 | Transfer of patent right |
Effective date of registration: 20101221 Address after: 8, building 2, block B, The Union Buildings, No. 310013, Bauhinia Road, Hangzhou, Xihu District, Zhejiang Patentee after: Zhejiang Tianfu Tianfu Intelligent Technology Co., Ltd. Address before: 310027 Hangzhou, Zhejiang Province, Xihu District, Zhejiang Road, No. 38, No. Patentee before: Zhejiang University |
|
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090304 Termination date: 20140116 |