CN112089491A - 电磁导航系统中磁场形变校准工具及方法 - Google Patents
电磁导航系统中磁场形变校准工具及方法 Download PDFInfo
- Publication number
- CN112089491A CN112089491A CN202010833891.4A CN202010833891A CN112089491A CN 112089491 A CN112089491 A CN 112089491A CN 202010833891 A CN202010833891 A CN 202010833891A CN 112089491 A CN112089491 A CN 112089491A
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- deformation
- coordinate
- navigation system
- electromagnetic navigation
- 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 28
- 238000013179 statistical model Methods 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 3
- 229920000122 acrylonitrile butadiene styrene Polymers 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 229910052755 nonmetal Inorganic materials 0.000 claims description 2
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 230000002452 interceptive effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2051—Electromagnetic tracking systems
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Surgery (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Robotics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明提供了一种电磁导航系统中磁场形变校准工具及方法,包括以下步骤:首先,在数据训练阶段将磁场发生器放置在磁场所使用的工作环境下,在坐标采集板上等距网格位置采集磁场坐标;得到三维等距网格线上所有点的磁场坐标值,根据网格线定义的物理坐标与所采集对应的磁场坐标求得磁场坐标的形变量;在不同的磁场工作环境下求得的形变量组合形成磁场形变统计模型;在使用运行阶段,采集一次坐标采集板上的局部磁场坐标获取局部形变量,通过得到稀疏表达得到当前的磁场坐标全局的形变量数值。本发明具有快速、准确和通用的优势,通过磁场形变数据的统计模型的稀疏求解,快速解决电磁导航易受到环境因素干扰的问题,使导航系统更加准确。
Description
技术领域
本发明涉及磁场形变校准技术领域,具体地,涉及一种电磁导航系统中磁场形变校准工具及方法。尤其地,涉及一种在医用手术导航系统下进行辅助准确定位的应用。
背景技术
随着科学技术的不断开发,微创精准介入定位手术越来越被人们重视,介入方法虽然好,但是定位难度十分大,因此,通过计算机辅助的方法能够帮助器械进行准确的空间定位。在实际中,辅助手术定位的方法主要有光学定位和电磁定位。光学定位的精度高,但是在实际使用的过程中容易受到环境中物体的遮挡影响,数据容易丢失,使得整个应用场景下的视野要求十分严格。而电磁定位通过电磁感应进行空间的方向和角度的定位,能够达到同级别的定位精度并且使用更加方便。
电磁定位由磁场发生器产生稳定的工作磁场,植入式的定位线圈能够在工作磁场范围能产生信号提供定位数据。然而在工作范围内的磁干扰物质能够产生感应涡流磁场,影响定位的精度。一般的电磁导航系统中往往会忽略掉磁场工作环境中的干扰物质对磁场稳定工作定位的影响,对于医疗手术导航系统这种需要高精度定位要求的使用场景来说,往往无法忽视这类干扰的影响。在手术室的环境中,磁场发生器支撑物,影像设备,比如CT床、C形臂等设备均可以引起磁场的干扰形变。
相关专利US973336中披露了一种对磁场形变补偿的方法,通过局部4个传感器进行小范围内修正,或通过光学定位的装置在磁场的范围内对线圈进行追踪,计算线圈在磁场范围之内的绝对误差量。此专利中的方法对局部或特定的磁场形变场景适用,不适用于广泛的磁场定位场景。实现一种在普遍电磁导航系统工作下的复杂电磁干扰手术环境的形变校准计算十分困难。
专利文献CN110544309A(申请号:201910688868.8)公开了一种基于大尺度网格模型的实时稀疏编辑方法,包括:构建由相同网格拓扑的模型组成的模型库;提取该模型的原始特征向量,将该原始特征向量转化为可以编码模型大尺度形变的第一特征向量,并形成该模型库的第一特征向量集合;对该第一特征向量集合进行全局稀疏优化,提取第二特征向量集合;确定目标模型上的控制顶点位置;以该目标模型的特征向量为该第二特征向量集合的一个线性组合向量,对该线性组合向量进行稀疏优化,得到编辑特征向量;根据该控制顶点位置和该编辑特征向量,得到该目标模型的形变模型。此专利披露了数字媒体技术领域和计算机图形学领域中用于网格模型库,通过局部控制点的形变生成符合模型库中变形特征的编辑后网格模型的方法。此技术面向的是形态网格点位置数据而非形变量数据,无法解决电磁导航系统中磁场形变量问题。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种电磁导航系统中磁场形变校准工具及方法。
根据本发明提供的电磁导航系统中磁场形变校准工具,包括:坐标采集板和磁场发生器,所述坐标采集板上有等距坐标采集孔和磁场发生装置卡槽,所述磁场发生器上带有螺纹松紧卡口。
优选的,通过调节螺纹松紧卡口,使磁场发生器在磁场发生装置卡槽上等距移动。
优选的,所述坐标采集板,使用的材料是非金属板,包括有机玻璃、PC板和ABS塑料。
根据本发明提供的电磁导航系统中磁场形变校准方法,包括:
步骤1:将磁场发生器放置在磁场中,在坐标采集板上等距网格位置采集磁场坐标;
步骤2:根据网格线定义的物理坐标与采集对应的磁场坐标求得磁场坐标的形变量;
步骤3:将不同磁场工作环境下求得的形变量组合形成磁场形变统计模型;
步骤4:在使用运行阶段,采集坐标采集板上的局部磁场坐标获取局部形变量,通过稀疏表示算法和磁场形变统计模型计算得到当前的磁场坐标全局的形变量数值,将磁场坐标数据与对应的形变量数值相加作为形变修正补偿,得到磁场坐标返回值进行磁场形变校准。
优选的,所述磁场形变统计模型中,磁场坐标值由插入锁定在坐标采集孔中依次等距移动采集,通过等距移动工装卡口得到等间隔三维物理坐标集合P,对应每个采集孔的磁场坐标集合M,磁场形变量为D=P-M。
优选的,在电磁导航系统的工作环境下,改变磁场发生器的位置和周围存在的磁场烦扰物质,重复使用坐标采集板获得多组磁场形变量形成磁场形变统计模型S={D1,D2,...,Dn},n表示采集的磁场形变量的组数。
优选的,所述磁场形变统计模型,采集到的磁场形变量是三维形变量,磁场形变统计模型中将所有采集到三维空间形变量拼接转为一维,多组空间形变量进行组合形成多维。
优选的,在使用运行阶段获取局部形变量的局部磁场为靠近电磁导航定位目标所处的感兴趣区域。
优选的,稀疏表示算法的公式为:
其中,d为运行使用阶段的局部形变量,s为形变统计模型中的对应形变量,w为稀疏表示系数,e为误差项,α和β分别为控制w和e的稀疏性正则化参数,使用同伦算法进行最优化求解得到稀疏表示系数w的最优值。
优选的,稀疏表示系数w为所有磁场形变统计模型S={D1,D2,...,Dn}中对应的线性组合系数,根据稀疏表示系数w与磁场形变统计模型相乘求得运行使用阶段的全局形变量T,计算公式为:T=S*w。
与现有技术相比,本发明具有如下的有益效果:本发明具有快速、准确和通用的优势,通过磁场形变数据的统计模型的稀疏求解,快速解决电磁导航易受到环境因素干扰的问题,使导航系统更加准确。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1是本发明的一个实施例的磁场形变校准计算方法流程图;
图2为磁场采集板结构示意图;
图3为磁场发生器结构示意图;
图4为磁场形变校准工具结构示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
实施例:
如图2-4,坐标采集板1具有等距坐标采集孔101,工装具有磁场发生装置卡槽102,通过调节磁场发生器的螺纹松紧卡口201,能使磁场发生器2在卡口102上等距移动。
本实施例使用NDIAurora磁场发生器,标准磁场有效范围为磁场发生器500*500*500mm的范围,电脑为Windows10,采用C++语言编程实现,坐标采集板使用有机玻璃材料按照符合NDIAurora磁场发生器规格的尺寸制作,坐标采集孔的间距为10mm,采集孔的数量为30*30。测试数据再不同的CT床边获得。
如图1,根据本发明提供的电磁导航系统中磁场形变校准方法,包括如下步骤:
测试阶段,将磁场发生器放置在磁场所使用的工作环境下,将坐标采集板放置在与磁场发生器卡槽150mm的位置,物理空间坐标原点为磁场发生器平面中心点为原点,采集板的水平方向为x轴,竖直方向为y轴,原理磁场发生器移动卡槽的方向为z轴。按照先横向后纵向的方式依次采集坐标,再将卡槽每次等距移动10mm采集一轮数据,得到总共为30*30*30维度的磁场坐标集合M和记录对应的物理空间坐标集合为P。将集合M和P中的(x,y,z)分量转为一维并依次拼接,形成81000*1维度的数据;
计算得到磁场形变量为D=P-M;
改变磁场发生器的位置与朝向,在不同的CT床周围存,重复使用坐标采集板获得多组磁场形变量形成磁场形变统计模型S={D1,D2,...,Dn};
使用运行阶段,距离在磁场发生器150mm至350mm的位置,200*200*200mm范围内的感兴趣区域随机抽样采集100个磁场坐标数据集合m和记录对应的物理坐标集合p,计算形变量数据d=p-m;
将对应物理坐标集合p的坐标值在磁场形变统计模型S中得到对应的局部磁场形变量集合s={d1,d2,...,dn};
使用稀疏表达方法,运用以下公式计算:
其中,d为运行使用阶段的局部形变量,s为形变统计模型中的对应形变量,w为稀疏表示系数,e为误差项,α和β分别为控制w和e的稀疏性正则化参数,优选地,使用α=5和β=0.2,使用同伦(Hotomopy)算法进行最优化求解得到稀疏组合系数w的最优值;
稀疏系数w为所有磁场形变统计模型S={D1,D2,...,Dn}中对应的线性组合系数,根据求得w可求得运行使用阶段的全局形变量T,计算公式为:
T=S*w。
在获取到全局的磁场形变量后,使用运行阶段通,在定位采集到磁场坐标数据后,将坐标数据与对应的形变类相加作为形变修正补偿,可以获取准确的磁场坐标返回值。
在本申请的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以,本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件,而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构;也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (10)
1.一种电磁导航系统中磁场形变校准工具,其特征在于,包括:坐标采集板和磁场发生器,所述坐标采集板上有等距坐标采集孔和磁场发生装置卡槽,所述磁场发生器上带有螺纹松紧卡口。
2.根据权利要求1所述的电磁导航系统中磁场形变校准工具,其特征在于,通过调节螺纹松紧卡口,使磁场发生器在磁场发生装置卡槽上等距移动。
3.根据权利要求1所述的电磁导航系统中磁场形变校准工具,其特征在于,所述坐标采集板,使用的材料是非金属板,包括有机玻璃、PC板和ABS塑料。
4.一种电磁导航系统中磁场形变校准方法,其特征在于,采用权利要求1-3中任一种或任多种所述的电磁导航系统中磁场形变校准工具,包括:
步骤1:将磁场发生器放置在磁场中,在坐标采集板上等距网格位置采集磁场坐标;
步骤2:根据网格线定义的物理坐标与采集对应的磁场坐标求得磁场坐标的形变量;
步骤3:将不同磁场工作环境下求得的形变量组合形成磁场形变统计模型;
步骤4:在使用运行阶段,采集坐标采集板上的局部磁场坐标获取局部形变量,通过稀疏表示算法和磁场形变统计模型计算得到当前的磁场坐标全局的形变量数值,将磁场坐标数据与对应的形变量数值相加作为形变修正补偿,得到磁场坐标返回值进行磁场形变校准。
5.根据权利要求4所述的电磁导航系统中磁场形变校准方法,其特征在于,所述磁场形变统计模型中,磁场坐标值由插入锁定在坐标采集孔中依次等距移动采集,通过等距移动工装卡口得到等间隔三维物理坐标集合P,对应每个采集孔的磁场坐标集合M,磁场形变量为D=P-M。
6.根据权利要求5所述的电磁导航系统中磁场形变校准方法,其特征在于,在电磁导航系统的工作环境下,改变磁场发生器的位置和周围存在的磁场烦扰物质,重复使用坐标采集板获得多组磁场形变量形成磁场形变统计模型S={D1,D2,...,Dn},n表示采集的磁场形变量的组数。
7.根据权利要求5所述的电磁导航系统中磁场形变校准方法,其特征在于,所述磁场形变统计模型,采集到的磁场形变量是三维形变量,磁场形变统计模型中将所有采集到三维空间形变量拼接转为一维,多组空间形变量进行组合形成多维。
8.根据权利要求4所述的电磁导航系统中磁场形变校准方法,其特征在于,在使用运行阶段获取局部形变量的局部磁场为靠近电磁导航定位目标所处的感兴趣区域。
10.根据权利要求9所述的电磁导航系统中磁场形变校准方法,其特征在于,稀疏表示系数w为所有磁场形变统计模型S={D1,D2,...,Dn}中对应的线性组合系数,根据稀疏表示系数w与磁场形变统计模型相乘求得运行使用阶段的全局形变量T,计算公式为:T=S*w。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010833891.4A CN112089491B (zh) | 2020-08-18 | 2020-08-18 | 电磁导航系统中磁场形变校准工具及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010833891.4A CN112089491B (zh) | 2020-08-18 | 2020-08-18 | 电磁导航系统中磁场形变校准工具及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112089491A true CN112089491A (zh) | 2020-12-18 |
CN112089491B CN112089491B (zh) | 2023-12-19 |
Family
ID=73754606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010833891.4A Active CN112089491B (zh) | 2020-08-18 | 2020-08-18 | 电磁导航系统中磁场形变校准工具及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112089491B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007113719A1 (en) * | 2006-03-31 | 2007-10-11 | Koninklijke Philips Electronics, N.V. | System for local error compensation in electromagnetic tracking systems |
US20120165656A1 (en) * | 2010-12-22 | 2012-06-28 | Avram Dan Montag | Compensation for magnetic disturbance due to fluoroscope |
CN105030331A (zh) * | 2015-04-24 | 2015-11-11 | 长春理工大学 | 位置传感器与三维腹腔镜摄像机标定装置及方法 |
CN108495585A (zh) * | 2016-01-26 | 2018-09-04 | 圣犹达医疗用品国际控股有限公司 | 磁定位系统中的磁场变形检测和校正 |
-
2020
- 2020-08-18 CN CN202010833891.4A patent/CN112089491B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007113719A1 (en) * | 2006-03-31 | 2007-10-11 | Koninklijke Philips Electronics, N.V. | System for local error compensation in electromagnetic tracking systems |
US20120165656A1 (en) * | 2010-12-22 | 2012-06-28 | Avram Dan Montag | Compensation for magnetic disturbance due to fluoroscope |
CN105030331A (zh) * | 2015-04-24 | 2015-11-11 | 长春理工大学 | 位置传感器与三维腹腔镜摄像机标定装置及方法 |
CN108495585A (zh) * | 2016-01-26 | 2018-09-04 | 圣犹达医疗用品国际控股有限公司 | 磁定位系统中的磁场变形检测和校正 |
Also Published As
Publication number | Publication date |
---|---|
CN112089491B (zh) | 2023-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11270465B2 (en) | Multiple camera calibration | |
CN105030331A (zh) | 位置传感器与三维腹腔镜摄像机标定装置及方法 | |
AU2012356327A1 (en) | Systems and methods for designing and generating devices using accuracy maps and stability analysis | |
CN1961335A (zh) | 用于对图像上的对象进行缩放的图像处理设备、成像系统、计算机程序、和方法 | |
Zhu et al. | Accuracy-and simplicity-oriented self-calibration approach for two-dimensional precision stages | |
Brecher et al. | Volumetric thermo-elastic machine tool behavior | |
Kwartowitz et al. | A novel technique for analysis of accuracy of magnetic tracking systems used in image guided surgery | |
Son et al. | Path planning of multi-patched freeform surfaces for laser scanning | |
CN112089491B (zh) | 电磁导航系统中磁场形变校准工具及方法 | |
Hu et al. | A holistic self-calibration algorithm for $ xy $ precision metrology systems | |
JP2007167374A (ja) | 医用画像処理装置 | |
Esmaeili et al. | An integrated geometric and hysteretic error model of a three axis machine tool and its identification with a 3D telescoping ball-bar | |
DE102021115214A1 (de) | Verzerrungskorrektur für elektromagnetische Felder mit Inside-Out-Tracking | |
Dorveaux et al. | On-the-field calibration of an array of sensors | |
Wang et al. | A fiducial-aided reconfigurable artefact for the estimation of volumetric errors of multi-axis ultra-precision machine tools | |
Sun et al. | Simulation of an optical-sensing technique for tracking surgical tools employed in computer-assisted interventions | |
Wang et al. | Space quantization between the object and image spaces of a microscopic stereovision system with a stereo light microscope | |
Hastedt | Monte-carlo-simulation in close-range photogrammetry | |
Pérez et al. | Lateral error compensation for stitching-free measurement with focus variation microscopy | |
Zhao | A high-precision 6D sensor based on optoelectronic measurement | |
JP2006139713A (ja) | 3次元物体位置検出装置、プログラム | |
Rudy et al. | The Off-Line Simulation on Measuring through Software PC-DMIS CAD++ V4. 3 | |
Bermejillo Seco | Automatic experimental characterization of a magnetic transport for low energy proton beams | |
Drégelyi-Kiss | Towards Traceable Dimensional Measurements by Micro Computed Tomography | |
Cao et al. | Dimensional modeling of single pass welding in wire and arc additive manufacturing |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |