多标记物夹具的识别和对准的制作方法

本发明涉及通过拓扑标测的图像变换。更特别地，本发明涉及将2维医学图像与3维标测图组合。
背景技术：
：用于心内膜表面的3维几何标测和重构的方法是本领域中已知的。例如，美国专利号5,738,096(其公开内容以引用方式并入本文)描述了基于使探头与心脏壁上的多个位置接触并且确定在位置中的每个位置处的探头的位置坐标，用于标测心内膜的方法。位置坐标被组合以形成心脏的至少一部分的标测图。基于检测身体内部探头的位置的商业电生理和物理标测系统目前是可商购获得的。在它们中，可从biosensewebster,inc.(3333diamondcanyonroad,diamondbar,ca91765)商购获得的3系统是用于局部电活动与导管位置的自动关联和标测的系统。在当前的心导管插入术系统中，操作医师常常必须同时观察两个不同的屏幕上的两个不同的图像：胸腔的2维荧光镜图像和心脏的3维标测图。例如，可以使用心脏中的导管末端的磁跟踪来生成此类3维标测图。荧光镜图像和3维标测图两者都可以示出导管，但是是从不同的角度和透视图示出。因为荧光镜的视图和3维视图之间缺乏自动对准和配位，所以需要医师在显示器之间来回切换他的或她的注意力，并且在心理上记住它们包含的不同的信息。用于对准解剖图像和电解剖标测图与通过不同模态获取的3维图像的一些现有的方法大体依赖于位置数据。标测导管被放置在感兴趣的器官诸如心脏中的多个已知位置处，并且记录位置坐标。在3维图像中标记或以其它方式记录了这些相同的位置。除被采取作为标测过程本身的一部分的动作之外，该技术通常需要系统的操作者花时间找到和标记用于对准的目的的期望的位置。在专利文献中用于自动对准荧光镜图像与3维标测图的各种方法是已知的。例如，在授予ben-haim等人的共同转让的美国专利号6,314,310中，描述了此类方法，该专利的公开内容以引用方式并入本文。最近，在共同转让的美国专利申请公布号2014/0114173中公开了不同的方法，该专利申请以引用方式并入本文。该文档论述了将胸腔的2维荧光镜图像放置为与心脏的3维标测图功能的电解剖标测图对准。坐标系对准模块包括射线不可透元件，以固定的预先确定的图案布置该射线不可透元件，并且响应于射线不可透元件生成荧光镜图像，该射线不可透元件被配置成将模块的位置限定在荧光镜参考坐标系中。模块还包括一个或多个连接部，该一个或多个连接部被配置成将模块固定地连接到在相对于焊盘预先确定的位置和取向处的磁场传输焊盘，以便表征由磁场传输焊盘限定的磁参考坐标系中的对准模块的位置。在pilu等人的美国专利号8,300,941中提出一种图像变换的方法，该专利涉及通过识别用于表面编码的图案的最佳网格假设来校正透视失真。方法包括：从编码的表面图案提取一组直线假设；按取向将直线假设聚类；对于每个聚类，提取一组线束假设；从成对的线束假设生成一组规则网格假设；以及确定最佳规则网格假设。技术实现要素：本发明的实施方案提供了用于夹具标记物识别的通用方法，该方法是灵活的以支持夹具取向和位置，其中在图像中没有必要整个夹具是可视的。采用恒定尺寸，而不是使用不同尺寸的标记物。根据本发明的实施方案提供了用于对准坐标系的方法，通过定位具有对于成像区域中的成像模态不透明的标记物的校准夹具来执行该方法。标记物包括具有具有相应的交比的至少两个共线标记物四元组。还通过利用成像模态产生标记物的图像、检测图像上的标记物、在图像上识别所检测的标记物的共线图像四元组、计算共线图像四元组的相应的图像交比来执行该方法。还通过以下步骤执行该方法：使根据预定义的标准具有匹配的标记物交比和图像交比的至少一个共线标记物四元组与至少一个共线图像四元组相关联来限定候选线；使用成对的候选线，利用图像执行校准夹具的相应对准；确定对准的残余；以及选择对准中的一个，其中该对准的该对候选线与合格的残余相关联，合格的残余小于预先确定的值。根据该方法的另一方面，候选线的至少一部分具有多个共线图像四元组。根据该方法的另外的方面，沿着单个轴线设定标记物。该方法的又一方面包括使用所选择的对准，使标记物重新投影到图像上，以及此后，通过将所重新投影的标记物中的其它标记物定位在图像上来建立新的对准。根据该方法的另外的方面，定位所重新投影的标记物中的其它标记物包括确定所重新投影的标记物与所检测的标记物之间的接近度。根据该方法的一个方面，定位所重新投影的标记物中的其它标记物包括根据所重新投影的标记物和所检测的标记物的均方根误差来记下所重新投影的标记物。根据该方法的另一方面，建立新的对准包括使在六自由度上重新投影标记物的残余误差最小化。根据该方法的附加的方面，标记物的尺寸是相同的。根据本发明的实施方案还提供了用于对准坐标系的方法，通过在由荧光镜成像器照射的区域中定位具有射线不可透标记物的校准夹具来执行该方法。标记物包括至少两个共线标记物四元组，该至少两个共线标记物四元组具有相应的标记物交比。还通过以下步骤执行该方法：利用荧光镜产生标记物的荧光镜图像；检测图像上的标记物；在荧光镜图像上识别所检测的标记物的共线图像四元组；计算共线图像四元组的相应的图像交比；通过使根据预定义的标准具有匹配的标记物交比和图像交比的至少一个共线标记物四元组与至少一个共线图像四元组相关联来限定候选线；使用成对的候选线，利用荧光镜图像执行校准夹具的相应对准；确定对准的残余；以及选择对准中的一个，其中该对准的该对候选线与合格的残余相关联，合格的残余小于预先确定的值。根据本发明的实施方案，还提供了用于对准坐标系的设备，包括具有射线不可透标记物的校准夹具，其中标记物包括具有相应的标记物交比的至少两个共线标记物四元组。设备还包括处理器；处理器可访问的存储器，该存储器在其中存储程序和数据对象，其中程序的执行致使处理器执行方法，该方法包括检测荧光镜图像上的标记物；在图像上识别所检测的标记物的共线图像四元组；计算共线图像四元组的相应的图像交比；通过使根据预定义的标准具有匹配的标记物交比和图像交比的至少一个共线标记物四元组与至少一个共线图像四元组相关联来限定候选线；使用成对的候选线，利用图像执行校准夹具的相应对准；确定对准的残余；以及选择对准中的一个，其中该对准的该对候选线与合格的残余相关联，合格的残余小于预先确定的值。附图说明为了更好地理解本发明，以举例的方式引用本发明的详细描述，本发明的详细描述应结合以下附图来阅读，附图中相似的元件被给定相似的参考编号，并且在附图中：图1为根据本发明的实施方案示出在校准阶段中荧光镜图像和磁性标测集成系统的示意图；图2为根据本发明的实施方案的在校准阶段中使用的螺旋形校准夹具的示意性透视图；图3为根据本发明的实施方案所应用的投影几何结构的方面的图示；图4是根据本发明的实施方案的与2维图像夹具对准的方法的流程图；图5是根据本发明的实施方案的到2维图像上的夹具的示例性投影；图6是根据本发明的实施方案的到2维图像上的夹具的示例性投影；图7是根据本发明的另选实施方案的到2维图像上的具有单条标记物线的夹具的投影；以及图8是根据本发明的另选实施方案的与2维图像夹具对准的方法的流程图。具体实施方式为了提供对本发明的各种原理的全面理解，在以下描述中阐述了许多具体细节。然而，对于本领域的技术人员将显而易见的是，并非所有这些细节都是实践本发明所必需的。在这种情况下，未详细示出熟知的电路、控制逻辑、以及用于常规算法和过程的计算机程序指令的细节，以免不必要地使一般概念模糊不清。以引用方式并入本专利申请的文献将视为本专利申请的整体部分，不同的是如果在这些并入的文献中定义的任何术语与在本说明书中明确或隐含地给出的定义在某种程度上相冲突，则应只考虑本说明书中的定义。在以下的描述中，描述了荧光镜图像。这仅以举例的方式而不是限制。本文中所描述的过程同样适用于由很多模态产生的2维图像，例如，磁共振图像、计算机断层扫描图像、光学图像和标准x射线图像。现在转向附图，初始参考图1，图1为根据本发明的实施方案示出在校准阶段中的荧光镜图像和磁性标测集成系统20的示意图。系统20将由磁跟踪系统22获取的身体器官的3维标测图与由荧光镜24获取的患者的2维荧光镜图像组合，因此，形成组合显示26，组合显示26在屏幕28上呈现给系统20的操作者。在用于图1中所示的系统20的校准阶段中，未呈现患者。在图5中所示的系统20的后续可操作阶段中，假设患者躺在系统20的手术台30上，并且磁跟踪系统22和荧光镜24获取患者的3维标测图和2维图像，如下面更详细描述的。通常，由荧光镜获取的2维图像是患者的胸部的2维图像，并且由磁跟踪系统标测的身体器官包括患者的心脏。在实践中，虽然系统22和荧光镜24通常可以被配置为具有单独控制单元的单独物理单元，但在本说明书中，为了简单，假设由单个控制单元32操作系统20。控制单元32包括操作荧光镜的荧光镜控制器34和操作磁跟踪系统的磁系统控制器36，并且该单元在系统处理器38特别是在屏幕28上生成组合显示26的处理器的整体控制下。处理器38通常包括利用合适的软件编程以用于执行下文所述功能的通用或嵌入式计算机处理器。因此，尽管处理器38可以有时被描述为包括多个单独的功能块，但是这些功能块并非必须是单独的物理实体，而是表示处理器可访问的存储器中所存储的不同的计算任务或数据对象。这些任务可在运行于单个处理器上或运行于多个处理器上的软件中执行。软件可在与计算机系统一起使用的多种已知非暂时介质诸如磁盘或硬盘驱动器或cd-rom中的任一种上实施。代码可以分布在此类介质上，或者可以通过网络从另一个计算机系统(未示出)的存储器或存储装置分布到处理器38。另选地或除此之外，处理器38可包括数字信号处理器或硬连线逻辑。处理器38通常包括用于数据库29的存储装置。荧光镜24包括x射线源40和检测器42，辐射器和检测器安装在本文中假设为包括圆的圆弧的c形臂44的相对的两端部上。c形臂44通常保持在l形臂46的下端45，l形臂在其上端附接到手术室的天花板，或者在其下端附接到手术室地板。c形臂44可以围绕水平枢轴pa旋转，在该附图中，水平枢轴pa处于纸的平面内。c形臂44还可以围绕c形臂轴ca旋转，c形臂轴ca垂直于纸的平面，并且c形臂轴ca穿过c形臂圆的中心。名义上，两个轴线相交于c形臂旋转中心(也称为等中心)处，并且彼此正交。系统20的操作者能够使用荧光镜控制器34调节围绕轴线pa和轴线ca的荧光镜24的旋转，以及荧光镜的其它几何参数。(如下文所更详细描述，校准阶段提供对荧光镜24的标称特性的调节，标称特性诸如上面所举例说明的特性。)为了对准坐标系，或荧光镜24的参照系，在系统20的校准阶段中，系统20使用一个或多个对准元件。给定的对准元件具有特性：可以在该两个坐标系中同时确定给定的对准元件的位置和取向。本发明的实施方案使用螺旋形校准夹具60作为对准元件。参考图2更详细地描述了螺旋形校准夹具60(本文中还称为夹具60)。两个坐标系(一个基于定位垫50且一个基于对准模块62)呈现在系统20中。现在参考图2，图2为根据本发明的实施方案的夹具60的示意性透视图。如将从以下描述显而易见的，在荧光镜系统中夹具60在可以确定其中夹具60的位置的范围内是“可视的”。夹具60被制成塑料圆筒70，塑料圆筒70可以使用附接的支撑物72被放置在手术台30上，使得圆筒的轴线平行于手术台。圆筒70具有嵌入在圆筒中的相同的金属球体74。金属球体74在检测器42处提供良好的对比度荧光镜图像。只要球体74尺寸是相同的，球体74就可以具有任何适宜的尺寸。相同的球体的使用允许夹具定位的灵活性，因为完全避免了由标记物的尺寸考虑所需要的约束。球体74以螺旋形图案76布置在圆筒内，螺旋形图案被配置成具有变化的周期长度。此外，在每个周期内，存在不同数量的球体。由连接螺旋形的球体74的虚线示出螺旋形图案。而且，有至少两组共线标记物，每组包括四个标记物，被称为四元组。夹具60中的标记物间的距离是已知的，并且可以被编码以使得能够实现唯一识别。在本发明的所公开的实施方案中，圆筒70具有约120mm的直径，并且螺旋形图案76的总体长度约200mm。在所公开的实施方案中，球体74被布置成限定六个周期，并且下面的表1给出在平行于圆筒70的轴线的方向上的每个周期的示例性长度，以及每个周期内球体74的示例性数量。通常，在每个周期内，球体74均匀分布。表1螺旋形周期长度[mm]周期中球体72的数量20725113015351940234527另外的球体74可以被添加到螺旋形图案76中的那些周期。例如，已经添加更多的球体74以形成球体的直线78，线78对应于螺旋形周期的终止端。不同长度的螺旋形图案周期连同每个周期中的不同数量的球体74使得荧光镜控制器34能够唯一地识别球体74中的每个。如上所述，将另外的球体添加到螺旋形图案还有利于荧光镜控制器34识别球体74。球体74的识别被用于对准2维图像与荧光镜，如下面所描述的。在一些实施方案中，圆筒70包括后向反射器90，其中可以定位可拆卸球92。通常，在圆筒70中在平行于手术台30的平面中对称地布置后向反射器，并且可以借助激光器使用后向反射器用于跟踪手术台30上夹具60的位置。在另选的实施方案中，大体类似于球体74的其它金属球体可以被添加到夹具60以另外帮助其定位。对准。现在参考图3，图3为示出根据本发明的实施方案所应用的投影几何结构的方面的图示。方法依赖于投影几何结构的两个特性：(1)空间中的线210被投影到图像中的线212(2)如果点a、点b、点c、点d是线212中的共线点的四元组，则交比被定义为(a,b；c,d)＝(ac*bd)/(bc*ad)。在投影下，共线点的四元组的交比是不变量。在图3中，点a、点b、点c、点d分别是线210中点a'、点b'、点c'、点d'的投影。因此，交比(a,b；c,d)和(a',b'；c',d')相等。根据本发明的原理，夹具与2维图像的初始对准需要最少两条线，每条线具有在图像中可视的至少四个共线标记物。现在参考图4，图4为根据本发明的实施方案的与2维图像夹具对准的方法的流程图。为了呈现清楚起见，在本文中的该流程图和其它流程图中的处理步骤以特别的线性序列示出。然而，将很显然，可并行地、异步地或以不同的次序执行这些步骤中的许多。本领域的技术人员还应当理解，另选地，过程可被表示为多个相互联系的状态或事件，例如，在状态图表中。而且，可能不需要所有示出的过程步骤来实施所述方法。在初始步骤214，确定荧光镜几何参数和荧光镜投影函数之间的关系。可以使用校准的模型(如果可用的话)进行该确定。另选地，可以使用标称、非校准的荧光镜模型。模型提供了足够的近似度，以使得在下面所描述的步骤中能够进行另外的精化。现在参见图5，图5是根据本发明的实施方案的到2维图像上的夹具60(图2)的示例性投影216。标记物218是夹具60上球体74的投影。回转到图4，在步骤220，在投影216上检测标记物218(图5)。可以由方法执行该步骤，方法诸如模板匹配、椭圆检测或边缘检测。这些是常规图像处理技术，并且在本文中不再另外论述这些。步骤220结果是一组已检测的标记物。接下来，在步骤222，在投影216上在步骤220中生成的一组标记物中识别共线标记物。可以通过使用已知的方法诸如随机采样一致性(ransac)或霍夫变换的线拟合来执行该步骤。步骤222产生包括是用于共线性的候选者的至少标记物的四元组的一组线假设(即，子组)。每个四元组定义交比。是用于共线性的候选者的所有四元组的列表定义特征。子组包包括其中一个或多个标记物缺失的所有的可能性，并且列表指定缺失的标记物。例如，在5个标记物的线中，有一个5标记物子组和5个4标记物子组。每个子组被存储指定该子组的成员的相应的索引。接下来，在步骤224，通过使在步骤222中制定的线假设的交比与每条线内的夹具上的已知的标记物间距离的交比匹配，评估在步骤222中制定的线假设。步骤224包括若干过程：在步骤226中，评估夹具中四个或更多个共线标记物的全部已知的子组，以便覆盖未检测的标记物或在图像外的标记物的情况。每个四元组定义交比。关于模拟其中一个或多个标记物缺失的所有可能的选项，创建标记物的子组。缺失的标记物可以是在各种位置，并且甚至可以是在图6中所示的视场的外部。对于每个子组，存在对所有可能的情况进行说明的“特征”。可以使用有效地实施的穷举法确定特征，通过在每个子组内所有四元组上的迭代，计算所有子组的此类标记物交比，从而创建子组特定的特征。交比特征被保存在参考数据库29(图1)中。例如，可以保存四个连续的标记物的所有组的交比。对于每个夹具设计，仅需要执行一次步骤226。结果可以应用于具有相同的设计的所有后续制造的夹具。在步骤228中，为与该线假设有关且在步骤222中识别的图像上的每个子组的标记物中的每个四元组计算每个线假设交比。每个四元组定义交比。关于模拟其中共线标记物中的任一个可能是假的内围层(即，其中存在标记物意外地位于图像上的线上但是标记物在夹具机械结构上不是共线的情况)的所有可能的选项，创建标记物的子组。线假设交比是以与在步骤226中获得的参考数据库中的特征相同的格式在图像交比数据库中保存的特征。参考数据库和图像交比数据库两者可以包括在数据库29(图1)中。在匹配步骤230中，进行在步骤228中获得的交比特征与在步骤226中获得的交比特征之间的匹配的搜索。目的是找到在图像中的共线标记物与在夹具中它们的对应的标记物之间的潜在的匹配。每个匹配是用于标记物的子组中的标记物的识别的选项。通过比较相等尺寸的子组执行匹配。通过欧几里得距离，比较在步骤228中的交比与在步骤226中的交比，以找到两组数的对应。另选地，可以使用距离的其它测量。在任何情况下，当发现图像中的子组和夹具中的子组比预定义的阈值距离更加靠近彼此时，宣告匹配。在步骤232中，根据标记物的数量，或另选地根据特征匹配的均方根误差(rmse)，对在步骤230中获得的匹配进行优先排序，创建与两条夹具线匹配的线假设的优先排序的列表。在步骤238中，执行初始夹具对准过程。评估来自步骤232中挑选的线假设的候选线的穷举成对选择。以根据关联性和相似性的次序评估线对，将较大的子组的匹配放在前面。将被评估的第一对是与夹具中的一条线匹配的图像中的最大子组和与另一条线匹配的图像中的最大子组。如果这不产生符合要求的初始对准，则渐进地测试较小的子组。根据最小rmse值，每个测试需要初始对准和批准或拒绝。其它考虑包括图像中标记物的位置。根据c形臂取向，可以期望一条线在第二条线的左边(或右边)。交比匹配质量的质量可以是另一个决定性标准。步骤238将一个或多个图像线候选者给予给夹具线中的每条(仅一条是真的匹配)。然后在以下步骤中考虑和评估与两条夹具线匹配的每对图像线候选者，直到达到有效的对准。对于效率，具有较长的特征(即，较长的子组)的候选者在前面。将被评估的第一对是与夹具中的一条线匹配的图像中的最大子组和与另一条线匹配的图像中的最大子组。除此之外，图像中标记物的位置是另一个考虑。用于评估算法的终止标准是期望的，例如，迭代的阈值数量。然而，在实践中，几乎总是在几个迭代内达到正确的匹配。使用每个此类对执行与图像的夹具对准，并且通过将标记物重新投影到图像上，并且测量在第一投影上所检测的标记物与在重新投影上所检测的标记物之间的距离来确定残余。接下来，在决定步骤234，确定初始夹具对准是否是合格的残余，即，残余小于用户配置的标准。如果在决定步骤234处确定是否定的，则控制返回步骤232，其中在步骤232中评估另一个线对。如果在决定步骤234处确定是肯定的，则控制前进至步骤236。一旦确定有效的初始对准，则在步骤236处，基于初始对准，识别剩余的所检测的标记物。这通过在步骤238中使用线对的初始对准将夹具重新投影到图像上来实现。接近度标准(即，所检测的标记物和重新投影的标记物之间的距离)(例如，rmse)被用于识别剩余的标记物。如果初始对准足够，则将以最小误差定位所有剩余的标记物。在最终步骤240，使用在步骤236中识别的所有标记物来执行最终夹具对准。这通过寻找产生标记物(在步骤236中所识别的)的重新投影的最小残余误差的六自由度上的夹具位置和取向的最小化函数来进行。现在参考图6，图6是根据本发明的实施方案的到2维图像上的夹具60(图2)的完全对准的投影242。现在由标记(例如，相邻的数字)唯一识别标记物218。线244、线246是在步骤238、步骤240(图4)中执行的成功的初始对准和后续完全对准中挑选的所使用的候选线组中的最佳对。另选实施方案。该实施方案类似于先前的实施方案。然而，步骤238(图4)被修改以执行五自由度而不是六自由度上的初始夹具对准。这使用单个挑选的线假设而不是两条线来进行。该选项对围绕单个轴对称的夹具有用，该夹具诸如在对称轴上具有标记物的杆。该实施方案节省时间和计算机资源，其中仅需要确定图像中所挑选的线的轴线。现在参考图7，图7是根据本发明的另选实施方案的到2维图像上的具有单条标记物线的夹具的完全对准的投影248。标记物之间的距离有意图地变化，以便使得能够有各种交比。现在参考图8，图8为根据本发明的另选实施方案的与2维图像夹具对准的方法的流程图。图8中的步骤与在图4的论述中所描述的那些步骤相同，直到步骤250。为简洁起见，此处不再重复细节。在步骤250中，执行彻底的夹具对准，而不是如在第一实施方案中所描述的初步对准。接下来，在决定步骤252处，确定在步骤250中执行的对准是否是可接受的。如果确定是否定的，则控制返回到步骤232。如果在决定步骤252处确定为肯定的，则控制前进到最终步骤254，在最终步骤254处过程终止。如在步骤236和最终步骤240(图4)中所描述的其它标记物的识别和第二对准是不必要的。由于步骤250的对准仅解决五自由度，所以线性夹具围绕其自己的轴线的旋转不会影响投影。实施例。在该实施例中，激光跟踪设备被用于验证用于螺旋形夹具的校准和对准的上面所描述的算法的精度和重复性。比较用于相对于荧光镜检查系统的不同的螺旋形夹具位置的校准和对准结果。验证荧光镜检查系统坐标系和激光跟踪器坐标系之间所计算的变换维持恒定。在不同的螺旋形夹具位置和取向处对螺旋形夹具拍摄各种校准序列和对准图像。还使用激光跟踪器测量每个螺旋形夹具位置。从利用激光跟踪器获得的螺旋形夹具位置测量，和用于螺旋形夹具坐标系中的该测量的后向反射器的已知的机械位置，可能计算从激光跟踪器坐标系到螺旋形夹具坐标系的变换。为每个螺旋形夹具位置计算laser2helix变换。每个校准和对准算法的结果包括从螺旋形坐标系到荧光镜坐标系的所计算的变换。为每个校准结果且为每个对准图像计算helix2fluoro变换。根据上面的变换，从荧光镜坐标系到激光跟踪器坐标系的变换可以计算如下：方程1：fluoro2laser＝(laser2helix)-1*(helix2fluoro)-1因此，为每个校准结果且为每个对准图像计算fluoro2laser变换。由于在整个该测试中荧光镜坐标系和激光跟踪器坐标系是恒定的，所以该变换应是相同的，与螺旋形夹具位置或校准/对准图像无关。从用于每个螺旋形夹具位置的参考helix2fluoro变换，使用上面的方程1计算参考fluoro2laser。上面所计算的变换被用于以下比较，以验证执行校准和对准算法的荧光镜检查校准工具的精度和重复性。校准算法为不同的校准结果计算的所有fluoro2laser变换之间的比较。使用不同的校准结果，在荧光镜坐标系中的相同的点到在各种c形臂位置处拍摄的相同图像中的投影之间的比较。对准算法为相同的螺旋形夹具位置拍摄的不同的对准图像计算的所有fluoro2laser变换之间的比较。越过所有螺旋形夹具位置，为在各种位置处以c形臂拍摄的所有对准图像计算的所有fluoro2laser变换之间的比较。备注：理想的是，上面的比较应该产生接近于零的值；实际值是算法重复性的测量。测试工具和辅助仪器。荧光镜检查系统。具有univu模块的carto系统。激光跟踪器(可从farotechnologies,inc(250technologyparklakemary,fl32746,usa)商购获得)和后向反射器。螺旋形夹具。软件。从microsoftvisual源控制)取得的分析代码。从microsoftvisualstudio源控制)取得的fctfluoro积分算法。测试过程。校准算法测试计划。螺旋形夹具被放置在荧光镜cor中的患者手术台上，并且采集校准图像序列，即，在以下荧光镜c形臂位置处获取图像：lao90；lao60；lao30；ap；rao30；rao60；rao90；cra20；cra40；cau20和cau40。为患者手术台上的若干稍微不同的螺旋形夹具位置重复数据采集，如在下面的结果节中的表2中所详细说明的。根据在测试概念和方法节中描述的方法，使用在工具和仪器节中指定的分析代码执行数据分析。对准算法测试计划。螺旋形夹具被放置在荧光镜cor中的患者手术台上，并且在各种c形臂角度处拍摄对准图像。为患者手术台上的若干不同的患者手术台/螺旋形夹具位置重复数据采集，如在下面的结果节中的表6中所详细说明的。根据在上面的测试概念和方法节中描述的方法，使用在上面的工具和仪器节中指定的分析代码执行数据分析。接受标准。校准算法：为不同的校准结果计算的所有fluoro2laser变换之间的比较将产生不超过1mm的平均值。使用不同的校准结果，在荧光镜坐标系中的相同的点到在各种c形臂位置处拍摄的相同图像中的投影之间的比较将产生不超过1mm的平均值。对准算法。为相同的螺旋形夹具位置拍摄的不同的对准图像计算的所有fluoro2laser变换之间的比较将产生不超过1mm的平均值。越过所有螺旋形夹具位置，为在ap位置处以c形臂拍摄的不同的对准图像计算的所有fluoro2laser变换之间的比较将产生不超过1mm的平均值。越过所有螺旋形夹具位置，为在lao90位置处以c形臂拍摄的不同的对准图像计算的所有fluoro2laser变换之间的比较将产生不超过1mm的平均值。为用于不同的螺旋形夹具位置的对准图像计算的参考fluoro2laser变换之间的比较将产生不超过1mm的平均值。校准算法测试结果。下面的表2详细说明用于所收集的不同校准序列图像的螺旋形夹具位置中的差别。表2.校准数据集合。用于所有校准的患者手术台位置是[-54,-865,930]，不同的是第一校准(t124248)，在第一校准(t124248)中用于所有校准的患者手术台位置是[-59,-864,930]。用于所有校准的荧光镜缩放设置是42cm。下面的表3示出用于不同的校准的所校准的参数值。表3.用于所执行的校准的所校准的参数值。表4示出为不同的校准结果计算的所有fluoro2laser变换之间的比较。值表示所比较的fluoro2laser变换之间的平移中的差值的范数，即，如下执行计算：方程2：fluoro2lasercomparison12＝(fluoro2laser1)-1*fluoro2laser2方程3：δ＝norm(fluoro2lasercomparison(1:3,4))。表4中示出δ的值。表4用于所有校准的fluoro2laser变换的比较。校准#1校准#2校准#3校准#4校准#20.48mm校准#30.58mm0.11mm校准#40.47mm0.08mm0.12mm校准#50.61mm0.14mm0.03mm0.15mm上面比较的平均值是0.28mm，其低于1mm的接受标准阈值，从而测试通过。表5使用不同的校准结果示出在荧光镜坐标系中的相同的点到在各种c形臂位置处拍摄的相同图像中的投影。表5中示出的值表示从像素转换成毫米的投影之间的平均和最大差值。表5荧光镜与图像投影的比较对于所有校准上面比较的平均值是0.18mm，其低于1mm的接受标准阈值。因而，测试通过。对准算法测试结果。表6详细说明用于所收集的不同对准图像的螺旋形夹具位置中的区别。表6对准数据收集计算用于每个对准图像的helix2fluoro对准结果与用于该螺旋形夹具位置的总体参考helix2fluoro之间的差值，并且该差值被显示在附图中，附图为每个螺旋形夹具位置示出平移矩阵t元素[tx,ty,tz]对c形臂角度。表7示出为相同的螺旋形夹具位置拍摄的不同的对准图像计算的所有fluoro2laser变换的比较。使用上面所示的方程2和方程3计算表7中所示的值，并且表7中所示的值表示为相同的螺旋形夹具位置拍摄的对准图像所计算的fluoro2laser变换之间的平均和最大差值，以及cdf95值，即，95％的所估计的较高阈值。对于图像对准之间的每个比较，仅为垂直于两个图像角度的方向计算δ。表7.用于在相同的螺旋形夹具位置处的所有对准图像的fluoro2laser变换的比较位置平均δ[mm]最大δ[mm]cdf9510.151.110.5320.060.210.1830.180.790.7440.160.430.3950.200.970.70上面比较的平均值是0.15mm，其低于1mm的接受标准阈值，从而测试通过。表8示出越过所有螺旋形夹具位置，为在ap位置处和在lao90位置处以c形臂拍摄的不同的对准图像计算的所有fluoro2laser变换之间的比较。使用上面所示的方程2和方程3计算表8中所示的值，并且表8中所示的值表示为相同的c形臂位置拍摄的对准图像计算的fluoro2laser变换之间的平均和最大差值，以及cdf95值，即，95％的所估计的较高阈值。对于图像对准之间的每个比较，仅为垂直于图像c形臂角度的方向计算δ(即，z轴δ不被考虑到ap图像比较中，并且x轴δ不被考虑到lao90图像比较中)。表8.用于在相同的c形臂角度处的对准图像的fluoro2laser变换的比较c形臂角度平均δ[mm]最大δ[mm]cdf95ap0.401.071.06lao900.090.230.18上面比较的平均值是0.4mm和0.09mm，其都低于1mm的接受标准阈值。因而，测试通过。表9示出为不同的螺旋形夹具位置的对准图像计算的参考fluoro2laser变换之间的比较。使用上面所示的方程2和方程3计算表9中所示的值。表9.用于不同的螺旋形夹具位置的参考fluoro2laser变换的比较上面比较的平均值是0.63mm，其低于1mm的接受标准阈值。因而，测试通过。本领域的技术人员应当理解，本发明并不限于上文中特别示出和描述的内容。相反，本发明的范围包含上文所述各种特征的组合与子组合两者，以及不在现有技术范围内的其变型和修改，所属领域的技术人员在阅读上述描述时应当想到这些变型和修改。当前第1页12