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]); %采用二维中值滤波函数medfilt2对受椒盐噪声干扰的图像滤波 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; % 二阶巴特沃斯(Butterworth)滤波器 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');
- m0_746783182023-09-01这个资源值得下载,资源内容详细全面,与描述一致,受益匪浅。
- m0_736411962022-11-16资源不错,内容挺好的,有一定的使用价值,值得借鉴,感谢分享。
- 粉丝: 408
- 资源: 5624
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- Surfer,线性函数
- MyBatis 的动态 SQL 是其核心特性之一.txt
- 时代的sdddsddsddsd
- 基于哈希链表的简单人员信息管理系统
- 其他类别JdonFramework开源框架 v5.1 Build20071025-jdonframework-5.1.rar
- 2001~2022年上市公司数字赋能指数.dta
- 2001~2022年上市公司数字赋能指数.xlsx
- 信息办公石大在线财务管理系统(含源码)-shidacaiwu.rar
- 信息办公电信计费系统完整代码-netctossconformity.rar
- matlab实现TD-SCDMA中初始同步捕捉DwPTS下行同步导频时隙的仿真.zip