### Python处理GEOTIFF图像裁剪并保持坐标信息的库和方法 在Python中，处理GeoTIFF图像时，可以使用多种库来实现裁剪操作，并确保地理空间坐标信息不丢失。以下是一些常用的库及其方法： #### 1. 使用Rasterio进行裁剪 Rasterio是一个强大的Python库，用于读写栅格地理空间数据。它支持多种栅格数据格式，并能够很好地保留地理空间坐标信息。以下是使用Rasterio裁剪GeoTIFF图像的示例代码： ```python import rasterio from rasterio.mask import mask from shapely.geometry import box import geopandas as gpd # 打开GeoTIFF文件 with rasterio.open("path/to/your/geotiff.tif") as src: # 定义裁剪区域（例如，一个矩形） bbox = box(minx, miny, maxx, maxy) # 替换为实际坐标值 geo = gpd.GeoDataFrame({'geometry': [bbox]}, crs=src.crs) # 裁剪图像 out_image, out_transform = mask(src, geo.geometry, crop=True) out_meta = src.meta # 更新元数据以匹配裁剪后的图像 out_meta.update({"driver": "GTiff", "height": out_image.shape[1], "width": out_image.shape[2], "transform": out_transform}) # 保存裁剪后的图像 with rasterio.open("path/to/cropped_geotiff.tif", "w", **out_meta) as dest: dest.write(out_image) ``` 此方法通过`mask`函数裁剪图像，并利用`out_transform`更新裁剪后的仿射变换矩阵，从而确保坐标信息正确[^2]。 #### 2. 使用GDAL进行裁剪 GDAL是一个广泛使用的开源地理空间数据处理库，支持多种栅格数据格式。虽然GDAL主要是一个C++库，但它提供了Python绑定（osgeo.gdal），使得在Python中使用GDAL变得非常方便。以下是使用GDAL裁剪GeoTIFF图像的示例代码： ```python from osgeo import gdal # 定义裁剪区域的边界（minX, minY, maxX, maxY） clip_bounds = (minX, minY, maxX, maxY) # 打开GeoTIFF文件 input_file = "path/to/your/geotiff.tif" output_file = "path/to/cropped_geotiff.tif" # 使用gdal.Translate进行裁剪 options = gdal.TranslateOptions(projWin=clip_bounds) gdal.Translate(output_file, input_file, options=options) ``` 此方法通过`gdal.Translate`函数指定裁剪区域，并自动保留地理空间坐标信息[^1]。 #### 3. 使用ArcPy进行裁剪 ArcPy是Esri提供的Python模块，专门用于与ArcGIS软件交互。它可以轻松裁剪GeoTIFF图像，并确保地理空间坐标信息不丢失。以下是使用ArcPy裁剪GeoTIFF图像的示例代码： ```python import arcpy # 输入和输出文件路径 input_raster = "path/to/your/geotiff.tif" output_raster = "path/to/cropped_geotiff.tif" # 定义裁剪区域（可以通过Shapefile或直接定义边界框） clip_feature = "path/to/clip_boundary.shp" # 或者使用arcpy.Extent对象 # 使用ExtractByMask工具进行裁剪 arcpy.CheckOutExtension("Spatial") arcpy.sa.ExtractByMask(input_raster, clip_feature).save(output_raster) ``` 此方法通过`ExtractByMask`工具裁剪图像，并自动保留地理空间坐标信息[^4]。 #### 注意事项 - 在裁剪过程中，确保输入的裁剪区域边界与原始图像的坐标系一致。 - 如果裁剪区域由矢量数据（如Shapefile）定义，则需要确保矢量数据与GeoTIFF图像具有相同的坐标参考系统（CRS）。如果不同，需先进行投影转换。 - 对于Rasterio和GDAL，裁剪后会自动生成新的仿射变换矩阵，从而确保地理空间坐标信息正确。