Landsat数据在USGS中无法下载Surface Reflectance产品的解决方法

本文介绍在USGS官网下载Landsat遥感影像数据时，出现报错信息，无法下载地表反射率产品（Surface Reflectance）的解决办法。

最近，利用这篇文章批量下载Landsat遥感影像的方法中遥感影像在USGS官网批量下载方法下载Landsat数据时，出现了如下报错：

Errors were found in your request.1 validation errors→Remove invalid product request 'sr' in oli8_collection.products - Available products: ['source_metadata', 'l1', 'toa', 'aq_refl', 'stats', 'pixel_qa']

如下图所示。

这一错误我仅仅在下载一年当中部分时段的Landsat 8数据时才会遇到，因此证明可能是部分时段的遥感影像数据缺失导致。经过不断探索，终于发现问题所在：官网中提到，部分时段由于对应的初始遥感数据或计算地表反射率的环境变量数据缺失，因此Landsat 7与Landsat 8部分时段的影像暂时没有地表反射率产品，如下图（如果大家下载影像的时间不在下表所列时间范围，但还是出现了类似问题，看本文末尾即可）。

因此，我们需要结合上述博客中的批量下载方法，对导入订单时所用到的相关信息加以修改。

例如，如果我需要下载Landsat 8在2019年11月至12月的数据，那么依据上表，我就需要将12月20日与21日的数据直接剔除掉。

然后将新的.txt文件重新导入，生成订单，即可不报错：

同时，还需要注意的是，有些时段的遥感影像并不是所有分幅都有问题，例如表中红色的部分。如果我们需要用对应时间段内的遥感数据，但是遥感影像所在的Path号不是表中的号，那么其实也不会有问题。

例如，如果我需要下载Landsat 8在2020年11月至12月的数据，那么依据上表，我就需要将11月01日至08日、11月12日至13日的数据直接全部剔除掉；同时需要对11月09日与11月14日的数据加以检查，尽管这两天的数据不用全部剔除，但如果其有表内对应列出的Path号，那么也需要剔除。

首先，剔除11月01日至08日、11月12日至13日的全部数据。

其次，剔除11月09日Path号为121、137的图像。

11月14日同理，剔除Path为124的，上面的140由于没有问题，因此不用剔除。

同样可以成功生成订单了。

接下来，我们来看看一个新的问题。在下载我国范围内2015年01月至2015年04月的Landsat 8遥感影像时，同样出现了本文开头所示的这一报错，导致不能成功生成订单；但是这一时间范围并未在本文前述的表格中所提及。那么是怎么回事呢？

我们在官网页面继续往下看，可以看到这句话：

也就是说，凡是Landsat 8遥感影像中只有OLI传感器数据（也就是不含TIRS传感器数据）的图像，也是无法计算地表反射率的。那么，上面2015年这个数据是不是因为这个原因而导致不能生成订单呢？我们来验证一下。

首先，将无法生成订单的2015年01月至2015年04月的Landsat 8遥感影像所对应的.csv文件打开，在Sensor Identifier这一列做筛选，如下图所示。

可以看到，其中是有OLI与OLI_TIRS两个类别的数据的。

那么，我们再来看一个没有问题的数据。这里我以我国2015年05月至2015年08月的Landsat 8数据为例，这个时间范围内的遥感影像数据是没有问题的（也就是说可以成功生成订单的）。打开其所对应的.csv文件，同样的在Sensor Identifier这一列做筛选，如下图所示。