CN114795183A - 一种磁电结合定位跟踪的方法和装置 - Google Patents

一种磁电结合定位跟踪的方法和装置 Download PDF

Info

Publication number
CN114795183A
CN114795183A CN202210453845.0A CN202210453845A CN114795183A CN 114795183 A CN114795183 A CN 114795183A CN 202210453845 A CN202210453845 A CN 202210453845A CN 114795183 A CN114795183 A CN 114795183A
Authority
CN
China
Prior art keywords
electric field
positioning
magnetic field
matching
coef
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.)
Pending
Application number
CN202210453845.0A
Other languages
English (en)
Inventor
邓立
史天才
朱显钊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Jinjiang Electronic Science and Technology Co Ltd
Original Assignee
Sichuan Jinjiang Electronic Science and Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan Jinjiang Electronic Science and Technology Co Ltd filed Critical Sichuan Jinjiang Electronic Science and Technology Co Ltd
Priority to CN202210453845.0A priority Critical patent/CN114795183A/zh
Publication of CN114795183A publication Critical patent/CN114795183A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/06—Devices, other than using radiation, for detecting or locating foreign bodies ; determining position of probes within or on the body of the patient
    • A61B5/061—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body
    • A61B5/062—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body using magnetic field
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/06—Devices, other than using radiation, for detecting or locating foreign bodies ; determining position of probes within or on the body of the patient
    • A61B5/061—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body
    • A61B5/063—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body using impedance measurements

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Human Computer Interaction (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明涉及磁电结合医疗定位导航领域，具体涉及一种磁电结合定位跟踪的方法和装置，包括磁场定位模块、电场定位模块以及磁电匹配定位模块，磁电匹配定位模块被配置为针对定位空间，基于磁场定位模块得到磁传感器的磁场位置数据，与电场定位模块得到电场阻抗数据建立匹配关系，通过匹配关系对电场定位进行校准；具体匹配方法包括通过对定位空间划分单元格，记录空间点位置对应的磁场位置数据和电场阻抗数据，计算该单元格的匹配关系，最后利用该单元格的匹配关系，经运算获得校准后的电场定位结果；本发明解决了电场非线性定位的问题，提高电场定位的准确性。

Description

一种磁电结合定位跟踪的方法和装置
技术领域
本发明涉及一种磁电结合医疗定位导航领域，特别是一种磁电结合定位跟踪的方法和装置。
背景技术
目前医疗定位导航有基于电阻抗的定位系统包括在患者体外的一对或多对身体表面电极(例如，贴片)，附接到患者身体的参考传感器(例如，另一贴片)，以及附接到医疗设备的一个或多个传感器(例如，电极)。通过在电极对上施加电流，测量在设备电极处感应的相应电压(即相对于参考传感器)，处理测量的电压来确定定位。
而基于磁场的定位系统通常包括一个或多个磁场发生器，该磁场发生器放置在患者床或操作环境的其它组件附近，还包括与医疗设备耦合的一个或多个磁场检测线圈，磁场发生器与医疗设备耦合，并且检测线圈可以附接到操作环境的组件或放置在该操作环境的组件附近。发生器在解剖区域中提供磁场，检测线圈产生感应信号，系统处理这些信号，以产生与线圈相关联的一个或多个定位信息。与基于电阻抗的系统不同，在坐标系是相对于应用了身体表面电极的患者的情况下，基于磁场的系统具有与患者无关的坐标系。
基于电阻抗和基于磁场二者的定位系统提供优点，例如，基于电阻抗的系统提供了在多个医疗设备上同时定位相对大量传感器的能力。然而，由于基于电阻抗的系统利用人体中的电流流动，因此此类系统可能受到电干扰，会导致基于位置测量呈现的几何形状和表示可相对于受检区域的实际图像会发生失真。另一方面，基于磁场的坐标系不取决于患者解剖结构的特性，并且通常提供提高的准确度。然而，基于磁场的定位系统通常仅限于跟踪相对较少的传感器。
发明内容
本发明的目的在于：针对现有技术存在的基于电场定位系统受到电干扰，导致位置测量不准确的问题，提供一种磁电结合定位跟踪装置和方法。
为了实现上述目的，本发明采用的技术方案为：
一种磁电结合定位跟踪装置，包括磁场定位模块、电场定位模块以及磁电匹配定位模块；
磁电匹配定位模块，被配置为在建立磁电匹配关系阶段，对定位空间进行单元格划分，基于磁场定位模块得到磁场位置数据，与电场定位模块得到电场阻抗数据，建立各个单元格内的匹配关系；在校准电场阶段，通过电场定位模块确定测量对象所在的单元格位置，利用该单元格内的匹配关系校准其电场定位。
作为本发明的优选方案，一种磁电结合定位跟踪方法，包括：
磁场定位，获取磁场位置数据，包括磁传感器的空间位置和摆放角度；
电场定位，获取电场阻抗数据；
建立匹配关系，对定位空间划分单元格，记录经过各个单元格内空间点位置的磁场位置数据和电场阻抗数据，分别计算各个单元格的匹配关系；
校准电场定位，利用单元格的匹配关系，校准该单元格内采集的电场阻抗数据，获得校准后的电场定位结果。
作为本发明的优选方案，一种磁电结合定位跟踪方法，利用毕奥-萨伐尔定律，获取所述磁场位置数据。
作为本发明的优选方案，一种磁电结合定位跟踪方法，磁场位置数据计算公式为：
Figure BDA0003619899560000031
Figure BDA0003619899560000032
Figure BDA0003619899560000033
Figure BDA0003619899560000034
Voli＝γ*(B(x,i)*cos(α)*cos(β)+B(y,i)*cos(α)*sin(β)+B(z,i)*sin(α))
其中，磁场生成器方向向量
Figure BDA0003619899560000035
(xi,yi,zi)为三维空间位置，i≥6；(x,y,z)为磁传感器的三维空间位置，(α，β)为磁传感器的极角和方位角，γ为增益系数；Voli为第i个磁场生成器生成磁场作用于磁传感器上产生的信号量。
作为本发明的优选方案，一种磁电结合定位跟踪方法，采用LM算法，求解磁传感器的空间位置和摆放角度。
作为本发明的优选方案，一种磁电结合定位跟踪方法，电场定位还包括以下步骤：
步骤S1，通过电极片上的磁传感器获得电极片的位置信息，筛选最优电极片的作为共地电极片，并记录到系统中，其最优的筛选评判方法表达式如下：
Figure BDA0003619899560000036
其中，f(n)表示选择第n个电极片作为共地电极片，n＝1,2,3,4,5,6；α、β、γ、μ为权重系数；
Figure BDA0003619899560000037
为第k个电极片与共地电极片之间连线单位向量
Figure BDA0003619899560000038
与第t个电极片与共地电极片之间连线单位向量
Figure BDA0003619899560000039
之间点乘结果；Dis(n,k)为第k个电极片与共地电极片之间距离；VAng表示
Figure BDA0003619899560000041
与
Figure BDA0003619899560000042
之间夹角组成集合的方差，k∈1,2,3,4,5,6，t∈1,2,3,4,5,6，k≠n,t≠n,k≠t；VDis表示其余电极片与共地电极片之间距离集的方差；
步骤S2，筛选出相对共地电极片正交的三个电极片，并记录到系统中，筛选评判方法表达式如下：
Figure BDA0003619899560000043
其中，nbest为共地电极片的编号，
Figure BDA0003619899560000044
表示
Figure BDA0003619899560000045
与
Figure BDA0003619899560000046
点乘结果，a,b,c为正交的三个电极片，满足a∈1,2,3,4,5,6，b∈1,2,3,4,5,6，c∈1,2,3,4,5,6，a≠b≠c≠nbest；
步骤S3，其余电极贴片相对于共地电极片采用分频发放方法或者分时发放方法发放激励，采集介入导管的电极相对于电极片的电场阻抗数据。
作为本发明的优选方案，一种磁电结合定位跟踪方法，匹配关系包括：
利用正交电极片组(a,b,c)采集的电场阻抗数据建立线性匹配关系：
Figure BDA0003619899560000047
利用正交电极片组(a,b,c)采集的电场阻抗数据建立非线性匹配关系：
Figure BDA0003619899560000048
利用正交电极片组(a,b,c)采集的两个位置或者两个时刻(m，n)的电场阻抗数据建立线性匹配关系：
Figure BDA0003619899560000049
利用所有激励下采集的电场阻抗数据建立线性匹配关系：
Figure BDA00036198995600000410
利用所有激励下采集的电场阻抗数据建立非线性匹配关系：
Figure BDA0003619899560000051
利用所有激励下采集的两个位置或者两个时刻(m，n)的电场阻抗数据建立线性匹配关系：
Figure BDA0003619899560000052
作为本发明的优选方案，一种磁电结合定位跟踪方法，利用匹配关系校准后直接输出电场定位结果，或者经加权组合处理后输出电场定位结果，加权组合处理的表达式为：
Out＝Coef1*Out1+Coef2*Out2+Coef3*Out3+Coef4*Out4+Coef5*Out5+Coef6*Out6
其中，Out1、Out2、Out3、Out4、Out5、Out6分别为匹配关系(1)、匹配关系(2)、匹配关系(3)、匹配关系(3)、匹配关系(5)、匹配关系(6)校准后的电场定位结果；Coef1、Coef2、Coef3、Coef4、Coef5、Coef6分别为电场定位结果Out1、Out2、Out3、Out4、Out5、Out6的权重系数。
作为本发明的优选方案，一种磁电结合定位跟踪方法，单元格内空间点位置的间距分布不小于0.5mm。
综上所述，由于采用了上述技术方案，本发明的有益效果是：通过对定位空间进行单元格划分，采集电场阻抗数据和磁场位置数据，分别计算各个单元格内电场阻抗数据和磁场位置数据的匹配关系，再对相应单元格内的电场定位进行校正，有效地提高了电场定位的准确度，同时利用各个单元格内的匹配关系去校准对应单元格内的电场定位，也提高了电场定位校准的精确度。
附图说明
图1是本发明的系统装置示意图。
图2是磁场发生器示意图。
图3是磁传感器定位示意图。
图4是电场定位示意图。
图5是磁场与电场匹配关系示意图。
图标：101-患者；102-电极片；103-介入导管；104-操控手柄；105-磁场驱动装置；106-磁场发生器；107-电场驱动装置；108-集线器；109-磁场定位求解模块；110-定位盒；111-电场定位求解模块；112-磁电匹配定位模块；104-操控手柄；104-操控手柄；104-操控手柄；104-操控手柄；104-操控手柄；104-操控手柄。
具体实施方式
下面结合附图，对本发明作详细的说明。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
实施例1
如图1所示，一种磁电结合定位跟踪装置包括磁场定位模块、电场定位模块以及磁电匹配定位模块112，100为系统示意图，配置有介入导管103，介入导管103的远端安装有磁传感器和电传感器，104为介入导管103的操控手柄，101为患者；定位盒110可连接多个磁场信息采集传感器或电场信息采集传感器。
磁电匹配定位模块112，被配置为在建立磁电匹配关系阶段，对定位空间进行单元格划分，基于磁场定位模块得到磁场位置数据，与电场定位模块得到电场阻抗数据，建立各个单元格内的匹配关系；在校准电场阶段，通过电场定位模块确定测量对象所在的单元格位置，利用该单元格内的匹配关系校准其电场定位。
磁场定位模块包括磁场驱动装置105、磁场发生器106和磁场定位求解模块109，磁场驱动装置105与磁场发生器106连接，驱动磁场生成器106产生磁场，磁场定位求解模块109用于获取磁传感器的磁场位置数据。
具体的，磁场发生器106包括至少6个磁场生成器，磁场生成器可分组布置，也可分散布置。如图2所示的磁场生成器分组布置示意图，202A、202B、202C、202D为磁场生成器组，每组磁场生成器组包括三个磁场生成器，如202A磁场生成器组包括205、206、207三个磁场生成器，彼此位置近似正交。
电场定位模块包括电极片102、电场驱动装置107和电场定位求解模块111，如图4所示，电极片102包括401、402、403、404、405、406贴敷于介入对象101的体表，用于施加激励电信号，电极片102经过集线器108与电场驱动装置107连接，电场驱动装置107在电极片之间施加恒流或者恒压电信号激励，电场定位求解模块111用于获取电场阻抗数据。
一种磁电结合定位跟踪的方法包括：
磁场定位，磁场驱动装置105驱动所述磁场发生器106产生交变磁场，作用于磁传感器，经磁场定位求解模块109得到磁场位置数据，包括磁传感器的空间位置和摆放角度；
磁传感器为介入导管103远端安装的磁传感器，该磁传感器放入患者体内，医疗领域中常见安装有磁传感器的装置包括导管、导丝、导引器、探针等，应用领域包括心脏介入治疗导航、肺支气管定位导航、肾动脉消融导航等。
具体的，磁场驱动装置105驱动磁场生成器的方式分为两种：一种是分频驱动，即磁场驱动装置105为每个磁场生成器调制互不相同频率的信号，产生交变磁场，交变磁场作用于磁传感器产生感应电流，通过磁场信息采集传感器获得感应电压，感应电压经解调处理获得每一个磁场生成器作用于磁传感器上的电压；另一种是分时驱动，即磁场驱动装置105分时驱动每个磁场生成器，通过磁场信息采集传感器采集磁传感器上对应的感应电压。
进一步的，获得每个磁场生成器作用于磁传感器上产生的感应电压后，可根据磁偶极子模型，求解磁传感器的空间位置和摆放角度P(x,y,z,α,β)。
图3是磁场定位跟踪原理图。同理，302A、302B、302C、302D为磁场生成器组，每组包括3个磁场生成器，已知其中一个磁场生成器的位置和摆放角度P(xi,yi,zi,αi,βi)，304为介入导管103远端安装的磁传感器。
因磁场生成器与磁传感器之间的距离远大于磁场生成器自身尺寸，可将二者看做磁偶极子，根据毕奥-萨伐尔定律(Biot-Savart Law)，根据磁场生成器302位置和摆放角度P(xi,yi,zi,αi,βi)，得到归一化的磁场生成器方向向量
Figure BDA0003619899560000084
Dir(x,i)＝cos(αi)*cos(βi) (1)
Dir(y,i)＝cos(αi)*sin(βi) (2)
Dir(z,i)＝sin(αi) (3)
Figure BDA0003619899560000081
其中，(xi,yi,zi)为磁场生成器302的三维空间位置，(αi,βi)为磁场生成器302的摆放角度，即俯仰角度αi和旋转角度βi，i≥6。
磁传感器到磁场生成器302的距离为：
Figure BDA0003619899560000082
第i个磁场生成器生成磁场作用于磁传感器上产生的信号量Voli，对应电极片405中解调输出结果：
Figure BDA0003619899560000083
Figure BDA0003619899560000091
Figure BDA0003619899560000092
Figure BDA0003619899560000093
Voli＝γ*(B(x,i)*cos(α)*cos(β)+B(y,i)*cos(α)*sin(β)+B(z,i)*sin(α)) (10)
其中，(x,y,z)为磁传感器的三维空间位置，(α，β)为磁传感器的俯仰角度和旋转角度，γ为增益系数；由于，磁传感器的空间位置和摆放角度P(x,y,z,α,β,γ)包含6个待解未知量，设定有12个磁场生成器，获得包含6个未知量的12个方程，联立组成超定方程组：
Figure BDA0003619899560000094
求解超定方程组，可以按照一定筛选准则挑选出其中部分或全部方程联立求解，方程组个数大于等于6个，常用求解方法是LM(Levenberg-Marquardt)算法或其改进型，本专利采用其改进型，可在10次左右迭代获得收敛。根据非线性模型求解问题。
电场定位，电场驱动装置107对电极片102施加电信号激励，利用电场定位求解模块111获得介入导管103的电极相对于电极片102的电场阻抗数据；
具体的，如图4所示电场定位示意图，401、402、403、404、405、406为贴敷于介入对象101的体表的用于施加激励电信号的贴片，电场定位模块包括激励发放最优方式求解、激励发放和阻抗采集，其中，激励发放最优方式求解包括最优共地贴片求解与近似最正交贴片组合求解，方法步骤如下：
步骤S1，最优共地贴片求解，旨在筛选最优激励共地电极片。筛选原则是其余体表贴片与共地体表贴片距离尽量远，两个体表贴片与共地贴片连线夹角尽量大，距离间差异尽量小，夹角间差异尽量小；通过电极片上的磁传感器获得电极片的位置信息，筛选最优电极片的作为共地电极片，如图4中的共地电极片402，其最优的筛选评判方法表达式如下：
Figure BDA0003619899560000101
其中，n＝1,2,3,4,5,6，表示选择第n个体表电极片作为共地体表电极片。α、β、γ、μ为权重系数，分别描述夹角、距离、夹角方差VAng、距离方差VDis的重要程度，一般情况下，α＝β＝γ＝μ＝0.25；
Figure BDA0003619899560000102
为第k个电极片与共地电极片n之间连线单位向量
Figure BDA0003619899560000103
与第t个电极片与共地电极片之间连线单位向量
Figure BDA0003619899560000104
之间的点乘结果；Dis(n,k)为第k个电极片与共地电极片n之间距离；VAng表示
Figure BDA0003619899560000105
与
Figure BDA0003619899560000106
之间夹角组成集合的方差，k∈1,2,3,4,5,6，t∈1,2,3,4,5,6，k≠n,t≠n,k≠t；VDis表示其余电极片与共地电极片之间距离集的方差。
Figure BDA0003619899560000107
通过以上处理，得到激励发放最优共地电极片编号为nbest，记录到系统中。
步骤S2，近似最正交电极片组合求解，旨在在确定最优共地电极片后，筛选出相对共地电极片最正交的三个贴片，筛选评判方法可通过以下表达式描述：
Figure BDA0003619899560000108
其中，nbest为最优共地电极片的编号，
Figure BDA0003619899560000109
表示
Figure BDA00036198995600001010
与
Figure BDA00036198995600001011
的点乘结果，a,b,c为从除共地电极片外的其余5个贴片中选择的三个贴片，即满足a∈1,2,3,4,5,6，b∈1,2,3,4,5,6，c∈1,2,3,4,5,6，a≠b≠c≠nbest。
(a,b,c)＝argmin f(a,b,c,nbest) (15)
通过以上处理，得到相对共地贴片nbest最正交的电极片组合(a,b,c)，并记录到系统中。
步骤S3，激励发放和阻抗采集，其余电极片相对于共地电极片，采用分频发放方式或者分时发放方式发放激励，采集介入导管的电极相对于电极片的电场阻抗数据。
其中，其余电极片代表除共地电极片之外的电极片，分频发放方式为其余5个电极片以不同频率相对共地电极片发放激励，分时发放方式为依次轮替对共地电极片发放激励。
进一步的，阻抗采集包括采集介入导管103的电极在每对激励发放下相对电极片的电场阻抗数据；分频发放方式下，需对采集的信息进行解调，获得介入导管103的电极在每对激励下的电场阻抗数据Ele(v1,v2,v3,v4,v5)；分时发放方式可直接采得介入导管103电极在每对激励下的电场阻抗数据Ele(v1,v2,v3,v4,v5)。
建立匹配关系，对定位空间划分单元格，单元格内存储介入导管103经过的空间点位置，记录各个单元格内的空间点位置对应的磁场位置数据和电场阻抗数据，分别计算该单元格的匹配关系；
具体的，如图5所示将定位空间划分成单元格，501为心脏，502为心腔空间内单元格划分，单元格划分可采用多种尺寸划分方式，不同尺寸大小代表单元格的不同分辨率，如划分为2.5×2.5×2.5，5×5×5，10×10×10，单位为mm，503为其中一个单元格，单元格内存储介入导管103经过的空间点位置504，首先进行数据采集阶段，判断介入导管103所在的单元格，获取介入导管103在该单元格内经过的空间点位置504，记录每个空间点位置对应的磁场位置数据和电场阻抗数据，要求数据点之间满足空间间距分布，即空间点位置彼此之间间距不小于0.5mm。
其次，当单元格内数据点个数满足求解匹配关系最低数据个数要求时，求解更新该单元格的匹配关系信息，可以理解为，根据各个单元格内记录的数据点个数，选择能够求解的匹配方式，求解更新该单元格的匹配关系，本发明涉及的几种建立磁场与电场匹配关系的方式，描述如下：
a.利用最正交的电极片组采集的信息建立线性匹配关系1，描述如下：
Figure BDA0003619899560000121
b.利用最正交的电极片组采集的处信息建立非线性匹配关系2，描述如下：
Figure BDA0003619899560000122
c.利用最正交的电极片组采集的两个位置(或两个时刻)m和n的信息建立线性匹配关系3，描述如下：
Figure BDA0003619899560000123
d.利用所有激励下采集的阻抗信息建立线性匹配关系4，描述如下：
Figure BDA0003619899560000124
e.利用所有激励下采集的阻抗信息建立非线性匹配关系5，描述如下：
Figure BDA0003619899560000125
f.利用所有激励下采集的两个位置(或两个时刻)m和n的阻抗信息建立线性匹配关系6，描述如下：
Figure BDA0003619899560000126
校准电场定位，即应用各个单元格的匹配关系对该单元格电场定位进行校准，采集介入导管103相对各电极的电场信息，经处理后得到各电极在每组激励下的电场阻抗数据Ele(v1,v2,v3,v4,v5)和最正交激励组发放下的电场阻抗数据Ele(va,vb,vc)，找到当前电场阻抗数据对应的单元格，利用该单元内建立的匹配关系，经运算最终获得校准后的电场定位结果，包括如下：
a.应用匹配关系1，电场定位结果Out1(x,y,z)为：
Out1(x,y,z)＝Ele(va,vb,vc)*R3x3 (22)
b.应用匹配关系2，电场定位结果Out2(x,y,z)为：
Out2(x,y,z)＝Ele(va,vb,vc,va*va,vb*vb,vc*vc)*R6x3 (23)
c.应用匹配关系3，电场定位结果Out3(x,y,z)为：
Figure BDA0003619899560000131
其中，
Figure BDA0003619899560000132
或者
Figure BDA0003619899560000133
Figure BDA0003619899560000134
为与当前测量的电场阻抗数据Ele(va,vb,vc)最近的单元格内记录的历史数据。
d.应用匹配关系4，电场定位结果Out4(x,y,z)为：
Out4(x,y,z)＝Ele(v1,v2,v3,v4,v5)*R5x3 (25)
e.应用匹配关系5，电场定位结果Out5(x,y,z)为：
Out5(x,y,z)＝Ele(v1,v2,v3,v4,v5,v1*v1,v2*v2,v3*v3,v4*v4,v5*v5)*R10x3 (26)
f.应用匹配关系6，电场定位结果Out6(x,y,z,x0,y0,z0)为：
Figure BDA0003619899560000135
其中，(x0,y0,z0)表示多余输出，是上一个位置或上一个时刻对应的校准结果；
Figure BDA0003619899560000141
或者
Figure BDA0003619899560000142
Figure BDA0003619899560000143
为与当前测量的电场阻抗数据
Ele(v1,v2,v3,v4,v5)最近的单元格内记录的历史数据。
需要注意的是，以上电场校正结果可以直接作为最终校准结果，也和任意通过加权组合处理作为最终结果：
Out＝Coef1*Out1+Coef2*Out2+Coef3*Out3+Coef4*Out4+Coef5*Out5+Coef6*Out6(28)
其中，Out1、Out2、Out3、Out4、Out5、Out6分别为匹配关系(1)、匹配关系(2)、匹配关系(3)、匹配关系(3)、匹配关系(5)、匹配关系(6)校准后的电场定位结果；Coef1、Coef2、Coef3、Coef4、Coef5、Coef6分别为电场定位结果Out1、Out2、Out3、Out4、Out5、Out6的权重系数。
综上所述，本发明利用同一空间位置采集的磁场位置数据和电场阻抗数据，对空间位置进行单元格划分，分别建立各个单元格的磁场与电场位置的匹配关系，并利用匹配关系去校正对应的单元内采集的电场阻抗位置，实现提高电场定位校准的精确度，解决了电场非线性定位问题，进而提高电场定位的准确性。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种磁电结合定位跟踪装置，其特征在于，包括磁场定位模块、电场定位模块以及磁电匹配定位模块(112)；
所述磁电匹配定位模块(112)，被配置为在建立磁电匹配关系阶段，对定位空间进行单元格划分，基于磁场定位模块得到磁场位置数据，与电场定位模块得到电场阻抗数据，建立各个单元格内的匹配关系；在校准电场阶段，通过电场定位模块确定测量对象所在的单元格位置，利用该单元格内的匹配关系校准其电场定位。
2.一种磁电结合定位跟踪装置的方法，其特征在于，包括：
磁场定位，获取所述磁场位置数据，包括磁传感器的空间位置和摆放角度；
电场定位，获取所述电场阻抗数据；
建立匹配关系，对定位空间划分单元格，记录经过各个单元格内空间点位置的磁场位置数据和电场阻抗数据，分别计算各个单元格的匹配关系；
校准电场定位，利用单元格的匹配关系，校准该单元格内采集的电场阻抗数据，获得校准后的电场定位结果。
3.根据权利要求2所述的一种磁电结合定位跟踪方法，其特征在于，利用毕奥-萨伐尔定律，获取所述磁场位置数据。
4.根据权利要求2所述的一种磁电结合定位跟踪方法，其特征在于，所述磁场位置数据计算公式为：
Figure FDA0003619899550000011
Figure FDA0003619899550000012
Figure FDA0003619899550000013
Figure FDA0003619899550000021
Voli＝γ*(B(x,i)*cos(α)*cos(β)+B(y,i)*cos(α)*sin(β)+B(z,i)*sin(α))
其中，磁场生成器方向向量
Figure FDA0003619899550000022
(xi,yi,zi)为三维空间位置，i≥6；(x,y,z)为磁传感器的三维空间位置，(α，β)为磁传感器的极角和方位角，γ为增益系数；Voli为第i个磁场生成器生成磁场作用于磁传感器上产生的信号量。
5.根据权利要求4所述的一种磁电结合定位跟踪方法，其特征在于，采用LM算法，求解磁传感器的空间位置和摆放角度。
6.根据权利要求2所述的一种磁电结合定位跟踪方法，其特征在于，所述电场定位还包括以下步骤：
步骤S1，筛选最优的电极片作为共地电极片，并记录到系统中，其最优的筛选评判方法表达式如下：
Figure FDA0003619899550000023
其中，f(n)表示选择第n个电极片作为共地电极片，n＝1,2,3,4,5,6；α、β、γ、μ为权重系数；
Figure FDA0003619899550000024
为第k个电极片与共地电极片之间连线单位向量
Figure FDA0003619899550000025
与第t个电极片与共地电极片之间连线单位向量
Figure FDA0003619899550000026
之间点乘结果；Dis(n,k)为第k个电极片与共地电极片之间距离；VAng表示
Figure FDA0003619899550000027
与
Figure FDA0003619899550000028
之间夹角组成集合的方差，k∈1,2,3,4,5,6，t∈1,2,3,4,5,6，k≠n,t≠n,k≠t；VDis表示其余电极片与共地电极片之间距离集的方差；
步骤S2，筛选出相对所述共地电极片正交的三个电极片，并记录到系统中，筛选评判方法表达式如下：
Figure FDA0003619899550000031
其中，nbest为共地电极片的编号，
Figure FDA0003619899550000032
表示
Figure FDA0003619899550000033
与
Figure FDA0003619899550000034
点乘结果，a,b,c为正交的三个电极片，满足a∈1,2,3,4,5,6，b∈1,2,3,4,5,6，c∈1,2,3,4,5,6，a≠b≠c≠nbest；
步骤S3，对所述其余电极片相对于所述共地电极片发放激励，采集介入导管的电极相对于电极片的电场阻抗数据。
7.根据权利要求2所述的一种磁电结合定位跟踪方法，其特征在于，所述匹配关系包括：
利用正交电极片组(a,b,c)采集的电场阻抗数据建立线性匹配关系(1)：
Figure FDA0003619899550000035
利用正交电极片组(a,b,c)采集的电场阻抗数据建立非线性匹配关系(2)：
Figure FDA0003619899550000036
利用正交电极片组(a,b,c)采集的两个位置或者两个时刻(m，n)的电场阻抗数据建立线性匹配关系(3)：
Figure FDA0003619899550000037
利用所有激励下采集的电场阻抗数据建立线性匹配关系(4)：
Figure FDA0003619899550000038
利用所有激励下采集的电场阻抗数据建立非线性匹配关系(5)：
Figure FDA0003619899550000039
利用所有激励下采集的两个位置或者两个时刻(m，n)的电场阻抗数据建立线性匹配关系(6)：
Figure FDA0003619899550000041
8.根据权利要求7所述的一种磁电结合定位跟踪方法，其特征在于，利用所述匹配关系校准后直接输出电场定位结果，或者经加权组合处理后输出电场定位结果，加权组合处理的表达式为：
Out＝Coef1*Out1+Coef2*Out2+Coef3*Out3+Coef4*Out4+Coef5*Out5+Coef6*Out6
其中，Out1、Out2、Out3、Out4、Out5、Out6分别为匹配关系(1)、匹配关系(2)、匹配关系(3)、匹配关系(3)、匹配关系(5)、匹配关系(6)校准后的电场定位结果；Coef1、Coef2、Coef3、Coef4、Coef5、Coef6分别为电场定位结果Out1、Out2、Out3、Out4、Out5、Out6的权重系数。
9.根据权利要求2所述的一种磁电结合定位跟踪方法，其特征在于，单元格内空间点位置的间距分布不小于0.5mm。
CN202210453845.0A 2022-04-27 2022-04-27 一种磁电结合定位跟踪的方法和装置 Pending CN114795183A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210453845.0A CN114795183A (zh) 2022-04-27 2022-04-27 一种磁电结合定位跟踪的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210453845.0A CN114795183A (zh) 2022-04-27 2022-04-27 一种磁电结合定位跟踪的方法和装置

Publications (1)

Publication Number Publication Date
CN114795183A true CN114795183A (zh) 2022-07-29

Family

ID=82508785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210453845.0A Pending CN114795183A (zh) 2022-04-27 2022-04-27 一种磁电结合定位跟踪的方法和装置

Country Status (1)

Country Link
CN (1) CN114795183A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117357237A (zh) * 2023-10-25 2024-01-09 四川锦江生命科技有限公司 一种穿刺消融针

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117357237A (zh) * 2023-10-25 2024-01-09 四川锦江生命科技有限公司 一种穿刺消融针

Similar Documents

Publication Publication Date Title
US10524692B2 (en) Respiratory effect reduction in catheter position sensing
JP6301051B2 (ja) 複数の追跡方法を用いたプローブ追跡
CN101248996B (zh) 用于电磁跟踪方法和系统的线圈排列
JP5202880B2 (ja) インピーダンス校正用基準カテーテル
JP5642369B2 (ja) 電流局在化追跡機
AU2005234715B2 (en) Current-based position sensing
CN101657153B (zh) 用于非均匀场校正的系统和方法
AU2004203488B2 (en) Detection of metal disturbance in a magnetic tracking system
JP5052994B2 (ja) 周波数外挿を利用した歪みのない位置追跡方法およびシステム
US20140221803A1 (en) Current localization tracker
US7902817B2 (en) Electromagnetic tracking method and system
CN109091134B (zh) 组合躯干背心以标测心脏电生理学特性
JP2016105803A (ja) インピーダンスおよび磁界の測定を使用するカテーテルのナビゲーション
CN101120877A (zh) 利用冗余测量的抗失真位置跟踪
CN114795183A (zh) 一种磁电结合定位跟踪的方法和装置
Gürsoy et al. Imaging artifacts in magnetic induction tomography caused by the structural incorrectness of the sensor model
EP3791785B1 (en) Magnetic field distortion detection and correction in a magnetic localization system
CN110013311B (zh) 使用主成分分析改善基于阻抗的位置跟踪性能
WO2016200967A1 (en) Non-invasive system and method of spatial localization of specific electrocardiac elements
US20220211292A1 (en) Magnetic field distortion detection and correction in a magnetic localization system
EP4129168A1 (en) System and method for automatic localization of the spatial position of electrodes on a conductive body

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: No. 5, Wuke East 3rd Road, Wuhou District, Chengdu, Sichuan 610000

Applicant after: Sichuan Jinjiang Electronic Medical Device Technology Co.,Ltd.

Address before: 610045 Sichuan city of Chengdu province Wuhou District three Vuko East Road No. 5

Applicant before: SICHUAN JINJIANG ELECTRONIC SCIENCE AND TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information