CN104008269B - 一种基于人工标记物的手术导航系统空间自动注册方法 - Google Patents
一种基于人工标记物的手术导航系统空间自动注册方法 Download PDFInfo
- Publication number
- CN104008269B CN104008269B CN201410134319.3A CN201410134319A CN104008269B CN 104008269 B CN104008269 B CN 104008269B CN 201410134319 A CN201410134319 A CN 201410134319A CN 104008269 B CN104008269 B CN 104008269B
- Authority
- CN
- China
- Prior art keywords
- space
- point
- sigma
- points
- class
- 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
Landscapes
- Image Analysis (AREA)
Abstract
一种基于人工标记物的手术导航系统空间自动注册方法,该方法的步骤为:(1)在图像引导空间,对术前采集的医学影像图像进行表面重建,获得人头模型;采用不同的尺度空间因子,基于顶点预测器的原理,对模型进行光顺,构建模型的多尺度空间;(2)在多尺度空间中,计算相邻尺度的高斯差分,遍历当前尺度和相邻尺度的一环邻域寻找极值点;(3)对极值点进行标记物形状相关聚类,定位出标记物中心点;(4)在病人空间,使用双目相机实时获取病人头部图片,基于角点检测提取出标记物中心点,基于双目立体视觉原理恢复标记物中心点三维坐标;(5)基于随机迭代一致性(Ransac)方法匹配两空间标记物,计算欧式变换矩阵,完成手术导航系统的空间注册。
Description
技术领域
本发明涉及计算机辅助手术领域,更具体地讲,涉及一种图像引导手术导航系统空间自动注册方法。
背景技术
把病人在空间中的实际位置和三维模型中的虚拟位置进行精确配准的过程,在手术导航中称为注册。空间注册的精度直接影响着整个导航系统的精度,是手术前准备阶段中最重要的一步。基于点匹配的刚体配准法是目前手术导航系统临床应用中使用最多的空间注册法,其基本思想是分别在病人空间和三维图像空间中选取对应点集,通过数学原理计算出点集之间的变换关系,再扩展到病人空间中任意点在三维图像空间上的对应位置。
目前基于点匹配的空间注册方法大都需要医生手动、顺序地在图像空间和图像引导空间选择人工标记物,这种方法获得的注册精度依赖于医生的经验,并且手动选取精度往往不稳定,通常需要多次选取然后求取平均值,增加了医生的工作量并延长了手术准备时间。
本发明提出了一种基于人工标记物的手术导航系统空间自动注册方法,该方法设计了十字结构的圆形人工标记物,并利用标记物的形状特征,构建模型的多尺度空间,检测位于标记物圆周上的特征点,从而定位出标记物中心,在病人空间,利用角点检测和双目立体视觉原理,自动识别相同的标记物中心点,进而完成注册。这种方法不需要人工干预,稳定高效,大大节约了术前准备时间,提高了手术效率。
发明内容
本发明的目的是:解决现有手术导航系统需要医生手动、顺序选择标记物的费时费力的问题,提供了一种基于人工标记物的手术导航系统空间自动注册方法,该方法可以有效地完成病人术前医学影像空间与实际术中空间的注册,满足手术导航系统的速度和精度要求。
为了实现上述目的,本发明的原理是:在图像引导空间,利用人工标记物和人头表面接触处不同的凹凸性,检测出曲率骤变的标记物圆周,根据几何特点计算标记物中心点。在病人空间,本发明设计了圆形十字标记物,以标记物中心点为圆点,以黑白两色填充,每90°变换颜色,这样十字标记物中心点为一稳定的角点。利于人工标记物的特殊设计,检测十字标记的中心角点。同时采用双目相机基于立体视觉的技术恢复这些标记物中心点的三维坐标,然后寻找两个空间的标记物对应关系,计算两空间的欧式变换矩阵,从而完成注册。
本发明所采取的技术方案是:首先,在图像引导空间,对术前采集的医学影像图像进行表面重建,获得人头模型;采用不同的尺度空间因子,基于顶点预测器的原理,对模型进行光顺,构建模型的多尺度空间;在多尺度空间中,计算相邻尺度的高斯差分,遍历当前尺度和相邻尺度的一环邻域寻找极值点;然后对极值点进行标记物形状相关聚类,定位出标记物中心点。在病人空间,使用双目相机实时获取病人头部图片,基于角点检测提取出标记物中心点,基于双目立体视觉原理恢复标记物中心点三维坐标。最后,基于随机迭代一致性(Ransac)方法匹配两空间标记物,计算欧式变换矩阵,完成空间注册。
本发明与现有技术相比的优点在于:充分利用了标记物的几何特征和灰度特征,并结合多尺度空间原理和立体视觉原理,自动定位出人工标记物,有效地提高了空间注册的精度和效率,减少了医生的工作量。
附图说明
图1为本发明基于人工标记物的手术导航系统空间自动注册方法的流程图;
图2为图像引导空间人工标记物定位的流程图。
图3为定位出标记物示意图
具体实施方法
下面结合附图对本发明作进一步说明。
本发明提出的基于人工标记物的手术导航系统空间自动注册方法的流程如图1所示,其步骤如下:
(1)如图2所示,在图像引导空间,采用不同的尺度空间因子,基于顶点预测器的原理,对模型进行光顺,构建模型的多尺度空间;在多尺度空间中,计算相邻尺度的高斯差分,遍历当前尺度和相邻尺度的一环邻域寻找极值点;对极值点进行标记物形状相关聚类,定位出标记物中心点(如图3所示):
(1.1)设原始的人头模型为M,p为M中的一个顶点。对于尺度空间因子σ,M经过以σ为参数的高斯函数光顺后的结果为Mσ,p在Mσ中对应的顶点是pσ,定义顶点p处的信号强度l(p,σ)为p在pσ法向方向nσ(p)上的位移:
l(p,σ)=<nσ(p),p-pσ>
(1.2)在以σ为因子的尺度空间中,估算p点对应顶点pσ的法向nσ(p)。(是数学符号,表示对于属于M的任意一点p),设N(p)是p的邻域,q是邻域内的一点,即nq是q点的法向量。定义两点空间距离权值为||q-p||,同时考虑到网格M可能是非均匀采样的曲面,引入q的Voronoi面积作为顶点的面积权值Aq,利用如下公式估算出nσ(p):
其中Gσ(x)是高斯函数:
(1.3)计算p点的信号强度l(p,σ)。已知nσ(p),利用p在预测器∏q下的预测信号强度S(∏q(p))的加权平均来计算l(p,σ),如下式:
其中p点在预测器∏q下的预测信号强度S(∏q(p))为:
S(∏q(p))=nσ(p)·(q-p)
(1.4)选取不同的尺度空间因子:σ1,σ2,….,σn,其中σi+1=kσi,i=1,2,…,n-1,分别做上述处理,就可以得到三维模型M在不同尺度空间中的表示L(M,σi),i=1,2,…,n;
(1.5)计算出高斯差分函数的值D(M,σi):
D(M,σi)=L(M,σi+1)-L(M,σi)
(1.6)在p点所在尺度的1环邻域,以及相邻尺度的对应顶点和该顶点的1环领域内,寻找具有最大和最小DOG值的顶点作特征点。
(1.7)从检测到的特征点中使用Kmeans++的方法选取k个(k大于实际粘贴的标记点数目)特征点做初始聚类中心,然后扫描所有特征点并计算它到每一个类中心的距离,并将该点加入到距离最近的类中。
(1.8)迭代循环所有的类,如果发现某两个类之间的距离(类间距离定义为两个类中心之间的距离)小于一个阈值d(例如20mm,约为标记点直径),则将这两个类合并,直至所有类间距离都大于这个阈值为止。遍历所得的聚类,如果类中包含的特征点数少于阈值t,则删除该类,否则保留该类。
(1.9)计算类内特征点到类中心的距离,如果距离在某一范围内(如7~13mm,约为标记点直径±3mm),则保留该特征点,否则删除。处理完所有类后,重复(1.7)、(1.8)。
(1.10)由于对应于标记物的类中的特征点都分布在标记物圆周上,任意取三个特征点,根据三角形三边垂直平分线的交点即三角形外接圆圆心的原理,计算外接圆圆心,圆心到三点中任一点的距离即外接圆半径。然后计算类内特征点到圆心的距离,若与半径差值小于一定阈值(标记物尺寸是已知的),则认为该点合格,统计合格的特征点的比例,如果大于2/3,则认为该圆心合格。
(2)在病人空间,使用双目相机实时获取病人头部图片,基于角点检测提取出标记物中心点,基于双目立体视觉原理恢复标记物中心点三维坐标。
(3)基于随机迭代一致性(Ransac)方法匹配病人空间和图像引导空间标记物,计算欧式变换矩阵,完成空间注册:
(3.1)假设在图像引导空间和病人空间得到的标记物中心点集分别是P、Q,从P、Q中分别任意选三个点,组成3个点对。
(3.2)根据这3个点对计算出刚体变换矩阵T。
(3.3)用矩阵T对P进行变换,遍历变换后P中的点,在Q中找其最近点,如果小于预定义的阈值,则认为(Pi,Qj)是正确匹配的点对,否则(Pi,Qj)是错误匹配的点对。正确匹配的点对记为局内点(Inlier),错误的点对记为局外点(Outlier)。
(3.4)如果局内点的数目达到一定数量(与总点对的比例达到一定阈值),对所有的局内点采用基于单位四元数的参数估计方法重新计算M,否则迭代(3.1)-(3.3)步。
本发明中未详细阐述的部分属于本领域技术人员的公知技术。
Claims (1)
1.一种基于人工标记物的手术导航系统空间自动注册方法,其特征在于包括如下步骤:
(1)在图像引导空间,采用不同的尺度空间因子,基于顶点预测器的原理,对通过术前采集的医学影像图像进行表面重建的人头模型进行光顺,构建模型的多尺度空间;在多尺度空间中,计算相邻尺度的高斯差分,遍历当前尺度和相邻尺度的一环邻域寻找极值点;对极值点进行标记物形状相关聚类,定位出标记物中心点;
(2)在病人空间,使用双目相机实时获取病人头部图片,基于角点检测提取出标记物中心点,基于双目立体视觉原理恢复标记物中心点三维坐标;
(3)基于随机迭代一致性方法匹配两空间标记物,计算欧式变换矩阵,完成空间注册;
所述步骤(1)中标记物中心点定位方法,包括步骤如下:
(2.1)设原始的三维模型为M,p为M中的一个顶点;对于尺度空间因子σ,M经过以σ为参数的高斯函数光顺后的结果为Mσ,p在Mσ中对应的顶点是pσ,定义顶点p处的信号强度l(p,σ)为p在pσ法向方向nσ(p)上的位移:
l(p,σ)=<nσ(p),p-pσ>
(2.2)在以σ为因子的尺度空间中,估算p点对应顶点pσ的法向nσ(p);设N(p)是p的领域,q是邻域内的一点,即nq是q点的法向量;定义两点空间距离权值为||q-p||,同时考虑到网格M可能是非均匀采样的曲面,引入q的Voronoi面积作为顶点的面积权值Aq,利用如下公式估算出nσ(p):
其中Gσ(x)是高斯函数:
(2.3)计算p点的信号强度l(p,σ);已知nσ(p),利用p在预测器∏q下的预测信号强度S(∏q(p))的加权平均来计算l(p,σ),如下式:
其中p点在预测器∏q下的预测信号强度S(∏q(p))为:
S(∏q(p))=nσ(p)·(q-p)
(2.4)选取不同的尺度空间因子:σ1,σ2,…,σn,其中σi+1=Aσi,i=1,2,…,n-1,分别做上述处理,就得到三维模型M在不同尺度空间中的表示L(M,σi),i=1,2,…,n;n为选取的不同尺度因子的个数;
(2.5)计算出高斯差分函数的值D(M,σi):
D(M,σi)=L(M,σi+1)-L(M,σi)
(2.6)在p点所在尺度的1环领域,以及相邻尺度的对应顶点和该顶点的1环领域内,寻找具有最大和最小DOG值的顶点作特征点;
(2.7)从检测到的特征点中使用Kmeans++的方法选取k个特征点做初始聚类中心,k大于实际粘贴的标记点数目,然后扫描所有特征点并计算它到每一个类中心的距离,并将该点加入到距离最近的类中;
(2.8)迭代循环所有的类,如果发现某两个类之间的距离小于一个阈值d,则将这两个类合并,直至所有类间距离都大于这个阈值为止;所述类间距离定义为两个类中心之间的距离;遍历所得的聚类,如果类中包含的特征点数少于阈值t,则删除该类,否则保留该类;
(2.9)计算类内特征点到类中心的距离,如果距离在某一阈值范围内,则保留该特征点,否则删除;处理完所有类后,再重复一遍(2.7)、(2.8)步骤;
(2.10)由于对应于标记物的类中的特征点都分布在标记物圆周上,任意取三个特征点,根据三角形三边垂直平分线的交点即三角形外接圆圆心的原理,计算外接圆圆心,圆心到三点中任一点的距离即外接圆半径;然后计算类内特征点到圆心的距离,若与半径差值小于一定阈值,则认为该点合格,统计合格的特征点的比例,如果大于2/3,则认为该圆心合格。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410134319.3A CN104008269B (zh) | 2014-04-03 | 2014-04-03 | 一种基于人工标记物的手术导航系统空间自动注册方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410134319.3A CN104008269B (zh) | 2014-04-03 | 2014-04-03 | 一种基于人工标记物的手术导航系统空间自动注册方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104008269A CN104008269A (zh) | 2014-08-27 |
CN104008269B true CN104008269B (zh) | 2017-05-03 |
Family
ID=51368921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410134319.3A Expired - Fee Related CN104008269B (zh) | 2014-04-03 | 2014-04-03 | 一种基于人工标记物的手术导航系统空间自动注册方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104008269B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106560163B (zh) * | 2015-09-30 | 2019-11-29 | 合肥美亚光电技术股份有限公司 | 手术导航系统及手术导航系统的配准方法 |
CN106127753B (zh) * | 2016-06-20 | 2019-07-30 | 中国科学院深圳先进技术研究院 | 一种外科手术中ct影像体表人工标记自动提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103325143A (zh) * | 2013-06-13 | 2013-09-25 | 华南理工大学 | 基于模型匹配的标记点自动注册方法 |
CN103417295A (zh) * | 2012-05-17 | 2013-12-04 | 中国科学院深圳先进技术研究院 | 手术导航系统及手术导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9649160B2 (en) * | 2012-08-14 | 2017-05-16 | OrthAlign, Inc. | Hip replacement navigation system and method |
-
2014
- 2014-04-03 CN CN201410134319.3A patent/CN104008269B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103417295A (zh) * | 2012-05-17 | 2013-12-04 | 中国科学院深圳先进技术研究院 | 手术导航系统及手术导航方法 |
CN103325143A (zh) * | 2013-06-13 | 2013-09-25 | 华南理工大学 | 基于模型匹配的标记点自动注册方法 |
Non-Patent Citations (3)
Title |
---|
Distinctive image features from scale-invariant keypoints;David G.Lowe;《International Journal of Computer Vision》;20040105;第3-4节 * |
手术导航中双目立体视觉系统的设计及关键技术的研究;曹莉;《中国优秀硕士学位论文全文数据库 信息科技辑》;20111115(第11期);第2.1,5.3节 * |
手术导航中的三位注册方法;向华;《空军医学杂志》;20110930;第27卷(第3期);149-152页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104008269A (zh) | 2014-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108520519B (zh) | 一种图像处理方法、装置及计算机可读存储介质 | |
US11229419B2 (en) | Method for processing 3D image data and 3D ultrasonic imaging method and system | |
CN108898606B (zh) | 医学图像的自动分割方法、系统、设备及存储介质 | |
Kim et al. | Machine-learning-based automatic identification of fetal abdominal circumference from ultrasound images | |
CN105741251B (zh) | 一种肝脏cta序列图像的血管分割方法 | |
CN107274399B (zh) | 一种基于Hession矩阵和三维形状指数的肺结节分割方法 | |
CN108052942B (zh) | 一种飞机飞行姿态的视觉图像识别方法 | |
CN103325143B (zh) | 基于模型匹配的标记点自动注册方法 | |
CN106890031B (zh) | 一种标记物识别及标记点定位方法及手术导航系统 | |
EP2840553A1 (en) | 3d model data generation device, method and program | |
CN109087708B (zh) | 用于斑块分割的模型训练方法、装置、设备及存储介质 | |
CN103914823B (zh) | 基于稀疏表示的快速精确非线性配准立体医学影像的方法 | |
CN106127753B (zh) | 一种外科手术中ct影像体表人工标记自动提取方法 | |
CN105662474B (zh) | 胎儿头围超声图像的自动检测方法及检测系统 | |
CN104143193B (zh) | 一种超声图像中目标对象分割方法、装置及系统 | |
US9042614B2 (en) | Image registration method | |
CN107507212B (zh) | 数字脑可视化方法、装置、计算设备及存储介质 | |
CN107610121B (zh) | 一种肝脏统计形状模型初始位姿设定方法 | |
CN105232161A (zh) | 一种手术机器人标志点识别定位方法 | |
CN104008269B (zh) | 一种基于人工标记物的手术导航系统空间自动注册方法 | |
Ma et al. | Random walk based segmentation for the prostate on 3D transrectal ultrasound images | |
CN109410257A (zh) | 多信息3d医学图像高精度配准方法 | |
CN114155193A (zh) | 一种基于特征强化的血管分割方法及装置 | |
CN107464243B (zh) | 一种主动脉瓣的定位方法、装置和设备 | |
WO2014106747A1 (en) | Methods and apparatus for image processing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170503 Termination date: 20210403 |