利用Landsat8进行穗帽变换初探

穗帽变换/ K-T变换（Kauth和Thomas最先提出）不是由严格的理论体系推导而产生的，而是基于对Landsat的大量MSS图像统计研究中提出的。穗帽变换是一种特殊的主成分分析，和主成分分析不同的是其转换系数是固定的，因此它独立于单个图像，后来通过研究得出TM和ETM+数据的穗帽变换。

Landsat8 OLI传感器目前还没有相应的穗帽变换模型，参照下表1 中OLI和ETM+波段对照，发现OLI的Band2~7和ETM+的6个多光谱波段范围相近，在没有出现最新OLI穗帽变换模型时，我们是否可以使用ETM+的模型应用于OLI？本文没有探讨这个问题及给出答案，本文假设是可以使用的情况下结束操作流程。

表1： OLI陆地成像仪和ETM+对照表

OLI陆地成像仪			ETM+
序号	波段(μm)	空间分辨率 (m)	序号	波段(μm)	空间分辨率 (m)
1	0.433–0.453	30
2	0.450–0.515	30	1	0.450–0.515	30
3	0.525–0.600	30	2	0.525–0.605	30
4	0.630–0.680	30	3	0.630–0.690	30
5	0.845–0.885	30	4	0.775–0.900	30
6	1.560–1.660	30	5	1.550–1.750	30
7	2.100–2.300	30	7	2.090–2.350	30

    下面在ENVI5.1中，简单介绍下使用ETM+的模型应用于OLI穗帽变换的操作步骤。

（1） 打开一个Landsat8 数据。

（2） 在Toolbox中，选择/Radiometric Correction/Radiometric Calibration，选择LC8多光谱数据。

（3） Calibration Type：Reflectnace，选择输出路径和文件名。单击OK执行定标

图：定标为大气表观反射率

 （4） 在主界面中，选择File->Save as，选择LC8大气表观反射率定标结果，单击Spectral Subset按钮，选择后面6个波段。 选择数据路径和文件名。

（5）          在Toolbox中，选择Transform/Tasseled Cap。在文件对话框中选择前面得到的6个波段大气表观反射率图像文件。

（6）          在Tasseled Cap Transform Parameters面板中，传感器类型选择Landsat 7 ETM。

（7）          选择输出路径和文件名。

图：Tasseled Cap Transform Parameters面板