CN103908346B - 一种高精度自动神经导航空间配准方法 - Google Patents
一种高精度自动神经导航空间配准方法 Download PDFInfo
- Publication number
- CN103908346B CN103908346B CN201210595045.9A CN201210595045A CN103908346B CN 103908346 B CN103908346 B CN 103908346B CN 201210595045 A CN201210595045 A CN 201210595045A CN 103908346 B CN103908346 B CN 103908346B
- Authority
- CN
- China
- Prior art keywords
- points
- point set
- registration
- point
- matched
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 239000000203 mixture Substances 0.000 claims description 19
- 230000009466 transformation Effects 0.000 claims description 9
- 210000005036 nerve Anatomy 0.000 claims description 6
- 238000013138 pruning Methods 0.000 claims description 6
- 230000001537 neural effect Effects 0.000 description 7
- 230000007547 defect Effects 0.000 description 6
- 238000013507 mapping Methods 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Image Processing (AREA)
Abstract
本发明属医学图像处理及应用领域,涉及一种高精度自动神经导航空间配准方法。本发明方法中采用基于高斯混合模型的全局优化点集配准方法,将基于激光扫描仪获得的点云与CT获得的点云得到两个点云之间的坐标变换实现初步配准;然后采用无序点自动配准的方法,将初步配准的结果与神经导航中通过CT重建出来的病人空间再次进行配准。使用结果表明,本发明所述的方法实施简单,精度可靠,便于临床应用,可集成在现有导航***中,从而大幅度提高导航***精度。
Description
技术领域
本发明属医学图像处理及应用领域,涉及一种神经导航空间配准方法;尤其是一种高精度自动神经导航空间配准方法。该方法能使手术导航在临床应用中更精确、实用和方便。
背景技术
神经导航***作为一种医学辅助定位设备,能够为医生提供实时的术中信息,有效弥补了传统神经外科手术方式的不足,提高了手术质量。从临床应用情况可以看出,神经导航***的核心功能是使用图像来对手术器械进行定位和引导,因此,手术真实场景与图像引导空间的定位和配准精度将直接影响神经导航***的性能。
目前,神经导航***中使用的空间配准方法主要基于点匹配和面匹配;所述点匹配方法使用较多,但其存在一些固有缺陷,如:需要专门为导航进行一次影像扫描;粘贴在头皮表面的人工标记物容易发生移位,从而引起较大的空间配准误差;标记点识别过程消耗时间较长。为了解决上述点匹配方法存在的问题,基于面匹配的空间配准方法得到关注。所述面配准使用病人固有特征进行空间配准,可以直接使用已有的图像进行导航,克服了点配准存在的缺陷,同时,因为使用了更多的信息进行空间配准,可以达到更高的配准精度。使用基于面匹配的空间配准方法时,医生首先在图像空间提取一个代表病人面部轮廓的表面点云,然后在病人空间使用激光扫描仪获取相同表面的另一点云,然后基于两个点云之间的相关性得到两个空间之间的坐标变换。
所述的点云之间的相关性,总体上可分为单点与单点的一一对应,单点与多点间的相关性,多点与多点间相关性即模糊相关性三个方面;其中,基于单点与单点一一对应进行点云配准的主流方法是迭代最近点(ICP:IterativeClosestPoint)算法及其各种改进算法;该类方法直观、简单且具有较高的匹配精度,但存在对初始位置敏感等缺陷。所述单点与多点间的相关性,也就是点与点间的对应关系通过一定的权值来表示,即模型点集中的每一个点与场景点集中所有点的加权组合相对应。然而,在临床应用中,单点间的一一对应和点与多点之间的相关性往往较难确定;为了解决上述问题,无需建立明确的点之间相关性,即基于模糊相关性的点集配准方法被提出。该类方法主要基于统计和测度理论,使用概率密度分布描述点集空间;然而,点集配准各类方法中存在的一个比较共性的问题——收敛到局部最优解而非全局最优解依然存在。
因此,目前需要一种自动神经导航空间配准方法;该方法能使手术导航在临床应用中更精确、实用和方便。
与本发明有关的参考文献有:
1G.Eggers,J.MuhlingandR.Marmulla,Image-to-patientregistrationtechniquesinheadsurgery,Int.J.Oral.Max.Surg.,2006,35:1081-1095
2P.J.BeslandN.D.McKay,Amethodforregistrationof3-Dshapes,IEEETransactionsonPatternAnalysisandMachineIntelligence,1992,14(2):239-256。
发明内容
本发明的目的在于提供一种自动神经导航空间配准方法;尤其是一种高精度自动神经导航空间配准方法。该方法能使手术导航在临床应用中更精确、实用和方便。
本发明采用基于高斯混合模型的全局优化点集配准方法,并结合无序点自动匹配方法进一步提高配准精度,即:通过空间定位仪使手术真实场景与图像引导空间处于同一空间,并使用激光扫描仪获取病人头部表面的点云;然后基于激光扫描仪获得的点云与CT获得的点云之间相关性得到两个点云之间的坐标变换实现初步配准;在神经导航过程中,使用无序点自动配准方法将初步配准的结果进行再次匹配,提高配准精度。
本发明所述的方法中,首先采用基于高斯混合模型的全局优化点集配准方法,将基于激光扫描仪获得的点云与CT获得的点云得到两个点云之间的坐标变换实现初步配准,然后采用无序点自动配准的方法,将初步配准的结果与神经导航中通过CT重建出来的病人空间再次进行配准;
其中,
所述采用基于高斯混合模型的全局优化点集配准方法,将基于激光扫描仪获得的点云与CT获得的点云得到两个点云之间的坐标变换实现初步配准,是用激光扫描仪获得的点云和CT获得的点云采用基于高斯混合模型的全局优化点集配准方法实现配准;
所述采用无序点自动配准的方法,将初步配准的结果与神经导航中通过CT重建出来的病人空间再次进行配准,是将初步配准结果与神经导航中通过CT重建出来的病人空间采用无序点自动配准的方法再次进行配准。
具体而言,本发明的自动神经导航空间配准方法,其特征在于,包括以下步骤:
(1)基于高斯混合模型的全局优化点集配准方法
传统的点集配准方法通常将点集中各点的位置视作独立特征,点集的整体特征往往被忽略,而该方法将进行配准的点集看作一个整体,将点集内点与点在整体空间和子空间的相关性作为点集的特征,实施配准,从而最终实现了基于全局优化的点集配准;本发明提出的方法引入了点与点间的相关性,不仅解决了局部收敛的缺陷,同时还提高了配准的精度;
通常,高斯混合模型的概率密度函数可表示为:
其中,
根据高斯混合模型的定义,本发明定义了点集所表示的数据空间中点与点之间的相关性,可表示为:
从而,对于点集P所表示的整个数据空间中所有点相关性总和,即该数据空间的全局特征可表示为:
对于3维点集P表示的3维数据空间经过投影变换,可形成3个1维数据空间;每个1维数据空间可使用混合高斯模型进行表示,其表示形式与(4)一致;
假设,需要进行配准的两个点集,分别是同一数据空间不同的采样结果,因此它们代表了同一的数据空间,其相应的投影变换形成的子空间也是一致的,对应子空间的点与点之间相关性总和也可认为是接近相等的;基于上述思想,刚体点集配准问题可以转变为最小化式(6)所表示的目标函数,Ml和Sl分别表示两个点集代表的对应子空间。
dRMSl=(gmm_correlation(Ml)-gmm_correlation(Sl))2(6)
在进行点集配准的过程中,假设表示整个数据空间的高斯混合模型各个混合分量是均匀分布的,即式(4)中的同时,各个混合分量的协方差矩阵是一致的,即式(4)中的∑i=∑j;
(2)无序点自动匹配方法
将待匹配的点集按照点与点之间的距离关系映射到新的点集空间,然后根据距离阈值与目标点集通过剪枝的方法获得最终的匹配关系;该方法大大提高了导航配准精度和效率;
其包括步骤为:
①获取手术真实场景中标记点做为目标点集,获取步骤(1)中初步配准结果中的特征点做为待匹配点集;
②从待匹配点集中取出四个点,对其中每一个点获得与其他三个点之间的距离关系并按照一定的顺序进行排列,形成一个新的点的空间描述;
③对四个点的新的空间描述做如下处理:每两个点对应元素之间做差值并取绝对值然后求和;取所有和值中的最小值;
④对待匹配点集中所有的四个点的组合重复上述步骤②和步骤③中的操作,获得步骤③中所有最小值中的最大值对应的四个点的组合;
⑤对目标点集做四个点的排列,每一组排列分别与上述步骤④中获得的四个点的组合做剪枝,逐步排除不符合阈值判定的目标点集的排列;
⑥最后获得与待匹配点集相对应的目标点集的一组排列。
本发明的自动神经导航空间配准方法与现有技术相比,具有以下优点:
(1)解决了基于点匹配的的空间配准弊端;
(2)本方法为一种全自动的图像病人空间配准方法,减少了医生的工作量;
(3)本方法为一种全局优化的点集配准方法,降低了初始位置对点集配准方法的影响;
(4)提高了配准效率和精度,减小了配准误差。
本发明所述的方法实施简单,精度可靠,便于临床应用,可集成在现有导航***中,从而大幅度提高导航***精度。
具体实施方式
实施例1
本发明的自动神经导航空间配准方法,包括以下步骤:
(1)基于高斯混合模型的全局优化点集配准方法
将进行配准的点集看作一个整体,将点集内点与点在整体空间和子空间的相关性作为点集的特征,实施配准,从而最终实现了基于全局优化的点集配准;本发明提出的方法引入了点与点间的相关性,不仅解决了局部收敛的缺陷,同时还提高了配准的精度;
通常,高斯混合模型的概率密度函数可表示为:
其中,
根据高斯混合模型的定义,本发明定义了点集所表示的数据空间中点与点之间的相关性,可表示为:
从而,对于点集P所表示的整个数据空间中所有点相关性总和,即该数据空间的全局特征可表示为:
对于3维点集P表示的3维数据空间经过投影变换,可形成3个1维数据空间;每个1维数据空间可使用混合高斯模型进行表示,其表示形式与(4)一致;
假设,需要进行配准的两个点集,分别是同一数据空间不同的采样结果,因此它们代表了同一的数据空间,其相应的投影变换形成的子空间也是一致的,对应子空间的点与点之间相关性总和也可认为是接近相等的;基于上述思想,刚体点集配准问题可以转变为最小化式(6)所表示的目标函数,Ml和Sl分别表示两个点集代表的对应子空间。
dRMSl=(gmm_correlation(Ml)-gmm_correlation(Sl))2(6)
在进行点集配准的过程中,假设表示整个数据空间的高斯混合模型各个混合分量是均匀分布的,即式(4)中的同时,各个混合分量的协方差矩阵是一致的,即式(4)中的∑i=∑j;
(2)无序点自动匹配方法
将待匹配的点集按照点与点之间的距离关系映射到新的点集空间,然后根据距离阈值与目标点集通过剪枝的方法获得最终的匹配关系;该方法大大提高了导航配准精度和效率;
其包括步骤为:
①获取手术真实场景中标记点做为目标点集,获取步骤(1)中初步配准结果中的特征点做为待匹配点集;
②从待匹配点集中取出四个点,对其中每一个点获得与其他三个点之间的距离关系并按照一定的顺序进行排列,形成一个新的点的空间描述;
③对四个点的新的空间描述做如下处理:每两个点对应元素之间做差值并取绝对值然后求和;取所有和值中的最小值;
④对待匹配点集中所有的四个点的组合重复上述步骤②和步骤③中的操作,获得步骤③中所有最小值中的最大值对应的四个点的组合;
⑤对目标点集做四个点的排列,每一组排列分别与上述步骤④中获得的四个点的组合做剪枝,逐步排除不符合阈值判定的目标点集的排列;
⑥最后获得与待匹配点集相对应的目标点集的一组排列。
上述实施例的结果表明,本发明解决了基于点匹配的的空间配准弊端;所述方法为一种全自动的图像病人空间配准方法,减少了医生的工作量;该方法还是一种全局优化的点集配准方法,降低了初始位置对点集配准方法的影响;本发明提高了配准效率和精度,减小了配准误差;此外,所述方法实施简单,精度可靠,便于临床应用,可集成在现有导航***中,从而大幅度提高导航***精度。
Claims (2)
1.一种高精度自动神经导航空间配准方法,其特征在于,包括步骤:
(1)采用基于高斯混合模型的全局优化点集配准方法,将基于激光扫描仪获得的点云与CT获得的点云得到两个点云之间的坐标变换实现初步配准;
(2)采用无序点自动配准的方法,将初步配准的结果与神经导航中通过CT重建出来的病人空间再次进行配准。
2.根据权利要求1所述的方法,其特征在于,所述步骤(2)中无序点自动匹配方法为,将待匹配的点集按照点与点之间的距离关系映射到新的点集空间,然后根据距离阈值与目标点集通过剪枝的方法获得最终的匹配关系,其包括步骤为:
①获取手术真实场景中标记点做为目标点集,获取所述步骤(1)中初步配准结果中的特征点做为待匹配点集;
②从待匹配点集中取出四个点,对其中每一个点获得与其他三个点之间的距离关系并按照一定的顺序进行排列,形成一个新的点的空间描述;
③对四个点的新的空间描述做如下处理:每两个点对应元素之间做差值并取绝对值然后求和;取所有和值中的最小值;
④对待匹配点集中所有的四个点的组合重复上述步骤②和步骤③中的操作,获得步骤③中所有最小值中的最大值对应的四个点的组合;
⑤对目标点集做四个点的排列,每一组排列分别与上述步骤④中获得的四个点的组合做剪枝,逐步排除不符合阈值判定的目标点集的排列;
⑥最后获得与待匹配点集相对应的目标点集的一组排列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210595045.9A CN103908346B (zh) | 2012-12-31 | 2012-12-31 | 一种高精度自动神经导航空间配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210595045.9A CN103908346B (zh) | 2012-12-31 | 2012-12-31 | 一种高精度自动神经导航空间配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103908346A CN103908346A (zh) | 2014-07-09 |
CN103908346B true CN103908346B (zh) | 2016-04-20 |
Family
ID=51034556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210595045.9A Active CN103908346B (zh) | 2012-12-31 | 2012-12-31 | 一种高精度自动神经导航空间配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103908346B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105701788A (zh) * | 2014-11-24 | 2016-06-22 | 复旦大学 | 一种基于全局优化的神经导航自动空间配准方法 |
US10709352B2 (en) * | 2015-10-27 | 2020-07-14 | Covidien Lp | Method of using lung airway carina locations to improve ENB registration |
CN105852971A (zh) * | 2016-05-04 | 2016-08-17 | 苏州点合医疗科技有限公司 | 基于骨骼三维中点云的配准导航方法 |
KR102604919B1 (ko) * | 2016-07-14 | 2023-11-23 | 인튜어티브 서지컬 오퍼레이션즈 인코포레이티드 | 컴퓨터 보조 원격조작 시스템에서의 2차 기구 제어 |
CN106344154B (zh) * | 2016-09-14 | 2018-11-09 | 大连理工大学 | 一种基于最大相关熵的手术器械尖端点的标定方法 |
CN109965979A (zh) * | 2017-12-27 | 2019-07-05 | 上海复旦数字医疗科技股份有限公司 | 一种无需标志点的稳健的神经导航自动注册方法 |
CN108324369B (zh) * | 2018-02-01 | 2019-11-22 | 艾瑞迈迪医疗科技(北京)有限公司 | 基于面的术中配准方法及神经导航设备 |
CN108836479B (zh) * | 2018-05-16 | 2020-01-24 | 山东大学 | 一种医学图像配准方法和手术导航*** |
CN109745120B (zh) * | 2018-12-24 | 2020-07-31 | 罗雄彪 | 一种影像配准转换参数优化方法及*** |
CN111260702B (zh) * | 2020-02-13 | 2022-05-17 | 北京航空航天大学 | 激光三维点云与ct三维点云配准方法 |
CN114241021A (zh) * | 2021-12-15 | 2022-03-25 | 商丘市第一人民医院 | 一种基于空间配准的脊柱外科导航方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975951A (zh) * | 2010-06-09 | 2011-02-16 | 北京理工大学 | 一种融合距离和图像信息的野外环境障碍检测方法 |
US8295547B1 (en) * | 2010-05-26 | 2012-10-23 | Exelis, Inc | Model-based feature tracking in 3-D and 2-D imagery |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6791542B2 (en) * | 2002-06-17 | 2004-09-14 | Mitsubishi Electric Research Laboratories, Inc. | Modeling 3D objects with opacity hulls |
US8884618B2 (en) * | 2009-02-25 | 2014-11-11 | Zimmer, Inc. | Method of generating a patient-specific bone shell |
-
2012
- 2012-12-31 CN CN201210595045.9A patent/CN103908346B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8295547B1 (en) * | 2010-05-26 | 2012-10-23 | Exelis, Inc | Model-based feature tracking in 3-D and 2-D imagery |
CN101975951A (zh) * | 2010-06-09 | 2011-02-16 | 北京理工大学 | 一种融合距离和图像信息的野外环境障碍检测方法 |
Non-Patent Citations (2)
Title |
---|
G.Eggers,J.MUhling,R.Marmulla.Image-to-patient registration techniques in headsurgery.《Int. J.Oral.Max.Surg.》.2006,第35卷1081*1095. * |
Paul J. Besl, Neil D. McKay.A method for registration of 3-D shapes.《IEEE Transactions on Pattern Analysis and Machine Intelligence》.1992,第14卷(第2期),239-256. * |
Also Published As
Publication number | Publication date |
---|---|
CN103908346A (zh) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103908346B (zh) | 一种高精度自动神经导航空间配准方法 | |
CN109785374B (zh) | 一种牙科增强现实手术导航的自动实时无标记图像配准方法 | |
CN103325143B (zh) | 基于模型匹配的标记点自动注册方法 | |
CN106934821B (zh) | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 | |
CN110533664B (zh) | 一种基于大样本数据驱动的颅神经自动分割方法 | |
Deng et al. | A novel skull registration based on global and local deformations for craniofacial reconstruction | |
CN103860268B (zh) | 一种标记点配准方法、装置及外科手术导航*** | |
CN109166087A (zh) | 医学图像的风格转换方法、装置、医学设备、影像***及存储介质 | |
CN104867104B (zh) | 基于xct图像非刚度配准的目标鼠解剖结构图谱获取方法 | |
CN105310776B (zh) | 一种基于子块的软组织表面变形追踪方法 | |
CN105205810B (zh) | 基于距离场融合的mr图像海马体分割方法 | |
CN104881873B (zh) | 一种用于复杂纤维束精确重构的多级调整混合加权稀疏成像方法 | |
CN109965979A (zh) | 一种无需标志点的稳健的神经导航自动注册方法 | |
CN111640143A (zh) | 一种基于PointNet的神经导航快速面配准方法及*** | |
CN107507189A (zh) | 基于随机森林与统计模型的小鼠ct图像肾脏分割方法 | |
CN115578320A (zh) | 一种骨科手术机器人全自动空间注册方法及*** | |
CN116650115A (zh) | 一种基于uwb标记点的骨科手术导航注册方法 | |
CN107330948B (zh) | 一种基于流行学习算法的fMRI数据二维可视化方法 | |
CN113274130A (zh) | 用于光学手术导航***的无标记手术注册方法 | |
CN115116113B (zh) | 一种光学导航方法 | |
CN113920179B (zh) | 一种基于标记点的多视觉2d-3d图像非刚性配准方法及*** | |
CN115690045A (zh) | 基于曲面断层片的牙周炎治疗前后骨增量定量评估方法 | |
CN113648060B (zh) | 超声引导的软组织形变跟踪方法、装置、存储介质及*** | |
LU500422B1 (en) | Highly accurate automatic neural navigation spatial registration method | |
CN116168097A (zh) | 构建cbct勾画模型和勾画cbct图像的方法、装置、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Song Zhijian Inventor after: Chen Xinrong Inventor after: An Yong Inventor after: Wang Manning Inventor after: Yao Demin Inventor before: Song Zhijian Inventor before: Chen Xinrong Inventor before: An Yong Inventor before: Wang Manning Inventor before: Yao Demin |
|
COR | Change of bibliographic data | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |