遥感新手必看：手把手教你从USGS EarthExplore下载Landsat 5数据（附云量筛选技巧）-平芜编程栈

遥感数据获取实战指南：从USGS EarthExplore高效下载Landsat 5数据

第一次接触遥感数据分析时，最令人头疼的往往不是算法或软件操作，而是如何获取高质量的基础数据。作为全球历史最悠久的对地观测计划之一，Landsat系列卫星数据因其免费开放、时间跨度长（Landsat 5运行时间达28年）、空间分辨率适中（30米）等特点，成为学术研究和行业应用的首选。本文将聚焦USGS EarthExplore平台，以Landsat 5 TM数据为例，手把手教你避开数据下载中的常见陷阱，特别是如何根据研究目的智能筛选云量，以及理解不同处理级别数据的核心差异。

1. 准备工作与环境配置

1.1 注册USGS EarthExplore账号

访问https://earthexplorer.usgs.gov/点击右上角"Login"按钮。对于首次使用的用户，需要先完成账号注册：

点击"Register"进入注册页面
填写基本信息（姓名、邮箱、机构等）
验证邮箱后设置密码
登录后建议在"Profile"中完善研究领域信息

注意：部分教育网用户可能会遇到访问缓慢的情况，建议尝试切换网络环境或使用国际网络出口

1.2 了解Landsat 5数据基本参数

在开始搜索前，需要明确几个关键参数：

参数名称	典型值	说明
空间分辨率	30m (可见光/近红外)	热红外波段分辨率为120m
时间分辨率	16天	同一区域的重访周期
幅宽	185km	单景影像的覆盖范围
波段数量	7个	B1-B7对应不同光谱范围

# 典型波段组合示例（植被分析常用） vegetation_bands = { 'Band 3': 'Green (0.52-0.60μm)', 'Band 4': 'Near Infrared (0.76-0.90μm)', 'Band 5': 'Shortwave Infrared 1 (1.55-1.75μm)' }

2. 数据检索全流程详解

2.1 区域选择的三种实用方法

在"Search Criteria"标签页中，有多个方式可以定义研究区域：

Path/Row编号法（最精确）
- 通过WRS-2网格系统查询目标区域的Path和Row值
- 适合已知具体研究区坐标的情况
地图绘制法（最直观）
- 点击"Polygon"工具手动绘制区域边界
- 支持导入KML/Shapefile文件
- 适用于不规则研究区域
地名搜索法（最快捷）
- 在"Address/Place"输入城市名或地名
- 系统自动匹配到相应坐标范围

提示：对于大范围研究区，建议分多个小区域下载，避免单次结果过多导致系统超时

2.2 时间范围与云量筛选策略

在"Date Range"设置中，时间跨度不宜过长（建议不超过6个月），否则会导致结果过多。云量筛选是获取可用数据的关键步骤：

城市变化检测：建议云量<10%
植被监测：可放宽至<30%（云阴影影响较小）
水体研究：需<5%（水面反射易受云层干扰）

# 云量筛选逻辑示例 def cloud_filter(scene_list, max_cloud_cover): return [scene for scene in scene_list if scene['cloud_cover'] <= max_cloud_cover]

实际操作中，可以先将云量阈值设为20%进行初筛，然后在结果页面根据缩略图二次筛选。

3. 数据产品类型深度解析

3.1 Collection 1与Collection 2的区别

USGS目前提供两个版本的数据产品：

特性	Collection 1	Collection 2
时间范围	1972-2021	1982-至今
几何校正	L1TP/L1GT	统一L1TP标准
辐射定标	需要手动处理	已优化
建议使用场景	历史数据分析	新项目首选

3.2 Level-1与Level-2的核心差异

Level-1数据是经过辐射定标和几何校正的原始数据，而Level-2产品额外包含：

大气校正结果（地表反射率）
地表温度产品（仅限热红外波段）
质量评估波段（QA）

对于初学者，建议优先选择Level-2数据，除非有特殊处理需求。下表对比了两种级别的典型应用：

研究领域	推荐级别	原因
植被指数计算	Level-2	直接使用地表反射率
城市热岛分析	Level-2	包含温度产品
辐射传输模拟	Level-1	需要原始辐射值
传感器验证	Level-1	保留完整辐射特性

4. 下载优化与常见问题解决

4.1 多波段下载技巧

虽然EarthExplore支持多波段选择，但下载时可能会遇到以下情况：

批量下载失败：网络不稳定时建议单波段下载
部分波段缺失：检查该景数据是否完整归档
下载速度慢：尝试非高峰时段（UTC时间凌晨2-5点）

推荐使用Bulk Download Application工具进行大批量下载管理：

# 安装USGS批量下载工具 wget https://dds.cr.usgs.gov/bulkdownloads/BDM_Setup.exe chmod +x BDM_Setup.exe ./BDM_Setup.exe

4.2 数据质量验证方法

下载完成后，建议通过以下步骤验证数据完整性：

检查各波段文件大小是否正常（通常应>30MB）
使用QGIS或ENVI快速预览影像
查看元数据文件（*_MTL.txt）中的关键参数：
- SUN_ELEVATION（太阳高度角）
- DATA_TYPE（数据类型）
- CLOUD_COVER（实际云量）

遇到问题时的排查顺序：网络连接 → 账号权限 → 数据可用性 → 软件兼容性。对于持续出现的问题，可以截图保存错误信息，通过USGS的"Contact Us"页面提交支持请求。

5. 进阶技巧与个性化设置

5.1 创建数据提醒服务

对于长期监测项目，可以设置自动提醒：

在"Search Criteria"中保存当前搜索条件
点击"Create Alert"设置通知频率
选择邮件或站内信通知方式

5.2 API自动化访问

对于需要定期获取数据的用户，可以考虑使用USGS的Machine-to-Machine接口：

import requests api_url = "https://m2m.cr.usgs.gov/api/api/json/stable/" auth_payload = { "username": "your_username", "password": "your_password" } # 获取访问令牌 response = requests.post(api_url + "login", json=auth_payload) token = response.json()['data']