from osgeo import gdal, ogr # 定义降水、土地类型和坡度的权重 weights = {'precipitation': 0.3, 'land_type': 0.3, 'slope': 0.4} # 打开土地类型矢量数据 land_type_dataset = ogr.Open('path/to/land_types.shp') # 栅格化降水矢量数据和土地类型矢量数据 gdal.Rasterize('precipitation

感谢你的代码分享！根据你提供的代码，我注意到一个问题是在清理变量时未定义 land_type_band。此外，还有一些改进的空间，例如更好地处理异常情况和错误处理。以下是根据你的代码进行了一些修改和改进后的示例：

from osgeo import gdal, ogr

# 定义降水、土地类型和坡度的权重
weights = {'precipitation': 0.3, 'land_type': 0.3, 'slope': 0.4}

# 打开土地类型矢量数据
land_type_dataset = ogr.Open('path/to/land_types.shp')

if land_type_dataset is None:
    print("无法打开土地类型矢量数据")
    exit()

# 栅格化降水矢量数据和土地类型矢量数据
gdal.Rasterize('precipitation_raster.tif', land_type_dataset, attribute='PRECIPITATION')
gdal.Rasterize('land_type_raster.tif', land_type_dataset, attribute='LAND_TYPE')

# 打开坡度栅格数据
slope_raster = gdal.Open('path/to/slope.tif')

if slope_raster is None:
    print("无法打开坡度栅格数据")
    exit()

# 获取栅格波段
precipitation_band = gdal.Open('precipitation_raster.tif').GetRasterBand(1)
land_type_band = gdal.Open('land_type_raster.tif').GetRasterBand(1)
slope_band = slope_raster.GetRasterBand(1)

# 创建输出数据集
driver = gdal.GetDriverByName('GTiff')
out_dataset = driver.Create('suitability.tif', slope_band.XSize, slope_band.YSize, 1, gdal.GDT_Float32)
out_band = out_dataset.GetRasterBand(1)

# 进行叠加分析
for i in range(precipitation_band.YSize):
    precipitation_row = precipitation_band.ReadAsArray(0, i, precipitation_band.XSize, 1)
    land_type_row = land_type_band.ReadAsArray(0, i, land_type_band.XSize, 1)
    slope_row = slope_band.ReadAsArray(0, i, slope_band.XSize, 1)
    
    # 假设每个栅格数据的值已经经过标准化处理，可以直接相乘
    suitability_row = weights['precipitation'] * precipitation_row + weights['land_type'] * land_type_row + weights['slope'] * slope_row
    out_band.WriteArray(suitability_row, 0, i)

# 设置输出数据集的地理变换和投影
out_dataset.SetGeoTransform(slope_raster.GetGeoTransform())
out_dataset.SetProjection(slope_raster.GetProjection())

# 清理
out_band.FlushCache()
out_dataset = None
precipitation_raster = None
land_type_raster = None
slope_raster = None

print("完成")

请确保替换代码中的路径为实际数据文件的路径，并根据需要进行其他适当的修改。同时，这只是一个示例，具体的操作和处理还需根据你的需求进行调整。

内容由零声教学AI助手提供，问题来源于学员提问