版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、数字图像处理报告非线性平滑中值滤波一. 实验目的1. 了解统计排序滤波器在数字图象处理中的作用;2. 理解中值滤波的原理、特点、适用对象;3. 掌握中值滤波的各种方法。二. 实验原理中值滤波是一种非线性信号处理方法,与其对应的中值滤波器是一种统计排序滤波器,用象素邻域各灰度级排序的中值来代替该象素的灰度级(用一个奇数点的移动窗口,将窗口中心点的值用窗口的中值代替),一定条件下,可以克服线性滤波器(如邻域平滑滤波等)所带来的图象细节模糊,而且对滤除脉冲干扰及图象扫描噪声最为有效(中值滤波的输出与输入噪声的密度分布有关,对随机噪声的抑制能力, 中值滤波比平均值滤波要差一些;但对脉冲干扰,特别是脉冲
2、宽度小于窗口一半、相距较远的窄脉冲干扰,中值滤波 的效果较好)。设有一个一维序列fi, f2,fn ,取窗口长度(点数)为 m (m为奇数),对其进行中 值滤波,就是从输入序列中相继抽出m个数fi-v,fi-i, fi, fi+i,fi+v (其中fi为窗口中心点值,v=(m-1)/ 2),再将这m个点按其数值大小排序,取其序号为中心点的那个数作为 滤波输出。用数学公式表示为:m 1y Medfiv, ,fi, fi v i N,v 丐 (4 i)这样处理的特点是:周期小于m/2 (窗口一半)的脉冲会被抑制,而周期大于m/2 (窗口一半)的脉冲得到保留,在窗口单调增加或单调减少的序列,中值滤波
3、输出信号仍保持输入信号不变(利用这个特点,可以使中值滤波既能去除图象中的噪声,又能保持图象中一些物体的边缘)。同理,二维中值滤波可由下式表示:Yj Mqd fj ( 4 2)A式中:A为窗口; fj为二维数据序列。这样处理的特点是:当在窗口,被噪声污染的象 素不超过滤波器区域的一半时,噪声会被抑制;当在窗口,被噪声污染的象素超过滤波器区域的一半时,噪声会得到保留。二维中值滤波的窗口形状和尺寸对滤波效果影响较大,不同的图象容和不同的应用要求,往往采用不同的窗口形状和尺寸。常用的二维中值滤波窗口有线状、十字形、方形等。 窗口尺寸一般先用 3X 3,再取5X 5逐渐增大, 直到滤波效果满意为止。就一
4、般经验来讲, 对于有缓变的较长轮廓线物体的图象,采用方形窗口为宜。 对于包含有尖顶物体的图象, 用十字形窗口; 而窗口大小则以不超过图象中最小有效物体的尺寸为宜。如果图象 中点、线、尖角细节较多,则不宜采用中值滤波。对一些容复杂的图象,可以使用复合型 中值滤波, 如中值滤波的线性组合, 即将几种窗口尺寸大小和形状不同的中值滤波器复合 使用, 只要各窗口都与中心对称, 滤波输出可保持几个方向上的边缘跳变, 而且跳变幅度 可调节;其线性组合方程如下:NYijak Med fij (43)k 1Ak式中:ak为不同中值滤波的系数,Ak为窗口。三. 实验容1. 分别对含有高斯噪声、瑞利噪声、椒盐噪声的
5、图象进行中值滤波,采用多种窗口形状(线形、十字形、方形),窗口尺寸可变(先用 3X 3,再取5X 5逐渐增大);2. 比较中值滤波对含有不同噪声的图象的平滑效果,讨论中值滤波最适用于平滑哪种噪声;3. 比较不同的窗口形状和窗口尺寸对整幅图象以及局部图象平滑效果的影响;4. 比较中值滤波和低通掩模法在平滑效果以及轮廓模糊程度上的差异;5. 尝试多种窗口线形组合的中值滤波。四. 实验程序清单主程序clear allI=(imread('cameraman.tif'); J1=imnoise(I,'salt & pepper',0.02);J2=imnoise
6、(I,'gaussian',0.02); J3=raylrnd(im2double(I);domain3ju=1 1 1%domain 需为 n*n 的窗口,形状可自己定义1 1 11 1 1;K3j1=ord(J1,5,domain3ju); %ord函数中需传递窗口,及窗口采中值num的位置, 此处为 5,K3j2=ord(J2,5,domain3ju); % 矩 形 窗 num=(n*n+1)/2, 十 字 窗 num=4*(n-1)/2)*(n-1)/2)+nK3j3=ord(J3,5,domain3ju); %字型窗 num=2*(n-1)+(n+1)/2,其余自己推
7、算 K301=LP(J1);x31=uint8(real(K301);K3=x31K401=LP(J2); x32=uint8(real( K401); K4=x32K5=LP(J3);figure(1) subplot(221),imshow(I) title('原图像') subplot(222),imshow(J1) title('加椒盐噪声') subplot(223),imshow(J2) title('加高斯噪声') subplot(224),imshow(J3) title('加瑞利噪声')figure(2)subp
8、lot(321),imshow(uint8( K3j1)title('加椒盐后滤波')subplot(322),imshow(uint8( K3j2)title('加高斯后滤波')subplot(323),imshow( K3j3)title('加瑞利后滤波')subplot(324),imshow( K3)title('加椒盐butterworth低通滤波')subplot(325),imshow( K4)title('加高斯butterworth低通滤波')subplot(326),imshow( K5)tit
9、le('加瑞利butterworth低通滤波')domain3shi=0 1 01 1 10 1 0;K3s1=ord(J1,7,domain3shi);K3s2=ord(J2,7,domain3shi);K3s3=ord(J3,7,domain3shi);figure(3)subplot(321),imshow(uint8( K3s1) title(' 加椒盐后滤波 ') subplot(322),imshow(uint8( K3s2) title(' 加高斯后滤波 ') subplot(323),imshow( K3s3) title(
10、9; 加瑞利后滤波 ') subplot(324),imshow( K3)title('加椒盐butterworth低通滤波')subplot(325),imshow( K4)title('加高斯butterworth低通滤波') subplot(326),imshow( K5)title('加瑞利butterworth低通滤波')domain3yi=0 0 01 1 10 0 0;K3y1=ord(J1,8,domain3yi); K3y2=ord(J2,8,domain3yi);K3y3=ord(J3,8,domain3yi); fi
11、gure(4)subplot(321),imshow(uint8( K3y1) title('加椒盐后滤波') subplot(322),imshow(uint8( K3y2) title('加高斯后滤波') subplot(323),imshow( K3y3)title('加瑞利后滤波') subplot(324),imshow( K3)title('加椒盐butterworth低通滤波') subplot(325),imshow( K4)title('加高斯butterworth低通滤波') subplot(3
12、26),imshow( K5)title('加瑞利butterworth低通滤波')domain5ju=1111 11 11111 11111 11111 1111 K5j1=ord(J1,13,domain5ju);K5j2=ord(J2,13,domain5ju);K5j3=ord(J3,13,domain5ju);figure(5)subplot(321),imshow(uint8( K5j1) title(' 加椒盐后滤波 ') subplot(322),imshow(uint8( K5j2) title(' 加高斯后滤波 ') subp
13、lot(323),imshow( K5j3) title(' 加瑞利后滤波 ') subplot(324),imshow( K3) title('加椒盐butterworth低通滤波') subplot(325),imshow( K4)title('加高斯butterworth低通滤波') subplot(326),imshow( K5)title('加瑞利butterworth低通滤波')domain5shi=0 0 1 0 0 0 0 1 0 0 1 1 1 1 1 0 0 1 0 0 0 0 1 0 0 ;K5s1=ord(
14、J1,21,domain5shi); K5s2=ord(J2,21,domain5shi); K5s3=ord(J3,21,domain5shi); figure(6) subplot(321),imshow(uint8( K5s1)title('加椒盐后滤波') subplot(322),imshow(uint8( K5s2) title('加高斯后滤波') subplot(323),imshow( K5s3) title('加瑞利后滤波') subplot(324),imshow( K3)title('加椒盐butterworth低通
15、滤波') subplot(325),imshow( K4)title('加高斯butterworth低通滤波') subplot(326),imshow( K5)title('加瑞利butterworth低通滤波')domain5yi=0 0 0 0 00 0 0 0 01 1 1 1 10 0 0 0 00 0 0 0 0 ;K5y1=ord(J1,23,domain5yi);K5y2=ord(J2,23,domain5yi);K5y3=ord(J3,23,domain5yi);figure(7)subplot(321),imshow(uint8( K
16、5y1)title(' 加椒盐后滤波 ') subplot(322),imshow(uint8( K5y2) title(' 加高斯后滤波 ') subplot(323),imshow( K5y3) title(' 加瑞利后滤波 ') subplot(324),imshow( K3)title('加椒盐butterworth低通滤波')subplot(325),imshow( K4)title('加高斯butterworth低通滤波') subplot(326),imshow( K5)title('加瑞利b
17、utterworth低通滤波')中值滤波function w=ord(J,num,domain) c1,k1=size(J);n,k2=size(domain);A=zeros(c1+n-1,k1+n-1);for i=1:c1%补零for j=1:k1A(i+(n-1)/2,j+(n-1)/2)=J(i,j);endendm=(n+1)/2;m0=(n-1)/2;B=zeros(n,n);for k=m:(c1+m-1)%窗口中心移动位置for l=m:(k1+m-1)for p1=-m0:m0%根据窗口产矩阵for q1=-m0:m0B(p1+m,q1+m)=A(k+p1,l+q1
18、); end endC=B.*domain;n2=n*n;D=reshape(C,1,n2) ;for x=1:n2-1%冒泡排序for y=1:n2-xif D(y)>D(y+1) temp=D(y); D(y)=D(y+1);D(y+1)=temp;endendendw(k-(n-1)/2,l-(n-1)/2)=D(num);endend低通滤波function x2=LP(I2)f=double(I2);g=fft2(f);g=fftshift(g);N1,N2=size(g);n=2;d0=50;n1=fix(N1/2);n2=fix(N2/2);for i=1:N1for j=
19、2:N2d=sqrt(i-n1F2+(j-n2)T);h=1心+0.414*(d/d0F(2* n);s1(i,j)=h*g(i,j);if(g(i,j)>50);s2(i,j)=0;elses2(i,j)=g(i,j);endendends1=ifftshift(s1);s2=ifftshift(s2);x2=ifft2(s1);五. 实验结果力口高Mh槳声力仃IK A'J d早产原图及加噪图像力口能挞后酒満力訂离E斤行淀漩力C 4rfl 就' butterworth-fH; 19'冏Jju 离 Hfr t»Li it*? r wo rt fi fkJ; i®. ?|f 匚盛丿MI 土福禾U bu tie rwortti fl£ i® 订E -披3*3矩形窗口力口嵩晰右滤漩<jn橄盐匚滤遮力n枫 盐
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年专业管材交易协议简明样本
- 环境保护与企业可持续发展紧密相连的关系考核试卷
- 六年级数学下册 第3单元达标测试卷(苏教版)
- 《基于诺丁斯关怀理论的临床护士人文关怀能力评价指标体系的构建与初步应用》
- 《中学生睡眠质量的影响因素及教育对策研究》
- 个人值班脱岗检讨书
- 生物质燃气的流动床气化技术研究考核试卷
- 天然气开采业的供应链风险与供应保障考核试卷
- 2024年楼梯间设计与施工一体化协议
- 淀粉黏附剂在造纸工业的应用考核试卷
- 工作岗位风险评估工作规程
- 智能存包柜(储物柜)产品技术说明书
- 電鍍技術資料大全
- 控轧控冷技术在钢材生产中的应用
- 连接器成本分析-B版
- 金融业不良资产处置法律服务方案书
- 基底节区解剖位置关系.ppt
- 电子商务购物平台的设计与开发论文
- xx大学成人高等教育校外教学点检查自查报告
- 《在政府教育工作督导评估反馈会上的表态发言》
- 六年级上册数学比的计算题
评论
0/150
提交评论