数字图像处理课程设计实验报告

数字图像处理课程设计

题 目：数字图像处理及Huufman(或小波变换)编码仿真实现

学生姓名：

学 院：信息工程学院

系 别：电子信息工程系

专 业：电子信息工程

班 级：电子

指导教师：

2012 年 月 日

内蒙古工业大学课程设计(论文)任务书

课程名称： 数字图像处理课程设计 学院： 信息工程 班级： ___

学生姓名： 学号： 指导教师：

一、题目

数字图像处理及Huufman(或小波变换)编码仿真实现二、目的与意义

数字图像处理课程设计是电子信息工程专业实践教学的一个重要环节，通过本课程设计使学生运用所学知识进行图像处理的实际训练；是对所学知识的提高和综合应用数字图像直方图的统计及绘制等；图像平滑算法实现及应用图像复原程序设计

数字图像直方图的统计及绘制等；基于Matlab的图像平滑算法实现及应用

图像复原程序设计

创建一个仿真运动/均值模糊PSF来模糊一幅图像(图像自选)。

针对退化设计出复原滤波器，对退化图像进行复原(复原的方法自定)。

对退化图像进行复原，显示复原前后图像，对复原结果进行分析，并评价复原算法。(ft)与逆变换(ift)、快速傅里叶变换(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);