一种从超声波图像中提取目标特征的图像处理方法与流程

本发明涉及一种从超声波图像中提取目标特征的图像处理方法，属于图像处理领域。

技术背景：

因为超声波图像在医学检测领域中扮演着十分重要的角色，所以超声波图像处理的精确度也对医学诊断更加至关重要。比如在孕早期，监测胚胎胎心的变化和周期是判断胎儿健康与否的一项重要指标。然而，由于在用超声波检测给孕妇胎儿检查胎心的过程中用血流影像上彩会能量过大，从而潜在地对胎儿健康有影响，因此本发明研究并提出了一种能从不上彩超声波图片中提取胎心所在特征区域并以此追踪胚胎胎心的跳动周期。由此，可以利用标准的算法统一地从超声波图片中获得胎心率的变化数据，替代人为选取胎心周期来计算胎心率的方法，消除人为误差。

具体来讲，由于心脏的周期性扩张和缩小，胎心区域在一段连续的超声波视频中的尺寸会周期性地变大变小。在孕早期，胚芽和卵黄囊连在一起在孕囊内部显现。而且由于孕早期胚芽很小，心脏在胚芽中所占面积比率很大。同时，在心跳周期中卵黄囊的尺寸并不会发生变化。因此，可以近似地使用整个胚芽和卵黄囊区域的尺寸发生的变化来代表胎心区域尺寸的变化。该发明着重讲述精确提取两个重要指标：整个胚芽和卵黄囊区域（胚胎区域）和孕囊区域的图像处理算法。

名词解释：

连通区域：复平面上的一个区域g，如果在其中任做一条简单闭曲线，而闭曲线的内部总属于g，就称g为单连通区域。一个区域如果不是单连通区域，就称为多连通区域。

胚胎区域：在该发明中特指卵黄囊和胚芽所占的整个区域。

上彩：指高清晰度的黑白b超再加上彩色多普勒，彩色多普勒超声一般是用自相关技术进行多普勒信号处理，把自相关技术获得的血流信号经彩色编码后实时地叠加在二维图像上，即形成彩色多普勒超声血流图像，由于采用的彩色多普勒超声，其辐射量大于黑白b超。

技术实现要素：

为解决上述问题，本发明公开了一种从超声波图像中提取目标特征的图像处理方法，主要包含几个部分：分割孕囊所在区域、从孕囊区域中分割出胚胎区域和计算胎心率。这几个主要部分里面分别包含了一系列的图像处理算法步骤，来最终实现精确地提取胚胎区域的目标。

为实现上述目的，本发明技术方案如下：

一种从超声波图片中提取目标特征的图像处理方法，其特征在于，包括如下步骤：

步骤一、读取超声波原图：将未上彩的超声波检查过程中录制的视频转化为超声波图像序列；

步骤二、精确分割孕囊区域：

a.从超声波图像序列的第一幅图开始，使用阈值方法将图片二值化：阈值是统一选取为th，低于阈值th的部分设定灰度值为0，高于阈值th的部分设定灰度值为1；二值化后的图像中，孕囊边缘有个缺口，即得到有缺口的孕囊区域图；

b.补齐缺口：首先使用凸包检测获取二值化以后孕囊区域的外轮廓线，然后填充轮廓线包围内部区域的空洞，得到凸包填充后的结果图；

c.对步骤a获得的有缺口的孕囊区域图和步骤b中获得的凸包填充后的结果图求差值，获得差值图；

d.对差值图使用距离变换的方法，计算差值图中各个连通区域的距离变化的结果图并给每个区域标上记号，找到最大内切圆所在的连通区域并提取出来，即为胚芽和卵黄囊所在的区域即目标胚胎区域；

步骤三、获得胎心率

统计得到的胚芽和卵黄囊区域内的像素个数后获得的总数作为胚胎区域的面积；

步骤四、计算获得胎心率的数值

对在整个超声波图像序列中的图片均进行步骤一到步骤三的同样处理，可得到完整超声波图像序列的图片中显示出来的胚胎区域的面积变化曲线，胚胎区域的面积变化即胚芽区域的面积变化；胚芽区域的面积变化即作为胎心的面积变化；根据变化曲线上相邻的两个峰值之间的距离，即计算获得胎心率的数值。

进一步的改进，

胎心率的计算方法如下所示：

fhr表示胎心率；fps表示图像序列的帧率；dm表示胚胎区域的面积变化曲线上相邻波峰之间的距离。

进一步的改进，所述步骤一的具体步骤如下：将从超声波检查过程中录制的.avi格式的视频转换为图片并存放在一个单独的文件夹中；通过输入命令框输入原图所在的文件夹路径，然后自动定位并读取所有设定图片格式的图片文件并统计图片文件总数。

进一步的改进，所述步骤二中的阈值方法中使用的阈值：th=130。

进一步的改进，还包括步骤五:步骤d中，提取胚胎区域的轮廓线，胚胎区域的轮廓线与步骤a中获得的有缺口的孕囊区域图进行叠加算法，获得补齐孕囊区域边界以后的完整孕囊结果图。

本发明优点：

不用上彩，且能够精确的从超声波图片中提取出重要的目标区域，并且不会引入多余的边界像素，从而提高了超声波图像的检测指标的准确度。

附图说明：

附图1：该发明中算法的流程图；

附图2：一张超声波原图；

附图3：原图二值化以后获得的孕囊区域图；

附图4：使用凸包检测来获得孕囊区域外部轮廓线并填充空洞后的结果图；

附图5：将附图4与附图3做差值来获得差值图；

附图6：在附图5上使用距离变换以后的结果图；

附图7：精确提取出来的胚胎区域结果图；

附图:8：提取胚胎区域轮廓线并和附图3叠加后获得精确孕囊轮廓图；

附图9：精确提取出来的孕囊区域结果图。

具体实施方式：

为了更具体地阐述该发明，以一张孕囊所在区域中存在白色噪音并且胚胎区域与孕囊边界相连接的超声波图像为例，从超声波图像处理过程中比较复杂的情况来详细地讲述本发明的实施方式。本发明的流程如附图1所示，通过以下四个部分的数据处理实现：

一、读取超声波原图

将从超声波检查过程中录制的.avi格式的视频流转换为图片并存放在一个单独的文件夹中。该算法通过输入命令框输入原图所在的文件夹路径，然后自动定位并读取所有特定图片格式的文件并统计文件总数。

二、分割孕囊区域

这个步骤中有一系列的图像处理算法小步骤来实现在超声波图像中精确的提取出孕囊所在的区域。子步骤分别如下：

a.从一个完整的超声波图像序列的第一幅图开始，用阈值方法将图片二值化，阈值是统一选取为130。低于阈值的部分设定灰度值为0，高于阈值部分设定灰度值为1。

b．使用凸包检测来获取二值化以后孕囊区域的外轮廓线，然后使用填充空洞的方法，填充满轮廓线内部所包含的范围。

c.对子步骤二a和子步骤二b中获得的结果求差值获得差值图，找到二者之间的有差异的区域。

d.对差值图使用距离变换的方法，获得差值图中连通区域的距离变换结果图。通过标记连通区域和寻找最大内切圆所在的连通区域来找到胚胎区域，也即二值化以后孕囊结果图中的缺口区域。

三、精确提取胚胎区域

统计提取出来的胚胎区域的像素个数总和，作为胚胎区域面积。

四、计算胎心率

对整个超声波图片序列的中的所有图片重复以上步骤一，二和三中的所有步骤，获得完整的胚胎区域面积的变化曲线。根据相邻的两个峰值之间的距离和胎心率fhr的计算公式，即可获得该序列图片中呈现出来的胎心率结果。

五、精确提取孕囊区域

提取步骤二d中找到的胚胎区域的轮廓线，并与步骤二a中获得的有精确边界线的孕囊区域图进行叠加算法，即可获得有完整且精确边界的孕囊区域轮廓图。然后就可以获得精确提取出来的孕囊区域结果图。孕囊区域结果图用于判断孕囊的形状，从而诊断孕囊时候存在病变等情况。

该算法可以使用统一的计算机图像处理步骤来处理不同图像质量、噪音多、边界连接等等复杂情况的超声波图像，并且精确地提取出目标区域。能够为监测孕期胎儿的重要健康指标作出贡献。

具体示例如下：

1.首先，读取一张超声波图片序列中的一张原图，如附图2所示。

2.根据原图中像素的灰度值使用阈值方法，将原图二值化为如附图3所示的黑白图。这一步骤中的阈值是经过了一系列测试多组图片获得一个最合适阈值，130。

3.从附图3中可见，由于孕囊中的胚芽和卵黄囊区域与边界相连接，二值化以后的孕囊区域有缺口。这一步骤中使用凸包检测功能获得孕囊区域的外轮廓线，并使用填充空洞算法来补全轮廓线内部的所有区域，结果如附图4所示。

4.将附图4与二值化的附图3使用求差算法，来获得差值图，如附图5.

5.从差值图-附图5中可以看到很多大小不等的连通区域。然后使用距离变换，获得差值图的距离变换结果图，如附图6所示。

6.在附图6中将各个连通区域进行标记，找到最大内切圆所在的连通区域（距离变换的结果图中内部像素和边界像素距离最大的连通区域）并将其提取出来，即为目标胚胎区域，结果图如附图7所示。

7.提取胚胎区域的轮廓线，然后与二值化的附图3进行叠加算法，补齐整个孕囊区域的轮廓线，结果如附图8所示。这样，即可提取获得边界完整且精确的孕囊结果图，如附图9所示。

从结果图中可见，目标区域-胚胎区域和孕囊区域，都被相当精确地提取了出来，边界完整且几乎没有带入多余的边界信息。

然后对整个超声波图像序列重复步骤1-6的操作，即可获得整个序列的胎心区域面积的变化曲线，根据胎心率计算公式，即可获得该图像序列中呈现出来的胎心率结果。

以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变换，均仍属于本发明技术方案的保护范围内。