获取shp文件，使用ArcMap进行栅格裁剪与波段合成

因为某种需要，需要处理重庆主城区的遥感影像然后进行监督分类。

什么数据都没有，只能自己下载。

第一步，先获取重庆市主城区的shp数据。

通过阿里的免费开放平台 来下载shp数据。

逐一选择重庆九区，复制JSON API打开后复制内容 或者直接复制“其他类型”，如下图所示。

然后将复制的内容粘贴到mapshaper中打开，这里以北培区为例，下载后导出为shp即可。

第二步，获取重庆市某一时刻的遥感影像

获取遥感影像免费的可以通过地理空间数据云、GEE(Google Earth Engine)等方式来获取，也可以付费购买。

这里我们通过地理空间数据云来下载重庆市的遥感影像，GEE的方式略过，有兴趣自己去了解，GEE可以在线处理影像，不局限于下载影像。

我们进入地理空间数据云，选中高级检索，然后进行各项设置，选择范围、云量等等。

最后点击检索，下载符合条件的影像。

第三步，使用envi对影像进行预处理操作。

这一步具体操作省略，Landsat8不需要几何校正，对影像进行辐射定标和大气校正即可。（不校正对后续处理有影响吗？欢迎评论区交流

第四步，使用Arcmap打开下载好的各区shp，将其合成为一个shp。

选择地理处理中的合并选项，将多个shp合成为1个。

第5步，使用ArcMap进行栅格裁剪。

选择数据管理工具中的栅格->栅格裁剪，

进行各项参数设置，勾选使用输入要素裁剪几何。

 （因为要裁剪的shp影像是多波段的，如果裁剪报错的话，在地理处理中，将并行处理因素设置为0然后保存重复操作即可）

这里我预先通过envi将下载的影像镶嵌到一块了，具体步骤略过了，输入的栅格是镶嵌后的栅格。

裁剪完成，查看最终结果即可。

第6步，波段合成

仍然在栅格处理下，选择波段合成，将栅格合成为你需要的数据格式即可。