CN107918925A - 磁跟踪系统与成像装置的配准 - Google Patents
磁跟踪系统与成像装置的配准 Download PDFInfo
- Publication number
- CN107918925A CN107918925A CN201710941123.9A CN201710941123A CN107918925A CN 107918925 A CN107918925 A CN 107918925A CN 201710941123 A CN201710941123 A CN 201710941123A CN 107918925 A CN107918925 A CN 107918925A
- Authority
- CN
- China
- Prior art keywords
- point
- voxel
- probe
- processor
- point set
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 16
- 239000000523 sample Substances 0.000 claims abstract description 67
- 238000007689 inspection Methods 0.000 claims abstract description 36
- 238000006243 chemical reaction Methods 0.000 claims abstract description 33
- 230000003834 intracellular effect Effects 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims description 29
- 238000013459 approach Methods 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 2
- 239000013598 vector Substances 0.000 description 21
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 description 9
- 210000001331 nose Anatomy 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000001815 facial effect Effects 0.000 description 5
- 210000003128 head Anatomy 0.000 description 5
- 210000003695 paranasal sinus Anatomy 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000006698 induction Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 241001442234 Cosa Species 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 230000002688 persistence Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明题为“磁跟踪系统与成像装置的配准”。本发明公开了描述的实施方案包括处理器,该处理器在探头被放置在受检者的身体上的相应位置处时接收来自所述探头的信号,并且从所述信号导出与所述相应位置对应的点集。然后该处理器通过以下方式将点集与身体的三维图像配准:针对所述点中的每个点,识别所述图像的一个或多个体素,所述一个或多个体素的相应表面法线被取向在所述探头保持在所述点对应的所述位置处的取向的给定角度内;以及计算将所述点中的每个点与针对所述点所识别的所述体素中选择的一个体素对齐的转换。随后,该处理器基于所述配准在用于获取所述图像的成像装置的坐标系中跟踪胞内抗体工具的位置。
Description
技术领域
本发明涉及医疗装置以及相关联的软件的领域,并且具体地讲,涉及医学规程的背景内的坐标系的自动配准。
背景技术
迭代最近点(ICP)算法通常被用来将两个点集或点“云”彼此配准。使用ICP,“源云”被转换成与“参考云”最佳匹配,使得两朵云之间的距离的量度最小化。在其最基本的具体实施中,ICP算法迭代地执行以下一系列步骤直至针对收敛的某些条件得到满足:
1)针对源云内的每个点,发现参考云中对应的点(其可被称为“最近邻”)。
2)计算使两个点云之间的总距离最小化的转换,该总距离为源点与对应参考点之间的相应距离(例如,计算为平方欧氏距离)的总和。
3)使用所计算出的转换来转换源云。
以引用方式并入本文的Rusu、Radu Bogdan等人的“Aligning point cloud viewsusing persistent feature histograms”(2008IEEE/RSJ智能机器人与系统国际会议,IEEE,2008)针对将点云数据视图对齐到一致的全球模型中的问题而调查永久点特征直方图的用途。给定噪点云的集合,所呈现的算法估计局部描述每个点的几何结构的一组稳固的16D特征。通过分析在不同标度下的特征的持久性,该算法提取最好地表征给定点云的最佳集合。所得的持续性特征被用于初始对齐算法来估计大约将输入数据集配准的刚性转换。该算法通过将数据集转换成其收敛槽来为迭代配准算法诸如ICP(迭代最近点)提供好的起始点。
发明内容
根据本发明的一些实施方案,提供了包括电接口和处理器的一种系统。当探头被放置在受检者的身体上的相应位置处时,该处理器被配置成经由该电接口接收来自该探头的多个信号。该处理器还被配置成由所述信号导出与所述相应位置对应的点集。该处理器还被配置成通过以下方式将点集与受检者的身体的三维图像配准:针对所述点中的每个点,识别所述图像的一个或多个体素,所述一个或多个体素的相应表面法线被取向在所述探头保持在所述点对应的所述位置处的取向的给定角度内;以及计算将所述点中的每个点与针对所述点识别的所述体素中选择的一个体素对齐的转换。随后,基于所述配准,该处理器还被配置成在用于获取所述图像的成像装置的坐标系中跟踪胞内抗体工具的位置。
在一些实施方案中,所述处理器被配置成初始以第一精度水平来将所述点集与所述图像对齐,并且所计算出的转换以大于所述第一精度水平的第二精度水平将所述点集与所述图像对齐。
在一些实施方案中,所述处理器被配置成通过迭代地执行所述识别和所述计算来将所述点集与所述图像配准。
在一些实施方案中,所述处理器被配置成通过使用迭代最近点(ICP)算法来迭代地执行所述识别和所述计算。
在一些实施方案中,该给定角度小于50度。
在一些实施方案中,所述转换使所述点集与所述体素中相应的选择的体素之间的距离度量最小化。
在一些实施方案中,所述距离度量基于(i)所述体素中所述选择的体素的相应的表面法线与(ii)所述探头保持的相应的取向之间的相应角度。
在一些实施方案中,所述处理器还被配置成从所述信号导出所述探头保持在所述位置中的每个位置处的取向。
在一些实施方案中,所述处理器还被配置成在所述图像上在所述图像中的与所跟踪到的位置对应的位置处覆盖所述胞内抗体工具的视觉表示。
根据本发明的一些实施方案,还提供了一种方法,该方法包括将点集与受检者的身体的三维图像配准,所述点集与由跟踪系统跟踪的探头所放置的所述身体上的相应位置对应。然后通过以下方式执行该配准:针对所述点中的每个点,识别所述图像的一个或多个体素,所述一个或多个体素的相应表面法线被取向在所述探头保持在所述点对应的所述位置处的取向的给定角度内;以及计算将所述点中的每个点与针对所述点识别的所述体素中选择的一个体素对齐的转换。随后,该方法还包括使用所述跟踪系统并基于所述配准在用于获取所述图像的成像装置的坐标系中跟踪胞内抗体工具的位置。
在一些实施方案中,所述点与所述受检者的面部上的相应位置对应。
在一些实施方案中,跟踪系统包括磁跟踪系统。
结合附图,从本发明实施方案的以下详细描述,将更全面地理解本发明,其中:
附图说明
图1为根据本发明的一些实施方案的磁跟踪系统的示意图;
图2为根据本发明的一些实施方案的由处理器执行的配准技术的流程图;并且
图3为根据本发明的一些实施方案的识别候选对应体素的示意图。
具体实施方式
概述
本发明的实施方案提供了用于在医学规程之前将受检者的三维图像(或“容积”)与用于在该规程期间跟踪操作工具的跟踪系统配准的改进的方法和系统。这允许操作医师在该规程期间探知该操作工具相对于受检者的解剖结构的位置。例如,基于该配准,该工具的表示可覆盖在该规程期间实时获取的图像上,于该图像中与该工具的当前位置对应的位置处。
在一些实施方案中,例如,在对受试者执行鼻窦手术之前,使用计算机断层扫描(CT)获取受试者头部的容积。然后将所获取的CT容积与磁跟踪系统配准,该磁跟踪系统用于跟踪将被插入受检者的鼻窦中的操作工具。为了执行该配准,由磁跟踪系统跟踪的探头在获取CT容积之前或之后穿过受检者的面部。在面部上的各个位置处,记录在跟踪系统的坐标系中被表示为点的探头的位置。然后通过将从探头获得的点集与CT容积配准来将磁跟踪系统与CT容积配准。
为了执行这种配准,可使用ICP算法。例如,探头坐标集可被指定为源云,使得ICP算法迭代地(i)将每个探头坐标与对应的最近邻CT容积体素配对,并且然后(ii)将使对应点对之间的总距离最小化的转换应用到探头坐标。经过多次迭代计算的转换的合成限定两个坐标系之间的配准。
使用ICP算法的挑战是有时该算法由于算法的一个或多个迭代处较差的最近邻选择而缓慢收敛。
本发明的实施方案通过使用从探头获得的表面几何信息来更好地执行最近邻选择而解决了这个问题。具体地讲,随着探头坐标的获取,探头被保持成使得探头被配准在法向量对面部的合适角度阈值内。(换句话说,该探头被保持为大约垂直于面部。)因此,通过记录探头在探头坐标中的每一个探头坐标处的取向,可获得法向量在探头坐标处的大约相应取向。然后,在迭代中的每一次迭代处,探头坐标中的任一个探头坐标的候选最近邻可限于法向量在探头取向的给定角度内的那些体素。因此,可过滤掉不适当的最近邻候选,从而加速算法的收敛。
例如,在不使用本文所述的技术的情况下,由于CT容量体素的接近度,紧邻受检者的鼻部右侧的探头坐标可与位于受检者的鼻部右侧的底部处的该体素配对。然而,使用本文所述的技术,该CT容量体素将从作为候选最近邻的考虑中去除,由于该CT容量体素处的法线大约垂直于探头坐标处的法线。
虽然本公开主要讨论了上述具体应用,即,在其中磁跟踪系统与CT容量配准的鼻窦手术,需要注意的是,本文所述的技术可被应用于在任何种类的手术的背景下加速配准,而不管使用的特定成像模式或跟踪系统。一般来讲,从探头获取的点可对应于受检者的面部上的相应位置,和/或受检者的身体的任何其他部分上的相应位置。此外,本文所述的技术可应用于加速任何种类的点集配准算法,其中上述点对应挑战呈现其自己,包括例如鲁棒点匹配和核相关。
系统说明
首先参见图1,该图为根据本发明一些实施方案的磁跟踪系统20的示意图。如以下进一步描述的,系统20可用于促进鼻窦手术的执行,其中胞内抗体操作工具被插入受检者28的鼻窦中。
系统20包括一个或多个磁场发生器32。在手术期间,发生器32产生相应的磁场,其在联接到操作工具的一个或多个磁传感器中诱导相应电压。系统20还包括处理器34,其被配置成控制磁场发生器32、基于磁传感器中诱导的电压探知工具的位置和取向,以及将跟踪系统的坐标系与其他坐标系配准,如下文更详细所述。
在鼻窦手术之前,受检者28通常躺在手术台30上,例如以仰卧位,使得受检者的头部靠近磁场发生器。随后,在磁场发生器产生相应磁场时,医师26用其远侧端部包括一个或多个(例如,三个)磁传感器24的探头22的远侧尖端37来接触受检者的面部的多个位置。在探头22和受检者皮肤之间的每个接触点处,远侧尖端37的位置和探头的取向由处理器34基于传感器中诱导的电压来探知。(换句话说,磁跟踪系统用于探知探头所保持在的位置和取向。)因此医师获取在跟踪系统的三维坐标系(其还可被称为探头的三维坐标系)中与受检者的面部上探头所放置的相应位置对应的点集44。通常,如上所述,医师使探头保持大约垂直于受检者的面部,同时获取点集44。
如图1所示,可以围绕受检者的头部的马蹄构型将发生器32布置在手术台30的表面上。另选地,发生器可位于手术台旁边或下方,和/或可被布置成任何其他合适的构型。
如图1所示,系统20通常包括具有一个或多个电接口38的控制台36。探头22和发生器32经由电接口38连接至处理器34。因此,随着探头被放置在受检者的面部上的相应位置处,经由电接口38,处理器可控制磁场发生器,以及从探头接收指示探头位置和取向的信号(例如,诱导的电压)。然后,处理器可从该信号导出点集44,以及相关联的探头取向。(为简单起见,可以这样说,处理器从探头接收点集44以及相关联的探头取向。)
在获取到点集之后,处理器34将该点集与受检者的头部的三维图像40配准,该图像例如可通过CT成像装置获取。通常,处理器首先以第一粗精度水平将点集与图像配准。例如,可通过将与明确位置例如受检者的鼻部的尖端对应的若干点与图像中的对应体素配准来执行初始配准。通常,医师手动地识别这些点和体素,之后处理器计算将所识别的点映射到所识别的体素的转换,并将该转换应用到点集。
接下来,处理器通过计算如下转换来精修配准,该转换以大于第一精度水平的第二精度水平将点集与容量配准。通常,处理器执行迭代求精,从而在每个迭代期间,处理器计算相继更精确的转换。例如,处理器可使用如上所述执行迭代配准的ICP算法。在一些实施方案中,处理器34在配准程序期间在显示器42上显示图像40和点44。
一般来讲,处理器34可体现为单个处理器或协作式联网或集群的处理器集。处理器34通常是编程化数字计算装置,该装置包括中央处理单元(CPU)、随机存取存储器(RAM)、非易失辅助存储器诸如硬盘驱动器或光盘驱动器、网络接口和/或外围装置。如本领域所公知的,将包括了软件程序的程序代码和/或数据加载到RAM中用于由CPU执行和处理,并且生成结果用于显示、输出、传输或存储。程序代码和/或数据可以电子形式例如经网络下载到计算机,或者,另选地或除此之外,程序代码和/或数据可提供和/或存储在诸如磁性存储器、光学存储器或电子存储器的非临时性有形介质上。此类程序代码和/或数据在提供给处理器之后,产生被配置成执行本文所述任务的机器或专用计算机。
另外现在参见图2,该图是根据本发明的一些实施方案的由处理器34执行的配准技术的流程图。
如上所述,处理器34首先在初始配准步骤46处执行点集与图像的初始配准。换句话讲,如上所述,处理器计算将点与图像中适当的对应体素粗略配准的转换,并将所计算出的转换应用到点。
接下来,在识别步骤48处,针对点集中的每个点,处理器识别一个或多个候选对应体素51,候选对应体素中的每一个体素将可能被选择为该点的对应体素,即最近邻体素。通常,处理器仅识别作为候选体素的位于点的给定距离内(在其当前转换位置处)并且其相应表面法线被取向在所述探头保持在所述点对应的所述位置处的取向的给定角度内的那些体素。(施加到点集的任何转换也施加到探头的取向向量,使得取向向量与表面法线之间的角度可被正确地计算。)如下文参考图3进一步描述,通过基于这些体素处的表面法线取向选择点的候选体素,处理器可避免选择不适当的对应体素。(为了探知表面法线,处理器可使用任何合适的算法,诸如在函数“surfnorm.”中所施用的算法)
随后,在转换计算步骤50处,处理器计算将每个点与该点的候选对应体素中的选择的一个体素配准的转换。在计算该转换时,处理器通常使点集与所选择的对应体素之间的距离度量最小化。例如,在ICP算法的一个具体实施中,假设有N个点,最小化的距离D为||TP-Q||2,其中(i)P为在每列中包含点中的相应一个点的齐次坐标的4xN矩阵,(ii)Q为在每列中包含对应体素的坐标的3xN矩阵,并且(iii)T为施加到P的3x4转换矩阵。换句话讲,处理器找到T和Q,其中Q中的每个体素均选自相关的坐标集,其使距离D最小化。(需要注意的是,任何其他合适的距离度量诸如“L1”或“L∞”距离可被最小化,另选地最小化至上面提到的“L2”距离。)
随后,在转换步骤52处,处理器将所计算出的转换施加到所述点。例如,在第i次迭代期间,假设发现转换矩阵Ti使点集与对应体素之间的距离最小化,则处理器可计算Pi+1=TiPi。然后,在检查步骤54处,处理器检查收敛是否已经达到。例如,处理器可检查距离D是否低于阈值;另选地,处理器可检查所计算出的转换是否已经将点移动了小于阈值距离。如果算法尚未收敛,则处理器执行识别步骤48、计算步骤50和转换步骤52的另一迭代。否则,点集与图像的配准结束。随后,如果在先前的步骤中尚未执行,则处理器计算将磁跟踪系统的坐标系与成像装置的坐标系配准的复合转换。例如,假设针对初始配准使用的初始转换T0以及在配准算法的k次迭代期间计算的转换T1…Tk,则复合转换Tc将为T0×T1×…Tk。
随后,在受检者手术期间,处理器从联接到胞内抗体工具的磁传感器接收诱导的电压。根据这些诱导的电压,处理器导出工具在磁跟踪系统的坐标系中的坐标。然后,处理器将复合转换施加到这些坐标,以便将这些坐标映射到用于获取图像的成像装置的坐标系。因此,处理器跟踪胞内抗体工具在成像装置的坐标系中的位置。换句话说,使用跟踪系统并基于该配准,在成像装置的坐标系中跟踪胞内抗体工具的位置;跟踪系统跟踪该工具在该跟踪系统的坐标系中的位置,并然后在配准期间计算的转换将该位置映射到成像装置的坐标系中。
通常,在手术期间,处理器进一步在由成像装置获取的图像上在所述图像中的与胞内抗体工具的所跟踪到的位置对应的位置处覆盖该工具的视觉表示。胞内抗体工具的表示可覆盖在图像40上或由成像装置获取的另一个图像上,诸如在手术期间实时获取的图像上。
现在参见图3,其为根据本发明的一些实施方案的识别候选对应体素的示意图。
在左侧,图3示出了受检者的面部的一部分。探头22接触紧邻受检者的鼻部的受检者的面部。由沿探头的中心纵向轴线延伸的取向向量58来指示探头的取向。通常,如上面所指出的,医师握住探头使得探头大约垂直于受检者的面部;因此,取向向量58大约与受检者的面部的探头接触该面部的位置60处的法向量配准。
在右侧,图3示出了受检者的面部的相同部分的三维图像40,该图像40包括多个体素62。如上所述,在识别步骤48中,处理器识别从探头获取的点中的每个点的候选对应体素。图3示出了在探头的坐标系中与位置60对应的特定点44的该过程。(假设一个或多个转换已被执行到点44,使得该点在该图像中处于粗略正确的位置。)
如上所述,在识别步骤48期间的挑战不是去识别作为候选对应体素的具有很小的可能性与点对应的那些体素。例如,图3指定靠近点44的两个特定体素:第一体素64和第二体素70。在第二体素70潜在地与点44适当匹配时,位于受检者的鼻部上的第一体素64明显为不适当匹配;因此,将第一体素64选择为点44的对应体素将有可能导致配准算法相对于选择第二体素70来说更加缓慢地收敛。
如上所述,处理器通过基于表面法线信息来执行候选对应体素的选择识别而解决了这一难题。例如,在识别点44的候选对应体素时,处理器可考虑第一体素64的法向量66和第二体素70的法向量68。具体地讲,处理器可以计算这些法向量与取向向量58之间的相应角度θ。在法向量68的情况下,该角度相对较小;另一方面,在法向量66的情况下,该角度相对较大。因此,处理器可将第二体素70而不是第一体素64识别为点44的候选对应体素。
假设医师握住探头使其大约垂直于受检者的面部,针对适当的候选对应体素,θ可能接近零。然而,处理器可为θ设定为探头的取向留有一定余地的阈值。例如,所述阈值可具有小于50度的任何合适的值,诸如45度。此类阈值为探头的取向留有一定余地,同时仍过滤掉不适当的候选体素。
上述实施方案将表面法线信息用作过滤器,以选择性地识别候选对应体素。在识别这些候选之后,例如基于通常距离最小化的目标来执行针对每个点的相应对应体素(即,最近邻)的选择,而不进一步考虑该表面法线信息。换句话说,上述实施方案使用识别步骤48而不是计算步骤50中的表面法线信息。
在其他实施方案中,表面法线信息另选地或除此之外在计算步骤50中使用,用于对应体素的选择和转换的计算。例如,在计算步骤50中最小化的距离可基于所选择的对应体素的相应表面法线与探头所保持的相应取向之间的相应角度。例如,距离度量可采取||TPi-Q||2+αQ(θ1,θ2,...θN)的形式,其中α为适当的加权系数,并且Q为将角度θ映射到标量值的“罚”函数,从而对于较高的θ值,标量值增大。(每个角度θk为探头针对第k个点的取向向量与第k个点的对应体素处的法向量之间的角度。)
示例伪码
在该示例性实施方案中,当以任何合适的编程语言实施时,以下伪码可用于在识别步骤48中识别针对点集中的点中的每一个的候选对应体素。以下实施例包括若干评论,每个都开始于双斜杠(“//”),其解释了伪码的功能。
//初始化将保持所有候选的向量
向量<vec4>skinSurfaceNeighbors;
//限定45度(45度为角度阈值)的余弦
实型常量COS45=0.707f;
//在所有小区近邻(在特定距离内的相邻体素)上循环
针对(inti=0;i<cellNeigbors.size();i++){
//得到体素
vec4pos=m_dAsset->GetVolumeCoords(cellNeigbors[i]);
//得到体素的法向量
vec4n=m_dAsset->GetNormal(pos);
n.Normalize();
//得到探头的取向向量
vec4normPointDir=cathPointingDir.Normalize();
//计算法向量与取向向量之间的角度的余弦
实型cosA=normPointDir.Dot(n).GetElement<0>();
//如果该体素位于皮肤表面上并且角度小于45度,则将该体素列为候选
如果(isSkinSurfaceCell(cellNeigbors[i])&&(cosA>=COS45))
skinSurfaceNeighbors.push_back(pos);
}
本领域技术人员应当理解,本发明不限于上文具体示出和描述的内容。相反,本发明实施方案的范围包括上文所述各种特征的组合与子组合两者,以及本领域的技术人员在阅读上述说明书时可能想到的未在现有技术范围内的变型和修改。以引用方式并入本专利申请的文献被视为本专利申请的整体部分,不同的是如果这些并入的文献中定义的任何术语与本说明书中明确或隐含地给出的定义相冲突,则应只考虑本说明书中的定义。
Claims (20)
1.一种系统,包括:
电接口;和
处理器,所述处理器被配置成:
当探头被放置在受检者的身体上的相应位置处时经由所述电接口接收来自所述探头的多个信号,
从所述信号导出与所述相应位置对应的点集,
通过以下方式来将所述点集与所述受检者的所述身体的三维图像配准:
针对所述点中的每个点,识别所述图像的一个或多个体素,所述一个或多个体素的相应表面法线被取向在所述探头保持在所述点对应的所述位置处的取向的给定角度内,以及
计算将所述点中的每个点与针对所述点所识别的所述体素中选择的一个体素对齐的转换,并且
随后,基于所述配准,在用于获取所述图像的成像装置的坐标系中跟踪胞内抗体工具的位置。
2.根据权利要求1所述的系统,其中所述处理器被配置成初始以第一精度水平来将所述点集与所述图像对齐,并且其中所计算出的转换以大于所述第一精度水平的第二精度水平将所述点集与所述图像对齐。
3.根据权利要求1所述的系统,其中所述处理器被配置成通过迭代地执行所述识别和所述计算来将所述点集与所述图像配准。
4.根据权利要求3所述的系统,其中所述处理器被配置成通过使用迭代最近点(ICP)算法来迭代地执行所述识别和所述计算。
5.根据权利要求1所述的系统,其中所述给定角度小于50度。
6.根据权利要求1所述的系统,其中所述转换使所述点集与所述体素中相应的选择的体素之间的距离度量最小化。
7.根据权利要求6所述的系统,其中所述距离度量基于(i)所述体素中所述选择的体素的相应的表面法线与(ii)所述探头保持的相应的取向之间的相应角度。
8.根据权利要求1所述的系统,其中所述处理器还被配置成从所述信号导出所述探头保持在所述位置中的每个位置处的取向。
9.根据权利要求1所述的系统,其中所述处理器还被配置成在所述图像上在所述图像中的与所跟踪到的位置对应的位置处覆盖所述胞内抗体工具的视觉表示。
10.一种方法,包括:
通过以下方式将点集与受检者的身体的三维图像配准,所述点集与由跟踪系统跟踪的探头所放置的所述身体上的相应位置对应:
针对所述点中的每个点,识别所述图像的一个或多个体素,所述一个或多个体素的相应表面法线被取向在所述探头保持在所述点对应的所述位置处的取向的给定角度内,以及
计算将所述点中的每个点与针对所述点所识别的所述体素中选择的一个体素对齐的转换;以及
随后,使用所述跟踪系统并基于所述配准,在用于获取所述图像的成像装置的坐标系中跟踪胞内抗体工具的位置。
11.根据权利要求10所述的方法,其中将所述点集与所述图像配准还包括初始以第一精度水平来将所述点集与所述图像对齐,并且其中所计算出的转换以大于所述第一精度水平的第二精度水平将所述点集与所述图像对齐。
12.根据权利要求10所述的方法,其中将所述点集与所述图像配准包括迭代地执行所述识别和所述计算。
13.根据权利要求12所述的方法,其中迭代地执行所述识别和所述计算包括使用迭代最近点(ICP)算法。
14.根据权利要求10所述的方法,其中所述给定角度小于50度。
15.根据权利要求10所述的方法,其中所述点与所述受检者的面部上的相应位置对应。
16.根据权利要求10所述的方法,其中所述转换使所述点集与所述体素中相应的选择的体素之间的距离度量最小化。
17.根据权利要求16所述的方法,其中所述距离度量基于(i)所述体素中所述选择的体素的相应的表面法线与(ii)所述探头保持的相应的取向之间的相应角度。
18.根据权利要求10所述的方法,其中所述跟踪系统包括磁跟踪系统。
19.根据权利要求10所述的方法,还包括使用所述跟踪系统探知所述探头保持在所述位置中的每一个位置处的所述取向。
20.根据权利要求10所述的方法,还包括在所述图像上在所述图像中的与所跟踪到的位置对应的位置处覆盖所述胞内抗体工具的视觉表示。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/290968 | 2016-10-11 | ||
US15/290,968 US10152786B2 (en) | 2016-10-11 | 2016-10-11 | Registration of a magnetic tracking system with an imaging device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107918925A true CN107918925A (zh) | 2018-04-17 |
CN107918925B CN107918925B (zh) | 2023-05-02 |
Family
ID=60201806
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710941123.9A Active CN107918925B (zh) | 2016-10-11 | 2017-10-11 | 磁跟踪系统与成像装置的配准 |
Country Status (7)
Country | Link |
---|---|
US (2) | US10152786B2 (zh) |
EP (1) | EP3309749B1 (zh) |
JP (1) | JP7046553B2 (zh) |
CN (1) | CN107918925B (zh) |
AU (1) | AU2017235896A1 (zh) |
CA (1) | CA2981027A1 (zh) |
IL (1) | IL254681B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111281539A (zh) * | 2018-12-06 | 2020-06-16 | 韦伯斯特生物官能(以色列)有限公司 | 用于执行面部配准的方法和设备 |
CN111759463A (zh) * | 2020-07-31 | 2020-10-13 | 南京普爱医疗设备股份有限公司 | 一种提高手术机械臂定位精度的方法 |
CN113660912A (zh) * | 2019-04-10 | 2021-11-16 | 博医来股份公司 | 受试者的相关表面点的采样方法 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11000336B2 (en) * | 2016-09-23 | 2021-05-11 | Koninklijke Philips N.V. | Visualization of an image object relating to an instrucment in an extracorporeal image |
US10152786B2 (en) | 2016-10-11 | 2018-12-11 | Biosense Webster (Israel) Ltd. | Registration of a magnetic tracking system with an imaging device |
US11653853B2 (en) | 2016-11-29 | 2023-05-23 | Biosense Webster (Israel) Ltd. | Visualization of distances to walls of anatomical cavities |
US10510171B2 (en) | 2016-11-29 | 2019-12-17 | Biosense Webster (Israel) Ltd. | Visualization of anatomical cavities |
US10452263B2 (en) * | 2017-09-13 | 2019-10-22 | Biosense Webster (Israel) Ltd. | Patient face as touchpad user interface |
CN116725679A (zh) * | 2022-08-12 | 2023-09-12 | 北京和华瑞博医疗科技有限公司 | 配准点确定以及配准方法、装置、设备、介质和程序产品 |
CN116993791B (zh) * | 2023-09-28 | 2023-12-19 | 真健康(北京)医疗科技有限公司 | 一种基于体表定位装置的医学图像配准方法及设备 |
CN117257459B (zh) * | 2023-11-22 | 2024-03-12 | 杭州先奥科技有限公司 | 一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统 |
Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050054900A1 (en) * | 2003-07-21 | 2005-03-10 | Vanderbilt University | Ophthalmic orbital surgery apparatus and method and image-guided navigation system |
CN101019765A (zh) * | 2007-03-29 | 2007-08-22 | 新奥博为技术有限公司 | 一种磁共振图像引导下的手术系统及手术导航方法 |
CN101051387A (zh) * | 2006-03-21 | 2007-10-10 | 韦伯斯特生物官能公司 | 使用局部加权拟合的图像配准 |
CN101133431A (zh) * | 2005-02-03 | 2008-02-27 | 布拉科成像S.P.A.公司 | 能够减少物体运动造成的成像伪影的生物医学图像配准方法和计算机程序产品 |
CN101229080A (zh) * | 2006-10-09 | 2008-07-30 | 韦伯斯特生物官能公司 | 采用器官外部解剖特征的器官图像的配准 |
CN101297321A (zh) * | 2005-10-25 | 2008-10-29 | 布拉科成像S.P.A.公司 | 配准图像的方法,用于执行配准图像的方法的算法,用于使用所述算法配准图像的程序以及处理生物医学图像以减少由对象移动导致的成像伪影的方法 |
US20080287787A1 (en) * | 2006-02-02 | 2008-11-20 | Frank Sauer | Line-based calibration of ultrasound transducer integrated with a pose sensor |
JP2010282437A (ja) * | 2009-06-04 | 2010-12-16 | Seiko Epson Corp | 物体認識装置、物体認識方法及びプログラム |
US20110107270A1 (en) * | 2009-10-30 | 2011-05-05 | Bai Wang | Treatment planning in a virtual environment |
CN102499701A (zh) * | 2011-09-29 | 2012-06-20 | 华中科技大学 | X射线和荧光双模式活体成像系统的几何校准方法 |
CN102736112A (zh) * | 2012-07-23 | 2012-10-17 | 复旦大学 | 可变旋转平面的电磁跟踪搜索方法及其系统 |
CN102970926A (zh) * | 2010-06-28 | 2013-03-13 | 皇家飞利浦电子股份有限公司 | Em校准的实时质量控制 |
CN103295455A (zh) * | 2013-06-19 | 2013-09-11 | 北京理工大学 | 基于ct影像模拟与定位的超声培训系统 |
CN103402453A (zh) * | 2011-03-03 | 2013-11-20 | 皇家飞利浦有限公司 | 用于导航系统的自动初始化和配准的系统和方法 |
CN104299260A (zh) * | 2014-09-10 | 2015-01-21 | 西南交通大学 | 一种基于sift和lbp的点云配准的接触网三维重建方法 |
CN104392426A (zh) * | 2014-10-23 | 2015-03-04 | 华中科技大学 | 一种自适应的无标志点三维点云自动拼接方法 |
US20150182191A1 (en) * | 2014-01-02 | 2015-07-02 | Metritrack, Inc. | System and method for tracking completeness of co-registered medical image data |
US20150217958A1 (en) * | 2014-01-31 | 2015-08-06 | Xerox Corporation | Systems and methods for implementing unique offsetting stacker registration using omni-directional wheels for set compiling in image forming devices |
CN105025803A (zh) * | 2013-02-28 | 2015-11-04 | 皇家飞利浦有限公司 | 从多个三维视图对大对象的分割 |
CN105069743A (zh) * | 2015-07-28 | 2015-11-18 | 中国科学院长春光学精密机械与物理研究所 | 探测器拼接实时图像配准的方法 |
CN105078573A (zh) * | 2014-05-11 | 2015-11-25 | 复旦大学 | 基于手持式扫描仪的神经导航空间配准方法 |
WO2016153427A1 (en) * | 2015-03-26 | 2016-09-29 | Nanyang Technological University | Photo-acoustic imaging apparatus and methods of operation |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5999840A (en) * | 1994-09-01 | 1999-12-07 | Massachusetts Institute Of Technology | System and method of registration of three-dimensional data sets |
JP3135068B2 (ja) * | 1994-09-15 | 2001-02-13 | ビジュアリゼイション テクノロジー インコーポレイテッド | 患者の頭部に固定された基準ユニットを用いた医学的応用に使用される位置追跡及び画像生成システム |
EP0951874A3 (en) | 1994-09-15 | 2000-06-14 | Visualization Technology, Inc. | Position tracking and imaging system for use in medical applications using a reference unit secured to a patients head |
US7366562B2 (en) | 2003-10-17 | 2008-04-29 | Medtronic Navigation, Inc. | Method and apparatus for surgical navigation |
US7505809B2 (en) | 2003-01-13 | 2009-03-17 | Mediguide Ltd. | Method and system for registering a first image with a second image relative to the body of a patient |
US7103399B2 (en) | 2003-09-08 | 2006-09-05 | Vanderbilt University | Apparatus and methods of cortical surface registration and deformation tracking for patient-to-image alignment in relation to image-guided surgery |
US7697748B2 (en) | 2004-07-06 | 2010-04-13 | Dimsdale Engineering, Llc | Method and apparatus for high resolution 3D imaging as a function of camera position, camera trajectory and range |
WO2006036842A2 (en) * | 2004-09-24 | 2006-04-06 | The University Of North Carolina At Chapel Hill | Methods, systems, and computer program products for hierarchical registration between a blood vessel and tissue surface model for a subject and blood vessel and tissue surface image for the subject |
JP4836122B2 (ja) | 2006-02-09 | 2011-12-14 | 国立大学法人浜松医科大学 | 手術支援装置、方法及びプログラム |
BRPI0917864A2 (pt) | 2008-08-15 | 2015-11-24 | Univ Brown | aparelho e método para estimativa da forma corporal |
GB0913930D0 (en) | 2009-08-07 | 2009-09-16 | Ucl Business Plc | Apparatus and method for registering two medical images |
FR2963693B1 (fr) | 2010-08-04 | 2013-05-03 | Medtech | Procede d'acquisition automatise et assiste de surfaces anatomiques |
EP2919194B1 (en) | 2011-03-09 | 2020-08-19 | Osaka University | Image data processing device and transcranial magnetic stimulation apparatus |
EP2720171B1 (en) | 2012-10-12 | 2015-04-08 | MVTec Software GmbH | Recognition and pose determination of 3D objects in multimodal scenes |
KR102114417B1 (ko) | 2013-06-11 | 2020-05-22 | 삼성메디슨 주식회사 | 영상 정합 방법 및 장치 |
JP5808446B2 (ja) | 2014-02-24 | 2015-11-10 | キヤノン株式会社 | 情報処理装置および情報処理方法 |
US10152786B2 (en) | 2016-10-11 | 2018-12-11 | Biosense Webster (Israel) Ltd. | Registration of a magnetic tracking system with an imaging device |
-
2016
- 2016-10-11 US US15/290,968 patent/US10152786B2/en active Active
-
2017
- 2017-09-25 IL IL254681A patent/IL254681B/en active IP Right Grant
- 2017-09-26 AU AU2017235896A patent/AU2017235896A1/en not_active Abandoned
- 2017-10-02 CA CA2981027A patent/CA2981027A1/en not_active Abandoned
- 2017-10-10 JP JP2017196651A patent/JP7046553B2/ja active Active
- 2017-10-10 EP EP17195644.4A patent/EP3309749B1/en active Active
- 2017-10-11 CN CN201710941123.9A patent/CN107918925B/zh active Active
-
2018
- 2018-11-28 US US16/202,314 patent/US10535136B2/en active Active
Patent Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050054900A1 (en) * | 2003-07-21 | 2005-03-10 | Vanderbilt University | Ophthalmic orbital surgery apparatus and method and image-guided navigation system |
CN101133431A (zh) * | 2005-02-03 | 2008-02-27 | 布拉科成像S.P.A.公司 | 能够减少物体运动造成的成像伪影的生物医学图像配准方法和计算机程序产品 |
CN101297321A (zh) * | 2005-10-25 | 2008-10-29 | 布拉科成像S.P.A.公司 | 配准图像的方法,用于执行配准图像的方法的算法,用于使用所述算法配准图像的程序以及处理生物医学图像以减少由对象移动导致的成像伪影的方法 |
US20080287787A1 (en) * | 2006-02-02 | 2008-11-20 | Frank Sauer | Line-based calibration of ultrasound transducer integrated with a pose sensor |
CN105574878A (zh) * | 2006-03-21 | 2016-05-11 | 韦伯斯特生物官能公司 | 使用局部加权拟合的图像配准 |
CN101051387A (zh) * | 2006-03-21 | 2007-10-10 | 韦伯斯特生物官能公司 | 使用局部加权拟合的图像配准 |
CN101229080A (zh) * | 2006-10-09 | 2008-07-30 | 韦伯斯特生物官能公司 | 采用器官外部解剖特征的器官图像的配准 |
CN101019765A (zh) * | 2007-03-29 | 2007-08-22 | 新奥博为技术有限公司 | 一种磁共振图像引导下的手术系统及手术导航方法 |
JP2010282437A (ja) * | 2009-06-04 | 2010-12-16 | Seiko Epson Corp | 物体認識装置、物体認識方法及びプログラム |
US20110107270A1 (en) * | 2009-10-30 | 2011-05-05 | Bai Wang | Treatment planning in a virtual environment |
CN102970926A (zh) * | 2010-06-28 | 2013-03-13 | 皇家飞利浦电子股份有限公司 | Em校准的实时质量控制 |
CN103402453A (zh) * | 2011-03-03 | 2013-11-20 | 皇家飞利浦有限公司 | 用于导航系统的自动初始化和配准的系统和方法 |
US20140193053A1 (en) * | 2011-03-03 | 2014-07-10 | Koninklijke Philips N.V. | System and method for automated initialization and registration of navigation system |
CN102499701A (zh) * | 2011-09-29 | 2012-06-20 | 华中科技大学 | X射线和荧光双模式活体成像系统的几何校准方法 |
CN102736112A (zh) * | 2012-07-23 | 2012-10-17 | 复旦大学 | 可变旋转平面的电磁跟踪搜索方法及其系统 |
CN105025803A (zh) * | 2013-02-28 | 2015-11-04 | 皇家飞利浦有限公司 | 从多个三维视图对大对象的分割 |
CN103295455A (zh) * | 2013-06-19 | 2013-09-11 | 北京理工大学 | 基于ct影像模拟与定位的超声培训系统 |
US20150182191A1 (en) * | 2014-01-02 | 2015-07-02 | Metritrack, Inc. | System and method for tracking completeness of co-registered medical image data |
US20150217958A1 (en) * | 2014-01-31 | 2015-08-06 | Xerox Corporation | Systems and methods for implementing unique offsetting stacker registration using omni-directional wheels for set compiling in image forming devices |
CN105078573A (zh) * | 2014-05-11 | 2015-11-25 | 复旦大学 | 基于手持式扫描仪的神经导航空间配准方法 |
CN104299260A (zh) * | 2014-09-10 | 2015-01-21 | 西南交通大学 | 一种基于sift和lbp的点云配准的接触网三维重建方法 |
CN104392426A (zh) * | 2014-10-23 | 2015-03-04 | 华中科技大学 | 一种自适应的无标志点三维点云自动拼接方法 |
WO2016153427A1 (en) * | 2015-03-26 | 2016-09-29 | Nanyang Technological University | Photo-acoustic imaging apparatus and methods of operation |
CN105069743A (zh) * | 2015-07-28 | 2015-11-18 | 中国科学院长春光学精密机械与物理研究所 | 探测器拼接实时图像配准的方法 |
Non-Patent Citations (2)
Title |
---|
李莉等: "基于角度和径向差分函数的抗平移、旋转以及尺度的图像配准技术", 《中国光学期刊网》 * |
鲁通等: "超声引导下微波消融治疗机器人系统的初步研究", 《中国医疗器械信息》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111281539A (zh) * | 2018-12-06 | 2020-06-16 | 韦伯斯特生物官能(以色列)有限公司 | 用于执行面部配准的方法和设备 |
CN113660912A (zh) * | 2019-04-10 | 2021-11-16 | 博医来股份公司 | 受试者的相关表面点的采样方法 |
CN111759463A (zh) * | 2020-07-31 | 2020-10-13 | 南京普爱医疗设备股份有限公司 | 一种提高手术机械臂定位精度的方法 |
Also Published As
Publication number | Publication date |
---|---|
EP3309749B1 (en) | 2021-06-30 |
IL254681A0 (en) | 2017-11-30 |
US10152786B2 (en) | 2018-12-11 |
JP7046553B2 (ja) | 2022-04-04 |
CA2981027A1 (en) | 2018-04-11 |
IL254681B (en) | 2020-01-30 |
JP2018061837A (ja) | 2018-04-19 |
CN107918925B (zh) | 2023-05-02 |
US20180101950A1 (en) | 2018-04-12 |
US20190164282A1 (en) | 2019-05-30 |
US10535136B2 (en) | 2020-01-14 |
EP3309749A1 (en) | 2018-04-18 |
AU2017235896A1 (en) | 2018-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107918925A (zh) | 磁跟踪系统与成像装置的配准 | |
US10296707B2 (en) | System and method for patient-specific image-based guidance of cardiac arrhythmia therapies | |
CN111388092B (zh) | 一种定位跟踪件、配准方法、存储介质及电子设备 | |
JP5334692B2 (ja) | 制約されたマージナル空間学習を使用して3d解剖学的オブジェクトを検出するための方法およびシステム | |
AU2015312327A1 (en) | Systems and methods for segmenting medical images based on anatomical landmark-based features | |
CN111161333B (zh) | 一种肝脏呼吸运动模型的预测方法、装置及存储介质 | |
CN113171062B (zh) | 识别寸、关、尺的方法、终端、系统、介质及计算机设备 | |
US20060122480A1 (en) | Segmenting occluded anatomical structures in medical images | |
US8437518B2 (en) | Registration of electroanatomical mapping points to corresponding image data | |
Xie et al. | Feature‐based rectal contour propagation from planning CT to cone beam CT | |
CN113112486B (zh) | 一种肿瘤的运动估计方法、装置、终端设备和存储介质 | |
CN105793893A (zh) | 解剖结构的基于模型的分割 | |
EP4156096A1 (en) | Method, device and system for automated processing of medical images to output alerts for detected dissimilarities | |
US20240074689A1 (en) | Volumetric lat map | |
CN116824173A (zh) | 医学图像的处理方法、设备及存储介质 | |
KR102084256B1 (ko) | 복수의 후보 대응점을 이용하는 영상 정합 장치 및 방법 | |
CN116583223A (zh) | 解剖结构的3d表示的引导式采集 | |
WO2021137115A1 (en) | Method and apparatus for registering live medical image with anatomical model | |
EP4233761A1 (en) | Method of automatic registration of a 3d model in a medical navigation system display device | |
Baluwala et al. | Automatic landmark detection using statistical shape modelling and template matching | |
WO2024089423A1 (en) | System and method for three-dimensional imaging | |
EP4221597A1 (en) | Methods and systems for tool tracking | |
CN114693750A (zh) | 一种图像配准方法及装置 | |
CN115423789A (zh) | 一种提取cbct图像中定位用金属球位置的方法及装置 | |
Leung et al. | Experimental validation of a 4D elastic registration algorithm |
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 |