数字图像采集与处理工程导论作业

《数字图像采集与处理工程导论作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、“数字图像采集与处理工程导论”实验报告姓名:***学号：0209****年级：2009级班级：0209**班学院：电子工程学院专业：信息对抗技术电话:用Matlab语言完成如下实验:(1)打开一幅BMP图像；(2)将局部区域的灰度值进行改变;(3)另存为一个新的BMP文件。实验过程:(1)选取用matlab打开:I=imread(9avril.bmp5)(1)使用命令1(150:300,0=300改变其部分区域值;用命令：image(I)杳看改变灰度后的结杲如25010015020025030035040045050100150200250300350

2、(1)使用命令imwrite(I,new.bmp,)将修改后的图像另存为名为new.bmp的文件，可见在当前工作口录下产生了名为new.bmp的图像文件。MATLAB编程实现傅里叶高通，低通滤波，给出算法原理及实验结果。程序:%%DOCUMENTTITLE%INTRODUCTORYTEXTfilename='avril.jpg1;Img=imread(filename);imshow(Img);%functionIdealHighFilter(Img,p,q,d)%%if〜isgray(Img)Img=rgb2gray(Img);%end;figure

3、(l);subplot(2,2,l);imshow(Img);title('原图像')；f=fftshift(fift2(Img));subplot(2,2,2);imshow(abs(f),[])；title('原图像的傅里叶频谱J;subplot(2,2,3);imshow(log(abs(f)),[]);title('原图像的傅里叶变换取对数后的频谱J;[ab]=size(f);a0=round(a/2);b0=round(b/2);fori=l:aforii=l:bd=10;p=0.2;q=0.5;dis=sqrt((i-aO)A2+(ii-b

4、O)A2);ifdis>=dh=l；elseh=0;end;s(i,ii)=(p+q*h)*f(i,ii);end;%end-SHUinsrea一(iff〔2(iff〔shift(s))))八Subpot(224)jmshow(s)2-o(^s5^^nn3^w)；end%funsionIdearOWFi一Or(lmg9d)%%if〜isgray(Img)%Imgurgb2gnly(Imgr%end-figure(2)八subpoi(22l)；imshow(Img)；ts:e(・WIH^_)；Irfftshift(fft2(lmg))；subpot(22

5、2)；imshow(abs(psule(^B^satffl#®i5;subpot(223)；imshow(log(abs(D)=)；titk(・^函^saJR三^^®xl^EA3Bi5-「ab〒size(D;aouround(巴2)八bouround(b、2)；forn-rraforiH-rrbdisHsqrt((7-ao)>2+(7:bo)>2)_ifdiscedhurre-sehHpend八s(Lii)uhwLii)；end八%end八s=uint8(real(ifft2(ifftshift(s))));subplot(2,2,4),imshow

6、(s),title(，高通滤波后的频谱');End实验结果：原图像原图像的傅里叶频谱原图像的傅里叶变换取对数后的频谱低通滤波后的频谱原图像的傅里叶变换取对数后的频谱高通滤波后的频谱找一个曝光不足的灰度(或彩色)图像，用matlab按照直方图均衡化的方法进行处理。程序：clc;clear;PS=rgb2gray(imread('3.jpg'));subplot(2,2,l)imshow(PS)title。曝光不足的图像)%绘制直方图[m,n]=size(PS);P=zeros(1,256);%预创建存放灰度IH现概率的向量fori=1:256P(i)=l

7、ength(find(PS==(i-1)))/(m*n);%计算灰度概率,存入GP相应位置endsubplot(2,2,2)bar(0:255,P,g)titleC原图像直方图)xlabelC灰度值)ylabel('出现概率')%直方图均衡化Sl=zeros(1,256);fori=1:256forj=l:iSl(i)=P(j)+Sl(i);endendS2=round(S1*255);fori=1:256Pl(i)=sum(P(find(S2==i)));%绘制直方图%计算Sk%将$1＜归到相近级的灰度%计算现冇每个灰度级出现的概率endsubplo

8、t(2,2,4)bar(0:255,Pl,TV)title(均衡化后的直方图Jxlabelf灰