CN106023140B - 一种术前三维图像与术中三维图像的无标志物配准方法 - Google Patents
一种术前三维图像与术中三维图像的无标志物配准方法 Download PDFInfo
- Publication number
- CN106023140B CN106023140B CN201610293632.0A CN201610293632A CN106023140B CN 106023140 B CN106023140 B CN 106023140B CN 201610293632 A CN201610293632 A CN 201610293632A CN 106023140 B CN106023140 B CN 106023140B
- Authority
- CN
- China
- Prior art keywords
- image
- locator
- preoperative
- mutual information
- benchmark
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Processing Or Creating Images (AREA)
Abstract
本发明涉及一种术前三维图像与术中三维图像的无标志物配准方法,包括以下步骤:1)设置一配准***,其包括***、数据计算模块和互信息图像配准模块;2)术前配置阶段,***获取术前三维图像的位置;3)判断***是否记录了基准图像及其位置;4)数据计算模块计算***与C型臂中心的误差,并将计算出误差的1.5‑2倍作为搜索边界,获取手术中三维图像并保存至搜索边界;5)判断手术中***获取三维图像位置是否成功;6)数据计算模块计算手术中三维图像位置与基准图像位置的偏移距离;7)由互信息配准模块对手术中三维图像和基准图像进行互信息精准配准,然后输出精准配准后的图像。本发明配准成本低、配准速度快。
Description
技术领域
本发明涉及一种无标志物配准方法,特别是关于一种医学领域中使用的术前三维图像与术中三维图像的无标志物配准方法。
背景技术
目前,随着CBCT、o-Arm等手术室三维影像设备的普及,术中获取三维影像成为可能,但此类设备多是移动布置的,移动前后的图像不在同一个坐标系下,如需同时观察就需要对其进行配准。一般配准方法主要按照有无患者标志物分为两大类,患者标志物通常是固定在患者被扫描部分的一个或多个标志,可以在图像上清楚识别,但在手术环境中通常难以安装患者标识物。常规的无标志物三维配准方法多是基于互信息或是基于高精度定位设备的。
基于互信息的配准方法,是利用多组图之间的自然特征进行搜索的,配准成功与否及配准速度受两幅图的初始位置影响,一旦初始位置不良进行互信息搜索时就会占用很长世间甚至无法配准。这个方法一般用于CT工作站上进行术前术后图像的分析比对,对时间效率要求不严格且操作医生多会人工调整图像的初始位置使其近似配准。但手术过程中时间紧,且由于消毒的原因医生不方便精细的操作来手工调整图像初始位置。
基于高精度定位设备是利用成像设备成的三维影像的空间位置与成像设备位置存在固定关系,使用光电或电磁定位设备跟踪成像设备的位置,进而计算出图像的位置的,此方法方便快捷,但配准精度受***本身精度限制,高精度***价格不菲且一般是光电原理的。光电***精度与定位标尺尺寸正相关,维持必要精度的标尺大小一般较大,导致较难被***有效范围覆盖,实际上缩小了成像设备可以配准的范围。
发明内容
针对上述问题,本发明的目的是提供一种术前三维图像与术中三维图像的无标志物配准方法,其配准成本低、配准速度快。
为实现上述目的,本发明采取以下技术方案:一种术前三维图像与术中三维图像的无标志物配准方法,其特征在于它包括以下步骤:1)设置一配准***,其包括***、数据计算模块和互信息图像配准模块;2)术前配置阶段,将***固定在CBCT的C型臂可旋转部,启动配准***,获取一个系列的三维图像,***获取术前三维图像的位置;3)判断***是否记录了基准图像及其位置,如果记录了基准图像及其位置则进行下一步,否则,记录术前三维图像及其位置作为基准图像和基准图像位置并且输出基准图像,进入下一步;4)数据计算模块计算***与C型臂中心的误差,并将计算出误差的1.5-2倍作为搜索边界,获取手术中三维图像并保存至搜索边界;5)判断手术中***获取三维图像位置是否成功,如果***定位成功则进行下一步,否则直接输出手术中的原三维图像,本次术前三维图像与术中三维图像的无标志物配准结束;6)数据计算模块计算手术中三维图像位置与基准图像位置的偏移距离,并将计算结果作为初始参数传入互信息配准模块;7)由互信息配准模块对手术中三维图像和基准图像进行互信息精准配准,然后输出精准配准后的图像。
优选地,所述步骤4)中,***与C型臂的中心的误差,根据三角函数计算误差e:
式中,L为***输出的坐标中心与C型臂的中心的距离,单位为mm;R为***最大角度误差。
优选地,所述步骤7)中,互信息配准模块的配准方法,其步骤如下:1)互信息配准模块包括图像移动器、图像评估器和图像移动优化器;2)获取基准图像和手术中三维图像,用图像移动器移动两组图像,使基准图像和手术中三维图像相同地方相互重叠;3)图像评估器对移动的两组图像效果进行评估;4)判断移动的两组图像效果是否达到了精度要求,如果达到要求的精度则直接输出配准图像,否则使用图像移动优化器对两组图像均进行优化,并返回步骤2)。
优选地,所述图像移动器采用欧拉角刚体移动器。
优选地,所述***是非高精度的***,所述***包括光学、电磁学和惯性***,所述***输出的结果能转化为表达位置的点(x,y,z)和相应的欧拉角(a,b,r)或转成4x4矩阵的形式。
本发明由于采取以上技术方案,其具有以下优点:1、本发明采用了术前三维图像与术中三维图像的无标志物配准方法,有效降低了对***精度的要求。2、本发明采用了术前三维图像与术中三维图像的无标志物配准方法,加快了基于互信息配准的速度。3、本发明采用了术前三维图像与术中三维图像的无标志物配准方法,在有限的范围搜索,能有效防止基于互信息的配准陷入局部伪最优解而返回错误信息。
附图说明
图1是本发明整体流程示意图;
图2是本发明互信息配准流程示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供一种术前三维图像与术中三维图像的无标志物配准方法,其包括以下步骤:
1)设置一配准***,其包括***、数据计算模块和互信息图像配准模块;
2)术前配置阶段,将***固定在CBCT的C型臂可旋转部,启动配准***,获取一个系列的三维图像,***获取术前三维图像的位置;
3)判断***是否记录了基准图像及其位置,如果记录了基准图像及其位置则进行下一步,否则,记录术前三维图像及其位置作为基准图像和基准图像位置并且输出基准图像,进入下一步;
4)数据计算模块计算***与C型臂中心的误差,并将计算出误差的1.5-2倍作为搜索边界,获取手术中三维图像并保存至搜索边界;
5)判断手术中***获取三维图像位置是否成功,如果***定位成功则进行下一步,否则直接输出手术中的原三维图像,本次术前三维图像与术中三维图像的无标志物配准结束;其中,***定位失败等异常处理全部静默进行,不会影响手术的正常操作;
6)数据计算模块计算手术中三维图像位置与基准图像位置的偏移距离,并将计算结果作为初始参数传入互信息配准模块;
7)由互信息配准模块对手术中三维图像和基准图像进行互信息精准配准,然后输出精准配准后的图像,完成术前三维图像与术中三维图像的无标志物配准。
上述步骤4)中,***与C型臂的中心的误差,根据三角函数计算误差e:
式中,L为***输出的坐标中心与C型臂的中心的距离,单位为mm;R为***最大角度误差。
上述步骤7)中,如图2所示,本发明还提供互信息配准模块的配准方法,其步骤如下:
1)互信息配准模块包括图像移动器、图像评估器和图像移动优化器;
2)获取基准图像和手术中三维图像,用图像移动器移动两组图像,使基准图像和手术中三维图像相同地方相互重叠;
3)图像评估器对移动的两组图像效果进行评估;
4)判断移动的两组图像效果是否达到了精度要求,如果达到要求的精度则直接输出配准图像,否则使用图像移动优化器对两组图像均进行优化,并返回步骤2)。
上述实施例中,图像移动器采用欧拉角刚体移动器。
上述各实施例中,***是非高精度的***,***包括但不限于光学、电磁学和惯性***,***输出的结果可以转化为表达位置的点(x,y,z)和相应的欧拉角(a,b,r)或转成4x4矩阵的形式。
上述各实施例仅用于说明本发明,各部件的结构、尺寸、设置位置及形状都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别部件进行的改进和等同变换,均不应排除在本发明的保护范围之外。
Claims (5)
1.一种术前三维图像与术中三维图像的无标志物配准方法,其特征在于:它包括以下步骤:
1)设置一配准***,其包括***、数据计算模块和互信息图像配准模块;
2)术前配置阶段,将***固定在CBCT的C型臂可旋转部,启动配准***,获取一个系列的三维图像,***获取术前三维图像的位置;
3)判断***是否记录了基准图像及其位置,如果记录了基准图像及其位置则进行下一步,否则,记录术前三维图像及其位置作为基准图像和基准图像位置并且输出基准图像,进入下一步;
4)数据计算模块计算***与C型臂中心的误差,并将计算出误差的1.5-2倍作为搜索边界,获取手术中三维图像并保存至搜索边界;
5)判断手术中***获取三维图像位置是否成功,如果***定位成功则进行下一步,否则直接输出手术中的原三维图像,本次术前三维图像与术中三维图像的无标志物配准结束;
6)数据计算模块计算手术中三维图像位置与基准图像位置的偏移距离,并将计算结果作为初始参数传入互信息图像配准模块;
7)由互信息图像配准模块对手术中三维图像和基准图像进行互信息精准配准,然后输出精准配准后的图像。
2.如权利要求1所述的一种术前三维图像与术中三维图像的无标志物配准方法,其特征在于:所述步骤4)中,***与C型臂的中心的误差,根据三角函数计算误差e:
式中,L为***输出的坐标中心与C型臂的中心的距离,单位为mm;R为***最大角度误差。
3.如权利要求1所述的一种术前三维图像与术中三维图像的无标志物配准方法,其特征在于:所述步骤7)中,互信息图像配准模块的配准方法,其步骤如下:
7.1)互信息图像配准模块包括图像移动器、图像评估器和图像移动优化器;
7.2)获取基准图像和手术中三维图像,用图像移动器移动两组图像,使基准图像和手术中三维图像相同地方相互重叠;
7.3)图像评估器对移动的两组图像效果进行评估;
7.4)判断移动的两组图像效果是否达到了精度要求,如果达到要求的精度则直接输出配准图像,否则使用图像移动优化器对两组图像均进行优化,并返回步骤7.2)。
4.如权利要求3所述的一种术前三维图像与术中三维图像的无标志物配准方法,其特征在于:所述图像移动器采用欧拉角刚体移动器。
5.如权利要求1-4任一项所述的一种术前三维图像与术中三维图像的无标志物配准方法,其特征在于:所述***是非高精度的***,所述***包括光学、电磁学和惯性***,所述***输出的结果能转化为表达位置的点(x,y,z)和相应的欧拉角(a,b,r)或转成4x4矩阵的形式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610293632.0A CN106023140B (zh) | 2016-05-05 | 2016-05-05 | 一种术前三维图像与术中三维图像的无标志物配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610293632.0A CN106023140B (zh) | 2016-05-05 | 2016-05-05 | 一种术前三维图像与术中三维图像的无标志物配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106023140A CN106023140A (zh) | 2016-10-12 |
CN106023140B true CN106023140B (zh) | 2018-12-21 |
Family
ID=57081797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610293632.0A Active CN106023140B (zh) | 2016-05-05 | 2016-05-05 | 一种术前三维图像与术中三维图像的无标志物配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106023140B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113573776A (zh) * | 2020-02-14 | 2021-10-29 | 西安大医集团股份有限公司 | 图像引导方法、装置、放疗设备和计算机存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6775404B1 (en) * | 1999-03-18 | 2004-08-10 | University Of Washington | Apparatus and method for interactive 3D registration of ultrasound and magnetic resonance images based on a magnetic position sensor |
CN103020976A (zh) * | 2012-12-31 | 2013-04-03 | 中国科学院合肥物质科学研究院 | 一种基于带权模糊互信息的三维医学图像配准方法及*** |
-
2016
- 2016-05-05 CN CN201610293632.0A patent/CN106023140B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6775404B1 (en) * | 1999-03-18 | 2004-08-10 | University Of Washington | Apparatus and method for interactive 3D registration of ultrasound and magnetic resonance images based on a magnetic position sensor |
CN103020976A (zh) * | 2012-12-31 | 2013-04-03 | 中国科学院合肥物质科学研究院 | 一种基于带权模糊互信息的三维医学图像配准方法及*** |
Non-Patent Citations (3)
Title |
---|
Fast 3D medical image rigid registration using the mutual information in equivalent meridian plane;Z T Lu等;《Journal of Photographic Science》;20131231;第58卷(第2期);第96-104页 * |
基于互信息的多分辨率三维脑图像配准方法;任海萍等;《基于互信息的多分辨率三维脑图像配准方法》;20021231;第19卷(第4期);第599-601、610页 * |
超声与磁共振图像三维配准方法研究;李婷婷;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160315(第03期);第1-51页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106023140A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107468350B (zh) | 一种三维图像专用标定器、手术定位***及定位方法 | |
CN107330926A (zh) | 一种术中导航***中的非标记医学图像配准***及方法 | |
EP2953569B1 (en) | Tracking apparatus for tracking an object with respect to a body | |
US11103314B2 (en) | Methods and devices for tracking objects by surgical navigation systems | |
US9173715B2 (en) | Ultrasound CT registration for positioning | |
EP2931159B1 (en) | Registration and navigation using a three-dimensional tracking sensor | |
EP3057064B1 (en) | Pose determination from a pattern of four leds | |
EP3116401B1 (en) | Anatomical image projection system | |
CN102266250A (zh) | 超声手术导航***及超声手术导航方法 | |
EP3212104B1 (en) | Hybrid navigation system for surgical interventions | |
Kyme et al. | Markerless motion estimation for motion-compensated clinical brain imaging | |
Liu et al. | On-demand calibration and evaluation for electromagnetically tracked laparoscope in augmented reality visualization | |
CN115038401A (zh) | 用于外部固定器的方法和装置 | |
CN106023140B (zh) | 一种术前三维图像与术中三维图像的无标志物配准方法 | |
CN203953815U (zh) | 辅助外科手术的术中导航*** | |
US11826107B2 (en) | Registration system for medical navigation and method of operation thereof | |
Xie et al. | Image-guided navigation system for minimally invasive total hip arthroplasty (MITHA) using an improved position-sensing marker | |
Haliburton et al. | A visual odometry base-tracking system for intraoperative C-arm guidance | |
TW201519856A (zh) | 手術定位方法 | |
US20240016550A1 (en) | Scanner for intraoperative application | |
WO2022198198A1 (en) | Method and system for non-contact patient registration in image‑guided surgery | |
CN114209433A (zh) | 一种手术机器人导航定位方法及装置 | |
KR102577964B1 (ko) | 간 수술을 위한 정렬 시스템 | |
Yang et al. | A novel craniotomy simulation system for evaluation of stereo-pair reconstruction fidelity and tracking | |
KR101544712B1 (ko) | 수술용 항법 시스템의 교정 방법 |
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 |