landsat7数据裁剪

### 使用GIS工具或遥感软件裁剪Landsat 7卫星影像数据

对于地理空间数据分析而言，裁剪 Landsat 7 卫星影像是一项常见的预处理操作。这可以通过多种 GIS 工具和遥感软件来完成，如 QGIS、ArcGIS 或者 ERDAS Imagine。

#### 方法一：使用QGIS进行裁剪

在 QGIS 中执行此任务相对简单直观：

1. 打开 QGIS 并加载要裁剪的 Landsat 7 影像文件。
2. 添加用于定义裁剪范围的矢量层（通常是面要素），也可以手动绘制一个多边形作为裁剪边界。
3. 转到 `Raster` -> `Extraction` -> `Clip Raster by Mask Layer` 功能菜单项。
4. 设置输入栅格图层为待裁剪的 Landsat 7 数据集；设置掩模图层为你之前准备好的矢量或多边形。
5. 配置输出选项并点击运行按钮即可获得经过裁剪后的结果[^1]。

import os
from qgis.core import QgsProject, QgsVectorLayer, QgsRasterLayer, QgsProcessingFeedback, QgsProcessingContext
from qgis.analysis import QgsNativeAlgorithms

context = QgsProcessingContext()
feedback = QgsProcessingFeedback()

# 加载栅格图层
raster_layer_path = 'path_to_landsat_image.tif'
raster_layer = QgsRasterLayer(raster_layer_path, "Landsat Image")

# 加载矢量图层(裁剪边界)
vector_layer_path = 'path_to_clip_boundary.shp'
vector_layer = QgsVectorLayer(vector_layer_path, "Boundary", "ogr")

# 运行裁剪算法
output_file = '/tmp/clipped_raster.tif'
processing.run("native:cliprasterbymasklayer",
               {'INPUT': raster_layer,
                'MASK': vector_layer,
                'OUTPUT': output_file},
               context=context,
               feedback=feedback)

print(f"Cropped image saved to {output_file}")

#### 方法二：利用 ArcGIS Pro 实现相同目的

如果倾向于使用 ESRI 的产品，则可以在 ArcGIS Pro 中按照如下方式操作：

1. 启动应用程序并将目标 Landsat 7 文件拖放到项目内。
2. 创建一个新的要素类或者导入现有的用来指定裁切区间的几何对象。
3. 寻找名为 “Extract by Mask”的工具，在其对话框里指明源影像以及遮罩形状。
4. 完成参数配置之后按下 OK 键等待处理结束，新生成的数据会自动加入地图视窗中显示出来[^2]。

以上两种方法均能有效地帮助用户根据特定需求调整 Landsat 7 图像的空间范围，从而更好地服务于后续的研究工作或是应用场合。