opencv图像中频域滤波去除噪声,使用c++代码

时间: 2023-11-14 17:22:26 浏览: 41
以下是使用C语言和OpenCV库实现频域滤波去除噪声的示例代码: ``` #include <opencv2/opencv.hpp> #include <iostream> using namespace std; using namespace cv; int main(int argc, char** argv) { // 读取图像 Mat img = imread("lena.jpg", IMREAD_GRAYSCALE); if (img.empty()) { cout << "读取图像失败!" << endl; return -1; } // 添加高斯噪声 Mat noisyImg; GaussianBlur(img, noisyImg, Size(7, 7), 10); imwrite("noisyImage.jpg", noisyImg); // 傅里叶变换 Mat planes[] = { Mat_<float>(noisyImg), Mat::zeros(noisyImg.size(), CV_32F) }; Mat complexImg; merge(planes, 2, complexImg); dft(complexImg, complexImg); // 显示频谱图 Mat magImg; magnitude(planes[0], planes[1], magImg); magImg += Scalar::all(1); log(magImg, magImg); normalize(magImg, magImg, 0, 1, NORM_MINMAX); imshow("频谱图", magImg); // 生成滤波器 Mat filter = Mat::zeros(noisyImg.size(), CV_32F); int cx = noisyImg.cols / 2; int cy = noisyImg.rows / 2; for (int i = 0; i < noisyImg.rows; i++) { for (int j = 0; j < noisyImg.cols; j++) { float d = sqrt(pow(i - cy, 2) + pow(j - cx, 2)); if (d < 30) { filter.at<float>(i, j) = 1; } } } // 显示滤波器 imshow("滤波器", filter); // 应用滤波器 Mat filteredImg; mulSpectrums(complexImg, filter, complexImg, 0); idft(complexImg, filteredImg, DFT_SCALE | DFT_REAL_OUTPUT); // 显示去噪后的图像 imshow("去噪后的图像", filteredImg); waitKey(); return 0; } ``` 以上代码实现了以下功能: 1. 读取图像; 2. 给图像添加高斯噪声; 3. 进行傅里叶变换,并显示频谱图; 4. 生成滤波器; 5. 应用滤波器; 6. 进行反傅里叶变换,得到去噪后的图像; 7. 显示去噪后的图像。 在代码中,我们使用`imread`函数读取图像,使用`GaussianBlur`函数添加高斯噪声。然后,我们将图像转换为频域表示,通过计算距离生成一个低通滤波器,将其应用于频域图像,再进行反傅里叶变换得到去噪后的图像。最后,我们使用`imshow`函数显示频谱图、滤波器和去噪后的图像,并使用`waitKey`函数等待用户按下键盘上的任意键。 注意:这里的代码只适用于灰度图像。如果要处理彩色图像,需要将其转换为YUV或HSV等颜色空间,对亮度通道进行处理,再转换回RGB。

相关推荐

最新推荐

recommend-type

opencv 图像滤波(均值,方框,高斯,中值)

主要介绍了opencv 图像滤波(均值,方框,高斯,中值),文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
recommend-type

opencv实现轮廓高斯滤波平滑

主要为大家详细介绍了opencv实现轮廓高斯滤波平滑,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
recommend-type

使用Python和OpenCV检测图像中的物体并将物体裁剪下来

主要介绍了使用Python和OpenCV检测图像中的物体并将物体裁剪下来,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
recommend-type

opencv3/C++图像边缘提取方式

今天小编就为大家分享一篇opencv3/C++图像边缘提取方式,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
recommend-type

opencv+python实现均值滤波

主要为大家详细介绍了opencv+python实现均值滤波,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

实现实时数据湖架构:Kafka与Hive集成

![实现实时数据湖架构:Kafka与Hive集成](https://img-blog.csdnimg.cn/img_convert/10eb2e6972b3b6086286fc64c0b3ee41.jpeg) # 1. 实时数据湖架构概述** 实时数据湖是一种现代数据管理架构,它允许企业以低延迟的方式收集、存储和处理大量数据。与传统数据仓库不同,实时数据湖不依赖于预先定义的模式,而是采用灵活的架构,可以处理各种数据类型和格式。这种架构为企业提供了以下优势: - **实时洞察:**实时数据湖允许企业访问最新的数据,从而做出更明智的决策。 - **数据民主化:**实时数据湖使各种利益相关者都可
recommend-type

解释minorization-maximization (MM) algorithm,并给出matlab代码编写的例子

Minorization-maximization (MM) algorithm是一种常用的优化算法,用于求解非凸问题或含有约束的优化问题。该算法的基本思想是通过构造一个凸下界函数来逼近原问题,然后通过求解凸下界函数的最优解来逼近原问题的最优解。具体步骤如下: 1. 初始化参数 $\theta_0$,设 $k=0$; 2. 构造一个凸下界函数 $Q(\theta|\theta_k)$,使其满足 $Q(\theta_k|\theta_k)=f(\theta_k)$; 3. 求解 $Q(\theta|\theta_k)$ 的最优值 $\theta_{k+1}=\arg\min_\theta Q(
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。