CN104665833B - 补偿由荧光检查仪产生的磁干扰 - Google Patents
补偿由荧光检查仪产生的磁干扰 Download PDFInfo
- Publication number
- CN104665833B CN104665833B CN201510029457.XA CN201510029457A CN104665833B CN 104665833 B CN104665833 B CN 104665833B CN 201510029457 A CN201510029457 A CN 201510029457A CN 104665833 B CN104665833 B CN 104665833B
- Authority
- CN
- China
- Prior art keywords
- camera
- locating pad
- fluoroscope
- orientation
- field
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 claims abstract description 57
- 238000003384 imaging method Methods 0.000 claims description 13
- 238000007689 inspection Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 abstract description 37
- 239000000523 sample Substances 0.000 abstract description 29
- 239000011159 matrix material Substances 0.000 description 22
- 238000013507 mapping Methods 0.000 description 18
- 238000010276 construction Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 11
- 238000005259 measurement Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 239000013598 vector Substances 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 230000005405 multipole Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 230000006978 adaptation Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 239000004020 conductor Substances 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002594 fluoroscopy Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- TVEXGJYMHHTVKP-UHFFFAOYSA-N 6-oxabicyclo[3.2.1]oct-3-en-7-one Chemical compound C1C2C(=O)OC1C=CC2 TVEXGJYMHHTVKP-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000036651 mood Effects 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012883 sequential measurement Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000004454 trace mineral analysis Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- 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
-
- 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7217—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/585—Calibration of detector units
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/587—Alignment of source unit to detector unit
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B2017/00681—Aspects not otherwise provided for
- A61B2017/00725—Calibration or performance testing
-
- 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
-
- 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/2055—Optical tracking systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/37—Surgical systems with images on a monitor during operation
- A61B2090/376—Surgical systems with images on a monitor during operation using X-rays, e.g. fluoroscopy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
- A61B2560/02—Operational features
- A61B2560/0223—Operational features of calibration, e.g. protocols for calibrating sensors
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/12—Arrangements for detecting or locating foreign bodies
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
- A61B6/487—Diagnostic techniques involving generating temporal series of image data involving fluoroscopy
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Signal Processing (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- High Energy & Nuclear Physics (AREA)
- Psychiatry (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Human Computer Interaction (AREA)
- Robotics (AREA)
- Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明涉及补偿由荧光检查仪产生的磁干扰。本发明公开了一种方法,所述方法包括使用多个磁发射器在某一区域中产生磁场,以及将场扰动元件引入所述区域中。所述方法包括表征每一个磁发射器在所述场扰动元件中的多幅图像,以及根据所述表征的图像计算所述区域中的反应磁场。所述方法还包括将探针定位在所述区域中,测量在所述探针处的被扰动的磁场,以及针对所测量的被扰动的磁场和所计算的反应磁场来确定所述探针的位置。
Description
技术领域
本发明总体涉及感测设置在活体内的物体的位置,具体地讲,涉及补偿对位置传感器造成影响的磁干扰。
背景技术
宽泛的医疗手术范围涉及在体内设置诸如传感器、管、导管、分配设备和植入物等物体。实时成像方法常常用来辅助医师在手术过程中观察物体及其周边环境。一些方法通过使用磁场的物体来跟踪物体。然而,磁场方面的干扰可能导致跟踪出错。
美国专利6,714,008描述了测定大物体磁场的梯度测量方法该专利被授予Holmes等人,其公开内容以引用方式并入本文。
美国专利申请2004/0034515描述了在电磁系统中估计位置和取向测量值的方法,该专利申请授予Bar Tal等人,其公开内容以引用方式并入本文。该方法包括使介于测量用模型与一个或多个测量值之间的差值最小化。
W.I.P.O(世界知识产权组织)专利公布WO/2004/006770描述了用于校准医学成像系统的方法,该专利公布授予Sati等人,其公开内容以引用方式并入本文。据说,针对C形臂成像装置的多个取向中的每一个取向,该方法都能够测定成像源的位置和(如果适用)磁场畸变。
美国专利申请2005/0107687描述了用于减少电磁跟踪器中磁场畸变的系统,该专利申请授予Anderson,其公开内容以引用方式并入本文。据说,该系统的某些实施例包括:跟踪分析单元,其用于分析仪器的跟踪行为;和跟踪修正单元,其用于补偿仪器的跟踪行为。
美国专利申请2007/0055125描述了电磁跟踪系统,该系统包括:场发生器和场传感器,其被布置用于分别产生和检测电磁场,该专利申请授予Anderson等人,其公开内容以引用方式并入本文。
美国专利申请2008/0183064描述了用于检测电磁(EM)场畸变的方法,该专利申请授予Chandonnet等人,其公开内容以引用方式并入本文。该方法包括:对设置在所关注体积内的传感器组件进行采样,以获得所关注体积内的电磁场的测量值;以及监测所述测量值,以检测所关注体积内的电磁场畸变。
美国专利6,147,480描述了用于跟踪物体的方法:由于响应能量场而引入的物品产生了干涉作用,在存在干涉作用的情况下,该方法在物体附近使用能量场,该专利被授予Osadchy等人,其公开内容以引用方式并入本文。
美国专利5,767,669描述了用于使用从固定位置产生的脉冲磁场来测定远程传感器的位置和取向的系统,该专利被授予Hansen等人,其公开内容以引用方式并入本文。由该系统单独感测和消减涡流畸变。
美国专利7,657,075描述了用于确定三维术前图像数据集的变换的方法,以获得三维图像数据与X射线成像系统的配准,该专利被授予Viswanathan,其公开内容以引用方式并入本文。
美国专利7,689,019描述了用于将物体的二维投影图像相对于同一物体的三维图像数据记录配准的方法和装置,该专利被授予Boese等人,其公开内容以引用方式并入本文。物体中包含的三维特征(其也可在三维图像中识别出)用符号重建。
美国专利申请2009/0010540描述了用于进行图像配准的方法,该专利申请授予Mullick等人,其公开内容以引用方式并入本文。该方法包括获得基准图像数据集和目标图像数据集,以及在基准图像数据集中限定用于所关注的区域的图像掩码。
美国专利申请2002/0172328描述了用于将仪器的空间坐标转换成其对应的X射线投影图像的方法,该专利申请授予Dekel,其公开内容以引用方式并入本文。据说,该方法的根据是:通过同时记录校准工具的空间坐标和X射线投影图像,将X射线束成像系统的坐标系与定位装置配准。
发明内容
本发明的实施例提供方法,该方法包括:
使用多个磁发射器在某一区域中产生磁场;
将场扰动元件引入该区域中;
表征每一个磁发射器在场扰动元件中的多幅图像;
根据表征的图像计算该区域中的反应磁场;
将探针定位在该区域中以及测量在该探针处的被扰动的磁场;以及
针对测量的被扰动的磁场和计算的反应磁场来测定探针的位置。
产生磁场的步骤通常包括在不存在场扰动元件的情况下测量磁场,以及在该区域中存在场扰动元件的情况下测量磁场。
在本发明所公开的实施例中,将场扰动元件引入该区域中的步骤包括测量场扰动元件相对于磁发射器限定的轴线的位置和取向。测量场扰动元件的位置和取向的步骤可包括针对反应磁场适应地计算位置和取向。
在本发明所公开的另一个实施例中,表征多幅图像的步骤包括假定多幅图像相对于彼此处于预定构型。预定构型可包括矩形,并且多幅图像可包括位于矩形的拐角和中心处的五幅图像。
在本发明所公开的又一个实施例中,计算反应磁场的步骤包括根据球谐函数展开来计算场。计算反应磁场的步骤可包括进行高达3阶的球谐函数展开。
作为另外一种选择或除此之外,根据球谐函数展开来计算场的步骤可包括将该展开配置成空间转移矩阵,该方法还可包括测定作为反应场矩阵的场扰动元件的性质,并且计算反应磁场的步骤可响应于空间转移矩阵与反应场矩阵的乘积。
通常,定位探针的步骤包括将探针定位在患者体内。
在一个实施例中,场扰动元件包括荧光检查仪的至少一部分。该方法可包括将荧光检查仪的荧光检查仪轴线与发射器所限定的轴线配准。配准荧光检查仪轴线的步骤可包括用固定地连接到发射器的摄像头来形成附接到荧光检查仪的基准标记的图像。引入场扰动元件的步骤可包括针对图像来测量场扰动元件的位置和取向。
根据本发明的实施例,还提供了设备,该设备包括:
多个磁发射器,其被配置用于在某一区域中产生磁场;
场扰动元件,其被引入所述区域中;和
处理器,其被配置用于:
表征每一个磁发射器在场扰动元件中的多幅图像,
根据表征的图像计算该区域中的反应磁场,
测量设置在该区域中的探针处的被扰动的磁场,以及
针对测量的被扰动的磁场和计算的反应磁场来测定探针的位置。
根据本发明的实施例,还提供了方法,该方法包括:
在定位垫上安装磁发射器,所述磁发射器被配置用于在患者体内产生磁场;
将定位垫摄像头以相应的固定取向附接到定位垫;
将可旋转摄像头连接到定位垫;
将基准标记附接到被配置用于使患者成像的荧光检查仪;
将荧光检查仪定位到不同位置中;以及
对每一个位置都执行如下操作:
将可旋转摄像头取向为已知取向,用可旋转摄像头和定位垫摄像头形成相应的基准标记的图像,以及
分析相应的图像,以将荧光检查仪的位置和取向与定位垫的轴线配准。
该方法可包括从定位垫移除可旋转摄像头,以及仅使用由定位垫摄像头形成的基准标记的图像来确定荧光检查仪相对于定位垫轴线的位置和取向。在一个实施例中,将定位垫摄像头附接到定位垫的步骤包括将可拆卸夹具附接到定位垫,以及通过用摄像头使夹具成像而将定位垫摄像头与相应的固定取向对齐。
根据本发明的实施例,还提供了设备,该设备包括:
定位垫;
磁发射器,其被配置用于在患者体内产生磁场,并且被安装在定位垫上;
定位垫摄像头,其以相应的固定取向被附接到定位垫;
可旋转摄像头,其被连接到定位垫;
荧光检查仪,其被配置用于使患者成像;
基准标记,其被附接到荧光检查仪;和
处理器,其被配置用于:
将荧光检查仪定位到不同位置中,并对每一个位置都执行如下操作:
将可旋转摄像头取向为已知取向,用可旋转摄像头和定位垫摄像头形成相应的基准标记的图像,以及
分析相应的图像,以将荧光检查仪的位置和取向与定位垫的轴线配准。
通过以下与附图结合在一起的本发明实施例的详细说明,将更全面地理解本发明。
附图说明
图1为根据本发明实施例的位置感测系统的示意性立体图;
图2为示出根据本发明实施例的探针的远端的示意性详细视图;
图3为根据本发明实施例的位置感测系统中使用的元件的示意图;
图4为根据本发明实施例测定探针远端的位置的工序的流程图;
图5为根据本发明实施例的可选位置感测系统的示意性立体图;
图6A和图6B为根据本发明实施例的可选位置感测系统的部分的示意图和系统校准中使用的元件的示意图;
图7为根据本发明实施例由控制单元在操作可选位置感测系统中执行的步骤的流程图;以及
图8A-8G示出了根据本发明实施例的矩阵元素。
具体实施方式
概述
本发明的实施例提供了用于对某区域中的磁场中产生的扰动进行补偿的方法。将扰动元件(通常为金属元件)引入到磁发射器所产生的场中会引起扰动。为了对扰动元件的存在进行补偿,本发明的实施例中使用的反应场模型假定每以个磁发射器都在元件中产生多幅发射器图像。该模型假定每一幅图像都产生相应的反应场,这些反应场总体起到扰动由发射器产生的场的作用。
每一幅图像通常都可表征为多极(即,双极、四极和/或更高阶的极)的组合。每一幅图像的特征还特别都取决于产生图像的发射器场。该模型通常通过假定场可由球谐函数展开表示,并根据图像的特征来计算来自多极图像中的每一幅的反应场。通常,通过在存在扰动元件以及在不存在扰动元件的情况下测量发射器发射的场来有效测定图像特征,这些特征可用空间转移矩阵和反应场矩阵的方式表示。
一旦图像被表征,设置在该区域中的探针就测量被扰动的磁场。然后,可针对测量的磁场和如上所述计算的反应场来计算探针在该区域中的位置。
系统说明
现在参见图1、图2和图3,图1为根据本发明实施例的位置感测系统20的示意性立体图,该系统被配置用于感测导管探针22的远端21的位置,图2为示出根据本发明实施例的远端21的示意性详细视图,图3为根据本发明实施例的系统20中使用的元件的示意图。通常由医务人员在系统20的操作阶段将探针22插入患者体腔或器官中,而在系统的校准阶段不提供探针22。为清楚起见,探针22在图1中用虚线表示。在医务人员(其可以操作位置感测系统20)施行手术过程中,患者通常平躺在手术台24上。为简单明了起见,患者和医务人员均未在图1中示出。
以举例的方式,在下文的描述中,假定探针22用于患者心室内进行的侵入式手术。将患者置于处于磁场中的手术台24上,该磁场是通过在患者下方设置定位垫26而产生的,该定位垫包含大致类似的交变磁场发射器线圈32A、32B、的构造28,这些线圈在本文也统称为磁场发射器线圈32。线圈32在区域30中产生其交变磁场,该区域在图1中以椭圆形示意性示出。图3示出了线圈32的示例性构造。
通常情况下以及如图3所示,构造28包括发射器线圈32,该发射器线圈成形为三组三轴线圈34、36和38,每一组三轴线圈都包括彼此正交的三个线圈。然而,除了构造28以外,磁场发射器线圈32的其他构造也是可行的,并且此类构造对于本领域普通技术人员而言是显而易见的。授予Govari的美国专利6,484,118描述了与构造28类似的线圈构造以及可用于系统20中的线圈32的其他构造,该专利被授予Govari,其公开内容以引用方式并入本文。
本文的描述假定远端21包括三个作为位于远端的电磁(EM)传感器46的部件的大致正交的线圈40、42和44。(远端通常包括其他元件,诸如图2所示的电极48。)线圈32产生的磁场在传感器46的线圈中根据线圈感测的场产生电信号。将来自传感器46线圈的电信号传送到控制单元50,该控制单元对这些信号进行分析,以便测定探针22的位置和取向的坐标。假定坐标参照一组正交的xyz轴线(如图3所示),这些轴线相对于定位垫26固定。
远端21中用于检测探针22位置和取向的线圈的其他构造在本领域是已知的。一种此类构造采用一个线圈,该线圈测量磁场的投影。本领域的普通技术人员在加以必要的变通后将能够对本描述进行修改,以阐释与线圈40、42和44例示的构造不同的线圈构造。
控制单元50包括处理器52,通常为具有适当信号处理电路的计算机。处理器使用存储器54,该存储器通常包括其内存储操作系统20的数据的易失性和非易失性数据存储器件。连接处理器以驱动控制台,该控制台可提供探针22的位置的视觉显示器56。
系统20包括荧光检查仪60,该荧光检查仪由荧光检查仪控制器61操作,并且能够产生手术台24上的患者的荧光透视图像。假定荧光检查仪控制器61为控制单元50的子单元。荧光检查仪60具有多个部分,包括准直X射线源62(在本文也称为准直器62)和检测器64。准直器和检测器通过另一部分即“C形臂”66连接在一起,“C形臂”66允许它们围绕两条轴线旋转,所述两条轴线即水平轴线68和垂直于纸面经过轴线68的轴线。C形臂也允许准直器和检测器在空间中平移,诸如沿平行于水平轴线的方向平移。在荧光检查仪围绕轴线68旋转的过程中,C形臂66使准直器和检测器保持彼此固定地对齐,并且保持彼此之间恒定的距离。通常,荧光检查仪60形成的图像可以用围绕轴线68旋转至任何取向的荧光检查仪形成,所述取向根据患者的需要和专业操作系统20的要求而选择。
通常,荧光检查仪60基本上在采用探针22的同一时间进行操作,以产生其成像。然而,荧光检查仪中靠近区域30的金属元件改变线圈32在该区域中产生的磁场。如果不补偿这些变化,则会造成测量的探针22的位置不准确。如本文所述,本发明的实施例在不论荧光检查仪相对于手术台24处于何种取向或位置的情况下,对荧光检查仪60引起的磁场变化进行了补偿,从而防止了测量的探针位置出现任何不准确。
通常,系统20包括其他元件,出于简洁起见,这些元件未在图中示出,并且下列描述会视需要提及。例如,系统20可以包括ECG监视器,该ECG监护仪被连接以接收来自一个或多个身体表面电极的信号,以便向控制单元50提供ECG同步信号。
图1的配置为示例性配置,仅仅是出于概念明晰目的而选择。在可供选择的实施例中,也可以使用任何其他合适的配置。处理器52通常在软件内编程,以执行本文所述功能。例如,该软件可以电子形式通过网络下载到处理器,或者作为另外一种选择或除此之外,该软件可以被提供和/或存储在非临时性有形介质(诸如磁性、光学或电子存储器)上。
系统20还包括能够标绘线圈32在区域30中产生的磁场的设备。在本发明的一个实施例中,标测器70用于标绘磁场,该标测器包括磁场检测器72的阵列,这些磁场检测器固定地安装在固体基部(诸如塑料板)的已知位置中。标测器70被配置用于使得其可以相对于定位垫26的已知预定位置和取向而设置在手术台24上。在一个实施例中,标测器70包括50个检测器。通常,检测器在标测器中被配置用于使得磁场和所有其不可忽略的梯度为可测量。在可供选择的实施例中,标测器具有78个分布在矩形框中的检测器,该矩形框的近似尺寸(高×宽×长)为等于150mm×250mm×250mm。
通常可使用来自检测器72的磁测量值,将标测器相对于上下左右大致居中地设置。
虽然在有关系统20操作的下列描述中假定标测器70用于标绘磁场,但应当理解标测器是用于测量区域30中的磁场的示例性系统,并且可以使用任何其他合适的系统,诸如可移动到已知位置中的一个或多个场检测器。此类可供选择的磁场标绘系统对于本领域内的普通技术人员而言将显而易见,并应认为涵盖于本发明的范围内。
检测器72可包括用于测量磁场大小和方向的任何便捷传感器,诸如大致类似于传感器46的霍尔探针或霍尔传感器。通常通过电缆74将来自检测器的读数传输到控制单元50,但也可使用任何其他方便的传输方法,诸如无线传输。
如下文更详细描述的,标测器70用于系统20的校准阶段,从而标测器及其连接电缆在图1中用虚线表示。在系统20处于其操作阶段时移除标测器及其电缆。
在一些实施例中,系统20包括一个或多个与传感器46大致类似的参考传感器76。如下文所述,传感器76控制单元50提供信号,以使该单元能够确定荧光检查仪60的位置和取向。通常,传感器76相对于手术台24固定,并便利地定位在手术台下方。作为另外一种选择,参考传感器76可固定到荧光检查仪60。
现在参见图4,该图为根据本发明实施例测定探针22远端21的位置的工序的流程图100。流程图100包括第一校准阶段102,然后是第二操作阶段104。该流程图的下列描述还涉及下列标题为“反应场模型”章节中的部分。
在校准阶段的第一标测步骤106中,将可扰动线圈32所产生的磁场的所有物体(本文称为场扰动器或场扰动元件)从区域30及其附近移除。此类物体包括荧光检查仪60。将标测器70以其相对于定位垫26的预定位置和取向设置在手术台24上,并启动发射器线圈32。控制单元50操控磁场检测器72,以便标测在区域30中未受扰动的磁场。
在第二标测步骤108中,在标测器70保持就位的情况下,将荧光检查仪60设置在区域30附近的相对于定位垫轴线的具体位置和取向(LO)。控制单元50可使用荧光检查仪控制器61测定具体位置和取向。作为另外一种选择,在包括参考传感器76的实施例中,控制单元50可使用来自传感器的信号测定位置和取向,如下文在标题为“荧光检查仪的适应定位”章节中所述。下文结合图5、图6A、图6B和图7描述了用于确定荧光检查仪已知位置和取向的另一可供选择的光学方法。
如下文所述,假定荧光检查仪60具有两个扰动元件:准直器62和检测器64。为了使下文针对步骤110描述的表征有效,选择了不同的位置和取向,以使得它们分成两组:第一组中准直器62远离区域30,以使得仅检测器64为扰动元件;第二组中检测器64远离区域30,以使得仅准直器62为扰动元件。
对于每一个已知位置和取向而言,控制单元50都操控磁场检测器72,以便标测区域30中的被扰动的磁场。
在多极表征步骤110中,控制单元假定由荧光检查仪60的每一个扰动元件引起的扰动都归因于相应的元件中磁图像源产生的反应场。在下文标题为“反应场模型”的章节中更详细描述了在步骤110中进行的分析。磁性源可视为发射器线圈32在扰动元件(即在两种元件中:准直器62和检测器64)中产生的图像。如果荧光检查仪60具有不同数量的元件,或者对于其他场扰动器而言,本领域内的普通技术人员将能够对该描述进行调整。
理论上,如果场扰动器为理想导体球并且发射器为力矩朝向球体中心的双极,那么球体产生的发射器图像为单个双极。为了容许存在与此类理论模型的差异(尤其包括:发射器不是理想双极,并且扰动器也不是理想导体球),本发明实施例假定在荧光检查仪的每一个扰动元件内,单个发射器线圈32都产生线圈的两幅或更多幅多极磁性源图像。
假定每一个多极都产生相应的反应场,并且多极可包括双极、四极和/或更高阶的极。考虑到反应场的图像源可能具有除双极之外的元件,控制单元50并非使用分立的极来分析数据,而是使用球谐函数分析。如章节“反应场模型”所阐释的,并且以举例的方式,本文假定每一个线圈32在准直器62中都产生五幅多极图像,并在检测器64中产生五幅多极图像。图像位置充当假定描述反应场的球谐函数的展开点。本领域内的普通技术人员将能够对本描述作出调整,以适于准直器和/或检测器中的其他数量的图像。
在步骤110结束时,通常为校准阶段结束时,控制单元已表征了每一个扰动元件所有的多个图像源。对于每一个扰动元件,该表征都包括定位图像源在扰动器中的位置,以及测定每一个扰动元件的图像源所产生的反应场矩阵[T反应]的元素。对于此处考虑的示例性实施例(其中准直器62和检测器64为扰动元件),测定了相应的矩阵[T反应]准直器和[T反应]检测器。
在一个实施例中,检测器具有尺寸为38cm×48cm的外部矩形框架。假定用于检测器的五个图像源位于距离矩形拐角大约10cm和位于矩形中心处,以及位于检测器表面下方大约10cm的位置。准直器的外部具有相似的矩形尺寸,并且准直器的五个图像源按照与检测器的那些图像源相似的方式设置。
在操作阶段,控制单元使用矩阵[T反应]准直器和[T反应]检测器来计算图像源产生的反应场。
在第一操作阶段的步骤112中,如果发射器线圈32在步骤110后已停用,则控制单元会启动发射器线圈32。
在扰动步骤114中,荧光检查仪60移动到已知位置和取向,即检测器位于区域30附近的已知位置和取向并且准直器也位于区域30附近的已知位置和取向,并且控制单元50配准已知位置和取向。可通过以上针对校准过程所述方法之一将已知位置和取向提供给控制单元50。
在初始位置测定步骤118中,控制单元测量从探针22远端21中的传感器产生的信号。作为第一逼近程度,控制单元假定由这些信号测定的场为未受扰动的场。控制单元使用得自第一标测步骤106的结果来测定远端的大致位置。
在场校正步骤120中,控制单元50使用方程(7)计算在前述步骤中测定的位置处的反应场。控制单元从测量的场中减去反应场,以得出在该位置处未受扰动的场的改进估值。
在改进的位置测定步骤122中,控制单元使用未受扰动的场的值和步骤120中测定的磁场的值,对远端位置作出改进的测定。
在可选的比较步骤124中,控制单元检查介于步骤122中测定的最后位置与倒数第二位置之间的差值是否小于预设值,预设值通常为大约0.1mm。如果差值大于或等于预设值,则控制单元返回到步骤122。如果差值小于预设值,则控制单元继续到最终步骤126。
应当理解,步骤120-124描述了一种可选的迭代循环,其中控制单元对远端位置执行逐渐精准的测定,同时还对该位置的未受扰动的场执行逐渐精准的测定。
在一些实施例中,对时序测量值执行大致类似于步骤120和124(但不包括步骤122)限定的循环。在这种情况下,可按时对每一个传感器的下一个测量值都实施反应场校正。在每16ms对测量值进行计算的一个此类实施例中,本发明人已经发现,对于固定传感器而言,需要进行三次或更少次迭代,从而得到传感器位置的真值。
在最终步骤126中,控制单元使用步骤122中得到的最终位置测定值作为远端的位置。
图5为根据本发明实施例的可选位置感测系统150的示意性立体图。图6A和图6B为根据本发明实施例的系统150的一部分的示意图和系统校准中使用的元件的示意图;除了下文所述的差异外,系统150的操作大致类似于系统20(图1-4)的操作,并且在系统20和150两者中用相同附图标号指示的元件通常在构造和操作方面大致类似。为简单明了起见,系统150的一些元件(存在于系统20中并在图1中示出)未在图5中示出。
在系统150中,一个或多个基准标记152固定到荧光检查仪60的元件。本文以举例的方式假定两个标记152附接到准直器62。此外,两个或更多个摄像头154固定地附接到定位垫26。摄像头154受控制单元50的控制,以使得处理器52能够接收并处理由摄像头生成的图像。摄像头固定到具有位置和取向的定位垫26,以使得由摄像头形成的图像包括基准标记152的图像。
如下文结合图7的流程图所详述,控制单元50使用基准标记152的摄像头图像来确定荧光检查仪60(包括其准直器62和检测器64)相对于定位垫26的位置和取向。为执行此测定操作,使摄像头的位置相对于定位垫配准。摄像头与定位垫的配准,以及基准标记的图像与荧光检查仪的位置和取向的关系的测定,均根据图7的流程图由控制单元50实现。
可使用装入定位垫中的已知位置中的可拆卸的预制校准夹具156来完成配准操作。在一些实施例中,夹具156可伸缩地铰接在铰链158和160处。通过使夹具156可伸缩,现场工程师可以简单有效地改进现有位置感测系统20中定位垫26中的夹具,如图6B所示。来自夹具156的元件的图像用于在定位垫中对准摄像头154,并用于校正可能由摄像头引起的任何光学误差。
控制单元50可使用荧光校准系统(FCS)162来确定基准标记的图像与荧光检查仪的位置和取向的关系,如图5和图6A所示。FCS 162包括第二对摄像头164,该对摄像头可以将其图像传输至控制单元。摄像头安装在支撑件166上,并且每一个摄像头都能够根据从控制单元50接收的命令独立地全景拍摄和倾斜。
图7为根据本发明实施例由操作系统150中控制单元50执行的步骤的流程图200。流程图200分为校准阶段202和操作阶段204。
在校准阶段的第一步骤206中,将基准标记152固定到荧光检查仪60。基准标记被布置为使得在如下文所述将摄像头连接到定位垫26时位于摄像头154的视野之内。通常,对将标记固定到其上的荧光检查仪的元件加以选择,以使得对于荧光检查仪的基本上所有的操作位置和取向,标记都位于摄像头154中的至少一个的视野之内。本文以举例的方式假定基准标记固定到准直器。在一些实施例中,在系统150中使用不止两个摄像头154,以使得对于荧光检查仪的所有位置和取向,基准标记152都位于摄像头中的至少一个的视野之内。
在摄像头安装步骤208中,摄像头154固定地附接到定位垫26,以使得摄像头位于已知的位置中并具有已知的取向。在一些实施例中,在附接之后,使用充当摄像头的校准物体的夹具156,将摄像头对齐到其已知位置和取向。由摄像头中的每一个形成的夹具的图像都可以使摄像头能够精确对准,以及提供使控制单元50能够校正摄像头成像中的缺陷的数据。作为另外一种选择,可采用用于将摄像头154相对于定位垫固定到已知位置和取向,以及用于校正摄像头的图像中的缺陷的任何其他简便方法。
在设置步骤210中,FCS 162相对于定位垫26固定地设置在已知位置中。FCS 162的定位可通过将FCS的支撑件166附接到手术台24(如图5所示),或附接到任何其他方便的物体(诸如定位垫26自身)来进行。设置FCS 162,使得安装在支撑件166上的摄像头164能够由控制单元旋转,以使得它们始终能够使基准标记152成像。
在数据测量步骤212中,使用荧光检查仪控制器61将荧光检查仪60移入多个不同位置中,即多个不同位置和取向。在每一个不同位置中,控制单元50都将信号传输至摄像头164,以使得摄像头位于已知全景拍摄和倾斜角度。控制单元50使摄像头全景拍摄和/或倾斜,以使得每一个摄像头164都使基准标记152成像。控制单元50记录每一个摄像头的全景拍摄和倾斜角度(由于控制单元发送全景拍摄和倾斜信号,因此它可感知每一个摄像头164的取向)以及由每一个摄像头形成的基准标记的图像。
此外,在荧光检查仪的每一个不同位置处,控制单元50都记录由每一个摄像头154生成的基准标记的图像。
在相关性步骤214中,控制单元50生成介于摄像头164的基准标记的图像和摄像头154的基准标记的图像之间的相关性。所述相关性允许控制单元配准两个系统的轴线,即荧光检查仪的轴线和定位垫的轴线。
在校准阶段的最后一个步骤216中,一旦进行了配准,FCS 162(包括其摄像头164)就被移除。
在操作阶段的操作步骤218中,启动摄像头154。在本实例中,在基准标记附接到准直器位置处,控制单元从摄像头接收基准标记的图像,并使用步骤214中确定的相关性对图像进行分析,以确定准直器相对于定位垫的位置和取向。由于荧光检查仪60的准直器和检测器彼此处于固定的物理关系,因此,控制单元可通过准直器的位置和取向确定检测器的位置和取向,具体由它们之间的距离而定。通常使用优化技术来确定该缺失值,该技术类似于下文所述的“荧光检查仪的适应定位”法。
已结合确定荧光检查仪相对于定位垫的位置和取向对流程图200进行了描述。应当理解,加以必要的修改后,也可以实施流程图中的各个步骤,来确定另一物体的位置和取向,所述物体通常为定位在定位垫附近并进行移动的、可能会干扰定位垫的发射器的物体。所有这类具体实施都被认为在本发明的范围内。
反应场模型
本发明实施例的反应场模型假定由设置在源磁场中的磁场扰动元件所导致的磁反应场是通过所述元件中的多个点源生成的。源磁场可由任何数量的场发射器生成。点源可视为场发射器的扰动元件中的图像。
为简单起见,对与本发明实施例有关的模型的描述进行了以下假设:
●扰动元件。荧光检查仪60包括通过C形臂66连接的准直器62和检测器64。本发明人已经发现,如果假定仅检测器和准直器为扰动元件,则会获得由荧光检查仪生成的反应场的良好结果,因此,下面的分析中假定仅有两个此类元件。然而,也可对其他扰动元件采用相同类型的分析,以使得本发明的范围包括任何数量的此类元件。
●发射器的数量和类型。为与构造28的三个三轴线圈相对应,假定有九个磁场发射器。发射器中的每一个都为近似的双极。然而,通常对场发射器的数量或构造没有任何限制。
●每一个扰动元件中的点源。以举例的方式,假定给定发射器在每一个扰动元件中成像时都生成五个点源,本文也称为展开点(因为下文提及的球谐函数是通过这些点展开的)。在其他实施例中,可以假定每个发射器的其他图像点源数量,并且用于不同扰动元件的源数量不必相等。在本发明的实施例中,假定与不同发射器对应的五个点源彼此全等。此外,假定该五个点源还位于矩形的中心处和拐角处。
●球谐函数。可将每一个点源都视为由双极、四极和/或更高阶的极组成。模型不是分析单独的分量,而是假定来自每一个点源的场都可通过球谐函数展开来表示。本发明人已经发现,展开至3阶会得到良好的结果,但应当理解,本发明的范围包括展开至更低阶或更高阶。
●位置测量。为简单起见,假定相对于单组轴线测量所有位置。以举例的方式,假定该单组轴线基于检测器64,其原点位于检测器的中心图像点源处,其中x轴和y轴由检测器中的点源的矩形的边限定。
一个发射器
首先考虑一个发射器以及由发射器在扰动元件(除非另外指明,否则在本文中假定为检测器64)中形成的一个展开点。为简单起见,假定所述一个展开点位于基于检测器的(x,y,z)坐标组的原点(0,0,0)处,那么离开原点的位置(x,y,z)处的磁反应场的表达式(表示为3元素列向量)为:
其中为表示为列向量的反应场,
[C0]为基于1-3阶球谐函数项、围绕原点展开(下文将给出[C]展开形式的表达式)的15元素列向量,并且
[T空间(x,y,z)]为得自[C0]的球谐函数项的3×15空间转移矩阵。
图8A-8G示出了根据本发明实施例的[T空间(x,y,z)]的元素中的每一个。
对于五个展开点(0,0,0)、(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)、(x4,y4,z4)而言,方程(1)变为:
其中[T空间(x,y,z)(5)]为3×75空间转移矩阵,并且
[Cs]为75元素列向量。
[T空间(x,y,z)(5)]通过将四个[T空间]矩阵(每一个矩阵距离原点的偏移量都为其相应位移的坐标)与[T空间(x,y,z)]连接在一起而形成。即,
[T空间(x,y,z)(5)]≡
[T空间(x,y,z)…T空间(x-x4,y-y4,z-z4)] (3)
[Cs]通过将5个列向量叠加在一起而形成,每一个列向量都基于围绕相应展开点展开的球谐函数项。即,
来自一个发射器的场
来自区域30中任意点处的任何给定发射器的场都可通过定位垫制造过程期间执行的校准步骤进行计算。为了表征发射器,考虑检测器64的区域中的多个点。本发明人已经发现,17个点已足够,但也可采用任何其他适合的数字。在17个点处的场可通过未受扰动的场的51元素列向量[B卫星]来表示,其中
根据以下方程,假定[Cs]随[B卫星]而变化:
[C5]=[T反应]·[B卫星] (6)
其中[T反应]为75×51反应矩阵,其中的元素表示检测器的反应特性。
一般情况
可通过上述推导,将用于在点(x,y,z)处的检测器反应磁场的表达式写为下面的方程(7)。该表达式假定检测器中有九个发射器和五个展开点。
[B反应(x,y,z)]=[T空间(x,y,z)(5)]·[T反应]·[B卫星(9)] (7)
其中[B反应(x,y,z)]为3×9反应场矩阵,该矩阵有9列,每一列都类似于方程(1)的反应场向量的列,
[B卫星(9)]为未受扰动的场的51×9矩阵,该矩阵具有类似于向量[B卫星](方程(5))的列,
并且对于方程(2)和(6)而言,[T空间(x,y,z)(5)]和[T反应]具有相同的量纲,而且元素符合方程(7)的要求。
应当理解,乘积[T空间(x,y,z)(5)]·[T反应]能有效地表征在检测器中形成的五幅图像。
方程(7)在系统20的校准阶段和操作阶段中均可使用,如下所述。
在一些实施例中,通过用[B反应(x,y,z)]和[T空间(x,y,z)(5)]除以||B(x,y,z)||来改写方程(7),而不是按原样使用方程(7),其中B(x,y,z)为来自在(x,y,z)处的发射器的未受扰动的场。本发明人已经发现,此类改写可为发射器场较弱的位置赋予更大的权重,并可以补偿该场的测量误差。
为简单起见,本文的具体实施方式假定按原样使用方程(7),并且本领域的普通技术人员加以必要的变通后将能够更改所述具体实施方式,以用于实施例实现对上述方程(7)的改写。
校准阶段
在校准阶段中,方程(7)用于获得[T反应]的值。
对于具体位置(x,y,z),[B反应(x,y,z)]的元素可通过如下方式取得:当检测器64位于对磁场有扰动的位置时和检测器被设定为不会扰动磁场时,获得这两个场之间的差值。
[T空间(5)(x,y,z)]的元素可通过图8中所示的方程取得,假定在适用时在所述方程中使用展开点(结合方程(2)提及)的位置。
因此,在方程(7)中,除[T反应]之外的所有项都是已知的。
可使用将为本领域的普通技术人员所熟知的方法,通过处理方程(7)取得[T反应]。一种用于取得[T反应]的此类方法是使[B反应(x,y,z)]和[T反应]向量化,并应用克罗内克积生成换算的方程:
其中上标T表示转置。
vec[T反应]通过将方程(10)的两边都乘以(的伪逆算子)取得。[T反应]通过将vec[T反应]分为合适的列长度而生成。
此类方法取得(x,y,z)的一个值和用于具体位置中的检测器(以及由此得到的展开点)的[T反应]。
在校准阶段中,假定标测器70使用N个传感器测量(x,y,z)的N个值。对于N个点(x1,y1,z1)、...、(xN,yN,zN)而言,可通过叠加方程的矩阵将方程(10)改写为:
方程(11)用于一个检测器位置,并可被改写为:
假定在校准阶段将检测器放入M个位置中,则可通过进一步叠加矩阵来改写方程(12),以形成方程(13):
就用于解出方程(10)的上述方法而言,可通过将方程的两边都乘以的伪逆算子而取得方程(13)中[T反应]的解。
在N=50且M=30的实施例中,矩阵具有40,000×4,000的近似量纲。如本领域所已知,可使用奇值分解来取得伪逆算子。在一个实施例中,通过如下方法使伪逆算子中的奇值数最小化:找到解的性质(即,解对扰动的校正程度)开始降级的点。通常,使用最多大约10-6的值。
通常执行校准阶段以取得用于每一个扰动元件的[T反应]。因此,对于其中假定仅准直器62和检测器64为扰动元件的荧光检查仪60而言,取得矩阵[T反应]准直器和矩阵[T反应]检测器。要取得上述矩阵[T反应]准直器,通常在更改荧光检查仪的位置的同时降低手术台24,以使得检测器的影响可以忽略不计。相似地,要取得上述矩阵[T反应]检测器,通常在更改荧光检查仪的位置的同时升高手术台24。
操作阶段
在操作阶段中,将在校准阶段中确定的[T反应]值代入方程(7)中,以取得[B反应(x,y,z)]的值。如果存在不止一个扰动元件,则假定反应场为来自每一个元件的[B反应(x,y,z)]的相应值的线性叠加。
荧光检查仪的适应定位
重新参考图1,本发明的实施例提供了用于相对于定位垫26的适应定位荧光检查仪60的方法。该方法假定通过任何方便的方法,完成基准荧光检查仪架与基准定位垫架的初步配准。此类方法通常使用荧光检查仪使手术台24上的物体成像,取得该物体离定位垫的线圈32的距离,并通过将图像和所述距离作比较来执行配准。该方法通常应用于荧光检查仪的多个不同位置,并使用已知的荧光检查仪参数,诸如准直器到检测器的距离,和荧光检查仪在使物体成像时的取向。通常在对系统20中的患者开始进行临床手术时进行初步配准。
在手术期间,通常在患者胸部上放置三个贴片电极,并且本文假定每一个贴片都包括类似于传感器76的传感器。虽然贴片传感器的位置通常相对于定位垫变化,但可用贴片传感器检测到的场来估算这些位置。所有五个传感器(每一个都具有三个线圈)均受到来自线圈32的九个传输的场以及来自准直器和检测器的图像场的作用,以使得在任意时间都有135个输入处理器52中的来自传感器的独立评估的测量信号(本文称为“Meas”)。应当理解,对于完整描述荧光检查仪及其元件的位置和取向所需的自由度数量而言,该信号数量提供一组高度冗余的值。此外,在传感器76处测量的未受扰动的场允许直接测量在传感器处的反应场。
因此,处理器52能够通过信号Meas来计算荧光检查仪的位置和取向以及检测器64和准直器62的位置。处理器计算以下代价函数:
其中M为由定位垫26生成的磁场与上述反应场模型的总和,
为传感器的位移,
为传感器的取向,
为荧光检查仪60的仰角,
θ荧光检查仪为荧光检查仪60的方位角,
d检测器为检测器64的距离,以及
d准直器为准直器62的距离,
并且其中所有上述变量均相对于定位垫的轴线和原点计算。
为了计算荧光检查仪及其元件的位置和取向,处理器优化代价函数C。该优化通常涉及计算C的一阶和更高阶导数,可使用用于本领域中已知的任何方法。
虽然上述具体实施方式提及的是荧光检查仪的元件为磁场扰动元件,但应当理解,本发明的范围也包括其他扰动元件。例如,探针22可包括扰动由构造28生成的磁场的金属元件,并且可以基本上类似于本文针对荧光检查仪元件所述的方式补偿这些元件所产生的扰动。此类金属元件可包括、但不限于手术台或其部件。
因此,应当理解,上述实施例是以举例的方式引用,并且本发明不限于上文所具体示出和描述的内容。更确切地说,本发明的范围包括上述各种特征的组合和子组合、以及本领域技术人员在阅读上述说明书时可能想到的并且现有技术中未公开的变型形式和修改形式。
Claims (6)
1.一种用于位置感测的方法,包括:
在定位垫上安装磁发射器,所述磁发射器被配置用于在患者体内产生磁场;
将定位垫摄像头以相应的固定取向附接到所述定位垫;
将可旋转摄像头连接到所述定位垫;
将基准标记附接到荧光检查仪,所述荧光检查仪被配置用于使所述患者成像;
将所述荧光检查仪定位到不同位置中;以及
对每一个位置都执行如下操作:
将所述可旋转摄像头取向为已知取向,用所述可旋转摄像头和所述定位垫摄像头形成所述基准标记的相应的图像,并且
分析所述相应的图像,以将所述荧光检查仪的位置和取向与所述定位垫的轴线配准。
2.根据权利要求1所述的方法,还包括从所述定位垫移除所述可旋转摄像头,以及仅使用由所述定位垫摄像头形成的所述基准标记的图像来确定所述荧光检查仪相对于所述定位垫轴线的所述位置和所述取向。
3.根据权利要求1所述的方法,其中将所述定位垫摄像头附接到所述定位垫的步骤包括将可拆卸夹具附接到所述定位垫,以及通过用所述摄像头使所述夹具成像而将所述定位垫摄像头与所述相应的固定取向对齐。
4.一种用于位置感测的设备,包括:
定位垫;
磁发射器,所述磁发射器安装在所述定位垫上,并且被配置用于在患者体内产生磁场;
定位垫摄像头,所述定位垫摄像头以相应的固定取向附接到所述定位垫;
可旋转摄像头,所述可旋转摄像头连接到所述定位垫;
荧光检查仪,所述荧光检查仪被配置用于使所述患者成像;
基准标记,所述基准标记附接到所述荧光检查仪;和
处理器,所述处理器被配置用于:
将所述荧光检查仪定位到不同位置中,并对每一个位置都执行如下操作:
将所述可旋转摄像头取向为已知取向,用所述可旋转摄像头和所述定位垫摄像头形成所述基准标记的相应的图像,以及
分析所述相应的图像,以将所述荧光检查仪的位置和取向与所述定位垫的轴线配准。
5.根据权利要求4所述的设备,其中所述处理器被配置用于当所述可旋转摄像头从所述定位垫移除时,仅使用所述定位垫摄像头形成的所述基准标记的图像来确定所述荧光检查仪相对于所述定位垫轴线的所述位置和所述取向。
6.根据权利要求4所述的设备,还包括可拆卸夹具,所述可拆卸夹具被配置用于附接到所述定位垫,并且其中所述定位垫摄像头被配置用于通过用所述摄像头使所述夹具成像而与所述相应的固定取向对齐。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/975,915 US8812079B2 (en) | 2010-12-22 | 2010-12-22 | Compensation for magnetic disturbance due to fluoroscope |
US12/975915 | 2010-12-22 | ||
CN201110461552.9A CN102525471B (zh) | 2010-12-22 | 2011-12-22 | 补偿由荧光检查仪产生的磁干扰 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110461552.9A Division CN102525471B (zh) | 2010-12-22 | 2011-12-22 | 补偿由荧光检查仪产生的磁干扰 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104665833A CN104665833A (zh) | 2015-06-03 |
CN104665833B true CN104665833B (zh) | 2017-03-29 |
Family
ID=45442953
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110461552.9A Active CN102525471B (zh) | 2010-12-22 | 2011-12-22 | 补偿由荧光检查仪产生的磁干扰 |
CN201510029457.XA Active CN104665833B (zh) | 2010-12-22 | 2011-12-22 | 补偿由荧光检查仪产生的磁干扰 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110461552.9A Active CN102525471B (zh) | 2010-12-22 | 2011-12-22 | 补偿由荧光检查仪产生的磁干扰 |
Country Status (8)
Country | Link |
---|---|
US (1) | US8812079B2 (zh) |
EP (3) | EP2468208B1 (zh) |
JP (2) | JP6000544B2 (zh) |
CN (2) | CN102525471B (zh) |
AU (1) | AU2011254069B2 (zh) |
CA (1) | CA2762262C (zh) |
ES (1) | ES2792111T3 (zh) |
IL (1) | IL216807B (zh) |
Families Citing this family (54)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8672837B2 (en) | 2010-06-24 | 2014-03-18 | Hansen Medical, Inc. | Methods and devices for controlling a shapeable medical device |
US9622714B2 (en) * | 2011-06-08 | 2017-04-18 | Siemens Medical Solutions Usa, Inc. | System and method for photographic determination of multichannel collimator channel pointing directions |
US8818486B2 (en) * | 2012-07-12 | 2014-08-26 | Biosense Webster (Israel) Ltd. | Position and orientation algorithm for a single axis sensor |
US9310446B2 (en) * | 2012-10-18 | 2016-04-12 | Analog Devices, Inc. | Magnetic field direction detector |
US10441236B2 (en) * | 2012-10-19 | 2019-10-15 | Biosense Webster (Israel) Ltd. | Integration between 3D maps and fluoroscopic images |
US9057600B2 (en) | 2013-03-13 | 2015-06-16 | Hansen Medical, Inc. | Reducing incremental measurement sensor error |
US9629595B2 (en) | 2013-03-15 | 2017-04-25 | Hansen Medical, Inc. | Systems and methods for localizing, tracking and/or controlling medical instruments |
US9271663B2 (en) | 2013-03-15 | 2016-03-01 | Hansen Medical, Inc. | Flexible instrument localization from both remote and elongation sensors |
US9014851B2 (en) | 2013-03-15 | 2015-04-21 | Hansen Medical, Inc. | Systems and methods for tracking robotically controlled medical instruments |
US11020016B2 (en) | 2013-05-30 | 2021-06-01 | Auris Health, Inc. | System and method for displaying anatomy and devices on a movable display |
US9186115B2 (en) * | 2013-07-01 | 2015-11-17 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for compensating for magnetic field during medical imaging |
US9696131B2 (en) * | 2013-12-24 | 2017-07-04 | Biosense Webster (Israel) Ltd. | Adaptive fluoroscope location for the application of field compensation |
US9638820B2 (en) * | 2014-03-03 | 2017-05-02 | Biosense Webster (Israel) Ltd. | Calibration jig for a flat location pad |
US20190004122A1 (en) * | 2014-04-25 | 2019-01-03 | Purdue Research Foundation | Wireless position sensing using magnetic field of single transmitter |
US10213133B2 (en) * | 2014-12-22 | 2019-02-26 | Biosense Webster (Israel) Ltd | Modeling of a magnetic field |
CA2919901A1 (en) | 2015-02-04 | 2016-08-04 | Hossein Sadjadi | Methods and apparatus for improved electromagnetic tracking and localization |
US10105117B2 (en) | 2015-02-13 | 2018-10-23 | Biosense Webster (Israel) Ltd. | Compensation for heart movement using coronary sinus catheter images |
US10307078B2 (en) | 2015-02-13 | 2019-06-04 | Biosense Webster (Israel) Ltd | Training of impedance based location system using registered catheter images |
JP7004574B2 (ja) * | 2015-06-12 | 2022-01-21 | コーニンクレッカ フィリップス エヌ ヴェ | 電磁装置の追跡 |
US9727963B2 (en) | 2015-09-18 | 2017-08-08 | Auris Surgical Robotics, Inc. | Navigation of tubular networks |
US10359458B2 (en) * | 2015-10-05 | 2019-07-23 | Covidien Lp | Systems and methods for automated mapping and accuracy-testing |
US10143526B2 (en) | 2015-11-30 | 2018-12-04 | Auris Health, Inc. | Robot-assisted driving systems and methods |
US11294003B2 (en) | 2016-03-23 | 2022-04-05 | Analog Devices International Unlimited Company | Magnetic field detector |
EP3432794B1 (en) * | 2016-05-03 | 2020-09-09 | St. Jude Medical International Holding S.à r.l. | Magnetic field distortion detection and correction in a magnetic localization system |
US10244926B2 (en) | 2016-12-28 | 2019-04-02 | Auris Health, Inc. | Detecting endolumenal buckling of flexible instruments |
JP7282685B2 (ja) | 2017-03-31 | 2023-05-29 | オーリス ヘルス インコーポレイテッド | 生理学的ノイズを補償する管腔ネットワークのナビゲーション用ロボットシステム |
US10022192B1 (en) | 2017-06-23 | 2018-07-17 | Auris Health, Inc. | Automatically-initialized robotic systems for navigation of luminal networks |
JP7330902B2 (ja) * | 2017-06-28 | 2023-08-22 | オーリス ヘルス インコーポレイテッド | 電磁歪み検出 |
CN110809452B (zh) | 2017-06-28 | 2023-05-23 | 奥瑞斯健康公司 | 电磁场发生器对准 |
US10739165B2 (en) | 2017-07-05 | 2020-08-11 | Analog Devices Global | Magnetic field sensor |
US10555778B2 (en) | 2017-10-13 | 2020-02-11 | Auris Health, Inc. | Image-based branch detection and mapping for navigation |
US11058493B2 (en) | 2017-10-13 | 2021-07-13 | Auris Health, Inc. | Robotic system configured for navigation path tracing |
EP3684562A4 (en) | 2017-12-14 | 2021-06-30 | Auris Health, Inc. | SYSTEM AND METHOD OF ESTIMATING THE LOCATION OF AN INSTRUMENT |
EP3684283A4 (en) | 2017-12-18 | 2021-07-14 | Auris Health, Inc. | METHODS AND SYSTEMS FOR MONITORING AND NAVIGATION OF INSTRUMENTS IN LUMINAL NETWORKS |
EP3773131A4 (en) | 2018-03-28 | 2021-12-15 | Auris Health, Inc. | LOCATION SENSOR RECORDING SYSTEMS AND METHODS |
US10827913B2 (en) | 2018-03-28 | 2020-11-10 | Auris Health, Inc. | Systems and methods for displaying estimated location of instrument |
US10905499B2 (en) | 2018-05-30 | 2021-02-02 | Auris Health, Inc. | Systems and methods for location sensor-based branch prediction |
EP3801189A4 (en) | 2018-05-31 | 2022-02-23 | Auris Health, Inc. | PATH-BASED NAVIGATION OF TUBULAR NETWORKS |
MX2020012897A (es) | 2018-05-31 | 2021-05-27 | Auris Health Inc | Sistemas roboticos y metodos para navegacion de la red luminal que detectan ruido fisiologico. |
WO2019232236A1 (en) | 2018-05-31 | 2019-12-05 | Auris Health, Inc. | Image-based airway analysis and mapping |
JP2022505955A (ja) * | 2018-11-18 | 2022-01-14 | トリグ メディカル リミテッド | 画像装置のための空間的位置合わせ方法 |
US20200237459A1 (en) | 2019-01-25 | 2020-07-30 | Biosense Webster (Israel) Ltd. | Flexible multi-coil tracking sensor |
US20200372409A1 (en) * | 2019-05-24 | 2020-11-26 | Boston Scientific Scimed Inc. | Electromagnetic distortion compensation for device tracking |
IT201900014559A1 (it) * | 2019-08-09 | 2021-02-09 | Mectron S P A | Assieme inserto con identificatore a radiofrequenza |
KR20220058918A (ko) | 2019-08-30 | 2022-05-10 | 아우리스 헬스, 인코포레이티드 | 기구 이미지 신뢰성 시스템 및 방법 |
CN114340542B (zh) | 2019-08-30 | 2023-07-21 | 奥瑞斯健康公司 | 用于位置传感器的基于权重的配准的系统和方法 |
WO2021044297A1 (en) | 2019-09-03 | 2021-03-11 | Auris Health, Inc. | Electromagnetic distortion detection and compensation |
CA3093980C (en) * | 2019-09-26 | 2023-10-31 | Ascension Technology Corporation | Reconfigurable transmitter array for electromagnetic tracking systems |
US11504023B2 (en) | 2019-12-16 | 2022-11-22 | Biosense Webster (Israel) Ltd. | Sparse calibration of magnetic field created by coils in metal-rich environment |
US11602372B2 (en) | 2019-12-31 | 2023-03-14 | Auris Health, Inc. | Alignment interfaces for percutaneous access |
JP2023508525A (ja) | 2019-12-31 | 2023-03-02 | オーリス ヘルス インコーポレイテッド | 経皮的アクセスのための位置合わせ技術 |
WO2021137072A1 (en) | 2019-12-31 | 2021-07-08 | Auris Health, Inc. | Anatomical feature identification and targeting |
CN112089491B (zh) * | 2020-08-18 | 2023-12-19 | 上海精劢医疗科技有限公司 | 电磁导航系统中磁场形变校准工具及方法 |
EP4019886B1 (en) * | 2020-12-28 | 2024-05-01 | Mitutoyo Corporation | Inductive position detection configuration for indicating a measurement device stylus position |
Family Cites Families (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5767669A (en) | 1996-06-14 | 1998-06-16 | Ascension Technology Corporation | Magnetic field position and orientation measurement system with dynamic eddy current rejection |
US6147480A (en) | 1997-10-23 | 2000-11-14 | Biosense, Inc. | Detection of metal disturbance |
US7366562B2 (en) * | 2003-10-17 | 2008-04-29 | Medtronic Navigation, Inc. | Method and apparatus for surgical navigation |
US6400139B1 (en) | 1999-11-01 | 2002-06-04 | Polhemus Inc. | Methods and apparatus for electromagnetic position and orientation tracking with distortion compensation |
US6484118B1 (en) | 2000-07-20 | 2002-11-19 | Biosense, Inc. | Electromagnetic position single axis system |
IL137520A (en) | 2000-07-25 | 2010-06-16 | Elbit Systems Ltd | Estimating position and orientation in electromagnetic systems |
CA2348135A1 (en) | 2001-05-17 | 2002-11-17 | Cedara Software Corp. | 3-d navigation for x-ray imaging system |
US6774624B2 (en) | 2002-03-27 | 2004-08-10 | Ge Medical Systems Global Technology Company, Llc | Magnetic tracking system |
US7065393B2 (en) | 2002-07-11 | 2006-06-20 | Cedara Software Corp. | Apparatus, system and method of calibrating medical imaging systems |
US6714008B1 (en) | 2002-07-29 | 2004-03-30 | The United States Of America As Represented By The Secretary Of The Navy | Gradiometric measurement methodology for determining magnetic fields of large objects |
CN100457034C (zh) * | 2002-12-11 | 2009-02-04 | 皇家飞利浦电子股份有限公司 | 具有校准装置的c形臂x射线装置 |
US7974680B2 (en) * | 2003-05-29 | 2011-07-05 | Biosense, Inc. | Hysteresis assessment for metal immunity |
CN1901835A (zh) * | 2003-11-14 | 2007-01-24 | 通用电气公司 | 减少电磁跟踪器中失真的系统和方法 |
US20050107687A1 (en) | 2003-11-14 | 2005-05-19 | Anderson Peter T. | System and method for distortion reduction in an electromagnetic tracker |
US7062391B2 (en) * | 2004-03-12 | 2006-06-13 | Vsm Medtech Systems Inc. | Motion compensation in biomagnetic measurement |
US20060229515A1 (en) * | 2004-11-17 | 2006-10-12 | The Regents Of The University Of California | Fiber optic evaluation of tissue modification |
US7657075B2 (en) | 2005-05-06 | 2010-02-02 | Stereotaxis, Inc. | Registration of three dimensional image data with X-ray imaging system |
DE102005023167B4 (de) | 2005-05-19 | 2008-01-03 | Siemens Ag | Verfahren und Vorrichtung zur Registrierung von 2D-Projektionsbildern relativ zu einem 3D-Bilddatensatz |
US8000772B2 (en) * | 2005-10-19 | 2011-08-16 | Biosense Webster, Inc. | Metal immunity in a reverse magnetic system |
US7532997B2 (en) | 2006-04-17 | 2009-05-12 | General Electric Company | Electromagnetic tracking using a discretized numerical field model |
US8040127B2 (en) * | 2006-08-15 | 2011-10-18 | General Electric Company | Multi-sensor distortion mapping method and system |
US10016148B2 (en) * | 2006-09-27 | 2018-07-10 | General Electric Company | Method and apparatus for correction of multiple EM sensor positions |
US7917337B2 (en) * | 2006-09-29 | 2011-03-29 | Optichron, Inc. | Adaptive self-linearization with separation filter |
US20080183064A1 (en) | 2007-01-30 | 2008-07-31 | General Electric Company | Multi-sensor distortion detection method and system |
US20080228065A1 (en) * | 2007-03-13 | 2008-09-18 | Viswanathan Raju R | System and Method for Registration of Localization and Imaging Systems for Navigational Control of Medical Devices |
DE102007019827A1 (de) * | 2007-04-26 | 2008-11-06 | Siemens Ag | System und Verfahren zur Ermittlung der Position eines Instruments |
WO2008154183A2 (en) | 2007-06-05 | 2008-12-18 | Ascension Technology Corporation | Systems and methods for compensating for large moving objects in magnetic-tracking environments |
US7995864B2 (en) | 2007-07-03 | 2011-08-09 | General Electric Company | Method and system for performing image registration |
US20090082665A1 (en) * | 2007-09-26 | 2009-03-26 | General Electric Company | System and method for tracking medical device |
FR2941861B1 (fr) * | 2009-02-12 | 2012-06-22 | Gen Electric | Procede de calibration de la position relative d'une table et d'un bras en c d'un systeme d'imagerie medical |
-
2010
- 2010-12-22 US US12/975,915 patent/US8812079B2/en active Active
-
2011
- 2011-12-07 IL IL216807A patent/IL216807B/en active IP Right Grant
- 2011-12-15 AU AU2011254069A patent/AU2011254069B2/en active Active
- 2011-12-15 CA CA2762262A patent/CA2762262C/en active Active
- 2011-12-21 JP JP2011279383A patent/JP6000544B2/ja active Active
- 2011-12-21 ES ES18177011T patent/ES2792111T3/es active Active
- 2011-12-21 EP EP11194809.7A patent/EP2468208B1/en active Active
- 2011-12-21 EP EP17176440.0A patent/EP3241495B1/en active Active
- 2011-12-21 EP EP18177011.6A patent/EP3403610B1/en active Active
- 2011-12-22 CN CN201110461552.9A patent/CN102525471B/zh active Active
- 2011-12-22 CN CN201510029457.XA patent/CN104665833B/zh active Active
-
2016
- 2016-08-31 JP JP2016168988A patent/JP6189509B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
EP3403610A1 (en) | 2018-11-21 |
EP3241495B1 (en) | 2018-06-13 |
AU2011254069B2 (en) | 2014-10-02 |
CN104665833A (zh) | 2015-06-03 |
IL216807A0 (en) | 2012-01-31 |
EP3403610B1 (en) | 2020-03-25 |
AU2011254069A1 (en) | 2012-07-12 |
EP3241495A1 (en) | 2017-11-08 |
EP2468208A3 (en) | 2012-09-05 |
CN102525471A (zh) | 2012-07-04 |
JP6189509B2 (ja) | 2017-08-30 |
US20120165656A1 (en) | 2012-06-28 |
CA2762262A1 (en) | 2012-06-22 |
CA2762262C (en) | 2019-07-23 |
ES2792111T3 (es) | 2020-11-10 |
EP2468208B1 (en) | 2017-10-25 |
EP2468208A2 (en) | 2012-06-27 |
IL216807B (en) | 2018-01-31 |
JP2012130703A (ja) | 2012-07-12 |
US8812079B2 (en) | 2014-08-19 |
JP2016195923A (ja) | 2016-11-24 |
CN102525471B (zh) | 2016-03-16 |
JP6000544B2 (ja) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104665833B (zh) | 补偿由荧光检查仪产生的磁干扰 | |
US20210022809A1 (en) | Robotic fluoroscopic navigation | |
CN103767683B (zh) | 3d标测图和荧光镜图像之间的整合 | |
US7097357B2 (en) | Method and system for improved correction of registration error in a fluoroscopic image | |
US10842461B2 (en) | Systems and methods of checking registrations for surgical systems | |
US20060025668A1 (en) | Operating table with embedded tracking technology | |
US8131342B2 (en) | Method and system for field mapping using integral methodology | |
JP2018110841A (ja) | 外科用システムのための位置合わせをチェックするシステム及び方法 | |
EP3295887B1 (en) | Robotic fluoroscopic navigation | |
US20240138927A1 (en) | Robotic fluoroscopic navigation | |
EP4289387A1 (en) | Surgical navigation system with flat panel registration fixture | |
US20240074819A1 (en) | Robotic fluoroscopic navigation | |
ES2679818T3 (es) | Compensación por perturbaciones magnéticas debido al fluoroscopio |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |