CN108470357B - 耦合呼吸时相的弹性配准方法 - Google Patents
耦合呼吸时相的弹性配准方法 Download PDFInfo
- Publication number
- CN108470357B CN108470357B CN201810258592.5A CN201810258592A CN108470357B CN 108470357 B CN108470357 B CN 108470357B CN 201810258592 A CN201810258592 A CN 201810258592A CN 108470357 B CN108470357 B CN 108470357B
- Authority
- CN
- China
- Prior art keywords
- image
- registration
- respiratory
- ray
- images
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- 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/10076—4D tomography; Time-sequential 3D tomography
-
- 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/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- 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
- G06T2207/30061—Lung
-
- 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
- G06T2207/30096—Tumor; Lesion
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种耦合呼吸时相的弹性配准方法,将治疗过程中X射线图像与4D CT图像进行结合,基于X射线透视成像提取呼吸信号,采用4D CT不同时相图像对应X射线透视图像,将呼吸时相信息耦合到治疗前和治疗过程中的弹性配准当中,根据相似性测度对配准结果进行实时更新,通过金字塔多分辨率配准同时调节邻域灰度范围与配准图像的范围,使配准范围呈缩小态势,能够更快速和精确得到弹性配准结果,使弹性配准达到既快又优,进一步提高医学图像的配准质量。
Description
技术领域
本发明涉及肿瘤放射治疗和手术导航过程中的医学图像配准和摆位验证领域,特别是涉及一种耦合呼吸时相的弹性配准方法。
背景技术
对于胸腹部肿瘤如肺癌、肝癌、胰腺癌等,肿瘤运动会受到呼吸运动的影响而出现位置偏差,特别是在治疗过程中,肿瘤与周围正常组织器官之间的相对位置会发生一定程度的改变,同时还存在解剖结构形变等误差,严重影响了人体内剂量分布,降低了治疗效果,甚至会带来多种放疗并发症的产生。医学图像配准是肿瘤靶区勾画和放射治疗前摆位验证的核心。医学图像配准按照空间变换关系可分为刚性配准与弹性配准两种。刚性配准结果取得的是图像整体平移、旋转或缩放尺度,无法实现器官在配准过程中的局部形变。弹性配准将图像看作是一个弹性物体,在内外力的相互作用下使图像向外拉伸、向内压缩、弯曲或扭转,不断迭代驱动浮动图像发生局部形变最终与参考图像对齐,适用于处理器官在配准过程中的局部形变。
胸腹部器官在治疗过程中受到心脏跳动、肠胃蠕动、呼吸运动等因素影响会产生局部形变。为了降低病人在整个治疗过程中胸腹部呼吸运动及肿瘤形态变化对放射治疗疗效的影响,需要在整个治疗过程中采集病人不同时期或同一时期不同呼吸时相的图像进行精确快速的配准以便快速的确定下一步治疗方案,因此弹性配准的精确性对胸腹部肿瘤的治疗精度具有重要意义。传统CT图像配准大多考虑图像与图像之间的静态配准,而很少考虑到呼吸运动控制场合下图像与图像之间的动态配准,因此亟需提供一种新型的耦合呼吸时相的弹性配准方法来解决上述问题。
发明内容
本发明所要解决的技术问题是提供一种耦合呼吸时相的弹性配准方法,能够显著提高弹性配准的速度和精确性。
为解决上述技术问题,本发明采用的一个技术方案是:提供一种耦合呼吸时相的弹性配准方法,包括以下步骤:
(1)基于X射线透视成像提取呼吸信号,将X射线的图像信息及不同呼吸相位图像的时序信息进行对应;获取并分析4D CT每个时相下不同位置的CT片之间对应关系;
(2)利用步骤(1)中采集得到的X射线透视图像和4D CT中各个时相的CT片之间的对应关系,将放射治疗过程中X射线透视图像与4D CT图像进行结合,利用4D CT提供的呼吸时序信息建立与X射线透视图像中呼吸时序信号的对应规则;
(3)利用步骤(2)中得到的耦合呼吸时相信息的X射线图像和4D CT图像,用于放射治疗前和放射治疗过程中的弹性配准当中,根据相似性测度对配准结果进行实时更新,通过金字塔多分辨率配准同时调节邻域灰度范围与配准图像的范围,确定出肿瘤位置的摆位误差。
在本发明一个较佳实施例中,在步骤(2)中,4D CT提供的呼吸时序信息与X射线透视图像中呼吸时序信号的对应规则的建立过程为:
对已有的4D CT图像进行图像滤波、图像层间插值,提取图像物理信息和呼吸时相信息;采集全呼吸周期下X射线透视图像得到多个时相的2D X射线图像,基于二维图像中的标记点信息提取呼吸信号,进而通过2D和4D图像的配准,建立4D CT中的呼吸时序信息与X射线图像中的呼吸时序信号的对应规则。
在本发明一个较佳实施例中,在步骤(3)中,进行相似性测度时,当待配准图像的灰度均方误差或梯度均方误差达到最小值时,配准结束并反馈待配准图像对应的呼吸时相,再将待配准图像对应呼吸时相上的图像取出叠加到配准图像上完成配准。
本发明的有益效果是:本发明所述弹性配准方法同时关心X射线的图像信息及不同呼吸相位图像的时序信息与图像信息,将治疗过程中X射线图像与4D CT图像进行结合,基于X射线透视成像提取呼吸信号,采用4D CT不同时相图像对应X射线透视图像,将呼吸时相信息耦合到治疗前和治疗过程中的弹性配准当中,根据相似性测度对配准结果进行实时更新,通过金字塔多分辨率配准同时调节邻域灰度范围与配准图像的范围,使配准范围呈缩小态势,能够更快速和精确得到弹性配准结果,使弹性配准达到既快又优,进一步提高医学图像的配准质量。
附图说明
图1是本发明耦合呼吸时相的弹性配准方法一较佳实施例的流程图。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
请参阅图1,本发明实施例包括:
采用一组肺癌病人在自由深呼吸模式下的4D肺部CT图像作为测试用例,采集到的图像的体积为512×512×640,分辨率为0.97×0.97×0.5mm3。首先采集4D肺部CT图像,同时记录体外标记的三维空间坐标;接着采用距离对应关联法建立4D肺部CT图像与呼吸信号的关系,分析4D CT每个时相下不同位置的CT片之间对应关系,通过标记物相同时间段内在4D肺部CT图像中和呼吸信号中的位移量的比较,确定4D肺部CT图像和呼吸信号的对应关系。
基于X射线透视成像提取呼吸信号,将X射线的图像信息及不同呼吸相位图像的时序信息进行对应;利用基于空间信息与灰度信息的图像弹性配准方法配准4D肺部CT图像与X射线透视图像,获得形变场T{T|i1,2,...,m},第i时刻配准控制点的形变位移设为D。将位于同一呼吸周期中的4D肺部CT图像与X射线图像进行匹配,在预先构建的序列影像与呼吸信号对应关系的基础上,利用4D肺部CT弹性配准获取每一个配准控制点对于时间的三维形变位移,基于4D CT图像生成一个周期局部控制点的邻域灰度值或邻域梯度值作为模板,建立邻域模板相似性测度函数,将待配准的耦合呼吸时相信息的X射线图像对应局部控制点邻域与4D CT图像的全周期邻域进行相似性匹配,引入邻域模板作为约束,使配准既快又优。当待配准图像的灰度均方误差或梯度均方误差达到最小值时,配准结束并反馈待配准图像对应的呼吸时相,再将待配准图像对应呼吸时相上的图像取出叠加到配准图像上完成配准。
由于参与相似性匹配的计算量较大,使用多分辨率金字塔策略作用于邻域模板,动态调整邻域模板尺寸大小,通过由粗到精的匹配,加速弹性配准的进程。
根据相似性测度对配准结果进行实时更新,通过金字塔多分辨率配准同时调节邻域灰度范围与配准图像的范围,最终确定出肿瘤位置的摆位误差。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种耦合呼吸时相的弹性配准方法,包括以下步骤:
(1)基于X射线透视成像提取呼吸信号,将X射线的图像信息及不同呼吸相位图像的时序信息进行对应,获取并分析4D CT中每个时相下不同位置的CT片之间对应关系得到呼吸时序信息;
(2)利用步骤(1)中采集得到的X射线透视图像和4D CT中各个时相的CT片之间的对应关系,将X射线透视图像与4D CT图像进行结合,利用4D CT提供的呼吸时序信息建立与X射线透视图像中呼吸时序信号的对应规则;
(3)利用步骤(2)中得到的耦合呼吸时相信息的X射线图像和4D CT图像进行弹性配准,根据相似性测度对配准结果进行实时更新,通过金字塔多分辨率配准同时调节邻域灰度范围与配准图像的范围,确定出肿瘤位置的摆位误差。
2.根据权利要求1所述的耦合呼吸时相的弹性配准方法,其特征在于,在步骤(2)中,4DCT提供的呼吸时序信息与X射线透视图像中呼吸时序信号的对应规则的建立过程为:
对已有的4D CT图像进行图像滤波、图像层间插值,提取图像物理信息和呼吸时相信息;采集全呼吸周期下X射线透视图像得到多个时相的2D X射线图像,基于二维图像中的标记点信息提取呼吸信号,进而通过2D和4D图像的配准,建立4D CT中的呼吸时序信息与X射线图像中的呼吸时序信号的对应规则。
3.根据权利要求1所述的耦合呼吸时相的弹性配准方法,其特征在于,在步骤(3)中,进行相似性测度时,当待配准图像的灰度均方误差或梯度均方误差达到最小值时,配准结束并反馈待配准图像对应的呼吸时相,再将待配准图像对应呼吸时相上的图像取出叠加到配准图像上完成配准。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810258592.5A CN108470357B (zh) | 2018-03-27 | 2018-03-27 | 耦合呼吸时相的弹性配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810258592.5A CN108470357B (zh) | 2018-03-27 | 2018-03-27 | 耦合呼吸时相的弹性配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108470357A CN108470357A (zh) | 2018-08-31 |
CN108470357B true CN108470357B (zh) | 2022-04-01 |
Family
ID=63264837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810258592.5A Active CN108470357B (zh) | 2018-03-27 | 2018-03-27 | 耦合呼吸时相的弹性配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108470357B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113592804B (zh) * | 2021-07-27 | 2024-05-07 | 东软医疗系统股份有限公司 | 图像处理方法、装置及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1672651A (zh) * | 2004-02-17 | 2005-09-28 | 西门子共同研究公司 | 在有呼吸运动的情况下用于放射治疗的患者定位系统和方法 |
CN101773395A (zh) * | 2009-12-31 | 2010-07-14 | 华中科技大学 | 一种从单臂x射线造影图中提取呼吸运动参数的方法 |
CN104268914A (zh) * | 2014-10-24 | 2015-01-07 | 山东师范大学 | 一种4d-ct不同时相序列图像的重建方法 |
CN104268895A (zh) * | 2014-10-24 | 2015-01-07 | 山东师范大学 | 一种联合空域和时域信息的4d-ct形变配准方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012042449A2 (en) * | 2010-09-30 | 2012-04-05 | Koninklijke Philips Electronics N.V. | Image and annotation display |
-
2018
- 2018-03-27 CN CN201810258592.5A patent/CN108470357B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1672651A (zh) * | 2004-02-17 | 2005-09-28 | 西门子共同研究公司 | 在有呼吸运动的情况下用于放射治疗的患者定位系统和方法 |
CN101773395A (zh) * | 2009-12-31 | 2010-07-14 | 华中科技大学 | 一种从单臂x射线造影图中提取呼吸运动参数的方法 |
CN104268914A (zh) * | 2014-10-24 | 2015-01-07 | 山东师范大学 | 一种4d-ct不同时相序列图像的重建方法 |
CN104268895A (zh) * | 2014-10-24 | 2015-01-07 | 山东师范大学 | 一种联合空域和时域信息的4d-ct形变配准方法 |
Non-Patent Citations (3)
Title |
---|
4D-CT imaging of a volume influenced by respiratory motion on multi-slice CT;Tinsu Pan等;《Med. Phys.》;20040229;全文 * |
基于4DCT的腹部器官呼吸运动分析;习勉等;《癌症》;20091231;全文 * |
基于透视成像的呼吸与肺部组织运动的关系研究;杨词慧等;《南昌航空大学学报:自然科学版》;20150930;第29卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108470357A (zh) | 2018-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10293182B2 (en) | Systems and methods for adaptive replanning based on multi-modality imaging | |
US20220047247A1 (en) | Apparatus and method for real-time tracking of tissue structures | |
Blackall et al. | MRI-based measurements of respiratory motion variability and assessment of imaging strategies for radiotherapy planning | |
McClelland et al. | Respiratory motion models: a review | |
Blackall et al. | Alignment of sparse freehand 3-D ultrasound with preoperative images of the liver using models of respiratory motion and deformation | |
EP2175931B1 (en) | Systems for compensating for changes in anatomy of radiotherapy patients | |
US7117026B2 (en) | Physiological model based non-rigid image registration | |
Ébastien Clippe et al. | Patient setup error measurement using 3D intensity-based image registration techniques | |
US20090306497A1 (en) | Calibration method for catheter tracking system using medical imaging data | |
US20080107312A1 (en) | Device and Process for Multimodal Registration of Images | |
AU2015360825A1 (en) | Motion management in MRI-guided Linac | |
Fu et al. | Xsight lung tracking system: a fiducial-less method for respiratory motion tracking | |
Penney et al. | Overview of an ultrasound to CT or MR registration system for use in thermal ablation of liver metastases | |
Alam et al. | Medical image registration: Classification, applications and issues | |
CN108470357B (zh) | 耦合呼吸时相的弹性配准方法 | |
KR20170065118A (ko) | 치료계획용 4d mdct 영상과 치료시 획득한 4d cbct 영상 간 영상 정합 및 종양 매칭을 이용한 방사선 치료시 종양 움직임 추적 방법 및 장치 | |
Preiswerk et al. | Robust tumour tracking from 2D imaging using a population-based statistical motion model | |
JP6692817B2 (ja) | 対象物体の変位を計算する方法及びシステム | |
Timinger et al. | Motion compensated coronary interventional navigation by means of diaphragm tracking and elastic motion models | |
Krueger et al. | Modality-integrated magnetic catheter tracking for x-ray vascular interventions | |
Wein et al. | Ultrasound based respiratory motion compensation in the abdomen | |
CN116128838A (zh) | 基于双目视觉的呼吸检测系统、4d-ct图像重建系统和方法 | |
KR102383753B1 (ko) | 의료영상 정합 방법 및 그 장치 | |
Hawkes et al. | Measuring and modeling soft tissue deformation for image guided interventions | |
Ranjbar | Simulating the breathing of lung cancer patients to estimate tumor motion and deformation at the time of radiation treatment |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: Room 406, block B, phase I, Zhongdan Ecological Life Science Industrial Park, No. 3-1, xinjinhu Road, Jiangbei new district, Nanjing City, Jiangsu Province Applicant after: Zhongke Chaojing (Nanjing) Technology Co., Ltd Address before: 230601 NO.560, Haiheng building, No.6 Cuiwei Road, Hefei Economic and Technological Development Zone, Anhui Province Applicant before: ZHONGKE CHAOJING (ANHUI) TECHNOLOGY Co.,Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |