mattlab《数字图像处理》实验报告02
更新时间:2024-03-02 22:26:01 阅读量: 综合文库 文档下载
GDOU-B-11-112
广东海洋大学学生实验报告书(学生用表)
实验名称 图像的空域与频域增强 学院(系) 学生姓名
课程名称 数字图像处理
实验地点
课程号 19242504 班级
实验日期 12.10.19
学号
专业
一、实验目的
1.熟练运用Matlab中的函数进行图像的空域增强和频域增强
2.能自己编程实现Matlab中的二维滤波函数filter2和中值滤波函数medfilt2的部分功能
二、实验内容 必做:
1. 对图像’eight.tif’分别添加椒盐噪声和高斯噪声,分别用中值滤波和空域低通滤波对噪声图像进行去噪。
2. 对图像’saturn.tif’分别添加椒盐噪声和高斯噪声,分别用理想低通滤波和巴特沃斯低通滤波对噪声图像进行去噪。
3. 使用梯度算子和空域高通滤波对图像'moon.tif'进行锐化。 选做:
4. 自己编写一个中值滤波函数mymedfilt2,实现Matlab中的medfilt2函数的功能(边缘采用复制的方式);对图像‘eight.tif’加上椒盐噪声,用自编的函数对其进行中值滤波。
三、程序清单与运行结果 1. 添加椒盐噪声和高斯噪声 I=imread('eight.tif');
J=imnoise(I,'salt & pepper',0.02); K=imnoise(I,'gaussian');
subplot(231),imshow(I);title('原图像');
subplot(232),imshow(J);title('椒盐噪声图像'); subplot(233),imshow(K);title('高斯噪声图像'); k1=medfilt2(J,[5 5]); k2=medfilt2(K,[5 5]);
subplot(234),imshow(k1);title('椒盐中值滤波'); subplot(235),imshow(k2);title('高斯中值滤波'); k1=filter2(fspecial('gaussian',3),K); k2=filter2(fspecial('gaussian',5),K); k3=filter2(fspecial('gaussian',7),K);
subplot(336);imshow(uint8(k1));title('3*3模板空域低通滤波'); subplot(337);imshow(uint8(k2));title('5*5模板空域低通滤波'); subplot(338);imshow(uint8(k3));title('7*7模板空域低通滤波');
2. I=imread('pout.tif');
J=imnoise(I,'salt & pepper',0.02); K=imnoise(I,'gaussian',0.02);
subplot(231),imshow(I);title('原图像');
subplot(232),imshow(J);title('椒盐噪声图像'); subplot(233),imshow(K);title('高斯噪声图像'); A=double(J); f1=fft2(J);
g1=fftshift(f1); [M,N]=size(f1); n=3; d0=20;
n1=floor(M/2); n2=floor(N/2);
for i=1:M for j=1:N
d=sqrt((i-n1)^2+(j-n2)^2); h=1/(1+(d/d0)^(2*n)); g1(i,j)=h*g1(i,j); end end
g1=ifftshift(g1);
g1=uint8(real(ifft2(g1))); subplot(234);
imshow(g1);title('巴特沃斯低通滤波去椒盐');
B=double(J); f2=fft2(J);
g2=fftshift(f1); [M,N]=size(f2); d0=40;
n1=floor(M/2); n2=floor(N/2); for i=1:M for j=1:N
d=sqrt((i-n1)^2+(j-n2)^2); if d<=d0 h=1; else h=0; end
g2(i,j)=h*g2(i,j); end end
g2=ifftshift(g2);
g2=uint8(real(ifft2(g2))); subplot(235);
imshow(g2);title('理想低通滤波去椒盐');
(没有图像’saturn.tif’,用这个代替)
3. 梯度算子进行锐化。 I=imread('moon.tif'); H=fspecial('sobel'); TH=filter2(H,I);
subplot(131),imshow(I); subplot(132),imshow(TH);
空域高通滤波对图像'moon.tif'进行锐化 I=imread('moon.tif');
subplot(121),imshow(I);title('原图像'); I=double(I);
[Gx,Gy]=gradient(I); G=sqrt(Gx.*Gx+Gy.*Gy); J=I;
K=find(G<=7); J(K)=0;
Q=find(G>=7); J(Q)=255;
subplot(122),imshow(J);title('拉普拉斯算子法');
4. h=3;
I=imread('eight.tif');
subplot(231),imshow(I);title('原图像'); noisy=imnoise(I,'salt & pepper',0.05); imshow(noisy);
k1=medfilt2(noisy,[h h]);
subplot(231),imshow(k1);title('3*3'); h=5;
k2=medfilt2(noisy,[h h]);
subplot(232),imshow(k2);title('5*5'); h=7;
k3=medfilt2(noisy,[h h]);
subplot(233),imshow(k3);title('7*7');
5.
四、实验的收获与体会
成绩 指导教师
日期
注:请用A4纸书写,不够另附纸。
第 页,共 页
正在阅读:
孔雀东南飞- 课后练习及答案01-19
化工企业节能与地区循环经济发展10-19
仁爱版九年级英语上册词组归纳复习05-23
15附件V 电梯驱动主机型式试验细则05-15
20以内退位减法-教案-完整版06-27
旅游局节能减排的工作报告02-25
2010年高考试题政治分类解析:第四单元认识社会与价值选择(必修4) - 图文12-21
- 多层物业服务方案
- (审判实务)习惯法与少数民族地区民间纠纷解决问题(孙 潋)
- 人教版新课标六年级下册语文全册教案
- 词语打卡
- photoshop实习报告
- 钢结构设计原理综合测试2
- 2014年期末练习题
- 高中数学中的逆向思维解题方法探讨
- 名师原创 全国通用2014-2015学年高二寒假作业 政治(一)Word版
- 北航《建筑结构检测鉴定与加固》在线作业三
- XX县卫生监督所工程建设项目可行性研究报告
- 小学四年级观察作文经典评语
- 浅谈110KV变电站电气一次设计-程泉焱(1)
- 安全员考试题库
- 国家电网公司变电运维管理规定(试行)
- 义务教育课程标准稿征求意见提纲
- 教学秘书面试技巧
- 钢结构工程施工组织设计
- 水利工程概论论文
- 09届九年级数学第四次模拟试卷
- 图像处理
- mattlab
- 实验
- 数字
- 报告
- 三年级语文第三单元教案 - 图文
- 学校食堂整改措施
- 初三数学下第2章《简单事件的概率》教材分析
- 计算机系教师队伍现状分析与建设对策
- 2017年中国乳胶床垫行业发展现状分析及市场供需预测报告(目录)
- 2011计算机一级考试复习 - Dreamweaver
- 职教新梦,筑巢引凤 - 打造公共实训基地 服务区域经济发展
- 李正双同志2014年度考核优秀材料(定稿)
- 人教版语文二年级《清澈的湖水》教学设计
- 如今家乡山连山音乐课教案
- 学前教育大学生实习心得体会
- 病案1
- 商务文书例文-全套
- 大学校园才艺展示落地大赛活动策划完整方案
- GENEX+Probe使用指导书 - 图文
- 在返乡农民工座谈会上的讲话
- 高中英语 必修4 Unit4《Body Language》 精品学案6
- 消防主机说明及使用规程
- 最新计量经济学参考文献大全
- 第二十五章&nbsp;&nbsp;概率初步(第二课时)