Python GDAL 栅格叠加运算：详细教程

要在Python中使用gdal实现栅格的叠加运算，你可以按照以下步骤进行操作：

1. 导入必要的库：

import gdal
import numpy as np

2. 打开要进行叠加运算的栅格文件：

dataset1 = gdal.Open('path/to/raster1.tif')
dataset2 = gdal.Open('path/to/raster2.tif')

3. 获取栅格的投影、地理转换和大小信息：

projection = dataset1.GetProjection()
geotransform = dataset1.GetGeoTransform()
rows, cols = dataset1.RasterYSize, dataset1.RasterXSize

4. 读取栅格数据为Numpy数组：

band1 = dataset1.GetRasterBand(1)
band2 = dataset2.GetRasterBand(1)
array1 = band1.ReadAsArray()
array2 = band2.ReadAsArray()

5. 进行栅格叠加运算：

result = array1 + array2

6. 创建输出文件：

driver = gdal.GetDriverByName('GTiff')
output = driver.Create('path/to/output.tif', cols, rows, 1, band1.DataType)
output.SetProjection(projection)
output.SetGeoTransform(geotransform)

7. 将结果写入输出文件：

output.GetRasterBand(1).WriteArray(result)
output.FlushCache()

8. 关闭文件：

dataset1 = None
dataset2 = None
output = None

这样，你就可以使用gdal实现栅格的叠加运算了。请注意，这只是一个简单的示例，实际操作中可能需要进行更多的错误处理和参数设置。