毕业论文-导航机器人系统标定实验

导航机器人系统标定实验 摘 要 三维超声成像系统以其形象直观的可视效果以及由此带来的临床价值， 受到了广大研究人员和医疗工作者的关注。 而基于磁定位系统的三维超声成像系统是近几年来三维超声成像的研究热点。 本文主要针对其实现过程中遇到的标定问题进行深入的探讨和研究。 在基于磁定位系统的三维超声成像系统中， 为确定磁定位系统的接收器与超声探头平面之间的精确定位关系， 必须进行磁定位系统的标定。 本文主要探讨了磁定位系统及标定模型中的各种坐标系， 研究了各种坐标系之间转换的实现方法。 针对三维超声系统， 在分析了标定模板的优缺点后， 采用了 “ N” 型模板标定方法及其相应的奇异值分解法。 标定算法能自动提取特征点， 提高了标定精度。 本文采用 “ N”型模型的几何特性来计算标定参数的方案， 不仅减少了超声束在远场发散对标定精度的影响， 而且标定方法简单， 提高了标定效率。 经过实验表明本文利用的这种计算标定参数的方案具有较高的准确性。在标定软件的设计方面，引入了 Matlab 工具，利用 Matlab 强大的数学计算功能及其与 Visual C+灵活的接口编程，设计了计算标定矩阵的软件，为快速、有效地进行磁定位系统的标定提供了有利的工具。 关键词 ： 三维超声，坐标系转换，标定，磁定位 导航机器人系统标定实验 Abstract 3D ultrasound imaging system has potential application in clinical evaluation of diseases.And 3D ultrasound imaging system based on magnetic orientation system is a researchfocus of 3D ultrasound imaging in recent years.In this dissertation,several key problems， including calibration and images capturing are in-depth discussed and investigated. A key requirement of 3D ultrasound imaging system based on magnetic orientation system is calibration:determining the position and orientation of the scan-head with respectto the receiver mounted on it .In this dissertation,coordinate systems in the magnetic orientation system and phantom are discussed and methods of coordinate systems conversion among them are presented,which provides the basis in theory for calibration accurately.For system of 3D ultrasound,a N-shape phantom and a least-squares based on singular value decomposition algorithm are presented in this dissertation.The method not only reduces the effect of the thickness and divergence of ultrasound beam,but also improves the test efficiency greatly.The experiment results show that the new methodhas good performances in precision and accuracy.In design of calibration program,a programming method using mixed programming with Visual C+and Matlab is presented,which makes full use of the advances of Matlabs strong computing ability and flexible interfere with Visual C+. Key words： 3D ultrasound; conversion of coordinate system;calibration;magnetic orientation system 导航机器人系统标定实验 目 录 第一章 前言 .1 1.1 超声引导微波消融治疗肝癌技术 . 1 1.2 肝肿瘤消融机器人系统原理 . 2 1.3 医疗机器人发展现状 . 4 1.4 相关技术研究 . 6 1.4.1 图像的 三维重建 . 6 1.4.2 重建中 的坐标系与坐标转换 . 7 1.4.3 标定模 板 . 8 1.4.4 奇异值 分解法 . 11 1.4.5 图像的 匹配 .12 1.5 本论文的研究内容 . 14 第二章 导航机器人标定实验原理、设备 .15 2.1 实验原理 . 15 2.2 实验设备 . 16 2.2.1 实验 仪器 .16 2.2.2“N” 型标定模板 .20 第三章 导航机器人标定实验 .22 3.1 探针标定 .22 3.1.1 实验步 骤 .23 3.1.2 实验数 据 .23 3.2 模板相对发射器标定 .25 3.2.1 实验步 骤 .26 3.2.2 实验数 据 .27 3.3 图像相对于接收器标定 .29 3.3.1 实验步 骤 .30 3.3.2 实验数 据 .32 3.4 误差分析 .35第四章 结论与展望 .37 导航机器人系统标定实验 4.1 结论 .37 4.2 对进一步研究的展望 .37 参考文献 .39 致 谢 .41 附 录 .42 声 明 .48 导航机器人系统标定实验 第一章 前 言 1.1 超声引导微波消融治疗肝癌技术 肝癌为常见恶性肿瘤之一， 目前每年全世界肝癌发生病人至少为一百万， 且肝脏为仅次于淋巴结的其它实体癌的好发转移部位1。 肝脏常为结直肠腺癌的转移部位。 原发性肝癌以及肝脏转移癌的治疗为目前困扰外科医生的难题之一， 全身及局部的化疗很难取得较好疗效， 相反， 其毒副反应严重的影响了病人的生存质量。 迄今为止， 为提高病人生存率及生活质量， 手术切除仍为治疗肝癌的最好方法。 手术切除肝细胞肝癌，结直肠癌以及其它原发肿瘤肝转移灶时，可提高 25%-30%病人的生存率， 对一些神经内分泌肿瘤肝转移病人， 手术切除肝脏转灶可明显缓解因过量激素分泌水平所至的症状， 由于很多神经内分泌肿瘤和结直肠癌肝转移的病人可得到治愈性切除， 肿物局限于肝脏并不是手术切除的绝对适应症， 肝癌的多中心起源， 肿物临近肝内大的血管及胆管， 并存的肝硬化而致的肝功能不全均不能保证根治性切除。 必须探寻出新的治疗方法。 现在应用较多的是局部消融技术， 包括经皮肤无水酒精注射、 微波凝固、 间质激光凝固、 冷冻及射频消融2。 射频消融治疗肝脏肿瘤于 1995 年由意大利的 Rossi率先应用于临床， 为当今最新的肝肿瘤导向治疗方法之一，通过产生热能而破坏肿瘤细胞，当肿瘤细胞加热超过 45-50 摄氏度时，细胞内蛋白质变性， 细胞膜双脂质膜溶解导致细胞膜破坏， 从而细胞结构改变， 射频消融技术通过使细胞内温度超过 60 摄氏度而达到破坏肿瘤细胞的目的3。因此 ,发展有效的非手术疗法是临床迫切需要解决的问题 .其中介入治疗被公认为是治疗不可切除的中、晚期肝癌的重要手段 .超声引导微波消融治疗肝癌技术就是介入治疗的一种典型代表 .超声引导微波消融治疗肝癌过程中病人无需开刀 ,在局部麻醉后 ,医生利用超声图像引导 ,将针型微波辐射天线经皮肤穿刺进入肿瘤中 ,利用微波加热将肿瘤消融 ,达到治疗效果。 超声引导微波消融治疗肝肿瘤是热消融的一种方式。 所谓肿瘤的热消融治疗方法是指在影像技术的导航下应用化学的或加热的方法直接作用于病灶实体肿瘤， 根除或毁坏肿瘤组织的治疗方法。 热消融治疗在临床应用的主要方法有： 射频、 微波、激光及高强聚焦超声等4。 在治疗过程中患者无需开刀， 在局部麻醉后， 医生利用- 1 - 导航机器人系统标定实验 超声图像导航将针型微波辐射天线经皮肤穿刺进入肿瘤中， 利用微波加热将肿瘤消融达到治疗效果。其工作原理如图 1-1 所示。 和传统的手术相比 ,微波消融治疗肝癌技术存在着明显的优点 :微创 ,安全 ,无化疗、 放疗的毒副作用 ,灭活效果好 ,对正常组织损伤小 .正是由于上述优点 ,在超声引导下的针对肿块的各种局部热消融的治疗方法在临床上普遍开展 ,显示出强劲的发展潜力。 图 1-1 超声引导微波消融技术 1.2 肝癌消融机器人系统原理 图 1-2 为基于超声影像导航的肝癌消融机器人系统原理图。 患者腹部肝脏病变部位经过 FreeHand（自由手臂扫描） 三维超声成像系统扫描后输入到计算机中的手术规划软件， 医生借助手术规划软件完成病灶及周边重要血管的三维重建， 建立术前三维模型坐标系、 设定靶点、 规划路径、 模拟热场。 手术时， 由于术前和术中两个阶段患者在发射器中的位置会发生变化， 因此需要通过肝脏血管的术前和术中脊线模型配准，将术前的模型映射到患者实体上使得术前的模型与病人实体一一对应。 最后， 再将术前的模型映射到机器人坐标系中， 同时将三维模型中设定的靶点和穿刺路径等手术规划信息映射到机器人坐标系中， 作为机器人位置和姿态的目标- 2 - 导航机器人系统标定实验 值， 控制机器人运动， 具体的坐标变换原理详见 1.4 节。 机器人运动完成后， 末端执行器上的导向装置可辅助医生对患者病灶进行定位， 对手术器械进行导引， 完成一系列的手术操作。 术后， 再次扫描患者图像， 评价手术效果。 根据该原理， 不难发现基于超声影像导航的肝癌消融机器人系统主要由手术规划软件、 定位装置、 超声仪、机器人及控制器组成5。 图 1-2 系统原理图 根据上述原理，典型的机器人辅助肝癌消融外科手术流程如下： (1) 手术前， 医生通过 FreeHand 三维超声成像系统， 得到患者肝脏病变组织、 肝脏血管的超声医学图像及其位姿信息。 (2) 把超声医学图像及其位姿信息输入到手术规划软件，该软件首先将非规则的超声医学图像重建成标准的断层图像。 然后医生借助手术规划软件完成病灶区域轮廓的勾勒， 再由手术规划软件重建出病灶区域的三维模型。 在此基础上， 医生进行靶点设定、路径规划、热场模拟等手术规划。 (3)手术时，医生通过 FreeHand 三维超声成像系统，得到患者术中肝脏血管的超声医学图像信息。 在此基础上， 通过提取术前、 术中血管的脊线模型完成术中患者实体和术前的三维模型的空间映射， 从而获得术前三维模型坐标系到机器人坐标系的映射。 (4)医生在手术导航软件中所做的手术规划将被映射到机器人坐标系， 该信息被用来控制机器人运动，使消融针按照预定的路线到达患者体内的病灶点。在此过程中，获取实时的 B 型超声图像以监控手术过程，防止消融针偏离预先指定的路径。 (5)消融针到达指定靶点后，医生对患者体内的病灶进行一系列的处理。 - 3 - 导航机器人系统标定实验 (6)手术后，医生通过 FreeHand 三维超声成像系统，得到患者肝脏组织超声医学图像通过计算机辅助分析判断手术效果。 左图 1-3 为常用射频电极， 此电极较适合于 B 超引导下操作。 右图 1-4 是改进的射频电极，较适合 CT 引导下操作。 图 1-3 B 超射频电极 图 1-4 CT 射频电极 1.3 医疗 机器人发展现状 当前医疗机器人系统一种可能的分类情况如图 1-5，该系统还可以用来确定未来几年外科手术机器人的发展趋势。 一个外科手术机器人系统至少需要包含以下几个系统： 规划系统 导航系统 安全系统 机器人进行外科手术时所需的一些特殊需求： 消毒 保证病人及医疗人员的安全 - 4 - 导航机器人系统标定实验 表 1.5 外科手术机器人的分类 手术类型 传统手术 微创手术 创内手术 自主系统 ROBODOC CASPAR Slereotaxis Inc MUSYCEMIL 交互式系统 Eye scalpelRinC AESOP MIAS Active Catheters 遥控操作系统 Mammotome PAKY da Vinci ZEUS MINOSC 交互类型 被动系统 PinPoint HALS (not robotic)(非机器人) Given Imaging (not robotic) (非机器人) 美国 Computer Motion公司开发了 Zeus微创外科机器人系统，该系统结构如图1-6 所示，适用于内窥镜微创手术。系统由三只置于手术床上的交互性机械臂、计算机控制器及医生控制台组成。 手术时， 医生坐在控制台前， 通过摄像机观察手术情况， 用语音控制内窥镜， 并通过操作手柄来控制手术工具。 其中一只机械臂采用声控的交互方式，操作内窥镜，另外两只机械臂则在医生控制下操作手术工具6。 图 1-6 zeus 微创外科机器人系统 从上述资料看来，医疗机器人系统的研究主要集中在西方发达国家。我国人口众多， 正在逐步成为一个医疗器械的生产大国和使用大国。 然而， 我国目前现代医疗的整体水平和高新技术的医疗器械研究与开发还十分落后。 随着我国社会文明- 5 - 导航机器人系统标定实验 和科学技术的进步， 人民生活水平的提高， 人们对疾病的诊断、 治疗、 预防以及卫生健康也提出了越来越高的要求。 我国政府有关部门已经开始关注医疗机器人的研究与开发。 国内清华大学、 北京航空航天大学、 上海交通大学、 哈尔滨工业大学等在医疗外科机器人方面开展了初步的研究工作， 已经取得了一些可喜成果。 上世纪90 年代中后期，北京航空航天大学机器人研究所与清华大学、解放军海军总医院合作， 在国内开展了远程脑外科机器人辅助手术系统的研究， 研制成功了具有临床实用性的机器人辅助脑外科立体定向手术系统 机器人技术在微创外科领域的开发与应用将从根本上改变 21 世纪外科手术的面貌， 推动医疗领域新的技术变革。 正如该领域中的著名学者美国 Taylor 教授的评价， 机器人和计算机技术对于外科手术的影响将如同计算机集成制造对于产品制造的影响一样广泛而深远。 1.4 相关技术研究 1.4.1 图像的三维重建 由于二维超声只能反映出人体三维结构的平面图像， 所以超声诊断医生需要通过观察不同方向的二维切片， 然后在大脑中对图像信息整合， 才能够想象出人体器官的三维立体情景。 这常常受到器官的内部结构、 切面成像等诸多因素的制约， 所以二维超声不能全面直观地显示出器官结构的真实形态。 而三维超声可以比较直观的显示复杂器官的立体形态。 例如在肝癌介入治疗中， 三维肿瘤模型能够给医生提供更加丰富的信息。 它能在计算机屏幕上， 用不同颜色显示不同解剖模型， 还可以从一些极为特殊的角度对器官进行观测，这有助于医生术前更加全面地了解病情。由于三维肿瘤模型提供了肿瘤在体内精确的空间位置以及三维形态， 可以用来进行手术规划和模拟， 病灶分析和测量， 进而提高体外超声治疗和超声导向介入治疗手术的准确度，因此被广泛引用到手术导航中来。 目前， 有三种方法获取三维超声7： 第一种方法是利用二维面阵探头发射体积超声束从而获得实时的三维空间数据的方法； 第二种方法是利用马达带动超声探头进行匀速地扫描从而获得三维空间的数据，该方式可以旋转扫描也可以扇形扫描；第三种方法是 FreeHand三维超声方法， 该方法利用附属在超声探头上的定位装置来记录当前二维 B型超声图像在三维空间中位置从而获得三维空间的标准断层图像数- 6 - 导航机器人系统标定实验 据，如 所示。大部分商业系统考虑到操作的方便性和实用性，选择马达带动超声探头进行均匀扫描获得三维超声， 但是 FreeHand三维超声技术具有体积小， 任意角度扫描，甚至可以进行不连续扫描的优点而成为手术导航过程中的首选。 图 1-7 三维超声扫描方式 应用 FreeHand 三维超声的前提是标定出超声平面图像坐标系和附属在超声探头定位装置坐标系之间的转换矩阵。 三维超声标定的一般思路是用粘贴有传感器的超声探头对一个空间位置及几何特性已知的模板做多次扫描，得到一组测量方程，然后通过求解标定方程的方法进行标定。 1.4.2 重建中的坐标系与坐标转换 在三维重建过程中，共涉及到 3 个笛卡儿坐标系。 1.空间绝对坐标系 T,由电磁定位系统中的发射器确定在重建过程中始终保持不变 2.超声图像平面坐标系 P 3.电磁定位系统中的接收器的坐标系 R 接收器固定在超声探头上因此接收器坐标系与超声平面坐标系有固定的转换关系。 首先需要确定接收器坐标系 R 与发射器坐标系 T 间的关系根据空间几何性质 从 R 坐标系统中一点 转换到 T 坐标系统中的位置pRPT可用公式 (1)表示 13333333 XTOXRXRTXTRPxMP += 其中 为 R 坐标系统的原点到 T 坐标系统位置矢量 , 为 R 坐标系统的坐标轴与 T 坐标系统的坐标轴之间的旋转矩阵。 在磁定位系统中， 已提供了 6个自由度信息，它们是 3 个旋转信息XRTO333XRTM , 以及 3 个平移信息 x, y ,z。 利用这 6 个自由度信息我们就可以确立 和 即为： 33XRTT13XTOR- 7 - 导航机器人系统标定实验 33XRTT =8+coscossincossinsincoscossinsincoscossinsinsincossinsinsincossincoscossinsinsincoscoscos=xyz 13XTORTPR点先经过旋转变换 ,再进行平移从而转换到 T 坐标系统的 PT点。 将位置偏移量和旋转矩阵合起来也可以写成一个 4 4 齐次坐标矩阵 。 RTTRTT = 10001333 XTOXRTRM于是 P 坐标系统的点 PR转换到 T 坐标系统中可以写成行齐次矩阵乘法的形式PT=RTT PR。我们进一步就可以推出超声平面坐标系上的点 PP转换到 T 坐标系统的关系 : PT=RTT PRT PR。其 中 是接收器坐标系与超声平面坐标系间的转换关系，即要通过校准得到的目标。 PRT由此可见，任意两坐标系 A,B 之间的转换关系均可以表示为： A=旋转矩阵 平移向 量+B 1.4.3 标定模板 校准一般是通过扫查一个已知几何尺寸的模型来进行的。 利用模型的几何特性和磁定位系统的测量来确定校准参数。 1.交叉线模型9一般的标定模型是交叉线标定模型，图 1-8a。两条交叉的棉线被固定在一个水槽中， 超声探头从不同的方向扫查得到交叉点的图像， 图 1-8b。 利 用交叉点相对于磁定位系统中发射器的位置不变的特性来计算标定参数。 为了计算方便， C 坐标系的原点并不是和重建的坐标一致，而是放在交叉点处。 - 8 - 导航机器人系统标定实验 图 1-8a 交叉线模型 图 1-8b 交叉线模型扫描图像 这种标定方法也可以采用其他的点物体 （比如针头或者是一个悬挂着的小铜珠）来代替交叉点。 这种标定方法的精度主要依靠准确定位交叉点的中心。 一般都是通过手动的方法来确定在 B 超图像中的位置，这使得标定的过程枯燥而且耗时。 2.三交线模型10另一种标定技术是利用三交线模型， 图 1-9a。 这种方法类似于交叉点的标定方法。 因为它也采用在一个水槽中扫描交叉线的方法， 但是它要求三条线必须精确定位在三个相互垂直的方法上，图 1-9b。 图 1-9a 三交线模型 图 1-9b 三交线模型扫描图像 这种方法的优势就是容易操作， 然而要求三条线必须严格垂直， 标定时也必须对于三条线分别扫查，可操作性较差。 3.单平面模型11- 9 - 导航机器人系统标定实验 当扫查交叉线和三交叉线模型时，我们可以观测到在模型装置的底端会在 B超图像上超声清晰连续的直线， 所以我们可以单独采用水槽底面来进行标定， 如图1-10a。标定的精度主要决定于水槽底面的平稳程度，图 1-10b。 图 1-10a 单平面模型 图 1-10b 单平面模型扫描图像 这种通过扫查水槽底面的标定方法比较简单， 不需要特别的标定模型。 而且由于采用平面法可以在超声图像上产生一条清晰的直线， 这条直线可以通过程序来自动识别。 4.剑桥模型12单平面标定方法由于受到超声束的厚度和光谱反射特性的影响， 因为当探头以倾斜角度扫查时，这种结果就会引入位置误差。于是英国学者 R.W.Prager 等人于1998 年设计了一种称为剑桥模型的标定装置，这个模型主要包括两个部分：一个夹住探头的夹子和一个加在两个圆盘中间的细铜条，图 1-11a。夹子限制细铜条只能在超声束的中间来回运动， 这样便有效的避免了声束厚度对标定精度的影响， 并且大大节约了测量时间，图 1-11b。但是与上面的方法相似它也涉及复杂设计和精密加工的问题。 并且， 不同的探头需要不同的夹具。 需要重新确定声束的中间位置。 - 10 - 导航机器人系统标定实验 图 1-11a 剑桥模型 图 1-11b 剑桥模型扫描图像 5.三棱锥标定模型13清华大学的刘洁等人设计的三棱锥装置 （如图 1-12） 在 一个塑料水槽中放置一个三棱锥模型通过反复扫查棱锥得到超声平面与棱锥的交点并计算出标定矩阵基磁定位系统的三维医学超声系统的图像采集与标定方法研究。 这种方法不对固定目标点进行扫查而是任意截取三棱锥的平面从而增强了实验的灵活性而且制作简单但精度不高。 图 1-12 三棱锥模型 1.4.4 奇异值分解法 设有一一对应的两点集 i=X x ， i=Y y 。其 中 ，ixiy 分别表示坐标系 X 和 Y中第 i 个点的三维位置坐标。 现需要对 和 Y 进行配准。 假设它们之间的转换关系为刚体变换，即： X- 11 - 导航机器人系统标定实验 ii= +y Rx p (1) 其中 R 是 旋转矩阵；33 p 是 13 的平移向量。 配准的目标是找到最优旋转矩阵 R 和平移向量 p ，使以下目标函数最小： 211(,) ( ( )NiiifN= +Rp y Rx p (2) 1.4.5 图像的匹配 在医疗机器人系统中， 为了按照三维模型坐标系中的手术规划引导机器人准确到达病灶靶点，必须将三维模型中的靶点位置通过空间映射传递到机器人空间中。因此， 求得三维模型坐标系与机器人坐标系之间的坐标映射关系是微创外科机器人系统的关键技术。该映射的求解涉及到术前模型和术中患者实体之间的图像配准。配准的主要任务是将术前的患者模型和术中患者实体信息统一在一个公共的坐标系中。 在统一坐标系的过程中， 一般是以矩阵变换的方式来进行坐标空间之间的转换。 它是坐标映射中最为关键的步骤之一， 其精度直接决定了坐标映射乃至系统的精度。 按照配准所基于的特征的不同来源， 可以分为基于体外特征的配准和基于体内特征的配准。 基于体外特征的配准分为有创配准和无创配准两类。 有创配准的包括神经外科手术过程中的立体定向框架13及植入标志点的方法14，立体定向框架方法需要用螺钉拧到患者的头骨上进行定位， 给患者带来了极大的痛苦。 因此， 研究者开始设计一些无创的配准方法，如利用模具、框架和牙托等方法。 Schad等人使用一个特殊的面具来固定框架，配准误差不超过 2mm15。利用体表标志点的方法也属于该类方法， 由于该方法操作简便， 因此在不同形式的医学图像配准和医疗机器人系统中得到了广泛的应用。 Vinas等人开发的无框架脑外科机器人就属于此类16。 骨科和神经外科手术过程中，由于骨骼不变形，使得术前获取的图像和术中患者实体之间很容易通过体表或骨骼上的标志点进行刚性配准， 从而可以较容易地按照术前的计划进行手术。 但是一些软组织容易变形、 移动， 必须通过体内的一些稳定的特征将术前获取的模型与手术中患者实体之间进行配准。 基于体内特征的配准主要有以下三类： 基于标记的配准、 基于分割的配准、 和基于像素的配准。 基于标记分为基于解剖知识的标记17（如： 利用人体特殊的解剖结构， 一般由人工直接描- 12 - 导航机器人系统标定实验 述） 和基于几何知识的标记18。 19（如： 运 用数学知识得到大量点、 线、 面的曲率、角度特征等） ； 基于分割的配准分为基于刚体变换的配准20（如： 头骨的表皮分割）和基于变形模型的配准21（心脏的表皮分割） 。 基于标记和分割的配准方法运算速度比较快（在分割完成之后），然而这两种方法的主要缺点在于： 1、需要对患者植入人工标记； 2、配准依赖于分割，而分割复杂度高，很容易出错，所以通常需要交互操作。 基于像素的配准分为部分运用像素灰度信息22（如： 对像素二值化处理后的基于矩的配准） 和全部运用像素灰度信息 （如： 基于傅立叶变换23或者熵的配准24） 。 基于像素灰度信息的配准直接对像素灰度进行操作， 与基于分割的配准方法相比， 虽然需要对应的点数增加了， 但是不再需要特征提取的步骤， 也就是说，不需要进行图像分割。 在肝脏介入治疗的计算辅助导航系统中， 术前模型和术中模型的配准也取得了很好的发展。 S.Eulenstein25在基于超声导航的肝脏手术培训系统中，首先将术前CT图像中的血管部分利用区域增长算法分割出来，然后利用骨骼化算法提出出血管的脊线。 术中的三维超声数据通过多谱勒能量图提取出脊线。 配准过程分为三个步骤： 第一步， 通过人工选择 3-4 个点， 用刚性配准实现粗配准； 第二步， 在粗配准的基础上采用迭代最近点算法 (Iterative Closed Point， ICP)再次进行刚性配准； 第三步，使用 B样条曲线进行非刚性配准。 P.Windyga26在基于超声导航的计算机辅助手术中采用了两步配准方法： 第一步粗配准采用肿瘤区域的力矩和主轴进行粗配准；第二步用基于血管标志点进行精配准。 Aylward27在导航系统采用了一种基于图像和特征的配准方法， 该方法只要求在术前图像中提取特征点， 然后与术中的图像直接进行配准， 无须再进行特征提取。 该方法避免了对术中图像的分割， 极大地提高了配准效率和精度。 由于配准的精度直接影像系统的精度，因此不少学者也对配准精度展开了研究， 目前对于配准精度的研究主要集中在基于标志点配准精度的研究。 Sibson首先给出了基于标志点配准误差的统计分布模型28。在此基础上， Fitzpatrick利用六个均匀分布的标志点实验平台进行配准精度的实验29。 Fitzpatrick还推导了基于标志点配准方法中靶点映射误差的公式30。 Pennec等人利用最小二乘原理对基于标志点配准获得的变换矩阵进行了误差分析31、 33。 - 13 - 导航机器人系统标定实验 1.5 本论文的研究内容 本文以基于超声影像导航的肝癌消融机器人为代表， 研究医疗机器人系统的标定问题。 标定问题的实质是坐标系之间的转换。 标定之所以困难是由于超声探头和接收器都是密封装置 ,从外部无法直接测量两者的空间关系。因此 ,通常借助于某种标定模型（“ N”型线框） ，通过多次 测量这个模型，并结合其已知的物理特性，计算出接收器与图像平面间的空间关系。 三维超声标定的一般思路是用固定有传感器的超声探头对一个空间位置及几何特性的已知的模板做多次扫描，得到测量方程，通过求解方程得到转换关系。 标定依据： FIRRTTPFITTTP =标定过程： 首先进行了探针标定， 之后进行模板标定， 最后进行了超声探头的标定， 得到了图像上一点 F 与模板坐标系 P 之间的转换关系。 分析了误差产生的可能原因，对进一步精确测量提供了理论依据。 - 14 - 导航机器人系统标定实验 第二章 导航机器人系统标定实验原理、设备 2.1 实验原理 要进行三维重建和模型配准， 必须标定出固定在超声探头上的接收器坐标系与超声平面坐标系之间的转换关系。 在实验前， 有必要了解一下标定系统中存在的 4组坐标系。 标定系统中存在着四组坐标系： 图像坐标系 I 、 磁定位仪接收器坐标系 R 、磁定位仪发射器坐标系 T 、 以及模板坐标系 ， 如图所示。 设模板上有一点 处于某幅图像的视野当中， 则点 在坐标系 和坐标系P FF P I 中的坐标转换关系以齐次坐标变换矩阵表示为： FIRRTTPFITTTP =式中： 01TFu vsusv= I 为点 在超声图像平面中的坐标， 、 为超声图像两个方向的分辨率， 一般情况下， 。 为 了避免超声探头深度调整而造成图像坐标原点位置变化，将图像坐标原点设置在扇形探头曲面的中点； 为图像坐标系FusvsusvsRITI 相对于磁定位仪接收器坐标系 R 的转换矩阵，包含 3 个旋转参数、 3 个平移参数。 TRT 为磁定位仪接收器坐标系 R 相对于磁定位仪发射器坐标系 T 的转换矩阵，由 磁 定 位 仪 读 数 求 解 得 到 ；P为发射器坐标系TT T 相对于模板坐标系 的转换矩阵，未知待求解。 P标定的目的是求解转换矩阵 ， ，以及超声图像分辨率 和 。标定的步骤是先标定 再标定 ，其中， 标定又分为两步：第一步求出探针针尖在磁定位仪发射器坐标；第二步标定模板空间与磁定位仪发射器的转换矩阵。 RITPTTusvsPTTRITPTT- 15 - 导航机器人系统标定实验 图 1-13 坐标系定义及转换 2.2.实验 设备 2.2.1 实验仪器 标定实验所需硬件包括： 磁定位仪 （发 射器固定位置） 、 B 超机、 计算机、水槽（灌满水）、探针（一根）、N 型 模板。实验平台如图 2-1 所示 图 2-1 实验平台 - 16 - 导航机器人系统标定实验 1）磁定位仪 本次实验所选用的磁定位仪是 NDI公司生产的介入治疗电磁定位系统 Aurora。NDI的电磁定位系统 Aurora具有无视线遮挡，精度高，抗金属干扰和易于集成等特点， 已成功应用于介入治疗， 内窥镜诊断， 和骨科手术的计算机辅助手术和治疗中，系统的超小型定位线圈可与多种介入工具结合， 在体内提供实时定位和引导， 如心导管， 活检穿刺针， 粒子放疗针， 消融针等工具的引导， 随着计算机辅助手术和治疗技术在临床的进一步推广， Aurora电磁定位技术的应用前景非常广阔34。 系统由磁场发生器，系统控制器，工具联结器和定位线圈组成，主要技术指标包括 : 测量范围 90 厘米 X 60 厘米 RMS 测量精度可达 0.7 毫米 小型 5D 定位线圈 0.55 毫米 直径 8 毫米 小型 6D 定位线圈 1.8 毫米 直径 8 毫米 抗金属对测量的影响 可同时跟踪 8 个定位线圈 定位工具开发软件 市场上还有一些其他品牌的磁定位仪，如美国Ascension Technology公司生产的产品Flock of Birds35。如图 2-2 所示。 图 2-2 Flock of Birds 磁定位仪 - 17 - 导航机器人系统标定实验 2） B 型超声诊断仪 本次试验所选用的 B型超声诊断仪是江苏中科电子生产的 ZK-3000B型B超36，该设备具有 高精度 DSC数字成像技术； ABS注塑一体式机壳； 智能化 TGC全程 8 段增益控制；先进的 SMT贴片工艺确保整机可靠性；独特的 4B显示功能，多方位显示病变信息。 性能介绍： 可配探头：电子凸阵、电子线阵、腔内探头 工作频率： 2.5MHz 7.5MHz （ 2.5MHz/3.5MHz/5.0MHz/6.5MHz/7.5MHz 宽频、变频、依据探头） 显示方式： B 型、 B/B 型、 B/M 型、 M 型 灰阶： 256 级 显示屏： 14 四种图像倍率： 1.0、 1.2、 1.5、 2.0 及局部放大功能 最大显示深度： 240mm 图像处理：边缘增强、帧相关、伽玛校正（ 8 种）、直方图、左右翻 转、上下翻转、极性翻转、上下提升功能，多达 512 幅以上的电影 回放功能、 32 幅图像存储功能 体位标记： 18 种带探头标记的体位标记 测量项目：距离、面积、周长、体积（椭圆法）、角度、心率、斜率、 时间间隔； GS、 CRL、 BPD、 HC、 FL 等测量胎儿的胎重、妊娠期、预产期 字符标记：自动显示日期、星期、时间、探头频率、帧频、病历编号、 医生姓名、医院名称、倍率、计测值、体位标记、探头标记（可转动）、 全屏幕字符编辑、字符注解、中英文显示等 标准配置：主机一台、 2.5/3.5/5.0MHz 宽频、变频电子凸阵探头一只 . 选购探头： 5.0/6.5/7.5MHz 宽频、变频电子 线阵探头、 5.0/6.5/7.5MHz 宽频、变频腔体探头 体积： 560mm（宽） 1180mm（高） 680mm（深） 重量：约 45kg 选购件：视频打印机 可连接：影像工作站 - 18 - 导航机器人系统标定实验 市场上还有一些其他品牌的B型超声诊断仪，比如江苏大为电子设备有限公司生产的 DW3102A数字化B超 （如图 2-3）37。 这两款产品在外形， 功能， 技术参数方面都相仿，从品牌角度出发还是选择了 ZK-3000B型B超机。 图 2-3 江苏大为 B 型超声诊断仪 3）水槽 实验选用的是一个长*宽*高： 402525(CM)玻璃制水槽。 充入介质为纯净水。有 2 个可控条件对实验结果有明显影响。 a.温度 超声仪的计算速度设定为超声波在人体内传播速度v,超声仪所直接测得的数据是超声波到达目标点的时间t,那么目标点位置的计算为v t38。 如果介质的传播速度与人体传播速度 1540m/s不一样,将导致目标点位置的计算不准确,因此需要- 19 - 导航机器人系统标定实验 使标定时水槽中介质与人体传播速度一致。50热水的传播速度为