《数字图像处理课程论文》.docVIP

滨江学院 课程论文（设计）课程名称 《数字图像处理》学 期 2014-2015（1） 院 系 电子工程系专 业 通信工程班 级 12通信3班学生姓名范勤考学 号 904二Ｏ一四 年 十二 月 二十 日考核内容 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~255A=(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~255A=(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(原始图像)I=rgb2gray(I);subplot(223)imshow(I);title(原始灰度图像)dct_I=dct2(I);A=abs(dct_I);A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255;subplot(224)imshow(A);title(离散余弦频谱图像)3. clear clc X=imre