[发明专利]机器人辅助外科手术空间映射方法

权利要求书

1、一种机器人辅助外科手术空间映射方法，所述的机器人辅助外科手术空间包括图像空间、机器人空间和手术时病人空间，

摄像头空间：是以双目视觉摄像头，上两个摄像头的中间点为原点，两个摄像头的连线为X方向，两个摄像头镜头所在的平面为XY平面，垂直于该平面的轴为Z轴，建立的坐标系；

图象空间：由以第一张CT/MR片的左下角为圆点，CT/MR的高度方向为Y轴，宽度方向为X轴，CT/MR序号增加的方向是Z轴，建立的坐标系；

机器人空间：以机器人上贴者的摄像头识别标志建立起来的XYZ坐标系，该标志贴在摄像头易于捕获的位置；

手术时病人空间：是指手术时病人所在的物理空间；

其特征在于该方法包括以下步骤：

(1)标识点的建立：利用球体的几何特性，将球体作为定位球，而将其球心作为标识点，标识点的个数为3个或者3个以上；

(2)建立图像空间中仿射坐标系：

(2.1)在具有定位球截面的CT/MR图像上的定位球的截面内部用鼠标选取一个点，同时设定检测区域，通过边缘检测方法找到所有定位球边界点；

所述的检测区域是以选取的点为中心、1.5倍定位球直径折算成像素点数为边长的矩形；

(2.2)采用圆的拟合方法，利用找出的定位球边界点求出定位球在该CT/MR截面上的圆心坐标值；

(2.3)根据已知的定位球的直径，算出定位球的球心到该张CT/MR片的距离ΔZ，而定位球球心的Z坐标等于圆心Z坐标加上或减去ΔZ，即Z_球心＝Z_圆心±ΔZ，求出定位球球心的坐标值；

所述的加上或减去ΔZ是通过判断前后两张CT/MR片在所述的矩形区域内CT/MR灰度总值的大小，灰度总值大的CT/MR片的方向就是定位球球心所在的方向，总值大时加ΔZ，总值小时减ΔZ；

(2.4)利用上述方法找出全部以定位球球心作为标识点的坐标，并以其中的一个标识点作为仿射坐标系的原点，该点与其他标识点的连线作为仿射坐标系的基向量，建立标识点在图像坐标系下的仿射坐标系，得到标识点的仿射坐标矩阵A_I。A_I是3×3的矩阵，矩阵的列是各个基向量；

(3)建立摄像头空间中仿射坐标系：

(3.1)在摄像头坐标系下捕捉到3个标识点及两个长短矢量在摄像头坐标系下的坐标值和矢量，叉乘出垂直于该3个标识点构成的标识面的矢量；

(3.2)通过预先将定位球置于摄像头可识别的标识下的某一位置，根据向量加法计算出定位球球心，从而获得贴在人体身上的其他标识点在摄像头下的坐标，以其中的一个标识点为坐标原点，建立标识点在摄像头空间中的仿射坐标系和仿射坐标矩阵A_C；

(4)机器人辅助外科手术空间之间的相互映射：

(4.1)贴在人体身上的标识点彼此没有相对运动，在摄像头坐标系及图像坐标系下，其仿射坐标矩阵完全相同，通过获得每个仿射坐标系的在其基础坐标系的仿射矩阵和原点坐标，实现空间之间的相互映射；

上述空间相互映射的实现过程如下：

式(1)和式(2)分别是摄像头坐标系和图像坐标系下任意一点在标识点仿射坐标系下的仿射坐标值：

x′y′z′=AC-1xC-xOCyC-yOCzC-zOC---(1)]]>

x′′y′′z′′=AI-1xI-xOIyI-yOIzI-zOI---(2)]]>

其中(x_C，y_C，z_C)、(x_I，y_I，z_I)分别为摄像头仿射坐标系下的任意一点和其在图像坐标系下的对应点，(x_OC，y_OC，z_OC)、(x_OI，y_OI，z_OI)分别为摄像头坐标系和图像坐标系下作为原点的标识点在其两个基础坐标系下的坐标值，(x′，y′，z′)为(x_C，y_C，z_C)以摄像头坐标系为基坐标系建立的仿射坐标值，(x″，y″，z″)为(x_I，y_I，z_I)以图像坐标系为基坐标系建立的仿射坐标值，(x′，y′，z′)与(x″，y″，z″)理论上完全相同(因为是这是不同坐标系下的同一个点)，合并式(1)和式(2)中的两个矩阵得式(3)：

xIyIzI=AIAc-1xc-xocyc-yoczc-zoc+xoIyoIzoI=AIAc-1xcyczc-AIAc-1xocyoczoc+xoIyoIzoI---(3)]]>

这样即建立了(x_C，y_C，z_C)与(x_I，y_I，z_I)的对应关系，其映射过程是这样的：首先获得手术时病人空间的点在摄像头空间坐标系下的点，将点的坐标值转换成以摄像头空间坐标系为基坐标系建立的仿射坐标值，该值也是以图像坐标系为基坐标系建立的仿射坐标值，再将该值转换到图像坐标系下，即完成了图像坐标系与手术时病人空间坐标系的映射，从而实现图像空间与手术时病人空间的相互映射；

(4.2)通过上述方法实现机器人空间与图像空间之间的映射，即在机器人基座上某固定位置放置贴在机器人平台上的标识，在摄像头坐标系下捕获到该标识的仿射矩阵A_RC，该标识的长短轴XZ、XY标识均是已知，将XZ和YX分别作为X轴和Y轴，X轴叉乘Y轴作为Z轴，X点作为原点，建立仿射坐标系得到仿射矩阵A_R，套用式(1)、式(2)得到式(4)：

xRCyRCZRC=ARCAR-1xR-xORyR-yORzR-zOR+xORCyORCzORC---(4)]]>

其中(x_R，y_R，z_R)、(x_OR，y_OR，z_OR)、(x_ORC，y_ORC，z_ORC)、(x_RC，y_RC，z_RC)分别为机器人坐标系下任意一点、机器人坐标系标识点原点、机器人坐标系标识点原点在摄像头坐标下的坐标、机器人坐标系下任意一点映射到摄像头坐标系下的坐标，再带入式(3)中，得到式(5)：

xRIyRIzRI=AIAc-1ARCAR-1xR-xORyR-yORzR-zOR+xORCyORCzORC-xoc-yoc-zoc+xoIyoIzoI---(5)]]>

其中，(x_RI，y_RI，z_RI)为机器人坐标系下的点映射到图像坐标系下的坐标，完成机器人空间与图像空间之间的相互映射，这个过程相对以上映射，增加了机器人坐标系下的坐标映射到摄像头坐标的过程。