CN109692041A - 一种针对覆盖软骨的骨表面配准方法 - Google Patents
一种针对覆盖软骨的骨表面配准方法 Download PDFInfo
- Publication number
- CN109692041A CN109692041A CN201910090929.0A CN201910090929A CN109692041A CN 109692041 A CN109692041 A CN 109692041A CN 201910090929 A CN201910090929 A CN 201910090929A CN 109692041 A CN109692041 A CN 109692041A
- Authority
- CN
- China
- Prior art keywords
- sample point
- bone
- bone surface
- point
- cartilage
- 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.)
- Granted
Links
- 210000000988 bone and bone Anatomy 0.000 title claims abstract description 69
- 210000000845 cartilage Anatomy 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 238000013507 mapping Methods 0.000 claims abstract description 5
- 239000000523 sample Substances 0.000 claims description 62
- 239000013074 reference sample Substances 0.000 claims description 6
- 238000002059 diagnostic imaging Methods 0.000 claims description 4
- 208000006735 Periostitis Diseases 0.000 description 3
- 210000003460 periosteum Anatomy 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 230000008569 process Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 210000000689 upper leg Anatomy 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
-
- 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
- A61B2034/2046—Tracking techniques
- A61B2034/2065—Tracking using image or pattern recognition
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Surgery (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Robotics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明的目的是提出一种针对覆盖软骨的骨表面配准方法,以提高配准精度,具体如下:分别获取待配准骨的骨表面样本点在影像坐标系下的第一骨表面点集S1和在世界坐标系下的第二骨表面点集S2;然后使用位置传感器获取待配准骨的辅助样本点在世界坐标系下的软骨表面点云S3;再配准第一骨表面点集S1与软骨表面点云S3,得到初始旋转矩阵R1和初始平移向量T1;并配准第二骨表面点集S2与软骨表面点云S3,得到修正旋转矩阵R2和修正平移向量T2;最后根据初始旋转矩阵R1、初始平移向量T1、修正旋转矩阵R2和修正平移向量T2,得到影像坐标系到世界坐标系的映射关系如下:旋转矩阵R=R2R1,平移向量T=R2T1+T2。
Description
技术领域
本发明属于医疗技术领域,具体涉及到一种骨科手术中,针对覆盖软骨的骨表面配准方法。
背景技术
配准技术是寻找最优空间变换的过程。手术导航系统广泛使用配准技术获得点集、线、面或者体之间的空间映射关系。配准主体则通常根据应用场景而变化,包括 但不限于点集之间、点集与曲面、高维特征之间、图像与图像等等。手术导航过程中, 通常需要借助位置传感器(光学、电磁定位系统)获取目标(器官、骨骼、肌肉组织 等)在世界坐标系中的坐标,借助医学成像(CT、MRI、超声)和重构技术获取目 标在影像坐标系下的坐标。
传统骨表面配准的方法如下:采集若干点在两个欧式空间的坐标,进行如下计算操作,即可获两个空间映射关系。计算过程表示如下:首先定义两个点集 Q=(q1,q2,...,qN)T和P=(p1,p2,...,pN)T,分别表示N个点在尺度相同的两个坐标系的对 应点集,再定义旋转矩阵R和平移向量T,表示两个欧式空间的变换关系,那么 Q=RP+T;
定义能量函数E,有其中,和分别表示为Q和P的样本平均;
通过奇异值分解(SVD),有
旋转矩阵R和平移向量T分别表示为R=VUT,
上述方法在多数情况下能够获得精准的配准结果,但是并不适用于股骨这类覆盖软骨的场景。一方面,在医学影像中,骨头边缘灰度特征明显,但软骨及骨膜厚度小 且灰度与软组织类似,因此,骨头边缘点易于获取,而骨膜与软骨边缘难以界定。借 助CT/MRI等医学影像往往可以提取骨表面坐标,对应地,借助位置传感器获取的则 是软骨/骨膜表面坐标。不难看出,两组点集所表示的并非同一目标,而是空间位置 接近的两个实体,这种差别会为精准配准引入误差;另一方面,从局部骨表面来看, 缺乏特征明显且能够精准定位的标记点,样本点的选择天然就带有误差,仅选取少量 所谓样本点对进行计算,是难以精确的。
发明内容
本发明的目的是提出一种针对覆盖软骨的骨表面配准方法,以提高配准精度。
本发明的针对覆盖软骨的骨表面配准方法包括如下步骤:
A、从待配准骨的医疗影像中提取骨表面的若干样本点在影像坐标系下的坐标,形成第一骨表面点集S1;
B、使用位置传感器探针获取待配准骨的骨表面的所述样本点在世界坐标系下的坐标,形成第二骨表面点集S2;
C、使用位置传感器获取待配准骨的辅助样本点在世界坐标系下的坐标,形成软骨表面点云S3;所述辅助样本点为待配准骨的软骨表面在所述样本点附近的预定数 量的点;
D、配准第一骨表面点集S1与软骨表面点云S3,得到初始旋转矩阵R1和初始平 移向量T1;
E、配准第二骨表面点集S2与软骨表面点云S3,得到修正旋转矩阵R2和修正平 移向量T2;
F、根据初始旋转矩阵R1、初始平移向量T1、修正旋转矩阵R2和修正平移向量 T2,得到影像坐标系到世界坐标系的映射关系如下:旋转矩阵R=R2R1,平移向量T= R2T1+T2。
上述A、B、C步骤中的获取相应点集(点云)的操作为现有技术中的常规操作, 此处不再赘述。
进一步地,由于B步骤中获取第二骨表面点集S2时,需要刺穿软骨或骨膜,出 于最大限度保护软组织的目的,建议点集S1、S2的规模较小,即所述A、B步骤中, 所述样本点的数量低于8个,且分布于感兴趣区域内。当然,有时为了提高准确性, 也可以适当扩大的点集S1、S2的规模,即所述A、B步骤中的样本点还包括位于所述 感兴趣区域外的参考样本点;所述C步骤中的辅助样本点还包括对应于参考样本点的 辅助参考样本点。
进一步地,为了将点集S1、S2中的稀疏样本都替换为密集点云,从而丰富软骨 表面的几何信息,提高配准的准确性,所述C步骤中,针对每个样本点,均在其周围 的软骨表面处设定若干个辅助样本点。由于C步骤无需伤害软组织,因此适当增多取 样点的数量,并不会对软骨造成伤害。
进一步地,为了提高点云S3与点集S1、S2的相关性,辅助样本点的取点原则如 下:所述C步骤中,针对同一个样本点的辅助样本点均匀分布于以该样本点为中心的 区域内。具体可以采用如下原则取点:所述C步骤中,针对同一个样本点的辅助样 本点位于以该样本点为圆心的若干个同心圆上。当然,还可以采用画三角形、画正方 形等中心对称图形(使对应的样本点位于该图形的中心位置处)的方法进行辅助样本 点的采集。
本发明将骨表面的样本点与软骨表面的辅助样本点结合起来,利用软骨表面的数量较多的辅助样本点点云与两种坐标系下的骨表面的数量较少的样本点点集分别进 行配准,相当于将点云粗配准和点集精配准相结合,既提高了配准精度,又能最大限 度地保护健康组织,具有配准抗噪能力强,鲁棒性能好的优点。
附图说明
图1是医学影像中的骨表面样本点的分布图。
图2是实体的骨表面样本点和辅助样本点的分布图。
具体实施方式
下面对照附图,通过对实施实例的描述,对本发明的具体实施方式如所涉及的各构件的形状、构造、各部分之间的相互位置及连接关系、各部分的作用及工作原理等 作进一步的详细说明。
实施例1:
本实施例提出了一种针对覆盖软骨的骨表面配准方法,以提高配准精度。
本实施例的针对覆盖软骨的骨表面配准方法包括如下步骤:
A、从待配准骨的医疗影像(CT模态)中提取骨表面的若干样本点在影像坐标系 下的坐标,形成第一骨表面点集S1;如图1所示,第一骨表面点集S1包括A、B、 C、D、E五个从股骨边缘上的样本点,采用的坐标系为影像坐标系;
B、使用位置传感器探针获取待配准骨的骨表面的所述样本点在世界坐标系下的坐标,形成第二骨表面点集S2;如图2所示,第二骨表面点集S2包括A′、B′、 C′、D′、E′五个从股骨边缘上的样本点,采用的坐标系为世界坐标系;其中A′、 B′、C′、D′、E′与A、B、C、D、E的位置分别对应;
C、使用位置传感器获取待配准骨的辅助样本点在世界坐标系下的坐标,形成软骨表面点云S3;所述辅助样本点为待配准骨的软骨表面在所述样本点附近的预定数 量的点;如图2所示,软骨表面点云S3包括A1′、A2′、A3′、A4′、B1′、B2′、B3′、B4′、C1′、C2′、C3′、C4′、D1′、D2′、D3′、D4′、E1′、E2′、E3′、E4′共20个辅助样本点,其中A1′、A2′、A3′、A4′均匀分布于样本点A的周围,具体来说,是位于以样本点A为圆心的若干个同心圆上;以同样的取点 原则来获取样本点B、C、D、E所对应的辅助样本点,此处不再赘述;
D、配准第一骨表面点集S1与软骨表面点云S3,得到初始旋转矩阵R1和初始平 移向量T1;
E、配准第二骨表面点集S2与软骨表面点云S3,得到修正旋转矩阵R2和修正平 移向量T2;
F、根据初始旋转矩阵R1、初始平移向量T1、修正旋转矩阵R2和修正平移向量 T2,得到影像坐标系到世界坐标系的映射关系如下:旋转矩阵R=R2R1,平移向量T= R2T1+T2。
本实施例在经典刚性配准方法基础上,增加了二次配准操作,提高了配准精度。
上面结合附图对本发明进行了示例性描述,显然本发明具体设计并不受上述方式的限制,只要采用了本发明的构思和技术方案进行的各种非实质性的改进,或未经改 进将本发明的构思和技术方案直接应用于其它场合的,均在本发明的保护范围之内。
Claims (6)
1.一种针对覆盖软骨的骨表面配准方法,其特征在于包括如下步骤:
A、从待配准骨的医疗影像中提取骨表面的若干样本点在影像坐标系下的坐标,形成第一骨表面点集S1;
B、使用位置传感器探针获取待配准骨的骨表面的所述样本点在世界坐标系下的坐标,形成第二骨表面点集S2;
C、使用位置传感器获取待配准骨的辅助样本点在世界坐标系下的坐标,形成软骨表面点云S3;所述辅助样本点为待配准骨的软骨表面在所述样本点附近的预定数量的点;
D、配准第一骨表面点集S1与软骨表面点云S3,得到初始旋转矩阵R1和初始平移向量T1;
E、配准第二骨表面点集S2与软骨表面点云S3,得到修正旋转矩阵R2和修正平移向量T2;
F、根据初始旋转矩阵R1、初始平移向量T1、修正旋转矩阵R2和修正平移向量T2,得到影像坐标系到世界坐标系的映射关系如下:旋转矩阵R=R2R1,平移向量T=R2T1+T2。
2.根据权利要求1所述的针对覆盖软骨的骨表面配准方法,其特征在于所述A、B步骤中,所述样本点的数量低于8个,且分布于感兴趣区域内。
3.根据权利要求2所述的针对覆盖软骨的骨表面配准方法,其特征在于所述A、B步骤中的样本点还包括位于所述感兴趣区域外的参考样本点;所述C步骤中的辅助样本点还包括对应于参考样本点的辅助参考样本点。
4.根据权利要求1所述的针对覆盖软骨的骨表面配准方法,其特征在于所述C步骤中,针对每个样本点,均在其周围的软骨表面处设定若干个辅助样本点。
5.根据权利要求4所述的针对覆盖软骨的骨表面配准方法,其特征在于所述C步骤中,针对同一个样本点的辅助样本点均匀分布于以该样本点为中心的区域内。
6.根据权利要求5所述的针对覆盖软骨的骨表面配准方法,其特征在于所述C步骤中,针对同一个样本点的辅助样本点位于以该样本点为圆心的若干个同心圆上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910090929.0A CN109692041B (zh) | 2019-01-30 | 2019-01-30 | 一种针对覆盖软骨的骨表面配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910090929.0A CN109692041B (zh) | 2019-01-30 | 2019-01-30 | 一种针对覆盖软骨的骨表面配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109692041A true CN109692041A (zh) | 2019-04-30 |
CN109692041B CN109692041B (zh) | 2021-09-21 |
Family
ID=66234651
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910090929.0A Active CN109692041B (zh) | 2019-01-30 | 2019-01-30 | 一种针对覆盖软骨的骨表面配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109692041B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111227932A (zh) * | 2020-02-19 | 2020-06-05 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
CN111297475A (zh) * | 2020-02-19 | 2020-06-19 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
CN112155732A (zh) * | 2020-09-29 | 2021-01-01 | 苏州微创畅行机器人有限公司 | 可读存储介质、骨建模配准系统及骨科手术系统 |
CN113081265A (zh) * | 2021-03-24 | 2021-07-09 | 重庆博仕康科技有限公司 | 手术导航空间配准方法、装置以及手术导航系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102933163A (zh) * | 2010-04-14 | 2013-02-13 | 史密夫和内修有限公司 | 用于基于患者的计算机辅助手术程序的系统和方法 |
WO2017204832A1 (en) * | 2016-05-27 | 2017-11-30 | Mako Surgical Corp. | Preoperative planning and associated intraoperative registration for a surgical system |
-
2019
- 2019-01-30 CN CN201910090929.0A patent/CN109692041B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102933163A (zh) * | 2010-04-14 | 2013-02-13 | 史密夫和内修有限公司 | 用于基于患者的计算机辅助手术程序的系统和方法 |
WO2017204832A1 (en) * | 2016-05-27 | 2017-11-30 | Mako Surgical Corp. | Preoperative planning and associated intraoperative registration for a surgical system |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111227932A (zh) * | 2020-02-19 | 2020-06-05 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
CN111297475A (zh) * | 2020-02-19 | 2020-06-19 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
CN111227932B (zh) * | 2020-02-19 | 2021-06-04 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
CN111297475B (zh) * | 2020-02-19 | 2021-06-18 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
WO2021164219A1 (zh) * | 2020-02-19 | 2021-08-26 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
WO2021164220A1 (zh) * | 2020-02-19 | 2021-08-26 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
US11426247B2 (en) | 2020-02-19 | 2022-08-30 | Suzhou MicroPort Orthobot Co., Ltd. | Method, system for registration of bone, and trackable element |
CN112155732A (zh) * | 2020-09-29 | 2021-01-01 | 苏州微创畅行机器人有限公司 | 可读存储介质、骨建模配准系统及骨科手术系统 |
WO2022068341A1 (zh) * | 2020-09-29 | 2022-04-07 | 苏州微创畅行机器人有限公司 | 可读存储介质、骨建模配准系统及骨科手术系统 |
CN112155732B (zh) * | 2020-09-29 | 2022-05-17 | 苏州微创畅行机器人有限公司 | 可读存储介质、骨建模配准系统及骨科手术系统 |
CN113081265A (zh) * | 2021-03-24 | 2021-07-09 | 重庆博仕康科技有限公司 | 手术导航空间配准方法、装置以及手术导航系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109692041B (zh) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10271822B2 (en) | Sensor coordinate calibration in an ultrasound system | |
CN109692041A (zh) | 一种针对覆盖软骨的骨表面配准方法 | |
EP2961322B1 (en) | Segmentation of large objects from multiple three-dimensional views | |
CN110163897B (zh) | 一种基于合成超声图像的多模态图像配准方法 | |
CN100456323C (zh) | 三维图像的快速配准方法 | |
CN106934821B (zh) | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 | |
CN111161333B (zh) | 一种肝脏呼吸运动模型的预测方法、装置及存储介质 | |
CN104021547A (zh) | 肺部 ct 的三维配准方法 | |
JP2016536035A (ja) | 変形可能レジストレーションのためのシステム及びモジュラーネットワーク | |
CN114155376A (zh) | 目标特征点提取方法、装置、计算机设备和存储介质 | |
Schumann et al. | An integrated system for 3D hip joint reconstruction from 2D X-rays: a preliminary validation study | |
Xie et al. | Feature‐based rectal contour propagation from planning CT to cone beam CT | |
CN113515875B (zh) | 基于多模态图像的骨生物力学建模方法、系统及装置 | |
Lorenz et al. | Automated abdominal plane and circumference estimation in 3D US for fetal screening | |
Maier-Hein et al. | Accounting for anisotropic noise in fine registration of time-of-flight range data with high-resolution surface data | |
CN108428245A (zh) | 基于自适应正则项的滑移图像配准方法 | |
CN116650115A (zh) | 一种基于uwb标记点的骨科手术导航注册方法 | |
CN115311258A (zh) | 一种spect平面图像中自动分割器官的方法和系统 | |
Haneishi et al. | Image synthesis using a mini gamma camera and stereo optical cameras | |
CN111166373B (zh) | 一种定位注册方法、装置和系统 | |
Yang et al. | A novel neurosurgery registration pipeline based on heat maps and anatomic facial feature points | |
Schumann et al. | Determination of pelvic orientation from ultrasound images using patch-SSMs and a hierarchical speed of sound compensation strategy | |
CN116012526B (zh) | 一种基于二维影像的三维ct影像病灶重建方法 | |
Dos Santos et al. | Minimally deformed correspondences between surfaces for intra-operative registration | |
Garg et al. | A computational approach to pre-align point cloud data for surface registration in image guided liver surgery |
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 | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: Room 305 and 306, building 8, No. 181, Wuchang Avenue, Wuchang Street, Yuhang District, Hangzhou, Zhejiang 311100 Patentee after: Hangzhou Jianjia Medical Technology Co.,Ltd. Address before: 310000 room 201-15, building 9, No. 181, Wuchang Avenue, Wuchang Street, Yuhang District, Hangzhou City, Zhejiang Province Patentee before: HANGZHOU JIANJIA ROBOT Co.,Ltd. |