一种人类胚胎心脏超声图像去噪方法技术

技术编号:21894451 阅读:94 留言:0更新日期:2019-08-17 15:30
本发明专利技术公开了一种人类胚胎心脏超声图像去噪方法,包括以下步骤:S1:获取具有时间序列和空间序列特征的超声图像数据集并选取中心图像,确定该中心图像的相邻图像;S2:将中心图像中的当前待处理像素标记为中心像素,计算该中心像素与相邻图像对应搜索域中每个像素的相似度,S3:根据相似度计算相邻图像对应的中心像素灰度值并对该中心像素灰度值进行取平均值操作得到中心像素的最终灰度值,对中心图像的每一个像素采用上述方法计算对应的最终灰度值,遍历整张中心图像后即获得去噪后的清晰图像。本方法在去噪过程中可以自由调整参数,在去噪效果与时间效率方面做出平衡,本方法不需要训练集的支撑,易于编程实现,算法复杂度较低。

A Noise Removal Method for Human Embryonic Heart Ultrasound Images

【技术实现步骤摘要】
一种人类胚胎心脏超声图像去噪方法
本专利技术涉及图像处理
,尤其涉及一种人类胚胎心脏超声图像去噪方法。
技术介绍
目前,超声技术是检测胎儿心脏是否健康的重要手段,但是由于胎儿的特殊性,超声需要额外穿过胎儿母亲的腹部脂肪才能针对胚胎心脏进行成像,这就导致超声图像往往会存在比其他情况下更多的伪影与噪声.针对这种情况,目前的解决方法有以下几种,首先是依据医疗人员的经验,利用未处理的超声图像进行诊断,这种方法对医疗人员的技术水平要求非常高,会消耗大量时间和精力,并且容易造成误诊;第二个是利用传统滤波方式去除噪声,这种方法虽然有一定的去噪效果,但是由于算法未针对超声图像进行优化,导致了图像中一些关键信息的丢失,不利于医生诊断;最后是利用机器学习的方法,使用大量的人类胚胎心脏的超声图像作为训练集,最后得到去噪效果,这种方法需要大量的数据作为训练集,而临床上这些数据很难获得,并且这种方法对硬件的要求很高,耗时也比较长。
技术实现思路
根据现有技术存在的问题,本专利技术公开了一种人类胚胎心脏超声图像去噪方法,包括以下步骤:S1:获取具有时间序列和空间序列特征的超声图像数据集并选取中心图像,确定该中心图像的相邻图像;S2:将中心图像中的当前待处理像素标记为中心像素,计算该中心像素与相邻图像对应搜索域中每个像素的相似度:设置中心像素对应的搜索域、计算搜索域对应累计方差,计算平均累计方差、计算搜索域中每个像素的邻域方差、计算搜索域中每个像素的邻域平均欧氏距离、计算搜索域中每个像素与中心像素的相似度;S3:根据相似度计算相邻图像对应的中心像素灰度值并对该中心像素灰度值进行取平均值操作得到中心像素的最终灰度值,对中心图像的每一个像素采用上述方法计算对应的最终灰度值,遍历整张中心图像后即获得去噪后的清晰图像。进一步的,所述中心图像的相邻图像选取方式为:将中心图像的时间序列当中相邻的两张图像和空间序列当中相邻的两张图像共4张作为相邻图像。进一步的,所述中心图像的每一个像素与相邻图像对应搜索域的相似度采用如下方式确认:S21:在中心图像中选取中心像素,将中心像素在相邻图像中相同位置的像素定义为目标像素,以目标像素为中心m×m像素范围的区域定义为搜索域,对于每一个搜索域当中的每一个像素,以该像素为中心的n×n像素范围定义为邻域;S22:计算相邻图像对应的搜索域中的像素P的累计方差:设像素P的灰度值为s,像素P对应的邻域中所有像素的灰度值分别表示为ti,i∈[1,n2],由此计算像素P在这一搜索域上的累计方差为:S23:对于中心像素对应的所有搜索域,均按照S22采用的方法计算累计方差,则该像素对应的平均累计方差为:S24:遍历整个搜索域,对搜索域中的每一个像素都按照S23采用的方法计算平均累计方差;S25:计算搜索域中每个像素的邻域方差:对于一个搜索域的一个像素P,由S23计算出其平均累计方差e,则像素P对应的邻域方差为:S26:计算邻域方差对应的高斯加权权重:设σ为高斯标准差,h为滤波系数,此时像素P对应的高斯加权权重表示为:S27:遍历整个搜索域,对于该搜索域当中的每一个像素按照S26所提出的方法计算高斯加权权重,将高斯加权权重表示为邻域平均欧氏距离,对于该搜索域当中的第j个像素所对应的高斯加权权重记为WGj,j∈[1,m2];S28:将一个搜索域中所有像素对应的高斯加权权重的和定义为归一化系数,由S27获取搜索域中第j个像素对应的高斯加权权重为WGj,则归一化系数表示为:S29:计算中心像素与搜索域中像素的相似度,其中中心像素和搜索域中的第j个像素的相似度表示为:S3中具体采用如下方式:S31:计算第k张相邻图像对应的中心像素的灰度值:设sj为此图像搜索域中的第j个像素的灰度值,设第j个像素对应的相似度为Wj,则此相邻图像对应的中心像素灰度值为:S32:计算中心像素的最终像素值:对于相邻图像采用S31提出的方法计算对应中心像素的灰度值,对相邻图像的对应中心像素的灰度值进行取平均操作即得到中心像素的最终像素值:S33:对于整张中心图像,遍历其所有的像素,利用上述方案得到对应的像素值最终得到整幅图像的去噪结果。由于采用了上述技术方案,本专利技术提供的一种人类胚胎心脏超声图像去噪方法,该方法基于时间序列和空间序列对人类胚胎心脏超声图像进行去噪处理,对质量较低的图像也能有较好的处理效果,对之后的医生诊断和三维重建等工作提供了比较好的前期准备工作。本方法在去噪过程中可以自由调整参数,在去噪效果与时间效率方面做出平衡,与基于机器学习的方法相比,本方法不需要训练集的支撑,易于编程实现,算法复杂度较低。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。图1为本专利技术中方法的流程图图2为本专利技术中输入图像的示意图图3为本专利技术所用的中心图像原图图4为本专利技术中心图像与搜索域示意图图5为本专利技术搜索域与邻域示意图图6为本专利技术的最终结果具体实施方式为使本专利技术的技术方案和优点更加清楚,下面结合本专利技术实施例中的附图,对本专利技术实施例中的技术方案进行清楚完整的描述:如图1所示的一种人类胚胎心脏超声图像去噪方法,具体包括以下步骤:S1:获取具有时间序列和空间序列特征的超声图像数据集并选取中心图像,确定该中心图像的相邻图像,如图2所示;S11:将病例数据包含的超声图像转化为灰度图像;S12:选择一幅图像为中心图像,中心图像原图如图3所示,获得其在时间序列当中相邻的两张图像和空间序列当中相邻的两张图像共4张作为相邻图像。S2:将中心图像中的当前待处理像素标记为中心像素,计算该中心像素与相邻图像对应搜索域中每个像素的相似度:设置中心像素对应的搜索域、计算搜索域对应累计方差,计算平均累计方差、计算搜索域中每个像素的邻域方差、计算搜索域中每个像素的邻域平均欧氏距离、计算搜索域中每个像素与中心像素的相似度。S21:确定搜索域与邻域:如图4所示,在中心图像当中选择一点作为中心像素,将其在4张相邻图像当中相同位置的像素称为目标像素,如图5所示,对于一个目标像素,以它为中心m×m像素范围的区域称为搜索域,对于每一个搜索域当中的每一个像素,以它为中心的n×n像素范围定义为邻域;对于一个搜索域中的一个像素,以它为中心的n×n像素范围定义为邻域,搜索域中的每一个像素都对应一个邻域。S22:计算累计方差:对于第k张相邻图像对应的搜索域中的一个像素P,P的灰度值表示为s,P对应的邻域中所有像素的灰度值分别表示为ti,i∈[1,n2],由此可计算P在这一搜索域上的累计方差为:S23:计算平均累计方差:对其他的搜索域,均按照S22采用的方法计算累计方差,则该像素对应的平均累计方差为:S24:遍历整个搜索域,对搜索域中的每一个像素都按照S23采用的方法计算平均累计方差;S25:计算每个像素的邻域方差:对于一个搜索域的一个像素P,由S23可以计算出其平均累计方差e,则像素P对应的邻域方差为:S26:计算邻域方差对应的高斯加权权重:设σ为高斯标准差,h为滤波系数,此时像素P对应的高斯加权权本文档来自技高网...

【技术保护点】
1.一种人类胚胎心脏超声图像去噪方法,其特征在于包括以下步骤:S1:获取具有时间序列和空间序列特征的超声图像数据集并选取中心图像,确定该中心图像的相邻图像;S2:将中心图像中的当前待处理像素标记为中心像素,计算该中心像素与相邻图像对应搜索域中每个像素的相似度:设置中心像素对应的搜索域、计算搜索域对应累计方差,计算平均累计方差、计算搜索域中每个像素的邻域方差、计算搜索域中每个像素的邻域平均欧氏距离、计算搜索域中每个像素与中心像素的相似度;S3:根据相似度计算相邻图像对应的中心像素灰度值并对该中心像素灰度值进行取平均值操作得到中心像素的最终灰度值,对中心图像的每一个像素采用上述方法计算对应的最终灰度值,遍历整张中心图像后即获得去噪后的清晰图像。

【技术特征摘要】
1.一种人类胚胎心脏超声图像去噪方法,其特征在于包括以下步骤:S1:获取具有时间序列和空间序列特征的超声图像数据集并选取中心图像,确定该中心图像的相邻图像;S2:将中心图像中的当前待处理像素标记为中心像素,计算该中心像素与相邻图像对应搜索域中每个像素的相似度:设置中心像素对应的搜索域、计算搜索域对应累计方差,计算平均累计方差、计算搜索域中每个像素的邻域方差、计算搜索域中每个像素的邻域平均欧氏距离、计算搜索域中每个像素与中心像素的相似度;S3:根据相似度计算相邻图像对应的中心像素灰度值并对该中心像素灰度值进行取平均值操作得到中心像素的最终灰度值,对中心图像的每一个像素采用上述方法计算对应的最终灰度值,遍历整张中心图像后即获得去噪后的清晰图像。2.根据权利要求1所述的一种人类胚胎心脏超声图像去噪方法,其特征还在于:所述中心图像的相邻图像选取方式为:将中心图像的时间序列当中相邻的两张图像和空间序列当中相邻的两张图像共4张作为相邻图像。3.根据权利要求2所述的一种人类胚胎心脏超声图像去噪方法,其特征还在于:所述中心图像的每一个像素与相邻图像对应搜索域的相似度采用如下方式确认:S21:在中心图像中选取中心像素,将中心像素在相邻图像中相同位置的像素定义为目标像素,以目标像素为中心m×m像素范围的区域定义为搜索域,对于每一个搜索域当中的每一个像素,以该像素为中心的n×n像素范围定义为邻域;S22:计算相邻图像对应的搜索域中的像素P的累计方差:设像素P的灰度值为s,像素P对应的邻域中所有像素的灰度值分别表示为ti,i∈[1,n2],由此计算像素P在这一搜索域上的累计方差为:S23:...

【专利技术属性】
技术研发人员:刘斌许钊
申请(专利权)人:大连理工大学
类型:发明
国别省市:辽宁,21

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1