CN111667515A - 基于流形正则化的血管3d/2d弹性配准方法及装置 - Google Patents
基于流形正则化的血管3d/2d弹性配准方法及装置 Download PDFInfo
- Publication number
- CN111667515A CN111667515A CN202010325644.3A CN202010325644A CN111667515A CN 111667515 A CN111667515 A CN 111667515A CN 202010325644 A CN202010325644 A CN 202010325644A CN 111667515 A CN111667515 A CN 111667515A
- Authority
- CN
- China
- Prior art keywords
- points
- vessel
- matrix
- elastic registration
- manifold
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 47
- 238000006073 displacement reaction Methods 0.000 claims abstract description 25
- 230000009466 transformation Effects 0.000 claims abstract description 24
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 239000008280 blood Substances 0.000 claims abstract description 6
- 210000004369 blood Anatomy 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 59
- 230000005489 elastic deformation Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000844 transformation Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 2
- 230000002792 vascular Effects 0.000 description 4
- 238000010968 computed tomography angiography Methods 0.000 description 3
- 210000003484 anatomy Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000002583 angiography Methods 0.000 description 1
- 210000000702 aorta abdominal Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000002586 coronary angiography Methods 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
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
-
- 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/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- 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
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
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
Description
技术领域
本发明涉及医学图像处理的技术领域,尤其涉及一种基于流形正则化的血管3D/2D弹性配准方法,以及基于流形正则化的血管3D/2D弹性配准装置。
背景技术
血管3D/2D配准目的是,在术前的3D图像和术中的2D图像中找到对应血管解剖结构的对应关系,并估计3D血管的刚性和弹性变换,使得3D和2D血管在对应结构上达到空间位置的一致。3D和2D血管分别从CTA(计算机断层扫描血管造影)和XRA(X射线血管造影)图像中提取,经3D/2D配准后,3D和2D血管在对应解剖结构上达到空间位置的一致。根据3D/2D配准任务的不同,对其进行任务分解,包括:初始化、刚性变换和弹性变换。
弹性变换中的形变模型可以分为两类:一类是基于大量训练数据统计的形变模型,如统计形状模型(SSM)、主动表观模型(AAM)和可形变模型(Morphable model);另一类是基于平滑先验假设的形变模型,如薄板样条模型(TPS)、B样条模型和相干点漂移模型(CPD)。
针对血管配准问题,许多研究者在以上这些模型的基础上做出了适用于血管模型的改进。Groher等人于2009年提出了一种基于单视角的血管形变模型,在构建的能量函数中对血管长度变化和非自然弯曲进行惩罚,通过求解目标函数对3D血管点位移的偏导构建目标函数梯度场,利用BFGS优化方法求得最优的弹性变换参数。2010年,Groher等人改进了方法,将配准问题表示为包含以图像信息响应为主的外能和以3D血管长度变化约束为主的内能的能量最小化问题,该方法直接利用2D造影图像的管状结构增强,无需对2D图像进行血管分割。Liao等人将所有血管点视为自由点,在腹主动脉血管的图模型的基础上,构建以距离测度、长度保持项和平滑正则项为主的相似性测度,使用多尺度优化策略求取血管点的位移值。Rivest-Henault等人使用了与类似的内能和外能构建的弹性配准能量函数,不同的是作者利用了血管分叉点处的形变约束来构建外能,并在双视图的冠脉造影图像上求解冠脉3D中心线的位移值。Liu等人利用2D TPS函数实现投影后的3D血管点到2D血管中心线的非刚性配准。Serradell等人针对心脏运动的先验知识,提出了一种合成的冠脉血管树形变模型。利用该先验随机的生成大量的血管中心线训练集,并使用SSM来构建血管形变模型。Baka等人则在真实的临床冠脉4D CTA数据上,提取出3D+t冠脉中心线,并构建真实冠脉的SSM。但是,这些改进的血管模型均只适用于特定的具有拓扑结构的点云结构。
发明内容
为克服现有技术的缺陷,本发明要解决的技术问题是提供了一种基于流形正则化的血管3D/2D弹性配准方法,其能够适用于任意具有拓扑结构的点云结构,可以得到最优的血管3D/2D弹性配准结果。
本发明的技术方案是:这种基于流形正则化的血管3D/2D弹性配准方法,其包括以下步骤:
(1)在已知血管3D和2D点的匹配关系基础上,表示3D点,表示2D点,i和j分别为点集的索引值;用 表示3D点集矩阵,其中为3D点坐标,用表示2D点集矩阵,其中 为2D点坐标,两个点集中对应点具有相同的索引值,K为点数量;
(4)获得血管3D/2D弹性配准的目标函数;
(5)基于梯度优化的算法对目标函数进行求解。
本发明提出一种基于流形正则化的血管弹性形变模型,该形变模型利用了血管点的位移在其血管拓扑上具有平滑约束的先验,因此该模型适用于任意具有拓扑结构的点云结构;采用基于梯度优化的算法对目标函数进行求解,可以得到最优的血管3D/2D弹性配准结果。
还提供了基于流形正则化的血管3D/2D弹性配准装置,其包括:
目标函数获得模块,其获得血管3D/2D弹性配准的目标函数;
求解模块,其基于梯度优化的算法对目标函数进行求解。
附图说明
图1示出了根据本发明的基于流形正则化的血管3D/2D弹性配准方法的流程图。
具体实施方式
如图1所示,这种基于流形正则化的血管3D/2D弹性配准方法,其包括以下步骤:
(1)在已知血管3D和2D点的匹配关系基础上,表示3D点,表示2D点,i和j分别为点集的索引值;用 表示3D点集矩阵,其中为3D点坐标,用表示2D点集矩阵,其中 为2D点坐标,两个点集中对应点具有相同的索引值,K为点数量;
(4)获得血管3D/2D弹性配准的目标函数;
(5)基于梯度优化的算法对目标函数进行求解。
本发明提出一种基于流形正则化的血管弹性形变模型,该形变模型利用了血管点的位移在其血管拓扑上具有平滑约束的先验,因此该模型适用于任意具有拓扑结构的点云结构;采用基于梯度优化的算法对目标函数进行求解,可以得到最优的血管3D/2D弹性配准结果。
优选地,所述步骤(2)中,采用高斯函数φ作为核函数,空间变换表示为公式(3)
优选地,所述步骤(4)中,血管3D/2D弹性配准的目标函数为公式(6)
优选地,所述步骤(5)中,目标函数的第二项Q2=λ1tr(WTUW)的梯度通过公式(8)计算
优选地,所述步骤(5)中,对于目标函数第三项Q3=λ2tr(τTLτ)的梯度通过公式(9)计算,τ=ΦTW-Y,
本领域普通技术人员可以理解,实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,包括上述实施例方法的各步骤,而所述的存储介质可以是:ROM/RAM、磁碟、光盘、存储卡等。因此,与本发明的方法相对应的,本发明还同时包括一种基于流形正则化的血管3D/2D弹性配准装置,该装置通常以与方法各步骤相对应的功能模块的形式表示。该装置包括:
目标函数获得模块,其获得血管3D/2D弹性配准的目标函数;
求解模块,其基于梯度优化的算法对目标函数进行求解。
以下更详细地说明本发明的内容。
在已知血管3D和2D点的匹配关系基础上,用 表示3D点集矩阵,其中为3D点坐标,用 表示2D点集矩阵,其中为2D点坐标,两个点集中对应点具有相同的索引值,K为点数量。弹性配准的目的是估计的空间变换使得3D点的变换可表示为τ(yk)为yk的位移。该问题是不适定的(ill-posed),因其具有无限数量的解。解决该问题的常用方法是使用在可再生核希尔伯特空间(Reproducing Kernel Hilbert Space,RKHS)中运行的正则化技术。令为一个RKHS,通过Tikhonov正则化来最小化目标函数可求得定义在的空间变换
其中k为正定的核函数 为任意非空集合。k满足可再生性,即对于任意与核函数k的內积等于其自身<f,k(ci,·)>k=f(ci,·),为核函数的输入参数。在此定义下RKHS的范数可用其內积表示 在本方法中,采用高斯函数φ作为核函数。因此空间变换可以表示为下式
可写成公式(3)
空间变换的范数项约束了由控制点C和系数矩阵W描述的空间形变场的平滑程度。而对于3D/2D配准,本方法对血管点的位移在其血管拓扑上进行约束,令矩阵A∈{0,1}K×K表示3D血管点拓扑的邻接矩阵,当yi和yj之间存在直接连接时aij=1,否则aij=0。
拓扑流形的约束可进一步限制空间变换在血管上的平滑性。通过观察对比,可发现拓扑流形约束同样可限制血管的长度变化。根据图拉普拉斯的定义,可用流形正则化表示,令D表示邻接矩阵A的度矩阵 则血管拓扑的拉普拉斯矩阵为L=D-A。令 表示3D点的位移矩阵,因此血管3D/2D弹性配准的目标函数可如下所示,
W为弹性形变模型的参数,λ1和λ2是权重系数。
目标函数的第二项Q2=λ1tr(WTUW)的梯度可直接通过矩阵偏导公式计算
对于目标函数第三项Q3=λ2tr(τTLτ),τ=ΦTW-Y,
以上所述,仅是本发明的较佳实施例,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属本发明技术方案的保护范围。
Claims (10)
2.根据权利要求1所述的基于流形正则化的血管3D/2D弹性配准方法,其特征在于:所述步骤(2)中,使用在可再生核希尔伯特空间RKHS中运行的正则化技术,令为一个RKHS,通过Tikhonov正则化来最小化目标函数,利用公式(1)求得定义在的空间变换
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010325644.3A CN111667515B (zh) | 2020-04-22 | 2020-04-22 | 基于流形正则化的血管3d/2d弹性配准方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010325644.3A CN111667515B (zh) | 2020-04-22 | 2020-04-22 | 基于流形正则化的血管3d/2d弹性配准方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111667515A true CN111667515A (zh) | 2020-09-15 |
CN111667515B CN111667515B (zh) | 2023-11-14 |
Family
ID=72382828
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010325644.3A Active CN111667515B (zh) | 2020-04-22 | 2020-04-22 | 基于流形正则化的血管3d/2d弹性配准方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111667515B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112465813A (zh) * | 2020-12-17 | 2021-03-09 | 北京工业大学 | 一种基于应力应变的血管内超声弹性分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170243361A1 (en) * | 2016-02-18 | 2017-08-24 | Roman Schaffert | 2D/3D Registration |
CN107392947A (zh) * | 2017-06-28 | 2017-11-24 | 西安电子科技大学 | 基于轮廓共面四点集的2d‑3d图像配准方法 |
CN109448031A (zh) * | 2018-09-28 | 2019-03-08 | 武汉大学 | 基于高斯场约束和流形正则化的图像配准方法及*** |
CN109993730A (zh) * | 2019-03-20 | 2019-07-09 | 北京理工大学 | 3d/2d血管配准方法及装置 |
-
2020
- 2020-04-22 CN CN202010325644.3A patent/CN111667515B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170243361A1 (en) * | 2016-02-18 | 2017-08-24 | Roman Schaffert | 2D/3D Registration |
CN107392947A (zh) * | 2017-06-28 | 2017-11-24 | 西安电子科技大学 | 基于轮廓共面四点集的2d‑3d图像配准方法 |
CN109448031A (zh) * | 2018-09-28 | 2019-03-08 | 武汉大学 | 基于高斯场约束和流形正则化的图像配准方法及*** |
CN109993730A (zh) * | 2019-03-20 | 2019-07-09 | 北京理工大学 | 3d/2d血管配准方法及装置 |
Non-Patent Citations (3)
Title |
---|
JIAYI MA 等: "Non-Rigid Point Set Registration with Robust Transformation Estimation under Manifold Regularization", 《THIRTY-FIRST AAAI CONFERENCE ON ARTI CIAL INTELLIGENCE》 * |
YINLONG LIU 等: "2D-3D Point Set Registration Based on Global Rotation Search", 《IEEE TRANSACTIONS ON IMAGE PROCESSING * |
李致勋;张英涛;唐降龙;: "基于节点拓扑一致性的2D/3D冠脉血管非刚性配准", 智能计算机与应用 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112465813A (zh) * | 2020-12-17 | 2021-03-09 | 北京工业大学 | 一种基于应力应变的血管内超声弹性分析方法 |
CN112465813B (zh) * | 2020-12-17 | 2024-05-28 | 北京工业大学 | 一种基于应力应变的血管内超声弹性分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111667515B (zh) | 2023-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lynch et al. | Segmentation of the left ventricle of the heart in 3-D+ t MRI data using an optimized nonrigid temporal model | |
KR101855369B1 (ko) | 의료 이미지들 및 임상학적 데이터로부터 생리학적 심장 측정치들을 추정하기 위한 시스템들 및 방법들 | |
McInerney et al. | Deformable models in medical image analysis | |
Frangi et al. | Automatic construction of multiple-object three-dimensional statistical shape models: Application to cardiac modeling | |
US10496729B2 (en) | Method and system for image-based estimation of multi-physics parameters and their uncertainty for patient-specific simulation of organ function | |
US10296809B2 (en) | Method and system for fast patient-specific cardiac electrophysiology simulations for therapy planning and guidance | |
Wang et al. | A review of deformation models in medical image registration | |
CN109993730B (zh) | 3d/2d血管配准方法及装置 | |
Chen et al. | Parametric shape representation by a deformable NURBS model for cardiac functional measurements | |
Heyde et al. | Anatomical image registration using volume conservation to assess cardiac deformation from 3D ultrasound recordings | |
Chen et al. | Shape registration with learned deformations for 3D shape reconstruction from sparse and incomplete point clouds | |
Ahn et al. | Unsupervised motion tracking of left ventricle in echocardiography | |
CN109961508A (zh) | 一种心脏静态三维模型重建方法 | |
CN111724395B (zh) | 心脏图像四维上下文分割方法、设备、存储介质及装置 | |
CN116703994B (zh) | 医学图像配准的方法、计算设备和计算机可读存储介质 | |
CN111667515B (zh) | 基于流形正则化的血管3d/2d弹性配准方法及装置 | |
CN107767409A (zh) | 基于高维表达的一致点漂移配准方法 | |
Goyal et al. | MRI image based patient specific computational model reconstruction of the left ventricle cavity and myocardium | |
Bhan et al. | Patient-specific cardiac computational modeling based on left ventricle segmentation from magnetic resonance images | |
Lötjönen et al. | Four-chamber 3-D statistical shape model from cardiac short-axis and long-axis MR images | |
CN113256693A (zh) | 基于K-means与正态分布变换的多视角配准方法 | |
Sermesant et al. | Biomechanical model construction from different modalities: Application to cardiac images | |
Clarysse et al. | Multi-modality cardiac imaging: processing and analysis | |
Krishnaswamy et al. | A novel 3D-to-3D Diffeomorphic registration algorithm with applications to left ventricle segmentation in MR and Ultrasound sequences | |
Mihalef et al. | Implementation of a patient-specific cardiac model |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |