一种用于股骨手术的图像配准系统及方法与流程

1.本发明涉及手术图像配准领域，特别是涉及一种用于股骨手术的图像配准系统及方法。


背景技术：

2.股骨颈需要承担人体大部分活动时的运动和支撑功能，其骨折的治疗效果将直接影响患者术后的生活质量。目前，股骨颈骨折已占全身骨折的3.6％，占髋部骨折的一半，多发于骨质疏松的中老年人，随着社会老龄化严重，其发病率逐年递增，手术导航系统辅助股骨颈骨折手术具有重要的意义。
3.目前市面上根据医学影像数据的不同股骨手术中常用的导航方式主要分为以下四种：1)基于ct数据的术前导航系统；2)基于c型臂x光机的术中成像导航系统；3)基于术中三维重建的导航系统；4)基于三维-二维图像配准的导航系统。
4.基于ct数据的术前导航系统由于没有术中信息的支持，无法进行术中的空间定位，基于c型臂x光机的术中成像导航系统又不够直观，术中三维重建系统价格高且普及率低。而基于三维-二维图像配准的导航系统可以有效避免以上缺点，难点在于对三维图像和二维图像的精确配准。


技术实现要素：

5.本发明的目的是提供一种用于股骨手术的图像配准系统及方法，提高了图像配准的精度，能够保证手术辅助导航的精度，进而改善了手术质量。
6.为实现上述目的，本发明提供了如下方案：
7.一种用于股骨手术的图像配准系统，包括：投影模块、预处理模块以及配准模块；
8.所述投影模块，用于对三维股骨ct图像进行不同角度的drr投影，得到若干张二维drr投影图；
9.所述预处理模块，用于对各二维drr投影图和各二维x光图像中的股骨图像进行预处理，生成股骨坐标系和股骨特征点；所述二维x光图像包括不同角度下对股骨部位拍摄的多张x光平面图像；
10.所述配准模块，用于根据所述二维drr投影图和所述二维x光图像的股骨坐标系和股骨特征点，从所述若干张二维drr投影图中确定与所述二维x光图像中股骨位姿相同的二维drr投影图，得到多张目标drr投影图。
11.可选地，所述股骨包括股骨头、股骨颈和股骨干，所述预处理模块包括股骨轮廓线生成单元、股骨头拟合单元、股骨颈特征线选取单元、股骨坐标系建立单元、股骨特征圆生成单元、股骨特征线选取单元和股骨特征点选取单元；
12.所述股骨轮廓线生成单元，用于自动提取股骨外轮廓，得到股骨轮廓线；
13.所述股骨头拟合单元，用于将股骨头拟合为一个圆，得到股骨头拟合圆；
14.所述股骨颈特征线选取单元，用于在股骨颈处的两条股骨轮廓线之间选取距离最
短的两点连线作为股骨颈特征线；
15.所述股骨坐标系建立单元，用于以所述股骨头拟合圆圆心和所述股骨颈特征线中点的连线为股骨坐标系的一个坐标轴，建立股骨坐标系；
16.所述股骨特征圆生成单元，用于以所述股骨头拟合圆半径为半径，在所述股骨头拟合圆圆心和所述股骨颈特征线中点的连线的延长线上选取两点作为圆心，分别生成第一特征圆和第二特征圆；所述第一特征圆与所述股骨头拟合圆相切，所述第一特征圆的圆心在所述第二特征圆的边缘线上；
17.所述股骨特征线选取单元，用于选取多条垂直于所述股骨头拟合圆圆心和所述股骨颈特征线中点连线的股骨特征线；
18.所述股骨特征点选取单元，用于取各股骨特征线与所述股骨轮廓线的交点为股骨特征点。
19.可选地，所述股骨特征线选取单元包括：第一股骨特征线选取子单元、第二股骨特征线选取子单元、第三股骨特征线选取子单元、第四股骨特征线选取子单元和第五股骨特征线选取子单元；
20.所述第一股骨特征线选取子单元，用于选取第一股骨特征线；所述第一股骨特征线为过所述股骨头拟合圆和所述第一特征圆切点的线段；
21.所述第二股骨特征线选取子单元，用于选取第二股骨特征线；所述第二股骨特征线为过所述第一特征圆圆心的线段；
22.所述第三股骨特征线选取子单元，用于选取第三股骨特征线；所述第三股骨特征线为过所述第二特征圆圆心的线段；
23.所述第四股骨特征线选取子单元，用于选取第四股骨特征线；所述第四股骨特征线为平行于所述第二股骨特征线，且为所述第二特征圆切线的线段；
24.所述第五股骨特征线选取子单元，用于选取第五股骨特征线；所述第五股骨特征线为过所述股骨头拟合圆圆心和所述股骨颈特征线中点连线的延长线与所述第三股骨特征线和所述第四股骨特征线交点连线的中点的线段。
25.可选地，所述配准模块包括：股骨坐标系重合单元和股骨特征点配准单元；
26.所述股骨坐标系重合单元，用于将所述二维drr投影图和所述二维x光图像的股骨坐标系重合；
27.所述股骨特征点配准单元，用于将所述二维drr投影图和所述二维x光图像的股骨特征点位置一一匹配。
28.可选地，所述预处理模块还包括：感兴趣区域选取单元；
29.所述感兴趣区域选取单元，用于在所述股骨轮廓线的股骨头轮廓线和两段股骨颈轮廓线的位置各选取一块由线段和圆弧组成的感兴趣区域。
30.另一方面，针对前述的图像配准系统，本发明还提供了一种用于股骨手术的图像配准方法，包括以下步骤：
31.对三维股骨ct图像进行不同角度的drr投影，得到若干张二维drr投影图；
32.对各二维drr投影图和各二维x光图像中的股骨图像进行预处理，生成股骨坐标系和股骨特征点；所述二维x光图像包括不同角度下对股骨部位拍摄的多张x光平面图像；
33.根据所述二维drr投影图和所述二维x光图像的股骨坐标系和股骨特征点，从所述
若干张二维drr投影图中确定与所述二维x光图像中股骨位姿相同的二维drr投影图，得到多张目标drr投影图。
34.可选地，所述股骨包括股骨头、股骨颈和股骨干，所述对各二维drr投影图和各二维x光图像中的股骨图像进行预处理，生成股骨坐标系和股骨特征点，具体包括：
35.自动提取股骨外轮廓，得到股骨轮廓线；
36.将股骨头拟合为一个圆，得到股骨头拟合圆；
37.在股骨颈处的两条股骨轮廓线之间选取距离最短的两点连线作为股骨颈特征线；
38.以所述股骨头拟合圆圆心和所述股骨颈特征线中点的连线为股骨坐标系的一个坐标轴，建立股骨坐标系；
39.以所述股骨头拟合圆半径为半径，在所述股骨头拟合圆圆心和所述股骨颈特征线中点的连线的延长线上选取两点作为圆心，分别生成第一特征圆和第二特征圆；所述第一特征圆与所述股骨头拟合圆相切，所述第一特征圆的圆心在所述第二特征圆的边缘线上；
40.生成多条垂直于所述股骨头拟合圆圆心和所述股骨颈特征线中点连线的股骨特征线；
41.取各股骨特征线与所述股骨轮廓线的交点为股骨特征点。
42.可选地，所述生成多条垂直于所述股骨头拟合圆圆心和所述股骨颈特征线中点连线的股骨特征线，具体包括：
43.选取过所述股骨头拟合圆和所述第一特征圆切点的线段为第一股骨特征线；
44.选取过所述第一特征圆圆心的线段为第二股骨特征线；
45.选取过所述第二特征圆圆心的线段为第三股骨特征线；
46.选取平行于所述第二股骨特征线，且为所述第二特征圆切线的线段为第四股骨特征线；
47.选取过所述股骨头拟合圆圆心和所述股骨颈特征线中点连线的延长线与所述第三股骨特征线和所述第四股骨特征线交点连线的中点的线段为第五股骨特征线。
48.可选地，所述根据所述二维drr投影图和所述二维x光图像的股骨坐标系和股骨特征点，从所述若干张二维drr投影图中确定与所述二维x光图像中股骨位姿相同的二维drr投影图，得到多张目标drr投影图，具体包括：
49.将所述二维drr投影图和所述二维x光图像的股骨坐标系重合；
50.将所述二维drr投影图和所述二维x光图像的股骨特征点位置一一匹配。
51.可选地，所述方法还包括：
52.在所述股骨轮廓线的股骨头轮廓线和两段股骨颈轮廓线的位置各选取一块由线段和圆弧组成的感兴趣区域；
53.对各所述感兴趣区域中的轮廓线进行识别，生成所述股骨头轮廓线和两段所述股骨颈轮廓线的拟合曲线；
54.根据股骨头处的股骨轮廓线的拟合曲线，自动生成股骨头拟合圆；
55.在股骨颈处的两段股骨轮廓线的拟合曲线之间，选取距离最短的两点连线作为股骨颈特征线。
56.根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供的一种用于股骨手术的图像配准系统及方法，包括：投影模块、预处理模块以及配准模块；投影模
块，用于对三维股骨ct图像进行不同角度的drr投影，得到若干张二维drr投影图；预处理模块，用于对各二维drr投影图和各二维x光图像中的股骨图像进行预处理，生成股骨坐标系和股骨特征点；二维x光图像包括不同角度下对股骨部位拍摄的多张x光平面图像；配准模块，用于根据二维drr投影图和二维x光图像的股骨坐标系和股骨特征点，从若干张二维drr投影图中确定与二维x光图像中股骨位姿相同的二维drr投影图，得到多张目标drr投影图。本发明在进行手术之前将三维股骨ct图像进行各种角度的drr投影，分别在二维drr投影图与二维x光图像上进行股骨坐标系和股骨特征点的生成和选取，根据股骨坐标系和股骨特征点，对二维drr投影图和二维x光图像进行配准，提高了图像配准的精度，能够准确得到与进行二维x光图像扫描时的股骨位姿相同的二维drr投影图，进而可以获得真实的手术路径以进行手术导航，保证手术辅助导航的精度，进而改善了手术质量。
附图说明
57.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
58.图1为本发明提供的一种用于股骨手术的图像配准系统的结构框图；
59.图2为本发明提供的图像配准系统中股骨特征线选取单元的结构框图；
60.图3为本发明提供的一种用于股骨手术的图像配准方法的流程图；
61.图4为本发明提供的图像配准方法中步骤s2的具体流程图；
62.图5为本发明提供的图像配准方法中步骤s26的具体流程图；
63.图6为本发明提供的图像配准方法中重建的三维股骨ct图像；
64.图7为本发明提供的图像配准方法中对二维drr投影图生成的股骨轮廓线示意图；
65.图8为本发明提供的图像配准方法中对二维drr投影图生成的股骨特征点示意图；
66.图9为本发明提供的图像配准方法中扫描的二维x光图像；
67.图10为本发明提供的图像配准方法中对二维x光图像生成的股骨轮廓线示意图；
68.图11为本发明提供的图像配准方法中对二维x光图像生成的股骨特征点示意图。
69.图中1：投影模块；2：预处理模块；3：配准模块；21：股骨轮廓线生成单元；22：股骨头拟合单元；23：股骨颈特征线选取单元；24：股骨坐标系建立单元；25：股骨特征圆生成单元；26：股骨特征线选取单元；27：股骨特征点选取单元；28：感兴趣区域选取单元；261：第一股骨特征线选取子单元；262：第二股骨特征线选取子单元；263：第三股骨特征线选取子单元；264：第四股骨特征线选取子单元；265：第五股骨特征线选取子单元；31：股骨坐标系重合单元；32：股骨特征点配准单元。
具体实施方式
70.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
71.本发明的目的是提供一种用于股骨手术的图像配准系统及方法，提高了图像配准的精度，能够保证手术辅助导航的精度，进而改善了手术质量。
72.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
73.实施例1：
74.本发明提供了一种用于股骨手术的图像配准系统，如图1～图2所示，本发明用于股骨手术的图像配准系统包括：投影模块1、预处理模块2以及配准模块3；
75.投影模块1，用于对三维股骨ct图像进行不同角度的drr投影，得到若干张二维drr投影图；
76.预处理模块2，用于对各二维drr投影图和各二维x光图像中的股骨图像进行预处理，生成股骨坐标系和股骨特征点；二维x光图像包括不同角度下对股骨部位拍摄的多张x光平面图像；
77.配准模块3，用于根据二维drr投影图和二维x光图像的股骨坐标系和股骨特征点，从若干张二维drr投影图中确定与二维x光图像中股骨位姿相同的二维drr投影图，得到多张目标drr投影图。
78.股骨包括股骨头、股骨颈和股骨干，本实施例中，预处理模块2包括股骨轮廓线生成单元21、股骨头拟合单元22、股骨颈特征线选取单元23、股骨坐标系建立单元24、股骨特征圆生成单元25、股骨特征线选取单元26和股骨特征点选取单元27；
79.股骨轮廓线生成单元21，用于自动提取股骨外轮廓，得到股骨轮廓线；
80.股骨头拟合单元22，用于将股骨头拟合为一个圆，得到股骨头拟合圆；
81.股骨颈特征线选取单元23，用于在股骨颈处的两条股骨轮廓线之间选取距离最短的两点连线作为股骨颈特征线；
82.股骨坐标系建立单元24，用于以股骨头拟合圆圆心和股骨颈特征线中点的连线为股骨坐标系的一个坐标轴，建立股骨坐标系；
83.股骨特征圆生成单元25，用于以股骨头拟合圆半径为半径，在股骨头拟合圆圆心和股骨颈特征线中点的连线的延长线上选取两点作为圆心，分别生成第一特征圆和第二特征圆；第一特征圆与股骨头拟合圆相切，第一特征圆的圆心在第二特征圆的边缘线上；
84.股骨特征线选取单元26，用于选取多条垂直于股骨头拟合圆圆心和股骨颈特征线中点连线的股骨特征线；
85.股骨特征点选取单元27，用于取各股骨特征线与股骨轮廓线的交点为股骨特征点。
86.本实施例中，股骨特征线选取单元26包括：第一股骨特征线选取子单元261、第二股骨特征线选取子单元262、第三股骨特征线选取子单元263、第四股骨特征线选取子单元264和第五股骨特征线选取子单元265；
87.第一股骨特征线选取子单元261，用于选取第一股骨特征线；第一股骨特征线为过股骨头拟合圆和第一特征圆切点的线段；
88.第二股骨特征线选取子单元262，用于选取第二股骨特征线；第二股骨特征线为过第一特征圆圆心的线段；
89.第三股骨特征线选取子单元263，用于选取第三股骨特征线；第三股骨特征线为过
第二特征圆圆心的线段；
90.第四股骨特征线选取子单元264，用于选取第四股骨特征线；第四股骨特征线为平行于第二股骨特征线，且为第二特征圆切线的线段；
91.第五股骨特征线选取子单元265，用于选取第五股骨特征线；第五股骨特征线为过股骨头拟合圆圆心和股骨颈特征线中点连线的延长线与第三股骨特征线和第四股骨特征线交点连线的中点的线段。
92.具体到本实施例中，配准模块3包括：股骨坐标系重合单元31和股骨特征点配准单元32；
93.股骨坐标系重合单元31，用于将二维drr投影图和二维x光图像的股骨坐标系重合；
94.股骨特征点配准单元32，用于将二维drr投影图和二维x光图像的股骨特征点位置一一匹配。
95.为了减少计算量，在本实施例中，预处理模块2还包括：感兴趣区域选取单元28；
96.感兴趣区域选取单元28，用于在股骨轮廓线的股骨头轮廓线和两段股骨颈轮廓线的位置各选取一块由线段和圆弧组成的感兴趣区域。由线段和圆弧组成的选择模板在股骨头和股骨颈处进行感兴趣区域选择，这样不需要对整个x片进行轮廓识别，从而达到了减少计算量的效果。
97.实施例2：
98.对应于实施例1的用于股骨手术的图像配准系统，本发明还提供了一种用于股骨手术的图像配准方法，如图3所示，图像配准方法包括以下步骤：
99.s1、对三维股骨ct图像进行不同角度的drr投影，得到若干张二维drr投影图；
100.s2、对各二维drr投影图和各二维x光图像中的股骨图像进行预处理，生成股骨坐标系和股骨特征点；二维x光图像包括不同角度下对股骨部位拍摄的多张x光平面图像；
101.s3、根据二维drr投影图和二维x光图像的股骨坐标系和股骨特征点，从若干张二维drr投影图中确定与二维x光图像中股骨位姿相同的二维drr投影图，得到多张目标drr投影图。步骤s3具体包括：
102.s31、将二维drr投影图和二维x光图像的股骨坐标系重合；
103.s32、将二维drr投影图和二维x光图像的股骨特征点位置一一匹配。
104.股骨包括股骨头、股骨颈和股骨干，ct或者x片只能显示股骨干一部分。为了进行ct图像和x光图像之间的精确配准，在本实施例中，如图4所示，步骤s2具体包括：
105.s21、自动提取股骨外轮廓，得到股骨轮廓线；
106.s22、将股骨头拟合为一个圆，得到股骨头拟合圆；
107.s23、在股骨颈处的两条股骨轮廓线之间选取距离最短的两点连线作为股骨颈特征线；
108.s24、以股骨头拟合圆圆心和股骨颈特征线中点的连线为股骨坐标系的一个坐标轴，建立股骨坐标系；
109.s25、以股骨头拟合圆半径为半径，在股骨头拟合圆圆心和股骨颈特征线中点的连线的延长线上选取两点作为圆心，分别生成第一特征圆和第二特征圆；第一特征圆与股骨头拟合圆相切，第一特征圆的圆心在第二特征圆的边缘线上；
110.s26、生成多条垂直于股骨头拟合圆圆心和股骨颈特征线中点连线的股骨特征线；
111.s27、取各股骨特征线与股骨轮廓线的交点为股骨特征点。
112.关于各股骨特征线的选取规则，在本实施例中，如图5所示，步骤s26具体包括：
113.s261、选取过股骨头拟合圆和第一特征圆切点的线段为第一股骨特征线；
114.s262、选取过第一特征圆圆心的线段为第二股骨特征线；
115.s263、选取过第二特征圆圆心的线段为第三股骨特征线；
116.s264、选取平行于第二股骨特征线，且为第二特征圆切线的线段为第四股骨特征线；
117.s265、选取过股骨头拟合圆圆心和股骨颈特征线中点连线的延长线与第三股骨特征线和第四股骨特征线交点连线的中点的线段为第五股骨特征线。
118.为了减少计算量，本实施例的图像配准方法还包括：
119.a1、在股骨轮廓线的股骨头轮廓线和两段股骨颈轮廓线的位置各选取一块由线段和圆弧组成的感兴趣区域；由线段和圆弧组成的选择模板在股骨头和股骨颈处进行感兴趣区域选择，这样不需要对整个x片进行轮廓识别，从而达到了减少计算量的效果。
120.a2、对各感兴趣区域中的轮廓线进行识别，生成股骨头轮廓线和两段股骨颈轮廓线的拟合曲线。
121.a3、根据股骨头处的股骨轮廓线的拟合曲线，自动生成股骨头拟合圆。
122.a4、在股骨颈处的两段股骨轮廓线的拟合曲线之间，选取距离最短的两点连线作为股骨颈特征线。
123.下面结合一个具体的例子，来对本发明提供的用于股骨手术的图像配准方法作进一步详细的说明，包括以下步骤：
124.1)在手术开始前，通过ct影像设备获取患者股骨图像数据，并处理重建生成三维股骨ct图像，如图6，在其上规划手术路径。
125.2)根据手术需求在侧位单一视角z方向上对三维股骨ct图像进行多次投影，如图7左侧所示，o为股骨头拟合圆的圆心，oy过股骨颈特征线交股骨干中心轴于o1，oz垂直于oy且在oo1y1平面内，ox轴符合右手法则。
126.3)绕x、y轴每间隔预设旋转度数，α、β均在1
°‑2°
范围内，从z轴投影生成多张drr图像，获取对应的位姿参数p＝(θ
x
,θy,θz,p
x
,py,pz)，θ
x
、θy、θz分别表示参考坐标系下旋转角度，p
x
、py、pz分别表示参考坐标系下沿三个主轴的平移量；类似的也在正位单一视角x方向上进行多次投影，如图7右侧所示，得到多张二维drr投影图。
127.4)将每张二维drr投影图进行处理，如图8，自动提取股骨轮廓线，自动生成股骨头拟合圆，圆心为o2，o3为股骨颈特征线的中点，在o2o3延长线上作第一特征圆o4和第二特征圆o5，选取五条股骨特征线交股骨轮廓线于股骨特征点a1、a2、a3、a4、a5、a6、a7、a8、a9和a10。
128.5)在患者躺倒手术床要进行手术时，获取c型臂x线机不同角度拍摄的二维x光图像，如图9，两张二维图像拍摄夹角为90
°
，图9左侧所示为正位x光图像，即患者平躺，垂直拍摄的x片，图9右侧所示为侧位x光图像，即患者平躺，水平拍摄的x片。
129.6)如图10所示，以正位x光图像为例，在股骨头处划定感兴趣区域c1自动生成股骨头拟合圆c2，或在区域c1内手动选取三点，生成股骨头拟合圆c2，圆心为o6。手动选取股骨
颈感兴趣区域c3和c5，自动生成股骨颈外轮廓拟合曲线c4和c6，手动选取股骨颈最窄处两条拟合曲线，计算两条拟合曲线的最小距离，或手动选取目测股骨颈最窄处的两点，以最小距离作为股骨颈特征线，以股骨颈特征线为直径生成股骨颈最小圆c7，圆心为o7。参考操作步骤4)，自动对二维x光图像生成股骨坐标系、股骨特征圆、股骨特征线和股骨特征点，如图11。
130.7)对股骨每张正侧位二维drr投影与正侧位二维x光图像进行自动配准，将二维drr投影图与二维x光图像的股骨坐标系重合，计算所有对应股骨特征点的距离，进行特征点相似度比较，在误差范围内如果有二维drr投影图和二维x光图像匹配，则获取该二维drr投影图的投影位姿参数，若配准结果不止一组，选取误差最小一组，完成股骨三维股骨ct图像与二维x光图像的配准。
131.8)根据二维drr投影图的投影位姿参数和手术规划路径信息，利用双平面法输出精确配准后的手术路径。
132.本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。