文档章节

Matlab数字图像处理之三:图像算术运算及几何变换

天蚕宝衣
 天蚕宝衣
发布于 2017/03/16 20:29
字数 1495
阅读 395
收藏 0

一、 实验目的

1. 掌握算术运算、几何变换的方法;

2. 了解图像算术运算、几何变换在数字图像处理中的初步应用;

3. 观察算法(变换)处理结果,体会图像算术运算和几何变换处理的过程和处理前后图像的变化。

二、 实验原理

1. 图像算术运算:两幅输入图像之间进行的点对点的加、减、乘、除运算后得到输出图像的过程。如果输入图像为A(x,y)和B(x,y),输出图像为C(x,y),则图像的算术运算有如下四种形式:

C(x,y) = A(x,y) + B(x,y)

C(x,y) = A(x,y) – B(x,y)

C(x,y) = A(x,y) * B(x,y)

C(x,y) = A(x,y) / B(x,y)

2. MATLAB图像处理工具箱中的算术运算函数

matlab图像算术运算及几何变换

3. 图像的算术运算函数使用以下截取规则使运算结果符合数据范围的要求:

超出数据范围的整型数据将被截取为数据范围的极值,分数结果将被四舍五入。例如,如果数据类型是uint8,那么大于255的结果(包括无穷大inf)将被设置为255。

三、 实验内容及步骤

1. 图像的加法运算

imadd函数调用格式:Z=imadd(X,Y)

实现:图像+图像,图像+常数

2. 噪声处理

Imnoise函数调用格式:Z=imnoise(X, type, parameter)

type:噪声的类型 parameter:噪声的参数

3. 图像的减法运算(差分方法)

用途:检测图像变化及运动物体,作为许多图像处理工作的准备步骤。

函数调用格式: Z = imsubtract(X,Y)

Z = imabsdiff(X,Y)

实现:图像-图像,图像-常数

4. 图像的乘法运算

用途:

实现掩模操作,即屏蔽掉图像的某些部分;

实现缩放:一幅图像乘以一个常数,如果使用的缩放因子大于1,那么将增强图像的亮度,如果因子小于1则会使图像变暗,效果优于简单添加像素偏移量。

由于时域的卷积或相关运算与频域的乘积运算对应,因此乘法运算有时也被作为一种技巧来实现卷积或相关处理。

immultiply函数调用格式:Z = immultiply(X,Y)

5. 图像的除法运算

用途:

校正成像设备的非线性影响,这在特殊形态的图像(如断层扫描等医学图像)处理中常常用到。

检测两幅图像间的区别——相应像素值的变化比率,而不是每个像素的绝对差异,因而称为比率变换。

imdivide函数的调用格式如下:Z = imdivide(X,Y)

6. 图像的四则运算

函数imlincomb按照双精度执行所有算术运算操作,而且仅对最好的输出结果进行截取,该函数的调用格式如下

Z = imlincomb(A,X,B,Y,C);

其中,Z=A*X+B*Y+C。

MATLAB会自动根据输入参数的个数判断需要进行的运算。例如:

Z = imlincomb(A,X,C) 计算Z=A*X+C:

Z = imlincomb(A,X,B,Y,) 计算Z=A*X+B*Y:

7. 图像的缩放

图像缩放函数imresize调用格式:

B = imresize(A,m)

B = imresize(A,m,method)

B = imresize(A,[mrows ncols],method)

B = imresize(…,method,n)

B = imresize(…,method,h(FIR滤波器))

Method:插值方法

nearest—最近邻插值(默认)

Bilinear—双线性插值

Bicubic—三次卷积插值

8. 图像的旋转

图像旋转函数imrotate调用格式:

B = imrotate(A,angle,method)
B = imrotate(A,angle,method,’crop’)

四、 实验结果

1. 图像的加法运算

实现:图像+图像

A=imread(‘rice.png’)

B=imread(‘testpat1.png’)

C=imadd(A,B)

subplot(1,3,1),imshow(A);subplot(1,3,2),imshow(B);subplot(1,3,3),imshow(C)

matlab图像算术运算及几何变换

实现:图像+/-常数

A=imread(‘cameraman.tif’)

B=A+50

C=A-50

subplot(1,3,1),imshow(A);subplot(1,3,2),imshow(B);subplot(1,3,3),imshow(C)

matlab图像算术运算及几何变换

2. 噪声处理

对图像rice.png分别加入高斯噪声、椒盐噪声和乘性噪声:

A=imread(‘rice.png’)

B=imnoise(A,’gaussian’,0.02)

C=imnoise(A,’salt & pepper’,0.02)

D=imnoise(A,’speckle’,0.02)

subplot(2,2,1),imshow(A) ,title(‘原图’);

subplot(2,2,2),imshow(B) ,title(‘高斯噪声’);

subplot(2,2,3),imshow(C) ,title(‘椒盐噪声’);

subplot(2,2,4),imshow(D) ,title(‘乘性噪声’);

matlab图像算术运算及几何变换

在原图上分别叠加密度为0.01~0.06的椒盐噪声后取平均

A=imread(‘pout.tif’);

B=imnoise(A,’salt & pepper’,0.01);

C=imnoise(A,’salt & pepper’,0.02);

D=imnoise(A,’salt & pepper’,0.03);

E=imnoise(A,’salt & pepper’,0.04);

F=imnoise(A,’salt & pepper’,0.05);

G=imnoise(A,’salt & pepper’,0.06);

H=0.1*B+0.2*C+0.1*D+0.2*E+0.2*F+0.2*G;

subplot(2,4,1),imshow(A) ,title(‘原图’);

subplot(2,4,2),imshow(B) ,title(‘椒盐0.01′);

subplot(2,4,3),imshow(C) ,title(‘椒盐0.02′);

subplot(2,4,4),imshow(D) ,title(‘椒盐0.03′);

subplot(2,4,5),imshow(E) ,title(‘椒盐0.04′);

subplot(2,4,6),imshow(F) ,title(‘椒盐0.05′);

subplot(2,4,7),imshow(G) ,title(‘椒盐0.06′);

subplot(2,4,8),imshow(H) ,title(‘椒盐平均’);

matlab图像算术运算及几何变换

3. 图像的减法运算(差分方法)

A=imread(‘cameraman.tif’)

B =imopen (A, strel(‘disk’,15))

C=imabsdiff(A,B)

subplot(1,3,1),imshow(A);subplot(1,3,2),imshow(B);subplot(1,3,3),imshow(C)

matlab图像算术运算及几何变换

4. 图像的乘法运算

A=imread(‘cameraman.tif’)

B=immultiply(A,0.5)

C=immultiply(A,1.2)

subplot(1,3,1),imshow(A),title(‘原图’);

subplot(1,3,2),imshow(B),title(’0.5′);

subplot(1,3,3),imshow(C),title(’1.2′)

matlab图像算术运算及几何变换

5. 图像的除法运算

A=imread(‘cameraman.tif’)

B =imopen (A, strel(‘disk’,15))

C=imabsdiff(A,B)

D= imdivide(A,C)

subplot(1,3,1),imshow(A),title(‘原图’);

subplot(1,3,2),imshow(C),title(‘减去背景的图像’);

subplot(1,3,3),imshow(D),title(‘原图和减去背景图像相除’)

matlab图像算术运算及几何变换

6. 图像的四则运算

A=imread(‘cameraman.tif’)

B =imopen (A, strel(‘disk’,15))

C=imabsdiff(A,B)

D= imdivide(A,C)

E= imlincomb(0.2,C,0.4,D,12)

subplot(1,3,1),imshow(C),title(‘C’);

subplot(1,3,2),imshow(D),title(‘D’);

subplot(1,3,3),imshow(E),title(’0.2*C+0.4*D+12′)

matlab图像算术运算及几何变换

7. 图像的缩放

A = imread(‘cameraman.tif’)

B = imresize(A,0.5,’nearest’)

C = imresize(A,0.8,’Bilinear’)

D = imresize(A,1.2,’Bicubic’)

figure,imshow(B),title(‘最近邻差值 0.5′);

figure,imshow(C),title(‘双线性插值 0.8′);

figure,imshow(D),title(‘三次卷积插值1.2′)

matlab图像算术运算及几何变换

8. 图像的旋转

A = imread(‘cameraman.tif’)

B = imrotate(A,30, ‘bilinear’,’crop’)

C = imrotate(A,30, ‘bilinear’,’loose’)

subplot(1,3,1),imshow(A),title(‘原图’);

subplot(1,3,2),imshow(B),title(‘crop’);

subplot(1,3,3),imshow(C),title(‘loose’)

matlab图像算术运算及几何变换

本文转载自:http://www.jishuzh.com/gis/matlab%E6%95%B0%E5%AD%97%E5%9B%BE%E5%83%8F%E5%A4%84%E7%90%86%E4%B9%8B%E4%

共有 人打赏支持
天蚕宝衣
粉丝 24
博文 240
码字总数 179371
作品 0
天津
私信 提问
[Python图像处理] 十.形态学之图像顶帽运算和黑帽运算

版权声明:本文为博主原创文章,转载请注明CSDN博客源地址!共同学习,一起进步~ https://blog.csdn.net/Eastmount/article/details/83692456 该系列文章是讲解Python OpenCV图像处理知识,前...

Eastmount
2018/11/04
0
0
【工具使用系列】关于 MATLAB 形态学,你需要知道的事

如何进行图像操作 二值形态学基本运算 膨胀和腐蚀 膨胀 腐蚀 膨胀与腐蚀的结合 膨胀与腐蚀的对偶性 开启和闭合 击中或击不中变换 二值形态学图像处理 噪声滤除 边界提取 对象编著 图像的特性...

AllenMoore
2018/01/27
0
0
[Python图像处理] 九.形态学之图像开运算、闭运算、梯度运算

版权声明:本文为博主原创文章,转载请注明CSDN博客源地址!共同学习,一起进步~ https://blog.csdn.net/Eastmount/article/details/83651172 该系列文章是讲解Python OpenCV图像处理知识,前...

Eastmount
2018/11/02
0
0
(三)matlab数字图像处理实验-图像灰度变换处理

前两节都是熟悉一下怎么在matlab底下对图片做一些操作,并没有什么卵用,这一节稍微有点卵用,灰度变换一般是图像处理的第一步。 数字图像处理实验1-9点击下列链接有源码和链接: matlab数字...

aninstein
2018/01/06
0
0
【工具使用系列】关于 MATLAB 图像处理工具箱,你需要知道的事

如何使用 MATLAB 图像处理工具箱 MATLAB 图像类型转换 图像文件的读写和查询 图像文件的显示 不同类型的图像显示方法 什么是 MATLAB 图像处理工具箱 图像 图像的概念 图像的分类 图像的表示 ...

AllenMoore
2018/01/28
9
0

没有更多内容

加载失败,请刷新页面

加载更多

Linux Wireshark普通用户启动使用方案

当系统安装好Wireshark后请正常启动是否可以进行正常使用,如果不行请参考下列指导 向系统添加一个用户组 sudo groupadd wireshark //如提示此组存在可跳过 将指定用户添加到这个组中 sudo...

CHONGCHEN
今天
1
0
CSS 选择器参考手册

CSS 选择器参考手册 选择器 描述 [attribute] 用于选取带有指定属性的元素。 [attribute=value] 用于选取带有指定属性和值的元素。 [attribute~=value] 用于选取属性值中包含指定词汇的元素。...

Jack088
今天
2
0
数据库篇一

数据库篇 第1章 数据库介绍 1.1 数据库概述  什么是数据库(DB:DataBase) 数据库就是存储数据的仓库,其本质是一个文件系统,数据按照特定的格式将数据存储起来,用户可以对数据库中的数据...

stars永恒
今天
4
0
Intellij IDEA中设置了jsp页面,但是在访问页面时却提示404

在Intellij IDEA中设置了spring boot的jsp页面,但是在访问时,却出现404,Not Found,经过查找资料后解决,步骤如下: 在Run/Debug Configurations面板中设置该程序的Working Directory选项...

uknow8692
昨天
4
0
day24:文档第五行增内容|每月1号压缩/etc/目录|过滤文本重复次数多的10个单词|人员分组|

1、在文本文档1.txt里第五行下面增加如下内容;两个方法; # This is a test file.# Test insert line into this file. 分析:给文档后增加内容,可以用sed 来搞定;也可以用while do done...

芬野de博客
昨天
4
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部