内蒙古工业大学数字图像处理.docVIP

学校代码： 10128学 号： 201220203053 学校代码： 10128 学 号： 201220203053 数字图像处理课程设计 题 目： 数字图像处理及算术编码（或DCT压缩编码）仿真实现 学生姓名： 陈浩 学 院： 信息工程学院 系 别： 电子信息工程系 专 业： 电子信息工程 班 级： 电子12-2 指导教师： 吕方 黎玉玲 2015 年 12月 28 日 数字图像处理课程设计 1、课程设计目的 通过本课程设计使学生了解数字图像的基本概念，掌握数字图像处理的基本内容，如图像点运算、几何变换、增强处理、图像复原、边缘检测以及图像压缩等的基本原理和Matlab实现方法。 通过本次课程设计，让学生掌握如何学习一门语言，如何进行资料查阅搜集，如何自己解决问题等方法，养成良好的学习习惯。扩展理论知识，培养学生的综合设计能力。 2、课程设计内容 2.1 图像处理基本功能 1）数字图像的变换：普通傅里叶变换(ft)与逆变换（ift）、快速傅里叶变换（fft）与逆变换(ifft)、离散余弦变换（DCT），小波变换。 2) 数字图像直方图的统计及绘制等； 2.2 图像处理综合功能 1）图像平滑算法程序设计： 2）DCT压缩（保留不同系数），要求显示原图像、压缩后图像的文件大小、压缩比或算术编码压缩 3、 课程设计的一般步骤 选题与搜集资料：选择课题，进行系统调查，搜集资料。 分析与设计：根据搜集的资料，进行功能分析，并对系统功能与模块划分等设计。 程序设计：运用掌握的语言，编写程序，实现所设计的功能。 调试与测试：自行调试程序，同学之间交叉测试程序，并记录测试情况。 验收与评分：指导教师对每个成员开发的程序进行综合验收，结合设计报告，根据课程设计成绩的评定方法，评出成绩。 1.普通傅里叶变换(ft)与逆变换（ift） 傅里叶变换是一个非常复杂的理论，我们在图像处理中集中关注于其傅里叶离散变换离散傅立叶变换(Discrete Fourier Transform)。 图象是由灰度（RGB）组成的二维离散数据矩阵，则对它进行傅立叶变换是离散的傅立叶变换。 对图像数据f(x,y)（x=0,1,… ,Ｍ-1； y=0,1,… ,N-1）。则其离散傅立叶变换定义可表示为： 式中，u=0,1,…, M-1;v= 0,1,…, N-1 其逆变换为 式中，x=0,1,…, M-1;y= 0,1,…, N-1 程序如下： t=imread(D:\chenhao.jpg); t=rgb2gray(t); figure,subplot(2,1,1) imshow(t); title(原灰度图象); [m,n]=size(t); M=zeros(m,m); N=zeros(n,n); for i=1:m for k=1:m M(i,k)=exp(2*pi*i*k*(-j)/m); end end for k=1:n for i=1:n N(k,i)=exp(2*pi*i*k*(-j)/n); end end t=double(t); T = M*t*N; Ni=N; Mi=M; Ti = Mi*T*Ni/(m*n); subplot(2,2,3),imshow(uint8(T)) title(傅里叶变换频谱); subplot(2,2,4),imshow(uint8(Ti)); title(傅里叶逆变换图像); 图像如下： 2.快速傅里叶变换（fft）与逆变换(ifft) 可分离性的优点是二维的傅立叶变换或逆变换由两个连续的一维傅立叶变换变换来实现，对于一个影像f(x,y)，可以先沿着其每一列求一维傅立叶变换，再对其每一行再求一维变换。 正变换 逆变换 由于二维的傅立叶变换具有可分离性，故只讨论一维快速傅立叶变换。 正变换 逆变换 由于计算机进行运算的时间主要取决于所用的乘法的次数。 按照上式进行一维离散由空间域向频率域傅立叶变换时，对于N个F（u)值，中的每一个都要进行N次运算，运算时间与N2成正比。 程序如下： I=imread(D:/chenhao.jpg); subplot(1,2,1),imshow(I); title(原始图像); f=rgb2gray(I); F1=fft2(f); subplot(1,2,2),imshow(F1); title(二维离散傅里叶变换后的频谱); figure,subplot(2,2,1),imshow(log(abs(F1)+1),[0 10]); title(图像的频谱图); F2=fftshift(F1);