%% 学习目标:一个图像处理的经典综合案例
clear;clc;close all
%% 读取图像
I=imread('1.jpg');
try
I=rgb2gray(I); %如果是RGB图像,则转化为灰度图
end
%% 添加噪声
I=imnoise(I,'salt & pepper',0.04); % 叠加密度为0.04的椒盐噪声
figure
imshow(I)
title('原图')
%% 图像去噪
I=medfilt2(I,[2 4]); %采用二维中值滤波函数对图像滤波
figure
imshow(I)
title('去噪之后的图像')
%% 图像锐化
[M,N]=size(I); % 计算图像尺寸
f=double(I); % 数据类型转换,MATLAB不支持图像的无符号整型的计算
g=fft2(f); % 二维傅立叶变换,得到频域信息
g=fftshift(g); % 0频率移到(M/2,N/2)处
n=2; % 二阶巴特沃斯滤波器
D0=3; % 巴特沃斯滤波器的D0,D0越大,保留的高频信号就越少(对于不同的图片,可以自行调节,以便得到好的效果)
for i=1:M
for j=1:N
D=sqrt((i-M/2)^2+(j-N/2)^2); %计算离0频率(M/2,N/2)的距离
h=1/(1+(D0/D)^(2*n)); %计算传递函数的值
output(i,j)=h*g(i,j); %高通滤波之后(i,j)处的值
end
end
result=ifftshift(output); %与g=fftshift(g)对应,还原回去
I=ifft2(result); %傅里叶反变换
I=uint8(real(I)); %得到滤波之后的图像(锐化之后的图像)
figure
imshow(I) %滤波后图像显示
title('锐化之后的图像')
%% 直方图均衡化处理
I=histeq(I); %得到直方图均衡化之后的图像
figure
imshow(I)
title('直方图均衡之后的图像')
%% 边缘检测
I1=edge(I,'prewitt');
I2=edge(I,'roberts');
I3=edge(I,'sobel');
figure
imshow(I1)
title('prewitt边缘检测之后的图像')
figure
imshow(I2)
title('roberts边缘检测之后的图像')
figure
imshow(I3)
title('sobel边缘检测之后的图像')
%%
%% 大仙QQ:1960009019
%% 在线教育微信公众号:大仙一品堂
%% 大鱼号:在线教育大仙一品堂
%% 一点资讯号:大仙一品堂
MATLAB图像处理:8 图像去噪 滤波 锐化 边缘检测程序.zip
版权申诉
23 浏览量
2023-05-26
09:51:21
上传
评论
收藏 43.42MB ZIP 举报
![avatar](https://profile-avatar.csdnimg.cn/2ab41b4b301b4cc6861912179f3798cc_qq_41127811.jpg!1)
![avatar-vip](https://csdnimg.cn/release/downloadcmsfe/public/img/user-vip.1c89f3c5.png)
JGiser
- 粉丝: 7642
- 资源: 5110
最新资源
- 图标生成工具 Windows 版本(支持.ico和.icns)IconBuilder v1.1
- 666662222277777
- 图标生成工具 Mac 版本(支持.ico和.icns)IconBuilder v1.1
- 求最大公约数(简单).cpp
- elasticsearch数据库下载、配置、使用案例
- springboot的概要介绍与分析
- C语言的概要介绍与分析
- 第一个较大的Android项目,基于Android平台的图书管理系统(Android studio).zip
- Cisco Packet Tracer 6.2 for Windows Instructor Version
- 使⽤pyIAST计算⽓体吸附选择性
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback-tip](https://img-home.csdnimg.cn/images/20220527035111.png)