matlabe基于残差学习的图像去噪_《利用组稀疏残差约束进行图像去噪》

2023-10-29

题目:《利用组稀疏残差约束进行图像去噪》

来源:Zha Z, Liu X, Zhou Z, et al. Image denoising via group sparsity residual constraint[C]. 2017 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP). IEEE, 2017: 1787-1791.

作者:Zhiyuan Zha1,3, Xin Liu2, Ziheng Zhou2, Xiaohua Huang2, Jingang Shi2, Zhenhong Shang3, Lan Tang1, Yechao Bai1, Qiong Wang1, Xinggan Zhang1,3

作者单位:1 南京大学电子科学与工程学院;2 机器视觉和信号分析中心,芬兰奥卢大学;3 昆明理工大学信息工程与自动化学院

1.摘要

组稀疏残差在各种图像复原任务中显示出巨大的潜力,对此本文作者提出了一种通过组稀疏残差约束进行图像去噪的先验模型。为了提高组稀疏图像的去噪性能,提出了组稀疏残差的概念,将图像去噪问题转化为降低组稀疏残差的问题。为了能够减少残差,本文作者首先通过对含噪图像进行估计来获得对原始图像的组稀疏系数的一些好的估计,然后将噪声图像的组稀疏系数集中到估计中。实验结果表明,本文算法的去噪结果不仅优于已有的BM3D、WNNM等图像处理算法,而且具有更高的效率。

2. 组稀疏残差约束建模

2.1 基于组稀疏表示

最近的研究表明[1-3],结构稀疏、组稀疏表示在图像去噪任务中能够获得更好地性能。具体地,对于一幅图像

被分割成
个尺寸为
的图像块
。这时,对于每个样本块
,在大小为
的窗口中寻找
个相似块生成一个组
,定义为
,其中
表示第
组中的第
个相似块(列向量)。相似于块稀疏表示模型,设计一个字典
,每个组
的稀疏表示可以表示为
,并通过求解下面的
范数最小化问题,

这时,整个图像就可以通过所有的组稀疏编码集合来表示。基于此,图像去噪问题就转换为利用组稀疏编码从噪声组中重建图像,表示为:

如此一来,一旦能够获得所有的组稀疏系数

,就能够通过
重建原图像

2.2 组稀疏残差约束

尽管组稀疏已经证实了在图像去噪问题上的有效性,但是由于噪声的影响,很难从噪声图像Y中估计出来真实的组稀疏系数B。因此,图像去噪的质量很大程度上取决于组稀疏残差的水平,即定义组稀疏系数A与组稀疏系数B的差值,

因此,为了减少组稀疏残差R,以及提升组稀疏系数A的精确度,作者提出了一种新的图像去噪模型,组稀疏残差约束(Group Sparsity Residual Constraint, GSRC)模型,因此能够得到下式:

2.3 GSRC算法

在上式中,处理估计B的值,也需要决定p的取值。因此作者通过一些实验来研究残差的统计特性。在这写实验中,Leave(树叶)图像被用来进行图像去噪的操作,其中原图X已经加入了方差为30的高斯白噪声。作者绘制出来了R的直方图以及R拟合的高斯,拉普拉斯以及超拉普拉斯分布,如图1(a)所示。为了更好地观察其分布情况,在图1(b)中画出在对数域的直方图。

图1 (a)R的分布图

图1(b)R在对数域的分布图

很明显,R的直方图能够更好地通过拉普拉斯分布来描述,因此使用L1范数来调整残差,利用下式对提出的模型进行求解,

通过上述分析,作者提出一种基于组稀疏残差约束的图像去噪算法,目前的去噪算法如下表所示。

3 实验结果

在这一节中,作者将本文提出的算法与现有的图像去噪算法BM3D,NCSR[4],AST-NLS[5],MSEPLL[6]进行对比。GSRC算法的参数设置如下,对相似块的搜索窗半径为30*30,同时

设置为0.0001,每一个图像块的大小设置为6*6,7*7,8*8,9*9。并且文章中的源代码能够在网站:
http://www.ee.oulu.fi/˜xliu/research/gsrc/gsrc.html.中下载得到。文中首先在图2中的12幅图像中进行实验,加入方差为30,40,50,100的高斯白噪声。

图2 所有测试图像

各种去噪方法的PSNR结果如表1所示,能够发现本文提出的GSRC具有较强的性能。它能够得到0.39-0.6dB、0.46-0.76dB、0.08-0.14dB、0.3-0.41dB和0.61-0.66dB的提升针对BM3D、NCSR、WNNM、AST-NLS和MSEPLL算法,这说明本文提出的GSRC算法能够得到更好地图像去噪结果。举例来看,针对图像foremanHouse而言,当噪声的标准差为100的时候,GSRC算法的PSNR比BM3D算法在数值上要高出1.24dB和1.37dB。

图3和图4分别显示了在噪声级别为40和100时各种去噪算法的结果图像。可以发现BM3D、NCSR、WNNM、AST-NLS和MSEPLL不仅生成了一些不良的伪像,而且丢失了一些细节。相比之下,GSRC能够更好地保留图像的边缘信息。由此来看,本文提出的GSRC算法能够获得更好地图像去噪结果。

图3 不同去噪方法的lin结果图(噪声标准差为40)(a)原图;(b)噪声图像;(c)BM3D(PSNR=29.53dB);(d)NCSR(PSNR=29.45dB);(e)WNNM(PSNR=29.80dB);(f)AST-NLS(PSNR=29.43dB);(g)MSE-PLL(PSNR=29.68dB);(h)GSRC(PSNR=29.98dB)

图4 不同去噪方法的lin结果图(噪声标准差为40)(a)原图;(b)噪声图像;(c)BM3D(PSNR=25.87dB);(d)NCSR(PSNR=25.57dB);(e)WNNM(PSNR=26.72dB);(f)AST-NLS(PSNR=26.92dB);(g)MSE-PLL(PSNR=25.99dB);(h)GSRC(PSNR=27.24dB)

效率是评估算法的另一个关键因素,因此作者比较了上文提到的图像去噪方法的速度。所有实验都是在Matlab 2012b环境下,在装有Intel(R)Core(TM)i3-4150、3.56Hz CPU和4GB内存的计算机上进行的。表2显示了各种去噪方法在测试图像上的运行时间。可以看出,与BM3D相比,本文提出的GSRC使用的计算时间比其他所有方法(除了BM3D)都要少,这可能是因为BM3D是通过C++实现的,而GSRC是在Matlab中实现的。

4 结论

在本文中,作者提出了一种新的先验模型,用于通过组稀疏残差约束进行图像去噪。为了提高组稀疏图像去噪的性能,提出了组稀疏残差的概念,从而将图像去噪问题转化为减少组稀疏残差的问题。为此,作者首先通过噪声图像的初步估计来获得对原始图像的组稀疏系数的一些好的估计,然后将噪声图像的组稀疏系数集中到估计中。实验结果表明,在PSNR和时间方面,该方法优于许多最新的去噪方法,例如BM3D和WNNM图像去噪算法。

5 思考

本文作者提出一种基于组稀疏残差约束的图像去噪算法。不同于传统的基于单一非局部自相关先验信息方法,2种非局部自相关先验信息被用于图像去噪,即噪声图像的非局部自相关先验信息以及自然图像的非局部自相关先验信息。为了提升图像去噪的性能,提出了组稀疏残差的概念,将图像去噪问题转化为如何减少组稀疏残差,由于自然图像组包含了许多非局部自相关先验信息,作者首先利用高斯混合模型从自然图像中学习原图像的组稀疏估计;其次,噪声图像的组稀疏系数被用来逼近原图像的组稀疏估计。实验结果表明提出的图像去噪算法相比于现有的图像去噪方法,可以获得更好的去噪性能。

6 参考文献

[1] Julien Mairal, Francis Bach, Jean Ponce, Guillermo Sapiro, and Andrew Zisserman, “Non-local sparse models for image restoration,” in 2009 IEEE 12th International Conference on Computer Vision. IEEE, 2009, pp. 2272–2279.

[2] Kostadin Dabov, Alessandro Foi, Vladimir Katkovnik, and Karen Egiazarian, “Image denoising by sparse 3-d transform domain collaborative filtering,” IEEE Transactions on image processing, vol. 16, no. 8, pp. 2080–2095, 2007.

[3] Hui Ji, Sibin Huang, Zuowei Shen, and Yuhong Xu, “Robust video restoration by joint sparse and low rank matrix approxi mation,” SIAM Journal on Imaging Sciences, vol. 4, no. 4, pp. 1122–1142, 2011.

[4] Weisheng Dong, Lei Zhang, Guangming Shi, and Xin Li, “Nonlocally centralized sparse representation for image restoration,” IEEE Transactions on Image Processing, vol. 22, no. 4, pp. 1620–1630, 2013.

[5] Hangfan Liu, Ruiqin Xiong, Jian Zhang, and Wen Gao, “Image denoising via adaptive soft-thresholding based on non-local samples,” in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, 2015, pp. 484–492.

[6] Vardan Papyan and Michael Elad, “Multi-scale patch-based image restoration,” IEEE Transactions on Image Processing, vol. 25, no. 1, pp. 249–261, 2016.

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

matlabe基于残差学习的图像去噪_《利用组稀疏残差约束进行图像去噪》 的相关文章

  • 【Proteus仿真】555组成的多谐振荡器电路

    Proteus仿真 555组成的多谐振荡器电路 Proteus仿真演示 多谐振荡器电路 多谐振荡器电路是一种矩形波产生电路 属于数字电路 三极管不工作在放大线性区 这种电路不需要外加触发信号便能连续地 周期性地自行产生矩形脉冲 该脉冲是由基
  • Stable Diffusion:ChatGPT与AI绘画,引领艺术的未来

    人工智能 AI 的快速发展正在为各个领域带来革命性的变化 其中包括艺术与创意领域 AI绘画是一种将人工智能技术与艺术创作相结合的新兴范式 通过深度学习和生成对抗网络 GAN 等技术 AI绘画可以生成各种富有创意和想象力的艺术作品 本文将探讨
  • python - __str__ 和 __repr__

    内建函数str 和repr representation 表达 表示 或反引号操作符 可以方便地以字符串的方式获取对象的内容 类型 数值属性等信息 str 函数得到的字符串可读性好 故被print调用 而repr 函数得到的字符串通常可以用
  • Docker+docker-compose+nginx部署已有项目

    项目背景 在异地服务器拷docker相关项目到新的服务器 具体操作 1 新服务器安装好docker 2 新服务器安装好docker compose 3 从老服务器拷贝镜像到新服务器 4 新服务器导入镜像 5 构建项目地址挂载目录 找到doc
  • 用U盘作启动盘装Windows10系统整套流程 纯净版(不用其他乱七八糟的软件)(macOS适用)

    简介 本人的电脑是MacBook Air 2014年版的 因为内存小而且文件杂乱 所以一下子都给格式化了 但是要用Mac自带的恢复系统的话需要连接校园网 连接校园网又需要打开网页输入账号和密码 我们学校的校园网是这样的 所以只能用U盘作为格
  • gqrx编译过程记录

    gqrx编译过程记录 目标 环境 编译 下载源代码 建立编译位置 修改CMakefile txt中的模块 编译安装 运行界面 没有更多 目标 在ubuntun 20 04桌面版编译gqrx 通过USRP 205mini实现收音机功能 环境
  • 【解决】docker容器怎么使用宿主机的IPv6地址

    在IPv4时代 我们对外访问都是端口映射 都没有公网IP 但是在IPv6时大家都有公网IP 可能需要容器地址和主机地址一致 可以在docker run时使用参数 network host 则此容器网络和宿主机一致 docker run ne
  • AQS详解

    AQS详解 文章目录 AQS详解 AQS简单介绍 AQS原理 AQS原理概览 AQS对资源的共享方式 AQS定义两种资源共享方式 Exclusive 独占 Share 共享 AQS底层使用了模板方法模式 Semaphore 信号量 Coun
  • 浅谈可重入锁

    一 可重入锁 递归锁 1 概念 同一个线程在外层方法获取锁的时候 再进入该线程的内层方法会自动获取锁 前提是 锁对象是同一个对象 不是因为之前已经获取过还没有释放而阻塞 2 java中的ReentrantLock和synchronied都是
  • 关于Gdi+和GdiplusStartup

    GDI 实际上是一组类的定义 封装了gdi 的几乎所有API 当然使用方法就要从这些 例子 里边寻找了 本文正是尝试用GDI 写一个纯SDK的程序 语言自然是我最喜欢的语言WIN32ASM 这个程序很简单 就是用GDI 画了一条直线 算是抛
  • HCIA-FusionCompute华为企业级虚拟化

    一 云计算 按需付费 集中资源对外提供服务 1 云本身没有资源 云是资源整合者 整合底层的所有计算机资源 cpu 内存 磁盘等 云计算是一种模型 它可以实现随时随地 随需应变地从可配置计算资源共享池中获取所需的资源 例如 网络 服务器 存储
  • BigDecimal 问题小结

    BigDecimal 加法 add 函数 乘法multiply 函数 除法divide 函数 绝对值abs 函数 减法subtract 函数 ROUND CEILING 向正无穷方向舍入 ROUND DOWN 向零方向舍入 ROUND FL
  • 【Redis】新增数据结构

    BitMap位图 Redis提供了Bitmaps这个 数据类型 可以实现对位的操作 1 Bitmaps本身不是一种数据类型 实际上它就是字符串 key value 但是它可以对字符串的位进行操作 2 Bitmaps单独提供了一套命令 所以在
  • RabbitMQ与SpringBoot整合实战

    SpringBoot整合RabbitMQ SpringBoot与RabbitMQ集成非常筒単 不需要做任何的额外设置只需要两步即可 step1 引入相关依赖 spring boot starter amqp step2 対applicati
  • PyTorch-02梯度下降Gradient Descent、回归案例、手写数字识别案例

    PyTorch 02梯度下降Gradient Descent 回归案例 手写数字识别案例 了解梯度下降 梯度下降是深度学习的精髓 整个deep learning是靠梯度下降所支撑的 可以求解一个非常难的函数 使用的方法就是梯度下降算法 求一
  • 群体智能优化算法--烟花算法(附优化参数的通用代码)

    烟花算法是由北京大学谭营教授提出了烟花算法 这是一种既简单又具有较强优化能力的算法 根据烟花爆炸的原理 每个烟花爆炸之后会选择最好的烟花作为下一次爆炸的烟花 而且在多个烟花爆炸的同时 每个烟花都是相互独立的 寻找最优爆炸烟花只在自身本身爆炸
  • C语言不同操作系统不同编译器,msvc mingw gcc cmake VS MSVC的理解

    编译器的编译有三步 1 源代码生成汇编码 2 汇编语言生成中间代码 obj类型 一个源文件一个 obj 每个源文件通常编译成一个对应的目标文件 obj或 o 但在某些情况下 多个源文件可以编译成一个目标文件 3 连接 在汇编里称Link 在
  • 高并发网络编程之epoll详解

    在linux 没有实现epoll事件驱动机制之前 我们一般选择用select或者poll等IO多路复用的方法来实现并发服务程序 在大数据 高并发 集群等一些名词唱得火热之年代 select和poll的用武之地越来越有限 风头已经被epoll
  • JavaScript DOM 编程艺术学习笔记(一):静态标记

    JavaScriptDOM编程艺术 学习笔记 一 静态标记 DOM DOM脚本程序设计 涵盖了使用任何一种支持DOM API的程序设计语言去处理任何一种标记文档的情况 DOM是程序设计语言和标记文档之间的接口 它将文档表示成一棵节点树 每个

随机推荐

  • iOS系统和XCode各版本发布日期

    本人收集了iOS系统和XCode各版本发布日期 供大家参考 发布或推送日期 版本编号 更新 2023年9月7日 iOS16 6 1 推出iOS16 6 1正式版 2023年7月24日 iOS16 6 推出iOS16 6正式版 2023年6月
  • 【sybase】linux环境安装sybase数据库

    目录 安装步骤 使用root创建用户 用户组 安装目录 安装数据库 1 切换sybase用户登录 将ase157 linuxx86 64 tgz上传到 home sybase目录下 2 解压 3 执行安装启动文件 4 创建数据库server
  • vs中出现error LNK2038: 检测到“_MSC_VER”的不匹配项问题

    参考博客 https blog csdn net shenmifangke article details 50395116 例如原先是2013版本的 现在换成2015版本的话 方法 在项目 解决方案资源管理器或者属性管理器里都行 右键属性
  • oracle的时间表达,Oracle中的日期类型

    1 SYSDATE 获取当前系统时间 select SYSDATE from dual 格式化日期 TO CHAR SYSDATE YY MM DD HH24 MI SS 或 TO DATE SYSDATE YY MM DD HH24 MI
  • Oauth2源码剖析——密码式+数据库存储

    访问样例 授权服务器源码剖析 TokenEndPoint java gt postAccessToken principal parameters String clientId getclientId principal 得到Author
  • ARM机器学习新平台Trillium

    Project Trillium 主要面向机器学习和神经网络市场 平台主要有 Arm ML 处理器 OD 处理器 以及 NN 软件这三板斧 顾名思义 ML 处理器主打机器学习 特别是激动计算应用 最高性能可达 4 6 万亿次每秒 能效超过
  • Matlab 归一化(normalization)/标准化 (standarization)

    数据规范中的归一化与标准化 A 归一化 vs 标准化 归一化 要把你需要处理的数据经过处理后 通过某种算法 限制在你需要的一定范围内 首先归一化是为了后面数据处理的方便 其次是保正程序运行时收敛加快 一般指将数据限制在 0 1 之间 把数变
  • 2020年计算机、信安推免总结

    这里写自定义目录标题 个人情况 夏令营 浙软和南软 四川大学网络空间安全学院 湖南大学信科院 西工大软件学院 北京交通大学软件学院 预推免 中科院网络信息中心 浙江大学软件学院 天津大学智算学部 华中科技大学网络空间安全学院 东南大学网络空
  • 添加votedisk

    1 添加votedisk 必须有一半以上的votedisk同时可用 clusterware才能正常工作 否则cluster立刻宕掉 所以最好votedisk保持单数个 添加和删除votedisk的操作非常危险 必须在停止数据库 停止asm
  • Unity3D暂停,继续游戏,重新开始,退出,以及 UnityEditor.EditorApplication打包后不会执行

    1 暂停游戏 Time timescale 0 2 继续游戏 Time timescale 1 3 重新开始 using UnityEngine SceneManagement SceneManager LoadScene 0 其中 0 为
  • Python 中的<>和!= 区别

    今天在编写Python MySQL 采集脚本过程中 需要使用到 不等于 表达方式 第一种写法 在Python2 6以前版本 不等于 if string atof func get item mysql status Qcache hits
  • C++/C++11中头文件algorithm的使用

  • Ubuntu20.04 操作系统切换中文 & 安装中文输入法

    网上有很多的帖子指导如何安装搜狗输入法 我也尝试过 需要安装一些指令 解决依赖很麻烦 很消耗精力 也找不到简便的办法 因此还是用系统自带的输入法吧 一 换源 这一步很重要 系统下载输入法时 能够从阿里云服务器下载大大提升速度 具体请参考拙作
  • 中国数字经济发展指数报告(2023年)

    报告 显示 数字经济具有高度发展驱动性 各地区数字经济发展指数与地区综合经济实力具有高度正相关性 数字经济对于促进地区经济发展的驱动能力强劲 关注公众号 互联互通社区 回复 SZJJ396 获取全部报告内容 从总指数结果来看 我国数字经济呈
  • DDR3相关概念

    1 BL8 突发长度 理解一下概念 突发传输 在通信领域中一般指在短时间内进行相对高带宽的数据传输 Burst 突发 是指在同一行中相邻的存储单元连续进行数据传输的方式 BL突发长度 连续传输的周期数就是突发长度 在突发传输的模式下 多个数
  • tcp长连接的双向加密

    一 历史和基础 做即时通讯项目 网游登录等业务时 我们通常需要一种能在线回推的机制 这时候客户端的tcp长连接就可能派上用场了 示意图 二 问题产生 曾经做过一个项目 之前老的加密方案为 1 客户端发一个rsa公钥G1给服务端 明文 自己留
  • PaddleOCR详解和识别图片中文字

    前言 paddleocr可以说是最近ocr的主 军 家对于paddleocr的认可度是相当的 特别是最近推出的轻量级模型 可以识别近80种语 并且效率是这三种ocr工具种最 的 相同的图 paddleocr只需要2秒左右 对于多场景的 语种
  • 计算机键盘字母记忆,键盘26个字母口诀是怎样的?

    01 七 Q 碗 W 鹅 E 肉 R 汤 T 已 Y 无 W 一 Y 我 W 婆 P 爱 A 死 S 豆 D 腐 F 羹 G 回 H 家 J 快 K 乐 L 自 Z 行 X 车 C 未 V 帮 B 你 N 买 M 众所周知 计算机键盘字母是
  • MongoDB3.0+ 添加用户命令

    关键命令 db createUser user coolead pwd coolead roles role dbOwner db mycoolead C Program Files MongoDB Server 3 2 bin gt mo
  • matlabe基于残差学习的图像去噪_《利用组稀疏残差约束进行图像去噪》

    题目 利用组稀疏残差约束进行图像去噪 来源 Zha Z Liu X Zhou Z et al Image denoising via group sparsity residual constraint C 2017 IEEE Intern