数字图像处理论文合编.doc

滨江学院 课程论文（设计） 课程名称 《数字图像处理》 学 期 2014-2015（1） 院 系 电子工程系 专 业 通信工程 班 级 12通信3班 学生姓名 范勤考 学 号 20122334904 二Ｏ一四 年 十二 月 二十 日 考核内容 1.在Matlab下读入一幅图像，对其灰度图作快速傅立叶变换； 2.在Matlab下读入一幅图像，对其灰度图作DCT变换； 3.在Matlab下读入一幅图像，对其作DWT分解。 4.读入一幅有明显明暗缺陷的灰度图像，分析其直方图特征，分别用分段灰度线性变换法、直方图均衡法和图像灰度调整法(imadjust)分别进行处理； 5.读入一幅灰度图像 ，分别给其加上“乘性噪声”、“椒盐噪声”，然后分别用“均值滤波”、“中值滤波”和“巴特沃斯滤波”对其做平滑处理； 6.读入一幅灰度图像，分别用“Sobel算子”、“Laplacian算子”、“梯形滤波器”对其做锐化处理； 7.读入一幅灰度图像，分别用“灰度级分层法”、“灰度变换法”和“频域伪彩色处理法”对其进行伪彩色增强； 8.读入一幅图像，对其进行模糊化，然后用“逆滤波法”对其进行复原处理。 要 求 1.独立完成各项内容； 2.记录每一内容实现的步骤； 3.编写FFT、DCT、DWT变换Matlab程序； 4.打印输出原始图像、FFT的频谱图、将频率平面坐标原点移至窗口中心的FFT频谱图、DCT频谱图以及1级小波分解图像； 5.编写“灰度增强”、“平滑”、“锐化”、“伪彩色增强”、“图像复原”的Matlab程序； 6.打印输出经“灰度增强”、“平滑”、“锐化”、“伪彩色增强”、“图像复原”处理的图像及各自的原始图像； 7.独立撰写课程论文和设计。 1. 1.读入一幅图像，对其灰度图作快速傅立叶变换。 用imread()函数读如工作目录下的“风光壁纸33.jpg”图像存于I矩阵中，用rgb2gray()函数将其转换成灰度图像； 用fft2()函数对其进行FFT变换，并将变换得到的傅立叶频谱存于fft_I矩阵中； 用fftshift()函数将傅立叶频谱坐标原点移至窗口中央并存于sfft_I矩阵中； 显示FFT频谱图、移动后的频谱图。 由于fft_I和sfft_I均为复数矩阵，须用abs()函数求其模，并根据矩阵元素的特点将其归一化到0~255之间； 用subplot()函数及imshow函数在同一窗口下显示原始图像、灰度图像、FFT频谱图及移动后的频谱图。 附：FLY.m文件 I=imread(C:\Users\BJXY\pics\DIPUM2E_International_Original_Book_Images\DIPUM2E_International_Original_Book_Images\DIPUM2E_International_Ed_CH06_Images\Fig0622(a).tif) subplot(221) imshow(I); title(原始图像) I=rgb2gray(I); subplot(222) imshow(I); title(灰度图像) fft_I=fft2(I); % 2-D快速傅立叶变换 A=abs(fft_I); %将频谱矩阵元素归一化到0~255 A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255; subplot(223) imshow(A); title(傅立叶频谱图像) sfft_I=fftshift(fft_I); % 傅立叶频谱平面中心移至窗口中心 A=abs(sfft_I); %将频谱矩阵元素归一化到0~255 A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255; subplot(224) imshow(A); title(原点移到中心的傅立叶频谱图像) 2I=imread(C:\Users\BJXY\pics\DIPUM2E_International_Original_Book_Images\DIPUM2E_International_Original_Book_Images\DIPUM2E_International_Ed_CH06_Images\Fig0621(a).tif) subplot(221) imshow(I); title(原始