IDL代码

具体的代码下载链接如下：
LT4/5
LE7_SLC_off
LE7_SLC_on
LC8

使用方法

文件夹准备

首先，你需要这样一个文件夹目录结构

其中results文件夹是批处理的输出文件夹

results下面那些文件夹分别放解压后的遥感资料，如果你只需要处理LT5，你就只建一个LT4_5文件夹，如果只要处理LC8就只建一个LC8文件夹，另外两个同理。

红框里面的文件夹，除了母文件夹concave_convex可以按你的需要重命名，其他文件夹都要按图中的名字来（要么你可以自己改脚本中的路径）。

数据准备

需要准备的数据有landsat遥感资料、对应的掩膜shpfile

遥感影像资料准备

LT4/5、lc8、LE7_SCL_on的影像资料直接解压到相应文件夹，然后直接运行对应脚本即可。

对于LE7_SCL_off，即LE7在2003年之后有条带的影像资料，需要先在ENVI平台里面对多光谱数据集和全色数据集gapfill去条带，在对应文件夹下生成*_Multi和*_Pan文件：

掩膜shapefile准备

定义投影

你必须要把shp面要素定义投影（可以在arcmap里面做）为对应遥感影像的投影坐标系，不然掩膜会覆盖不到你的遥感影像。如我这里的投影坐标系是WGS84-UTM-49N,那我就必须把123040.shp这个文件定义投影为WGS84-UTM-49N。

掩膜文件名

shp文件的名字也要和遥感影像的条带号和行列号对应，这样脚本便可以识别文件夹下还没处理的对应的遥感影像。比如，我想要处理123040和124039两个区域的水体，但是它们都在同一个文件夹下，需要通过掩膜的名字来确定要提取的是对应的区域：

我这里用了123040和124039两个掩膜

示例

首先准备解压遥感数据到对应文件夹：

点击运行LC8.pro之后，首先会让你选择掩膜
这里要处理的数据是123040的，故选择掩膜123040.shp（掩膜的位置任意，只要你找得到）：

选择掩膜之后，会让你选择你要输出的水体shp文件的文件夹（同样的随便你放哪）：

然后等程序运行即可，运行后会在results下产生一个文件夹：

里面就是水体：

因为考虑到需要在arcmap里面打开这些水体，故脚本同时也把这些批处理的水体保存在了另外一个文件夹下，即我们之前选择的那个输出水体的文件夹：

脚本没有覆写功能，所以当你需要重新处理一个遥感，你需要在results里面把对应的文件夹删去。

下一步

我现在做遥感处理选择的方法是gee+python，不再写idl和arcpy的脚本了，参见我主页置顶文章

【IDL】使用IDL对Landsat系列进行批量水体提取相关推荐

Landsat系列卫星地表温度批量反演代码（大气校正法）
针对Landsat5.Landsat7.Landsat8的热红外波段反演地表温度的代码,可以批量进行温度的反演,但需要有前期的一些准备,包括大气校正参数的获取.可见光波段的大气校正等,以及文件夹的准备 ...
Landsat系列卫星介绍及影像下载
1.LandSat系列卫星介绍: 1.Landsat系列卫星概述: 美国NASA的陆地卫星(Landsat)计划从1972年7月23日以来,已发射8颗(第6颗发射失败).目前Landsat1-4均相继 ...
利用USGS网站下载Landsat系列影像（个人总结的小技巧，适用新手小白）
从USGS网站下载Landsat系列影像 USGS(United States Geological Survey, 美国地质勘探局)是美国内政部所属的科学研究机构,其主要负责自然灾害.地质.资源.地 ...
Landsat系列数据、GF1号数据和哨兵数据下载流程
今天,主要跟大家分享一些常用的中高分辨率遥感数据获取的途径及方法. (1)USGS的两个影像下载网站之一:https://glovis.usgs.gov/app,(直接搜索glovis即可找到该网站 ...
GDAL库扩展Landsat系列MTL文件格式支持
Landsat系列卫星提供的数据,一般都是每个波段一个tif文件,然后外加一个MTL.txt的元数据文件,使用gdal可以直接打开每个波段的tif文件,但是有时候想在打开tif数据的同时能够自动读取M ...
Landsat系列数据级别详解
Landsat系列数据级别详解转载自此文:https://www.cnblogs.com/icydengyw/p/12056211.html 一.Landsat Collection 1 Lands ...
AI Earth ——开发者模式案例2：Landsat系列影像数据去云
Landsat 系列数据去云¶ AIE平台中提供的 Landsat 影像均为 USGS Collection 2 中的数据,该集合使用 QA_PIXEL 波段对云.雪进行描述.因此,可通过位运算和逻辑 ...
Landsat系列卫星WRS条带号Path Row分布介绍与对照图
WRS,即Worldwide Reference System,是Landsat系列卫星全球影像标记符号系统,用以区分全球各区域对应的Landsat系列卫星影像编号:其用"Path&qu ...
自然语言处理系列篇——关键词智能提取
自然语言处理系列篇--关键词智能提取:https://zhuanlan.zhihu.com/p/25889937

【IDL】使用IDL对Landsat系列进行批量水体提取

目录