Forward and inverse problem for cardiac magnetic field and electric potential using two boundary element methods

This paper discusses the forward and inverse problem for cardiac magnetic fields and electric potentials.A torso-heart model established by boundary element method(BEM) is used for studying the distributions of cardiac magnetic fields and electric potentials.Because node-to-node and triangle-to-triangle BEM can lead to discrepant field distributions,their properties and influences are compared.Then based on constructed torso-heart model and supposed current source functional model-current dipole array,the magnetic and electric imaging by optimal constrained linear inverse method are applied at the same time.Through figure and reconstructing parameter comparison,though the magnetic current dipole array imaging possesses better reconstructing effect,however node-to-node BEM and triangleto-triangle BEM make little difference to magnetic and electric imaging.     

Cardiac magnetic source imaging based on current multipole model

It is widely accepted that the heart current source can be reduced into a current multipole. By adopting three linear inverse methods, the cardiac magnetic imaging is achieved in this article based on the current multipole model expanded to the first order terms. This magnetic imaging is realized in a reconstruction plane in the centre of human heart, where the current dipole array is employed to represent realistic cardiac current distribution. The current multipole as testing source generates magnetic fields in the measuring plane, serving as inputs of cardiac magnetic inverse problem. In the heart-torso model constructed by boundary element method, the current multipole magnetic field distribution is compared with that in the homogeneous infinite space, and also with the single current dipole magnetic field distribution. Then the minimum-norm least-squares (MNLS) method, the optimal weighted pseudoinverse method (OWPIM), and the optimal constrained linear inverse method (OCLIM) are selected as the algorithms for inverse computation based on current multipole model innovatively, and the imaging effects of these three inverse methods are compared. Besides, two reconstructing parameters, residual and mean residual, are also discussed, and their trends under MNLS, OWPIM and OCLIM each as a function of SNR are obtained and compared.     

基于心磁信号的心脏电流偶极子阵列成像及相关性质的研究

超导量子干涉器（SQUID）能探测到微弱的心脏磁场信号. 通过对所得的心磁信号进行分析,可为许多心脏疾病的诊断提供依据. 利用心磁信号,采用极小范数最小二乘法（MNLS）对心脏的电流偶极子阵列进行重建,从而实现了对心脏内部等效电流源的成像. 在使用MNLS进行电流偶极子阵列反演重建的过程中,反演所需的心磁信号,分别由单电流偶极子和电流多极子作为激发源模拟得到,以及由SQUID实际测量得到. 同时,对不同心磁信号反演得到的电流偶极子的分布规律进行了分析. 此外,还给模拟的人体外心磁信号施加了均匀噪声和随机噪声,研究不同信噪比的均匀噪声和随机噪声对电流偶极子阵列重建的影响. 关键词： 心磁信号 超导量子干涉器 电流偶极子阵列重建 极小范数最小二乘法     

基于心磁信号的心脏电流偶极子阵列成像及相关性质的研究

超导量子干涉器（SQUID）能探测到微弱的心脏磁场信号. 通过对所得的心磁信号进行分析,可为许多心脏疾病的诊断提供依据. 利用心磁信号,采用极小范数最小二乘法（MNLS）对心脏的电流偶极子阵列进行重建,从而实现了对心脏内部等效电流源的成像. 在使用MNLS进行电流偶极子阵列反演重建的过程中,反演所需的心磁信号,分别由单电流偶极子和电流多极子作为激发源模拟得到,以及由SQUID实际测量得到. 同时,对不同心磁信号反演得到的电流偶极子的分布规律进行了分析. 此外,还给模拟的人体外心磁信号施加了均匀噪声和随机噪声,研究不同信噪比的均匀噪声和随机噪声对电流偶极子阵列重建的影响.     

基于贪婪稀疏方法的心脏磁场源重构

心脏磁场源重构, 即通过人体胸腔表面的磁场阵列测量信号反演产生该磁场的电流源分布, 是一种无创地研究心脏电活动的方法. 本文用线性化方程描述测量磁场与心脏内部电流源的关系, 并通过一种贪婪优化的方法得到了源重构的稀疏解. 通过近似正交化过程和改变迭代算法中原子的选择方式, 降低了优化算法的复杂度, 可在保证源重构精度的情况下, 快速搜索源空间中强度比较大的位置. 文中通过一组正常人的心磁测量数据, 说明了源重构方法的有效性. 这组数据的研究结果表明, 强度大于65%的重构电流源的空间移动轨迹与心脏除极和复极电活动的传导过程基本符合. 其QRS 和ST-T 段的实测心磁图与重构电流源产生的磁场图的平均拟合优度分别为99.36%和99.78%. 关键词： 心磁图 逆问题 稀疏解 正交匹配追踪     

一种新的求解脑磁逆问题的搜索方法

给出一种新的求解真实头模型下脑磁逆问题的搜索方法.通过不同位置的源的相互关系，由上一个搜索源的计算结果通过简单计算，直接得到下一个搜索源的结果，避免了繁琐耗时的边界元积分方法，简化了求解过程，提高了求解速度. 关键词： 脑磁图 逆问题 搜索方法     

心脏磁场分布电流源重构及其精度分析

用测量到的心脏磁场信号重构其电流源是一种无创揭示心脏电活动的方法. 心脏电活动的时空信息可视化, 将有助于心脏功能的研究和心脏疾病诊断. 本文通过仿真实验研究了一种减时窗波束形成器对分布时变电流源的重构能力, 以及源重构结果与心室兴奋传播的关系. 采用元胞自动机模拟心室的兴奋传播, 产生分布随时间变化的电流源, 并用边界元法构建了一个心脏-躯干模型, 模拟体电导的作用. 仿真结果表明, 这种减时窗波束形成器能够重构分布时变电流源, 并达到较好的精度. 将该方法和无穷大均匀介质导联矩阵用于一例正常人的心脏磁场信号分析, 其结果可以反映心室兴奋传播的基本特征. 关键词： 心磁图 源重构 边界元 波束形成器     

用于心脏电活动成像的空间滤波器输出噪声抑制方法

利用人体表面测量的心脏磁场数据无创成像心脏电活动,需要解决的关键问题是提高其重建分布电流源偶极矩强度的分辨率.本文在最小方差波束成形(MVB)方法的基础上,提出了一种可抑制空间滤波器输出噪声功率增益(SONG)的波束成形方法,目的是通过构造一种新的滤波器权矩阵,约束空间滤波器的输出噪声功率增益,提高重建分布电流源偶极矩强度的分辨率,即分布电流源空间谱估计的源分辨能力,从而增强心脏电活动磁成像的分辨率.文中给出了电流源重建的理论分析和仿真结果;比较了该方法与MVB方法的差别;并给出了两个健康人36通道心脏磁场数据的电活动成像.结果表明,SONG方法分辨电流源的能力较强,能够观察到心脏电磁场信号R峰时刻健康人的心室电活动较强,心房电活动较弱等特征.     

多腔体心脏磁场模型的研究与应用

利用核磁共振图像(MRI)中提取的人体和心脏边界,根据边界元方法(BEM)建立了一个考虑左、右心房和心室的多腔体心脏磁场模型.分析了用该模型得到的36通道心脏磁场数据和特定时刻的磁场图.并在此基础上,研究了完全性右束支传导阻滞(CRBBB)和完全性左束支传导阻滞(CLBBB)病人ST-T段的心脏电活动.结果显示,用移动单电流偶极子模拟的单束支电兴奋传导所产生的磁场图与用超导量子干涉器(SQUID)测量的CRBBB/CLBBB病人数据绘制的心脏复极时的心磁图(MCG)十分相似.结果表明,该多腔体心脏BEM模型可用于CLBBB/CRBBB病人心脏磁场逆问题的研究.此外,文中给出了两个评价指标:测量平面上多腔体与单腔体的心脏磁场强度极大值之比,以及两种模型的36个测量点上磁场强度均方根之比.分析表明,多腔体心脏模型更贴近人体心脏的实际情况.该模型中心脏组织电导率参数的取值,以及等效电流偶极子的位置和个数决定了磁场的强度和分布.     

The deuteron and its magnetic moment the effects of non-local nuclear forces and of unitary transformations

The deuteron theory is renewed and completed for non-local Hamiltonians. A method to calculate deuteron wave functions is given. The formula for the current density and the magnetic dipole moment of a non-local potential is derived, and applied to the GCT, GT, and HJ-potential. The effects of a unitary transformation on the deuteron data — electric quadrupole moment, magnetic dipole moment, and radius — are shown. For the numerical calculations the Ristig-Kistler transformation and the HJ-potential were used. It is shown that an overall fit cannot be reached.     

Electromagnetic field perturbation by an internal void in a conducting cylinder excited by a wire loop

A thin prolate spheroidal void in an infinite conducting circular cylinder is used to model an internal flaw in a wire rope. The rope is excited by an electric ring current which is a model for a thin solenoid or multi-turn wire loop. The anomalous external fields are computed from the induced electric and magnetic dipole moments of the void. For this type of excitation, the induced axial magnetic dipole moment is the dominant contributor to the scattered field. The results have application to nondestructive testing of wire ropes.     

Scattering by a slender void in a homogeneous conducting wire rope

A thin prolate spheroidal void in an infinite conducting circular cylinder is used to model a broken strand in a wire rope. The rope is excited by an azimuthal magnetic line current which is a model for a thin toroidal coil. The anomalous external fields are computed from the induced electric and magnetic dipole moments of the void. The results have applications to nondestructive testing of wire ropes.     

Causal dynamic MRI reconstruction via nuclear norm minimization

This work addresses the problem of online reconstruction of dynamic magnetic resonance images (MRI). The proposed method reconstructs the difference between the images of previous and current time frames. This difference image is modeled as a rank deficient matrix and is solved from the partially sampled k-space data via nuclear norm minimization. Our proposed method has been compared against state-of-the-art offline and online reconstruction methods. Our method has similar reconstruction accuracy as the offline method and significantly higher accuracy compared to the online technique. It is about an order of magnitude faster than the online technique compared against. Our experimental data consisted of dynamic MRI data that were collected at 6 to 7 frames per second and having resolutions of 128×128 and 256×256 pixels per frame. Experimental evaluation indicates that our proposed method is capable of reconstructing 128×128 images at the rate of 4 frames per second and 256×256 images at the rate of 2 frames per second. The previous online method requires about 3.75 s for reconstructing each image. The improvement in reconstruction speed is clearly discernible. Moreover, our method has a reconstruction error that is about half that of the previous online method.     

托卡马克中等离子体的平衡重建

本文描述了用外部磁测量信号,和等离子体总电流,以及角向β_p的测量值,重建托卡马克中等高子体半衡的方法。对圆截面、非圆截面和偏滤器位形;对不同的电流分布分别作了重建计算。数值研究了测量误差、探针布置对重建结果的影响。     

Phase reconstruction via simulated tempering

A method for reconstructing phase images is described. Phase reconstruction is formulated as an optimization problem, finding the ground state of a locally constrained, finite size, spin-L Ising model under a non-uniform magnetic field. The optimization problem is solved by simulated tempering. Experiments on synthetic phase images show the effectiveness of the proposed algorithm.     

高能等离子体入射对磁场膨胀影响的模拟研究（已撤稿）

采用三维模型，使用混合网格质点法对等离子体入射偶极子磁场产生的磁场膨胀进行数值模拟.在模拟中考虑了高能等离子体注入两种不同类型磁场的情况：等离子体注入没有背景磁场的偶极子磁场和等离子体注入有背景磁场的偶极子磁场.研究表明背景磁场的存在不仅改变了粒子的分布，还改变了磁场膨胀的程度.还研究了注入的高能等离子体的速度对磁场膨胀的影响，结果表明入射的高能等离子体速度越大，磁场膨胀的程度就越大.对于低的入射速度，入射粒子在偶极子磁场中的回旋半径与偶极子磁场的特征长度相比较小，粒子被磁场束缚，对偶极子磁场的影响可以忽 关键词： 网格质点法 磁场膨胀 偶极子磁场     

Calculation of magnetic field for 4 T fast cycling superconducting dipole magnet

The problem of optimization of the two-dimensional magnetic field distribution of a dipole magnet with a field of 4 T and aperture diameter of 100–110 mm for a fast cycling synchrotron is considered. A singlelayer winding with a small number of turns is made from a hollow NbTi superconducting cable with an operating current up to 30 kA. The mathematical method providing optimization of the higher harmonics amplitude of the magnetic field by varying the position of current windings is described. The results of calculation of two-dimensional magnetic fields of the superconducting magnet are given.     

A multi-energy-level lattice Boltzmann model for Maxwell’s equations without sources

A lattice Boltzmann model for the Maxwell’s equations without sources is proposed by taking separate sets of distribution functions for the electric and magnetic fields and using the higher-order moment method. The higher-order moment method is based on the solution of a series of partial differential equations obtained by multi-energy-level techniques, multi-scale techniques, and Chapman–Enskog expansion, As numerical examples, some classical electromagnetic phenomena, such as the electric field and equipotential lines around an electrostatic dipole, the electric and magnetic fields around oscillating dipoles are given. These numerical results agree well with classical ones.     

真实头模型中的多电流偶极子脑磁源定位

给出了真实头模型下多电流偶极子定位的通用公式,并对多种源的情况进行了计算机模拟,逆问题的求解采用全局优化与局域优化相结合的优化方法,在较强的噪声背景下得到了较好的定位结果.其中,在偶极子源较少时,采用了球对称导体模型与真实头模型相结合的混合模型,进一步提高了求解速度 关键词： 脑磁技术 真实头模型 多偶极子脑磁源 优化方法     

Phase correction-based singularity function analysis for partial k-space reconstruction

Partial k-space acquisition is a conventional method in magnetic resonance imaging (MRI) for reducing imaging time while maintaining image quality. In this field, image reconstruction from partial k-space is a key issue. This paper proposes an approach fundamentally different from traditional techniques for reconstructing magnetic resonance (MR) images from partial k-space. It uses a so-called singularity function analysis (SFA) model based on phase correction. With such a reconstruction approach, some nonacquired negative spatial frequencies are first recovered by means of phase correction and Hermitian symmetry property, and then the other nonacquired negative and/or positive spatial frequencies are estimated using the mathematical SFA model. The method is particularly suitable for asymmetrical partial k-space acquisition owing to its ability of overcoming reconstruction limitations due to k-space truncations. The performance of this approach is evaluated using both simulated and real MR brain images, and compared with existing techniques. The results demonstrate that the proposed SFA based on phase correction achieves higher image quality than the initial SFA or the projection-onto-convex sets (POCS) method.