数字图像处理课程设计汇总.docx
《数字图像处理课程设计汇总.docx》由会员分享，可在线阅读，更多相关《数字图像处理课程设计汇总.docx（23页珍藏版）》请在冰豆网上搜索。
数字图像处理课程设计汇总
学校代码：
10128
学号：
200920203066
数字图像处理课程设计
题目：
数字图像处理及Huufman（或小波变换）编码仿真实现
学生姓名：
肖文杰
学院：
信息工程学院
系别：
电子信息工程系
专业：
电子信息工程
班级：
电子09-1班
指导教师：
吕方王红霞
2012年12月17日
内蒙古工业大学课程设计（论文）任务书
课程名称：
数字图像处理课程设计学院：
信息工程班级：
电子09-1___
学生姓名：
肖文杰___学号：
200920203066_指导教师：
吕芳王红霞
一、题目
数字图像处理及Huufman（或小波变换）编码仿真实现
二、目的与意义
数字图像处理课程设计是电子信息工程专业实践教学的一个重要环节，通过本课程设计使学生运用所学知识进行图像处理的实际训练；是对所学知识的提高和综合应用。
通过综合设计，使得学生进一步理解数字图像处理的基本概念、基本原理和基本方法（如掌握图像增强算法、图像复原以及图像压缩方法等），扩展理论知识，培养学生的综合设计能力。
三、要求（包括原始数据、技术参数、设计要求、图纸量、工作量要求等）
1、预习：
按要求学生可根据自己的情况预习或熟悉所用的语言，搜集资料。
2、分析与设计：
根据选定任务及搜集的资料设计实现方法，确定图像处理算法。
3、程序设计：
运用掌握的语言，编写程序，实现所设计的功能，需要在程序书写时做适当的注释。
4、调试与测试：
自行调试程序，同学之间交叉测试程序，并记录测试情况。
5、验收与评分：
指导教师对每个学生的程序进行综合验收，结合设计报告，根据课程设计成绩评定方法评出成绩。
四、工作内容、进度安排
1、图像处理基本功能（第一周）
普通傅里叶变换（ft）与逆变换（ift）、快速傅里叶变换（fft）与逆变换（ifft）、离散余弦变换（DCT），小波变换；数字图像直方图的统计及绘制等；图像平滑算法实现及应用。
2、图像处理综合功能（第二周）
1）图像复原程序设计
2）给定a,b,c,d概率，进行huffman编码，要求显示原图像、压缩后图像的文件大小、压缩比
五、主要参考文献
1、《图像工程》上册章毓晋清华大学出版社
2、Matlab相关书籍或资料
审核意见
系（教研室）主任（签字）
指导教师下达时间2012年12月14日
指导教师签字：
_______________
数字图像处理课程设计
1、课程设计的目的与意义
通过本课程设计使学生了解数字图像的基本概念，掌握数字图像处理的基本内容，如图像点运算、几何变换、增强处理、图像复原、边缘检测以及图像压缩等的基本原理和Matlab实现方法。
通过本次课程设计，让学生掌握如何学习一门语言，如何进行资料查阅搜集，如何自己解决问题等方法，养成良好的学习习惯。
扩展理论知识，培养学生的综合设计能力。
2、课程设计内容
2.1图像处理基本功能
1）数字图像的变换：
普通傅里叶变换（ft）与逆变换（ift）、快速傅里叶变换（fft）与逆变换（ifft）、离散余弦变换（DCT），小波变换。
2）数字图像直方图的统计及绘制等；
3）基于Matlab的图像平滑算法实现及应用
2.2图像处理综合功能
1）图像复原程序设计
●创建一个仿真运动/均值模糊PSF来模糊一幅图像（图像自选）。
●针对退化设计出复原滤波器，对退化图像进行复原（复原的方法自定）。
●对退化图像进行复原，显示复原前后图像，对复原结果进行分析，并评价复原算法。
2）给定a,b,c,d概率，进行huffman编码，要求显示原图像、压缩后图像的文件大小、压缩比；或采用小波变换进行编码
3、课程设计背景与基本原理
3.1课程设计背景
数字图像处理（DigitalImageProcessing）又称为计算机图像处理，它是指将图像信号转换成数字信号并利用计算机对其进行处理的过程。
是通过计算机对图像进行去除噪声、增强、复原、分割、提取特征等处理的方法和技术。
3.2课程设计基本原理
3.2.1傅里叶变换
傅里叶变换是可分离和正交变换中的一个特例，对图像的傅里叶变换将图像从图像空间变换到频率空间，从而可利用傅里叶频谱特性进行图像处理。
从20世纪60年代傅里叶变换的快速算法提出来以后，傅里叶变换在信号处理和图像处理中都得到了广泛的使用。
3.2.2离散余弦变换
离散余弦变换（discretecosinetransform,DCT）是一种可分离和正交变换并且是对称的。
它与傅里叶变换也有密切的联系，近年得到了广泛应用，特别是在图像压缩领域。
1-D离散余弦变换和其反变换由以下两式定义：
u=0,1,...,N-1
x=0,1,...,N-1
其中a（u）为归一化加权系数，由下式定义：
2-DdDCT对由下面两式定义：
u,v=0,1,...,N-1
x,y=0,1,...,N-1
3.2.3小波变换
小波变换的基础：
序列展开、缩放函数和小波函数。
小波变换和Gabor变换相同的事在时间上和频率上都具有局部化的能力，时间窗口长度与频率窗口长度的乘积是常数。
它与Gabor变换不同的是时间窗口尺寸与频率窗口尺寸互成反变化，对不同的图像都有较好的效果。
2-D小波变换的实现
因为缩放函数和小波函数都是可分离的，所以可对f（x,y）的行先进行1-D变换再对结果进行变换。
下图是小波变换的方框图：
列
小波变换的结果是将图像进行了分解，这种分解是从高尺度向低尺度进行的。
小波分解的结果是将图像划分成了子图像的集合。
在第一级小波分解时，原始图像被划分成了一个低频子图像LL和三个高频子图像HH，LH,和HL的集合。
在第二级小波分解时，低频子图像LL继续被划分成了LL的一个低频子图像和三个高频子图像的集合，而原来第一级分解得到的三个高频子图像不变。
上述分解过程可以这样继续下去，得到越来越多的子图像。
3.2.4哈夫曼编码
哈夫曼编码是消除编码冗余最常用的技术。
当对信源符号逐个编码时，哈夫曼编码能给出最短的码子。
根据无失真编码定理哈夫曼编码方式对固定阶数的信源是最优的。
4、源代码及处理结果
4.1图像处理基本功能
普通傅里叶变换（ft）与逆变换（ift）、快速傅里叶变换（fft）与逆变换（ifft）、离散余弦变换（DCT），小波变换；数字图像直方图的统计及绘制等；图像平滑算法实现及应用。
4.1.1正逆傅里叶变换
傅里叶变换和逆变换
clc;
img=imread（'D:
/xwj1.jpg'）;
subplot（2,2,1）,imshow（img）;
title（'原图'）;
f=rgb2gray（img）;
subplot（2,2,2）,imshow（f）;
b=dct2（a）;
imshow（log（abs（b））,[]）,colormap（jet（64））,colorbar;
title（'DCT变换结果'）;
subplot（224）;
b（abs（b）<10）=0;
%idct
c=idct2（b）/255;
imshow（c）;
title（'IDCT变换结果'）
实验结果：
小波变换
i=imread（'D:
/xwj1.jpg'）;
x=rgb2gray（i）;%真彩色图像转化为灰度图像
[cA,cH,cV,cD]=dwt2（x,'db1'）;%对图像进行单层分解，小波为bd1
A=upcoef2（'a',cA,'db1',1）;
H=upcoef2（'h',cH,'db1',1）;
V=upcoef2（'v',cV,'db1',1）;
D=upcoef2（'d',cD,'db1',1）;%图像编码
figure
（1）
subplot（221）;image（wcodemat（A,192））;
title（'近似细节系数'）;
subplot（222）;image（wcodemat（H,192））;
title（'水平细节系数'）;
subplot（223）;image（wcodemat（V,192））;
title（'垂直细节系数'）;
subplot（224）;image（wcodemat（D,192））;
title（'对角细节系数'）;
figure
（2）
d=idwt2（cA,cH,cV,cD,'db1'）;%对分解的细节系数执行单层重构，小波为db1
imshow（d,[]）;
实验结果：
数字图像直方图的统计及绘制
clc;
Y=imread（'D:
/xwj1.jpg'）;
length（size（Y））==3
s=rgb2gray（Y）;
figure;
subplot（321）;
imshow（Y）;
title（'原图'）;
Y=rgb2gray（Y）;
subplot（323）;
imshow（Y）;
title（'原始图像'）;
[J,T]=histeq（Y）;
subplot（324）;
imshow（J）;
title（'增强图像'）;
subplot（325）;
imhist（Y,64）;
title（'原始图像直方图'）;
subplot（326）;
imhist（J,64）;
title（'均衡化图像直方图'）;
实验结果：
平滑算法
M=imread（'D:
/xwj1.jpg'）;
I=rgb2gray（M）;%彩色图像变灰度图像
J=imnoise（I,'salt&pepper',0.02）;
subplot（221）,imshow（I）
title（'原图像'）
subplot（222）,imshow（J）
title（'添加椒盐噪声图像'）
K1=filter2（fspecial（'average',3）,J）/255;%应用3*3邻域窗口法
subplot（223）,imshow（K1）
title（'3x3窗的邻域平均滤波图像'）
K2=filter2（fspecial（'average',7）,J）/255;%应用7*7邻域窗口法
subplot（224）,imshow（K2）
title（'7x7窗的邻域平均滤波图像'）
实验结果：
4.2图像处理综合功能
1）图像复原程序设计
●创建一个仿真运动/均值模糊PSF来模糊一幅图像（图像自选）。
●针对退化设计出复原滤波器，对退化图像进行复原（复原的方法自定）。
●对退化图像进行复原，显示复原前后图像，对复原结果进行分析，并评价复原算法。
2）给定a,b,c,d概率，进行huffman编码，要求显示原图像、压缩后图像的文件大小、压缩比；或采用小波变换进行编码
运动模糊及图像复原
I=imread（'D:
/xwj1.jpg'）;
subplot（221）;imshow（I）;
title（'原图像'）;
PSF=fspecial（'motion',12,10）;
MF=imfilter（I,PSF,'circular'）;
noise=imnoise（zeros（size（I））,'gaussian',0,0.001）;
MFN=imadd（MF,im2uint8（noise））;
subplot（222）;imshow（MFN）;
title（'运动模糊图像'）;
Sn=abs（fft2（（MFN））.^2）;
Sf=abs（（fft2（I））.^2）;
NCORR=（real（ifft2（Sn）））;
ICORR=
升级会员 