数字图像处理课程论文

滨江学院

课程论文（设计）课程名称《数字图像处理》

学期 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_Origina l_Book_Images\DIPUM2E_International_Original_Book_Images\D IPUM2E_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_Or iginal_Book_Images\DIPUM2E_International_Original_Book_Ima ges\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=imread('C:\Users\BJXY\pics\DIPUM2E_International_Origina