资源描述:
《数字图像处理.docx》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、内蒙古科技大学本科生课程设计论文数字图像处理课设学生姓名:于鹏飞学号:专业:电子信息工程班级:13信息二班指导教师:张明2016年5月实验一、计算图象统计参数及直方图变换编写程序:1.读取图像(文件名为‘cameraman.tif’);2.最大值、最小值、均值3.直方图(要求不能调用imhist函数,只能用该函数来验证自编直方图程序块的正确性)。4.计算并分析图像直方图(imhist)5.利用直接灰度变换法对图像进行灰度变换(imadjust)6.直方图均衡化处理及显示(histeq)程序f=imread('C:UsersYSSADesktop数字图像课设于鹏飞图片came
2、raman.tif');%求最大值、最小值、均值max(f(:))%最大值min(f(:))%最小值mean(f(:))%均值%显示其直方图[m,n]=size(f);%图像大小X=zeros(256,1);%灰度级范围fori=1:m%对应的统计位加1forj=1:nx=f(i,j);X(x)=X(x)+1;endend;%计算并分析图像直方图(imhist)g=imhist(f,256);%利用直接灰度变换法对图像进行灰度变换(imadjust)g2=imadjust(f,[01],[10]);%对原图像进行直方图均衡化处理g3=histeq(f);subplot(2,3,1),im
3、show(f);title('原始图像');subplot(2,3,2);bar(X);title('原图像直方图')subplot(2,3,3);bar(g);title('验证原图像直方图')subplot(2,3,4),imshow(g2);title('灰度变换');subplot(2,3,6),imshow(g3);title('直方图均衡化');照片实验二傅里叶变换及空间域和频率域滤波1、利用MATLAB实现数字图像的傅里叶变换读入并显示图“test1.tif”,作该图的二维FFT变换F,将其直流分量移到频谱中心F1(fftshift()函数),计算其实部RR、虚部II,用两
4、种方法计算幅值A1=abs(F1)和A2=sqrt(RR.^2+II.^2),分别显示A1和A2,并加以比较。2、近似冲击函数二维傅里叶变换A=zeros(99,99);A(49:51,49:51)=1;作A的二维傅里叶变换B,将B直流分量移到频谱中心B1,分别用函数imshow和mesh显示A和B1模的对数(log(1+abs(B1)))。程序g=imread('C:UsersYSSADesktop数字图像课设于鹏飞图片test1.tif');%对图像进行傅立叶变换G=fft2(g);%将其直流分量移到频谱中心F=fftshift(G);%实部RR、虚部IIRR=real
5、(F);II=imag(F);%显示模A1和A2,并加以比较A1=abs(F);A2=sqrt(RR.^2+II^2);aa1=log(1+A1);aa2=log(1+A2);%生成一个99*99的零矩阵A=zeros(99,99);A(49:51,49:51)=1;%B=fft2(A);%fftshift就是对换数据的左右两边B1=fftshift(B);%显示A和B1模的对数(log(1+abs(B1)))figure;subplot(2,4,1),imshow(g);title('原图像');subplot(2,4,2),imshow(G);title('傅立叶变换');subpl
6、ot(2,4,3),imshow(aa1,[]);title('a1abs');subplot(2,4,4),imshow(aa2,[]);title('a2sqrt');subplot(2,4,5),imshow(A);title('冲击函数imsow');subplot(2,4,6),mesh(A);title('冲击函数mesh');subplot(2,4,7),imshow(log(1+abs(B1)));title('冲击函数频域模mesh');subplot(2,4,8),mesh(log(1+abs(B1)));title('冲击函数频域模mesh');实验2_1-2图片3
7、、空间滤波与频域滤波将图“test2.tif”分别进行空间与频域滤波。空间滤波:用fspecial产生9×9,标准差为2的高斯滤波器w,用函数imfilter对f进行空间滤波,得到滤波后的图像fi1。频域滤波:将上述高斯滤波器w用函数freqz2获得其256×256频域形式W,用W在频域滤波图像f(注意W的直流分量在频谱中心),得到Ff,求其傅里叶逆变换(ifft2),得到滤波后的图像fi2。用函数imshow显示w,用函数mesh