基于流形正则化的血管3D/2D弹性配准方法及装置制造方法及图纸

基于流形正则化的血管3D/2D弹性配准方法及装置，基于流形正则化的血管3D/2D弹性配准方法，其能够适用于任意具有拓扑结构的点云结构，可以得到最优的血管3D/2D弹性配准结果。方法包括：(1)在已知血管3D和2D点的匹配关系

【技术实现步骤摘要】
基于流形正则化的血管3D/2D弹性配准方法及装置
本专利技术涉及医学图像处理的
，尤其涉及一种基于流形正则化的血管3D/2D弹性配准方法，以及基于流形正则化的血管3D/2D弹性配准装置。
技术介绍
血管3D/2D配准目的是，在术前的3D图像和术中的2D图像中找到对应血管解剖结构的对应关系，并估计3D血管的刚性和弹性变换，使得3D和2D血管在对应结构上达到空间位置的一致。3D和2D血管分别从CTA(计算机断层扫描血管造影)和XRA(X射线血管造影)图像中提取，经3D/2D配准后，3D和2D血管在对应解剖结构上达到空间位置的一致。根据3D/2D配准任务的不同，对其进行任务分解，包括：初始化、刚性变换和弹性变换。弹性变换中的形变模型可以分为两类：一类是基于大量训练数据统计的形变模型，如统计形状模型(SSM)、主动表观模型(AAM)和可形变模型(Morphablemodel)；另一类是基于平滑先验假设的形变模型，如薄板样条模型(TPS)、B样条模型和相干点漂移模型(CPD)。针对血管配准问题，许多研究者在以上这些模型的基础上做出了适用于血管模型的改进。Groher等人于2009年提出了一种基于单视角的血管形变模型，在构建的能量函数中对血管长度变化和非自然弯曲进行惩罚，通过求解目标函数对3D血管点位移的偏导构建目标函数梯度场，利用BFGS优化方法求得最优的弹性变换参数。2010年，Groher等人改进了方法，将配准问题表示为包含以图像信息响应为主的外能和以3D血管长度变化约束为主的内能的能量最小化问题，该方法直接利用2D造影图像的管状结构增强，无需对2D图像进行血管分割。Liao等人将所有血管点视为自由点，在腹主动脉血管的图模型的基础上，构建以距离测度、长度保持项和平滑正则项为主的相似性测度，使用多尺度优化策略求取血管点的位移值。Rivest-Henault等人使用了与类似的内能和外能构建的弹性配准能量函数，不同的是作者利用了血管分叉点处的形变约束来构建外能，并在双视图的冠脉造影图像上求解冠脉3D中心线的位移值。Liu等人利用2DTPS函数实现投影后的3D血管点到2D血管中心线的非刚性配准。Serradell等人针对心脏运动的先验知识，提出了一种合成的冠脉血管树形变模型。利用该先验随机的生成大量的血管中心线训练集，并使用SSM来构建血管形变模型。Baka等人则在真实的临床冠脉4DCTA数据上，提取出3D+t冠脉中心线，并构建真实冠脉的SSM。但是，这些改进的血管模型均只适用于特定的具有拓扑结构的点云结构。
技术实现思路
为克服现有技术的缺陷，本专利技术要解决的技术问题是提供了一种基于流形正则化的血管3D/2D弹性配准方法，其能够适用于任意具有拓扑结构的点云结构，可以得到最优的血管3D/2D弹性配准结果。本专利技术的技术方案是：这种基于流形正则化的血管3D/2D弹性配准方法，其包括以下步骤：(1)在已知血管3D和2D点的匹配关系基础上，表示3D点，表示2D点，i和j分别为点集的索引值；用表示3D点集矩阵，其中为3D点坐标，用表示2D点集矩阵，其中为2D点坐标，两个点集中对应点具有相同的索引值，K为点数量；(2)构建空间变换使得3D点的变换表示为τ(yk)为yk的位移；(3)在血管拓扑流形上的约束描述为在拓扑上相邻的点的位移具有一致性，用流形正则化表示；(4)获得血管3D/2D弹性配准的目标函数；(5)基于梯度优化的算法对目标函数进行求解。本专利技术提出一种基于流形正则化的血管弹性形变模型，该形变模型利用了血管点的位移在其血管拓扑上具有平滑约束的先验，因此该模型适用于任意具有拓扑结构的点云结构；采用基于梯度优化的算法对目标函数进行求解，可以得到最优的血管3D/2D弹性配准结果。还提供了基于流形正则化的血管3D/2D弹性配准装置，其包括：初始化模块，其在已知血管3D和2D点的匹配关系基础上，用表示3D点集矩阵，其中为3D点坐标，用表示2D点集矩阵，其中为2D点坐标，两个点集中对应点具有相同的索引值，K为点数量；构建模块，其构建空间变换使得3D点的变换表示为τ(yk)为yk的位移；流形正则化模块，其在血管拓扑流形上的约束描述为在拓扑上相邻的点的位移具有一致性，用流形正则化表示；目标函数获得模块，其获得血管3D/2D弹性配准的目标函数；求解模块，其基于梯度优化的算法对目标函数进行求解。附图说明图1示出了根据本专利技术的基于流形正则化的血管3D/2D弹性配准方法的流程图。具体实施方式如图1所示，这种基于流形正则化的血管3D/2D弹性配准方法，其包括以下步骤：(1)在已知血管3D和2D点的匹配关系基础上，表示3D点，表示2D点，i和j分别为点集的索引值；用表示3D点集矩阵，其中为3D点坐标，用表示2D点集矩阵，其中为2D点坐标，两个点集中对应点具有相同的索引值，K为点数量；(2)构建空间变换使得3D点的变换表示为τ(yk)为yk的位移；(为作用符，等同于)(3)在血管拓扑流形上的约束描述为在拓扑上相邻的点的位移具有一致性，用流形正则化表示；(4)获得血管3D/2D弹性配准的目标函数；(5)基于梯度优化的算法对目标函数进行求解。本专利技术提出一种基于流形正则化的血管弹性形变模型，该形变模型利用了血管点的位移在其血管拓扑上具有平滑约束的先验，因此该模型适用于任意具有拓扑结构的点云结构；采用基于梯度优化的算法对目标函数进行求解，可以得到最优的血管3D/2D弹性配准结果。优选地，所述步骤(2)中，使用在可再生核希尔伯特空间RKHS中运行的正则化技术，令为一个RKHS，通过Tikhonov正则化来最小化目标函数，利用公式(1)求得定义在的空间变换其中第一项为距离误差，第二项则限制了空间变换的复杂度，λ为权衡两项贡献的正则化系数，为RKHS的范数；对于定义在RKHS的空间变换利用公式(2)表达为任意可再生核函数的线性组合，其中k为正定的核函数为任意非空集合，k满足可再生性，对于任意与核函数k的內积等于其自身<f,k(ci,·)>k＝f(ci,·)，为核函数的输入参数；在此定义下RKHS的范数用其內积表示优选地，所述步骤(2)中，采用高斯函数φ作为核函数，空间变换表示为公式(3)其中为核函数φ输入参数，则为输入的控制点矩阵，C为控制点数量，核函数的形式为β为高斯半径，为空间变换系数；则为系数矩阵，Φ(·,k)为矩阵的第k列，Φ的各元素φik＝φ(ci,yk)。优选地，所述步骤(2)中，空间变换的范数用公式(4)计算，其中矩阵的各元素为uij＝φ(ci,cj)，tr(·)为矩阵的迹，空间变换的范数项约束了由控制点C和系数矩阵W描述的空间形变场的平滑程度。优选地，所述步骤(3)中，在血管拓扑流形上的约束用公式(5)表示优选地，所述步骤(4本文档来自技高网...

【技术保护点】
1.基于流形正则化的血管3D/2D弹性配准方法，其特征在于：其包括以下步骤：/n(1)在已知血管3D和2D点的匹配关系

【技术特征摘要】
1.基于流形正则化的血管3D/2D弹性配准方法，其特征在于：其包括以下步骤：
(1)在已知血管3D和2D点的匹配关系基础上，表示3D点，表示2D点，i和j分别为点集的索引值；用表示3D点集矩阵，其中为3D点坐标，用表示2D点集矩阵，其中为2D点坐标，两个点集中对应点具有相同的索引值，K为点数量；
(2)构建空间变换使得3D点的变换表示为τ(yk)为yk的位移；
(3)在血管拓扑流形上的约束描述为在拓扑上相邻的点的位移具有一致性，用流形正则化表示；
(4)获得血管3D/2D弹性配准的目标函数；
(5)基于梯度优化的算法对目标函数进行求解。


2.根据权利要求1所述的基于流形正则化的血管3D/2D弹性配准方法，其特征在于：所述步骤(2)中，使用在可再生核希尔伯特空间RKHS中运行的正则化技术，令为一个RKHS，通过Tikhonov正则化来最小化目标函数，利用公式(1)求得定义在的空间变换



其中第一项为距离误差，第二项则限制了空间变换的复杂度，λ为权衡两项贡献的正则化系数，为RKHS的范数；
对于定义在RKHS的空间变换利用公式(2)表达为任意可再生核函数的线性组合，



其中k为正定的核函数为任意非空集合，k满足可再生性，对于任意与核函数k的內积等于其自身<f,k(ci,·)>k＝f(ci,·)，为核函数的输入参数；在此定义下RKHS的范数用其內积表示


3.根据权利要求2所述的基于流形正则化的血管3D/2D弹性配准方法，其特征在于：所述步骤(2)中，采用高斯函数φ作为核函数，空间变换表示为公式(3)



其中为核函数φ输入参数，则为输入的控制点矩阵，C为控制点数量，核函数的形式为β为高斯半径，为空间变换系数；

则为系数矩阵，Φ(·,k)为矩阵的第k列，Φ的各元素φik＝φ(ci,yk)。


4.根据权利要求3所述的基于流形正则化的血管3D/2D弹性配准方法，其特征在于：所述步骤(2)中，空间变换的范数用公式(4)计算，



其...

【专利技术属性】
技术研发人员：杨健，朱建军，艾丹妮，范敬凡，王涌天，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京;11