实验内容:已知某蚀变矿物在Landsat 7 ETM传感器中的 第三波段表现为强反射,第一波段表现为强吸收,请在图像上,增强该矿物的信息,并进行密度分割。同理 提取羟基信息
(在第五波段表现为强反射 第七波段表现为强吸收)。
步骤:1.先进行波段运算,用b3/b1,值大就表明信息强。点击主菜单点击basic tools>band math 输入运算公式,点击OK出现如图1-1所示对话框,定义如图。
然后进行假彩色密度分割,得到如图1-2所示图像。经过处理白色部分是矿物分布的区域。如图所示。
2.进行主成分分析,点击主菜单transform>principal components>forward pc rotation>compute new statistics and rotate选择1345波段进行主成分分析。出现如图2-1所示对话框,进行设置。并查看统计数据如图2-2所示.
在统计表中的Eigenvector(特征向量表明了各波段对图像贡献的大小)选项中我们观察band3一栏中band2为正band1为负。这符合我们的要求,强反射正值越大,强吸收负值数值越大,效果越好。那么就选择band3这一行即波段pc3,对它进行统计分析,点击主菜单basic tools>statistic>compute statistic,找出标准差。如图2-3所示。标准差为7.992046(stdev)。然后对pc3波段进行密度分割,最小值为2倍标准差,即16,最大值保留原有的值。如图2-4所示。为输出后图像与密度分割对话框。
图2-4红色区域就是蚀变矿物信息的分布
3.提取羟基信息
步骤:同上。先进行波段运算b5/b7,得到如图3-1所示图像
然后进行主成分分析用波段1457,统计结果如下,图3-2所示。
我们可以看到PC4符合要求,对其进行统计分析,得到标准差为3.5226如图3-3所示。
最后对PC4进行密度分割,最小值为标准差2倍,得到图像如图3-4所示。下载本文