CN103860268A - 一种标记点配准方法、装置及外科手术导航系统 - Google Patents

Abstract

本发明适用于医学图像处理及应用技术领域，提供了一种标记点配准的方法、装置及外科手术导航系统，所述方法包括：根据多次获取的各个物理空间标记点的物理坐标值和所述物理坐标值的平均值计算出各个物理空间标记点的空间定位误差；根据所述空间定位误差计算出各个物理空间标记点的权重矩阵；根据各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值，结合所述权重矩阵，在保证标记点配准误差FRE最小的前提下，用迭代法求解平移向量t和旋转矩阵R。本发明，可以充分考虑光学定位系统空间定位误差的各向异性性质，最终提高整个外科手术导航系统的导航精度。

Description

一种标记点配准方法、装置及外科手术导航系统

技术领域

本发明属于医学图像处理及应用技术领域，尤其涉及一种标记点配准方法、装置及外科手术导航系统。

背景技术

外科手术导航是医生进行微创外科手术的必要辅助技术，它可以结合术前和术中图像，为外科医生提供更加丰富的病人解剖信息。手术过程中，医生可以通过该外科手术导航实时监测手术器械与病人手术区域的空间位置。光学定位系统，如，加拿大NDI公司生产的Polaris系统，由于定位设备自身无法克服的定位误差，即光学定位相机与所跟踪对象的连线所在轴的方向上产生的误差是其他两个正交方向的三到五倍。换言之，通过光学定位系统所获取的手术器械的空间位置信息具有各向异性的性质。光学定位系统所获取的空间位置信息参与手术导航的配准过程时，空间位置信息自身的误差会累积到手术导航配准结果中，导致最终的外科手术导航系统误差加大。

配准是外科手术导航系统中至关紧要的一个步骤，它实现病人物理空间和图像空间坐标系的统一。目前所用的外科手术导航系统，术前采集的图像在术中必须通过配准，计算病人物理空间到图像空间的变换矩阵。一般的做法是在病人体表上粘贴有限数目的标记点（或通过识别具有一定特征的解剖标志点），然后利用光学定位系统来获取这些标记点的物理空间坐标，利用最小二乘法，求解图像空间对应部位的标记点与光学定位系统所获取的物理空间坐标之间的变换矩阵，最终实现图像空间和病人物理空间的统一。

现有技术也提到了外科手术导航中的三种误差，分别为标记点定位误差(Fiducial Localization Error,FLE)、标记点配准误差(Fiducial Registration Error,FRE)以及目标点配准误差(Target Registration Error,TRE)。标记点定位误差是指在图像空间或者病人物理空间选取用于配准的标记点时，标记点的真实位置与选取所得的坐标之间的距离；标记点配准误差是指在执行配准过程以后，图像空间与物理空间对齐的基础上，用于配准的两个空间坐标系对应标记点距离的均方根值；目标配准误差是指执行配准过程以后，图像空间与物理空间对齐的基础上，两个空间坐标系对应的非配准标记点距离的均方根值。因此，通过光学定位系统获取到的标记点空间坐标值，是计算手术导航配准变换矩阵的基础，估算出标记点定位误差，并考虑光学定位系统定位误差的各向异性性质，有利于提高整个手术导航系统的导航精度。

发明内容

本发明实施例提供了一种标记点配准的方法、装置及外科手术导航系统，旨在解决现有技术由于光学定位系统自身的定位误差的各向异性性质，导致外科手术导航系统的导航精度不高的问题。

一方面，提供一种标记点配准方法，所述方法包括：

分别多次获取各个物理空间标记点的物理坐标值和各个图像空间标记点的图像坐标值；

分别计算各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值；

根据多次获取的各个物理空间标记点的物理坐标值和所述物理坐标值的平均值计算出各个物理空间标记点的空间定位误差；

根据所述空间定位误差计算出各个物理空间标记点的权重矩阵；

根据各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值，结合所述权重矩阵，在保证标记点配准误差FRE最小的前提下，用迭代法求解平移向量t和旋转矩阵R。

另一方面，提供一种标记点配准装置，所述装置包括：

第一坐标获取单元，用于分别多次获取各个物理空间标记点的物理坐标值和各个图像空间标记点的图像坐标值；

第二坐标计算单元，用于分别计算各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值；

空间定位误差计算单元，用于根据多次获取的各个物理空间标记点的物理坐标值和所述物理坐标值的平均值计算出各个物理空间标记点的空间定位误差；

权重矩阵计算单元，用于根据所述空间定位误差计算出各个物理空间标记点的权重矩阵；

配准单元，用于根据各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值，结合所述权重矩阵，在保证标记点配准误差FRE最小的前提下，用迭代法求解平移向量t和旋转矩阵R。

再一方面，提供一种外科手术导航系统，所述外科手术导航系统包括如上所述的标记点配准装置。

在本发明实施例，在进行外科手术导航配准时，首先计算选取的物理空间上各个标记点的空间定位误差，然后根据该空间定位误差来计算各个标记点的权重矩阵，最后根据该权重矩阵来进行导航配准。由于选取的标记点是不在同一平面上的不共面的标记点，因此可以充分考虑光学定位系统空间定位误差的各向异性性质，最终提高整个外科手术导航系统的导航精度。

附图说明

图1是本发明实施例一提供的标记点配准方法的实现流程图；

图2是本发明实施例二提供的标记点配准装置的结构框图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

在本发明实施例中，在进行外科手术导航配准时，首先计算选取的物理空间上各个标记点的空间定位误差，然后根据该空间定位误差来计算各个标记点的权重矩阵，最后根据该权重矩阵来进行导航配准。

以下结合具体实施例对本发明的实现进行详细描述：

实施例一

图1示出了本发明实施例一提供的标记点配准方法的实现流程，详述如下：

在步骤S101中，分别多次获取各个物理空间标记点的物理坐标值和各个图像空间标记点的图像坐标值。

在本实施例中，在手术区域的表面贴放若干可以在CT或MR等成像方式显影的标记点，如鱼肝油颗粒，然后进行成像，得到一组带标记点的图像。使用外科手术导航系统选取图像空间的标记点，将所述图像空间标记点的图像坐标值记为{p_i}，其中，i=1,2,…,N，N≥3，N表示选取的标记点的个数，在本实施例中，N的值最小为3，表示最少选取3个图像空间标记点，并要求这些图像空间标记点尽量不在同一平面上。

由于标记点的体积大小比一般的CT或MR图像的体素要大很多，因此在图像空间上不考虑定位标记点时所引起的误差，但为了使图像空间标记点的坐标值更准确，选取每个图像空间标记点的次数不应少于5次，且所选取的标记点分布要尽可能地分散，不能共线，然后计算5次选取的标记点的坐标值的平均值，将该平均值作为该标记点图像空间的真实坐标值，用于后续的配准过程的计算。

使用光学定位系统选取物理空间的标记点，将该物理空间标记点的物理坐标值记为{f_i}，其中，i=1,2,…,N，N≥3，其中N表示选取的标记点的个数，在本实施例中，N的值最小为3，表示最少选取3个物理空间标记点，并要求这些物理空间标记点尽量不在同一平面上。

另外，利用光学定位系统选取物理空间标记点的物理坐标值时，对于每一个标记点的物理坐标值分别获取n次(n≥5)，而且要与图像空间的标记点一一对应，每个标记点可以得到n个物理坐标值。对于某个特定的标记点i，按x,y,z三个分量，将n个物理坐标值分成三组，分别记为[x_i1,x_i2,...,x_in]、[y_i1,y_i2,...,y_in]和[z_i1,z_i2,…,z_in]。然后计算n次选取的标记点的坐标值的平均值，将该平均值作为该标记点物理空间的真实坐标值，用于后续的配准过程的计算。

在步骤S102中，分别计算各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值。

在本实施例中，先将多次获取到的标记点i的n个物理坐标值按坐标轴的方向划分为三个分量[x_i1,x_i2,...,x_in]、[y_i1,y_i2,...,y_in]和[z_i1,z_i2,...,z_in]后，可以通过下述公式计算得到标记点i的物理坐标值在x、y、z三个分量上的平均值：

${\stackrel{\‾}{x}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{x}_{\mathrm{ij}}$

${\stackrel{\‾}{y}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{y}_{\mathrm{ij}}$

${\stackrel{\‾}{z}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{z}_{\mathrm{ij}}$

其中，x_ij表示标记点i在x方向上的第j个物理坐标值，y_ij表示标记点i在y方向上的第j个物理坐标值，z_ij表示标记点i在z方向上的第j个物理坐标值，n为获取同一标记点的物理坐标值的次数。

在步骤S103中，根据多次获取的各个物理空间标记点的物理坐标值和所述物理坐标值的平均值计算出各个物理空间标记点的空间定位误差。

在本实施例中，通过下述公式计算得到物理空间标记点i的n个物理坐标值的空间定位误差：

${\mathrm{\Δx}}_{i1}=x_{i1}-{\stackrel{\‾}{x}}_i,{\mathrm{\Δx}}_{i2}=x_{i2}-{\stackrel{\‾}{x}}_i,...,{\mathrm{\Δx}}_{\mathrm{in}}=x_{\mathrm{in}}-{\stackrel{\‾}{x}}_i$

${\mathrm{\Δy}}_{i1}=y_{i1}-{\stackrel{\‾}{y}}_i,{\mathrm{\Δy}}_{i2}=y_{i2}-{\stackrel{\‾}{y}}_i,...,{\mathrm{\Δy}}_{\mathrm{in}}=y_{\mathrm{in}}-{\stackrel{\‾}{y}}_i$

${\mathrm{\Δz}}_{i1}=z_{i1}-{\stackrel{\‾}{z}}_i,{\mathrm{\Δz}}_{i2}=z_{i2}-{\stackrel{\‾}{z}}_i,...,{\mathrm{\Δz}}_{\mathrm{in}}=z_{\mathrm{in}}-{\stackrel{\‾}{z}}_i$

其中，Δx_i1至Δx_in表示物理空间标记点i在x方向上的n个物理坐标值的空间定位误差，Δy_i1至Δy_in表示物理空间标记点i在y方向上的n个物理坐标值的空间定位误差，Δz_i1至Δz_in表示物理空间标记点i在z方向上的n个物理坐标值的空间定位误差。

在步骤S104中，根据所述空间定位误差计算出各个物理空间标记点的权重矩阵。

在本实施例中，通过以下步骤计算各个物理空间标记点的权重矩阵：

步骤1、计算出各个物理空间标记点的空间定位误差在各个方向上的方差分量以及各个物理空间标记点的空间定位误差的协方差矩阵。

在本实施例中，假设每个标记点的空间定位误差是正态分布的，期望值为0，并且方差可以分解为沿x,y,z方向的相互独立的三个分量，记为σ_iα，α=1,2,3，分别代表x,y,z三个方向。其中x方向的方差分量由数据Δx_i1，Δx_i2，……Δx_in得到，y方向的方差分量由数据Δy_i1，Δy_i2，……Δy_in得到，z方向的方差分量由数据Δz_i1，Δz_i2，……Δz_in得到，即：

σ_i1=(Δx_i1 ²+Δx_i2 ²+……+Δx_in ²)/n

σ_i2=(Δy_i1 ²+Δy_i2 ²+……+Δy_in ²)/n

σ_i3=(Δz_i1 ²+Δz_i2 ²+……+Δz_in ²)/n

在本实施例中，物理空间标记点的空间定位误差的协方差矩阵为满足：

${\mathrm{\Σ}}_i^{\mathrm{fle}}=\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\left[\begin{array}{ccc}{\mathrm{\Δx}}_{\mathrm{ij}}{\mathrm{\Δx}}_{\mathrm{ij}}& {\mathrm{\Δx}}_{\mathrm{ij}}{\mathrm{\Δy}}_{\mathrm{ij}}& {\mathrm{\Δx}}_{\mathrm{ij}}{\mathrm{\Δz}}_{\mathrm{ij}}\\ {\mathrm{\Δy}}_{\mathrm{ij}}{\mathrm{\Δx}}_{\mathrm{ij}}& {\mathrm{\Δy}}_{\mathrm{ij}}{\mathrm{\Δy}}_{\mathrm{ij}}& {\mathrm{\Δy}}_{\mathrm{ij}}{\mathrm{\Δz}}_{\mathrm{ij}}\\ {\mathrm{\Δz}}_{\mathrm{ij}}{\mathrm{\Δx}}_{\mathrm{ij}}& {\mathrm{\Δz}}_{\mathrm{ij}}{\mathrm{\Δy}}_{\mathrm{ij}}& {\mathrm{\Δz}}_{\mathrm{ij}}{\mathrm{\Δz}}_{\mathrm{ij}}\end{array}\right]$

步骤2、对所述协方差矩阵进行奇异值分解，得到各个物理空间标记点的空间定位误差的主轴。

在本实施中，对

进行分解，即使

满足

其中，分解得到的V_i是标记点i的FLE的主轴。

步骤3、根据所述主轴以及各个物理空间标记点的空间定位误差在各个方向上的方差分量，计算出各个物理空间标记点的权重矩阵。

在本实施例中，通过如下公式计算得到各个物理空间标记点的权重矩阵W_i，W_i满足：

$W_i={(V_i*\mathrm{diag}({{\mathrm{\σ}}_{i1}}^{-2},{{\mathrm{\σ}}_{i2}}^{-2},{{\mathrm{\σ}}_{i3}}^{-2}){*V}_i^T)}^{-1/2}$

在步骤S105中，根据各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值，结合所述权重矩阵，在保证标记点配准误差FRE最小的前提下，用迭代法求解三维平移向量t和旋转矩阵R。

在本实施例中，外科手术导航配准的目的就是经过三维平移向量t和旋转矩阵R的变换后，使物理空间坐标系和图像空间坐标系上对应点坐标值的距离d_i的均方根值最小，并考虑光学定位系统获取标记点空间位置时所产生误差的各向异性的性质。

即使得下述公式的值要求最小：

${\mathrm{FRE}}^2=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\left|W_i(R{\stackrel{\‾}{f}}_i+t-{\stackrel{\‾}{p}}_i)\right|}^2$

其中，

表示标记点i的物理坐标值的平均值，

表示标记点i的图像坐标值的平均值。

本实施例，在进行外科手术导航配准时，首先计算选取的物理空间上各个标记点的空间定位误差，然后根据该空间定位误差来计算各个标记点的权重矩阵，最后根据该权重矩阵来进行导航配准。由于选取的标记点是不在同一平面上的不共面的标记点，因此可以充分考虑光学定位系统空间定位误差的各向异性性质，最终提高整个外科手术导航系统的导航精度。另外，本实施例是在原来手术导航选取标记点空间坐标的过程中，增加对每个标记点选取的次数，空间定位误差的估算过程由算法自动完成，算法的实时性和实用性强，操作简单，应用范围广，适用于神经外科导航、骨科手术导航及腹部等软组织的手术导航中。

实施例二

图2示出了本发明实施例二提供的标记点配准装置的具体结构框图，为了便于说明，仅示出了与本发明实施例相关的部分。该标记点配准装置可以是外科手术导航系统中的一个软件单元、硬件单元或者软硬件结合的单元。在本实施例中，该标记点配准装置2包括：第一坐标获取单元21、第二坐标计算单元22、空间定位误差计算单元23、权重矩阵计算单元24和配准单元25。

其中，第一坐标获取单元21，用于分别多次获取各个物理空间标记点的物理坐标值和各个图像空间标记点的图像坐标值；

第二坐标计算单元22，用于分别计算各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值；

空间定位误差计算单元23，用于根据多次获取的各个物理空间标记点的物理坐标值和所述物理坐标值的平均值计算出各个物理空间标记点的空间定位误差；

权重矩阵计算单元24，用于根据所述空间定位误差计算出各个物理空间标记点的权重矩阵；

配准单元25，用于根据各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值，结合所述权重矩阵，在保证标记点配准误差FRE最小的前提下，用迭代法求解平移向量t和旋转矩阵R。

进一步地，所述权重矩阵计算单元24包括：方差计算模块、主轴计算模块和权重矩阵计算模块。

其中，方差计算模块，用于计算出各个物理空间标记点的空间定位误差在各个方向上的方差分量以及各个物理空间标记点的空间定位误差的协方差矩阵；

主轴计算模块，用于对所述协方差矩阵进行奇异值分解，得到各个物理空间标记点的空间定位误差的主轴；

权重矩阵计算模块，用于根据所述主轴以及各个物理空间标记点的空间定位误差在各个方向上的方差分量，计算出各个物理空间标记点的权重矩阵。

进一步地，所述第二坐标计算单元22包括：坐标值划分模块和平均值计算模块。

其中，坐标值划分模块，用于将多次获取的各个物理空间标记点的物理坐标值按坐标轴的方向划分为三个分量，分别记为x_ij、y_ij、z_ij；

平均值计算模块，用于计算每个物理坐标值在所述三个分量上的平均值，分别记为

具体计算公式如下：

${\stackrel{\‾}{x}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{x}_{\mathrm{ij}}$

${\stackrel{\‾}{y}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{y}_{\mathrm{ij}}$

${\stackrel{\‾}{z}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{z}_{\mathrm{ij}}$

其中，x_ij表示标记点i在x方向上的第j个物理坐标值，y_ij表示标记点i在y方向上的第j个物理坐标值，z_ij表示标记点i在z方向上的第j个物理坐标值，n为获取同一标记点的物理坐标值的次数。

本发明实施例提供的标记点配准装置可以应用在前述对应的方法实施例一中，详情参见上述实施例一的描述，在此不再赘述。

值得注意的是，上述系统实施例中，所包括的各个单元只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。

另外，本领域普通技术人员可以理解实现上述各实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，相应的程序可以存储于一计算机可读取存储介质中，所述的存储介质，如ROM/RAM、磁盘或光盘等。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种标记点配准方法，其特征在于，所述方法包括：

分别多次获取各个物理空间标记点的物理坐标值和各个图像空间标记点的图像坐标值；

分别计算各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值；

根据多次获取的各个物理空间标记点的物理坐标值和所述物理坐标值的平均值计算出各个物理空间标记点的空间定位误差；

根据所述空间定位误差计算出各个物理空间标记点的权重矩阵；

根据各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值，结合所述权重矩阵，在保证标记点配准误差FRE最小的前提下，用迭代法求解平移向量t和旋转矩阵R。

2.如权利要求1所述的方法，其特征在于，所述根据所述空间定位误差计算出各个物理空间标记点的权重矩阵具体包括：

计算出各个物理空间标记点的空间定位误差在各个方向上的方差分量以及各个物理空间标记点的空间定位误差的协方差矩阵；

对所述协方差矩阵进行奇异值分解，得到各个物理空间标记点的空间定位误差的主轴；

根据所述主轴以及各个物理空间标记点的空间定位误差在各个方向上的方差分量，计算出各个物理空间标记点的权重矩阵。

3.如权利要求1所述的方法，其特征在于，所述计算各个物理空间标记点的物理坐标值的平均值具体包括：

将多次获取的各个物理空间标记点的物理坐标值按坐标轴的方向划分为三个分量，分别记为x_ij、y_ij、z_ij；

计算每个物理坐标值在所述三个分量上的平均值，分别记为

具体计算公式如下：

${\stackrel{\‾}{x}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{x}_{\mathrm{ij}}$

${\stackrel{\‾}{y}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{y}_{\mathrm{ij}}$

${\stackrel{\‾}{z}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{z}_{\mathrm{ij}}$

其中，x_ij表示标记点i在x方向上的第j个物理坐标值，y_ij表示标记点i在y方向上的第j个物理坐标值，z_ij表示标记点i在z方向上的第j个物理坐标值，n为获取同一标记点的物理坐标值的次数。

4.一种标记点配准装置，其特征在于，所述装置包括：

第一坐标获取单元，用于分别多次获取各个物理空间标记点的物理坐标值和各个图像空间标记点的图像坐标值；

第二坐标计算单元，用于分别计算各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值；

空间定位误差计算单元，用于根据多次获取的各个物理空间标记点的物理坐标值和所述物理坐标值的平均值计算出各个物理空间标记点的空间定位误差；

权重矩阵计算单元，用于根据所述空间定位误差计算出各个物理空间标记点的权重矩阵；

配准单元，用于根据各个物理空间标记点的物理坐标值的平均值和各个图像空间标记点的图像坐标值的平均值，结合所述权重矩阵，在保证标记点配准误差FRE最小的前提下，用迭代法求解平移向量t和旋转矩阵R。

5.如权利要求4所述的装置，其特征在于，所述权重矩阵计算单元包括：

方差计算模块，用于计算出各个物理空间标记点的空间定位误差在各个方向上的方差分量以及各个物理空间标记点的空间定位误差的协方差矩阵；

主轴计算模块，用于对所述协方差矩阵进行奇异值分解，得到各个物理空间标记点的空间定位误差的主轴；

权重矩阵计算模块，用于根据所述主轴以及各个物理空间标记点的空间定位误差在各个方向上的方差分量，计算出各个物理空间标记点的权重矩阵。

6.如权利要求4所述的装置，其特征在于，所述第二坐标计算单元包括：

坐标值划分模块，用于将多次获取的各个物理空间标记点的物理坐标值按坐标轴的方向划分为三个分量，分别记为x_ij、y_ij、z_ij；

平均值计算模块，用于计算每个物理坐标值在所述三个分量上的平均值，分别记为

具体计算公式如下：

${\stackrel{\‾}{x}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{x}_{\mathrm{ij}}$

${\stackrel{\‾}{y}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{y}_{\mathrm{ij}}$

${\stackrel{\‾}{z}}_i=\frac{1}{n}\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{z}_{\mathrm{ij}}$

其中，x_ij表示标记点i在x方向上的第j个物理坐标值，y_ij表示标记点i在y方向上的第j个物理坐标值，z_ij表示标记点i在z方向上的第j个物理坐标值，n为获取同一标记点的物理坐标值的次数。

7.一种外科手术导航系统，其特征在于，所述系统包括权利要求4至6任一项所述的标记点配准装置。

Images