摘要
模糊的超声图像经过极端通道先验去模糊后不够稀疏,导致极端通道稀疏约束可能不存在。因此,为了充分利用图像通道信息,通过增强去模糊后超声图像的稀疏性,提出一种二次稀疏极端通道先验盲超声图像去模糊算法。首先,给出了相关的理论证明和实验说明二次稀疏极端通道先验用于约束模糊超声图像的可行性;然后,充分利用暗通道和亮通道的先验信息,在交替迭代过程中采用半二次分裂方法估计中间图像和模糊核;最后,用傅里叶变换求得最终的清晰图像和模糊核。在超声图像集上的实验结果表明,本文提出的二次稀疏极端通道先验盲超声图像去模糊算法和其他超声图像去模糊方法相比更具可行性和优越性。
医学超声成像技术具有动态、三维观察器官病变、价格低廉、无需手术、对人体无害以及实时性强等优点。近年来,超声设备在临床医学中得到广泛应用,但在超声成像中会由于一些不可控因素导致扫描器官运动而引起图像模糊。随着数字图像去模糊技术的快速发展,如何利用数字图像处理技术来提高医学超声图像的质量逐渐引起人们的关注。文献[
在实际应用中,医学超声成像系统的退化过程难以精确地描述出来,所以通常采用盲图像去模糊算法。由于超声成像系统在成像过程中存在固有的问题,只能得到畸变的图像。为便于描述,本文将去模糊后的图像称为清晰图像,将图像暗通道的再次变暗称为二次稀疏暗通道,将图像亮通道的再次变亮称为二次稀疏亮通道。由于超声图像经过二次稀疏后,模糊超声图像与去模糊超声图像的差异将更加明显。本文提出了二次稀疏极端通道先验盲超声图像去模糊算法,用以提高极端先验知识的精度。本文通过均方误差(Mean squared error, MSE)来测量模糊核的精度,通过峰值信噪比(Peak signal to noise ratio, PSNR)和结构相似性指数(Structural similarity index, SSIM)对去模糊结果进行了客观评
本文主要贡献如下:(1)提出了图像二次稀疏暗亮通道理论;(2)在图像二次稀疏暗亮通道理论基础上提出二次稀疏极端通道先验;(3)将二次稀疏极端通道先验应用于模糊超声图像以展示本文所提方法去模糊的可行性和优异性。
模糊图像的数学表达式可以表示为
(1) |
式中:(为图像的范围);表示模糊图像;表示清晰图像;(为图像的范围);为模糊核;为噪声;为卷积算
定义1 对于图像而言,卷积定义为
(2) |
式中:表示模糊核的范围,其大小为;[•]表示舍入算
定义2 图像的暗通道定义为
(3) |
式中:和表示图像的像素点;是以为中心的图像模块;为第种颜色通道。如果为灰度图像,则有。
定义3 图像的二次稀疏暗通道定义为
(4) |

图1 模糊超声图像及其两种稀疏暗通道
Fig.1 Blurry ultrasound image and two kinds of sparse dark channels

图2 两种暗通道的直方图对比
Fig.2 Histogram comparisons of two kinds of dark channels
定义4 图像的亮通道定义为
(5) |
如果为灰度图
定义5 图像的二次稀疏亮通道定义为
(6) |

图3 模糊超声图像及其两种稀疏亮通道
Fig.3 Blurry ultrasound image and two kinds of sparse bright channels

图4 两种亮通道的直方图对比
Fig.4 Histogram comparisons of two kinds of bright channels
将提出的二次稀疏极端通道先验项(
(8) |
式中:表示清晰图像;表示模糊图像;表示模糊核;和分别是模糊核先验与图像的二次稀疏极端通道先验;和为权重参数。
利用范数和范数分别度量二次稀疏暗通道和二次稀疏亮通道的稀疏性。同时本文使用范数作为数据保真度函数,可使去模糊图像和模糊核的卷积输出与模糊图像相似。对于模糊核的先验信息,通过快速傅里叶变换得到模糊核的范
(9) |
式中、、和为权重参数。
在文献[

图5 本文算法流程图
Fig.5 Flow chart of the proposed algorithm
由于模型(9)中清晰图像和模糊核未知,很难直接求出其解,因此本文采用了半二次分裂方
(10) |
(11) |
为了解决子问题(10)中正则项难以计算的问题,本文使用变量分离的最小化方法。通过引入图像梯度的辅助变量,二次稀疏暗通道与二次稀疏亮通道的梯度先验辅助变量和[
(12) |
式中、、、、和为权重参数。容易看出,当这些权重参数趋近无穷大时,
(13) |
式中与为非线性运算。为了便于优化,本文将转变为以保持一致性。下面通过构造一个映射矩阵来代替运算。令,,满足
(14) |
由的构造
(15) |
用快速傅里叶变换求解
(16) |
式中:;和分别表示图像的水平和垂直微分算子;与分别表示快速傅里叶变换和逆快速傅里叶变换。
给定,对于、和的子问题分别用以下方法求解
(17) |
(18) |
(19) |
数据来自布尔诺理工大学信号处理实验室。数据库包含10名不同质量(平均质量76.5±9.7 kg)的志愿者(平均年龄27.5±3.5岁)的颈总动脉(Common carotid artery, CCA)图像。图像数据库中包含84张CCA纵切面B超图像。
为了说明本文所提算法的可行性和优异性,本节将该算法运行在CCA模糊超声图像数据集上进行测试,并与文献[

图6 模糊超声图像去模糊结果对比
Fig.6 Comparisons of deblurring results for blurry ultrasound images
为进一步说明本文所提算法对模糊超声图像进行去模糊的优异性,本文将4种算法去模糊结果进行定量评价,评价结果如

图7 4种算法在CCA数据集上去模糊结果评价指标直方图
Fig.7 Evaluation index histogram of deblurring results of four algorithms on CCA dataset
本文提出了二次稀疏极端通道先验盲超声图像去模糊方法。通过实验得知,二次稀疏模糊超声图像与二次稀疏清晰超声图像之间的差异更加明显。因此,将图像的二次稀疏暗通道和图像的二次稀疏亮通道相结合,得到二次稀疏极端通道先验,提高了极端通道先验知识的准确性。本文算法充分利用了暗亮通道先验信息,对利用暗亮通道信息进行模糊超声图像去模糊具有重要意义。由于亮通道和暗通道的二次稀疏性运算增加,使得本文算法运行时间增加。接下来将针对如何缩短本文算法运行时间展开研究。
参考文献
胡双年,李艳艳,牛玉俊.基于模糊核稀疏先验正则化图像盲复原处理[J].安徽大学学报(自然科学版), 2019, 43(6): 43-47. [百度学术]
HU Shuangnian, LI Yanyan, NIU Yujun. Image blind restoration processing based on blur kernel sparse prior regularization [J]. Journal of Anhui University (Natural Sciences), 2019, 43(6): 43-47. [百度学术]
孙成阳,毛家发,胡亚红,等.基于RGB通道自适应的维纳滤波图像去模糊研究[J].计算机测量与控制,2019,27(11): 215-219, 238. [百度学术]
SUN Chengyang, MAO Jiafa, HU Yahong, et al. Research on adaptive Wiener filtering image deblurring based on RGB channel[J]. Computer Measurement & Control, 2019, 27(11): 215-219, 238. [百度学术]
余孝源, 谢巍, 陈定权,等. 一种基于改进的暗通道先验的运动模糊核估计方法[J]. 控制与决策, 2020, 35(7): 134-140. [百度学术]
YU Xiaoyuan, XIE Wei, CHEN Dingquan, et al. A motion blur kernel estimation method based on improved dark channel prior [J]. Control and Decision, 2020, 35(7): 134-140. [百度学术]
肖宿. 基于稀疏先验正则化的图像盲去模糊[J]. 绥化学院学报, 2019, 39(2): 143-147. [百度学术]
XIAO Su. Blind image deblurring based on sparse prior regularization[J]. Journal of Suihua University, 2019, 39(2): 143-147. [百度学术]
YEOH W S, ZHANG C S, CHEN M. An iterative constrained least squares filter for ultrasound image deconvolution[C]//Proceedings of 2006 International Conference of the IEEE Engineering in Medicine and Biology Society. [S.l.]: IEEE, 2006: 2758-2761. [百度学术]
OLOFSSON T, WENNERSTROM E. Sparse deconvolution of B-scan images[J]. IEEE Transactions on Ultrasonic Ferroelectrics and Frequency Control, 2007, 54(8): 1634-1640. [百度学术]
SHIN H C, PRAGER R. Sensitivity to point spread function parameters in medical ultrasound image deconvolution[J]. Ultrasonics, 2009, 49(2): 344-357. [百度学术]
PALLADINI A, TESTONI N, MARCHI L D, et al. A reduced complexity estimation algorithm for ultrasound images deblurring[J]. Computer Methods and Programs in Biomedicine, 2009, 95(2): S4-S11. [百度学术]
ZHAO N N, WEI Q, BASARAB A, et al. Single image super-resolution of medical ultrasound images using a fast algorithm[C]//Proceedings of 2016 IEEE 13th International Symposium on Biomedical Imaging (ISBI). [S.l.]: IEEE, 2016: 473-476. [百度学术]
涂春梅, 陈国彬, 刘超. 暗像素先验的模糊图像盲复原方法[J]. 计算机工程与应用, 2020, 56(10): 213-219. [百度学术]
TU Chunmei, CHEN Guobin, LIU Chao. Blind restoration method of blurred image with dark pixel prior[J]. Computer Engineering and Application, 2020, 56(10): 213-219. [百度学术]
PAN J S, SUN D Q, PFISTER H, et al. Blind image deblurring using dark channel prior[C]//Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. [S.l.]: IEEE, 2016: 1628-1636. [百度学术]
YAN Y Y, REN W Q, GUO Y F, et al. Image deblurring via extreme channels prior[C]//Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. [S.l.]: IEEE, 2017: 4003-4011. [百度学术]
YU X Y, XIE W. Single image blind deblurring based on salient edge-structures and elastic-net regularization[J]. Journal of Mathematical Imaging and Vision, 2020. DOI: 10.1007/s10851-020-00949-6. [百度学术]
LIU Q H, SUN L P, SHAO Z G. Compound L0 regularization method for image blind motion deblurring[J]. Journal of Electronic Imaging, 2016, 25(5): 053013.1-053013.15. [百度学术]
WEN F, YING R D, LIU Y P, et al. A simple local minimal intensity prior and an improved algorithm for blind image deblurring[J]. IEEE Transactions on Circuits and Systems for Video Technology, 2020(99): 1. [百度学术]
CHEN X, ZHU Q. A fast FFT-based iterative algorithm for image deblurring with anti-reflective boundary conditions[C]// Proceedings of the 2018 International Conference on Network, Communication, Computer Engineering (NCCE 2018). [S.l.]: Atlantis Press, 2018. [百度学术]
HEHN L, TILLEY S, PFEIFFER F, et al. Blind deconvolution in model-based iterative reconstruction for CT using a normalized sparsity measure[J]. Physics in Medicine & Biology, 2019, 64(21): 215010. [百度学术]
HU Z, CHO S, WANG J, et al. Deblurring low-light images with light streaks[C]//Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. [S.l.]: IEEE, 2014: 3382-3389. [百度学术]
SHU Q L, WU C S, ZHONG Q X, et al. Alternating minimization algorithm for hybrid regularized variational image dehazing[J]. Optik, 2019, 185: 943-956. [百度学术]
CHENG K H, DU J, ZHOU H X, et al. Image super-resolution based on half quadratic splitting[J]. Infrared Physics & Technology, 2020: 103193. [百度学术]
ZHANG H P, YUAN B, DONG B, et al. No-reference blurry image quality assessment by structural similarity index[J]. Applied Sciences, 2018, 8(10): 2003. [百度学术]