CN106023140B - 一种术前三维图像与术中三维图像的无标志物配准方法 - Google Patents

一种术前三维图像与术中三维图像的无标志物配准方法 Download PDF

Info

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
Application number
CN201610293632.0A
Other languages
English (en)
Other versions
CN106023140A (zh
Inventor
李寅岩
徐进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Tinavi Medical Technology Co Ltd
Original Assignee
Beijing Tinavi Medical Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Tinavi Medical Technology Co Ltd filed Critical Beijing Tinavi Medical Technology Co Ltd
Priority to CN201610293632.0A priority Critical patent/CN106023140B/zh
Publication of CN106023140A publication Critical patent/CN106023140A/zh
Application granted granted Critical
Publication of CN106023140B publication Critical patent/CN106023140B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

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矩阵的形式。
CN201610293632.0A 2016-05-05 2016-05-05 一种术前三维图像与术中三维图像的无标志物配准方法 Active CN106023140B (zh)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113573776A (zh) * 2020-02-14 2021-10-29 西安大医集团股份有限公司 图像引导方法、装置、放疗设备和计算机存储介质

Citations (2)

* Cited by examiner, † Cited by third party
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 中国科学院合肥物质科学研究院 一种基于带权模糊互信息的三维医学图像配准方法及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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) 一种术中导航系统中的非标记医学图像配准系统及方法
US20210346101A1 (en) Methods and devices for tracking objects by surgical navigation systems
EP2723268B1 (en) Ultrasound ct registration for positioning
CN104434313B (zh) 一种腹部外科手术导航方法及系统
EP2931159B1 (en) Registration and navigation using a three-dimensional tracking sensor
EP3116401B1 (en) Anatomical image projection system
EP3057064B1 (en) Pose determination from a pattern of four leds
CN102266250A (zh) 超声手术导航系统及超声手术导航方法
EP3212104B1 (en) Hybrid navigation system for surgical interventions
CN115038401A (zh) 用于外部固定器的方法和装置
CN106023140B (zh) 一种术前三维图像与术中三维图像的无标志物配准方法
CN203953815U (zh) 辅助外科手术的术中导航系统
US11826107B2 (en) Registration system for medical navigation and method of operation thereof
TW201519856A (zh) 手術定位方法
Gard et al. Image-based measurement by instrument tip tracking for tympanoplasty using digital surgical microscopy
US20240016550A1 (en) Scanner for intraoperative application
KR102577964B1 (ko) 간 수술을 위한 정렬 시스템
Yang et al. A novel craniotomy simulation system for evaluation of stereo-pair reconstruction fidelity and tracking
KR101544712B1 (ko) 수술용 항법 시스템의 교정 방법
Liu et al. A novel marker tracking method based on extended Kalman filter for multi-camera optical tracking systems
Khosravi et al. One-step needle pose estimation for ultrasound guided biopsies
CN104809717B (zh) 非刚体医学图像兴趣点标识、定位方法和装置、医疗设备
CN117357217A (zh) 一种经皮穿刺机械臂的系统标定方法
Maitland et al. A Video Based Tracker for use in Computer Aided Surgery.

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