一种处理人体骨骼肌肉层各向异性导电性的数值方法

重点讨论了应用电磁场数值计算方法求解心电图正问题和逆问题中人体骨骼肌肉层的各向异性导电性的处理方法．文中应用有限元和边界元结合的方法,构造了一个包含不同的纤维方向组合的骨骼肌肉层的三维胸腔模型,并在此模型下引入了局部坐标内的各向异性导电率向整体坐标转换的方法．据此进行的模拟计算结果以图像的方式清晰地展示了各向异性导电性对体表电位图的影响．     

人体骨骼肌肉层对心电传导和体表电位标测的影响

应用电磁场数值计算方法求解三维心电图正问题,重点考察了人体骨骼肌肉层的各向异性导电性对心电传导的影响。推导了描述人类心电场的有限元模型和边界元模型,以及二者结合的形式,并在一个包含异性导电肌肉层的三维人体模下进行了模拟计算。     

人体骨骼肌肉层对心电传导和体表电位标测的影响

应用电磁场数值计算方法求解三维心电图正问题,重点考察了人体骨骼肌肉层的各向异性导电性对心电传导的影响。推导了描述人类心电场的有限元模型和边界元模型,以及二者结合的形式,并在一个包含各向异性导电性肌肉层的三维人体模型下进行了模拟计算。讨论了骨骼肌肉层在不同的肌肉纤维方向组合下的胸腔模型,以及纤维的各向异性导电性对体表电位图的影响     

基于神经网络的复介电常数重建算法

提出将一种神经网络的方法应用于微波断层成象,重建复的介电常数图象。微波断层成象发展了二十多年,仍然处于实验阶段,这是因为必须求解逆散射问题,它具有非适应性和非线性。为了改善其非适应性及非线性,引入适当的正则化过程并且应用神经网络的方法来求解这个问题。数值模拟实例证明了此方法的有效性。     

逆代谢工程——基于有用表型的定向遗传工程策略

代谢工程分为推理性代谢工程和逆代谢工程,逆代谢工程是避开对代谢网络充分认识的一种全新的遗传工程。本简要介绍了逆代谢工程的策略,并以实例详细分析了逆代谢工程的具体应用,最后,分析了逆代谢工程的问题和发展前景。     

体表电位标测系统及心电逆问题

一、前言 心电学的研究分正问题研究和逆问题研究,前者研究体表心电是如何产生和形成的;后者研究的是如何从体表电位推断心脏状态。广义地讲,目前临床医生根据无损地从体表测得的12导联心电图或心向量图来作诊断也是心电逆问题求解。只是这种对心脏状态的推断依据的是书本知识和医生的临床经验,而且一般只能作定性的推断。     

狗的正常心电图

本文概要叙述了正常狗的心电图研究的进展及其研究方法,阐明了正常狗心电图的特征,叙述了狗的心率、心电轴和心电图各波的波形、时间、电压等正常数值的范围,同时还指明了影响正常狗心电图的各种有关因素,为某些医学基础的实验研究提供方法和依据。     

园林植物三维数字模型的构建与应用探索

利用三维数字模型进行辅助设计、效果展示、施工指导和后期管护等已经成为风景园林行业发展趋势。植物是风景园林的关键组成要素,通过探讨园林植物三维数字模型构建与应用的必要性,以栾树（Koelreuteria paniculata）为例,利用三维植物建模软件PlantFactory对其进行三维数字模型构建探索,构建一个包含15种典型园林树种的小型三维数字模型库。最后,通过对园林植物三维数字模型在植物形态认知、植物景观设计表现和构建风景园林信息模型（LIM）3方面的应用,探索其适用性和有效性。     

逆重复伪随机序列信号相关分析法在探测蝇复眼电反应中的应用

伪白噪声测试信号的相关分析法,对于黑箱辨识问题来说,是较为合适的实验手段。隨着这种方法应用的扩展,特别是应用于生理系统,提出了如何进一步选择合适的输入信号的问题。逆重复伪随机二进制序列是较为合适的。本文介绍了逆重复伪随机序列信号相关分析法在探测光引起的蝇复眼反应中的应用,描述了信号的产生,测量并计算了模拟输出,它和实测结果甚为符合,说明这种方法对预测反应具有较高的精确性     

三维几何形态学概述及其在昆虫学中的应用

长期以来二维（two-dimensional, 2D）数据是几何形态学（geometric morphometrics）分析的最主要的数据类型,在推动几何形态学的发展过程中起到了奠基性的作用,并也解决了很多重大的科学问题,展示了几何形态学强大的科学计算能力与问题解决能力。但有些特殊的科学问题或者特殊的形态结构,无法通过二维数据完美解决,亟需大规模、大尺度三维（three-dimensional, 3D）数据的支持,这对几何形态学的三维化发展提出需求。更重要的是,随着三维数据获取成本的日渐降低,大量三维数据涌现出来。因此,三维几何形态学应运而生。本文对三维几何形态学的原理及其应用进行了概述,重点探讨了三维几何形态学与二维几何形态学的异同点,并对前者的两个发展阶段（少量样本的形态模拟与准定量比较及大量样本的定量比较）进行了概述,评价了四维数据和有限元等方法的应用,指出了该方法在昆虫学领域的发展潜力,最后对该方法在样本量增加、硬件提升、数据分辨率提高、新算法的开发、分析结果的呈现及3D打印等方面的发展趋势进行了展望。     

Computational issues of importance to the inverse recovery of epicardial potentials in a realistic heart-torso geometry

In vitro data from a realistic-geometry electrolytic tank were used to demonstrate the consequences of computational issues critical to the ill-posed inverse problem in electrocardiography. The boundary element method was used to discretize the relationship between the body surface potentials and epicardial cage potentials. Variants of Tikhonov regularization were used to stabilize the inversion of the body surface potentials in order to reconstruct the epicardial surface potentials. The computational issues investigated were (1) computation of the regularization parameter; (2) effects of inaccuracy in locating the position of the heart; and (3) incorporation of a priori information on the properties of epicardial potentials into the regularization methodology. Two methods were suggested by which a priori information could be incorporated into the regularization formulation: (1) use of an estimate of the epicardial potential distribution everywhere on the surface and (2) use of regional bounds on the excursion of the potential. Results indicate that the a posteriori technique called CRESO, developed by Colli Franzone and coworkers, most consistently derives the regularization parameter closest to the optimal parameter for this experimental situation. The sensitivity of the inverse computation in a realistic-geometry torso to inaccuracies in estimating heart position are consistent with results from the eccentric spheres model; errors of 1 cm are well tolerated, but errors of 2 cm or greater result in a loss of position and amplitude information. Finally, estimates and bounds based on accurate, known information successfully lower the relative error associated with the inverse and have the potential to significantly enhance the amplitude and feature position information obtainable from the inverse-reconstructed epicardial potential map.     

Inverse electrocardiographic transformations: dependence on the number of epicardial regions and body surface data points

The inverse problem of electrocardiography, the computation of epicardial potentials from body surface potentials, is influenced by the desired resolution on the epicardium, the number of recording points on the body surface, and the method of limiting the inversion process. To examine the role of these variables in the computation of the inverse transform, Tikhonov's zero-order regularization and singular value decomposition (SVD) have been used to invert the forward transfer matrix. The inverses have been compared in a data-independent manner using the resolution and the noise amplification as endpoints. Sets of 32, 50, 192, and 384 leads were chosen as sets of body surface data, and 26, 50, 74, and 98 regions were chosen to represent the epicardium.The resolution and noise were both improved by using a greater number of electrodes on the body surface. When 60% of the singular values are retained, the results show a trade-off between noise and resolution, with typical maximal epicardial noise levels of less than 0.5% of maximum epicardial potentials for 26 epicardial regions, 2.5% for 50 epicardial regions, 7.5% for 74 epicardial regions, and 50% for 98 epicardial regions. As the number of epicardial regions is increased, the regularization technique effectively fixes the noise amplification but markedly decreases the resolution, whereas SVD results in an increase in noise and a moderate decrease in resolution. Overall the regularization technique performs slightly better than SVD in the noise-resolution relationship.There is a region at the posterior of the heart that was poorly resolved regardless of the number of regions chosen. The variance of the resolution was such as to suggest the use of variable-size epicardial regions based on the resolution.     

人体躯干模型中肺的存在对体表电位分布的影响

在所建三维人体躯干模型的基础上,给出了如何应用边界元方法对非均匀人体场域进行求解。在设定心外膜电位分别呈现为单偶极子和双偶极子时,求出相应非均匀场域中的体表电位分布,并将它们与相同情况下均匀场域的体表电位分布进行分析比较。结果表明：躯干模型中肺的存在虽然对体表电位中极值的大小和位置没有太大影响,但却较大程序地改变了整个体表电位的分布状况,具体地说,就是由于肺的存在使得体表电位值较均匀时的相对误差高     

Comparison and testing of least-squares time domain inverse solutions in electrocardiography

The use of several mathematical methods for estimating epicardial ECG potentials from arrays of body surface potentials has been reported in the literature; most of these methods are based on least-squares reconstruction principles and operate in the time-space domain. In this paper we introduce a general Bayesian maximum a posteriori (MAP) framework for time domain inverse solutions in the presence of noise. The two most popular previously applied least-squares methods, constrained (regularized) least-squares and low-rank approximation through the singular value decomposition, are placed in this framework, each of them requiring the a priori knowledge of a ‘regularization parameter’, which defines the degree of smoothing to be applied to the inversion. Results of simulations using these two methods are presented; they compare the ability of each method to reconstruct epicardial potentials. We used the geometric configuration of the torso and internal organs of an individual subject as reconstructed from CT scans. The accuracy of each method at each epicardial location was tested as a function of measurement noise, the size and shape of the subarray of torso sensors, and the regularization parameter. We paid particular attention to an assessment of the potential of these methods for clinical use by testing the effect of using compact, small-size subarrays of torso potentials while maintaining a high degree of resolution on the epicardium.     

Electrotonic load triggers remodeling of repolarizing current Ito in ventricle

A change in activation sequence electrically remodels ventricular myocardium, causing persistent changes in repolarizing currents (T-wave memory). However, the underlying mechanism for triggering activation sequence-dependent remodeling is unknown. Optical action potentials were mapped with high resolution from the epicardial surface of the arterially perfused canine wedge preparation (n = 23) during 30 min of baseline endocardial stimulation, followed by 40 min of epicardial stimulation, and, finally, restoration of endocardial stimulation. Immediately after the change from endocardial to epicardial stimulation, phase 1 notch amplitude of epicardial cells was attenuated by 74 +/- 8% (P < 0.001) compared with baseline and continued to diminish during the period of epicardial pacing, suggesting progressive remodeling of the transient outward current (Ito). When endocardial pacing was restored, notch amplitude did not immediately recover but remained attenuated by 23 +/- 10% (P < 0.001), also consistent with a remodeling effect. Peak Ito current measured from isolated epicardial myocytes changed by 12 +/- 4% (P < 0.025), providing direct evidence for Ito remodeling occurring on a surprisingly short time scale. The mechanism for triggering remodeling of Ito was a significant reduction (by 14 +/- 4%, P < 0.001) of upstroke amplitude in epicardial cells during epicardial stimulation. Reduction in upstroke amplitude during epicardial pacing was explained by electrotonic load on epicardial cells by fully repolarized downstream endocardial cells. These data suggest a novel mechanism for triggering electrical remodeling in the ventricle. Electrotonic load imposed by a change in activation sequence reduces upstroke amplitude, which, in turn, attenuates Ito according to its known voltage-dependent properties, triggering downregulation of current.     

Cardiac electric field on the epicardial and body surfaces during the period of depolarization and repolarization of ventricular myocardium in pikes

Body surface and ventricular epicardial potential distributions during the electrocardiographic QRST interval were studied in pikes with the aid of potential mapping. The earliest epicardial activation was observed at the posterior base near the atrioventricular orifice. The areas of the earliest repolarization were found at the apex and the posterior base, whereas the area of the latest repolarization was detected at the anterior base. In the initial period of the QRS, the minimum was developed in the middle third of the right lateral body surface, and the maximum in the middle third of the ventral body surface. The body surface potential distribution during the ST-Twas characterized by the clear-cut negative potential zone in the cranial ventral area with the rest of the body surface having positive potentials, a pattern being largely unchanged throughout the period of the T-wave. The ventricular epicardial repolarization sequence differed from the activation sequence. The ventricular epicardial depolarization and repolarization sequences as well as epicardial potential distributions are expressed in the cardiac electric field on the body surface during the QRS and ST-T complexes.     

Use of endogenous NADH fluorescence for real-time in situ visualization of epicardial radiofrequency ablation lesions and gaps

Radiofrequency ablation (RFA) aims to produce lesions that interrupt reentrant circuits or block the spread of electrical activation from sites of abnormal activity. Today, there are limited means for real-time visualization of cardiac muscle tissue injury during RFA procedures. We hypothesized that the fluorescence of endogenous NADH could be used as a marker of cardiac muscle injury during epicardial RFA procedures. Studies were conducted in blood-free and blood-perfused hearts from healthy adult Sprague-Dawley rats and New Zealand rabbits. Radiofrequency was applied to the epicardial surface of the heart using a 4-mm standard blazer ablation catheter. A dual camera optical mapping system was used to monitor NADH fluorescence upon ultraviolet illumination of the epicardial surface and to record optical action potentials using the voltage-sensitive probe RH237. Epicardial lesions were seen as areas of low NADH fluorescence. The lesions appeared immediately after ablation and remained stable for several hours. Real-time monitoring of NADH fluorescence allowed visualization of viable tissue between the RFA lesions. Dual recordings of NADH and epicardial electrical activity linked the gaps between lesions to postablation reentries. We found that the fluorescence of endogenous NADH aids the visualization of injured epicardial tissue caused by RFA. This was true for both blood-free and blood-perfused preparations. Gaps between NADH-negative regions revealed unablated tissue, which may promote postablation reentry or provide pathways for the conduction of abnormal electrical activity.     

A cylindrical model for studying subendocardial ischaemia in the left ventricle

In this paper a mathematical model of a left ventricle with a cylindrical geometry is presented with the aim of gaining a better understanding of the relationship between subendocardial ischaemia and ST depression. The model is formulated as an infinite cylinder and takes into account the full bidomain nature of cardiac tissue, as well as fibre rotation. A detailed solution method (based on Fourier series, Fourier transforms and a one dimensional finite difference scheme) for the governing equations for electric potential in the tissue and the blood is also presented. The model presented is used to study the effect increasing subendocardial ischaemia has on the epicardial potential distribution as well as the effects of changing the bidomain conductivity values. The epicardial potential distributions obtained with this cylindrical geometry are compared with results obtained using a previously published slab model. Results of the simulations presented show that the morphologies of the epicardial potential distributions are similar between the two geometries, with the main difference being that the cylindrical model predicts slightly higher potentials.     

Catheter ablation of recurrent polymorphic tachycardia: Use of sodium channel blockade to organize the tachycardia: A case report

A 55 year old male presented with recurrent implantable cardioverter defibrillator (ICD) shocks due to polymorphic ventricular tachycardia (PMVT). He had undergone prior catheter ablation for VT three years ago. During the prior attempt he underwent voltage guided substrate ablation. With programmed ventricular extrastimulation (PVES), PMVT was repeatedly induced requiring DC shock. Intravenous procainamide was administered and PVES was repeated which induced sustained monomorphic ventricular tachycardia (MMVT). This VT had pseudo delta waves with maximum deflection index of 0.68, suggestive of epicardial origin. Activation mapping was performed epicardially. Presystolic potentials were recorded in mid anterolateral wall of left ventricular epicardial region. Radiofrequency (RF) ablation at this site terminated the VT. Post ablation there was no inducible tachycardia and patient is free of arrhythmias during 2 years of follow-up.     

The effect of boundary conditions on epicardial potential distributions

This study presents a comparison of semi-analytical and numerical solution techniques for solving the passive bidomain equation in simple tissue geometries containing a region of subendocardial ischaemia. When the semi-analytical solution is based on Fourier transforms, recovering the solution from the frequency domain via fast Fourier transforms imposes a periodic boundary condition on the solution of the partial differential equation. On the other hand, the numerical solution uses an insulation boundary condition. When these techniques are applied to calculate the epicardial surface potentials, both yield a three well potential distribution which is identical if fibre rotation within the tissue is ignored. However, when fibre rotation is included, the resulting three-well distribution rotates, but through different angles, depending on the solution method. A quantitative comparison between the semi-analytical and numerical solutiontechniques is presented in terms of the effect fibre rotation has on the rotation of the epicardial potential distribution. It turns out that the Fourier transform approach predicts a larger rotation of the epicardial potential distribution than the numerical solution. The conclusion from this study is that it is not always possible to use analytical or semi-analytical solutions to check the accuracy of numerical solution procedures. For the problem considered here, this checking is only possible when it is assumed that there is no fibre rotation through the tissue.