MATLAB中deconvwnr函数用法

news/2024/7/21 5:09:06 标签: 算法, matlab, 图像处理

目录

语法

说明

示例

使用 Wiener 滤波对图像进行去模糊处理


        deconvwnr函数的功能是使用 Wiener 滤波对图像进行去模糊处理。

语法

J = deconvwnr(I,psf,nsr)
J = deconvwnr(I,psf,ncorr,icorr)
J = deconvwnr(I,psf)

说明

J = deconvwnr(I,psf,nsr) 使用 Wiener 滤波算法对图像 I 进行反卷积,从而返回去模糊后的图像 J。psf 是对 I 进行卷积的点扩散函数 (PSF)。nsr 是加性噪声的噪信功率比。在估计图像与真实图像之间的最小均方误差意义上,该算法是最优的。

J = deconvwnr(I,psf,ncorr,icorr) 对图像 I 进行反卷积,其中 ncorr 是噪声的自相关函数,icorr 是原始图像的自相关函数。

J = deconvwnr(I,psf) 使用 Wiener 滤波算法对图像 I 进行反卷积,无估计噪声。在不含噪情况下,Wiener 滤波等效于理想的逆滤波。

示例

使用 Wiener 滤波对图像进行去模糊处理

        将图像读入工作区并显示它。

I = im2double(imread('cameraman.tif'));
imshow(I);
title('Original Image (courtesy of MIT)');

如图所示:

仿真运动模糊。

LEN = 21;
THETA = 11;
PSF = fspecial('motion', LEN, THETA);
blurred = imfilter(I, PSF, 'conv', 'circular');
figure, imshow(blurred)

如图所示:

仿真加性噪声。

noise_mean = 0;
noise_var = 0.0001;
blurred_noisy = imnoise(blurred, 'gaussian', ...
                        noise_mean, noise_var);
figure, imshow(blurred_noisy)
title('Simulate Blur and Noise')

如图所示:

        尝试在假设没有噪声的情况下进行还原。

estimated_nsr = 0;
wnr2 = deconvwnr(blurred_noisy, PSF, estimated_nsr);
figure, imshow(wnr2)
title('Restoration of Blurred, Noisy Image Using NSR = 0')

        如图所示:

尝试使用更好的噪信功率比估计值进行还原。

estimated_nsr = noise_var / var(I(:));
wnr3 = deconvwnr(blurred_noisy, PSF, estimated_nsr);
figure, imshow(wnr3)
title('Restoration of Blurred, Noisy Image Using Estimated NSR');

如图所示:

参数说明

I — 模糊图像

        模糊图像,指定为任意维度的数值数组。

psf — 点扩散函数

        点扩散函数,指定为数值数组。

nsr — 噪信比

        噪信比,指定为非负标量或与图像 I 大小相同的数值数组。如果 nsr 是数组,则它表示频谱域。为 nsr 指定 0 等效于创建一个理想的逆滤波器。

ncorr — 噪声的自相关函数

噪声的自相关函数,指定为任意大小或维度(不超过原始图像)的数值数组。

  • ​如果 ncorr 的维度与图像 I 的维度匹配,则值对应于每个维度内的自相关。

  • ​如果 ncorr 是向量,psf 也是向量,则 ncorr 中的值表示第一个维度的自相关函数。

  • 如果 ncorr 是向量,psf 是数组,则一维自相关函数通过对称性外插至 psf 的所有非单一维度。

  • 如果 ncorr 是标量,则值表示噪声的功率。

icorr — 图像的自相关函数

图像的自相关函数,指定为任意大小或维度(不超过原始图像)的数值数组。

  • ​如果 icorr 的维度与图像 I 的维度匹配,则值对应于每个维度内的自相关。

  • ​如果 icorr 是向量,psf 也是向量,则 icorr 中的值表示第一个维度的自相关函数。

  • 如果 icorr 是向量,psf 是数组,则一维自相关函数通过对称性外插至 psf 的所有非单一维度。

  • 如果 icorr 是标量,则值表示图像的功率。

J — 去模糊后的图像

去模糊后的图像,以数值数组形式返回。J 与 I 的数据类型相同。

提示

  • ​输出图像 J 可能出现算法中使用的离散傅里叶变换引入的振铃效应。要减少振铃效应,请在调用 deconvwnr 之前使用 I = edgetaper(I,psf)。

参考

[1] Gonzalez, R. C., and R. E. Woods. Digital Image Processing. Addison-Wesley Publishing Company, Inc., 1992.


http://www.niftyadmin.cn/n/5170215.html

相关文章

目标检测YOLO实战应用案例100讲-基于多目标追踪的交通场景异常检测

目录 前言 国内外研究现状 目标检测 目标追踪 异常检测

【Java、MongoDB】程序控制非关系数据库

步骤&#xff1a; &#xff08;1&#xff09;连接 连接字符串 &#xff08;2&#xff09;CRUD 类与接口 解析 &#xff08;3&#xff09;maven管理方法 依赖 <dependency><groupId>org.mongodb</groupId><artifactId>mongodb-driver-legacy<…

全文索引前缀索引

什么是全文索引&#xff1f; 全文索引是一种用于支持文本数据搜索的索引类型&#xff0c;它可以帮助数据库系统高效地执行全文搜索操作&#xff0c;包括模糊匹配、关键词搜索等。在MySQL中&#xff0c;全文索引通常用于对文本类型的字段进行高效的搜索。 全文索引与传统的B-t…

LeetCode(3)删除有序数组中的重复项【数组/字符串】【简单】

目录 1.题目2.答案3.提交结果截图 链接&#xff1a; 26. 删除有序数组中的重复项 1.题目 给你一个 非严格递增排列 的数组 nums &#xff0c;请你** 原地** 删除重复出现的元素&#xff0c;使每个元素 只出现一次 &#xff0c;返回删除后数组的新长度。元素的 相对顺序 应该保…

uniapp 本身就是一个坑,里面还有无数的小坑

uniapp 本身就是一个坑&#xff0c;里面还有无数的小坑 uniapp 用起来就感觉是个很无章法的框架&#xff0c;并且里面会有无数你意想不到的小坑&#xff0c;坑中套坑。 还不如 Capacitor 好用&#xff0c;至少它没有那么多坑&#xff0c;一次性就能成功。 遇到的那些坑 初始…

群晖 DSM 7.0 Synology Photos IOS | Android 客户端下载

安卓版本低-官方说明&#xff1a;https://www.synology.com/zh-tw/dsm/7.1/software_spec/synology_photos Synology Photos Android 客户端下载大全&#xff1a; https://archive.synology.cn/download/Mobile 官网-Synology Photos | 群暉科技 Synology Inc. 全新的 DSM…

万物社用户运营工具:无代码开发下的电商平台和CRM集成

简介&#xff1a;万物社与集简云的引领式连接 万物社&#xff0c;隶属于厦门头号云信息科技有限公司&#xff0c;是一家专注于互联网和相关服务的企业。在日常的业务运营中&#xff0c;万物社通过与集简云的无代码集成&#xff0c;实现了业务流程的自动化和智能化&#xff0c;…

React 递归手写流程图展示树形数据

需求 根据树的数据结构画出流程图展示&#xff0c;支持新增前一级、后一级、同级以及删除功能&#xff08;便于标记节点&#xff0c;把节点数据当作label展示出来了&#xff0c;实际业务中跟据情况处理&#xff09; 文件结构 初始数据 [{"ticketTemplateCode": &…