第十讲高光谱遥感图象混合象元分析.docx
- 文档编号:12177898
- 上传时间:2023-04-17
- 格式:DOCX
- 页数:13
- 大小:143.11KB
第十讲高光谱遥感图象混合象元分析.docx
《第十讲高光谱遥感图象混合象元分析.docx》由会员分享,可在线阅读,更多相关《第十讲高光谱遥感图象混合象元分析.docx(13页珍藏版)》请在冰豆网上搜索。
第十讲高光谱遥感图象混合象元分析
第十讲高光谱遥感图象混合象元分析
混合象元的概念:
遥感器所获取的地面反射或发射光谱信号是以象元为单位记录的。
它是象元
所对应的地表物质光谱信号的综合。
图象中每个象元所对应的地表,往往包含不
同的覆盖类型,他们有着不同的光谱响应特征。
而每个象元则仅用一个信号记录这些“异质”成分。
若该象元仅包含一种类型,则为纯象元(purepixel),也称为
端元(endmembe),它所记录的正是该类型的光谱响应特征或光谱信号;若该象元包含不止一种土地覆盖类型,则成为混合象元(mixedpixel),它记录的是所
对应的不同土地覆盖类型光谱响应特征的综合。
由于传感器的空间分辨力限制以
及自然界地物的复杂多样性,混合像元普遍存在于遥感图象中。
混合象元模型
光谱混合形式上可以分为致密式(intrinsic)、聚合式(aggregate和整合式(area)三种情形(如图),本质上分可以分为线性混合和非线性混合两种模式。
线性混合模型假定到达传感器的光子只与一种物质发生作用[rast,1991];当混合
元素尺寸小,入射光子与多于一种以上的物质发生作用时,导致非线性混合[smiths,1985;Mustard,1987]。
IntrinsicAggregateAreal
i.线性混合模型
通常情况下,高光谱图象中每个象元都可以近似认为是图象中各个端元的线
性混合象元:
N
p八qein=Ecn
(1)
im
"Ci=1⑵
i=1
Q 其中N为端元数,p为图象中任意一L维光谱向量(L为图象波段数),E=[2…eN]为LN矩阵,其中的每列均为端元向量。 c二(cc…Cn亍为系数向量,Ci表示象元p中端元ei所占的比例,n为误差项。 在误差项n很小的情况下,满足 (1)、 (2)和(3)的所有点的集合正好构成一个高维空间的凸集,这些端元则坐落于这个凸面单形体的顶点。 以两个波段 三个端元为例来说明它们之间的几何关系(图2).从图2可以看出,端元a,b,c分别位于三角形体的顶点,三角形内部的点则对应着图象中的混合象元•这样,提取高光谱图象的端元问题就转化为求单形体的顶点的问题. 端兀a端兀b 波段i 图2两个波段三个端元的散点图在空间上具有明显的三角形结构 2.非线性混合模型 三.端元提取 1.PPI 当把特征空间中的所有散点往一个单位向量U上投影时,端元就会投影到u的两侧,而混合象元则会投影到中部。 基于这个思想,可以让图象在n个随机的单位向量上投影,并且记下每个象元被投影到端点的次数,即为纯象元指数(PPI). 当然,被投影到随机向量端点的次数越多,说明此象元为纯粹象元的可能性越大。 2.N-FINDR N-FINDR算法正是通过求最大单形体的体积而得到各个端元的,其体积公式 如下: 1 (n-1)! V(E)-abs(E)(5)其中ei为表征第i个端元的列向量,V是由ene? …,e”这 N个端元所构成的单形体的体积,为行列式运算符.由于用到了求行列式的运算,所以要求E必须为方阵,这样向量ei的维数必须为N-1,但原始的高光谱数据往往是不满足这个条件的,于是需要先对原始数据进行降维处理,这也正是 N-FINDR算法可能引起偏差(比如’忽视'小目标)的原因所在,同时也是此算法的不足之处 3.距离法 首先,以图象中所有象元的平均向量eo为初始值;从图象中找出距离e0最远的点即为第一个端元e;距离e最远的点e? 即为第二个端元,同时记下此距离为d (1);然后找出距离e1与e2所构成直线最远的点e3,即与这两个顶点端元围成的三角形面积最大的点,此即为第三个端元,同时记下此距离为d (2);然后,再找出距这三个端元0、e? 和e3所构成的平面最远的点e4,即与这三个端元围成的三棱锥体积最大的点,同时记下此距离为d(3)。 依此类推,可以找到图象中的所有端兀e1,e2,…,eN(其中ei=(e耳,…,eu),L为波段数), 上述算法的关键在于如何得到高维空间中一个点到一个超平面的距离。 假设现在已经得到端元e1,e2,…,el,p是图象中的任一象元(光谱向量),我们通过施密特正交化得到p到e1,e2,…,ei所构成的超平面的单位垂直向量u,过点p以u为方向可以得到L维空间的一个直线方程;过e1,e2/',ei中的任何一点以u为法方向可以得到L维空间的一个超平面方程;上述直线方程与超平面方程在L维空间 有唯一的交点,可以验证,p到此交点的距离既是p到e1,e2,…,ei所构成的超平面的距离,其中距离ei,e2,…,ei所构成的超平面最大的象元即为第I1个端元ei1,下面给出具体的数学描述: 在L维特征空间中,过点p方向为u的直线方程为: y-p二ku(6) uIy-ei)=0(7) 过点ei,以u为法向量的超平面方程为: 由(6)、(7)可以得到,两个方程的交点为: (8 于是,L维特征空间中,过点p方向为u的直线与以过任何一个端元edi=1,…,丨), 以u为法向量的超平面均交于相同的一点y0。 并且点p到点y0的距离即为点p到ei,e2,…,ei所构成的超平面的距离 4.体积法 蔡聪明把低维空间平行多面体的体积成功得推广到了高维,得到了m维平行 多面体的体积为: Vm=V(vi,v2,…,vm)二det(AmA)=、G(vi,v2,,vm)(9) 其中vi(i=1,2,…,m)为m维空间中的列向量,Am=(v1,v2,…,vm), (a)在m=1时,即得到一维平行多面体的体积就是线段的长度: V(v1)的体积-Vi的长度二Vi二v1v1 (b)m=2时,可得到二维平行多面体即平行四边形的面积: V(Vi,v2)的体积二由向量vi与v2所决定的平行四边形的面积 (c)m=3时,即可得到三维平行多面体即平行六面体的体积: V(vi,v2,v3)的体积=由向量v1,v2和v3所决定的平行六面体的体积 -G(v1,v2,v3)- 由于m维单形体包含于m维平行多面体之中,并且它们的体积有如下对应关 系: 1- VmVm(11) m! 其中Vm为m维单形体体积,V;为m维平行多面体体积。 显然,当m=1时,一维单形体与平行多面体均为直线,两者体积相同;当m=2时,二维单形体为三角 形,其体积为相应平行四边行体积的一半;当m=3时,三维单形体为三棱锥, 其体积为相应的平行六面体体积的1。 可以证明,在高维空间中同样有上述关系 6 成立。 上述体积公式中单形体均由向量组成,对于高光谱图象,象元在其特征空间中都是以散点的形式存在的,图象中的每个象元对应于特征空间中的一个点,假设P1,P2,…,Pn是图象中的门个象元(光谱向量),令 An*=(Pn-P1,Pn-卩2,…,Pn~卩2),则在特征空间中,以这n个点为顶点的门-1维单形体体积为: Vn=V(Pn-P1,Pn-P2,…,Pn-Pnd)=“I«AnJAn斗(12) (n-1)! Vl 其中•为行列式运算符,由于AiTAi(i=1,…,n-1)—定为方阵,所以式(12)对于任何维数的高光谱数据都是成立的,因而可直接把本公式应用于原始高光谱数据的 端元提取,而不需要对原始数据进行降维处理,这样就避免了N-FINDR算法中因数据降维而可能带来的偏差。 当n=2时,式(12)给出了象素》与p2之间的距 离;当n=3时,A2A2为由象素pi,p2和p3为顶点确定的平行四边形的面积, 再乘以系数-即是以象素pi,p2和p3为顶点的三角形的面积;当n=4时,.A3A32 1 为由象素Pi,p2,p3和p4为顶点确定的平行六面体的体积,再乘以系数即是以象 6 —即是以象(n-1)! 素p1,P2,p3和p4为顶点的三棱锥的体积;当n>4时,A爲AnJ为由象素 p1,p2,…,pn为顶点确定的高维平行多面体的体积,再乘以系数 素p1,p2,…,pn为顶点的高维凸面单行体的体积。 价。 下面给出证明: 由于: T p1 T p2一p1 pn一p1 =sqrt( L(pN-pi) T p2一p1pN T p2-P1p2p1 二sqrt( T pNp1p2-p1 =y/A pN-pj=Vn 因而: V(E"^abs(E)"(P2一p1^彳一p1, 目标提取就是根据目标的特性(比如吸收特征,光谱曲线等)在抑制背景的同时 从图象中得到该目标的分布情况。 我们的目标提取算法基于以下事实: 在高维特征空间中,每一个端元都游离于其它所有的端元构成的超平面之外,且是距离超 平面最远的点。 这样,对于任何一个端元,我们都可以得到一个最佳的投影方向(即此端元到其它端元所构成的超平面的垂线方向),在这个方向上投影,将得到此端元和别的端元的最佳分离效果,相应地,我们可以得到各种地物的成分图。 下面以两个波段三个端元时的情况为例来说明我们的算法,如下图3,端元A,B,C 分别位于三角形的三个顶点,点D,E,F是它们在各自对应边上的垂足•于是我们就可以得到三个单位矢量: 其中,直线段AD,BE,CF这里均看作是矢量。 这三个单位矢量就是我们所要得到端元投影向量,在其上投影就可以得到地物与背景的最佳分离效果。 比如,在图 3中,在直线AD方向上的投影将使得端元A所对应的地物与端元B和端元C所组成的背景得以最大的区分,而且端元B与端元C所对应的地物在此方向上 没有任何区别。 ih 波段i 个端元到另外所有端元所构成的超平面的垂线方向 我们利用施密特正交化很好 图3端元投影向量示意图 地解决了这一问题。 假设(e「e2,…,eN)是从图象中求出的所有端元,下面我们给 出从图象中提取出端元eN所对应目标的具体算法(对其它端元同 理): =e? _ei 1 b2二e3(eie? ) 2 =b2-Ci CiCi cbbNACN_2c CNi=bN」_CN_2_ Ci NJ-Ci CN_2CN_2 CNJ -JCnaCni 然后把图象中的所有象元投影到IN就得到端元eN所对应的地物在整个图象中的 分布情况。 2.最小二乘 线性混合模型一般可分为三种情形: 公式(i)为第一种情形,为无约束的线性混合模型,加上如下约束条件 (2)则为部分约束混合模型,再加上约束条件(3)则为全约束混合模型。 线性解混就是在已知所有端元的情况下求出它们图象的各个象元中所占的比例,从而得到与反应每个端元在图象中分布情况的比例系数图。 利用最小二乘法可以得到方程(i)的无约束解: (i5) 再加上 (2)可以得到部分约束的最小二乘解: 其中I为N阶单位矩阵,I为分量均为i的N维列向量 3.OSP 4.几何法 以两个波段、三个端元为例来阐明我们算法的原理,如图4,像元P是以端元 A,B,C为顶点的三角行内部的一个点,则此像元中端元A,B,C对应的地物的含量分别为: SpBC c1: Sabc A 波段i 图4二维情况下混合像元中各端元比例的几何示意图 其中Sabc为三角形ABC的面积,Spbc,Spac和Spab也分别是相应三角形的面积。 我们把上述结论推广到了高维空间,并且证明了在高维空间中对于凸面单形体仍有上述规律成立• 五.混合象元分析中代数与几何的关系 图像空间 特征空间
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第十 光谱 遥感 图象 混合 分析