matlab数字图像处理课程设计报告,数字图像处理课程设计实验报告.doc
数字图像处理课程设计实验报告
数字图像处理课程设计
题 目：数字图像处理及Huufman(或小波变换)编码仿真实现
学生姓名：
学 院：信息工程学院
系 别：电子信息工程系
专 业：电子信息工程
班 级：电子
指导教师：
2012 年 月 日
内蒙古工业大学课程设计(论文)任务书
课程名称： 数字图像处理课程设计 学院： 信息工程 班级： ___
学生姓名： 学号： 指导教师：
一、题目
数字图像处理及Huufman(或小波变换)编码仿真实现二、目的与意义
数字图像处理课程设计是电子信息工程专业实践教学的主要环节之一。该课程的主要目的是让学生运用已掌握的知识进行实际的图像处理训练。具体包括了对相关知识体系的深入学习以及掌握数字图像直方图统计与绘制的技术能力培养。此外还包含基于数值积分的方法实现平滑算法的设计 并在此基础上完成了基于插值法的图像复原程序的设计与实现
对数字图像直方图进行统计分析并生成分布曲线；采用Matlab开发基于图像平滑算法的去噪系统
图像复原程序设计
创建一个仿真运动/均值模糊PSF来模糊一幅图像(图像自选)。
针对退化设计出复原滤波器，对退化图像进行复原(复原的方法自定)。
针对退化图像进行修复，并展示修复前后的对比图；同时评估修复效果并分析该算法的性能。(FT及其反变换)，快速傅里叶变换(FFT)与其反变换(IFFT)，离散余弦变换(DCT)，小波变换；数字图像直方图的统计与绘制等；设计并实现一种基于边缘检测的图像平滑算法及其在实际中的应用
4.1.1正逆傅里叶变换
img=imread('wxm.jpg');
subplot(2,2,1),imshow(img);
title('原图');
f=rgb2gray(img);
subplot(2,2,2),imshow(f);
title('傅里叶变换');
%二维傅里叶变换
F=fft2(f); %======傅里叶变换
FS=fftshift(F);%======将变换的频率图
%频谱
S=log(1+abs(FS));
subplot(2,2,3);imshow(S,[])
title('直接变换频谱图');;
%% 二维傅里叶逆变换
fr=real(ifft2(ifftshift(FS)));%===频域的图反变
ret=im2uint8(mat2gray(fr)); %===取其灰度图
subplot(2,2,4),imshow(ret);
title('逆傅里叶变换');
处理结果：
4.1.2小波变换
i= imread('wxm.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
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('对角细节系数');
d=idwt2(cA,cH,cV,cD,'db1'); %对分解的细节系数执行单层重构，小波为db1
imshow(d,[]);
处理结果：
4.1.3直方图
I=imread('wxm.jpg');
length(size(I))==3
g=rgb2gray(I);
subplot(1,2,1)
imshow(I);
title('原图');
subplot(1,2,2)
I=rgb2gray(I);
subplot(2,2,1);imshow(I);title('原始图象');
[J,T] = histeq(I);
subplot(1,2,2);imhist(I,64);title('原始图象直方图');
处理结果：
4.1.4离散余弦
clc;
RGB= imread('wxm.jpg');
I = rgb2gray(RGB);
figure;
subplot(2, 2, 1); imshow(I); title('原图像');
J = dct2(I); % 离散余弦变换
J1 = fftshift(log(abs(J)));
subplot(222);imshow(J1);