【RS】如何从USGS上下载LANDSAT数据

### 如何从USGS下载Landsat卫星数据 为了自动化下载特定日期的Landsat影像，可以利用Python编写程序来完成此操作[^1]。下面是一个简单的流程说明以及相应的代码片段。 #### 安装必要的库 首先需要安装一些依赖项，比如`earthengine-api`用于访问Google Earth Engine平台上的公开遥感数据集；还有`requests`模块方便发起HTTP请求： ```bash pip install earthengine-api requests ``` #### 初始化Earth Engine并登录账户 接着初始化EE环境并与个人账号关联起来，在命令行输入如下指令完成授权过程（只需执行一次）: ```bash import ee ee.Authenticate() ee.Initialize() ``` #### 构建查询条件筛选目标场景 定义函数接收参数构建过滤器对象，进而选取符合条件的数据条目集合： ```python def get_landsat_image(date, region): # 设置时间范围 start_date = date + '-01' end_date = date + '-31' # 加载Landsat 8 OLI/TIRS集合 landsat_collection = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2') # 应用时间和地理边界约束 filtered = landsat_collection.filterDate(start_date, end_date)\ .filterBounds(region) return filtered.first() if not filtered.isEmpty() else None ``` 此处假设每个月份内只有一景满足需求，因此直接取第一个匹配的结果作为输出。 #### 下载选定的图像文件 最后一步就是把选好的影像保存到本地磁盘上去了。这里给出一段示范性的代码用来展示具体做法： ```python import os from urllib.parse import urlparse import requests def download_file(url, output_dir='.'): parsed_url = urlparse(url) filename = os.path.basename(parsed_url.path) filepath = os.path.join(output_dir, filename) response = requests.get(url, stream=True) total_length = int(response.headers.get('content-length')) with open(filepath, 'wb') as f: for chunk in response.iter_content(chunk_size=1024): if chunk: f.write(chunk) print(f'Downloaded {filename} successfully.') return filepath ``` 以上便是整个工作流的大致框架结构描述。值得注意的是实际应用过程中可能还会涉及到更多细节方面的考量因素，例如异常情况处理机制的设计等。