数字图像处理实验结果.docx
- 文档编号:23367739
- 上传时间:2023-05-16
- 格式:DOCX
- 页数:17
- 大小:90.13KB
数字图像处理实验结果.docx
《数字图像处理实验结果.docx》由会员分享,可在线阅读,更多相关《数字图像处理实验结果.docx(17页珍藏版)》请在冰豆网上搜索。
数字图像处理实验结果
注意:
打印后实验三第一个代码要解释代码
试验一
1.1灰度图像到索引图像转换
RGB=imread('2.jpg');
figure(3);
imshow(RGB);
figure
(1);
[RGB1,map1]=gray2ind(RGB,128);
imshow(RGB1,map1);
figure
(2);
[RGB2,map2]=gray2ind(RGB,16);
imshow(RGB2,map2);
imwrite(RGB1,map1,'3.bmp');
imwrite(RGB2,map2,'4.bmp');
info=imfinfo('3.bmp')
1.1.2索引到灰度转换
参考代码
[X,map]=imread('3.bmp');
figure
(1);
imshow(X);
figure
(2);
I=ind2gray(X,map);
imshow(I);
imwrite(I,'new.jpg');
info=imfinfo('new.jpg')
1.2RGB图像到索引图像的转换
RGB=imread('1.jpg');
figure
(1);
imshow(RGB);
figure
(2);
y=rgb2ind(RGB,256);
imwrite(y,'1s.bmp');
imshow(y);
info=imfinfo('1s.bmp')
1.2.2索引转换RGB
参考代码及结果
[X,MAP]=imread('1s.bmp');
figure
(1);
imshow(X);
figure
(2);
RGB=ind2rgb(X,MAP)
imwrite(RGB,'RGB.jpg');
imshow(RGB);
info=imfinfo('RGB.jpg')
1.3.1RGB图像到二值图像的转换
RGB=imread('1.jpg');%将图像“1.jpg”的数据赋给变量RGB
figure(3);%设定显示窗口3
imshow(RGB);%显示原“1.jpg”文件
bw1=im2bw(RGB,0.5);%阀值为0.5进行二值转换
bw2=im2bw(RGB,0.6);%阀值为0.6进行二值转换
figure
(1);%设定显示窗口1
imshow(bw1);%显示索引图像bw1
figure
(2);%设定显示窗口2
imshow(bw2);%显示索引图像bw2
1.3.2选取灰度图像进行二值转换
参考代码如下(阈值要修改)
RGB=imread('3.bmp');
figure(3);
imshow(RGB);
bw1=im2bw(RGB,0.25);bw2=im2bw(RGB,0.35);
figure
(1);
imshow(bw1);
figure
(2);
imshow(bw2);
2.1卷积实验
A=magic(3);
B=ones(3);
A(8,8)=0;
B(8,8)=0;
C=ifft2(fft2(A).*fft2(B));
C=C(1:
5,1:
5);
C=real(C)
C=
8.00009.000015.00007.00006.0000
11.000017.000030.000019.000013.0000
15.000030.000045.000030.000015.0000
7.000021.000030.000023.00009.0000
4.000013.000015.000011.00002.0000
参考代码
C=conv2(A,B);
C=C(1:
5,1:
5)
C=
891576
1117301913
1530453015
72130239
41315112
2.2对一副灰度图片进行二维傅立叶变换
参考代码
figure
(1);
RGB=imread('2.jpg');
imshow(RGB);
figure
(2);
B=fftshift(fft2(RGB));
imshow(log(abs(B)),[]),colormap(jet(64)),colorbar;
2.3创建大小为128*128,背景为黑色,中间开出一个8×8白色的窗口图像
参考代码
clear
x=zeros(128,128);
x(60:
68,60:
68)=1;
imshow(x)
imwrite(x,'new.bmp')
figure
(1);
RGB=imread('new.bmp');
imshow(RGB);
figure
(2);
B=fftshift(ifft2(RGB));
imshow(log(abs(B)),[]),colormap(jet(64)),colorbar;
3.1进行时域灰度变换
X1=imread('image.tif');
figure,imshow(X1)
f0=0;g0=0;
f1=70;g1=30;
f2=180;g2=230;
f3=255;g3=255;
r1=(g1-g0)/(f1-f0);
b1=g0-r1*f0;
r2=(g2-g1)/(f2-f1);
b2=g1-r2*f1;
r3=(g3-g2)/(f3-f2);
b3=g2-r3*f2;
[m,n]=size(X1);
X2=double(X1);
fori=1:
m
forj=1:
n
f=X2(i,j);
g(i,j)=0;
if(f>=0)&(f<=f1)
g(i,j)=r1*f+b1;
elseif(f>=f1)&(f<=f2)
g(i,j)=r2*f+b2;
elseif(f>=f2)&(f<=f3)
g(i,j)=r3*f+b3;
end
end
end
figure,imshow(mat2gray(g))
3.2实现Butterworth低通滤波器的参考代码
I1=imread('3.bmp');
figure,imshow(I1)
I2=imnoise(I1,'salt');
figure,imshow(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);
fori=1:
N1
forj=1:
N2
d=sqrt((i-n1)^2+(j-n2)^2);
h=1/(1+(d/d0)^(2*n));
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
X2=ifft2(result);
X3=uint8(real(X2));
figure,imshow(X3)
4.1均值滤波
3*3滤波代码
I=imread('2.jpg');
I=imnoise(I,'gaussian',0,0.02);%添加高斯噪声
figure
(1)
imshow(I);
h=[111
111
111];%产生滤波模板,不同模板对应的滤波形式也不同
h=h/25;%滤波模板归一化
J=conv2(I,h);%滤波并显示图像
figure
(2);
imshow(J,[])
5*5,滤波代码
I=imread('2.jpg');
I=imnoise(I,'gaussian',0,0.02);%添加高斯噪声
figure
(1)
imshow(I);
h=[11111
11111
11111
11111
11111];%产生滤波模板,不同模板对应的滤波形式也不同
h=h/25;%滤波模板归一化
J=conv2(I,h);%滤波并显示图像
figure
(2);
imshow(J,[])
7*7滤波参考代码
I=imread('2.jpg');
I=imnoise(I,'gaussian',0,0.02);%添加高斯噪声
figure
(1)
imshow(I);
h=[1111111
1111111
1111111
1111111
1111111
1111111
1111111];%产生滤波模板,不同模板对应的滤波形式也不同
h=h/49;%滤波模板归一化
J=conv2(I,h);%滤波并显示图像
figure
(2);
imshow(J,[])
4.2中值滤波
参考代码
I=imread('2.jpg');%将图像的数据赋给变量I
J=imnoise(I,'gaussian',0.02);%加入噪声
K=medfilt2(J);%进行中值滤波
subplot(211),imshow(J);%分别显示原图像和变换后的图像
subplot(212),imshow(K)
4.3加椒盐噪声后中值滤波与均值滤波对比
参考代码
I=imread('b.jpg');
N=rgb2gray(I);
subplot(2,2,1),imshow(N);title('灰度图')
J=imnoise(J,'salt&pepper',0.005);
K=medfilt2(J);
subplot(2,2,2),imshow(J);title('加椒盐噪声图')
subplot(2,2,3),imshow(K);title('中值滤波图')
h=[111
111
111];
h=h/9;
M=conv2(N,h);
subplot(2,2,4),imshow(M);title('均值滤波图')
处理图像结果
加高斯噪声后中值滤波与均值滤波对比
参考代码
I=imread('b.jpg');
N=rgb2gray(I);
subplot(2,2,1),imshow(N);title('灰度图')
J=imnoise(J,'gaussian',0,0.02);
K=medfilt2(J);
subplot(2,2,2),imshow(J);title('加高斯噪声图')
subplot(2,2,3),imshow(K);title('中值滤波图')
h=[11
11]
h=h/4;
M=conv2(N,h);
subplot(2,2,4),imshow(M);title('均值滤波图')
4.3边缘检测
I=imread('2.jpg');
subplot(411),imshow(I);title('原图')
BW1=edge(I,'Roberts');subplot(412),imshow(BW1);title('Roberts')
BW2=edge(I,'sobel');subplot(413),imshow(BW2);title('sobel')
BW3=edge(I,'Prewitt')
subplot(414),imshow(BW3);title('Prewitt')
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字图像 处理 实验 结果