CN111815527A - 基于Weibull分布的混合高阶变分超声图像去噪方法 - Google Patents
基于Weibull分布的混合高阶变分超声图像去噪方法 Download PDFInfo
- Publication number
- CN111815527A CN111815527A CN202010550145.4A CN202010550145A CN111815527A CN 111815527 A CN111815527 A CN 111815527A CN 202010550145 A CN202010550145 A CN 202010550145A CN 111815527 A CN111815527 A CN 111815527A
- Authority
- CN
- China
- Prior art keywords
- image
- order
- ultrasonic image
- function
- noise
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 6
- 230000004927 fusion Effects 0.000 claims abstract description 4
- 238000002604 ultrasonography Methods 0.000 claims description 16
- 238000003384 imaging method Methods 0.000 claims description 8
- 238000000926 separation method Methods 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000000007 visual effect Effects 0.000 abstract description 4
- 230000002708 enhancing effect Effects 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 210000003205 muscle Anatomy 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 238000002592 echocardiography Methods 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000012631 diagnostic technique Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 210000002249 digestive system Anatomy 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 210000003754 fetus Anatomy 0.000 description 1
- 230000000762 glandular Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000011551 log transformation method Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 230000002485 urinary effect Effects 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
Abstract
基于Weibull分布的混合高阶变分超声图像去噪方法,对超声图像进行噪声分布拟合,采用Weibull概率密度函数,对选定的图像局部均匀区域进行灰度直方图拟合,并对拟合分布参数采用最大似然法进行估计;根据拟合的噪声分布,推导对数真实无噪图像的似然函数,将似然函数进行负的对数变换,从而得到最小化能量函数的数据保真项;将对数真实无噪图像一阶正则化表达式和二阶正则化表达式进行加权融合,构建最小化能量函数的混合高阶正则化项;将数据保真项和混合高阶正则化项形成最小化能量函数模型,采用Split‑Bregman迭代方法进行快速求解,迭代收敛后,利用指数变换得到去噪超声图像。本发明方法能有效降低散斑噪声的同时更好地保留图像的细节和边缘信息,从而增强超声图像的视觉判读性。
Description
技术领域
本发明涉及图像处理技术领域,具体涉及一种基于Weibull分布的混合高阶变分超声图像去噪方法。
背景技术
目前,超声影像已成为临床多种疾病诊断的首选方法,并已发展成为一种能提供多种参数的系列诊断技术。超声成像方法常用来判断脏器的位置、大小、形态,确定病灶的范围和物理性质,提供一些腺体组织的解剖图,鉴别胎儿的正常与异常,在妇产科及心血管系统、消化系统、泌尿系统的应用十分广泛。超声成像采用回声原理,由仪器探头向人体发射一束超声进入体内,并进行线形、扇形或其它形式的扫描,遇到不同声阻抗的两种组织的交界面,即有超声反射回来,由探头接收后,经过信号放大或信息处理,形成人体的断层图像。除了反射回波外,人体组织的微小结构作为散射体使声波产生散射后形成的散射回波也会被接收,但这部分回波相互干涉通常产生所谓的散斑噪声。由于该类噪声的存在,使得图像分辨率和对比度严重下降,影响计算机辅助的临床诊断,并为后续的组织或病灶体分割、识别、配准或三维重建等造成极大困难。因此超声图像散斑噪声去除一直是医学成像领域研究的热点问题。
多年来,国内外学者致力于超声图像的去噪研究,提出了多种不同类型的去噪方法,包括局部空域滤波、各项异性扩散滤波、非局部均值滤波、总变分方法、多尺度方法以及同态滤波方法等。这些方法,有的去噪效果不够理想,有的去噪算法复杂度较高。因此,寻求一种算法简单有效且去噪效果较理想的方案,是目前超声图像实时处理的必要条件,也是去噪算法移植到实际医疗设备的先决条件。
发明内容
本发明提出了一种基于Weibull分布的混合高阶变分超声图像去噪方法,通过拟合超声图像均匀区域得到整体噪声分布参数,以此构建最小化能量函数的数据保真项;并设计自适应参数以及混合高阶变分得到正则化项。两者共同作用,能有效降低散斑噪声的同时更好地保留图像的细节和边缘信息,从而增强超声图像的视觉判读性。
本发明采取的技术方案为:
基于Weibull分布的混合高阶变分超声图像去噪方法,包括以下步骤:
步骤1:对超声图像进行噪声分布拟合,采用Weibull概率密度函数,对选定的图像局部均匀区域进行灰度直方图拟合,并对拟合分布参数采用最大似然法进行估计;
步骤2:根据拟合的噪声分布,推导对数真实无噪图像的似然函数,将似然函数进行负的对数变换,从而得到最小化能量函数的数据保真项;
步骤3:将对数真实无噪图像一阶正则化表达式和二阶正则化表达式进行加权融合,构建最小化能量函数的混合高阶正则化项,且权值能根据图像边缘特征自适应调整;
步骤4:将步骤2的数据保真项和步骤3的混合高阶正则化项形成最小化能量函数模型,采用Split-Bregman迭代方法进行快速求解,迭代收敛后,得到对数去噪超声图像,最后通过指数变换得到去噪超声图像。
本发明一种基于Weibull分布的混合高阶变分超声图像去噪方法,技术效果如下:
1)、本发明利用图像的噪声分布参数构建数据保真项,该项可根据超声图像的噪声强度自适应改变;采用混合高阶变分项并设计自适应参数,在有效去除噪声时能更精确地描述边缘信息以及微小细节信息,同时避免一阶变分计算时产生的阶梯块效应,且引入更少的参数以减少算法复杂度。综合两项的优势使去噪算法能够在均匀区域有效降低散斑噪声,在细节区域有效保留或增强对比度。
2)、本发明通过图像自身灰度统计、自适应参数设计、高阶变分混合的形式对图像进行滤波,引入更少的参数降低算法难度。且有效降低噪声以及一阶变分框架阶梯块效应,更好地保留图像细节和边缘信息,从而增强超声图像的视觉判读性。
3)、本发明方法采用Split-Bregman迭代求解算法,收敛速度快,适宜进行硬件固化。
附图说明
图1为本发明设计实现流程图。
图2为本发明实施例中的待去噪超声图。
图3为本发明实施例中的去噪结果图。
图4为本发明方法对比三幅不同噪声强度测试图像的去噪结果的峰值信噪比(Peak Signal to Noise Ratio,PSNR)随数据保真项权系数α的变化曲线图。
图5为本发明方法对比三幅不同噪声强度测试图像去噪结果的峰值信噪比PSNR随调节系数θ1的变化曲线图。
图6为本发明方法对比三幅不同噪声强度测试图像去噪结果的峰值信噪比PSNR随调节系数θ2的变化曲线图。
图7为去噪结果对比的局部区域示意图。
图8(1)为原始超声图像的局部区域;
图8(2)为Exp超声图像去噪方法的处理结果图;
图8(3)为MIDAL超声图像去噪方法的处理结果图;
图8(4)为NTV超声图像去噪方法的处理结果图;
图8(5)为Huang-Yang超声图像去噪方法的处理结果图;
图8(6)为TM超声图像去噪方法的处理结果图;
图8(7)为本发明WHHVM超声图像去噪方法的处理结果图。
具体实施方式
基于Weibull分布的混合高阶变分超声图像去噪方法,包括以下步骤:
步骤1:对超声图像进行噪声分布拟合,采用Weibull概率密度函数,对选定的图像局部均匀区域进行灰度直方图拟合,并对拟合分布参数采用最大似然法进行估计。
所述步骤1中,真实无噪图像v,受到乘性噪声n的干扰,生成观察到的超声图像f,且满足超声图像成像模型f=vn,其中,乘性噪声n在图像的局部灰度均匀区域近似为Weibull分布,即噪声的概率密度函数满足:
其中,exp()为指数函数;λ和k分别为Weibull分布的尺度参数和形状参数,选取超声图像的局部灰度均匀区域,获得直方图拟合的样本数据,利用最大似然法对尺度参数λ和形状参数k进行估计。
步骤2:根据拟合的噪声分布,推导对数真实无噪图像的似然函数,将似然函数进行负的对数变换,从而得到最小化能量函数的数据保真项。
所述步骤2中,对超声图像成像模型f=vn,进行对数变换可得logf=logv+logn。令g=logf,u=logv和w=logn,则有g=u+w。根据乘性噪声n的概率密度函数pN(n)推导对数变换后的噪声w满足:
进一步,求得对数变换后的对数真实无噪图像u的似然函数为:
将似然函数进行负的对数变换,且忽略常数部分,可得:
-logpG|U(g|u)=-k(g-u)+λ-kek(g-u);
假定图像区域Ω内的所有像素g(x,y)∈Ω满足独立同分布,则能构建最小化能量函数的数据保真项如下:
步骤3:将对数真实无噪图像一阶正则化表达式和二阶正则化表达式进行加权融合,构建最小化能量函数的混合高阶正则化项,且权值能根据图像边缘特征自适应调整。
其中,η为控制函数平滑程度的阈值,取值范围为0.1~0.9,本发明实施例中取值为0.65。步骤4:将步骤2的数据保真项和步骤3的混合高阶正则化项形成最小化能量函数模型,采用Split-Bregman迭代方法进行快速求解,迭代收敛后,得到对数去噪超声图像,最后通过指数变换得到去噪超声图像。
所述步骤4中,将步骤2的数据保真项和步骤3的混合高阶正则化项组合,引入数据保真项权系数α,得到最小化能量函数模型为:
采用Split-Bregman迭代方法,对最小化能量函数模型进行求解,具体过程如下:
其中,θ1、θ2是算法求解时引入的调节系数。
b):为求u的最优解,首先将函数对u求偏导,得出欧拉-拉格朗日方程:
令偏导为0,整理后有:
其中,div(·)、div2(·)分别为一、二阶散度算子。
根据图像偏微分以及网格离散差分运算,有:
其中,上标“+”代表前向差分格式,上标“-”代表后向差分格式。
将上述离散化表达式,代入偏微分方程中,利用Gauss-Seidel迭代算法,可以得出:
c):根据b)中得到的结果,更新步骤3中的自适应加权系数γ1、γ2。
d):求解最小化能量函数模型中分离变量q,s,参考步骤b),分别将函数对q,s求偏导并置零:
采用阈值收缩算法求得:
e):得到qm+1,sm+1后,即可对迭代参量进行更新:
f):得到第m+1次迭代图像um+1后,根据前后两次迭代的滤波图像均方差小于给定的阈值,即||um+1-um||≤ξ来判定迭代收敛,阈值ξ典型取值为0.01;迭代收敛的图像um+1为对数变换后的对数去噪超声图像,最后利用指数变换得到去噪超声图像v,即v=exp(um+1)。
图2为本发明实施例所用的肌肉组织超声图像,可以看到由于散斑噪声的影响,肌肉组织的细节比较模糊。图3为本发明实施例中的肌肉组织超声图像去噪结果,可以看到去噪后的图像有效减弱了散斑噪声,同时保留了肌肉组织结构以及纹理细节,视觉上获得有效增强。
由图4~图6可以看出本发明方法在三种不同强度噪声情况下,模型中的数据保真项权系数α对模型影响较大,噪声越强α取值应该适当增大以增大数据保真项权重;调节系数θ1、θ2对模型影响较小。
图8(2)~图8(7)为不同超声图像去噪方法的处理结果对比图:包括Exp(Exponential Model,指数模型)、MIDAL(Multiplicative Image Denoising byAugmented Lagrangian,增广拉格朗日乘性图像去噪)、NTV(Nakagami Total Variation,Nakagami变分去噪)、Huang-Yang(Huang-Yang变分模型)、TM(Taylor’s expansionapproximation Model,泰勒展开近似模型)以及本发明的WHHVM(Weibull-based HybridHigh-order Variational Model,)基于Weibull分布的混合高阶变分模型方法。由图8(2)~图8(7)可以看出:MIDAL、Huang-Yang和TM的滤波结果仍有少量的散斑噪声残留,Exp和NTV处理出现了伴生的伪影现象,而本发明WHHVM方法去噪同时保留了完整的细节,具有最好的视觉效果。
Claims (6)
1.基于Weibull分布的混合高阶变分超声图像去噪方法,其特征在于包括以下步骤:
步骤1:对超声图像进行噪声分布拟合,采用Weibull概率密度函数,对选定的图像局部均匀区域进行灰度直方图拟合,并对拟合分布参数采用最大似然法进行估计;
步骤2:根据拟合的噪声分布,推导对数真实无噪图像的似然函数,将似然函数进行负的对数变换,从而得到最小化能量函数的数据保真项;
步骤3:将对数真实无噪图像一阶正则化表达式和二阶正则化表达式进行加权融合,构建最小化能量函数的混合高阶正则化项,且权值能根据图像边缘特征自适应调整;
步骤4:将步骤2的数据保真项和步骤3的混合高阶正则化项形成最小化能量函数模型,采用Split-Bregman迭代方法进行快速求解,迭代收敛后,得到对数去噪超声图像,最后通过指数变换得到去噪超声图像。
6.根据权利要求5所述基于Weibull分布的混合高阶变分超声图像去噪方法,其特征在于:采用Split-Bregman迭代方法,对最小化能量函数模型进行求解,具体过程如下:
其中,θ1、θ2是算法求解时引入的调节系数;
b):为求u的最优解,首先将函数对u求偏导,得出欧拉-拉格朗日方程:
令偏导为0,整理后有:
其中,div(·)、div2(·)分别为一、二阶散度算子;
根据图像偏微分以及网格离散差分运算,有:
其中,上标“+”代表前向差分格式,上标“-”代表后向差分格式;
将上述离散化表达式,代入偏微分方程中,利用Gauss-Seidel迭代算法,可以得出:
c):根据b)中得到的结果,更新步骤3中的自适应加权系数γ1、γ2;
d):求解最小化能量函数模型中分离变量q,s,参考步骤b),分别将函数对q,s求偏导并置零:
采用阈值收缩算法求得:
e):得到qm+1,sm+1后,即可对迭代参量进行更新:
f):得到第m+1次迭代图像um+1后,根据前后两次迭代的滤波图像均方差小于给定的阈值,即||um+1-um||≤ξ来判定迭代收敛,阈值ξ典型取值为0.01;迭代收敛的图像um+1为对数变换后的对数去噪超声图像,最后利用指数变换得到去噪超声图像v,即v=exp(um+1)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010550145.4A CN111815527B (zh) | 2020-06-16 | 2020-06-16 | 基于Weibull分布的混合高阶变分超声图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010550145.4A CN111815527B (zh) | 2020-06-16 | 2020-06-16 | 基于Weibull分布的混合高阶变分超声图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111815527A true CN111815527A (zh) | 2020-10-23 |
CN111815527B CN111815527B (zh) | 2023-10-27 |
Family
ID=72845240
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010550145.4A Active CN111815527B (zh) | 2020-06-16 | 2020-06-16 | 基于Weibull分布的混合高阶变分超声图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111815527B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112561820A (zh) * | 2020-12-17 | 2021-03-26 | 三峡大学 | 一种适用于超声图像去噪的自适应加权混合总变分方法 |
CN112767272A (zh) * | 2021-01-20 | 2021-05-07 | 南京信息工程大学 | 一种权值自适应混合阶全变分图像去噪算法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160225169A1 (en) * | 2013-09-30 | 2016-08-04 | Koninklijke Philips N.V. | Method for local adjustment of regularization parameters for image quality optimization in fully 3d iterative ct reconstruction |
CN107369139A (zh) * | 2017-06-27 | 2017-11-21 | 五邑大学 | 一种基于混合数据拟合与加权全变差的噪声图像去模糊方法 |
CN107705265A (zh) * | 2017-10-11 | 2018-02-16 | 青岛大学 | 一种基于总曲率的sar图像变分去噪方法 |
CN109767400A (zh) * | 2019-01-14 | 2019-05-17 | 三峡大学 | 一种导向三边滤波的超声图像散斑噪声去除方法 |
-
2020
- 2020-06-16 CN CN202010550145.4A patent/CN111815527B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160225169A1 (en) * | 2013-09-30 | 2016-08-04 | Koninklijke Philips N.V. | Method for local adjustment of regularization parameters for image quality optimization in fully 3d iterative ct reconstruction |
CN107369139A (zh) * | 2017-06-27 | 2017-11-21 | 五邑大学 | 一种基于混合数据拟合与加权全变差的噪声图像去模糊方法 |
CN107705265A (zh) * | 2017-10-11 | 2018-02-16 | 青岛大学 | 一种基于总曲率的sar图像变分去噪方法 |
CN109767400A (zh) * | 2019-01-14 | 2019-05-17 | 三峡大学 | 一种导向三边滤波的超声图像散斑噪声去除方法 |
Non-Patent Citations (2)
Title |
---|
张宏群;陈小晴;陶兴龙;: "基于混合阶偏微分方程的超声图像降噪", 信息技术, no. 01 * |
毛天祺;刘伟;黄洁;赵拥军;: "二进小波增强与边缘局部信息FCM的SAR图像变化检测", 信号处理, no. 01 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112561820A (zh) * | 2020-12-17 | 2021-03-26 | 三峡大学 | 一种适用于超声图像去噪的自适应加权混合总变分方法 |
CN112561820B (zh) * | 2020-12-17 | 2023-10-27 | 三峡大学 | 一种适用于超声图像去噪的自适应加权混合总变分方法 |
CN112767272A (zh) * | 2021-01-20 | 2021-05-07 | 南京信息工程大学 | 一种权值自适应混合阶全变分图像去噪算法 |
Also Published As
Publication number | Publication date |
---|---|
CN111815527B (zh) | 2023-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109978778B (zh) | 基于残差学习的卷积神经网络医学ct图像去噪方法 | |
CN109767400B (zh) | 一种导向三边滤波的超声图像散斑噪声去除方法 | |
Xiao et al. | Segmentation of ultrasound B-mode images with intensity inhomogeneity correction | |
Zhou et al. | High spatial–temporal resolution reconstruction of plane-wave ultrasound images with a multichannel multiscale convolutional neural network | |
EP1690230B1 (en) | Automatic multi-dimensional intravascular ultrasound image segmentation method | |
CN109961411B (zh) | 非下采样剪切波变换医学ct图像去噪方法 | |
US20060184021A1 (en) | Method of improving the quality of a three-dimensional ultrasound doppler image | |
CN109598680B (zh) | 基于快速非局部均值和tv-l1模型的剪切波变换医学ct图像去噪方法 | |
CN109003232B (zh) | 基于频域尺度平滑Shearlet的医学MRI图像去噪方法 | |
CN111815527B (zh) | 基于Weibull分布的混合高阶变分超声图像去噪方法 | |
Long et al. | Incoherent clutter suppression using lag-one coherence | |
Yao et al. | Multi-view 3D echocardiography compounding based on feature consistency | |
CN109003280B (zh) | 一种双通道血管内超声影像的血管中内膜分割方法 | |
CN110610528A (zh) | 基于模型的双约束光声断层图像重建方法 | |
CN112348794A (zh) | 基于注意力增强u型网络的超声乳腺肿瘤自动分割方法 | |
Seabra et al. | A 3-D ultrasound-based framework to characterize the echo morphology of carotid plaques | |
Zhang et al. | Despeckling Methods for Medical Ultrasound Images | |
CN111340805A (zh) | B超图像处理装置、脂肪肝b超图像处理装置及b超系统 | |
CN117115452A (zh) | 可控的医学超声图像去噪方法、系统及计算机存储介质 | |
Yu et al. | Generative adversarial networks with dense connection for optical coherence tomography images denoising | |
CN112561820B (zh) | 一种适用于超声图像去噪的自适应加权混合总变分方法 | |
CN109584322B (zh) | 基于频域方向平滑Shearlet医学PET图像去噪方法 | |
Hasan et al. | Blind deconvolution of ultrasound images using $ l_ {1} $-norm-constrained block-based damped variable step-size multichannel LMS algorithm | |
CN109035156B (zh) | 基于dnst的医学ct图像去噪方法 | |
Yao et al. | Spatial compounding of large numbers of multi-view 3D echocardiography images using feature consistency |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |