数字图像作业详解.doc

滨江学院 课程论文（设计） 课程名称 《数字图像处理》 学 期 2014-2015（1） 院 系 电子工程系 专 业 通信工程 班 级 12通信（3）班 学生姓名 岳巍 学 号 20122334931 二Ｏ一四 年 十二 月 二十 日 考核内容 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.独立撰写课程论文和设计。 论文设计内容 在Matlab下读入一幅图像，对其灰度图作快速傅立叶变换； I=imread( C:\Users\Administrator\Desktop195852438137_2.jpg); imshow(I); 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(原点移到中心的傅立叶频谱图像) 在Matlab下读入一幅图像，对其灰度图作DCT变换； I=imread(E:\454545\4444.jpg); imshow(I); subplot(2,2,1) imshow(I); title(原始图像) I=rgb2gray(I); subplot(2,2,2) imshow(I); title(灰度图像) dct_I=dct2(I); A=abs(dct_I); A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255; subplot(223) imshow(A); title(离散余弦频谱图像) 在Matlab下读入一幅图像，对其作DWT分解 clear clc X= imread(E:\454545\4444.jpg); subplot(231) imshow(X) title(原始图像) X=rgb2gray(X); subplot(234) imshow(X) title(原始灰度图像) %对图象用2D小波进行一次分解 [c,s]=wavedec2(X,1,haar); %提取特征分量和细节分量 a1=wrcoef2(a,c,s,haar); h1=wrcoef2(h,c,s,haar); v1=wrcoef2(v,c,s,haa