数字图像第一次作业.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）
升级会员 