python转换tif数据投影

### 使用Python将TIF数据的投影进行转换在GIS领域，使用GDAL或Rasterio库可以轻松实现TIF数据的投影转换。以下是基于GDAL和Rasterio两种方法的代码示例。 #### 方法一：使用GDAL进行投影转换 GDAL提供了强大的功能来处理栅格数据，包括投影转换。以下是一个完整的代码示例[^1]： ```python from osgeo import gdal def convert_projection(input_tif, output_tif, new_crs): # 打开输入文件 dataset = gdal.Open(input_tif) # 定义目标CRS（如EPSG:4326） target_crs = f'EPSG:{new_crs}' # 使用gdal.Warp进行投影转换 gdal.Warp(output_tif, dataset, dstSRS=target_crs) # 示例调用 input_tif = "./Image/input.tif" output_tif = "./Image/output_converted.tif" new_crs = 4326 # 目标坐标系 EPSG代码 convert_projection(input_tif, output_tif, new_crs) ``` 此代码中，`gdal.Warp`函数是关键，它不仅支持投影转换，还可以用于裁剪、重采样等操作。 --- #### 方法二：使用Rasterio进行投影转换 Rasterio是另一个流行的Python库，专门用于处理地理空间栅格数据。以下是使用Rasterio实现投影转换的代码示例[^2]： ```python import rasterio from rasterio.warp import calculate_default_transform, reproject, Resampling def convert_projection_rasterio(input_tif, output_tif, dst_crs): with rasterio.open(input_tif) as src: transform, width, height = calculate_default_transform( src.crs, dst_crs, src.width, src.height, *src.bounds) kwargs = src.meta.copy() kwargs.update({ 'crs': dst_crs, 'transform': transform, 'width': width, 'height': height }) with rasterio.open(output_tif, 'w', **kwargs) as dst: for i in range(1, src.count + 1): reproject( source=rasterio.band(src, i), destination=rasterio.band(dst, i), src_transform=src.transform, src_crs=src.crs, dst_transform=transform, dst_crs=dst_crs, resampling=Resampling.nearest) # 示例调用 input_tif = "./Image/input.tif" output_tif = "./Image/output_converted_rasterio.tif" dst_crs = "EPSG:4326" # 目标坐标系 convert_projection_rasterio(input_tif, output_tif, dst_crs) ``` 在上述代码中，`calculate_default_transform`用于计算新的变换矩阵，而`reproject`则完成实际的投影转换操作。 --- ### 注意事项 - 在使用GDAL或Rasterio时，请确保已正确安装相关依赖项。 - 投影转换的目标坐标系通常以`EPSG`代码表示，例如`EPSG:4326`表示WGS84坐标系。 - 如果输入数据未定义坐标系，可能需要手动指定原始坐标系，否则可能导致错误[^3]。 ---

创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考

下一篇随便写一个python代码