### 使用Python库裁剪遥感影像为正方形 为了实现将遥感影像裁剪为正方形区域的目标，可以利用 `GDAL` 或 `Rasterio` 这些强大的 Python 库。以下是具体方法： #### 方法一：基于 GDAL 的解决方案 通过设置目标窗口的左上角坐标以及宽度和高度参数，能够精确控制裁剪范围并确保其为正方形。 ```python from osgeo import gdal def clip_raster_to_square(input_file, output_file, top_left_x, top_left_y, size): """ 裁剪栅格数据到指定大小的正方形。 参数: input_file (str): 输入的 TIFF 文件路径。 output_file (str): 输出的裁剪后的 TIFF 文件路径。 top_left_x (float): 正方形左上角 X 坐标（投影单位）。 top_left_y (float): 正方形左上角 Y 坐标（投影单位）。 size (int): 正方形边长（以像元数表示）。 """ options = f"-projwin {top_left_x} {top_left_y} {top_left_x + size} {top_left_y - size}" command = f'gdal_translate {options} "{input_file}" "{output_file}"' gdal.Translate(output_file, input_file, projWin=[top_left_x, top_left_y, top_left_x + size, top_left_y - size]) # 示例调用 clip_raster_to_square('input.tif', 'output_square.tif', 100000, 200000, 500) # 边长为500像素 ``` 上述代码片段展示了如何定义一个函数来完成此操作，并提供了具体的命令行选项解释[^1]。 #### 方法二：基于 Rasterio 的解决方案 Rasterio 提供了一种更直观的方式处理地理空间数据集。它允许读取、写入 GeoTIFF 并支持多种格式转换等功能。 ```python import rasterio from rasterio.windows import Window def clip_with_rasterio(input_path, output_path, row_start, col_start, square_size): with rasterio.open(input_path) as src: window = Window(col_start, row_start, square_size, square_size) kwargs = src.meta.copy() data = src.read(window=window) # 更新元数据中的宽高信息 kwargs.update({ "height": data.shape[1], "width": data.shape[2], "transform": rasterio.windows.transform(window, src.transform), }) with rasterio.open(output_path, "w", **kwargs) as dst: dst.write(data) # 示例调用 clip_with_rasterio('example.tif', 'clipped_example.tif', 100, 100, 200) # 切割位置从(100,100)，尺寸为200*200 ``` 这段脚本说明了怎样打开原始图像文件，创建一个新的子图对象作为视窗传给源对象用于提取所需部分的数据数组；接着调整输出文件的相关属性以便保存新生成的小块图像[^2]。 注意，在实际应用过程中可能还需要考虑其他因素比如边界条件判断等细节问题。