CN109300147A - 血管二维中心线和三维中心线配准的方法及装置 - Google Patents

血管二维中心线和三维中心线配准的方法及装置 Download PDF

Info

Publication number
CN109300147A
CN109300147A CN201811082383.6A CN201811082383A CN109300147A CN 109300147 A CN109300147 A CN 109300147A CN 201811082383 A CN201811082383 A CN 201811082383A CN 109300147 A CN109300147 A CN 109300147A
Authority
CN
China
Prior art keywords
dimensional
center line
dimensional center
centerline
blood vessel
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
CN201811082383.6A
Other languages
English (en)
Other versions
CN109300147B (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 Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201811082383.6A priority Critical patent/CN109300147B/zh
Publication of CN109300147A publication Critical patent/CN109300147A/zh
Application granted granted Critical
Publication of CN109300147B publication Critical patent/CN109300147B/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
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G06T7/344—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • 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
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Graphics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Image Analysis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明实施例提供一种血管二维中心线和三维中心线配准的方法及装置，所述方法包括：从目标血管的二维图像中提取所述目标血管的二维中心线的特征，从所述目标血管的三维图像中提取所述目标血管的三维中心线的特征；对所述二维中心线的特征和所述三维中心线的特征进行建模，获取所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型，根据所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型构建目标函数；基于粒子群优化算法对所述目标函数进行优化，根据优化后的所述目标函数对所述二维中心线和所述三维中心线进行配准。本发明实施例提高了配准的准确性。

Description

血管二维中心线和三维中心线配准的方法及装置
技术领域
本发明实施例属于图像处理技术领域，更具体地，涉及一种血管二维中心线和三维中心线配准的方法及装置。
背景技术
血管疾病严重威胁人类健康，相关的病症包括血管狭窄、血管硬化以及动脉粥样硬化。二维X-ray(X射线)造影图像和对应的三维CT(Computed Tomography，计算机断层扫描)造影图像从不同角度提供血管疾病的细节。因此术前三维CT图像和二维造影图像的融合在基于图像引导的血管微创手术中极为重要。
为了对血管的三维CT图像和二维造影图像进行融合，首先需要对三维CT图像中血管的中心线和二维造影图像中血管的中心线进行配准。二维和三维血管中心线配准的方法大多为基于点云的配准方法进行相应的改变，通过使用高斯混合模型对点云的位置信息进行建模，通过使用最大似然函数或者最小二乘的框架实现目标函数的建立，然后基于梯度下降的方法对目标函数进行优化得到最优结果。
上述基于梯度下降的方法容易受到噪声的影响，即容易陷入局部极小值，使得配准结果发生较大偏差。
发明内容
为克服上述现有的血管二维中心线和三维中心线配准方法配准结果不精确的问题或者至少部分地解决上述问题，本发明实施例提供一种血管二维中心线和三维中心线配准的方法及装置。
根据本发明实施例的第一方面，提供一种血管二维中心线和三维中心线配准的方法，包括：
从目标血管的二维图像中提取所述目标血管的二维中心线的特征，从所述目标血管的三维图像中提取所述目标血管的三维中心线的特征；
对所述二维中心线的特征和所述三维中心线的特征进行建模，获取所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型，根据所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型构建目标函数；
基于粒子群优化算法对所述目标函数进行优化，根据优化后的所述目标函数对所述二维中心线和所述三维中心线进行配准。
根据本发明实施例第二方面提供一种血管二维中心线和三维中心线配准的装置，包括：
提取模块，用于从目标血管的二维图像中提取所述目标血管的二维中心线的特征，从所述目标血管的三维图像中提取所述目标血管的三维中心线的特征；
构建模块，用于对所述二维中心线的特征和所述三维中心线的特征进行建模，获取所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型，根据所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型构建目标函数；
优化模块，用于基于粒子群优化算法对所述目标函数进行优化，根据优化后的所述目标函数对所述二维中心线和所述三维中心线进行配准。
根据本发明实施例的第三个方面，还提供一种电子设备，包括：
至少一个处理器；以及
与所述处理器通信连接的至少一个存储器，其中：
所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的血管二维中心线和三维中心线配准的方法。
根据本发明实施例的第四个方面，还提供一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的血管二维中心线和三维中心线配准的方法。
本发明实施例提供一种血管二维中心线和三维中心线配准的方法及装置，该方法为实现不同维度的血管中心线的配准，先提取血管二维中心线和三维中心线的特征，并对特征进行建模，然后根据特征模型构建目标函数，随后使用粒子群优化算法对目标函数进行优化和目标函数中的特征模型参数再估计，根据优化后的目标函数对不同维度的血管中心线进行配准，从而实现血管二维中心线和三维中心线的自动配准，使用粒子群的启发式算法可以有效的跳出优化目标函数过程中的局部极小值，得到全局最优值，提高了配准的准确性。
附图说明
图1为本发明实施例提供的血管二维中心线和三维中心线配准的方法整体流程示意图；
图2为本发明实施例提供的血管二维中心线和三维中心线配准的装置整体结构示意图；
图3为本发明实施例提供的电子设备整体结构示意图。
具体实施方式
下面结合附图和实施例，对本发明实施例的具体实施方式作进一步详细描述。以下实施例用于说明本发明实施例，但不用来限制本发明实施例的范围。
在本发明实施例的一个实施例中提供一种血管二维中心线和三维中心线配准的方法，图1为本发明实施例提供的血管二维中心线和三维中心线配准的方法整体流程示意图，该方法包括：S101，从目标血管的二维图像中提取目标血管的二维中心线的特征，从目标血管的三维图像中提取目标血管的三维中心线的特征；
其中，目标血管为需要进行二维中心线和三维中心线配准的血管。二维中心线为从目标血管的二维图像中提取出的中心线，其位置为二维坐标。三维中心线为从目标血管的三维图像中提取出的中心线，其位置为三维坐标。进一步从二维中心线中提取出二维中心线的一个或多个特征，从三维中心线中提取出三维中心线的一个或多个特征。其中，二维中心线的特征与三维中心线的特征一一对应，例如二维中心线的位置特征对应有三维中心线的位置特征。本发明实施例不限于二维中心线和三维中心线的提取方法，也不限于二维中心线和三维中心线的各特征的提取方法。
S102，对二维中心线的特征和三维中心线的特征进行建模，获取二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型，根据二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型构建目标函数；
对二维中心线的每个特征和三维中心线的每个特征进行建模，即用模型表示二维中心线的每个特征和三维中心线的每个特征。二维中心线的特征对应的特征模型是指用于表示二维中心线的每个特征的模型。三维中心线的特征对应的特征模型是指用于表示三维中心线的每个特征的模型。根据二维中心线的特征对应的特征模型和三维中心线的特征对应的模型构建目标函数。
S103，基于粒子群优化算法对目标函数进行优化，根据优化后的目标函数对二维中心线和三维中心线进行配准。
其中，粒子群优化算法首先初始化一群随机粒子，即随机解，然后通过迭代找到最优解。在每次迭代中，粒子通过跟踪两个极值进行自我更新，直到达到结束条件。一个极值为粒子本身所找到的最优解，这个解叫做个体极值，另一个极值为整个种群目前找到的最优解，这个极值为全局极值。本发明实施例基于粒子群优化算法对目标函数进行优化和目标函数中特征模型参数的再估计。在目标函数的优化过程中，建立最大期望的迭代数学框架，将粒子群的优化算法融入数学框架之中，在一定搜索区域内通过随机分布粒子进行最优值搜索，在得到本轮的最优值之后，根据期望最大化的方法计算出特征模型的估计参数，直到达到结束条件。经过迭代优化之后获取目标函数最终的最优值，在此情况之下会得到对应的旋转矩阵R和平移矩阵T，根据旋转矩阵R和平移矩阵T对血管二维中心线和血管三维中心线进行配准。其中，血管二维中心线和三维中心线配准是指通过旋转和平移的变换方式使得血管二维中心线和三维中心线的物理坐标得到最大程度的对应。
本发明实施例为实现不同维度的血管中心线的配准，先提取血管二维中心线和三维中心线的特征，并对特征进行建模，然后根据特征模型构建目标函数，随后使用粒子群优化算法对目标函数进行优化和目标函数中的特征模型参数再估计，根据优化后的目标函数对不同维度的血管中心线进行配准，从而实现血管二维中心线和三维中心线的自动配准，使用粒子群的启发式算法可以有效的跳出优化目标函数过程中的局部极小值，得到全局最优值，提高了配准的准确性。
在上述实施例的基础上，本实施例中二维中心线的特征包括方向、位置和半径，三维中心线的特征包括方向、位置和半径。
本发明实施例提取目标血管二维中心线的方向、位置和半径，以及目标血管三维中心线的方向、位置和半径，通过提取基于方向、位置和半径三个特征，使用三个特征对应的特征模型构建目标函数。本发明实施例基于血管中心线的多个特征进行配准，使得目标函数能够准确的衡量不同维度中心线之间的配准程度，提高了配准的精确度。
在上述实施例的基础上，本实施例中从目标血管的二维图像中提取目标血管的二维中心线的特征，从目标血管的三维图像中提取目标血管的三维中心线的特征的步骤具体包括：从目标血管的二维图像中提取目标血管的二维中心线的位置和半径，从目标血管的三维图像中提取目标血管的三维中心线的位置和半径；根据二维中心线中连续两个点的位置获取二维中心线的方向，根据三维中心线中连续两个点的位置获取三维中心线的方向。
具体地，从目标血管的二维中心线中提取出二维中心线的位置，即二维坐标。从目标血管的三维中心线中提取出三维中心线的位置，即三维坐标。从而获取二维中心线和三维中心线上各点的半径，并对得到的不同维度中心线上各点的半径进行邻域高斯平滑。对二维中心线上连续两个点的位置求差，获取二维中心线的方向。对三维中心线上连续两个点的位置求差，获取三维中心线的方向。
在上述实施例的基础上，本实施例中对二维中心线的特征和三维中心线的特征进行建模，获取二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型的步骤具体包括：基于高斯混合模型对二维中心线的位置和三维中心线的位置进行建模，获取二维中心线的位置对应的特征模型和三维中心线的位置对应的特征模型；基于修正的Fisher混合模型对二维中心线的方向和三维中心线的方向进行建模，获取二维中心线的方向对应的特征模型和三维中心线的方向对应的特征模型；基于高斯分布对二维中心线的半径和三维中心线的半径进行建模，获取二维中心线的半径对应的特征模型和三维中心线的半径对应的特征模型。
具体地，在对血管二维中心线和三维中心线建模的过程中，使用高斯混合模型对二维中心线的位置和三维中心线的位置进行建模获取相应的特征模型，使用修正的Fisher混合模型对二维中心线的方向和三维中心线的方向进行建模获取相应的特征模型，使用高斯分布模型对二维中心线的半径和三维中心线的半径进行建模获取相应的特征模型。其中，高斯混合模型是指将多个不同权值的高斯分布叠加而获取的模型。修正的Fisher混合模型的公式如下：
其中，I0表示0阶贝塞尔函数，k为方向的权重，g和v为两个方向。
在上述实施例的基础上，本实施例中根据二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型构建目标函数的步骤具体包括：根据二维中心线的各特征对应的特征模型和三维中心线的各特征对应的特征模型，基于最大似然函数方法构建目标函数。
其中，最大似然函数提供了一种给定观察数据来估计模型参数的方法，首先假设所得的样本服从某一分布，目标是估计出分布中的参数，当得到这一组样本的概率最大时获取到模型的参数值。将二维中心线的各特征对应的特征模型和三维中心线的各特征对应的特征模型融合至最大似然函数框架中构建目标函数。
在上述实施例的基础上，本实施例中通过以下公式根据二维中心线的各特征对应的特征模型和三维中心线的各特征对应的特征模型，基于最大似然函数方法构建目标函数：
其中，N为二维中心线中的总点数，M为三维中心线中像素的总点数，当m＝M+1时，当m≠M+1时，w表示三维中心线中离群点所占比重，a表示三维中心线中点集的投影面积，An表示二维中心线的第n个点，Bm表示三维中心线中的第m个点，p(An,Bm|m)表示An与Bm之间的相似度，表示二维中心线中第n个点的方向向量的转置，表示三维中心线中第m个点的方向向量，表示二维中心线中第n个点的位置，表示三维中心线中第m个点的位置，表示二维中心线中第n个点的半径，表示三维中心线中第m个点的半径，κ,σd,σr分别为κ表示方向的权重，σd表示位置的权重，σr表示半径的权重，表示三维中心线中第m个点的投影，表示三维中心线中第m个点形变后的投影，表示三维中心线的形变，I0表示0阶贝塞尔函数。
在上述各实施例的基础上，本实施例中基于粒子群优化算法对目标函数进行优化的步骤具体包括：根据预设搜索区域中对目标函数进行优化的各粒子的位置，计算出各粒子对应的目标函数值，将最大的目标函数值作为目标函数的当前最优值；根据当前最优值，对二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型进行更新；根据当前最优值和上次迭代中目标函数的最优值，对各粒子的位置进行更新，将更新位置后的粒子作为下一次对目标函数进行优化的粒子。
其中，预设搜索区域为预先设定的搜索区域，通过使用粒子群优化算法对目标函数进行最优值搜寻，根据每次迭代获取的目标函数最优值更新粒子所在的位置，直到目标函数最优值基本不发生改变或者达到预设迭代次数。在使用粒子群优化算法对目标函数进行优化时，首先初始化搜索区域、粒子数量、粒子位置、迭代次数、目标函数最优值的最小变化量、二维中心线的特征对应的特征模型的初始参数和三维中心线的特征对应的特征模型的初始参数。其中，目标函数最优值的最小变化量为相邻两次迭代获取的目标函数最优值之间差值的阈值。该阈值为预先设定的。然后在预设搜索区域中各粒子的位置，计算出各粒子对应的目标函数值，将最大的目标函数值作为目标函数的当前最优值。根据当前最优值，对二维中心线的特征对应的特征模型的初始参数和三维中心线的特征对应的特征模型的初始参数进行更新。判断当前最优值和上次迭代获取的目标函数的最优值之间的差值是否小于最小变化量，若大于最小变化量，则更新各粒子的位置，使用更新位置后的各粒子迭代执行上述对目标函数优化的步骤，直到当前最优值和上次迭代获取的目标函数的最优值之间的差值小于或等于最小变化量。
在本发明实施例的另一个实施例中提供一种血管二维中心线和三维中心线配准的装置，该装置用于实现前述各实施例中的方法。因此，在前述血管二维中心线和三维中心线配准的方法的各实施例中的描述和定义，可以用于本发明实施例中各个执行模块的理解。图2为本发明实施例提供的血管二维中心线和三维中心线配准的装置整体结构示意图，该装置包括提取模块201、构建模块202和优化模块203；其中：
提取模块201用于从目标血管的二维图像中提取目标血管的二维中心线的特征，从目标血管的三维图像中提取目标血管的三维中心线的特征；
其中，目标血管为需要进行二维中心线和三维中心线配准的血管。二维中心线为从目标血管的二维图像中提取出的中心线，其位置为二维坐标。三维中心线为从目标血管的三维图像中提取出的中心线，其位置为三维坐标。提取模块201从二维中心线中提取出二维中心线的一个或多个特征，从三维中心线中提取出三维中心线的一个或多个特征。其中，二维中心线的特征与三维中心线的特征一一对应，例如二维中心线的位置特征对应有三维中心线的位置特征。本发明实施例不限于二维中心线和三维中心线的提取方法，也不限于二维中心线和三维中心线的各特征的提取方法。
构建模块202用于对二维中心线的特征和三维中心线的特征进行建模，获取二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型，根据二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型构建目标函数；
构建模块202对二维中心线的每个特征和三维中心线的每个特征进行建模，即用模型表示二维中心线的每个特征和三维中心线的每个特征。二维中心线的特征对应的特征模型是指用于表示二维中心线的每个特征的模型。三维中心线的特征对应的特征模型是指用于表示三维中心线的每个特征的模型。根据二维中心线的特征对应的特征模型和三维中心线的特征对应的模型构建目标函数。
优化模块203用于基于粒子群优化算法对目标函数进行优化，根据优化后的目标函数对二维中心线和三维中心线进行配准。
其中，粒子群优化算法首先初始化一群随机粒子，即随机解，然后通过迭代找到最优解。在每次迭代中，粒子通过跟踪两个极值进行自我更新，直到达到结束条件。一个极值为粒子本身所找到的最优解，这个解叫做个体极值，另一个极值为整个种群目前找到的最优解，这个极值为全局极值。优化模块203基于粒子群优化算法对目标函数进行优化和目标函数中特征模型参数的再估计。在目标函数的优化过程中，建立最大期望的迭代数学框架，将粒子群的优化算法融入数学框架之中，在一定搜索区域内通过随机分布粒子进行最优值搜索，在得到本轮的最优值之后，根据期望最大化的方法计算出特征模型的估计参数，直到达到结束条件。根据优化后的目标函数对血管二维中心线和三维中心线进行配准。其中，血管二维中心线和三维中心线配准是指通过X轴偏移、Y轴偏移和Z轴偏移的变换方式使得血管二维中心线和三维中心线的物理坐标得到最大程度的对应。
本发明实施例为实现不同维度的血管中心线的配准，先提取血管二维中心线和三维中心线的特征，并对特征进行建模，然后根据特征模型构建目标函数，随后使用粒子群优化算法对目标函数进行优化和目标函数中的特征模型参数再估计，根据优化后的目标函数对不同维度的血管中心线进行配准，从而实现血管二维中心线和三维中心线的自动配准，使用粒子群的启发式算法可以有效的跳出优化目标函数过程中的局部极小值，得到全局最优值，提高了配准的准确性。
在上述实施例的基础上，本实施例中二维中心线的特征包括方向、位置和半径，三维中心线的特征包括方向、位置和半径。
在上述实施例的基础上，本实施例中提取模块具体用于：从目标血管的二维图像中提取目标血管的二维中心线的位置和半径，从目标血管的三维图像中提取目标血管的三维中心线的位置和半径；根据二维中心线中连续两个点的位置获取二维中心线的方向，根据三维中心线中连续两个点的位置获取三维中心线的方向。
在上述实施例的基础上，本实施例中构建模块具体用于：基于高斯混合模型对二维中心线的位置和三维中心线的位置进行建模，获取二维中心线的位置对应的特征模型和三维中心线的位置对应的特征模型；基于修正的Fisher混合模型对二维中心线的方向和三维中心线的方向进行建模，获取二维中心线的方向对应的特征模型和三维中心线的方向对应的特征模型；基于高斯分布对二维中心线的半径和三维中心线的半径进行建模，获取二维中心线的半径对应的特征模型和三维中心线的半径对应的特征模型。
在上述实施例的基础上，本实施例中构建模块具体用于：根据二维中心线的各特征对应的特征模型和三维中心线的各特征对应的特征模型，基于最大似然函数方法构建目标函数。
在上述各实施例的基础上，本实施例中构建模块具体通过以下公式根据二维中心线的各特征对应的特征模型和三维中心线的各特征对应的特征模型，基于最大似然函数方法构建目标函数：
其中，N为二维中心线中的总点数，M为三维中心线中像素的总点数，当m＝M+1时，当m≠M+1时，w表示三维中心线中离群点所占比重，a表示三维中心线中点集的投影面积，An表示二维中心线的第n个点，Bm表示三维中心线中的第m个点，p(An,Bm|m)表示An与Bm之间的相似度，表示二维中心线中第n个点的方向向量的转置，表示三维中心线中第m个点的方向向量，表示二维中心线中第n个点的位置，表示三维中心线中第m个点的位置，表示二维中心线中第n个点的半径，表示三维中心线中第m个点的半径，κ,σd,σr分别为κ表示方向的权重，σd表示位置的权重，σr表示半径的权重，表示三维中心线中第m个点的投影，表示三维中心线中第m个点形变后的投影，表示三维中心线的形变，I0表示0阶贝塞尔函数。
在上述各实施例的基础上，本实施例中优化模块具体用于：根据预设搜索区域中对目标函数进行优化的各粒子的位置，计算出各粒子对应的目标函数值，将最大的目标函数值作为目标函数的当前最优值；根据当前最优值，对二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型进行更新；根据当前最优值和上次迭代中目标函数的最优值，对各粒子的位置进行更新，将更新位置后的粒子作为下一次对目标函数进行优化的粒子。
本实施例提供一种电子设备，图3为本发明实施例提供的电子设备整体结构示意图，该设备包括：至少一个处理器301、至少一个存储器302和总线303；其中，
处理器301和存储器302通过总线303完成相互间的通信；
存储器302存储有可被处理器301执行的程序指令，处理器调用程序指令能够执行上述各方法实施例所提供的方法，例如包括：从目标血管的二维图像中提取目标血管的二维中心线的特征，从目标血管的三维图像中提取目标血管的三维中心线的特征；对二维中心线的特征和三维中心线的特征进行建模，获取二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型，根据二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型构建目标函数；基于粒子群优化算法对目标函数进行优化，根据优化后的目标函数对二维中心线和三维中心线进行配准。
本实施例提供一种非暂态计算机可读存储介质，非暂态计算机可读存储介质存储计算机指令，计算机指令使计算机执行上述各方法实施例所提供的方法，例如包括：从目标血管的二维图像中提取目标血管的二维中心线的特征，从目标血管的三维图像中提取目标血管的三维中心线的特征；对二维中心线的特征和三维中心线的特征进行建模，获取二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型，根据二维中心线的特征对应的特征模型和三维中心线的特征对应的特征模型构建目标函数；基于粒子群优化算法对目标函数进行优化，根据优化后的目标函数对二维中心线和三维中心线进行配准。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的电子设备实施例仅仅是示意性的，其中作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分的方法。
最后，本申请的方法仅为较佳的实施方案，并非用于限定本发明实施例的保护范围。凡在本发明实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明实施例的保护范围之内。

Claims (10)

1.一种血管二维中心线和三维中心线配准的方法，其特征在于，包括：
从目标血管的二维图像中提取所述目标血管的二维中心线的特征，从目标血管的三维图像中提取所述目标血管的三维中心线的特征；
对所述二维中心线的特征和所述三维中心线的特征进行建模，获取所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型，根据所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型构建目标函数；
基于粒子群优化算法对所述目标函数进行优化，根据优化后的所述目标函数对所述二维中心线和所述三维中心线进行配准。
2.根据权利要求1所述的方法，其特征在于，所述二维中心线的特征包括方向、位置和半径，所述三维中心线的特征包括方向、位置和半径。
3.根据权利要求2所述的方法，其特征在于，从目标血管的二维图像中提取所述目标血管的二维中心线的特征，从所述目标血管的三维图像中提取所述目标血管的三维中心线的特征的步骤具体包括：
从目标血管的二维图像中提取所述目标血管的二维中心线的位置和半径，从所述目标血管的三维图像中提取所述目标血管的三维中心线的位置和半径；
根据所述二维中心线中连续两个点的位置获取所述二维中心线的方向，根据所述三维中心线中连续两个点的位置获取所述三维中心线的方向。
4.根据权利要求2所述的方法，其特征在于，对所述二维中心线的特征和所述三维中心线的特征进行建模，获取所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型的步骤具体包括：
基于高斯混合模型对所述二维中心线的位置和所述三维中心线的位置进行建模，获取所述二维中心线的位置对应的特征模型和所述三维中心线的位置对应的特征模型；
基于修正的Fisher混合模型对所述二维中心线的方向和所述三维中心线的方向进行建模，获取所述二维中心线的方向对应的特征模型和所述三维中心线的方向对应的特征模型；
基于高斯分布对所述二维中心线的半径和所述三维中心线的半径进行建模，获取所述二维中心线的半径对应的特征模型和所述三维中心线的半径对应的特征模型。
5.根据权利要求2所述的方法，其特征在于，根据所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型构建目标函数的步骤具体包括：
根据所述二维中心线的各特征对应的特征模型和所述三维中心线的各特征对应的特征模型，基于最大似然函数方法构建目标函数。
6.根据权利要求5所述的方法，其特征在于，通过以下公式根据所述二维中心线的各特征对应的特征模型和所述三维中心线的各特征对应的特征模型，基于最大似然函数方法构建目标函数：
其中，N为所述二维中心线中的总点数，M为所述三维中心线中像素的总点数，当m＝M+1时，当m≠M+1时，w表示所述三维中心线中离群点所占比重，a表示所述三维中心线中点集的投影面积，An表示所述二维中心线的第n个点，Bm表示所述三维中心线中的第m个点，p(An,Bm|m)表示An与Bm之间的相似度，表示所述二维中心线中第n个点的方向向量的转置，表示所述三维中心线中第m个点的方向向量，表示所述二维中心线中第n个点的位置，表示所述三维中心线中第m个点的位置，表示所述二维中心线中第n个点的半径，表示所述三维中心线中第m个点的半径，κ,σd,σr分别为κ表示方向的权重，σd表示位置的权重，σr表示半径的权重，表示所述三维中心线中第m个点的投影，表示所述三维中心线中第m个点形变后的投影，表示所述三维中心线的形变，I0表示0阶贝塞尔函数。
7.根据权利要求1-6任一所述的方法，其特征在于，基于粒子群优化算法对所述目标函数进行优化的步骤具体包括：
根据预设搜索区域中对目标函数进行优化的各粒子的位置，计算出各粒子对应的目标函数值，将最大的所述目标函数值作为所述目标函数的当前最优值；
根据所述当前最优值，对所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型进行更新；
根据所述当前最优值和上次迭代中所述目标函数的最优值，对各所述粒子的位置进行更新，将更新位置后的所述粒子作为下一次对目标函数进行优化的粒子。
8.一种血管二维中心线和三维中心线配准的装置，其特征在于，包括：
提取模块，用于从目标血管的二维图像中提取所述目标血管的二维中心线的特征，从所述目标血管的三维图像中提取所述目标血管的三维中心线的特征；
构建模块，用于对所述二维中心线的特征和所述三维中心线的特征进行建模，获取所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型，根据所述二维中心线的特征对应的特征模型和所述三维中心线的特征对应的特征模型构建目标函数；
优化模块，用于基于粒子群优化算法对所述目标函数进行优化，根据优化后的所述目标函数对所述二维中心线和所述三维中心线进行配准。
9.一种电子设备，其特征在于，包括：
至少一个处理器、至少一个存储器和总线；其中，
所述处理器和存储器通过所述总线完成相互间的通信；
所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行如权利要求1至7任一所述的方法。
10.一种非暂态计算机可读存储介质，其特征在于，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行如权利要求1至7任一所述的方法。
CN201811082383.6A 2018-09-17 2018-09-17 血管二维中心线和三维中心线配准的方法及装置 Active CN109300147B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811082383.6A CN109300147B (zh) 2018-09-17 2018-09-17 血管二维中心线和三维中心线配准的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811082383.6A CN109300147B (zh) 2018-09-17 2018-09-17 血管二维中心线和三维中心线配准的方法及装置

Publications (2)

Publication Number Publication Date
CN109300147A true CN109300147A (zh) 2019-02-01
CN109300147B CN109300147B (zh) 2020-10-16

Family

ID=65163237

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811082383.6A Active CN109300147B (zh) 2018-09-17 2018-09-17 血管二维中心线和三维中心线配准的方法及装置

Country Status (1)

Country Link
CN (1) CN109300147B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109993730A (zh) * 2019-03-20 2019-07-09 北京理工大学 3d/2d血管配准方法及装置
CN110288637A (zh) * 2019-06-13 2019-09-27 北京理工大学 多角度dsa造影图像血管匹配方法及装置
CN110840561A (zh) * 2019-11-11 2020-02-28 郑健青 一种基于人工智能与图论算法的手术导航辅助系统
CN111476791A (zh) * 2019-05-21 2020-07-31 北京昆仑医云科技有限公司 图像处理方法、图像处理装置和非暂时性计算机可读介质
CN111612778A (zh) * 2020-05-26 2020-09-01 上海交通大学 一种术前cta与术中x光冠状动脉配准方法
CN113763543A (zh) * 2021-09-17 2021-12-07 北京理工大学 基于三维体素化结构的血管重构方法、评价方法及系统
CN117495667A (zh) * 2023-12-29 2024-02-02 柏意慧心(杭州)网络科技有限公司 用于确定血管中心线的方法、计算设备和介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060036167A1 (en) * 2004-07-03 2006-02-16 Shina Systems Ltd. Vascular image processing
CN104361626A (zh) * 2014-09-29 2015-02-18 北京理工大学 基于混合匹配策略的皮下静脉三维重建方法
CN104835112A (zh) * 2015-05-07 2015-08-12 厦门大学 一种肝脏多相期ct图像融合方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060036167A1 (en) * 2004-07-03 2006-02-16 Shina Systems Ltd. Vascular image processing
CN104361626A (zh) * 2014-09-29 2015-02-18 北京理工大学 基于混合匹配策略的皮下静脉三维重建方法
CN104835112A (zh) * 2015-05-07 2015-08-12 厦门大学 一种肝脏多相期ct图像融合方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JULIEN JOMIER等: "3D/2D Model-to-Image Registration Applied to TIPS Surgery", 《MED IMAGE COMPUT COMPUT ASSIST INTERV.》 *
SONGBO LIU等: "A 3D/2D registration of the coronary arteries based on tree topologyconsistency matching", 《BIOMEDICAL SIGNAL PROCESSING AND CONTROL》 *
丛伟建: "应用于肝脏介入诊疗的超声图像分析与导航方法研究", 《中国博士学位论文全文数据库_信息科技辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109993730A (zh) * 2019-03-20 2019-07-09 北京理工大学 3d/2d血管配准方法及装置
CN111476791A (zh) * 2019-05-21 2020-07-31 北京昆仑医云科技有限公司 图像处理方法、图像处理装置和非暂时性计算机可读介质
CN111476791B (zh) * 2019-05-21 2024-02-09 科亚医疗科技股份有限公司 图像处理方法、图像处理装置和非暂时性计算机可读介质
CN110288637A (zh) * 2019-06-13 2019-09-27 北京理工大学 多角度dsa造影图像血管匹配方法及装置
CN110288637B (zh) * 2019-06-13 2021-09-21 北京理工大学 多角度dsa造影图像血管匹配方法及装置
CN110840561A (zh) * 2019-11-11 2020-02-28 郑健青 一种基于人工智能与图论算法的手术导航辅助系统
CN111612778A (zh) * 2020-05-26 2020-09-01 上海交通大学 一种术前cta与术中x光冠状动脉配准方法
CN113763543A (zh) * 2021-09-17 2021-12-07 北京理工大学 基于三维体素化结构的血管重构方法、评价方法及系统
CN113763543B (zh) * 2021-09-17 2024-07-09 北京理工大学 基于三维体素化结构的血管重构方法、评价方法及系统
CN117495667A (zh) * 2023-12-29 2024-02-02 柏意慧心(杭州)网络科技有限公司 用于确定血管中心线的方法、计算设备和介质
CN117495667B (zh) * 2023-12-29 2024-03-22 柏意慧心(杭州)网络科技有限公司 用于确定血管中心线的方法、计算设备和介质

Also Published As

Publication number Publication date
CN109300147B (zh) 2020-10-16

Similar Documents

Publication Publication Date Title
CN109300147B (zh) 血管二维中心线和三维中心线配准的方法及装置
EP3365869B1 (en) System and method for image registration in medical imaging system
JP6433149B2 (ja) 姿勢推定装置、姿勢推定方法およびプログラム
Baka et al. Oriented Gaussian mixture models for nonrigid 2D/3D coronary artery registration
CN107133946B (zh) 医学图像处理方法、装置及设备
JP6505124B2 (ja) 適応性放射線治療における自動輪郭抽出システム及び方法
US20110026794A1 (en) Deformable 2D-3D Registration of Structure
JP2020503603A (ja) 解剖学的画像処理における確率的セグメンテーションのためのシステム及び方法
CN111095354A (zh) 经改进的3-d血管树表面重构
CN109993730B (zh) 3d/2d血管配准方法及装置
JP2018530401A5 (zh)
EP2688040B1 (en) 3D registration of tomographic models using particularly chosen points
CN109102490A (zh) 自动图像注册质量评估
US10762648B2 (en) Image processing apparatus, image processing method, image processing system, and program
US9984311B2 (en) Method and system for image segmentation using a directed graph
CN106133789B (zh) 用于分割感兴趣区域的图像处理装置和方法
CN106600624B (zh) 基于粒子群的粒子滤波视频目标跟踪方法
US11631177B2 (en) Machine learning device, estimation device, non-transitory computer readable medium, and learned model
CN113870215B (zh) 中线提取方法及装置
CN111862046A (zh) 一种心脏冠脉剪影中导管位置判别系统和方法
Fahmy et al. Myocardial segmentation using contour-constrained optical flow tracking
CN110751722A (zh) 同时定位建图方法及装置
Oh et al. A novel method to quantify and compare anatomical shape: application in cervix cancer radiotherapy
Aranda et al. A flocking based method for brain tractography
CN114787867A (zh) 用于医学图像配准的器官变形补偿

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