预览加载中,请您耐心等待几秒...
1/9
2/9
3/9
4/9
5/9
6/9
7/9
8/9
9/9

在线预览结束,喜欢就下载吧,查找使用更方便

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

齐齐哈尔大学 通信与电子工程学院 实验教师预作报告 实验课名称:通信系统仿真设计实验 项目名称:电视机模型建模与仿真(2) 教学班级:通信10级 教师姓名:房汉雄 实验地点:6J510 实验日期:2013.03.16 实验十二图像处理的仿真 实验目的及要求 1.了解图像增强的目的及意义,加深对图像增强的感性认识,巩固所学的图像增强的理论知识和相关算法。 2.熟练掌握直方图均衡化和直方图规定化的计算过程。 3.熟练掌握空域滤波中常用的平滑和锐化滤波器。 4.熟练掌握低通和高通滤波器的使用方法,明确不同性质的滤波器对图像的影响和作用。 二、实验环境 计算机,带Windows操作系统 三、实验原理 1.直方图均衡化 直方图均衡化处理是以累积分布函数为基础的直方图修改方法。直方图均衡化的目的是将原始图像的直方图变为均衡分布的形式,即将一幅已知灰度概率密度分布的图像,经过某种变换,变成一幅具有均匀灰度概率密度分布的新图像。相关的直方图均衡化的计算过程请查看教材。 2.均值滤波器和中值滤波器 邻域平均法将一个像元及其邻域中的所有像元的平均值赋给输出图像中相应的像元,从而达到平滑的目的,又称为均值滤波。以3*3模板为例,将图1(b)中的所有赋予相同的值,求得r后将其值除以9,再赋给相应的像元。 r(a)(b)(c) 图1空域模板滤波示意图 邻域滤波法虽然可以平滑图像,但在消除噪声的同时,会使图像中的一些细节变得模糊。中值滤波则在消除噪声的同时还能保持图像中的细节部分,防止边缘模糊,与邻域平均法不同,中值滤波是一种非线性滤波,它将邻域内所有像元值从小到大排序,取中间值作为中心像元的输出值,以3*3模板为例,将图4.1(a)中的~进行排序,将中间值赋给r。中值滤波与均值滤波不同,它不是通过对邻域内的所有像元求平均值来消除噪声的,而是让与周围像元灰度值的差比较大的像元改取近似于周围像元灰度值的值,从而达到消除噪声的目的。 3.锐化滤波器 图像的平滑处理会使图像的边缘纹理信息受到损失,图像变得比较模糊。锐化滤波则可以增强图像中物体的边缘轮廓并减弱图像的低频分量,使得除边缘以外的像元的灰度值趋向于零。常用的锐化滤波方法有梯度法和拉普拉斯算子法。 (1)梯度法 对于一幅图像,它的梯度可以表示为: 和分别表示在某个位置,沿x方向和y方向的灰度变化率(一阶偏导数),是一个矢量,它指向的最大变化率方向。对于数字图像用微分运算不方便,故一般采用差分形式,常用的梯度差分法为罗伯特梯度算子,可以表示为: 图4.6给出了罗伯特算子的模板。 (2)拉普拉斯算子 拉普拉斯算子法采用二阶偏导数,其定义为: 对于离散数字图像,在某个像元(i,j)处的拉普拉斯算子可采用如下差分形式: 其相应的模板如图2所示。 10 010-1-10罗伯特算子 0-10-1-1-1-14-1-18-10-10-1-1-1拉普拉斯算子 图2算子模板 4.频域滤波增强 频域增强技术是在图像的频率域空间对图像进行滤波,其主要步骤如下: (a)、对原始图像进行傅里叶变换得到; (b)、将与传递函数进行乘法运算得到; (c)、将进行傅里叶反变换得到增强图像。 图像从空间域变换到频率域后,其低频分量对应图像中灰度值变化比较缓慢的区域。低通滤波是指保留低频分量,而通过滤波器函数减弱或抑制高频分量的过程。因此低通滤波与空域中的平滑滤波器一样可以消除图像中的随机噪声,减弱边缘效应,起到平滑图像的作用。高通滤波是指保留高频分量,而减弱或抑制低频分量的过程。图像频率域的高频分量表征了图像的边缘及其它灰度变化较快的区域,因此高通滤波可以增强图像的边缘,起到锐化图像的作用。 四、实验内容及步骤: 1.编制一个对图像进行直方图均衡化处理的程序,增强的结果要能实时在屏幕上显示。 x=imread('w.bmp');x=rgb2gray(x); y=histeq(x); figure('name','直方图均衡化','NumberTitle','off','MenuBar','none'); subplot(1,2,1);imshow(x);title('原图像'); subplot(1,2,2);imshow(y);title('直方图均衡化后图像'); 2.编制程序实现m×n的均值滤波器和中值滤波器,注意边界像元的处理,设计一种简便的对图像边界像元进行处理的方法。 %%邻域平均法 %读取原始图像并显示 f0=imread('w.bmp'); %加入高斯噪声并显示 x=imnoise(f0,'gaussian',0,0.05);x=rgb2gray(x); %给定三种不同的滤波窗口 h1=[1/91/91/9;1/91/91/9;1/91