CN110123451B - 应用于无标记点光学手术导航系统的病人表面注册方法 - Google Patents
应用于无标记点光学手术导航系统的病人表面注册方法 Download PDFInfo
- Publication number
- CN110123451B CN110123451B CN201910306492.XA CN201910306492A CN110123451B CN 110123451 B CN110123451 B CN 110123451B CN 201910306492 A CN201910306492 A CN 201910306492A CN 110123451 B CN110123451 B CN 110123451B
- Authority
- CN
- China
- Prior art keywords
- point set
- matching
- point
- patient
- points
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000003287 optical effect Effects 0.000 title claims abstract description 17
- 230000009466 transformation Effects 0.000 claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 239000004576 sand Substances 0.000 claims description 2
- 238000012216 screening Methods 0.000 abstract 2
- 230000001678 irradiating effect Effects 0.000 abstract 1
- 239000003550 marker Substances 0.000 description 4
- 208000014674 injury Diseases 0.000 description 3
- 210000003484 anatomy Anatomy 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008733 trauma Effects 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
Abstract
本发明公开了一种应用于无标记点光学手术导航系统的病人表面注册方法,包括步骤:1)获取病人图像空间表面数据:原始点集与简化点集;2)获取病人实际空间表面数据:激光点集与简化点集;3)通过网格化和八方位匹配的方法筛选出最优匹配的粗匹配变换矩阵集合;4)在上一步的基础上求取精匹配的匹配矩阵集合和对应的匹配误差。根据误差筛选出精匹配步骤的最佳变换矩阵。5)根据3)和4)可求取总的最佳变换矩阵,完成空间注册。本发明使用红外激光灯照射病人表面获取实际空间的数据空间注册,采用网格化进行寻优,并将待匹配点集以八种不同的方位与目标点集进行匹配,这两步操作避免了匹配结果陷入局部最优的情况,实现了自动注册的功能。
Description
技术领域
本发明涉及光学手术导航的技术领域,尤其是指一种应用于无标记点光学手术导航系统的病人表面注册方法。
背景技术
手术导航系统在临床上越来越广泛的使用,成为外科手术治疗精准化和微创化的重要发展方向。在手术精度、手术消耗时间、手术创伤、术后疗效等方面与对应的传统外科手术相比有长足的改进。在术中导航时,系统可以通过可视化手术工具和患者特征结构的在实际空间中的位置来帮助医生进行手术,从而提高手术精度和减少对病人的伤害。导航系统的关键步骤是实际空间与图像空间的变换关系,这一步被称为注册。注册方法的优劣直接影响手术导航精度和导航的方便性。
目前手术导航中注册方法多采用标记点注册方法。标记点的种类主要包括骨植入螺钉标记点、解剖标记点和粘贴于皮肤表面的标记点。骨植入螺钉的标记点精度最高,但是需要术前植入器械,会对患者造成额外的创伤和疼痛。解剖标记点利用人体比较显著的解剖特征作为标记点,需要医生在术前提取图像空间的解剖标记点,以及在术中提取患者空间的解剖标记点,存在一定的人为操作误差。粘贴于皮肤表面的标记点精度较高并且操作简单,是最常用的一种标记点。但是在使用过程中存在标记点光线易被遮挡、标记点易掉落以及会由于皮肤的位移形变产生移动等问题,从而影响到导航精度。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提出了一种应用于无标记点光学手术导航系统的病人表面注册方法,该方法有效解决标记点光线易被遮挡,标记点易掉落以及标记点粘贴部位受限等问题,避免了手动提取标记点产生误差的情况,同时,注册过程操作简单,注册时间短并且注册精度可达到临床要求。
为实现上述目的,本发明所提供的技术方案为:应用于无标记点光学手术导航系统的病人表面注册方法,包括以下步骤:
1)导入病人CT图像序列,对CT图像序列提取外表面轮廓,记为点集P,将点集P简化为点集P';
2)用红外激光笔扫描病人表面,光学导航仪实时记录病人表面反射的红外光点的三维坐标,将这些坐标保存为原始激光点集Q,将点集Q简化得到点集Q';
3)分别将点集P'作为目标点集,点集Q'作为原始点集进行分网格初始化粗匹配,得到的候选粗匹配变换矩阵的集合M={Mi|i=1,…,k},其中k为网格的个数;
4)结合步骤3)中粗匹配求得的k种粗匹配变换矩阵结果,将点集P作为匹配的目标点集,将点集Q作为匹配的原始点集进行精匹配,求得精匹配中的最佳匹配矩阵M1s和对应的粗匹配变换矩阵Ms,0<s<k+1;
5)根据步骤4)中的结果得到表面注册的实际空间和图像空间的变换矩阵为Mat=M1s*Ms。
在步骤1)、步骤2)中,点集简化方法是:设原始点集为S,简化后的点集为S',简化边长为l,将点集S用边长为l的正方体进行划分,对于每个正方体,求取其中心点坐标,若在该正方体内部有属于点集S的点,则找到这些点中距离正方体中心点坐标最近的点并保存得到简化后的点集S'。
在步骤2)中,使用的红外激光笔包含两束激光,一束是波长为650nm的可见光,用于指示红外光的位置;另一束是波长为980nm的不可见红外光,将其照射到人体表面,红外光学定位仪获取实际空间人体表面信息。
在步骤3)中,将点集P'作为目标点集和点集Q'作为原始点集进行分网格初始化粗匹配求取候选初始位置对应的变换矩阵的集合,包括以下步骤:
3.1)将点集P'用边长为50mm的正方体网格进行划分,得到n个正方体网格,网格中心点保存为点集C={ci|i=1,…,n},求取点集Q'的中心点cq,对点集C中每一个点ci,分别设ci为点集C的中心点,将点集Q'和点集C进行八方位匹配,当ci为点集C的中心点时,点集C和点集Q'的最佳变换矩阵为Mi,将矩阵Mi保存到集合M={Mi|i=1,…,k};
3.2)对于步骤3.1)中提到的设置两点集的八方位匹配,设A为原始点集,B为目标点集,其方式是:平移点集A和点集B,使平移后的点集A'和点集B'的中心点都与坐标系原点重合;旋转点集A'和点集B',使旋转后的点集A1”和点集B”的拟平面法向量与坐标系z轴平行,点集B”为目标点集B的匹配初始位置,点集A1”为原始点集A的第一种匹配初始位置;分别将点集A1”绕z轴旋转90°、180°和270°得到点集A2”、A3”、A4”作为第二、三、四种初始位置;将点集A1”绕y轴旋转180°得到点集A5”作为第五种匹配初始位置;分别将点集A5”绕z轴旋转90°、180°和270°得到点集A6”、A7”、A8”作为第六、七、八种初始位置;分别将点集A1”~A8”与点集B”进行ICP匹配,得到八组变换矩阵和对应的误差,比较误差值,最小的误差对应的变换矩阵则为该八方位匹配步骤的最佳变换矩阵;
在步骤4)中,分别用步骤3.1)中求出来的变换矩阵集合M={Mi|i=1,…,k}中的每一个矩阵Mi对点集P中的坐标进行变换得到点集Pi,将变换后的点集Pi分别与点集Q进行匹配,求得k个匹配误差e={ei|i=1,…,k}和变换矩阵M1={M1i|i=1,…,k},找到最小误差值对应的序号s,0<s<k+1,则精匹配结果求得的变换矩阵为M1s。
本发明与现有技术相比,具有如下优点与有益效果:
1、本发明首次提出了通过将目标点集网格化初始位置匹配和八方位匹配的方式进行匹配寻优,避免在匹配过程中陷入局部最优的情况,从而得到全局最优匹配的结果。
2、本发明通过将患者在图像空间和实际空间的表面信息匹配起来,从而完成空间注册,这一过程中,不需要通过医生手动确定图像空间和实际空间的特征点来完成粗匹配,简化空间注册的操作。
附图说明
图1为本发明方法的流程示意图。
图2为简化的患者图像空间表面点集。
图3为患者实际空间表面激光点集。
图4为将患者图像空间表面点集进行正方体网格化。
图5为八方位匹配中激光点集的八种位置情况。
图6为图像空间点集和实际空间激光点集的最终匹配结果。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
如图1所示,本实施例所提供的应用于无标记点光学手术导航系统的病人表面注册方法,包括以下步骤:
步骤S101,患者术前三维医学图像表面数据的获取。首先是获得CT扫描医学影像序列,然后根据影响序列重建出患者三维解剖结构。用移动立方体算法(MC)从三维解剖结构中提取患者表面数据得到患者原始表面点集P。然后再将P均匀网格简化为点集P'。
步骤S102,患者实际空间表面数据的获取。首先用红外激光笔扫描人体表面,光学定位仪记录患者表面反射的激光点坐标保存为原始激光点集。对记录到的原始激光点集进行去噪得到患者表面激光点集Q。此处使用的红外激光笔包含两束激光,一束是波长为650nm的可见光,用于指示红外光的位置;另一束是波长为980nm的不可见红外光,将其照射到人体表面,红外光学定位仪获取实际空间人体表面信息。然后再将点集Q进行均匀网格简化得到点集Q'。
步骤S103,粗匹配。首先以特定边长的正方体将点集P'进行网格化(参见图4所示),保存这些网格的中心点C={ci|i=1,…,n},作为下一步粗匹配中点集P'的初始位置中心点。根据这n种初始位置中心点,设置点集P'的n种匹配初始位置与点集Q'进行八方位匹配(参见图5所示),从而求得候选的最佳变换矩阵的集合。其中,将点集P'作为目标点集和点集Q'作为原始点集进行分网格初始化粗匹配求取候选初始位置对应的变换矩阵的集合,包括以下步骤:
步骤S1031,将点集P'用边长为50mm的正方体网格进行划分,得到n个正方体网格,网格中心点保存为点集C={ci|i=1,…,n},求取点集Q'的中心点cq,对点集C中每一个点ci,分别设ci为点集C的中心点,将点集Q'和点集C进行八方位匹配,当ci为点集C的中心点时,点集C和点集Q'的最佳变换矩阵为Mi,将矩阵Mi保存到集合M={Mi|i=1,…,k};
步骤S1032,对于步骤S1031中提到的设置两点集的八方位匹配,设A为原始点集,B为目标点集,其方式是:平移点集A和点集B,使平移后的点集A'和点集B'的中心点都与坐标系原点重合;旋转点集A'和点集B',使旋转后的点集A1”和点集B”的拟平面法向量与坐标系z轴平行,点集B”为目标点集B的匹配初始位置,点集A1”为原始点集A的第一种匹配初始位置;分别将点集A1”绕z轴旋转90°、180°和270°得到点集A2”、A3”、A4”作为第二、三、四种初始位置;将点集A1”绕y轴旋转180°得到点集A5”作为第五种匹配初始位置;分别将点集A5”绕z轴旋转90°、180°和270°得到点集A6”、A7”、A8”作为第六、七、八种初始位置;分别将点集A1”~A8”与点集B”进行ICP匹配,得到八组变换矩阵和对应的误差,比较误差值,最小的误差对应的变换矩阵则为该八方位匹配步骤的最佳变换矩阵。
步骤S104,精匹配。分别用步骤S1031中求出来的变换矩阵集合M={Mi|i=1,…,k}中的每一个矩阵Mi对点集P中的坐标进行变换得到点集Pi,将变换后的点集Pi分别与点集Q进行匹配,求得k个匹配误差e={ei|i=1,…,k}和变换矩阵M1={M1i|i=1,…,k},找到最小误差值对应的序号s,0<s<k+1,则精匹配结果求得的变换矩阵为M1s。
步骤S105,求取总的变换矩阵。根据S103和S104求得的最佳匹配结果对应的粗匹配步骤和精匹配步骤的变换矩阵,可以求得最终总的最佳变换矩阵为Mat=M1s*Ms。将激光点集P用最佳变换矩阵Mat进行变换,可以得到点集P对应到图像空间中的位置。如图6所示,为最终匹配结果,其中灰度较暗,范围较大的是头部在图像空间中形成的点集,在该点集表面的灰度较亮,范围较小的是激光点集。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (3)
1.应用于无标记点光学手术导航系统的病人表面注册方法,其特征在于,包括以下步骤:
1)导入病人CT图像序列,对CT图像序列提取外表面轮廓,记为点集P,将点集P简化为点集P';
2)用红外激光笔扫描病人表面,光学导航仪实时记录病人表面反射的红外光点的三维坐标,将这些坐标保存为原始激光点集Q,将点集Q简化得到点集Q';
3)分别将点集P'作为目标点集,点集Q'作为原始点集进行分网格初始化粗匹配,得到的候选粗匹配变换矩阵的集合M={Mi|i=1,…,k},其中k为网格的个数,包括以下步骤:
3.1)将点集P'用边长为50mm的正方体网格进行划分,得到n个正方体网格,网格中心点保存为点集C={ci|i=1,…,n},求取点集Q'的中心点cq,对点集C中每一个点ci,分别设ci为点集C的中心点,将点集Q'和点集C进行八方位匹配,当ci为点集C的中心点时,点集C和点集Q'的最佳变换矩阵为Mi,将矩阵Mi保存到集合M={Mi|i=1,…,k};
3.2)对于步骤3.1)中提到的设置两点集的八方位匹配,设A为原始点集,B为目标点集,其方式是:平移点集A和点集B,使平移后的点集A'和点集B'的中心点都与坐标系原点重合;旋转点集A'和点集B',使旋转后的点集A1”和点集B”的拟平面法向量与坐标系z轴平行,点集B”为目标点集B的匹配初始位置,点集A1”为原始点集A的第一种匹配初始位置;分别将点集A1”绕z轴旋转90°、180°和270°得到点集A2”、A3”、A4”作为第二、三、四种初始位置;将点集A1”绕y轴旋转180°得到点集A5”作为第五种匹配初始位置;分别将点集A5”绕z轴旋转90°、180°和270°得到点集A6”、A7”、A8”作为第六、七、八种初始位置;分别将点集A1”~A8”与点集B”进行ICP匹配,得到八组变换矩阵和对应的误差,比较误差值,最小的误差对应的变换矩阵则为该八方位匹配步骤的最佳变换矩阵;
4)结合步骤3)中粗匹配求得的k种粗匹配变换矩阵结果,将点集P作为匹配的目标点集,将点集Q作为匹配的原始点集进行精匹配,求得精匹配中的最佳匹配矩阵M1s和对应的粗匹配变换矩阵Ms,0<s<k+1,具体如下:
分别用步骤3.1)中求出来的变换矩阵集合M={Mi|i=1,…,k}中的每一个矩阵Mi对点集P中的坐标进行变换得到点集Pi,将变换后的点集Pi分别与点集Q进行匹配,求得k个匹配误差e={ei|i=1,…,k}和变换矩阵M1={M1i|i=1,…,k},找到最小误差值对应的序号s,0<s<k+1,则精匹配步骤求得的变换矩阵为M1s;
5)根据步骤4)中的结果得到表面注册的实际空间和图像空间的变换矩阵为Mat=M1s*Ms。
2.根据权利要求1所述的应用于无标记点光学手术导航系统的病人表面注册方法,其特征在于:在步骤1)、步骤2)中,点集简化方法是:设原始点集为S,简化后的点集为S',简化边长为l,将点集S用边长为l的正方体进行划分,对于每个正方体,求取其中心点坐标,若在该正方体内部有属于点集S的点,则找到这些点中距离正方体中心点坐标最近的点并保存得到简化后的点集S'。
3.根据权利要求1所述的应用于无标记点光学手术导航系统的病人表面注册方法,其特征在于:在步骤2)中,使用的红外激光笔包含两束激光,一束是波长为650nm的可见光,用于指示红外光的位置;另一束是波长为980nm的不可见红外光,将其照射到人体表面,红外光学定位仪获取实际空间人体表面信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306492.XA CN110123451B (zh) | 2019-04-17 | 2019-04-17 | 应用于无标记点光学手术导航系统的病人表面注册方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306492.XA CN110123451B (zh) | 2019-04-17 | 2019-04-17 | 应用于无标记点光学手术导航系统的病人表面注册方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110123451A CN110123451A (zh) | 2019-08-16 |
CN110123451B true CN110123451B (zh) | 2020-07-28 |
Family
ID=67570039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910306492.XA Expired - Fee Related CN110123451B (zh) | 2019-04-17 | 2019-04-17 | 应用于无标记点光学手术导航系统的病人表面注册方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110123451B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112837391B (zh) * | 2021-03-04 | 2022-02-18 | 北京柏惠维康科技有限公司 | 坐标转换关系获取方法、装置、电子设备及存储介质 |
CN113274130A (zh) * | 2021-05-14 | 2021-08-20 | 上海大学 | 用于光学手术导航系统的无标记手术注册方法 |
CN114081626B (zh) * | 2021-10-21 | 2023-06-20 | 华南理工大学 | 基于光学定位仪的兼容少量标记点丢失的工具识别方法 |
CN117372661B (zh) * | 2023-12-07 | 2024-03-12 | 华科精准(北京)医疗科技有限公司 | 一种手术导航系统、手术机器人系统及注册方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
-
2019
- 2019-04-17 CN CN201910306492.XA patent/CN110123451B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110123451A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110123451B (zh) | 应用于无标记点光学手术导航系统的病人表面注册方法 | |
CN109785374B (zh) | 一种牙科增强现实手术导航的自动实时无标记图像配准方法 | |
CN107951538B (zh) | 联合骨切除与钛板定位的3d打印腓骨重建颌骨手术工具的制造方法 | |
JP4821940B2 (ja) | 画像照合装置及びこれを用いた患者位置決め装置 | |
CN109615642B (zh) | 一种放疗计划中的危及器官自动勾画方法及装置 | |
CN107392995B (zh) | 机械轴导航系统中的人体下肢配准系统 | |
US8463004B2 (en) | Determining shaft and femur neck axes and three-dimensional reconstruction | |
CN114119549B (zh) | 一种多模态医学图像三维点云配准优化方法 | |
CN111161241B (zh) | 一种肝脏图像识别方法、电子设备及存储介质 | |
CN109493943B (zh) | 一种结合光学手术导航的三维可视化头皮开颅定位方法 | |
US20210374452A1 (en) | Method and device for image processing, and elecrtonic equipment | |
CN109754396A (zh) | 图像的配准方法、装置、计算机设备和存储介质 | |
CN109994188A (zh) | 基于ndi的神经外科导航配准测试方法和系统 | |
CN112907642A (zh) | 术前ct或核磁图像与其术中对应病灶精确配准重合的方法、系统、存储介质及设备 | |
CN104867104A (zh) | 基于xct图像非刚度配准的目标鼠解剖结构图谱获取方法 | |
CN106327479A (zh) | 血管造影中介下先心病术中血管辨识的装置及方法 | |
CN109965979A (zh) | 一种无需标志点的稳健的神经导航自动注册方法 | |
CN114821031A (zh) | 基于c臂机的术中图像匹配方法、装置和系统 | |
CN115578320A (zh) | 一种骨科手术机器人全自动空间注册方法及系统 | |
CN113633377B (zh) | 一种面向胫骨高位截骨手术的胫骨优化配准系统与方法 | |
CN115526929A (zh) | 基于图像的配准方法和装置 | |
CN114283188A (zh) | 一种针对骨科手术机器人的定位配准方法 | |
CN113274130A (zh) | 用于光学手术导航系统的无标记手术注册方法 | |
JP2007534416A5 (zh) | ||
CN116883471B (zh) | 面向胸腹部经皮穿刺的线结构光无接触点云配准方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200728 |
|
CF01 | Termination of patent right due to non-payment of annual fee |