www.gusucode.com > 关于海航matlab和lingo的训练题 > 图像处理.m

    clear all
g0=imread('z1.bmp');
%g0 = g0(:,:,2);  %三维转二维

figure(1);imshow(g0) ;title('原图')     %如图 5
g1=imnoise(g0,'salt & pepper',0.2)
g1=im2double(g1);
figure(2);imshow(g1);title('加入椒盐噪声')    %如图 6

h1=fspecial('gaussian',4,0.3)
g2=filter2(h1,g1,'same')
figure(3);imshow(g2);title('进行高斯滤波')    %如图 7

h2=fspecial('sobel')
g3=filter2(h2,g1,'same')
figure(4);imshow(g3);title('进行sobel滤波')    %如图   8

h3=fspecial('prewitt')
g4=filter2(h3,g1,'same')
figure(5);imshow(g4);title('进行prewitt滤波')    %如图 9

h4=fspecial('laplacian',0.5);
g5=filter2(h4,g1,'same');
figure(6);imshow(g5);title('进行拉普拉斯滤波');    %如图 10

h5=fspecial('log',4,0.3);
g6=filter2(h5,g1,'same');figure(7);
imshow(g6);title('进行高斯拉普拉斯滤波');    %如图 11

h6=fspecial('average');
g7=filter2(h6,g1,'same');figure(8);
imshow(g7);title('进行均值滤波');    %如图 12

h7=fspecial('unsharp',0.3);
g8=filter2(h7,g1,'same');
figure(9);imshow(g8);title('进行模糊滤波');    %如图 13

h8=[0 -1 0;-1 5 -1;0 -1 0];
g9=filter2(h8,g1,'same');
figure(10);imshow(g9);title('进行高通高斯滤波');    %如图 14

h9=g1;g10=medfilt2(h9);
figure(11);imshow(g10);title('进行中值滤波');    %如图 15