CN109692041B - 一种针对覆盖软骨的骨表面配准方法 - Google Patents
一种针对覆盖软骨的骨表面配准方法 Download PDFInfo
- Publication number
- CN109692041B CN109692041B CN201910090929.0A CN201910090929A CN109692041B CN 109692041 B CN109692041 B CN 109692041B CN 201910090929 A CN201910090929 A CN 201910090929A CN 109692041 B CN109692041 B CN 109692041B
- Authority
- CN
- China
- Prior art keywords
- bone
- bone surface
- cartilage
- sample points
- coordinate system
- 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
- 210000000988 bone and bone Anatomy 0.000 title claims abstract description 65
- 210000000845 cartilage Anatomy 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000013519 translation Methods 0.000 claims abstract description 20
- 238000013507 mapping Methods 0.000 claims abstract description 6
- 239000000523 sample Substances 0.000 claims description 61
- 239000013074 reference sample Substances 0.000 claims description 6
- 208000006735 Periostitis Diseases 0.000 description 4
- 210000003460 periosteum Anatomy 0.000 description 4
- 210000004872 soft tissue Anatomy 0.000 description 3
- 210000000689 upper leg Anatomy 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000002604 ultrasonography Methods 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
-
- 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;
上述方法在多数情况下能够获得精准的配准结果,但是并不适用于股骨这类覆盖软骨的场景。一方面,在医学影像中,骨头边缘灰度特征明显,但软骨及骨膜厚度小且灰度与软组织类似,因此,骨头边缘点易于获取,而骨膜与软骨边缘难以界定。借助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 CN109692041A (zh) | 2019-04-30 |
CN109692041B true 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) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111227932B (zh) * | 2020-02-19 | 2021-06-04 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
CN111297475B (zh) * | 2020-02-19 | 2021-06-18 | 苏州微创畅行机器人有限公司 | 骨注册方法、骨注册系统、骨注册控制装置及可跟踪元件 |
CN112155732B (zh) * | 2020-09-29 | 2022-05-17 | 苏州微创畅行机器人有限公司 | 可读存储介质、骨建模配准系统及骨科手术系统 |
CN113081265B (zh) * | 2021-03-24 | 2022-11-15 | 重庆博仕康科技有限公司 | 手术导航空间配准方法、装置以及手术导航系统 |
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 |
Also Published As
Publication number | Publication date |
---|---|
CN109692041A (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109692041B (zh) | 一种针对覆盖软骨的骨表面配准方法 | |
US10561403B2 (en) | Sensor coordinate calibration in an ultrasound system | |
EP2961322B1 (en) | Segmentation of large objects from multiple three-dimensional views | |
CN106934821B (zh) | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 | |
US11257241B2 (en) | System and method for component positioning by registering a 3D patient model to an intra-operative image | |
JP2003265408A (ja) | 内視鏡誘導装置および方法 | |
EP3463097B1 (en) | Correcting probe induced deformation in an ultrasound fusing imaging system | |
CN112509022A (zh) | 一种术前三维影像与术中透视图像的无标定物配准方法 | |
CN105046644B (zh) | 基于线性相关性的超声与ct图像配准方法和系统 | |
CN109949349A (zh) | 一种多模态三维图像的配准及融合显示方法 | |
CN103337065A (zh) | 小鼠三维ct图像的非刚性配准方法 | |
CN106296613B (zh) | 一种基于dr机器的双能量减影方法 | |
CN116650115A (zh) | 一种基于uwb标记点的骨科手术导航注册方法 | |
CN113515875B (zh) | 基于多模态图像的骨生物力学建模方法、系统及装置 | |
EP3137917A1 (en) | High resolution intraoperative mri images | |
CN111640127B (zh) | 一种用于骨科的精准临床诊断导航方法 | |
CN108804861B (zh) | 一种具有真实力反馈的脊柱微创手术培训系统及方法 | |
Boussaid et al. | 3D Model-based reconstruction of the proximal femur from low-dose biplanar x-ray images | |
CN115137342A (zh) | 基于深度学习的图像采集医学系统和方法 | |
Hohlmann et al. | The interleaved partial active shape model (IPASM) search algorithm–towards 3D ultrasound-based bone surface reconstruction | |
CN105678794A (zh) | 一种基于图像处理的瘢痕治疗效果的评估方法 | |
Zhang et al. | Performance analysis of active shape reconstruction of fractured, incomplete skulls | |
US20060111630A1 (en) | Method of tomographic imaging | |
Dekomien et al. | System architecture for intraoperative ultrasound registration in image-based medical navigation | |
CN111932443A (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 | ||
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. |