CN105976372A - 一种术前三维影像与术中透视图像的无标定物配准方法 - Google Patents

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

Info

Publication number
CN105976372A
CN105976372A CN201610293634.XA CN201610293634A CN105976372A CN 105976372 A CN105976372 A CN 105976372A CN 201610293634 A CN201610293634 A CN 201610293634A CN 105976372 A CN105976372 A CN 105976372A
Authority
CN
China
Prior art keywords
image
images
drr
ray
dimensional
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
Application number
CN201610293634.XA
Other languages
English (en)
Other versions
CN105976372B (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 CN201610293634.XA priority Critical patent/CN105976372B/zh
Publication of CN105976372A publication Critical patent/CN105976372A/zh
Application granted granted Critical
Publication of CN105976372B publication Critical patent/CN105976372B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • G06T2207/10124Digitally reconstructed radiograph [DRR]

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明涉及一种术前三维影像与术中透视图像的无标定物配准方法,它包括以下步骤:1)手术开始前,CT或MRI获取一系列三维图像,根据手术需求生成N个DRR图像,得到每个图像的位姿坐标W;2)手术中,C型臂上的X射线机获取不同角度的X‑Ray透视图像;3)对DRR图像i和X‑Ray透视图像均做归一化处理;4)采用非线性直方图匹配算法;5)对DRR图像i与X‑Ray透射图像进行秩相关相似度比较,如果两组图像匹配,获取X‑Ray透射图像位姿参数,进入下一步;否则提取DRR图像i+1信息,返回步骤3);6)基于多分辨率机制和梯度差异相关性为术前三维影像和X‑Ray透视图像精确配准,输出精确配准后的图像。本发明降低了病人手术中所受的伤害。

Description

一种术前三维影像与术中透视图像的无标定物配准方法
技术领域
本发明涉及一种无标定物配准方法,特别是关于一种医学领域中使用的术前三维影像与术中透视图像的无标定物配准方法。
背景技术
目前,在常见医疗机器人的手术导航定位过程中,尤其在脊柱手术中,传统方法是在脊柱脊突上植入标定物,以建立机器人坐标系与手术空间坐标系之间的映射关系,并在此基础上进行手术路径规划和定位。该方法能够减少X射线透视次数,从而降低对病人及医生的危害;但术前准备繁琐,且需要植入标定物,病人所受损伤较大。
一种方法是利用2D/3D图像配准算法的图像导航的放射手术方法和系统,其利用了分层和迭代的2D/3D配准算法求出内层面和外层面的变换参数。此配准算法的缺点是准确度较低,速度较慢。另一种方法是使用DRR技术对MRI和X-Ray图像配准的算法和系统,其特点在于先对MRI图像做提取分割,由分割后的图像生成DRR用于配准,缺点是算法需要事先训练,从而获取特征用于分割,受限于训练的集合大小及范围,而且成本较大。
发明内容
针对上述问题,本发明的目的是提供一种术前三维影像与术中透视图像的无标定物配准方法,其操作简单,速度快,准确度高。
为实现上述目的,本发明采取以下技术方案:一种术前三维影像与术中透视图像的无标定物配准方法,其特征在于它包括以下步骤:1)在手术开始前,由CT或MRI获取一系列三维图像,根据手术需求,在冠状位,矢状位或横断位的视角方向上,从0°开始每间隔预设度数根据CT或MRI提供的三维图像生成N个DRR图像,获取对应的位姿参数;2)手术中,C型臂上的X射线机获取不同角度的X-Ray透视图像;3)对DRR图像i和X-Ray透视图像上的每一个像素均做二维高斯加权归一化处理,每一个像素原始灰度值减去高斯加权均值,然后除以高斯加权均方差从而获取归一化后的像素灰度值;4)采用非线性直方图匹配算法,校正DRR图像i与X-Ray透射图像的灰度差异,选择灰度值大于投影图像灰度均值的区域进行直方图匹配,消除背景像素对直方图的影响;5)对DRR图像i与X-Ray透射图像进行秩相关相似度比较,如果两组图像匹配,获取X-Ray透射图像位姿参数,并进入下一步;否则提取DRR图像i+1信息,并返回步骤3);6)基于多分辨率机制和梯度差异相关性为术前三维影像和X-Ray透视图像精确配准,输出精确配准后的图像,完成术前三维影像与术中透视图像的无标定物配准。
优选地,所述步骤1)中,从0°开始每间隔预设度数优选为1°或2°根据CT/MRI提供的三维图像生成N个DRR图像。
优选地,所述步骤1)中,所述DRR图像的生成采用基于GPU加速的光影透射算法中的cuberille算法,同时使用稀疏抽样提高图像性能。
优选地,所述步骤2)中,C型臂上的X射线机获取X-Ray透视图像至少获取两张。
优选地,所述步骤5)中,所述秩相关作为相似度测量的配准算法,秩相关MSRC函数如下:
其中,表示每个像素的秩方差,是随机选出的像素数量,其数值远小于图像所有像素数量。
本发明由于采取以上技术方案,其具有以下优点:1、本发明采用了基于GPU的DRR稀疏抽样生成算法,保证了初始参数的准确度和快速获取。2、本发明采用了术前三维影像与术中透视图像的无标定物配准方法,降低了病人所受的损伤。3、本发明采用了随机秩相关作为图像相似度的衡量标准,提高了结果的可靠性和准确性。
附图说明
图1是本发明从不同视野角度对术前CT或MRI图像进行数字投影重建示意图;
图2是本发明X射线源围绕病人旋转移动获取N个投影图像示意图;
图3是本发明整体流程示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1~图3所示,本发明提供一种术前三维影像与术中透视图像的无标定物配准方法,其包括以下步骤:
1)在手术开始前,由CT或MRI获取一系列三维图像,根据手术需求,在冠状位、矢状位或横断位的视角方向上,从0°开始每间隔预设度数根据CT或MRI提供的三维图像生成N个放射成像图像(DRR图像),获取对应的位姿参数W;其中,预设度数优选为1°或2°;
W=(θxyz,X,Y,Z) (1)
式中,θx、θy、θz表示旋转方向,X、Y、Z表示在坐标系各个方向下的平移量。
2)手术中,C型臂上的X射线机获取不同角度的X-Ray透视图像;其中,透视图像至少获取两张;
3)对DRR图像i(0<i<N,且i为整数)和X-Ray透视图像上的每一个像素均做二维高斯加权归一化处理,每一个像素原始灰度值减去高斯加权均值,然后除以高斯加权均方差从而获取归一化后的像素灰度值;
4)采用非线性直方图匹配算法,校正DRR图像i与X-Ray透射图像的灰度差异,选择灰度值大于投影图像灰度均值的区域进行直方图匹配,消除背景像素对直方图的影响;
5)对DRR图像i与X-Ray透射图像进行秩相关相似度比较,如果两组图像匹配,获取X-Ray透射图像位姿参数,并进入下一步;否则提取DRR图像i+1信息,并返回步骤3);
6)基于多分辨率机制和梯度差异相关性为术前三维影像和X-Ray透视图像精确配准,输出精确配准后的图像,完成术前三维影像与术中透视图像的无标定物配准。
上述步骤1)中,DRR图像的生成采用基于GPU加速的光影透射算法中的cuberille算法,使图像的边界尽可能的小,同时使用稀疏抽样提高图像性能。
上述步骤5)中,秩相关作为相似度测量的配准算法,秩相关MSRC函数如下:
其中,表示每个像素的秩方差,是随机选出的像素数量,其数值远小于图像所有像素数量。
上述步骤6)中,多分辨率机制用于加速配准进程,同时也可以在能量函数最小化过程中避免陷入局部最小。
上述各实施例仅用于说明本发明,各部件的结构、尺寸、设置位置及形状都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别部件进行的改进和等同变换,均不应排除在本发明的保护范围之外。

Claims (5)

1.一种术前三维影像与术中透视图像的无标定物配准方法,其特征在于:它包括以下步骤:
1)在手术开始前,由CT或MRI获取一系列三维图像,根据手术需求,在冠状位,矢状位或横断位的视角方向上,从0°开始每间隔预设度数根据CT或MRI提供的三维图像生成N个DRR图像,获取对应的位姿参数;
2)手术中,C型臂上的X射线机获取不同角度的X-Ray透视图像;
3)对DRR图像i和X-Ray透视图像上的每一个像素均做二维高斯加权归一化处理,每一个像素原始灰度值减去高斯加权均值,然后除以高斯加权均方差从而获取归一化后的像素灰度值;
4)采用非线性直方图匹配算法,校正DRR图像i与X-Ray透射图像的灰度差异,选择灰度值大于投影图像灰度均值的区域进行直方图匹配,消除背景像素对直方图的影响;
5)对DRR图像i与X-Ray透射图像进行秩相关相似度比较,如果两组图像匹配,获取X-Ray透射图像位姿参数,并进入下一步;否则提取DRR图像i+1信息,并返回步骤3);
6)基于多分辨率机制和梯度差异相关性为术前三维影像和X-Ray透视图像精确配准,输出精确配准后的图像,完成术前三维影像与术中透视图像的无标定物配准。
2.如权利要求1所述的一种术前三维影像与术中透视图像的无标定物配准方法,其特征在于:所述步骤1)中,从0°开始每间隔预设度数优选为1°或2°根据CT/MRI提供的三维图像生成N个DRR图像。
3.如权利要求1所述的一种术前三维影像与术中透视图像的无标定物配准方法,其特征在于:所述步骤1)中,所述DRR图像的生成采用基于GPU加速的光影透射算法中的cuberille算法,同时使用稀疏抽样提高图像性能。
4.如权利要求1所述的一种术前三维影像与术中透视图像的无标定物配准方法,其特征在于:所述步骤2)中,C型臂上的X射线机获取X-Ray透视图像至少获取两张。
5.如权利要求1所述的一种术前三维影像与术中透视图像的无标定物配准方法,其特征在于:所述步骤5)中,所述秩相关作为相似度测量的配准算法,秩相关MSRC函数如下:
其中,表示每个像素的秩方差,是随机选出的像素数量,其数值远小于图像所有像素数量。
CN201610293634.XA 2016-05-05 2016-05-05 一种术前三维影像与术中透视图像的无标定物配准方法 Active CN105976372B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610293634.XA CN105976372B (zh) 2016-05-05 2016-05-05 一种术前三维影像与术中透视图像的无标定物配准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610293634.XA CN105976372B (zh) 2016-05-05 2016-05-05 一种术前三维影像与术中透视图像的无标定物配准方法

Publications (2)

Publication Number Publication Date
CN105976372A true CN105976372A (zh) 2016-09-28
CN105976372B CN105976372B (zh) 2019-03-05

Family

ID=56991327

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610293634.XA Active CN105976372B (zh) 2016-05-05 2016-05-05 一种术前三维影像与术中透视图像的无标定物配准方法

Country Status (1)

Country Link
CN (1) CN105976372B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109700528A (zh) * 2019-02-27 2019-05-03 江苏霆升科技有限公司 一种实时构建心脏三维模型方法及装置
CN109978927A (zh) * 2019-03-12 2019-07-05 上海嘉奥信息科技发展有限公司 图像配准误差的测量装置及测量方法和系统
CN109993792A (zh) * 2019-03-04 2019-07-09 杭州三坛医疗科技有限公司 投影方法、装置及系统和可读存储介质
CN110148160A (zh) * 2019-05-22 2019-08-20 合肥中科离子医学技术装备有限公司 一种正交x射线影像快速2d-3d医学图像配准方法
CN111723836A (zh) * 2019-03-21 2020-09-29 杭州三坛医疗科技有限公司 图像相似度的计算方法及装置、电子设备、存储介质
CN112472293A (zh) * 2020-12-15 2021-03-12 山东威高医疗科技有限公司 一种术前三维影像与术中透视图像的配准方法
CN112509022A (zh) * 2020-12-17 2021-03-16 安徽埃克索医疗机器人有限公司 一种术前三维影像与术中透视图像的无标定物配准方法
CN112837358A (zh) * 2021-03-25 2021-05-25 中科超精(南京)科技有限公司 一种多透视图耦合的台面物体定位方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102147919A (zh) * 2010-02-10 2011-08-10 昆明医学院第一附属医院 一种校正术前三维图像的术中配准方法和装置
US8036441B2 (en) * 2003-04-25 2011-10-11 Medtronic Navigation, Inc. Method and apparatus for performing 2D to 3D registration

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8036441B2 (en) * 2003-04-25 2011-10-11 Medtronic Navigation, Inc. Method and apparatus for performing 2D to 3D registration
CN102147919A (zh) * 2010-02-10 2011-08-10 昆明医学院第一附属医院 一种校正术前三维图像的术中配准方法和装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHEN-RUI CHOU等: "2D/3D Image Registration using Regression Learning", 《COMPUTER VISION & IMAGE UNDERSTANDING》 *
WOLFGANG BIRKFELLNER等: "Stochastic rank correlation: A robust merit function for 2D/3D registration of image data obtained at different energies", 《MEDICAL PHYSICS》 *
张薇等: "基于灰度的二维/三维图像配准方法及其在骨科导航手术中的实现", 《中国医学影像技术》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109700528A (zh) * 2019-02-27 2019-05-03 江苏霆升科技有限公司 一种实时构建心脏三维模型方法及装置
CN109700528B (zh) * 2019-02-27 2021-02-02 江苏霆升科技有限公司 一种实时构建心脏三维模型方法及装置
CN109993792A (zh) * 2019-03-04 2019-07-09 杭州三坛医疗科技有限公司 投影方法、装置及系统和可读存储介质
CN109978927A (zh) * 2019-03-12 2019-07-05 上海嘉奥信息科技发展有限公司 图像配准误差的测量装置及测量方法和系统
CN111723836A (zh) * 2019-03-21 2020-09-29 杭州三坛医疗科技有限公司 图像相似度的计算方法及装置、电子设备、存储介质
CN110148160A (zh) * 2019-05-22 2019-08-20 合肥中科离子医学技术装备有限公司 一种正交x射线影像快速2d-3d医学图像配准方法
CN112472293A (zh) * 2020-12-15 2021-03-12 山东威高医疗科技有限公司 一种术前三维影像与术中透视图像的配准方法
CN112472293B (zh) * 2020-12-15 2022-10-21 山东威高医疗科技有限公司 一种术前三维影像与术中透视图像的配准方法
CN112509022A (zh) * 2020-12-17 2021-03-16 安徽埃克索医疗机器人有限公司 一种术前三维影像与术中透视图像的无标定物配准方法
CN112837358A (zh) * 2021-03-25 2021-05-25 中科超精(南京)科技有限公司 一种多透视图耦合的台面物体定位方法

Also Published As

Publication number Publication date
CN105976372B (zh) 2019-03-05

Similar Documents

Publication Publication Date Title
CN105976372B (zh) 一种术前三维影像与术中透视图像的无标定物配准方法
JP5906015B2 (ja) 特徴に基づいた2次元/3次元画像のレジストレーション
US9240046B2 (en) Method and system to assist 2D-3D image registration
JP2019526124A (ja) 三次元表面の画像を再構築するための方法、装置およびシステム
US10154820B2 (en) Method and apparatus for navigating CT scan with a marker
EP2849630B1 (en) Virtual fiducial markers
Uneri et al. 3D–2D registration for surgical guidance: effect of projection view angles on registration accuracy
CN109152566B (zh) 校正超声融合成像系统中的探头引起的变形
JP2010246883A (ja) 患者位置決めシステム
JP6349278B2 (ja) 放射線撮像装置、画像処理方法及びプログラム
US20210174502A1 (en) Method for metal artifact avoidance in x-ray imaging
TWI836493B (zh) 註冊二維影像資料組與感興趣部位的三維影像資料組的方法及導航系統
CN111093505B (zh) 放射线拍摄装置以及图像处理方法
CN111757701B (zh) 用于检测解剖区域的系统和方法
WO2009090572A2 (en) Image registration alignment metric
CN112509022A (zh) 一种术前三维影像与术中透视图像的无标定物配准方法
TW201821030A (zh) 提供光學對準特徵的牙齒影像收集裝置以及相關系統及方法
US9254106B2 (en) Method for completing a medical image data set
JP6145870B2 (ja) 画像表示装置および方法、並びにプログラム
EP4301229B1 (en) Image-based planning of tomographic scan
JP5016231B2 (ja) 撮像の幾何学的パラメータを決定する方法及び装置
US8787649B2 (en) System and method for processing an X-ray image of an organ
KR20160031794A (ko) 병변 검출 장치 및 그 방법
Zheng et al. Point similarity measures based on MRF modeling of difference images for spline-based 2D-3D rigid registration of X-ray fluoroscopy to CT images
CN114066947B (zh) 图像配准方法和图像配准装置

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