CN107945226A - 一种医学影像解剖结构的中心线确定方法及装置 - Google Patents

一种医学影像解剖结构的中心线确定方法及装置 Download PDF

Info

Publication number
CN107945226A
CN107945226A CN201711091221.4A CN201711091221A CN107945226A CN 107945226 A CN107945226 A CN 107945226A CN 201711091221 A CN201711091221 A CN 201711091221A CN 107945226 A CN107945226 A CN 107945226A
Authority
CN
China
Prior art keywords
target
image
point
center
estimation
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
CN201711091221.4A
Other languages
English (en)
Other versions
CN107945226B (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.)
Neusoft Medical Systems Co Ltd
Original Assignee
Neusoft Medical Systems 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 Neusoft Medical Systems Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN201711091221.4A priority Critical patent/CN107945226B/zh
Publication of CN107945226A publication Critical patent/CN107945226A/zh
Application granted granted Critical
Publication of CN107945226B publication Critical patent/CN107945226B/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
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20104Interactive definition of region of interest [ROI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geometry (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Analysis (AREA)

Abstract

本发明实施例公开了一种医学影像解剖结构的中心线确定方法,采用多级分辨率的分类器估算出目标中心点,首先使用低分辨率的分类器进行估算,通过估算后获得的位置偏移量来确定出目标中心点,而后采用高分辨率的分类器进行估算,此次估算时ROI图像的中心点为前次估算后获得的目标中心点,从而,逐渐获得靠近中心线的目标中心点,通过这些目标中心点精确地确定出中心线。该方法中,采用不同分辨率的分类器逐步确定目标中心点,兼顾中心线确定的速度和精度,可以适用于不同结构特点的解剖结构的中心线的确定,且对图像不均匀及噪声等干扰因素不敏感,中心线确定的准确性和稳定性高。

Description

一种医学影像解剖结构的中心线确定方法及装置
技术领域
本发明涉及医学影像的图像处理技术领域,尤其涉及一种医学影像解剖结构的中心线确定方法及装置。
背景技术
医学影像扫描是医学中常用的辅助诊断手段,在利用医学影像进行诊断的过程中,通常会需要在影像上标记出某些特定解剖结构的中心线位置,以便于为进一步的影像扫描、分析或诊断提供位置参照和依据。
目前,中心线位置的确定方法主要有细化法和关键点法,细化法是通过细化算法定义一定大小的模板,并在计算过程中使得该模板沿着扫描方向移动按照预设规则删除或保留目标图像的像素,最终获得单像素宽的连通区域,从而获得中心线,该方法原理简单且计算量小,但对解剖结构的形状有要求,仅适合线条类的简单结构。关键点法是在目标图像上选定多个关键点,在目标图像上识别关键点,或者在目标图像内分割出解剖结构,而后在分割出的解剖结构上识别关键点,最后以这些关键点的连线作为中心线,该方法对于目标图像及解剖结构的影响要求较高,当目标图像质量较低、关键点不明确或解剖结构轮廓不明确时,难以保证中心线的识别精度。
发明内容
本发明提供了一种医学影像解剖结构的中心线确定方法及装置,适用于不同的解剖结构,中心线识别的准确性和稳定性高。
本发明提供了一种医学影像解剖结构的中心线确定方法,预先设置有多个不同预设分辨率的分类器,所述分类器用于以目标图像上的ROI图像为输入,估算出ROI图像的中心点相对于目标图像中目标区域的中心线的位置偏移量;所述方法包括:
获得扫描图像,扫描图像具有目标区域;
进行中心点的估算,所述中心点的估算包括:
进行位置偏移量的估算以及目标中心点的确定,所述位置偏移量的估算包括:以目标图像上确定的多个ROI图像为输入,利用分类器估算出ROI图像的中心点相对于目标区域中心线的位置偏移量,目标图像为所利用的分类器的分辨率下的扫描图像;所述目标中心点的确定包括:通过位置偏移量,获得ROI图像的中心点在目标区域上对应的目标中心点;
以本次估算获得的目标中心点为下次估算中ROI图像的中心点,重复进行中心点的估算的步骤,直到到达预设的估算次数,其中,每次中心点的估算中采用的分类器的分辨率依次增大,最后一次中心点的估算中采用的分类器的分辨率为扫描图像的分辨率;
根据最后一次估算获得的目标中心点,获得扫描图像目标区域的中心线。
可选地,在所述中心点的估算的步骤中,对位置偏移量估算以及目标中心点的确定的次数大于或者等于1次,且本次修正步骤中获得的目标中心点作为下次修正步骤中ROI图像的中心点。
可选地,所述根据最后一次估算获得的目标中心点,获得扫描图像目标区域的中心线,包括:
对最后一次中心点的估算中获得的目标中心点进行直线拟合,以获得扫描图像目标区域的中心线。
可选地,所述预设分辨率的分类器的确定方法包括:
获得预设分辨率的扫描图像样本,扫描图像样本的目标区域上标记有中心线;
在所述扫描图像样本上确定ROI图像样本,并记录ROI图像样本中心点相对于目标区域中心线的位置偏移量;
通过所述ROI图像样本和位置偏移量进行分类器模型的训练,以获得所述预设分辨率的分类器。
可选地,所述分类器模型包括特征提取器和特征分类器,所述特征提取器由深度卷积神经网络构成,所述特征分类器由全连接层构成。
一种医学影像解剖结构的中心线确定装置,所述装置包括:
多个不同预设分辨率的分类器,所述分类器用于以目标图像上的ROI图像为输入,估算出ROI图像的中心点相对于目标图像中目标区域的中心线的位置偏移量;
扫描图像获取单元,用于获得扫描图像,扫描图像具有目标区域;
中心点估算单元,用于进行中心点的估算,所述中心点的估算包括:
进行位置偏移量的估算以及目标中心点的确定,所述位置偏移量的估算包括:以目标图像上确定的多个ROI图像为输入,利用分类器估算出ROI图像的中心点相对于目标区域中心线的位置偏移量,目标图像为所利用的分类器的分辨率下的扫描图像;所述目标中心点的确定包括:通过位置偏移量,获得ROI图像的中心点在目标区域上对应的目标中心点;
以本次估算获得的目标中心点为下次估算中ROI图像的中心点,重复进行中心点的估算的步骤,直到到达预设的估算次数,其中,每次中心点的估算中采用的分类器的分辨率依次增大,最后一次中心点的估算中采用的分类器的分辨率为扫描图像的分辨率;中心线确定单元,用于根据最后一次估算获得的目标中心点,获得扫描图像目标区域的中心线。
可选地,在所述中心点的估算的步骤中,执行多次进行位置偏移量的估算以及目标中心点的确定的修正步骤,且本次修正步骤中获得的目标中心点作为下次修正步骤中ROI图像的中心点。
可选地,所述中心线确定单元中,对最后一次中心点的估算中获得的目标中心点进行直线拟合,以获得扫描图像目标区域的中心线。
可选地,还包括:分类器确定单元,用于获得预设分辨率的扫描图像样本,扫描图像样本的目标区域上标记有中心线;在所述扫描图像样本上确定ROI图像样本,并记录ROI图像样本中心点相对于目标区域中心线的位置偏移量;通过所述ROI图像样本和位置偏移量进行分类器模型的训练,以获得所述预设分辨率的分类器。
可选地,所述分类器模型包括特征提取器和特征分类器,所述特征提取器由深度卷积神经网络构成,所述特征分类器由全连接层构成。
本发明实施例提供的医学影像解剖结构的中心线确定方法,采用多级分辨率的分类器估算出目标中心点,首先使用低分辨率的分类器进行估算,通过估算后获得的位置偏移量来确定出目标中心点,而后采用高分辨率的分类器进行估算,此次估算时ROI图像的中心点为前次估算后获得的目标中心点,从而,逐渐获得靠近中心线的目标中心点,通过这些目标中心点精确地确定出中心线。该方法中,采用不同分辨率的分类器逐步确定目标中心点,兼顾中心线确定的速度和精度,可以适用于不同结构特点的解剖结构的中心线的确定,且对图像不均匀及噪声等干扰因素不敏感,中心线确定的准确性和稳定性高。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为根据本发明实施例的医学影像解剖结构的中心线确定方法的流程示意图;
图2为根据本发明实施例的医学影像解剖结构的中心线确定方法中扫描图像的采样点示意图;
图3为根据本发明实施例的医学影像解剖结构的中心线确定方法中分类器的结构示意图;
图4为采用本发明实施例的方法确定中心线的各过程的图像示意图;
图5为采用本发明实施例的方法在不同质量的扫描图像中确定的中心线示意图;
图6为根据本发明实施例的医学影像解剖结构的中心线确定装置的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本发明提供了一种医学影像解剖结构的中心线确定方法,通过预先设置的不同分辨率的分类器进行中心点位置的估算,分类器是以目标图像上的ROI图像为输入,通过分类器估算出ROI图像的中心点相对于扫描图像中目标区域的中心线的位置偏移量。具体的,参考图1所示,该方法包括:
S01,获得扫描图像,扫描图像具有目标区域;
S02,进行中心点的估算,所述中心点的估算包括:
进行S021位置偏移量的估算以及S022目标中心点的确定,所述位置偏移量的估算包括:以目标图像上确定的多个ROI图像为输入,利用分类器估算出ROI图像的中心点相对于目标区域中心线的位置偏移量,目标图像为所利用的分类器的分辨率下的扫描图像;所述目标中心点的确定包括:通过位置偏移量,获得ROI图像的中心点在目标区域上的目标中心点;
以目标中心点为下次估算中ROI图像的中心点,重复进行中心点的估算的步骤,直到到达预设的估算次数,其中,每次中心点的估算中采用的分类器的分辨率依次增大,最后一次中心点的估算中采用的分类器的分辨率为扫描图像的分辨率;
S03,根据最后一次中心点的估算获得的目标中心点,获得扫描图像目标区域的中心线。
在该方法中,通过多次估算来获得扫描图像目标区域的中心点,每次估算采用不同分辨率的分类器进行,且先采用低分辨率的分类器进行估算,通过估算后获得的位置偏移量来确定出目标中心点,而后采用高分辨率的分类器进行估算,每次估算获得的目标中心点作为下一次估算时的ROI图像的中心点,从而逐渐获得靠近中心线的目标中心点,通过这些目标中心点精确地确定出中心线。
在本发明的实施例中,每次估算中采用某一分辨率的分类器,即每一分类器具有特定的分辨率,也就是说,该分类器适用于特定分辨率下的目标图像的中心点的确定。该目标图像可以为经过扫描重建后直接获得的原始扫描图像,该原始扫描图像具有较高的分辨率,目标图像也可以为原始扫描图像经过降采样处理之后的图像,降采样之后的图像较原始扫描图像具有更低的分辨率,该扫描图像可以为经过任意类型的医学扫描成像后获得的影像,医学扫描例如可以为核磁共振扫描成像或CT扫描成像等。目标图像中的目标区域即为需要确定中心线的解剖结构所在的区域。
本发明中的分类器用于根据目标图像上的ROI(Region Of Interest,感兴趣区域),估算出该ROI图像的中心点相对于目标区域的中心线的位置偏移量,也就是说,分类器的输入为目标图像上的ROI区域,输出为该ROI区域相对于目标区域的中心线的位置偏移量,通过该偏移量就可以确定该ROI区域相对应的中心点位置。
具体的,可以采用如下方法确定预定分辨率的分类器。
在S101,获得预设分辨率的扫描图像样本,扫描图像样本的目标区域上标记有中心线。
扫描图像样本为大量已经标记有中心线的扫描影像,这些扫描图像样本为包含同一类型的解剖结构的影像且这些图像样本的分辨率是相同的。扫描图像样本的分辨率可以为W1*H1*S1,其中,W1为图像宽度,H1为图像高度,S1为层数。
在S102,在所述扫描图像样本上确定ROI图像样本,并记录ROI图像样本中心点相对于目标区域中心线的偏移量。
ROI图像样本为在扫描图像样本上确定的区域,可以随机地在扫描图像样本上确定分辨率为w1*h1*s1的图像区域,作为ROI图像样本,采样时,可以先确定采样点,而后以采样点为中心,确定w1*h1*s1的图像区域作为ROI图像样本。更优地,为了提高样本数据的有效性,在选取ROI区域时,将完全为背景区域的排除在选取范围之外,此外,为了提高中心点估算的准确性,可以提高中心线处的采样密度。
确定ROI图像样本之后,需要记录ROI图像样本中心点相对于目标区域中心线的偏移量,偏移量作为训练标签,用于模型训练。对于不同的解剖结构,偏移量可以是一维或多维的,在当解剖结构的方向相对固定时,偏移量可以为一维,参考图2所示,该示例中头部的大致方向是相对固定的,可以仅记录ROI图像样本中心点沿水平方向上的位置偏移量dx。当然,在解剖结构的方向不确定时,可以记录多个方向上的位置偏离量,例如dx、dy和dz等沿不同方向上的偏移量。在本发明实施例中,位置偏移量可以为偏移的实际距离,也可以为进行归一化之后的距离系数,例如可以为:dx=2dx/W1,dy=2dy/H1,dz=2dz/S1。
在S103,通过所述ROI图像样本和位置偏移量进行分类器模型的训练,以获得所述预设分辨率的分类器。
分类器模型可以采用合适的模型,本发明优选实施例中,参考图3所示,分类器模型包括特征提取器和特征分类器,所述特征提取器由深度卷积神经网络构成,该网络包括若干个模块,每个模块都有卷积层和池化层组成,所述特征分类器由全连接层构成。在具体训练时,将ROI图像样本为输入,并以位置偏移量作为标签,对分类器模型的参数进行修正,经过多次迭代,在输出结果收敛之后,则完成训练,此时,输入任意的ROI图像得到的位置偏移量都与标签一致或接近。
这样,就可以获得预设分辨率的分类器,在不同的应用中,需要不同分辨率的分类器,可以通过对不同分辨率的扫描图像样本分别进行训练,获得相应分辨率的分类器,低分辨率的样本可以通过高分辨率的样本降采样后获得。在一个具体的示例中,参考图2所示,用于训练分类器的扫描图像样本为头部核磁共振定位扫描图像,样本的分辨率W1*H1*S1为512*512*3,ROI图像样本的分辨率为64*64*1,通过上述的方法进行训练之后,可以获得分辨率为512*512*3、解剖结构为头部的分类器。同样地,还需要分辨率为256*256*3的头部的分类器,则可以将分辨率为512*512*3的扫描图像样本进行降采样,得到分辨率为256*256*3的扫描图像样本,或者其他途径获得该低分辨率的样本,通过上述方法进行训练之后,即可以获得分辨率为256*256*3、解剖结构为头部的分类器。
以上对本发明实施例的分类器以及分类器的确定方法进行了详细的描述。基于上述的分类器,进行医学影像解剖结构的中心线的确定。为了便于理解,以下将以预设有两个不同分辨率的分类器模型的实施例对该方法进行详细的描述,即第一分辨率的分类器和第二分辨率的分类器,其中,第二分辨率即为扫描图像的分辨率,为高分辨率,第一分辨率为低分辨率,即低于扫描图像的分辨率。
在步骤S201,获得扫描图像,扫描图像具有目标区域。
该扫描图像为进行医学影像扫描之后获得的图像,最终的中心线确定在该扫描图像上,扫描图像上的目标区域为需要确定出中心线的区域,即需要标定中心线的解剖结构区域。
在步骤S202,进行第一次中心点的估算。
一次中心点的估算包括:
S2021,以第一分辨率下的扫描图像为第一目标图像,在所述第一目标图像上确定出多个ROI图像,并以确定出的ROI图像为输入,利用第一分辨率的分类器估算出ROI图像的中心点相对于第一目标图像的目标区域中心线的位置偏移量P(1,m)。
S2022,通过位置偏移量P(1,m),获得ROI图像中心点S(1,m)在目标区域上对应的目标中心点S(1,m)’,m从1至M,M≥1。
根据需要,在本发明实施例中,可以执行一次或多次步骤S2021和S2022,为了便于描述,将步骤S2021和步骤S2022记做修正步骤,在执行多次修正步骤时,本次修正步骤中获得的目标中心点S(1,m)’作为下次修正步骤中ROI图像的中心点S(1,m+1),M为一次估算中执行修正的总次数,可以根据需要确定修正的次数M,而最后一次修正步骤获得的目标中心点S(1,M)’则作为下一次中心点估算时的ROI图像的中心点S(2,1)。
在具体的实施例中,第一分辨率低于扫描图像的分辨率,可以将扫描图像进行降采样,从而获得第一分辨率下的扫描图像,为了便于描述记做第一目标图像。可以通过先确定多个采样点,而后分别以采样点为中心点确定ROI图像,采样点即为ROI图像的中心点S(1,1)。通常地,采样点为多个,可以通过随机的方式或者指定的方式在第一目标图像上选择采样点,优选地,可以在目标区域的两侧或轮廓处确定采样点,在第一目标图像上确定出多个采样点后,以这些采样点为中心点,可以分别确定出预设大小的ROI图像,预设大小的ROI图像的分辨率例如可以为64*64*1。
而后,以这些ROI图像为输入,利用第一分辨率的分类器进行估算,可以获得ROI图像的中心点相对于目标区域中心线的位置偏移量P(1,1),也就是ROI图像中心点相对于目标区域中心线的位置关系,那么,将ROI图像的中心点通过该位置偏移量修正P(1,1)之后,就可以得到目标区域上对应的一个中心点,为了便于描述,记做目标中心点S(1,1)’。在具体的应用中,位置偏移量可以为正数或负数,可以通过将当前的ROI图像中心点与位置偏移量相加,从而获得目标中心点。若仅执行一次修正步骤,则以目标中心点S(1,1)’作为下一次中心点估算时的ROI图像的中心点S(2,1)。若执行多次修正步骤,则以目标中心点S(1,1)’作为下次修正步骤中ROI图像的中心点S(1,2),并返回并执行步骤S2021和S2022,直到执行到预设的修正的次数M,获得目标中心点S(1,M)’。
S203,进行第二次中心点的估算。
第二次中心点的估算包括:
S2031,以第二分辨率下的扫描图像为第二目标图像,在所述第二目标图像上确定出多个ROI图像,并以确定出的ROI图像为输入,利用第二分辨率的分类器估算出ROI图像的中心点相对于第二目标图像的目标区域中心线的位置偏移量P(2,m)。
S2032,通过位置偏移量P(2,m),获得ROI图像中心点S(2,m)在目标区域上对应的目标中心点S(2,m)’,m从1至M,M≥1。
该次中心点的估算中的目标图像为第二分辨率下的扫描图像,当仅执行两次中心点的估算时,第二分辨率即为扫描图像的分辨率,若执行两次以上的真心点的估算,则该第二分辨率高于第一分辨率且低于扫描图像的分辨率。
同第一次中心点估算中的描述,根据需要,在第二次中心点的估算中,可以执行一次或多次步骤S2031和S2032的修正步骤,在执行多次修正步骤时,本次修正步骤中获得的目标中心点S(2,m)’作为下次修正步骤中ROI图像的中心点S(2,m+1),M为一次估算中执行修正的总次数,可以根据需要确定修正的次数M。该次估算中,以前次估算中获得的目标中心点S(1,M)’作为首次修正步骤中ROI图像的中心点S(2,1),以这些中心点为中心可以在第二目标图像上分别确定出预设大小的ROI图像,预设大小的ROI图像的分辨率例如可以为64*64*1。
在具体的实施例中,以该中心点S(2,1)在第二目标图像上确定出多个ROI图像,以这些ROI图像为输入,利用第二分辨率的分类器进行估算,可以获得ROI图像的中心点相对于目标区域中心线的位置偏移量P(2,1),将ROI图像的中心点通过该位置偏移量修正P(2,1)之后,就可以得到目标区域上对应的一个中心点,为了便于描述,记做目标中心点S(2,1)’。若仅执行一次修正步骤,则以目标中心点S(2,1)’作为最终的目标中心点。若执行多次修正步骤,则以目标中心点S(2,1)’作为下次修正步骤中ROI图像的中心点S(2,2),并返回并执行步骤S2031和S2032,直到执行到预设的修正的次数M,获得目标中心点S(2,M)’。
在该实施例中,仅执行两次中心点的估算,第二次中心点估算中获得的目标中心点S(2,M)’则为最终的中心点,用于确定中心线。
可以理解的是,此处仅以2次估算为例进行说明,本发明并不限于此,根据不同的需要,还可以进行更多次的估算,例如3次或4次估算,之后进行的估算比前次估算采用更高分辨率的分类器,以便获得更为精确的中心点位置。在其他实施例中,可以执行更多次的中心点的估算的步骤,每次估算中,以本次估算中获得的目标中心点作为下次中心点的估算中的ROI图像的中心点S(i,1),i从1至I,I≥2,I为中心点的估算的总次数,重复进行中心点的估算的步骤,直到到达预设的估算次数I,获得目标中心点S(I,M)’。
最后,S204,根据目标中心点S(I,M)’,获得目标区域的中心线。
用于中心线确定的目标中心点为最后一次估算中的最后一次修正后获得的,这些目标中心点经过多次估算修正之后,已经较为精确的确定出了目标区域的中心点位置,通过这些中心点即可以精确地确定出中心线。具体的,可以对这些目标中心点进行直线拟合,从而获得目标区域的中心线。
为了便于理解,以下结合附图对一个具体的示例进行说明。参考图4(A)所示,为低分辨率的目标图像,目标图像中的脑部为目标区域,在目标区域的两侧边缘处设定了多个采样点,为初始的采样点,通过以这些采样点为中心点,在第一目标图像上确定出多个ROI图像。以图4(A)中的初始采样点为中心点,确定出的ROI图像为输入,利用低分辨率的分类器进行第一次中心点的估算,该次估算中仅进行一次修正步骤,参考图4(B)所示,获得此次估算的目标中心点,该目标中心点作为第二次中心点的估算时的ROI图像的中心点。在该示例中,第二次中心点的估算中也仅执行一次位置修正,以图4(B)中的目标中心点为ROI图像的中心点,确定出的ROI图像为输入,利用高分辨率的分类器进行中心点的估算之后,参考图4(C)所示,获得此次估算的目标中心点,这些目标中心点可以用于中心线的确认。参考图4(D)所示,为将图4(C)中的目标中心点进行直线拟合之后,获得的中心线。
以上对本发明实施例进行了详细的描述,该方法中,采用不同分辨率的分类器逐步确定目标中心点,兼顾中心线确定的速度和精度,可以适用于不同结构特点的解剖结构的中心线的确定,且对图像不均匀及噪声等干扰因素不敏感,中心线确定的准确性和稳定性高。参考图5所示,采用本发明实施例的方法在不同质量的扫描图像中确定的中心线示意图,其中图(A)、(B)、(C)为头部核磁共振定位的影像,这些影像的质量较低,图像的对比度低且脑中线的位置并没有明显的分界,通过本发明的方法仍可以准确地确定出脑部中心线。
此外,本发明还提供了实现上述方法的医学影像解剖结构的中心线确定装置,参考图6所示,所述装置包括:
多个不同预设分辨率的分类器,所述分类器用于以目标图像上的ROI图像为输入,估算出ROI图像的中心点相对于目标图像中目标区域的中心线的位置偏移量;
扫描图像获取单元,用于获得扫描图像,扫描图像具有目标区域;
中心点估算单元,用于进行中心点的估算,所述中心点的估算包括:
进行位置偏移量的估算以及目标中心点的确定,所述位置偏移量的估算包括:以目标图像上确定的多个ROI图像为输入,利用分类器估算出ROI图像的中心点相对于目标区域中心线的位置偏移量,目标图像为所利用的分类器的分辨率下的扫描图像;所述目标中心点的确定包括:通过位置偏移量,获得ROI图像的中心点在目标区域上对应的目标中心点;
以本次估算获得的目标中心点为下次估算中ROI图像的中心点,重复进行中心点的估算的步骤,直到到达预设的估算次数,其中,每次中心点的估算中采用的分类器的分辨率依次增大,最后一次中心点的估算中采用的分类器的分辨率为扫描图像的分辨率;中心线确定单元,用于根据最后一次估算获得的目标中心点,获得扫描图像目标区域的中心线。
进一步地,在所述中心点的估算的步骤中,执行多次进行位置偏移量的估算以及目标中心点的确定的修正步骤,且本次修正步骤中获得的目标中心点作为下次修正步骤中ROI图像的中心点。
进一步地,所述中心线确定单元中,对最后一次中心点的估算中获得的目标中心点进行直线拟合,以获得扫描图像目标区域的中心线。
进一步地,还包括:分类器确定单元,用于获得预设分辨率的扫描图像样本,扫描图像样本的目标区域上标记有中心线;在所述扫描图像样本上确定ROI图像样本,并记录ROI图像样本中心点相对于目标区域中心线的位置偏移量;通过所述ROI图像样本和位置偏移量进行分类器模型的训练,以获得所述预设分辨率的分类器。
进一步地,所述分类器模型包括特征提取器和特征分类器,所述特征提取器由深度卷积神经网络构成,所述特征分类器由全连接层构成。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的系统实施例仅仅是示意性的,其中所述作为分离部件说明的模块或单元可以是或者也可以不是物理上分开的,作为模块或单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

Claims (10)

1.一种医学影像解剖结构的中心线确定方法,其特征在于,预先设置有多个不同预设分辨率的分类器,所述分类器用于以目标图像上的ROI图像为输入,估算出ROI图像的中心点相对于目标图像中目标区域的中心线的位置偏移量;所述方法包括:
获得扫描图像,扫描图像具有目标区域;
进行中心点的估算,所述中心点的估算包括:
进行位置偏移量的估算以及目标中心点的确定,所述位置偏移量的估算包括:以目标图像上确定的多个ROI图像为输入,利用分类器估算出ROI图像的中心点相对于目标区域中心线的位置偏移量,目标图像为所利用的分类器的分辨率下的扫描图像;所述目标中心点的确定包括:通过位置偏移量,获得ROI图像的中心点在目标区域上对应的目标中心点;
以本次估算获得的目标中心点为下次估算中ROI图像的中心点,重复进行中心点的估算的步骤,直到到达预设的估算次数,其中,每次中心点的估算中采用的分类器的分辨率依次增大,最后一次中心点的估算中采用的分类器的分辨率为扫描图像的分辨率;
根据最后一次估算获得的目标中心点,获得扫描图像目标区域的中心线。
2.根据权利要求1所述的方法,其特征在于,在所述中心点的估算的步骤中,对位置偏移量估算以及目标中心点的确定的次数大于或者等于1次,且本次修正步骤中获得的目标中心点作为下次修正步骤中ROI图像的中心点。
3.根据权利要求1所述的方法,其特征在于,所述根据最后一次估算获得的目标中心点,获得扫描图像目标区域的中心线,包括:
对最后一次中心点的估算中获得的目标中心点进行直线拟合,以获得扫描图像目标区域的中心线。
4.根据权利要求1-3中任一项所述的方法,其特征在于,所述预设分辨率的分类器的确定方法包括:
获得预设分辨率的扫描图像样本,扫描图像样本的目标区域上标记有中心线;
在所述扫描图像样本上确定ROI图像样本,并记录ROI图像样本中心点相对于目标区域中心线的位置偏移量;
通过所述ROI图像样本和位置偏移量进行分类器模型的训练,以获得所述预设分辨率的分类器。
5.根据权利要求4所述的方法,其特征在于,所述分类器模型包括特征提取器和特征分类器,所述特征提取器由深度卷积神经网络构成,所述特征分类器由全连接层构成。
6.一种医学影像解剖结构的中心线确定装置,其特征在于,所述装置包括:
多个不同预设分辨率的分类器,所述分类器用于以目标图像上的ROI图像为输入,估算出ROI图像的中心点相对于目标图像中目标区域的中心线的位置偏移量;
扫描图像获取单元,用于获得扫描图像,扫描图像具有目标区域;
中心点估算单元,用于进行中心点的估算,所述中心点的估算包括:
进行位置偏移量的估算以及目标中心点的确定,所述位置偏移量的估算包括:以目标图像上确定的多个ROI图像为输入,利用分类器估算出ROI图像的中心点相对于目标区域中心线的位置偏移量,目标图像为所利用的分类器的分辨率下的扫描图像;所述目标中心点的确定包括:通过位置偏移量,获得ROI图像的中心点在目标区域上对应的目标中心点;
以本次估算获得的目标中心点为下次估算中ROI图像的中心点,重复进行中心点的估算的步骤,直到到达预设的估算次数,其中,每次中心点的估算中采用的分类器的分辨率依次增大,最后一次中心点的估算中采用的分类器的分辨率为扫描图像的分辨率;中心线确定单元,用于根据最后一次估算获得的目标中心点,获得扫描图像目标区域的中心线。
7.根据权利要求6所述的装置,其特征在于,在所述中心点的估算的步骤中,执行多次进行位置偏移量的估算以及目标中心点的确定的修正步骤,且本次修正步骤中获得的目标中心点作为下次修正步骤中ROI图像的中心点。
8.根据权利要求6所述的装置,其特征在于,所述中心线确定单元中,对最后一次中心点的估算中获得的目标中心点进行直线拟合,以获得扫描图像目标区域的中心线。
9.根据权利要求6-8中任一项所述的装置,其特征在于,还包括:分类器确定单元,用于获得预设分辨率的扫描图像样本,扫描图像样本的目标区域上标记有中心线;在所述扫描图像样本上确定ROI图像样本,并记录ROI图像样本中心点相对于目标区域中心线的位置偏移量;通过所述ROI图像样本和位置偏移量进行分类器模型的训练,以获得所述预设分辨率的分类器。
10.根据权利要求9所述的装置,其特征在于,所述分类器模型包括特征提取器和特征分类器,所述特征提取器由深度卷积神经网络构成,所述特征分类器由全连接层构成。
CN201711091221.4A 2017-11-08 2017-11-08 一种医学影像解剖结构的中心线确定方法及装置 Active CN107945226B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711091221.4A CN107945226B (zh) 2017-11-08 2017-11-08 一种医学影像解剖结构的中心线确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711091221.4A CN107945226B (zh) 2017-11-08 2017-11-08 一种医学影像解剖结构的中心线确定方法及装置

Publications (2)

Publication Number Publication Date
CN107945226A true CN107945226A (zh) 2018-04-20
CN107945226B CN107945226B (zh) 2020-04-28

Family

ID=61934501

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711091221.4A Active CN107945226B (zh) 2017-11-08 2017-11-08 一种医学影像解剖结构的中心线确定方法及装置

Country Status (1)

Country Link
CN (1) CN107945226B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110473172A (zh) * 2019-07-24 2019-11-19 上海联影智能医疗科技有限公司 医学图像解剖中线确定方法、计算机设备和存储介质
CN112785591A (zh) * 2021-03-05 2021-05-11 杭州健培科技有限公司 一种ct影像中肋骨骨折的检测与分割方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103218800A (zh) * 2011-09-27 2013-07-24 西门子公司 用于自动肋骨中心线提取的方法和系统
US20170258433A1 (en) * 2016-03-10 2017-09-14 Siemens Healthcare Gmbh Method and System for Extracting Centerline Representation of Vascular Structures in Medical Images Via Optimal Paths in Computational Flow Fields

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103218800A (zh) * 2011-09-27 2013-07-24 西门子公司 用于自动肋骨中心线提取的方法和系统
US20170258433A1 (en) * 2016-03-10 2017-09-14 Siemens Healthcare Gmbh Method and System for Extracting Centerline Representation of Vascular Structures in Medical Images Via Optimal Paths in Computational Flow Fields

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
AMOS SIRONI ET.AL: "Multiscale Centerline Detection by Learning a Scale-Space Distance Transform", 《COMPUTER VISION & PATTERN RECOGNITION》 *
SHUANGLING WANG ET.AL: "Hierarchical retinal blood vessel segmentation based on feature and ensemble learning", 《NEUROCOMPUTING》 *
朱承璋 等: "彩色眼底图像视网膜血管分割方法研究进展", 《计算机辅助设计与图形学学报》 *
秦志光 等: "基于多模态卷积神经网络的脑血管提取方法研究", 《电子科技大学学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110473172A (zh) * 2019-07-24 2019-11-19 上海联影智能医疗科技有限公司 医学图像解剖中线确定方法、计算机设备和存储介质
CN110473172B (zh) * 2019-07-24 2022-07-05 上海联影智能医疗科技有限公司 医学图像解剖中线确定方法、计算机设备和存储介质
CN112785591A (zh) * 2021-03-05 2021-05-11 杭州健培科技有限公司 一种ct影像中肋骨骨折的检测与分割方法及装置

Also Published As

Publication number Publication date
CN107945226B (zh) 2020-04-28

Similar Documents

Publication Publication Date Title
CN1765330B (zh) 超声血管测量设备和方法
CN103325143B (zh) 基于模型匹配的标记点自动注册方法
DE102010000274A1 (de) Verfahren und Vorrichtung zur automatischen Identifizierung von Bildansichten in einem 3D-Datensatz
CN106934821A (zh) 一种基于icp算法和b样条的锥形束ct和ct图像配准方法
CN107644420A (zh) 基于中心线提取的血管图像分割方法、核磁共振成像系统
CN106056537A (zh) 一种医学图像拼接方法及装置
CN109978871B (zh) 融合概率型和确定型纤维束追踪的纤维束筛选方法
CN109949349A (zh) 一种多模态三维图像的配准及融合显示方法
EP2549929B1 (en) Systems and methods for improved tractographic processing
CN102592137A (zh) 多模态图像配准方法及基于多模态图像配准的手术导航方法
CN105894508B (zh) 一种医学图像自动定位质量的评估方法
DE102014204467A1 (de) Übertragen von validierten CAD-Trainingsdaten auf geänderte MR-Kontraste
CN104182757A (zh) 获取被测目标实际覆盖面积的方法及装置
CN107945226A (zh) 一种医学影像解剖结构的中心线确定方法及装置
CN107093190A (zh) 一种基于多图谱标签融合的医学图像配准算法
JP5154751B2 (ja) 医用画像処理装置
CN101241597B (zh) 使用时间减影技术对间隔变化的视觉增强
CN101241598B (zh) 显示参考图像与浮动图像的减影图像的方法及设备
CN104915989A (zh) 基于ct影像的血管三维分割方法
CN116823663A (zh) 基于神经网络的血管图像处理方法和装置
Cao et al. Registration of medical images using an interpolated closest point transform: method and validation
Lu et al. Estimation of Transformation Parameters Between Centre‐Line Vector Road Maps and High Resolution Satellite Images
DE102018119542A1 (de) Magnetfeldverzerrungsberechnungsvorrichtung, -verfahren und -programm
US20220139003A1 (en) Methods and apparatus for mri reconstruction and data acquisition
CN110276762A (zh) 一种多b值扩散加权腹部磁共振成像的呼吸运动全自动校正方法

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

Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Applicant after: DongSoft Medical System Co., Ltd.

Address before: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Applicant before: Dongruan Medical Systems Co., Ltd., Shenyang

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant