没有合适的资源?快使用搜索试试~ 我知道了~
实验05 MATLAB图像变换.doc
试读
37页
需积分: 0 0 下载量 62 浏览量
更新于2023-09-19
收藏 1.42MB DOC 举报
实验五主要探讨的是MATLAB中的图像变换,这是图像处理领域中的一个重要概念。图像变换是指通过对图像的像素矩阵进行操作,将其转换成另一种形式,以揭示或提取图像的特定信息,如频率特性,这对于图像的压缩、增强、复原等处理至关重要。
图像变换必须满足几个关键条件。第一,变换应该是可逆的,这意味着通过逆变换,可以将变换后的图像恢复到原始状态,保持原有的信息不变。第二,变换应该有助于后续的图像处理,提供更易于分析的图像特征或简化计算步骤。第三,变换的算法应尽可能简单,最好存在快速算法,因为图像处理通常涉及大量像素的计算,速度快的算法能显著提高效率。尤其是对于方形图像且尺寸为2的幂次方的图像,往往可以采用快速算法,如傅立叶变换中的快速傅立叶变换(FFT)。
图像变换的数学表达通常采用二维离散变换的形式,包括正变换和逆变换。正变换将原始图像(f)转换为频域表示(F),而逆变换则将频域表示恢复为图像表示。具体公式如下:
1. 正变换:\( F(u, v) = \sum_{x=0}^{N_x-1}\sum_{y=0}^{N_y-1} f(x, y)g(u - x, v - y) \)
2. 逆变换:\( f(x, y) = \frac{1}{N_xN_y}\sum_{u=0}^{N_u-1}\sum_{v=0}^{N_v-1} F(u, v)h(u - x, v - y) \)
其中,\( g \) 和 \( h \) 分别是正变换核和逆变换核,\( N_x \) 和 \( N_y \) 是图像的宽度和高度。
当变换核满足特定条件时,例如 \( g(u - x, v - y) = g(u)g(v) \) 和 \( h(u - x, v - y) = h(u)h(v) \),则称变换核是可分离的,这意味着二维变换可以分解为两个一维变换,先对行进行变换,然后对列进行变换,这样可以大大提高计算效率。
在MATLAB中,可以使用内置函数实现这些图像变换,例如`fft2`用于二维傅立叶变换,`ifft2`用于二维逆傅立叶变换。理解并掌握这些基本的图像变换及其性质,对于在MATLAB中进行图像处理和分析是十分必要的。通过实践实验,学生能够更好地理解和应用这些理论知识,提高图像处理技能。
实验五
实验五
图像变换
图像变换
一、概述
一、概述
1
1
、变换条件
、变换条件
图
图
像
像
可
可
以
以
看
看
作
作
是
是
一
一
个
个
矩
矩
阵
阵
,
,
所
所
谓
谓
图
图
像
像
变
变
换
换
,
,
就
就
是
是
通
通
过
过
变
变
换
换
矩
矩
阵
阵
,
,
将
将
图
图
像矩阵变换成另一个矩阵。变换后的矩阵能得到某些图像的信息。
像矩阵变换成另一个矩阵。变换后的矩阵能得到某些图像的信息。
通
通
常
常
,
,
变
变
换
换
后
后
的
的
图
图
像
像
能
能
体
体
现
现
图
图
像
像
的
的
频
频
率
率
特
特
征
征
,
,
可
可
以
以
用
用
于
于
图
图
像
像
的
的
数
数
据
据
压
压
缩和各种处理。
缩和各种处理。
图像变换必需满足一下三个条件:
图像变换必需满足一下三个条件:
(1)
(1)
变换是可逆的。
变换是可逆的。
变
变
换
换
后
后
的
的
图
图
像
像
能
能
保
保
持
持
原
原
始
始
图
图
像
像
的
的
信
信
息
息
,
,
可
可
以
以
通
通
过
过
逆
逆
变
变
换
换
矩
矩
阵
阵
把
把
图
图
像
像
真
真
实复原。
实复原。
(2)
(2)
变换后能给图像的进一步运算带来方便。
变换后能给图像的进一步运算带来方便。
也
也
就
就
是
是
说
说
,
,
图
图
像
像
的
的
变
变
换
换
具
具
有
有
一
一
定
定
的
的
含
含
义
义
,
,
变
变
换
换
后
后
的
的
图
图
像
像
要
要
么
么
体
体
现
现
图
图
像
像
的某些特征,要么在数据上带来某些方便的处理。
的某些特征,要么在数据上带来某些方便的处理。
(3)
(3)
变换的算法简单,最好有快速算法。
变换的算法简单,最好有快速算法。
图
图
像
像
的
的
变
变
换
换
通
通
常
常
要
要
经
经
过
过
两
两
次
次
矩
矩
阵
阵
乘
乘
法
法
的
的
运
运
算
算
,
,
运
运
算
算
的
的
速
速
度
度
关
关
系
系
到
到
图
图
像
像
变
变
换
换
的
的
好
好
坏
坏
。
。
大
大
多
多
数
数
图
图
像
像
变
变
换
换
,
,
要
要
求
求
图
图
像
像
是
是
方
方
阵
阵
,
,
且
且
行
行
列
列
数
数
是
是
2
2
的
的
幂
幂
次
次
方
方
才有快速算法。
才有快速算法。
2
2
、变换公式
、变换公式
图像变换都是二维的离散变换,通用公式由下列两式给出:
图像变换都是二维的离散变换,通用公式由下列两式给出:
正变换:
正变换:
� � � � � �
� �
�
�
�
�
�
1
0
1
0
,,,,,
N
x
N
y
vuyxgyxfvuF
;
;
逆变换:
逆变换:
� � � � � �
� �
�
�
�
�
�
1
0
1
0
,,,,,
N
u
N
v
vuyxhvuFyxf
。
。
其中:
其中:
1,...,2,1,0,1,...2,1,0, ���� NyxNvu ;
,
,
� �
vuyxg ,,,
-称为正变换核,
-称为正变换核,
� �
vuyxh ,,,
-称为逆变换核。
-称为逆变换核。
假如:
假如:
� � � � � �
vyguxgvuyxg ,,,,,
21
�
,
,
� � � � � �
vyhuxhvuyxh ,,,,,
21
�
则称变换核是可分离的。如果:
则称变换核是可分离的。如果:
� � � �
vyguxg ,,
21
�
,
,
� � � �
vyhuxh ,,
22
�
则称变换是加法对称的,这时,正变换可表示为:
则称变换是加法对称的,这时,正变换可表示为:
� � � � � � � �
� � � � � � � �
� �
� �
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
1
0
1
0
1
0
1
0
,,,,
,,,,
N
x
N
y
N
x
N
y
uxgvygyxfvuF
vyguxgyxfvuF
即二维的变换可分离成两次的一维变换:
即二维的变换可分离成两次的一维变换:
先对行进行变换:
先对行进行变换:
� � � � � �
�
�
�
�
1
0
,,,
N
y
vygyxfvxF
,
,
再对列进行变换:
再对列进行变换:
� � � � � �
�
�
�
�
1
0
,,,
N
x
uxgvxFvuF
。
。
用矩阵来表示:
用矩阵来表示:
y
x
u
v
y
x
x
v
u
v
� �
yxf ,
可表示成
可表示成
� � � � � �
� � � � � �
� � � � � �
�
�
�
�
�
�
�
�
�
�
�
�
����
�
�
�
1,1...1,10,1
............
1,1...1,10,1
1,0...1,00,0
NNfNfNf
Nfff
Nfff
f
� �
vyg ,
可表示成
可表示成
� � � � � �
� � � � � �
� � � � � �
�
�
�
�
�
�
�
�
�
�
�
�
����
�
�
�
1,1...1,10,1
............
1,1...1,10,1
1,0...1,00,0
NNgNgNg
Nggg
Nggg
G
� �
vxF ,
可表示成
可表示成
� � � � � �
� � � � � �
� � � � � �
�
�
�
�
�
�
�
�
�
�
�
�
����
�
�
�
1,1...1,10,1
............
1,1...1,10,1
1,0...1,00,0
111
111
111
1
NNFNFNF
NFFF
NFFF
F
显然有:
显然有:
� � � � � �
�
�
�
�
1
0
,,,
N
y
vygyxfvxF
可表示成:
可表示成:
fGF �
1
若
若
� �
uxg ,
可表示成
可表示成
� � � � � �
� � � � � �
� � � � � �
�
�
�
�
�
�
�
�
�
�
�
�
����
�
�
�
1,1...1,10,1
.............
1,1...1,10,1
1,0...1,00,0
NNgNgNg
Nggg
Nggg
G
� �
vuF ,
可表示成
可表示成
� � � � � �
� � � � � �
� � � � � �
�
�
�
�
�
�
�
�
�
�
�
�
����
�
�
�
1,1...1,10,1
............
1,1...1,10,1
1,0...1,00,0
NNFNFNF
NFFF
NFFF
F
显然有:
显然有:
� � � � � �
�
�
�
�
1
0
,,,
N
x
uxgvxFvuF
可表示成:
可表示成:
1
FGF
T
�
则有:
则有:
fGGF
T
�
同样反变换有:
同样反变换有:
FHHf
T
�
显然,若
显然,若
IHG �
,则:
,则:
ffIIGHfGHHfGGHFHHf
TTTTT
����� )()()(
即,要使变换可逆,必需要求变换是正交变换。
即,要使变换可逆,必需要求变换是正交变换。
所
所
以
以
,
,
图
图
像
像
变
变
换
换
,
,
就
就
是
是
找
找
出
出
这
这
样
样
的
的
变
变
换
换
矩
矩
阵
阵
,
,
产
产
生
生
正
正
交
交
变
变
换
换
。
。
下
下
面
面
介
介
绍傅立叶变换和余弦变换。
绍傅立叶变换和余弦变换。
二、图像的傅立叶变换
二、图像的傅立叶变换
1
1
、变换表达式
、变换表达式
正变换:
正变换:
� � � � � �
� �
�
�
�
�
�
�
�
�
�
�
���
1
0
1
0
2
exp,,
N
x
N
y
vyux
N
j
yxfvuF
�
逆变换:
逆变换:
� � � � � �
� �
�
�
�
�
�
�
�
�
�
�
��
1
0
1
0
2
2
exp,
1
,
N
u
N
v
vyux
N
j
vuF
N
yxf
�
其中:
其中:
1,...,2,1,0,,, �� Nvuyx
若令
若令
�
�
�
�
�
�
��
N
j
w
N
�
2
exp
,显然傅立叶变换的正变换核有:
,显然傅立叶变换的正变换核有:
� �
vy
N
ux
N
wwvyguxg ��),(,
,正变换矩阵为:
,正变换矩阵为:
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�������
����
����
)1()1(1)1(0)1(
)1(11101
)1(01000
...
............
...
...
NN
N
N
N
N
N
N
NNN
N
NNN
www
www
www
G
逆变换核
逆变换核
yv
N
xu
N
ww
N
vyhuxh
��
�
2
1
),(),(
,逆变换矩阵为:
,逆变换矩阵为:
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
����������
�������
�������
)1()1(1)1(0)1(
)1(11101
)1(01000
...
............
...
...
1
NN
N
N
N
N
N
N
NNN
N
NNN
www
www
www
N
H
。
。
下面例子,
下面例子,
以
以
N=4
N=4
为例:
为例:
【例】变换核
【例】变换核
w4
w4
=
=
inline('exp(-2*j*pi/4*x*u)','x','u');
inline('exp(-2*j*pi/4*x*u)','x','u');
%inline
%inline
定义一个函数
定义一个函数
G=[w4(0,0),w4(0,1),w4(0,2),w4(0,3);...
G=[w4(0,0),w4(0,1),w4(0,2),w4(0,3);...
w4(1,0),w4(1,1),w4(1,2),w4(1,3);...
w4(1,0),w4(1,1),w4(1,2),w4(1,3);...
w4(2,0),w4(2,1),w4(2,2),w4(2,3);...
w4(2,0),w4(2,1),w4(2,2),w4(2,3);...
w4(3,0),w4(3,1),w4(3,2),w4(3,3)];
w4(3,0),w4(3,1),w4(3,2),w4(3,3)];
H=[w4(-0,0),w4(-0,1),w4(-0,2),w4(-0,3);...
H=[w4(-0,0),w4(-0,1),w4(-0,2),w4(-0,3);...
w4(-1,0),w4(-1,1),w4(-1,2),w4(-1,3);...
w4(-1,0),w4(-1,1),w4(-1,2),w4(-1,3);...
w4(-2,0),w4(-2,1),w4(-2,2),w4(-2,3);...
w4(-2,0),w4(-2,1),w4(-2,2),w4(-2,3);...
w4(-3,0),w4(-3,1),w4(-3,2),w4(-3,3)];
w4(-3,0),w4(-3,1),w4(-3,2),w4(-3,3)];
H=H/4;
H=H/4;
G,H,G*H
G,H,G*H
G
G
=
=
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
1.0000
0.0000
0.0000
-
-
1.0000i
1.0000i
-1.0000
-1.0000
-
-
0.0000i
0.0000i
-0.0000
-0.0000
+
+
1.0000i
1.0000i
1.0000
1.0000
-1.0000
-1.0000
-
-
0.0000i
0.0000i
1.0000
1.0000
+
+
0.0000i
0.0000i
-1.0000
-1.0000
-
-
0.0000i
0.0000i
1.0000
1.0000
-0.0000
-0.0000
+
+
1.0000i
1.0000i
-1.0000
-1.0000
-
-
0.0000i
0.0000i
0.0000
0.0000
-
-
1.0000i
1.0000i
H
H
=
=
0.2500
0.2500
0.2500
0.2500
0.2500
0.2500
0.2500
0.2500
0.2500
0.2500
0.0000
0.0000
+
+
0.2500i
0.2500i
-0.2500
-0.2500
+
+
0.0000i
0.0000i
-0.0000
-0.0000
-
-
0.2500i
0.2500i
0.2500
0.2500
-0.2500
-0.2500
+
+
0.0000i
0.0000i
0.2500
0.2500
-
-
0.0000i
0.0000i
-0.2500
-0.2500
+
+
0.0000i
0.0000i
0.2500
0.2500
-0.0000
-0.0000
-
-
0.2500i
0.2500i
-0.2500
-0.2500
+
+
0.0000i
0.0000i
0.0000
0.0000
+
+
0.2500i
0.2500i
ans
ans
=
=
1.0000
1.0000
-0.0000
-0.0000
+
+
0.0000i
0.0000i
0
0
+
+
0.0000i
0.0000i
0.0000
0.0000
+
+
0.0000i
0.0000i
-0.0000
-0.0000
-
-
0.0000i
0.0000i
1.0000
1.0000
-0.0000
-0.0000
0
0
+
+
0.0000i
0.0000i
0
0
-
-
0.0000i
0.0000i
-0.0000
-0.0000
-
-
0.0000i
0.0000i
1.0000
1.0000
+
+
0.0000i
0.0000i
-0.0000
-0.0000
+
+
0.0000i
0.0000i
0.0000
0.0000
-
-
0.0000i
0.0000i
0
0
-
-
0.0000i
0.0000i
-0.0000
-0.0000
-
-
0.0000i
0.0000i
1.0000
1.0000
【例】傅立叶变换
【例】傅立叶变换
w4
w4
=
=
inline('exp(-2*j*pi/4*x*u)','x','u');
inline('exp(-2*j*pi/4*x*u)','x','u');
G=[w4(0,0),w4(0,1),w4(0,2),w4(0,3);...
G=[w4(0,0),w4(0,1),w4(0,2),w4(0,3);...
w4(1,0),w4(1,1),w4(1,2),w4(1,3);...
w4(1,0),w4(1,1),w4(1,2),w4(1,3);...
w4(2,0),w4(2,1),w4(2,2),w4(2,3);...
w4(2,0),w4(2,1),w4(2,2),w4(2,3);...
w4(3,0),w4(3,1),w4(3,2),w4(3,3)];
w4(3,0),w4(3,1),w4(3,2),w4(3,3)];
H=[w4(-0,0),w4(-0,1),w4(-0,2),w4(-0,3);...
H=[w4(-0,0),w4(-0,1),w4(-0,2),w4(-0,3);...
w4(-1,0),w4(-1,1),w4(-1,2),w4(-1,3);...
w4(-1,0),w4(-1,1),w4(-1,2),w4(-1,3);...
w4(-2,0),w4(-2,1),w4(-2,2),w4(-2,3);...
w4(-2,0),w4(-2,1),w4(-2,2),w4(-2,3);...
w4(-3,0),w4(-3,1),w4(-3,2),w4(-3,3)];
w4(-3,0),w4(-3,1),w4(-3,2),w4(-3,3)];
H=H/4;
H=H/4;
f=[1
f=[1
2
2
3
3
4;5
4;5
6
6
7
7
8;9
8;9
10
10
11
11
12;13
12;13
14
14
15
15
16];
16];
F=G'*f*G;
F=G'*f*G;
f1=H'*F*H;
f1=H'*F*H;
f,F,f1
f,F,f1
f
f
=
=
1
1
2
2
3
3
4
4
5
5
6
6
7
7
8
8
9
9
10
10
11
11
12
12
13
13
14
14
15
15
16
16
F
F
=
=
1.0e+002
1.0e+002
*
*
1.3600
1.3600
-0.0800
-0.0800
+
+
0.0800i
0.0800i
-0.0800
-0.0800
-
-
0.0000i
0.0000i
-0.0800
-0.0800
-
-
0.0800i
0.0800i
-0.3200
-0.3200
-
-
0.3200i
0.3200i
0.0000
0.0000
+
+
0.0000i
0.0000i
-0.0000
-0.0000
+
+
0.0000i
0.0000i
-0.0000
-0.0000
-
-
0.0000i
0.0000i
-0.3200
-0.3200
+
+
0.0000i
0.0000i
0.0000
0.0000
+
+
0.0000i
0.0000i
0.0000
0.0000
+
+
0.0000i
0.0000i
-0.0000
-0.0000
+
+
0.0000i
0.0000i
-0.3200
-0.3200
+
+
0.3200i
0.3200i
0.0000
0.0000
+
+
0.0000i
0.0000i
0.0000
0.0000
+
+
0.0000i
0.0000i
0.0000
0.0000
+
+
0.0000i
0.0000i
f1
f1
=
=
1.0000
1.0000
+
+
0.0000i
0.0000i
2.0000
2.0000
+
+
0.0000i
0.0000i
3.0000
3.0000
+
+
0.0000i
0.0000i
4.0000
4.0000
+
+
0.0000i
0.0000i
5.0000
5.0000
-
-
0.0000i
0.0000i
6.0000
6.0000
+
+
0.0000i
0.0000i
7.0000
7.0000
+
+
0.0000i
0.0000i
8.0000
8.0000
+
+
0.0000i
0.0000i
9.0000
9.0000
-
-
0.0000i
0.0000i
10.0000
10.0000
-
-
0.0000i
0.0000i
11.0000
11.0000
+
+
0.0000i
0.0000i
12.0000
12.0000
+
+
0.0000i
0.0000i
剩余36页未读,继续阅读
资源推荐
资源评论
2024-04-19 上传
193 浏览量
2021-09-14 上传
2021-09-14 上传
2022-07-05 上传
5星 · 资源好评率100%
145 浏览量
189 浏览量
2022-07-07 上传
5星 · 资源好评率100%
200 浏览量
117 浏览量
151 浏览量
2022-07-05 上传
2024-04-20 上传
2021-10-31 上传
181 浏览量
2020-07-12 上传
5星 · 资源好评率100%
118 浏览量
资源评论
BTCKing
- 粉丝: 712
- 资源: 35
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- C#上位机开发源码 上位机项目源代码 采用基于RS485通讯总线的ModbusRtu协议,支持用户权限管理、sqlite数据库、实时曲线、历史曲线、历史报表、导出Excel、主界面布局可调带记忆等功能
- 基于plc的污水处理,组态王动画仿真,带PLC源代码,组态王源代码,图纸,IO地址分配
- Rhino(犀牛)插件ladybug-tools-1-8-0
- 三相10Kw光伏并网逆变器 包含全套理图 PCB 源代码
- MATLAB代码:考虑P2G和碳捕集设备的热电联供综合能源系统优化调度模型 关键词:碳捕集 综合能源系统 电转气P2G 热电联产 低碳调度 参考文档:Modeling and Optimiza
- 永磁同步直线电机仿真实例,仿真教学 maxwell16.0版本 12槽11极 包括图中模型以及一个仿真设置要点word文档教程
- 基于mpx+vue+node.js的双端网盘系统的设计与实现源代码全套技术资料.zip
- welearn刷时长版本v3.0.bat
- 前端分析-2023071100789-y5
- 前端分析-2023071100789
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功