专利名称：基于弹性成像的图像处理方法及系统的制作方法
基于弹性成像的图像处理方法及系统技术领域：
本发明涉及图像处理技术，特别是涉及一种基于弹性成像的图像处理方法及系统。背景技术：
随着科学技术的进步，弹性成像技术在医疗卫生领域得到了迅速的发展。弹性成像技术能够得到组织机械力学的信息，丰富了模态成像内容。弹性成像，是对组织施加一个外部激励，在压力-形变等物理规律作用下，组织将产生一个响应，根据该响应得到弹性成像图。
弹性成像反映的是组织整体的弹性概况，目前尚没有对于弹性成像的定量后处理，或对弹性成像的局部图像区域(或者是重点关注区域)处理的相关研究，在图像数据的精确解析和利用方面尚有进一步提升的空间。
发明内容
基于此，有必要提供一种能够对弹性成像的局部区域进行处理且更加精确的基于弹性成像的处理方法。
另外，还有必要提供一种能够对弹性成像的局部区域进行处理且更加精确的基于弹性成像的处理系统。
一种基于弹性成像的图像处理方法，包括以下步骤获取弹性图像；从所述弹性图像中选取关注区域，并从所述关注区域提取灰度分布图；对所述灰度分布图的曲线进行平滑处理；通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区。
优选地，所述对所述灰度分布图的曲线进行平滑处理的步骤为按照所述灰度分布图曲线上的选取点分别提取相邻点，所述相邻点包括预设数量的上相邻点和下相邻点； 计算获得上相邻点和下相邻点的灰度平均值，并分别与所述选取点的灰度值做差；判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值，是，则所述选取点的灰度值与下相邻点的灰度平均值相同。
优选地，所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤为获取所述灰度分布图曲线上点的斜率值；将所述斜率值与预设的阈值范围进行比较，判断所述斜率值是否超过所述阈值范围，是，则所述斜率值对应的所述灰度分布图曲线的区域为对照区，否，则为过渡区。
优选地，所述阈值范围为0 0.5。
优选地，所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤之后还包括计算所述过渡区的面积和所述对照区的面积；所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较，判断是否大于所述面积阈值，是，则所述过渡区为重点关注区。
一种基于弹性成像的图像处理系统，包括弹性图像模块，用于加载弹性图像；提取模块，用于从所述弹性图像中选取关注区域，并从所述关注区域提取灰度分布图；平滑处理模块，用于对所述灰度分布图的曲线进行平滑处理；划分模块，用于通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区。
优选地，所属平滑处理模块包括选取单元，用于按照所述灰度分布图曲线上的选取点分别提取相邻点，所述相邻点包括预设数量的上相邻点和下相邻点；计算单元，用于计算获得上相邻点和下相邻点的灰度平均值，并分别与所述选取点的灰度值做差；判断单元， 用于判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值，是，则所述选取点的灰度值与下相邻点的灰度平均值相同。
优选地，所述划分模块包括斜率值获取单元，用于获取所述灰度分布图曲线上点的斜率值；比较单元，用于将所述斜率值与预设的阈值范围进行比较，判断所述斜率值是否超过所述阈值范围，是，则所述斜率值对应的所述灰度分布图曲线的区域为对照区，否，则为过渡区。
优选地，所述阈值范围为0 0.5。
优选地，还包括面积计算模块，用于计算所述过渡区的面积和所述对照区的面积；面积判断模块，用于所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较，判断是否大于所述面积阈值，是，则所述过渡区为重点关注区。
通过基于弹性成像的图像处理方法及系统在弹性图像中选取关注区域(局部区域)，并提取灰度分布图，然后根据灰度分布图曲线上点的斜率信息，获得过渡区及对照区， 得到丰富的图像定量数据，可以从繁多复杂的信息中有目的或有针对性的获取信息，进而获得更加精确的信息。
图1为一实施例中基于弹性成像的图像处理方法的流程图2为一实施例中步骤S30的具体流程图3为举例说明示意图4为一实施例中步骤S40的具体流程图5为另一实施例的方法流程图6为测试组织的实物图7为测试组织的沿指示线的灰度分布度；
图8为测试组织的磁共振Tl加权的图9为测试组织的磁共振Tl加权图沿指示线的灰度分布图10为测试组织的磁共振弹性图像；
图11为测试组织的磁共振弹性图像的灰度分布图12为另一组测试组织的实物图13为另一组测试组织的沿指示线的灰度分布度；
图14为另一组测试组织的磁共振Tl加权的图15为另一组测试组织的磁共振Tl加权图沿指示线的灰度分布图16为另一组测试组织的磁共振弹性图像；
图17为另一组测试组织的磁共振弹性图像的灰度分布图18为附图17进行平滑处理的图19为基于弹性成像的图像处理系统的逻辑图20为一实施例中基于弹性成像的图像处理系统的逻辑图。
具体实施方式
为了解决目前还没有对弹性成像定量的后处理，或对弹性成像的局部图像区域 (或者是重点关注区域)处理的相关研究，在图像数据的解析方面的缺陷，提供了一种基于弹性成像的图像处理方法，详细方案如下
结合附图1，一种基于弹性成像的图像处理方法，具体步骤如下
步骤SlO 获取弹性图像。弹性成像技术是一种较为成熟的成像技术，根据该技术可以获得弹性图像，该弹性图像可以基于MRI (磁共振)或超声弹性成像获取的图像。
步骤S20 从弹性图像中选取关注区域，并从关注区域提取灰度分布图。具体地， 在弹性图像中在关注区域上，选取一“关注线”，沿该“关注线”从弹性图像中提取灰度分布图。
步骤S30 对灰度分布图的曲线进行平滑处理。对获得的灰度分布图多数为锯齿状曲线，不易于从该灰度分布图获取数据。通过对灰度分布图的曲线进行平滑处理后，易于对该关注区域的曲线进行处理，得到丰富的图像数据。
在一实施例中，结合附图2，步骤S30具体为
步骤S31 按照灰度分布图曲线上的选取点分别提取相邻点，相邻点包括预设数量的上相邻点和下相邻点。具体地，在灰度分布曲线上任意处获得选取点，同时在该选取点相邻点处选取预设数量的上相邻点和下相邻点，该预设数量可以为2个或3个。
步骤S32 计算获得上相邻点和下相邻点的灰度平均值，并分别与选取点的灰度值做差。具体地，在灰度分布图曲线上获得预设数量的上相邻点和下相邻点的灰度值，同时根据算数平均法计算出上相邻点和下相邻点的灰度平均值。上相邻点和下相邻点的灰度值与选取点的灰度值做差，获得差值。
步骤S33 判断上相邻点与选取点的差值的绝对值是否大于下相邻点与选取点的差值的绝对值，是，则选取点的灰度值与下相邻点的灰度平均值相同。具体地，获取上相邻点的差值与选取点的差值的绝对值，以及下相邻点的差值与选取点的差值的绝对值，并判断绝对值大小。是，则选取点的灰度值与下相邻点的灰度平均值相同；否，则选取点的灰度值与上相邻点的灰度平均值相同。
结合具体例子说明，选取点的相邻点的预设数量为2个，结合附图3，具体如下
在灰度分布图曲线上随机的获得选取点，具体的选取点为A，其2个上相邻点为 A1和A2, 2个下相邻点为L1和A_2 ；所对应的上相邻点的灰度值为F1和F2, 2个下相邻点的灰度值为F—和F_2，A的灰度值为F。上相邻点对应的灰度平均值为Fs = (FJF2)/2 ；下相邻点对应的灰度平均值为I7x = (F_i+F_2)/2。若F与Fs的差值的绝对值大于F与!^x的差值的绝对值，则选取点A所对应的F变换为Fx，即与下相邻点的!^x值相同。
步骤S40 通过对灰度分布图曲线上点的斜率划分得到过渡区及对照区。具体地， 获取灰度分布图曲线上点的斜率值，根据斜率值划分区域数据分析区域，即过渡区及对照区，进一步的从弹性图像中的关注区域获取丰富的图像数据。
在一实施例中，结合附图4，步骤S40具体为
步骤S41 获取灰度分布图曲线上点的斜率值。具体地，在灰度分布曲线的点上获取对应的点的斜率。
步骤S42 将斜率值与预设的阈值范围进行比较，判断斜率值是否超过阈值范围， 是，则斜率值对应的灰度分布图曲线的区域为对照区，否，则为过渡区。具体地，将斜率值与预设的阈值范围进行比较，在其它实施例中，该阈值范围为0 0. 5，若超过该阈值范围，则该斜率值对应的灰度分布图曲线与横坐标(距离)所包围的区域为对照区；否，则为过渡区。
通过本方案在弹性图像中选取关注区域，并提取灰度分布图，然后根据灰度分布图曲线上点的斜率信息，获得过渡区域及对照区域。可见，对弹性图像的关注区域进行图像处理，得到了丰富的图像数据，而不仅仅是简单的提供弹性图像。
在另一实施例中，结合附图5，还包括如下步骤
步骤S50 计算过渡区的面积和对照区的面积。具体地，采用微积分方法、填补发或估算法分别得到过渡区的面积和对照区的面积。
步骤S60 过渡区的面积与对照区的面积的比值与预设的面积阈值进行比较，判断是否大于面积阈值，是，则过渡区为重点关注区。具体地，根据渡区的面积与对照区的面积的比值，若超过面积阈值，例如为1/3，则确认该过渡区为重点关注区域，为进一步获取信息提供准确的指引；若没有超过面积阈值，则为非重点关注区。
结合具体附图6 11的实验数据，说明本申请的特点，具体地
附图6为测试组织的实物图，附图7为该测试组织的沿指示线的灰度分布度；
附图8为测试组织的磁共振Tl加权的图，附图9为测试组织的磁共振Tl加权图沿指示线的灰度分布附图10为测试组织的磁共振弹性图像，附图11为测试组织的磁共振弹性图像的灰度分布图。
清楚的看到实物图的灰度分布图为锯齿状，相应的磁共振图的灰度分布图为小幅度波动状，趋于均一，进一步提取组织特征信息的难度较大；而磁共振弹性图像的灰度分布图为有一定斜率趋势的曲线，信息丰富。
同时，实物图的灰度分布图、磁共振的灰度分布图以及磁共振弹性图像的灰度分布图都不相一致，而实物图和磁共振的灰度分布图仅是从图像本身反映对实物或图像的认知，而无法获取丰富的数据。磁共振弹性图像的灰度分布图从另一个角度进行分析，通过对该灰度分布图的处理获取丰富的数据。
结合另一组附图12 17的实验数据，进一步说明本申请的优势，具体地
附图12为另一组测试组织的实物图，附图13为另一组测试组织的沿指示线的灰度分布度；
附图14为另一组测试组织的磁共振Tl加权的图，附图15为另一组测试组织的磁共振Tl加权图沿指示线的灰度分布附图16为另一组测试组织的磁共振弹性图像，附图17为另一组测试组织的磁共振弹性图像的灰度分布图。
实物图的灰度分布图和磁共振常规图像的灰度分布图与磁共振弹性图像的灰度分布相差较大。本申请方案从弹性成像中所获取的信息较为丰富准确，可对关注区域进一步分析。
现例举一具体方案，结合附图18，划分过渡区和对照区，以及辨别重点关注区和非重点关注区。
附图18为附图17进行平滑处理的图，由此可见，能获得较为光滑的曲线图。根据斜率值在0 0. 5范围内所对应的灰度分布图曲线的区域为对照区，不在该范围内的为过渡区。
同时，每个区域所对应的面积为Si、S2和S3。可见，S1/S2和S3/S2大于1/3，则属于重点关注区。
基于上述基于弹性成像的图像处理方法，还有必要建立基于弹性成像的图像处理系统，具体如下
结合附图19，一种基于弹性成像的图像处理系统，具体步骤如下
弹性图像模块10，用于加载弹性图像。弹性成像技术是一种较为成熟的成像技术，根据该技术可以获得弹性图像并加载，该弹性图像可以是磁共振弹性成像或超声弹性图像。
提取模块20，用于从弹性图像中选取关注区域，并从关注区域提取灰度分布图。具体地，在弹性图像中在关注区域上，选取一“关注线”，沿该“关注线”从弹性图像中提取灰度分布图。
平滑处理模块30，用于对灰度分布图的曲线进行平滑处理。对获得的灰度分布图多数为锯齿状曲线，不易于从该灰度分布图获取数据。通过对灰度分布图的曲线进行平滑处理后，易于对该关注区域的曲线进行处理，得到丰富的图像数据。
在一实施例中，平滑处理模块30包括
选取单元31，用于按照灰度分布图曲线上的选取点分别提取相邻点，相邻点包括预设数量的上相邻点和下相邻点。具体地，在灰度分布曲线上任意处获得选取点，同时在该选取点相邻点处选取预设数量的上相邻点和下相邻点，该预设数量可以为2个或3个。
计算单元32，用于计算获得上相邻点和下相邻点的灰度平均值，并分别与选取点的灰度值做差。具体地，在灰度分布图曲线上获得预设数量的上相邻点和下相邻点的灰度值，同时根据算数平均法计算出上相邻点和下相邻点的灰度平均值。上相邻点和下相邻点的灰度值与选取点的灰度值做差，获得差值。
判断单元33，用于判断上相邻点与选取点的差值的绝对值是否大于下相邻点与选取点的差值的绝对值，是，则选取点的灰度值与下相邻点的灰度平均值相同。具体地，获取上相邻点的差值与选取点的差值的绝对值，以及下相邻点的差值与选取点的差值的绝对值，并判断绝对值大小。是，则选取点的灰度值与下相邻点的灰度平均值相同；否，则选取点的灰度值与上相邻点的灰度平均值相同。
结合具体例子说明，选取点的相邻点的预设数量为2个，具体如下
在灰度分布图曲线上随机的获得选取点，具体的选取点为A，其2个上相邻点为 A1和A2, 2个下相邻点为L1和A_2 ；所对应的上相邻点的灰度值为F1和F2, 2个下相邻点的灰度值为F—和F_2，A的灰度值为F。上相邻点对应的灰度平均值为Fs = (FJF2)/2 ；下相邻点对应的灰度平均值为I7x = (F_i+F_2)/2。若F与Fs的差值的绝对值大于F与!^x的差值的绝对值，则选取点A所对应的F变换为Fx，即与下相邻点的!^x值相同。
划分模块40，用于通过对灰度分布图曲线上点的斜率划分得到过渡区及对照区。 具体地，获取灰度分布图曲线上点的斜率值，根据斜率值划分区域数据分析区域，即过渡区及对照区，进一步的从弹性图像中的关注区域获取丰富的图像数据。
在一实施例中，划分模块40包括
斜率值获取单元41，获取灰度分布图曲线上点的斜率值。具体地，在灰度分布曲线的点上获取对应的点的斜率。
比较单元42，用于将斜率值与预设的阈值范围进行比较，判断斜率值是否超过阈值范围，是，则斜率值对应的灰度分布图曲线的区域为对照区，否，则为过渡区。具体地，将斜率值与预设的阈值范围进行比较，在其它实施例中，该阈值范围为O 0. 5，若超过该阈值范围，则该斜率值对应的灰度分布图曲线与横坐标(距离)所包围的区域为对照区；否， 则为过渡区。
通过本方案在弹性图像中选取关注区域，并提取灰度分布图，然后根据灰度分布图曲线上点的斜率信息，获得过渡区域及对照区域。可见，对弹性图像的关注区域进行图像处理，得到了丰富的图像数据，而不仅仅是简单的提供弹性图像。
在另一实施例中，结合附图20，还包括
面积计算模块，用于计算过渡区的面积和对照区的面积。具体地，采用微积分方法、填补发或估算法分别得到过渡区的面积和对照区的面积。
面积判断模块，用于过渡区的面积与对照区的面积的比值与预设的面积阈值进行比较，判断是否大于面积阈值，是，则关注区域为重点关注区。具体地，根据渡区的面积与对照区的面积的比值，若超过面积阈值，例如为1/3，则确认该关注区域为重点关注区域，为进一步获取信息提供准确的指引；若没有超过面积阈值，则为非重点关注区。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
权利要求
1.一种基于弹性成像的图像处理方法，其特征在于，包括以下步骤获取弹性图像；从所述弹性图像中选取关注区域，并从所述关注区域提取灰度分布图；对所述灰度分布图的曲线进行平滑处理；通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区。
2.根据权利要求1所述的基于弹性成像的图像处理方法，其特征在于，所述对所述灰度分布图的曲线进行平滑处理的步骤为按照所述灰度分布图曲线上的选取点分别提取相邻点，所述相邻点包括预设数量的上相邻点和下相邻点；计算获得上相邻点和下相邻点的灰度平均值，并分别与所述选取点的灰度值做差；判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值，是，则所述选取点的灰度值与下相邻点的灰度平均值相同。
3.根据权利要求2所述的基于弹性成像的图像处理方法，其特征在于，所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤为获取所述灰度分布图曲线上点的斜率值；将所述斜率值与预设的阈值范围进行比较，判断所述斜率值是否超过所述阈值范围， 是，则所述斜率值对应的所述灰度分布图曲线的区域为对照区，否，则为过渡区。
4.根据权利要求3所述的基于弹性成像的图像处理方法，其特征在于，所述阈值范围为 0. 5。
5.根据权利要求3所述的基于弹性成像的图像处理方法，其特征在于，所述通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区的步骤之后还包括计算所述过渡区的面积和所述对照区的面积；所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较，判断是否大于所述面积阈值，是，则所述过渡区为重点关注区。
6.一种基于弹性成像的图像处理系统，其特征在于，包括弹性图像模块，用于加载弹性图像；提取模块，用于从所述弹性图像中选取关注区域，并从所述关注区域提取灰度分布图；平滑处理模块，用于对所述灰度分布图的曲线进行平滑处理；划分模块，用于通过对所述灰度分布图曲线上点的斜率划分得到过渡区及对照区。
7.根据权利要求6所述的基于弹性成像的图像处理系统，其特征在于，所属平滑处理模块包括选取单元，用于按照所述灰度分布图曲线上的选取点分别提取相邻点，所述相邻点包括预设数量的上相邻点和下相邻点；计算单元，用于计算获得上相邻点和下相邻点的灰度平均值，并分别与所述选取点的灰度值做差；判断单元，用于判断所述上相邻点与所述选取点的差值的绝对值是否大于所述下相邻点与所述选取点的差值的绝对值，是，则所述选取点的灰度值与下相邻点的灰度平均值相同。
8.根据权利要求7所述的基于弹性成像的图像处理系统，其特征在于，所述划分模块包括斜率值获取单元，用于获取所述灰度分布图曲线上点的斜率值； 比较单元，用于将所述斜率值与预设的阈值范围进行比较，判断所述斜率值是否超过所述阈值范围，是，则所述斜率值对应的所述灰度分布图曲线的区域为对照区，否，则为过渡区。
9.根据权利要求8所述的基于弹性成像的图像处理系统，其特征在于，所述阈值范围为 0. 5。
10.根据权利要求7所述的基于弹性成像的图像处理系统，其特征在于，还包括 面积计算模块，用于计算所述过渡区的面积和所述对照区的面积；面积判断模块，用于所述过渡区的面积与所述对照区的面积的比值与预设的面积阈值进行比较，判断是否大于所述面积阈值，是，则所述过渡区为重点关注区。
全文摘要
一种基于弹性成像的图像处理方法及系统，包括以下步骤获取弹性图像；从所述弹性图像中选取关注区域，并从所述关注区域提取灰度分布图；对所述灰度分布图的曲线进行平滑处理；通过对所述灰度分布图曲线上点的斜率划分得到过渡区及辨别区。通过对弹性图像的关注区域进行过渡区提取，得到图像数据的定量信息，可以从繁多复杂的信息中有目的或有针对性的获取特征。同时，也可以根据丰富的数据进行更加精确的数据分析。
文档编号G06T7/00GK102521823SQ20111035038
公开日2012年6月27日 申请日期2011年11月8日 优先权日2011年11月8日
发明者刘新, 张丽娟, 郑海荣, 钟耀祖, 黄文慧 申请人:中国科学院深圳先进技术研究院