1、第 三 章 基 于 MATLABD 的 模 拟 噪 声 生 成MATLAB 图 像 处 理 工 具 箱 提 供 的 噪 声 添 加 函 数 imnoise,它 可 以 对 图 添加 一 些 典 型 的 噪 声 。 其 语 法 是 ：J=imnoise(I,type)J=imnoise(I,type,parameters)其 功 能 是 ： 返 回 对 原 图 像 I添 加 典 型 噪 声 的 图 像 J,参 数 type和parameters用 于 确 定 噪 声 的 类 型 和 相 应 的 参 数 。三 种 典 型 的 噪 声 ：type=gaussian时 ,为 高 斯 噪 声 ；type
2、=salt %读 取 图 像I=rgb2gray(i);I1=imnoise(I,gaussian); %加 高 斯 噪 声I2=imnoise(I,salt %加 椒 盐 噪 声I3=imnoise(I,speckle); %加 乘 性 噪 声subplot(221),imshow(I); %显 示 图 像 Ititle(原 图 像 );subplot(222),imshow(I1);title(受 高 斯 噪 声 污 染 的 图 像 );subplot(223),imshow(I2);title(受 椒 盐 噪 声 污 染 的 图 像 );subplot(224),imshow(I3);t
3、itle(受 乘 性 噪 声 污 染 的 图 像 );处 理 结 果 如 图 3.1 所 示 ：图 3.1 图像添加噪声处理结果第 四 章 均 值 滤 波 处 理 方 法4.2 均 值 滤 波 法 对 图 像 的 处 理在 MATLAB 图 像 处 理 工 具 箱 中 ， 提 供 了 imfilter 函 数 用 于 实 现 均 值 滤波 ， imfilter 的 语 法 格 式 为 ：B=imfilter(A,H)其 功 能 是 ， 用 H 模 板 对 图 像 A 进 行 均 值 滤 波 ， 模 板 为 ：取 H1,程 序 如 下 ：i=imread(king.jpg); %读 取 图 像1
4、91 11251HI=rgb2gray(i);I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=imnoise(I,speckle);H1=ones(3,3)/9; %33领 域 模 板J=imfilter(I,H1); %领 域 平 均J1=imfilter(I1,H1);J2=imfilter(I2,H1);J3=imfilter(I3,H1);subplot(221),imshow(J);title(原 图 像 滤 波 后 );subplot(222),imshow(J1);title(高 斯 污 染 图 像 滤 波 后 );subplot(223)
5、,imshow(J2);title(椒 盐 污 染 图 像 滤 波 后 );subplot(224),imshow(J3);title(乘 法 污 染 图 像 滤 波 后 );运 行 结 果 如 图 4.3图 4.3 33领 域 模 板 处 理 结 果取 H2， 程 序 如 下 ：i=imread(king.jpg); %读 取 图 像I=rgb2gray(i);I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=imnoise(I,speckle);H2= ones(5,5)/25; %55领 域 模 板J=imfilter(I,H2); %领 域 平 均
6、J1=imfilter(I1,H2);J2=imfilter(I2,H2);J3=imfilter(I3,H2);subplot(221),imshow(J);title(原 图 像 滤 波 后 );subplot(222),imshow(J1);title(高 斯 污 染 图 像 滤 波 后 );subplot(223),imshow(J2);title(椒 盐 污 染 图 像 滤 波 后 );subplot(224),imshow(J3);title(乘 法 污 染 图 像 滤 波 后 ); 运 行 结 果 如 图 4.4：图 4.4 55领 域 模 板 处 理 结 果第 五 章 中 值
7、滤 波 处 理 方 法5.2 中 值 滤 波 法 对 图 像 的 处 理MATLAB图 像 处 理 工 具 箱 提 供 了 medfilt2函 数 用 于 中 值 滤 波 。 其 语 法 格式 为 ：B=medfilt2(A) 其 功 能 为 ： 用 33的 滤 波 窗 口 对 图 像 A进 行 中值 滤 波 ；B=medfilt2(A,m,n) 其 功 能 是 ： 用 大 小 为 mn的 窗 口 对 图 像 A进 行中 值 滤 波 ；B=medfilt2(A,indexed,.) 其 功 能 为 ： 对 索 引 图 像 A进 行 中 值 滤波 ；用 33中 值 滤 波 模 板 和 55中 值
8、 滤 波 模 板 对 图 像 进 行 处 理 ， 程 序 如 下 ：i=imread(king.jpg); I=rgb2gray(i);I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=imnoise(I,speckle);J1=medfilt2(I1,3,3); %33中 值 滤 波 模 板J2=medfilt2(I2,3,3);J3=medfilt2(I3,3,3);J4=medfilt2(I1,5,5); %55中 值 滤 波 模 板J5=medfilt2(I2,5,5);J6=medfilt2(I3,5,5);figure,subplot(121)
9、,imshow(J1);title(高 斯 3*3中 值 滤 波 后 );subplot(122),imshow(J2);title(椒 盐 3*3中 值 滤 波 后 );figure,subplot(121),imshow(J3);title(乘 法 3*3中 值 滤 波 后 );subplot(122),imshow(J4);title(高 斯 5*5中 值 滤 波 后 );figure,subplot(121),imshow(J5);title(椒 盐 5*5中 值 滤 波 后 );subplot(122),imshow(J6);title(乘 法 5*5中 值 滤 波 后 );处 理
10、结 果 如 图 5.1图 5.1 33中 值 滤 波 模 板 和 55中 值 滤 波 模 板 对 图 像 进 行 处 理第 六 章 频 域 低 通 滤 波 法用 理 想 低 通 滤 波 器 对 图 像 处 理 ， 程 序 如 下 （ 以 受 高 斯 噪 声 污 染 图 像 为例 ） ：i=imread(king.jpg);I=rgb2gray(i);I1=imnoise(I,gaussian);s=fftshift(fft2(I1);a,b=size(s);a0=round(a/2);b0=round(b/2);d0=10;%这 里 d取 10和 50作 比 较for i=1:a for j=
11、1:bd=sqrt(i-a0)2+(j-b0)2);if d =d1,h=0;else h=1-(d-d0)/(d1-d0);end;end;s(i,j)=h*s(i,j);end;end;s=uint8(real(ifft2(ifftshift(s);subplot(121),imshow(I1);title(受 高 斯 噪 声 污 染 的 图 像 );subplot(122),imshow(s);title(高 斯 污 染 图 像 处 理 后 );处 理 结 果 如 图 6.9所 示图 6.9 梯 形 低 通 滤 波 器 对 图 像 的 处 理 结 果6.5 构 建 二 维 滤 波 器 对
12、 图 像 的 处 理首 先 构 建 二 维 滤 波 器 d， 传 递 函 数 如 图 6.10所 示 ;f1,f2=freqspace(25,meshgrid);Hd=zeros(25,25);d=sqrt(f1.2+f2.2)0.5; %0.5为 截 止 半 径 大 小Hd(d)=1; h=fsamp2(Hd);figure,freqz2(h,64,64);然 后 用 所 构 建 的 二 维 滤 波 器 对 以 上 图 像 进 行 滤 波i=imread(king.jpg); I=rgb2gray(i);I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=
13、imnoise(I,speckle);J=imfilter(I,h,replicate);J1=imfilter(I1,h,replicate);J2=imfilter(I2,h,replicate);J3=imfilter(I3,h,replicate);subplot(221),imshow(J);title(原 图 像 滤 波 后 );subplot(222),imshow(J1);title(高 斯 污 染 图 像 滤 波 后 );subplot(223),imshow(J2);title(椒 盐 污 染 图 像 污 染 后 );subplot(224),imshow(J3);title(乘 法 污 染 图 像 滤 波 后 );其 运 行 结 果 如 下 图 6.11所 示 :图 6.10 用 频 率 采 样 法 构 建 的 二 维 滤 波 器参 考 文 献1 曹 茂 永 ， 数 字 图 像 处 理 .北 京 大 学 出 版 社 ， 20072 唐 向 宏 ， MATLAB 及 在 电 子 信 息 类 课 程 中 的 应 用 .电 子 工 业 出 版 社 ，20103 杨帆，等数字图像处理与分析北京航空航天大学出版社，2007.4http:/ ， 百 度 文 库
