基于MATLAB课程设计-图像处理完整版

摘要

21世纪是一个充满信息的时代,图像作为人类感知世界的视觉基础,是人类获取信息、表达信息和传递信息的重要手段。图像处理,是用计算机对图像进行分析,以达到所需结果的技术。又称影像处理。基本内容 图像处理一般指数字图像处理。数字图像是指用数字摄像机、扫描仪等设备经过采样和数字化得到的一个大的二维数组,该数组的元素称为像素,其值为一整数,称为灰度值。图像处理技术的主要内容包括图像压缩,增强和复原,匹配、描述和识别3个部分。 常见的处理有图像数字化、图像编码、图像增强、图像复原、图像分割和图像分析等。图像处理一般指数字图像处理。所谓数字图像处理[7]就是利用计算机对图像信息进行加工以满足人的视觉心理或者应用需求的行为。实质上是一段能够被计算机还原显示和输出为一幅图像的数字码。

关键词:DCT变换;图像压缩;真色彩增强;平滑;锐化;直方图均衡;

灰度变换;滤波;M文件的使用

目录

摘要………………………………………………………………I

1 概述……………………………………………………………II

2 课程设计任务及要求...............................III

  设计任务

  设计要求

3 系统设计原理.....................................Ⅳ

   DCT图像压缩原理

   真彩色增强

      平滑

      锐化

   灰度变换(直方图均衡化)

   图像滤波

      中值滤波器

      维纳滤波器

4 程序代码及实验结果与分析..............................Ⅵ

   DCT图像压缩

      程序代码

      实验结果

      结果分析

   真彩色增强

      平滑程序代码

      实验结果

      结果分析

      锐化程序代码

      实验结果

      结果分析

   灰度变换(直方图均衡化)

   程序代码

      实验结果

      结果分析

   图像滤波

      程序代码

      实验结果

      结果分析

5 收获体会.........................................Ⅶ

6 参考文献........................................Ⅷ

概述

MATLAB作为一种矩阵语言,进行数字图像处理是非常方便的。从理论上讲,数字图像是一个二维的整数阵列。掌握简单的MATLAB图像增强技术,从而对图像处理有进一步的了解。

课程设计任务及要求

设计任务

  1. 在图像的变换和压缩中,常常用到离散余弦变换(DCT)。DCT变换用于图像的压缩实例。请在测试图像中验证你的结论。
  2. 请编程实现图像的真彩色增强。
  3. 通过直方图均衡化的方法实现图像的灰度变换,在测试图像中验证你的结论,分析程序结果。
  4. 使用常用的滤波器对数字图像进行处理。如均值滤波器、中值滤波器、维纳滤波器等。

设计要求

在报告中给出本实验的实验结果和分析。

    1. 程序代码附有必要的注释。
    2. 报告中对实验的基本算法原理要进行简要叙述。

系统设计原理

DCT图像压缩原理

DCT变换是最小均方误差条件下得出的次最佳正交变换,且已获得广泛应用,并成为许多图像编码国际标准的核心。JPEG图像格式的压缩算法采用的就是DCT变换,DCT变换的变换核为余弦函数,计算速度较快,有利于图像压缩和其他处理。在编码过程中,JPEG算法首先将RGB分量转化为亮分量和色差分量,然后将图像分解为8*8的像素块,,对这个8*8块进行二维离散余弦变换,每个块就产生了64个DCT系数,其中一个是直流(DC),它表示了8*8输入矩阵全部值的平均数,其余63个系数为交流(AC)系数,接下来对DCT系数进行量化,最后将量化的DCT系数进行编码,就形成了压缩后的图像格式。在解码过程中,先对已编码的量化的系数进行解码,然后求逆量化并利用二维DCT反变换把DCT系数转化为8*8样本像块,最将反变换后的块组合成一幅图像。这样就完成了图像的压缩和解压过程。

离散余弦变换DCT的MATLAB实现有两种方 法,一种是基于FFR的快速算法,这是通过MATLAB工具箱提供的DCT2函数实现的;另一种是DCT变换是矩阵方法。变换矩阵方法非常适合做8*8或16 *16的图像块的DCT变换,工具箱提供了dctmtx函数来计算变换矩阵。

真彩色增强

真彩色增强主要是针对伪彩色增强而言的。图像的色彩增强技术主要分为为彩色增强和真彩色增强两种,这两种方法在原理上存在着本质的区别。伪彩色增强时对原灰度图像中不同灰度值区域分别付于不同的颜色,使人能够更明白的区分不同的灰度级。由于原始图像事实上是没有颜色的,所以称这种人工赋予的颜色为伪彩色,伪彩色增强实质上只是一个图像的着色过程是一种灰度到彩色的映射技术。真彩色增强则是对原始图像本身具有的颜色进行调节,是一个彩色到彩色的映射过程。

平滑

在图像中,通过相邻点的相互平均可以去掉一些突然变化的点,从而滤掉一定的噪声,达到平滑的目的,使图片看起来更柔和,颜色更均匀更清晰。

锐化

图像平滑往往使图像中的边界、轮廓变得模糊,为了减少这类不利效果的影响,需要用图像锐化技术使图像的边缘变的清晰。图像锐化处理的目的是为了使图像的边缘、轮廓线以及图像的细节变的清晰,经过平滑的图像变得模糊的根本原因是因为图像受到了平均或积分运算,因此可以对其进行逆运算(如微分运算)就可以使图像变的清晰。从频率域来考虑,图像模糊的实质是因为其高频分量被衰减,因此可以用高通滤波器来使图像清晰。在水下图像的增强处理中除了去噪,对比度扩展外,有时候还需要加强图像中景物的边缘和轮廓。而边缘和轮廓常常位于图像中灰度突变的地方,因而可以直观地想到用灰度的差分对边缘和轮廓进行提取。

灰度变换(直方图均衡化)

   直方图均衡化的基本思想是把原始图的直方图变换为均匀分布的形式,这样就增加了象素灰度值的动态范围从而可达到增强图像整体对比度的效果。设原始图像在(x,y)处的灰度为f,而改变后的图像为g,则对图像增强的方法可表述为将在(x,y)处的灰度f映射为g。在灰度直方图均衡化处理中对图像的映射函数可定义为:g = EQ (f),这个映射函数EQ(f)必须满足两个条件(其中L为图像的灰度级数):

    (1)EQ(f)在0≤f≤L-1范围内是一个单值单增函数。这是为了保证增强处理没有打乱原始图像的灰度排列次序,原图各灰度级在变换后仍保持从黑到白(或从白到黑)的排列。

(2)对于0≤f≤L-1有0≤g≤L-1,这个条件保证了变换前后灰度值动态范围的一致性。

图像滤波

  滤波器是一种选频装置可以使信号中特定的频率成分通过而极大地衰减其它频率成分,可以滤除干扰噪声。在数字图像处理中,常常会遇到图像中混杂有许多的噪声。因此,在进行图像处理中,有时要先进行祛除噪声的工作。最常用的祛除噪声的方法是用滤波器进行滤波处理。MATLAB的图像处理工具箱里也设计了许多的滤波器。如均值滤波器、中值滤波器、维纳滤波器等

维纳滤波器

   维纳滤波器(Wiener filter)是由数学家维纳(Rorbert Wiener)提出的一种以最小平方为最优准则的线性滤波器。在一定的约束条件下,其输出与一给定函数(通常称为期望输出)的差的平方达到最小,通过数学运算最终可变为一个托布利兹方程的求解问题。维纳滤波器又被称为最小二乘滤波器或最小平方滤波器,目前是基本的滤波方法之一。维纳滤波是利用平稳随机过程的相关特性和频谱特性对混有噪声的信号进行滤波的方法。

中值滤波器:

   中值滤波是一种非线性数字滤波器技术,经常用于去除图像或者其它信号中的噪声。这个设计思想就是检查输入信号中的采样并判断它是否代表了信号,使用奇数个采样组成的观察窗实现这项功能。观察窗口中的数值进行排序,位于观察窗中间的中值作为输出。然后,丢弃最早的值,取得新的采样,重复上面的计算过程。 中值滤波是图像处理中的一个常用步骤,它对于斑点噪声(en:speckle noise)和椒盐噪声(en:salt-and-pepper noise)来说尤其有用。保存边缘的特性使它在不希望出现边缘模糊的场合也很有用。

程序代码及实验结果与分析

DCT图像压缩

程序代码:

I=imread('压缩前');%读取原始图片

I=rgb2gray(I);%把图像变成灰度图像

I=im2double(I);%把图像数据类型转换为双精度浮点类型

T=dctmtx(8);%计算离散余弦变换

B=blkproc(I,[8 8],'P1*x*P2',T,T');%对图像I的每个不同8*8数据块应用矩阵式’P1*X *P2’进行处理,其中P1=T,P2=T'

mask=[1  1  1  1  0  0  0  0

      1  1  1  0  0  0  0  0

      1  1  0  0  0  0  0  0

      1  0  0  0  0  0  0  0

      0  0  0  0  0  0  0  0

      0  0  0  0  0  0  0  0

      0  0  0  0  0  0  0  0

      0  0  0  0  0  0  0  0];%二进制值掩码,用来压缩DCT的系数;可通过改变矩阵来改变压缩程度

B2=blkproc(B,[8 8],'P1.*x',mask);%只保留DCT变换的10个系数

I2=blkproc(B2,[8 8],'P1*x*P2',T',T);%逆DCT变换,用来重构图像

subplot(1,2,1);

imshow(I);title('压缩前');%显示原始图像

subplot(1,2,2);

imshow(I2);title('压缩后');%显示压缩图像

imwrite(I2,'D:\我的文档\My Pictures\压缩后')

实验结果:

结果分析:

   尽管抛弃了85%的DCT系数,压缩后占用空间减少了约75%,但是处理后的图片与原图的失真并不大,压缩图片的画质令人满意。由此也验证了JPEG格式由于占用空间小,图片质量高被广泛认可。

真彩色增强

平滑程序代码:

I=imread('');

subplot(1,2,1);imshow(I);title('原图');

h1=fspecial('average',3);%用均值滤波器是领域值平均,以达到平滑目的

I1=imfilter(I,h1);

subplot(1,2,2);imshow(I1);title('平滑滤波真彩色增强');

实验结果:

结果分析:

处理后的图片较处理前的图片看起来更柔和,色彩更均匀。

锐化程序代码:

I=imread('');

g=[0 0 0;0 2 0;0 0 0];

w=fspecial('laplacian',;

H=im2double(I);

J=imfilter(H,g,'replicate');

subplot(1,2,1);

imshow(H);title('原图');

subplot(1,2,2);imshow(J);title('锐化真彩色增强图像');

实验结果

结果分析:

可以很明显的看出处理后的图片轮廓更清晰,边界也更明显

灰度变换(直方图均衡化):

程序代码

I=imread('压缩前');

I=rgb2gray(I);

subplot(2,2,1);imshow(I);title('原图');

subplot(2,2,2);imhist(I);title('原直方图');

J=histeq(I);%自动直方图均衡化,均衡后只有64个灰度级

subplot(2,2,4);imhist(J);title('灰度变换直方图');

subplot(2,2,3);imshow(J);title('灰度变换图');%颜色分布相对均衡

实验结果:

结果分析:

   对比处理前后的直方图可知,原直方图中在灰度级100左右的值出现的平率很低可近似忽略而经过histeq函数后直方图的分布从0~200出现的频率相近,基本实现直方图的均匀分布;从图片上则体现为变换后的图片颜色分布均匀没有特别黑或特别白的地方。另如果将histeq中的n值改为200(原图的灰度级),直方图的分布则接近为原图(如下图),此一来更好的说明了直方图可以明显直观地改变图像的灰度级。

图像滤波

程序代码

I=imread('压缩前');

I=rgb2gray(I);

I=im2double(I);

J=imnoise(I,'gaussian',0,; %加入高斯噪声

%进行均值滤波

h=fspecial('average',3);   %fspecial函数用于产生预定义滤波器

I2=uint8(round(filter2(h,I)));  %filter2函数用于图像滤波,h是滤波参数(均值)

%进行中值滤波

I3=medfilt2(J);  %medfilt2函数用于图像的中值滤波

%进行维纳滤波

I4=wiener2(J);%进行一次维纳滤波

I5=wiener2(I4);%进行二次维纳滤波

subplot(2,3,1),imshow(I),title('原图象')

subplot(2,3,2),imshow(J),title('加噪声图象')

subplot(2,3,3),imshow(I2),title('均值滤波后图象')

subplot(2,3,4),imshow(I3),title('中值滤波后图象')

subplot(2,3,5),imshow(I4),title('维纳滤波后图象')

subplot(2,3,6),imshow(I5),title('两次维纳滤波后图象')

实验结果:

结果分析

从图上对比可看出中值滤波器和维纳滤波器滤除噪声的能力差不多,但中值滤波器对于图片亮度的改变相对于维纳滤波器的小。可以看出维纳滤波器不仅滤除了噪声而且也降低了图片的亮度。二次滤波后噪声更小了。

参考文献

[1]  张兆礼,赵春晖,梅晓丹.现代图像处理技术及MATLAB实现.北京:人民      邮电出版社,2001

[2]  孙兆林.MATLAB图像处理.北京:清华大学出版社,2002

[3]  朱习军,隋思涟,张宾,刘尊年.MATLAB在信号与图像处理中的应用.北京,      电子工业出版社

[4]  胡晓军、 徐飞 .MATLAB应用图像处理(第二版).西安电子科技大学出版

[5]  高展宏、 徐文波. 基于MATLAB的图像处理案例教程 清华大学出版社

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.xdnf.cn/news/4415.html

如若内容造成侵权/违法违规/事实不符,请联系一条长河网进行投诉反馈,一经查实,立即删除!

相关文章

基于51单片机的步进电机定时控制proteus仿真

地址: https://pan.baidu.com/s/1_F5bfyS_e_eKSblnja7RqA 提取码:1234 仿真图: 芯片/模块的特点: AT89C52/AT89C51简介: AT89C52/AT89C51是一款经典的8位单片机,是意法半导体(STMicroelectro…

【云原生开发】如何通过client-go来操作K8S集群

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,…

idea 配置自动导入设置

配置自动导入设置: 通过访问 Settings→Editor→General→Auto Import,勾选 Add unambiguous imports on the fly 和 Optimize imports on the fly 选项,可以实现自动添加明确的导入和快速优化导入,即自动删除无用的导入

【文本情感分析识别】Python+SVM算法+模型训练+文本分类+文本情感分析

一、介绍 使用Python作为开发语言,基于文本数据集(一个积极的xls文本格式和一个消极的xls文本格式文件),使用Word2vec对文本进行处理。通过支持向量机SVM算法训练情绪分类模型。实现对文本消极情感和文本积极情感的识别。并基于D…

iptables面试题

1、详述iptales工作流程以及规则过滤顺序? iptables过滤的规则顺序是由上至下,若出现相同的匹配规则则遵循由上至下的顺序 2、iptables的几个表以及每个表对应链的作用? Iptables有四表五链 Filter表 : Filter表是iptables中使用的默认表…

【青牛科技】GC2803:白色家电与安防领域中 ULN2803 的卓越替代者

在当今科技飞速发展的时代,电子元器件在各个领域都扮演着至关重要的角色。在白色家电和安防等产品的电路设计中,驱动芯片的选择尤为关键。传统的 ULN2803 曾是广泛应用的一款芯片,但如今,芯麦 GC2803 的出现为这些领域带来了新的选…

变脸、看完毒液3,我把自己的脸变毒液了,视频有点惊悚!请谨慎观看

大家好,我是Shelly,一个专注于输出AI工具和科技前沿内容的AI应用教练,体验过300款以上的AI应用工具。关注科技及大模型领域对社会的影响10年。关注我一起驾驭AI工具,拥抱AI时代的到来。 AI工具集1:大厂AI工具【共23款…

数据分析:微生物功能差异分析之Maaslin2

文章目录 介绍加载R包数据链接导入数据数据预处理Maaslin2检验输出结果画图系统信息介绍 Maaslin2(Microbiome Multivariable Associations with Linear Models)是一种用于微生物组学数据的多变量关联分析工具。它的原理和特点如下: 多变量关联分析: Maaslin2旨在高效确定…

电脑软件:推荐四款非常好用的电脑磁盘分析工具

一、WizTree WizTree 是一款Windows下磁盘空间分析工具。它可以快速扫描并分析你的电脑硬盘驱动器中文件和文件夹,并以可视化块状的方式展示哪些文件和文件夹使用的磁盘空间最多。这样你就可以很方便找到占用磁盘空间大的文件。 WizTree的特点 ● 磁盘空间利…

6. STM32之TIM实验--编码器接口()--(实验5:PWM驱动直流电机)

这篇文章是通用定时器的最后一章节,也就是编码器接口,主要是用来进行对精确测量旋转角度或速度的负载进行精确控制。 STM32 编码器模式详解-CSDN博客 STM32——编码器测速原理及STM32编码器模式_龙邱512编码器stm32历程-CSDN博客 代码可根据这个进行编…

如何批量创建文件夹并命名?6个一键批量创建的方法

如何批量创建文件夹并命名?在快节奏的现代工作环境中,时间成为了最宝贵的资源。面对海量的数据与文件,如何高效地组织与管理,成为了提升工作效率的关键。为了节省时间,批量创建文件夹并命名,成为了一项至关…

人工智能之人脸识别(face_recognition)

文章目录 face_recognition 介绍主要功能**与opencv联系联系检测人脸切割人脸提取人物关键特征计算人脸的欧几里得距离计算人脸匹配程度总结 face_recognition 介绍 face_recognition 介绍 face_recognition 是一个非常流行的 Python 库,专门用于人脸识别任务。它基…

Java学习路线:JUL日志系统(二)使用Properties配置文件

目录 认识properties 使用properties编写日志配置文件 认识properties 之前的学习中,我们学习了使用XML配置文件,但是XML的读取实在有些麻烦。那有没有更简单的方式来配置文件呢? 答案是:使用Properties配置文件 在这里了解pro…

【数据库】elasticsearch

1、架构 es会为每个索引创建一定数量的主分片和副本分片。 分片(Shard): 将索引数据分割成多个部分,每个部分都是一个独立的索引。 主要目的是实现数据的分布式存储和并行处理,从而提高系统的扩展性和性能。 在创建索…

Qt——常用控件

前言:本篇文章,将分享Qt中常用的,具有代表性的一些控件。 一.按钮类控件 在前边的文章中我们也多次分享过PushButton按钮,但Qt中并非只提供这一种按钮。 在Qt中,QPushButton并非直接继承自QWidget,QAbstr…

硬件基础06 滤波器——无源、有源(含Filter Solutions、Filter Pro、MATLAB Fdatool)

推荐设计RC无源滤波器使用,数字滤波器可以使用MATLAB,有源滤波器使用Filter Pro。 一、Filter Solutions 1、软件资源及安装教程如下 FilterSolutions14.1.rar资源 (1)、双击FS14_1_0.exe进行安装;不要安装到中文路…

js WebAPI黑马笔记(万字速通)

此笔记来自于黑马程序员,pink老师yyds 复习: splice() 方法用于添加或删除数组中的元素。 注意: 这种方法会改变原始数组。 删除数组: splice(起始位置, 删除的个数) 比如:1 let arr [red, green, b…

docker安装低版本的jenkins-2.346.3,在线安装对应版本插件失败的解决方法

提示:写完文章后,目录可以自动生成,如何生成可参考右边的帮助文档 文章目录 前言一、网上最多的默认解决方法1、jenkins界面配置清华源2、替换default.json文件 二、解决低版本Jenkins在线安装插件问题1.手动下载插件并导入2.低版本jenkins在…

spark-on-k8s 介绍

spark-on-k8s 介绍 摘要 最近一段时间都在做与spark相关的项目,主要是与最近今年比较火的隐私计算相结合,主要是在机密计算领域使用spark做大数据分析、SQL等业务,从中也了解到了一些spark的知识,现在做一个简单的总结&#xff…

React教程(详细版)

React教程(详细版) 1,简介 1.1 概念 react是一个渲染html界面的一个js库,类似于vue,但是更加灵活,写法也比较像原生js,之前我们写出一个完成的是分为html,js,css&…