基于分级搜索的大视场自适应壁滤波方法及系统与流程

1.本发明涉及医学图像处理领域，特别涉及一种基于分级搜索的大视场自适应壁滤波方法及系统。


背景技术：

2.超快超声平面波成像技术能提供高帧频的大视场同步回波数据，为大视场血流图像提供了更多有效原始信息，普通滤波算法在处理大视场数据方面有诸多不足。
3.现有壁滤波方法大致分三种：(1)传统时间域滤波器每次滤波只能处理单点时间域信息，通过逐点滤波拼凑出二维图像，在大视场图像处理上存在效率低抗干扰力差等问题，速度慢，滤波效果差。(2)新型的奇异值分解滤波器(svd滤波器，singular value decompostion)能够一次输入全部原始数据，提高了处理效率，但对整幅图像采取“一刀切”的滤波截止设置，面对大视场图像中非均匀分布的有效信息存在局部过度或欠缺处理。(3)改进的分区域svd滤波器，将整幅图画分为许多小区域各自进行特异性处理，提高了滤波质量，但牺牲了svd滤波器高效的优势。
4.所以，现在需要一种更可靠的方案。


技术实现要素：

5.本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于分级搜索的大视场自适应壁滤波方法及系统。本发明能够综合全域svd高效和分区域svd滤波质量高的优点，对整个观察区域可快速获取高质量的血流图像。
6.为实现上述目的，本发明采用的技术方案是：一种基于分级搜索的大视场自适应壁滤波方法，该方法用于对超声图像进行壁滤波处理，该方法先采用大窗(一级窗)将整幅图像分为若干区域，进行滤波并评估各区域的血流信息含量，将血流信息含量高于阈值的区域提取出来用细化的小窗(二级窗)再进行局部滤波处理，每个小窗根据对应区域的噪声水平采用最优化的特异性滤波参数进行处理，最后进行图像重构，输出完整的壁滤波处理后的超声血流图像。
7.优选的是，该方法包括以下步骤：
8.1)获取多帧时间连续的二维平面波超声回波信号，经iq解调得到大视场连续帧二维平面波解调数据；
9.2)一级窗评估：对整幅图像以大窗划分为若干区域，同步评估各区域的血流信息含量，将血流信息含量高于阈值的区域提取出来；
10.3)二级窗细化：计算并设置合适的小窗，将步骤2)提取出来的血流信息含量高于阈值的区域划分成一个个相邻区域之间相互有重叠的小窗；
11.4)每个小窗独立进行自适应壁滤波处理；
12.5)图像重构：将各个小窗独立自适应滤波后输出的结果，叠加到全图相应位置，通过统计及加权平均获得归一化图像，输出重构的完整血流图像。
13.优选的是，单个小窗进行自适应壁滤波处理的步骤为：进行小窗数据时空域重构，得到分解矩阵，通过奇异值分解提取特征值，绘制特征曲线，利用特征曲线计算最优化的滤波参数，采用最优化的滤波参数对该小窗数据进行滤波，重构图像并输出结果。
14.本发明还提供一种基于分级搜索的大视场自适应壁滤波系统，其包括：
15.iq解调模块，其对获取的多帧时间连续的二维平面波超声回波信号进行iq解调，得到大视场连续帧二维平面波解调数据；
16.一级窗评估模块，其对整幅图像以大窗划分为若干区域，同步评估各区域的血流信息含量，将血流信息含量高于阈值的区域提取出来；
17.二级窗细化模块，其计算并设置合适的小窗，将一级窗评估模块提取出来的血流信息含量高于阈值的区域划分成一个个相邻区域之间相互有重叠的小窗；
18.自适应壁滤波处理模块，其对每个小窗独立进行自适应壁滤波处理；
19.以及图像重构模块，其将自适应壁滤波处理模块得到的各个小窗独立自适应滤波后输出的结果叠加到全图相应位置，通过统计及加权平均获得归一化图像，输出重构的完整血流图像。
20.本发明还提供一种存储介质，其上存储有计算机程序，该程序被执行时用于实现如上所述的方法。
21.本发明还提供一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上所述的方法。
22.本发明的有益效果是：本发明提供的基于分级搜索的大视场自适应壁滤波方法能够适用于大容量数据的高效壁滤波方法，一方面，本发明通过分级搜索策略，只在有需要的部位分配算力，提高了运算效率；另一方面，本发明中滤波小窗虽然重叠，但彼此滤波过程间并无迭代，因此多个独立小窗可通过多线程同步计算，大大节约运算时间；
23.本发明的壁滤波方法能够很好的抑制杂波，一方面，通过细分滤波区域来降低滤波难度，单次滤波区域缩小使空间域信息相似度提高，有效信息更集中，更容易将杂波信号从中剔除；另一方面，针对每个独立区域的数据特性设计最优化滤波参数，大大提高了滤波质量。
附图说明
24.图1为本发明的基于分级搜索的大视场自适应壁滤波方法的流程图。
具体实施方式
25.下面结合实施例对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
26.应当理解，本文所使用的诸如“具有”、“包含”以及“包括”术语并不排除一个或多个其它元件或其组合的存在或添加。
27.实施例1
28.本实施例的一种基于分级搜索的大视场自适应壁滤波方法，该方法用于对超声图像进行壁滤波处理，该方法先采用大窗(一级窗)将整幅图像分为若干区域，进行滤波并评
估各区域的血流信息含量，将血流信息含量高于阈值的区域(有效信息丰富、图像复杂)提取出来用细化的小窗(二级窗)再进行局部滤波处理，每个小窗根据对应区域的噪声水平采用最优化的特异性滤波参数进行处理，最后进行图像重构，输出完整的壁滤波处理后的超声血流图像。
29.参照图1，本实施例中的基于分级搜索的大视场自适应壁滤波方法包括以下步骤：
30.1)通过阵列超声探头(包括凸阵和线阵)采集多帧时间连续的二维平面波超声回波信号，经iq解调得到大视场连续帧二维平面波解调数据；
31.2)一级窗评估：对整幅图像以大窗划分为若干区域，同步评估各区域的血流信息含量，将血流信息含量高于阈值的区域提取出来；
32.3)二级窗细化：计算并设置合适的小窗，将步骤2)提取出来的血流信息含量高于阈值的区域划分成一个个相邻区域之间相互有重叠的小窗；
33.4)每个小窗独立进行自适应壁滤波处理：
34.进行小窗数据时空域重构，得到分解矩阵，通过奇异值分解提取特征值，绘制特征曲线，利用特征曲线计算最优化的滤波参数，采用最优化的滤波参数对该小窗数据进行滤波，重构图像并输出结果。每个小窗的数据独立滤波，由于各窗互不干扰，该过程可同步进行提高处理速度。
35.5)图像重构：将各个小窗独立自适应滤波后输出的结果，叠加到全图相应位置，通过统计及加权平均获得归一化图像，输出重构的完整血流图像。
36.本发明适合处理大视场成像带来的大容量数据，能够快速提供高质量血流图像。平面波技术的发展及硬件设备的升级，使高帧频大视场二维同步超声信号采集成为现实，二维同步的声信号携带了大量空间方面的信息，同时高重复频率的连续帧信号又携带了大量时间方面的信息，处理如此大量的数据，需要兼顾速度和滤波质量的壁滤波方法。本发明提供的算法，利用分区域独立自适应滤波保证滤波质量，又因为各区域无迭代过程可通过同步运算节约运算时间，能够快速处理大视场数据提供高质量血流图像。
37.大视场图像中成像目标分布不均匀性导致其中必然存在有效信息含量高的核心区域(在血流图中表现为血管丰富区域)，和几乎不含有效信息的边缘区域(在血流图中表现为无血流区域)，无区别的小窗搜索算法会在边缘区浪费大量算力。本发明采用分级搜索策略，以粗略搜索确定核心区，再以精确搜索更精准地提取有效信息，可以在保证滤波质量的前提下，大大提高运算效率。
38.本发明中，各小窗数据进行独立处理，根据特征曲线分析能量分布情况，计算最优化滤波参数，保证对每个小窗进行最优化参数滤波，从而提高了滤波质量。
39.实施例2
40.本实施例提供一种基于分级搜索的大视场自适应壁滤波系统，其利用实施例1的方法进行超声图像的壁滤波处理，该系统包括：
41.iq解调模块，其对获取的多帧时间连续的二维平面波超声回波信号进行iq解调，得到大视场连续帧二维平面波解调数据；
42.一级窗评估模块，其对整幅图像以大窗划分为若干区域，同步评估各区域的血流信息含量，将血流信息含量高于阈值的区域提取出来；
43.二级窗细化模块，其计算并设置合适的小窗，将一级窗评估模块提取出来的血流
信息含量高于阈值的区域划分成一个个相邻区域之间相互有重叠的小窗；
44.自适应壁滤波处理模块，其对每个小窗独立进行自适应壁滤波处理；
45.以及图像重构模块，其将自适应壁滤波处理模块得到的各个小窗独立自适应滤波后输出的结果叠加到全图相应位置，通过统计及加权平均获得归一化图像，输出重构的完整血流图像。
46.本实施例还提供一种存储介质，其上存储有计算机程序，该程序被执行时用于实现实施例1的方法。
47.本实施例还提供一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现实施例1的方法。
48.尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节。