当前位置:文档之家› 多光谱变换

多光谱变换


3-红
TM
1-蓝
4-近红外
Correlation Band 1 Band 2 Band 3 Band 4 Band 1 1.0000 0.9568 0.9380 0.0130
2-绿 Band 2 0.9568 1.0000 0.9793 0.1187
Band 3 0.9380 0.9793 1.0000 0.0803 Band 4 0.0130 0.1187 0.0803 1.0000 Band 5 0.8070 0.8486 0.8807 0.0925 Band 6 0.8454 0.8691 0.8995 -0.0703
Inv
1
4
PC
2
5
3
6
真彩色:TM321 假彩色
TM432
TM741
TM PC123
缨帽变换
缨帽变换
• 缨帽变换是Kauth和Thomas通过分析陆地卫星 MSS图像反映农作物和植被生长过程的多时相 数据,提出的一种经验性的、固定核的多波段 图像的线性组合变换 ,又称K-T变换
缨帽变换要点
y1 11,12,,1n x1
y2
21,22
,,2n
x2
ym
m1
,
m2
,
,mn
xn
多光谱变换
• 主要方法
– 主成分变换,Principal Components Transformation
– 缨帽变换 ,Tasseled Cap Transformation
主成分变换
• 即主成分分析(PCA)、卡夫林-列夫变换(KL变换) , 是根据各波段之间的协方差或相关系数构建的一种 正交线性变换方法
• 变换的目的就是把原来多波段图像中的有用信息集 中到数目尽可能少的新的主成分图像中,并使这些 主成分图像之间互不相关,也就是说各个主成分包 含的信息内容是不重叠的,从而大大减少总的数据 量并使图像信息得到增加
• K-T变换可表示为
uRxc
• 其中x为图像象元矢量,R为变换矩阵,c为避免u中 出现负值而定义的常向量。不同传感器的数据,R不 同
TM数据
0.30370.27930.47430.55850.50820.1863
0.28480.24350.54360.72430.08400.180 0
R 00..1852049200..1098743900..4332972300..3045086000..27011122 00..247567 83
主成分变换的性质和特点
• K-L变换是一种线性变换,而且是当取Y的前p (p<m)个主成分经反变换而恢复的图像X′和原图 像X在均方误差最小意义上的最佳正交变换。它具 有以下性质和特点
– (1)由于是正交线性变换,所以变换前后的方差总和不 变,变换只是把原来的方差(方差大说明信息量大)不 等量的再分配到新的主成分图像中
加法运算
• 指两幅或多幅同样大小的图像对应象元的灰 度值相加
F(x,y)af1(x,y)f2(x,y) fn(x,y)
通a常 1 n
• 可用于削弱图像的加性噪声
差值运算
• 指两幅同样大小的图像对应象元的灰度 值相减
F ( x ,y ) f 1 ( x ,y ) f 2 ( x ,y ) b
若以各特征方量为列构成矩阵,即
U矩阵满足:UTU=UUT=I(单位矩阵),则U矩阵是正交 矩阵。 U矩阵的转置矩阵即为所求的K-L变换的变换矩阵T。
T=
主成分变换则:
y1 11,12,,1m x1

y2
21,22
,,2m
x2
TM5/TM7抑制地形阴影
可以消除或减弱地形阴 影、云影影响和植被干 扰以及显示隐伏构造等
植被指数
• 植被指数(VI, Vegetation Indices)是基
于植被在红光波段(R)的强吸收和近红外波段 (IR)的强反射,通过这两个波段影像的比值或 线性组合来表达植被信息(如叶面指数、植被 覆盖度、生物量等)的参考量 – 差值植被指数 DVIRR – 比值植被指数 RVI IRR – 归一化差值植被指数 ND (IVR R I)(IR R ) – 绿度植被指数:KT变换中的绿度
TM432
差值植被指数 DVI
TM432
比值植被指数 RVI
TM432
归一化差值植被指数 NDVI
植被指数
常用传感器的红外和红光波段
传感器 红外波段 红光波段
TM
4
3
MSS
7(4)
5(2)
QuickBird
4
3
多光谱变换
多光谱变换
• 多光谱遥感影像,特别是陆地卫星的 TM等传感器,波段多,信息量大
多光谱变换
图像运算
图像运算
• 图像运算是指对两幅或两幅以上的输入图像
的对应像元逐个地进行和、差、积、商的四 则运算,以产生有增强效果的图像
• 图像运算是一种比较简单和有效的增强 处理,是遥感图像将增强处理中常用的 一种方法
• 进行代数运算处理的遥感图像必须进行 大气辐射校正、噪音抑制、几何配准等 预处理
• 用于动态监测、运动目标检测和跟踪及目标 识别等
差值运算
• 应用
– 增强在两个不同波段间亮度差异大的目标 – 检测同一区域内目标的运动变化
TM4-TM3
TM41999.9-TM42001.9
比值运算
• 两个波段对应像元的灰度值之比或几个波段 组合的对应像元灰度值之比
– 简单比值运算
F (x ,y ) f2 (x ,y )f1 (x ,y )
3-红
多光谱变换
• 多光谱变换方法可通过函数变换, 达到保留主要信息、降低数据量, 增强有用信息的目的
• 变换的本质是对遥感图像实行线性 变换,使多光谱空间的坐标系按一 定规律进行旋转
多光谱空间
• 多光谱空间就是一个n维坐标系,每一个坐标轴代 表一个波段,坐标值为亮度值,坐标系内的每一个 点代表一个像元
主成分变换
矩阵
其中,m和n分别为波段数(或称变量数)和每幅图像中 的像元数;矩阵中每一行矢量表示一个波段的图像 对于线性变换Y=TX,如果变换矩阵T是正交矩阵,并且 它是由原始图像数据矩阵X的协方差矩阵S的特征向量所 组成,则称此变换为主成分变换
主成分变换的具体过程
1. 根据原始图像数据矩阵X,求出它的协方差矩阵S。
• 像元点在坐标系中的位置可以表示成一个n维向量, 其中每个分量xi表示该点在第i个坐标轴上的投影, 即亮度值
x1
X
x2
x1, x2 ,, xn T
xn
多光谱变换
• 通过线性变换,以获得新的波段数据,实现在降低 数据冗余、压缩数据量的同时保留主要信息、增强 有用信息的目的
• 本质是对多光谱遥感图像实行线性变换,使多光谱 空间的坐标系按一定规律进行旋转
– 组合比值运算
F(x,y)F1(x,y) F2(x,y)
ba00[[gf00((xx,,yy))]] ba11[[gf11((xx,,yy))]] bann[[gfnn((xx,,yy))]]
比值运算
• 应用
– 可以扩大不同地物的光谱差异,区分在单波段 中容易发生混淆的地物
– 可以消除或减弱地形阴影、云影影响和植被干 扰以及显示隐伏构造等
ym
m1
,
m
2
,,mm
xm
式中Y矩阵的行向量
为第j主成分 。
经过K-L变换后,得到一组(m个)新的变量(即Y 的各个行向量),它们依次被称为第一主成分、第 二主成分、…第m主成分。这时若将Y矩阵的各行 恢复为二维图像时,即可以得到m个主成分图像。
主成分变换的应用
• 数据压缩
– 以TM影像为例,共有7个波段处理起来数据量很 大。进行K-L变换后,第一,或前二或前三个主 分量已包含了绝大多数的地物信息,足够分析使 用,同时数据量却大大地减少了。应用中常常只 取前三个主分量作假彩色合成,实现数据压缩, 也可作为分类前的特征选择。
• 数据量太大,在图像处理计算时,也常 常耗费大量的机时、占据大量的磁盘空 间。
• 实际上,一些波段的遥感数据之间都有 不同程度的相关性,存在着数据冗余。
TM
1-蓝
2-绿
Correlation Band 1 Band 2 Band 3 Band 1 1.0000 0.9568 0.9380 Band 2 0.9568 1.0000 0.9793 Band 3 0.9380 0.9793 1.0000 Band 4 0.0130 0.1187 0.0803 Band 5 0.8070 0.8486 0.8807 Band 6 0.8454 0.8691 0.8995
0.28480.24350.5436 0.7243 0.0840 0.1800
0.1509 0.1973 0.3273 0.3406 0.71120.4573
R
0.82420.08490.43920.05800.2012 0.2768
0.32800.05490.1075 0.18550.4357 0.8085
0.32800.05490.10750.18550.43570.8085
0.10840.90220.41200.05730.02510.0238
缨帽变换的物理意义
• MSS数据为4波段数据
x[x1 x2 x3 x4]T
uR xc
• 变换后,
u[u1 u2 u3 u4]T
• u1为“亮度(Brightness)”特征,反映地物总体辐射强度,土 壤亮度;
S是一个实对称矩阵
Xi
xi
1 n
n k1
xik
(即为第i个波段的均值)
主成分变换的具体过程
相关主题