Python遥感图像处理应用篇(四)：python如何使用numpy读取遥感图像光谱值

1.数据准备

1.1 影像数据选取

这里选取一景OLI8（Landsat8）数据作为测试数据，该数据已经进行过辐射定标和大气校正。该数据显示如下

数据信息：轨道号12340  时间20190817  测试数据下载见本文结尾处。

由于一景数据过大，处理太耗费时间，这里裁剪一个规则矩形区域作为测试区。

1.2 绘制规则矩形

在ArcCatlog里面新建shapefile图层，投影选择和影像一致（对于较新版本实际上也可以随意设置，裁剪功能会自动转换投影进行裁剪），一般选择一致，这样更好一些。

选择矩形工具进行要素采集，即绘制一个规则矩形，这里是绘制一个水平的规则矩形。具体方法可以参考：Arcgis应用（二）使用editor工具绘制水平矩形

1.3 使用clip工具裁剪影像作为测试数据

使用以下第三个工具