数字图像第一次作业.docx
《数字图像第一次作业.docx》由会员分享,可在线阅读,更多相关《数字图像第一次作业.docx(22页珍藏版)》请在冰豆网上搜索。
数字图像第一次作业
数字图像处理上机
(一)
1.将一幅彩色图像转换成灰度图像
>s=imread('a2.jpg');
>>A=rgb2gray(s);
>>imshow(A)
>>imshow(s)
原图
Gray图
2.将一幅灰度图像转换成二值图像
image=imread('a1.jpg');
tt=graythresh(image);
BW=im2bw(image,tt);
imshow(BW)
3.对一幅图像求反
>>s=imread('a2.jpg');
A=rgb2gray(s);
I=double(A)
O=255-I;
O=uint8(O);
imshow(O)
4.对一幅图像进行线性灰度变换
5.s=imread('a2.jpg');
6.A=rgb2gray(s);
7.I=double(A);
8.[m,n]=size(I);
9.
10.fori=1:
m
11.forj=1:
n
12.ifI(i,j)<=30
13.I(i,j)=I(i,j);
14.elseifI(i,j)<=150
15.I(i,j)=(200-30)/(150-30)*(I(i,j)-30)+30;
16.else
17.I(i,j)=(255-200)/(255-150)*(I(i,j)-150)+200;
18.end
19.end
20.end
21.end
22.O=uint8(I);
23.imshow(O)
24.
5.将一幅灰度图像进行傅里叶变换,并求其傅里叶频谱、对数傅里叶频谱以及
其傅里叶频谱
s=imread('a2.jpg');
A=rgb2gray(s);
I=double(A);
[m,n]=size(I);
O=fft2(I);
imshow(abs(O),[]);
其傅里叶中心普
s=imread('a2.jpg');
A=rgb2gray(s);
I=double(A);
[m,n]=size(I);
O=fftshift(fft2(I));
imshow(abs(O),[]);
对数傅里叶频谱
s=imread('a2.jpg');
A=rgb2gray(s);
I=double(A);
[m,n]=size(I);
O=fftshift(fft2(I));
imshow(log(abs(O)),[]);
6.将一幅灰度图像进行离散余弦变换并求其频谱,最后请说明离散余弦变换的重要应用
s=imread('a2.jpg');
I=rgb2gray(s);
O=dct2(I);
imshow(log(abs(O)),[]);
可用于图像压缩。
7.将一幅图像分别进行直方图均衡化和规定化操作
s=imread('a2.jpg');
I=rgb2gray(s);
figure
imshow(I)
figure
imhist(I);
II=histeq(I);
figure
imhist(II);
figure
imshow(II)
s=imread('a2.jpg');
I=rgb2gray(s);
[m,n]=size(I);
fori=1:
8:
257
count(i)=i;
end
N=histeq(I,count);
figure
imshow(N)
figure
imhist(N);
8.选择两幅图像分别进行“加”、“减”、“乘”和“除”操作
c1=imread('c1.png');
c2=imread('c2.png');
figure
imshow(c1)
figure
imshow(c2)
d1=c1+c2;
figure
imshow(d1)
d2=c1-c2;
figure
imshow(d2)
d3=c1.*c2;
figure
imshow(d3)
d4=c1./c2;
figure
imshow(d4)
选择两幅图像分别进行“与”、“或”和“非”操作
c1=imread('c1.png');
c2=imread('c2.png');
figure
imshow(c1)
figure
imshow(c2)
d1=bitor(c1,c2);
figure
imshow(d1)
d2=bitand(c1,c2);
figure
imshow(d2)
c2=rgb2gray(c2);
[m,n]=size(c2);
p=ones(m,n)*255;
fori=1:
m
forj=1:
n
d(i,j)=bitxor(p(i,j),c2(i,j));
end
end
figure
imshow(d)