一种基于图像金字塔的鲁棒自适应帧相关方法及系统技术方案

本发明专利技术属于超声成像技术领域，具体涉及一种基于图像金字塔的鲁棒自适应帧相关方法及系统，其中，涉及的一种基于图像金字塔的鲁棒自适应帧相关方法，包括以下步骤：S1、根据噪声图像构建金字塔；S2、根据金字塔的顶层图像计算帧相关系数模板；S3、根据系数模板进行尺度变换；S4、根据尺度变换后的系数模板实现对噪声图像的帧相关处理；本发明专利技术通过金字塔顶层图像进行帧相关系数的计算，能够大大减小算法的复杂度；金字塔每层图像在构建过程中均会进行卷积核的平滑，以使图像金字塔具有良好的去噪效果。

【技术实现步骤摘要】
一种基于图像金字塔的鲁棒自适应帧相关方法及系统
本专利技术属于超声成像
，具体涉及一种基于图像金字塔的鲁棒自适应帧相关方法及系统。
技术介绍
医学超声是一种低成本、非侵入、实时性的临床检查手段，然而实际的医学超声图像中往往存在较大的噪声，噪声的存在会降低图像的信噪比和对比分辨率等，最终降低了图像质量。超声图像中的噪声主要来源于两类，一类是内部电路、热效应等引起的随机噪声，另一类是超声回波相干散射形成的散斑噪声。现有的去噪方法主要分为三类，第一类是基于复合的方法，如频率复合，空间复合和时间复合。第二类是基于噪声模型的方法，该方法的缺点在于噪声特征过于复杂，噪声模型的构建过于理想化。最后一类是基于后处理的方法，具体可以分为单尺度和多尺度两种方法。帧相关即时间复合方法可以在不显著增加系统复杂度的基础上实现较好的去噪效果。帧相关的基本模型为：yt＝ayt-1+(1-a)xt，式中的x、y分别表示帧相关前、后的图像，下标表示时刻，α为帧相关系数。现有的帧相关方法主要分为非自适应和自适应两类，非自适应即模型中的α为固定值，自适应即α是根据不同时刻的图像实时计算出来的，具体可以基于两帧或多帧进行计算。现有的自适应帧相关方法是在原图像尺寸上进行计算的，而原图在去噪之前含有较多噪声，从而使得相关性的计算不够准确，最终导致帧相关的效果不佳；现有的解决方法是在帧相关之前对图像进行空间平滑、数字滤波等预处理，但该方法的缺点是难以在去除噪声和保留图像特征之间达到很好的平衡；为此，有必要对此进行改进，以克服实际应用中的不足。
技术实现思路
基于现有技术中存在的上述缺点和不足，本专利技术的目的之一是至少解决现有技术中存在的上述问题之一或多个，换言之，本专利技术的目的之一是提供满足前述需求之一或多个的一种基于图像金字塔的鲁棒自适应帧相关方法及系统。为了达到上述专利技术目的，本专利技术采用以下技术方案：一种基于图像金字塔的鲁棒自适应帧相关方法，包括以下步骤：S1、根据噪声图像构建金字塔；S2、根据金字塔的顶层图像计算帧相关系数模板；S3、根据系数模板进行尺度变换；S4、根据尺度变换后的系数模板实现对噪声图像的帧相关处理。作为优选方案，所述步骤S1具体为：S11、设置待分解的金字塔层数L，并将原始噪声图像I0定义为第0层；S12、计算卷积核kernel，并对kernel进行归一化处理；x＝-3σ,-3σ+1,...,3σ-1,3σ其中，σ为高斯核的标准差；S13、初始化i＝0，将金字塔第i层图像Ii与kernel进行卷积，即Ii＝conv(Ii,kernel)；S14、提取Ii第x行y列作为金字塔第i+1层的图像Ii+1；依次计算金字塔下一层的图像，直至完成金字塔L层的计算，以完成整个图像金字塔的构建；S15、根据前一帧和当前帧分别进行对应图像金字塔的构建，并将前一帧和当前帧的金字塔顶层图像分别记为IPt-1和IPt。作为优选方案，所述步骤S2具体为：S21、根据IPt-1和IPt计算当前像素的相关度C；S22、根据相关度计算帧相关系数a，并依次计算下一像素对应的帧相关系数，以得到所有像素的帧相关系数；S23、将所有像素的帧相关系数构成系数模板A，其中，A映射于0-1之间。作为优选方案，所述步骤S21中通过平均绝对差、均方差、均方根误差、相关系数以及相似度函数计算相关度。作为优选方案，所述步骤S22中通过递减函数实现相关度C与帧相关系数a之间的转换，所述递减函数包括线性函数、分段线性函数、指数函数和对数函数。作为优选方案，所述步骤S3具体为：S31、对系数模板A进行去噪处理；S32、将去噪后的系数模板变换为金字塔底层图像的尺寸，变换后的系数模板为A’。作为优选方案，所述步骤S4具体为：对变换后的系数模板A’进行帧相关处理，帧相关公式为：其中，上标加f的表示帧相关之后的图像，下标t表示时刻。本专利技术还提供一种基于图像金字塔的鲁棒自适应帧相关系统，包括：构建模块，用于根据噪声图像构建金字塔；第一计算模块，用于根据金字塔的顶层图像计算帧相关系数模板；变换模块，用于根据系数模板进行尺度变换；处理模块，用于根据尺度变换后的系数模板实现对噪声图像的帧相关处理。作为优选方案，所述第一计算模块包括：第二计算模块，用于根据IPt-1和IPt计算当前像素的相关度C；第三计算模块，用于根据相关度计算帧相关系数a，并依次计算下一像素对应的帧相关系数，以得到所有像素的帧相关系数。作为优选方案，所述第二计算模块通过平均绝对差、均方差、均方根误差、相关系数以及相似度函数计算相关度。本专利技术与现有技术相比，有益效果是：本专利技术通过金字塔顶层图像进行帧相关系数的计算，能够大大减小算法的复杂度。本专利技术的金字塔每层图像在构建过程中均会进行卷积核的平滑，以使图像金字塔具有良好的去噪效果。附图说明图1是本专利技术实施例一的一种基于图像金字塔的鲁棒自适应帧相关方法的流程图；图2是本专利技术实施例二的一种基于图像金字塔的鲁棒自适应帧相关系统的结构图。具体实施方式为了更清楚地说明本专利技术实施例，下面将对照附图说明本专利技术的具体实施方式。显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，并获得其他的实施方式。实施例一：本实施例提供一种基于图像金字塔的鲁棒自适应帧相关方法，如图1所示，包括以下步骤：S1、根据噪声图像构建金字塔；S2、根据金字塔的顶层图像计算帧相关系数模板；S3、根据系数模板进行尺度变换；S4、对尺度变换后的噪声图像进行帧相关处理。具体地，步骤S1根据噪声图像构建金字塔，具体为：S11、设置待分解的金字塔层数L＝3，原始噪声图像I0被定义为第0层。S12、计算高斯卷积核kernel，计算公式如下：其中，表示向上取整，取σ＝1.5；对kernel进行归一化，即使得kernel所有元素之和等于1；S13、初始化i＝0，将第i层的图像Ii与kernel进行卷积，用公式表示为：Ii＝conv(Ii,kernel)。S14、提取Ii第x行y列作为第i+1层的图像Ii+1；如果Ii的行数M为偶数，则x＝2p,p＝0,1,2,...,M/2-1，否则同理如果Ii的列数N为偶数，则y＝2q,q＝0,1,2,...,M/2-1，否则符号表示向下取整；依次计算金字塔下一层的图像，直至完成金字塔L层的计算，以完成图像金字塔的构建；S15、根据前一帧和当前帧分别进行对应图像金字塔的构建，对当前帧和前一帧图像分别进行上述步骤S13至S14完成相对本文档来自技高网...

【技术保护点】
1.一种基于图像金字塔的鲁棒自适应帧相关方法，其特征在于，包括以下步骤：/nS1、根据噪声图像构建金字塔；/nS2、根据金字塔的顶层图像计算帧相关系数模板；/nS3、根据系数模板进行尺度变换；/nS4、根据尺度变换后的系数模板实现对噪声图像的帧相关处理。/n

【技术特征摘要】
1.一种基于图像金字塔的鲁棒自适应帧相关方法，其特征在于，包括以下步骤：
S1、根据噪声图像构建金字塔；
S2、根据金字塔的顶层图像计算帧相关系数模板；
S3、根据系数模板进行尺度变换；
S4、根据尺度变换后的系数模板实现对噪声图像的帧相关处理。


2.根据权利要求1所述的一种基于图像金字塔的鲁棒自适应帧相关方法，其特征在于，所述步骤S1具体为：
S11、设置待分解的金字塔层数L，并将原始噪声图像I0定义为第0层；
S12、计算卷积核kernel，并对kernel进行归一化处理；



其中，σ为高斯核的标准差；
S13、初始化i＝0，将金字塔第i层图像Ii与kernel进行卷积，即Ii＝conv(Ii,kernel)；
S14、提取Ii第x行y列作为金字塔第i+1层的图像Ii+1；依次计算金字塔下一层的图像，直至完成金字塔L层的计算，以完成整个图像金字塔的构建；
S15、根据前一帧和当前帧分别进行对应图像金字塔的构建，并将前一帧和当前帧的金字塔顶层图像分别记为IPt-1和IPt。


3.根据权利要求2所述的一种基于图像金字塔的鲁棒自适应帧相关方法，其特征在于，所述步骤S2具体为：
S21、根据IPt-1和IPt计算当前像素的相关度C；
S22、根据相关度计算帧相关系数a，并依次计算下一像素对应的帧相关系数，以得到所有像素的帧相关系数；
S23、将所有像素的帧相关系数构成系数模板A，其中，A映射于0-1之间。


4.根据权利要求3所述的一种基于图像金字塔的鲁棒自适应帧相关方法，其特征在于，所述步骤S21中通过平均绝对差、均方差、均方根误差、相关系数以及相似度函数计算相关度。
<...

【专利技术属性】
技术研发人员：肖梦楠，
申请(专利权)人：聚融医疗科技杭州有限公司，
类型：发明
国别省市：浙江;33