CN102521823B - 基于弹性成像的图像处理方法及系统 - Google Patents
基于弹性成像的图像处理方法及系统 Download PDFInfo
- Publication number
- CN102521823B CN102521823B CN2011103503850A CN201110350385A CN102521823B CN 102521823 B CN102521823 B CN 102521823B CN 2011103503850 A CN2011103503850 A CN 2011103503850A CN 201110350385 A CN201110350385 A CN 201110350385A CN 102521823 B CN102521823 B CN 102521823B
- Authority
- CN
- China
- Prior art keywords
- intensity profile
- zone
- curve
- point
- consecutive point
- 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
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Image Analysis (AREA)
Abstract
一种基于弹性成像的图像处理方法及系统,包括以下步骤:获取弹性图像;从所述弹性图像中选取关注区域,并从所述关注区域提取灰度分布图;对所述灰度分布图的曲线进行平滑处理;通过对所述灰度分布图曲线上点的斜率划分得到过渡区及辨别区。通过对弹性图像的关注区域进行过渡区提取,得到图像数据的定量信息,可以从繁多复杂的信息中有目的或有针对性的获取特征。同时,也可以根据丰富的数据进行更加精确的数据分析。
Description
【技术领域】
本发明涉及图像处理技术,特别是涉及一种基于弹性成像的图像处理方法及系统。
【背景技术】
随着科学技术的进步,弹性成像技术在医疗卫生领域得到了迅速的发展。弹性成像技术能够得到组织机械力学的信息,丰富了模态成像内容。弹性成像,是对组织施加一个外部激励,在压力-形变等物理规律作用下,组织将产生一个响应,根据该响应得到弹性成像图。
弹性成像反映的是组织整体的弹性概况,目前尚没有对于弹性成像的定量后处理,或对弹性成像的局部图像区域(或者是重点关注区域)处理的相关研究,在图像数据的精确解析和利用方面尚有进一步提升的空间。
【发明内容】
基于此,有必要提供一种能够对弹性成像的局部区域进行处理且更加精确的基于弹性成像的处理方法。
另外,还有必要提供一种能够对弹性成像的局部区域进行处理且更加精确的基于弹性成像的处理系统。
一种基于弹性成像的图像处理方法,包括以下步骤:获取弹性图像;从所述弹性图像中选取关注区域,并从所述关注区域提取灰度分布图;对所述灰度分布图的曲线进行平滑处理;通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区。
优选地,所述对所述灰度分布图的曲线进行平滑处理的步骤为:按照所述灰度分布图曲线上的选取点分别提取相邻点,所述相邻点包括预设数量的上相邻点和下相邻点;计算获得上相邻点和下相邻点的灰度平均值,并分别与所述选取点的灰度值做差;判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值,是,则所述选取点的灰度值与下相邻点的灰度平均值相同。
优选地,所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤为:获取所述灰度分布图曲线上点的斜率值;将所述斜率值与预设的阈值范围进行比较,判断所述斜率值是否超过所述阈值范围,是,则所述斜率值对应的所述灰度分布图曲线的区域为对照区,否,则为过渡区。
优选地,所述阈值范围为0~0.5。
优选地,所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤之后还包括:计算所述过渡区的面积和所述对照区的面积;所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较,判断是否大于所述面积阈值,是,则所述过渡区为重点关注区。
一种基于弹性成像的图像处理系统,包括:弹性图像模块,用于加载弹性图像;提取模块,用于从所述弹性图像中选取关注区域,并从所述关注区域提取灰度分布图;平滑处理模块,用于对所述灰度分布图的曲线进行平滑处理;划分模块,用于通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区。
优选地,所属平滑处理模块包括:选取单元,用于按照所述灰度分布图曲线上的选取点分别提取相邻点,所述相邻点包括预设数量的上相邻点和下相邻点;计算单元,用于计算获得上相邻点和下相邻点的灰度平均值,并分别与所述选取点的灰度值做差;判断单元,用于判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值,是,则所述选取点的灰度值与下相邻点的灰度平均值相同。
优选地,所述划分模块包括:斜率值获取单元,用于获取所述灰度分布图曲线上点的斜率值;比较单元,用于将所述斜率值与预设的阈值范围进行比较,判断所述斜率值是否超过所述阈值范围,是,则所述斜率值对应的所述灰度分布图曲线的区域为对照区,否,则为过渡区。
优选地,所述阈值范围为0~0.5。
优选地,还包括:面积计算模块,用于计算所述过渡区的面积和所述对照区的面积;面积判断模块,用于所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较,判断是否大于所述面积阈值,是,则所述过渡区为重点关注区。
通过基于弹性成像的图像处理方法及系统在弹性图像中选取关注区域(局部区域),并提取灰度分布图,然后根据灰度分布图曲线上点的斜率信息,获得过渡区及对照区,得到丰富的图像定量数据,可以从繁多复杂的信息中有目的或有针对性的获取信息,进而获得更加精确的信息。
【附图说明】
图1为一实施例中基于弹性成像的图像处理方法的流程图;
图2为一实施例中步骤S30的具体流程图;
图3为举例说明示意图;
图4为一实施例中步骤S40的具体流程图;
图5为另一实施例的方法流程图;
图6为测试组织的实物图;
图7为测试组织的沿指示线的灰度分布度;
图8为测试组织的磁共振T1加权的图;
图9为测试组织的磁共振T1加权图沿指示线的灰度分布图;
图10为测试组织的磁共振弹性图像;
图11为测试组织的磁共振弹性图像的灰度分布图;
图12为另一组测试组织的实物图;
图13为另一组测试组织的沿指示线的灰度分布度;
图14为另一组测试组织的磁共振T1加权的图;
图15为另一组测试组织的磁共振T1加权图沿指示线的灰度分布图;
图16为另一组测试组织的磁共振弹性图像;
图17为另一组测试组织的磁共振弹性图像的灰度分布图;
图18为附图17进行平滑处理的图;
图19为基于弹性成像的图像处理系统的逻辑图;
图20为一实施例中基于弹性成像的图像处理系统的逻辑图。
【具体实施方式】
为了解决目前还没有对弹性成像定量的后处理,或对弹性成像的局部图像区域(或者是重点关注区域)处理的相关研究,在图像数据的解析方面的缺陷,提供了一种基于弹性成像的图像处理方法,详细方案如下:
结合附图1,一种基于弹性成像的图像处理方法,具体步骤如下:
步骤S10:获取弹性图像。弹性成像技术是一种较为成熟的成像技术,根据该技术可以获得弹性图像,该弹性图像可以基于MRI(磁共振)或超声弹性成像获取的图像。
步骤S20:从弹性图像中选取关注区域,并从关注区域提取灰度分布图。具体地,在弹性图像中在关注区域上,选取一“关注线”,沿该“关注线”从弹性图像中提取灰度分布图。
步骤S30:对灰度分布图的曲线进行平滑处理。对获得的灰度分布图多数为锯齿状曲线,不易于从该灰度分布图获取数据。通过对灰度分布图的曲线进行平滑处理后,易于对该关注区域的曲线进行处理,得到丰富的图像数据。
在一实施例中,结合附图2,步骤S30具体为:
步骤S31:按照灰度分布图曲线上的选取点分别提取相邻点,相邻点包括预设数量的上相邻点和下相邻点。具体地,在灰度分布曲线上任意处获得选取点,同时在该选取点相邻点处选取预设数量的上相邻点和下相邻点,该预设数量可以为2个或3个。
步骤S32:计算获得上相邻点和下相邻点的灰度平均值,并分别与选取点的灰度值做差。具体地,在灰度分布图曲线上获得预设数量的上相邻点和下相邻点的灰度值,同时根据算数平均法计算出上相邻点和下相邻点的灰度平均值。上相邻点和下相邻点的灰度值与选取点的灰度值做差,获得差值。
步骤S33:判断上相邻点与选取点的差值的绝对值是否大于下相邻点与选取点的差值的绝对值,是,则选取点的灰度值与下相邻点的灰度平均值相同。具体地,获取上相邻点的差值与选取点的差值的绝对值,以及下相邻点的差值与选取点的差值的绝对值,并判断绝对值大小。是,则选取点的灰度值与下相邻点的灰度平均值相同;否,则选取点的灰度值与上相邻点的灰度平均值相同。
结合具体例子说明,选取点的相邻点的预设数量为2个,结合附图3,具体如下:
在灰度分布图曲线上随机的获得选取点,具体的:选取点为A,其2个上相邻点为A1和A2,2个下相邻点为A-1和A-2;所对应的上相邻点的灰度值为F1和F2,2个下相邻点的灰度值为F-1和F-2,A的灰度值为F。上相邻点对应的灰度平均值为FS=(F1+F2)/2;下相邻点对应的灰度平均值为FX=(F-1+F-2)/2。若F与FS的差值的绝对值大于F与FX的差值的绝对值,则选取点A所对应的F变换为FX,即与下相邻点的FX值相同。
步骤S40:通过对灰度分布图曲线上点的斜率划分得到过渡区及对照区。具体地,获取灰度分布图曲线上点的斜率值,根据斜率值划分区域数据分析区域,即过渡区及对照区,进一步的从弹性图像中的关注区域获取丰富的图像数据。
在一实施例中,结合附图4,步骤S40具体为:
步骤S41:获取灰度分布图曲线上点的斜率值。具体地,在灰度分布曲线的点上获取对应的点的斜率。
步骤S42:将斜率值与预设的阈值范围进行比较,判断斜率值是否超过阈值范围,是,则斜率值对应的灰度分布图曲线的区域为对照区,否,则为过渡区。具体地,将斜率值与预设的阈值范围进行比较,在其它实施例中,该阈值范围为0~0.5,若超过该阈值范围,则该斜率值对应的灰度分布图曲线与横坐标(距离)所包围的区域为对照区;否,则为过渡区。
通过本方案在弹性图像中选取关注区域,并提取灰度分布图,然后根据灰度分布图曲线上点的斜率信息,获得过渡区域及对照区域。可见,对弹性图像的关注区域进行图像处理,得到了丰富的图像数据,而不仅仅是简单的提供弹性图像。
在另一实施例中,结合附图5,还包括如下步骤:
步骤S50:计算过渡区的面积和对照区的面积。具体地,采用微积分方法、填补发或估算法分别得到过渡区的面积和对照区的面积。
步骤S60:过渡区的面积与对照区的面积的比值与预设的面积阈值进行比较,判断是否大于面积阈值,是,则过渡区为重点关注区。具体地,根据渡区的面积与对照区的面积的比值,若超过面积阈值,例如为1/3,则确认该过渡区为重点关注区域,为进一步获取信息提供准确的指引;若没有超过面积阈值,则为非重点关注区。
结合具体附图6~11的实验数据,说明本申请的特点,具体地:
附图6为测试组织的实物图,附图7为该测试组织的沿指示线的灰度分布度;
附图8为测试组织的磁共振T1加权的图,附图9为测试组织的磁共振T1加权图沿指示线的灰度分布图;
附图10为测试组织的磁共振弹性图像,附图11为测试组织的磁共振弹性图像的灰度分布图。
清楚的看到实物图的灰度分布图为锯齿状,相应的磁共振图的灰度分布图为小幅度波动状,趋于均一,进一步提取组织特征信息的难度较大;而磁共振弹性图像的灰度分布图为有一定斜率趋势的曲线,信息丰富。
同时,实物图的灰度分布图、磁共振的灰度分布图以及磁共振弹性图像的灰度分布图都不相一致,而实物图和磁共振的灰度分布图仅是从图像本身反映对实物或图像的认知,而无法获取丰富的数据。磁共振弹性图像的灰度分布图从另一个角度进行分析,通过对该灰度分布图的处理获取丰富的数据。
结合另一组附图12~17的实验数据,进一步说明本申请的优势,具体地:
附图12为另一组测试组织的实物图,附图13为另一组测试组织的沿指示线的灰度分布度;
附图14为另一组测试组织的磁共振T1加权的图,附图15为另一组测试组织的磁共振T1加权图沿指示线的灰度分布图;
附图16为另一组测试组织的磁共振弹性图像,附图17为另一组测试组织的磁共振弹性图像的灰度分布图。
实物图的灰度分布图和磁共振常规图像的灰度分布图与磁共振弹性图像的灰度分布相差较大。本申请方案从弹性成像中所获取的信息较为丰富准确,可对关注区域进一步分析。
现例举一具体方案,结合附图18,划分过渡区和对照区,以及辨别重点关注区和非重点关注区。
附图18为附图17进行平滑处理的图,由此可见,能获得较为光滑的曲线图。根据斜率值在0~0.5范围内所对应的灰度分布图曲线的区域为对照区,不在该范围内的为过渡区。
同时,每个区域所对应的面积为S1、S2和S3。可见,S1/S2和S3/S2大于1/3,则属于重点关注区。
基于上述基于弹性成像的图像处理方法,还有必要建立基于弹性成像的图像处理系统,具体如下:
结合附图19,一种基于弹性成像的图像处理系统,具体步骤如下:
弹性图像模块10,用于加载弹性图像。弹性成像技术是一种较为成熟的成像技术,根据该技术可以获得弹性图像并加载,该弹性图像可以是磁共振弹性成像或超声弹性图像。
提取模块20,用于从弹性图像中选取关注区域,并从关注区域提取灰度分布图。具体地,在弹性图像中在关注区域上,选取一“关注线”,沿该“关注线”从弹性图像中提取灰度分布图。
平滑处理模块30,用于对灰度分布图的曲线进行平滑处理。对获得的灰度分布图多数为锯齿状曲线,不易于从该灰度分布图获取数据。通过对灰度分布图的曲线进行平滑处理后,易于对该关注区域的曲线进行处理,得到丰富的图像数据。
在一实施例中,平滑处理模块30包括:
选取单元31,用于按照灰度分布图曲线上的选取点分别提取相邻点,相邻点包括预设数量的上相邻点和下相邻点。具体地,在灰度分布曲线上任意处获得选取点,同时在该选取点相邻点处选取预设数量的上相邻点和下相邻点,该预设数量可以为2个或3个。
计算单元32,用于计算获得上相邻点和下相邻点的灰度平均值,并分别与选取点的灰度值做差。具体地,在灰度分布图曲线上获得预设数量的上相邻点和下相邻点的灰度值,同时根据算数平均法计算出上相邻点和下相邻点的灰度平均值。上相邻点和下相邻点的灰度值与选取点的灰度值做差,获得差值。
判断单元33,用于判断上相邻点与选取点的差值的绝对值是否大于下相邻点与选取点的差值的绝对值,是,则选取点的灰度值与下相邻点的灰度平均值相同。具体地,获取上相邻点的差值与选取点的差值的绝对值,以及下相邻点的差值与选取点的差值的绝对值,并判断绝对值大小。是,则选取点的灰度值与下相邻点的灰度平均值相同;否,则选取点的灰度值与上相邻点的灰度平均值相同。
结合具体例子说明,选取点的相邻点的预设数量为2个,具体如下:
在灰度分布图曲线上随机的获得选取点,具体的:选取点为A,其2个上相邻点为A1和A2,2个下相邻点为A-1和A-2;所对应的上相邻点的灰度值为F1和F2,2个下相邻点的灰度值为F-1和F-2,A的灰度值为F。上相邻点对应的灰度平均值为FS=(F1+F2)/2;下相邻点对应的灰度平均值为FX=(F-1+F-2)/2。若F与FS的差值的绝对值大于F与FX的差值的绝对值,则选取点A所对应的F变换为FX,即与下相邻点的FX值相同。
划分模块40,用于通过对灰度分布图曲线上点的斜率划分得到过渡区及对照区。具体地,获取灰度分布图曲线上点的斜率值,根据斜率值划分区域数据分析区域,即过渡区及对照区,进一步的从弹性图像中的关注区域获取丰富的图像数据。
在一实施例中,划分模块40包括:
斜率值获取单元41,获取灰度分布图曲线上点的斜率值。具体地,在灰度分布曲线的点上获取对应的点的斜率。
比较单元42,用于将斜率值与预设的阈值范围进行比较,判断斜率值是否超过阈值范围,是,则斜率值对应的灰度分布图曲线的区域为对照区,否,则为过渡区。具体地,将斜率值与预设的阈值范围进行比较,在其它实施例中,该阈值范围为0~0.5,若超过该阈值范围,则该斜率值对应的灰度分布图曲线与横坐标(距离)所包围的区域为对照区;否,则为过渡区。
通过本方案在弹性图像中选取关注区域,并提取灰度分布图,然后根据灰度分布图曲线上点的斜率信息,获得过渡区域及对照区域。可见,对弹性图像的关注区域进行图像处理,得到了丰富的图像数据,而不仅仅是简单的提供弹性图像。
在另一实施例中,结合附图20,还包括:
面积计算模块,用于计算过渡区的面积和对照区的面积。具体地,采用微积分方法、填补发或估算法分别得到过渡区的面积和对照区的面积。
面积判断模块,用于过渡区的面积与对照区的面积的比值与预设的面积阈值进行比较,判断是否大于面积阈值,是,则关注区域为重点关注区。具体地,根据渡区的面积与对照区的面积的比值,若超过面积阈值,例如为1/3,则确认该关注区域为重点关注区域,为进一步获取信息提供准确的指引;若没有超过面积阈值,则为非重点关注区。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (8)
1.一种基于弹性成像的图像处理方法,其特征在于,包括以下步骤:
获取弹性图像;
从所述弹性图像中选取关注区域,并从所述关注区域提取灰度分布图;
对所述灰度分布图的曲线进行平滑处理;
通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区;
所述对所述灰度分布图的曲线进行平滑处理的步骤为:
按照所述灰度分布图曲线上的选取点分别提取相邻点,所述相邻点包括预设数量的上相邻点和下相邻点;
计算获得上相邻点和下相邻点的灰度平均值,并分别与所述选取点的灰度值做差;
判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值,是,则所述选取点的灰度值与下相邻点的灰度平均值相同。
2.根据权利要求1所述的基于弹性成像的图像处理方法,其特征在于,所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤为:
获取所述灰度分布图曲线上点的斜率值;
将所述斜率值与预设的阈值范围进行比较,判断所述斜率值是否超过所述阈值范围,是,则所述斜率值对应的所述灰度分布图曲线的区域为对照区,否,则为过渡区。
3.根据权利要求2所述的基于弹性成像的图像处理方法,其特征在于,所述阈值范围为0~0.5。
4.根据权利要求2所述的基于弹性成像的图像处理方法,其特征在于,所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤之后还包括:
计算所述过渡区的面积和所述对照区的面积;
所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较,判断是否大于所述面积阈值,是,则所述过渡区为重点关注区。
5.一种基于弹性成像的图像处理系统,其特征在于,包括:
弹性图像模块,用于加载弹性图像;
提取模块,用于从所述弹性图像中选取关注区域,并从所述关注区域提取灰度分布图;
平滑处理模块,用于对所述灰度分布图的曲线进行平滑处理;
划分模块,用于通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区;
所述平滑处理模块包括:
选取单元,用于按照所述灰度分布图曲线上的选取点分别提取相邻点,所述相邻点包括预设数量的上相邻点和下相邻点;
计算单元,用于计算获得上相邻点和下相邻点的灰度平均值,并分别与所述选取点的灰度值做差;
判断单元,用于判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值,是,则所述选取点的灰度值与下相邻点的灰度平均值相同。
6.根据权利要求5所述的基于弹性成像的图像处理系统,其特征在于,所述划分模块包括:
斜率值获取单元,用于获取所述灰度分布图曲线上点的斜率值;
比较单元,用于将所述斜率值与预设的阈值范围进行比较,判断所述斜率值是否超过所述阈值范围,是,则所述斜率值对应的所述灰度分布图曲线的区域为对照区,否,则为过渡区。
7.根据权利要求6所述的基于弹性成像的图像处理系统,其特征在于,所述阈值范围为0~0.5。
8.根据权利要求5所述的基于弹性成像的图像处理系统,其特征在于,还包括:
面积计算模块,用于计算所述过渡区的面积和所述对照区的面积;
面积判断模块,用于所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较,判断是否大于所述面积阈值,是,则所述过渡区为重点关注区。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103503850A CN102521823B (zh) | 2011-11-08 | 2011-11-08 | 基于弹性成像的图像处理方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103503850A CN102521823B (zh) | 2011-11-08 | 2011-11-08 | 基于弹性成像的图像处理方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102521823A CN102521823A (zh) | 2012-06-27 |
CN102521823B true CN102521823B (zh) | 2013-12-11 |
Family
ID=46292727
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011103503850A Active CN102521823B (zh) | 2011-11-08 | 2011-11-08 | 基于弹性成像的图像处理方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102521823B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109919918A (zh) * | 2019-02-21 | 2019-06-21 | 清华大学 | 弹性成像的控制方法及装置、计算机设备及可读存储介质 |
CN113030816B (zh) * | 2021-04-15 | 2022-03-11 | 上海交通大学 | 提高磁共振弹性成像模量计算分辨率的方法、系统及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013503A (zh) * | 2007-01-26 | 2007-08-08 | 清华大学 | 一种医学图像中腹部器官分割方法 |
CN101226155A (zh) * | 2007-12-21 | 2008-07-23 | 中国人民解放军第八一医院 | 智能化肺癌早期细胞病理图像识别处理方法 |
CN101276218A (zh) * | 2008-04-29 | 2008-10-01 | 上海交通大学 | 针对高实时性图像采集与处理的机器视觉系统 |
-
2011
- 2011-11-08 CN CN2011103503850A patent/CN102521823B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013503A (zh) * | 2007-01-26 | 2007-08-08 | 清华大学 | 一种医学图像中腹部器官分割方法 |
CN101226155A (zh) * | 2007-12-21 | 2008-07-23 | 中国人民解放军第八一医院 | 智能化肺癌早期细胞病理图像识别处理方法 |
CN101276218A (zh) * | 2008-04-29 | 2008-10-01 | 上海交通大学 | 针对高实时性图像采集与处理的机器视觉系统 |
Non-Patent Citations (4)
Title |
---|
冠脉造影图像的血管识别方法研究;李春林;《中国优秀硕士学位论文全文数据库 信息科技辑》;20110531(第5期);1-50 * |
基于直方图的感兴趣区图像边缘检测;张立保,李冬玲,余先川,王鹏飞,蔡蕾;《强激光与粒子束》;20100831;第22卷(第8期);1847-1850 * |
张立保,李冬玲,余先川,王鹏飞,蔡蕾.基于直方图的感兴趣区图像边缘检测.《强激光与粒子束》.2010,第22卷(第8期),1847-1850. |
李春林.冠脉造影图像的血管识别方法研究.《中国优秀硕士学位论文全文数据库 信息科技辑》.2011,(第5期),1-50. |
Also Published As
Publication number | Publication date |
---|---|
CN102521823A (zh) | 2012-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9445744B2 (en) | Methods, systems, and devices for spine centrum extraction and intervertebral disk dividing | |
CN104574298B (zh) | 一种基于互信息的多b值扩散权重图像的降噪方法 | |
Alperin et al. | PUBS: pulsatility‐based segmentation of lumens conducting non‐steady flow | |
CN103985099B (zh) | 一种弥散张量磁共振图像张量域非局部均值去噪方法 | |
Brecheisen et al. | Parameter sensitivity visualization for DTI fiber tracking | |
CN101669828A (zh) | 基于pet/ct图像纹理特征的肺部恶性肿瘤与良性结节检测系统 | |
CN107292884A (zh) | 一种识别mri图像中水肿和血肿的方法及装置 | |
CN103597370A (zh) | 空间编码的相位对比磁共振成像 | |
CN110391016B (zh) | 一种动态对比增强磁共振图像的分析方法 | |
CN102274051A (zh) | 一种超声图像膀胱容积自动测量方法及系统 | |
Xing et al. | A multi-scale contrast-based image quality assessment model for multi-exposure image fusion | |
CN1918481A (zh) | 高角分辨扩散加权mri | |
KR20090088404A (ko) | 의료 이미징 시스템 | |
CN108550145A (zh) | 一种sar图像质量评估方法和装置 | |
CN102521823B (zh) | 基于弹性成像的图像处理方法及系统 | |
Martin et al. | Contrast-to-noise ratio analysis of microscopic diffusion anisotropy indices in q-space trajectory imaging | |
CN102436551A (zh) | 基于目标跟踪的计算机辅助胃癌诊断方法 | |
CN103985110A (zh) | 数字图像中伪影或非均匀亮度的定位及显示方法 | |
CN103810712A (zh) | 一种能谱ct图像质量的评价方法 | |
US20100274119A1 (en) | Blood flow dynamic analysis apparatus, method of blood flow dynamic analysis, and magnetic resonance imaging system | |
CN110458834A (zh) | 一种乳腺肿瘤图像处理系统、方法及装置 | |
Lu et al. | An improved method for CT/MRI image fusion on bandelets transform domain | |
CN109242866A (zh) | 基于扩散磁共振图像的乳腺肿瘤自动辅助检测系统 | |
Konukoglu et al. | Monitoring slowly evolving tumors | |
CN102193076B (zh) | 用于确定背景相位变化的方法和磁共振装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |