Arc gis 栅格计算器加权叠加结果有白色区域怎么回事呢

Question

求解答，根据得出的权重在栅格计算器里对每个评价因子栅格图层进行叠加，得出的最终结果有白值，我感觉这部分白色区域像是我评价因子里面夜间灯光➕人口密度➕河网密度，请问是这三个图层本身就有问题吗

Answer 1

引自免费微信小程序：皆我百晓生

在 ArcGIS 中，当使用栅格计算器进行加权叠加时，出现白色区域可能是由于以下几个原因：

1. 数据范围不一致：确保所有输入栅格数据具有相同的地理空间参考（坐标系统）、投影和分辨率。
2. 权重设置错误：检查权重是否正确赋给各个栅格图层，并且它们都在有效范围内（通常为0到1之间）。
3. 数据溢出：如果加权之后的总和超过了栅格数据类型的最大值（例如，8位灰度图像的最大值为255），则超出部分可能会显示为白色。
4. 数据缺失或NoData值处理不当：确认所有栅格图层中的NoData值是否被正确处理，比如设置为0或其他合适的值。

解决这个问题的步骤可以分为以下几个部分：

步骤一：确保数据一致性

首先，你需要确保所有输入栅格的数据格式、投影和分辨率都是一致的。这可以通过 ArcGIS 的 "Raster Calculator" 或 "Raster Project" 工具实现。以下是使用 "Raster Calculator" 调整分辨率的示例代码：

# 假设你的三个栅格分别为raster_light, raster_population, 和 raster_river_density
from arcgis.raster.functions import *
from arcgis.geometry import SpatialReference

def project_raster(raster, target_res):
    sr = raster.spatial_reference
    new_sr = sr.projectAs(target_res)
    return arcpy.Resample_management(raster, "new_" + raster.name, str(target_res.MetersPerUnit), cell_size=target_res.MetersPerUnit)

target_projection = # 设置目标投影
target_resolution = # 设置目标分辨率 (例如，100米)

raster_light_projected = project_raster(raster_light, target_projection)
raster_population_projected = project_raster(raster_population, target_projection)
raster_river_density_projected = project_raster(raster_river_density, target_projection)

步骤二：检查权重

确保你的权重设置正确，没有负数或超过1的值。权重应该在0到1之间。这里是计算加权叠加的代码示例：

weighted_light = raster_light_projected * weight_light  # 假设weight_light是夜灯的权重
weighted_population = raster_population_projected * weight_population  # 假设weight_population是人口的权重
weighted_river_density = raster_river_density_projected * weight_river_density  # 假设weight_river_density是河流密度的权重

# 检查权重是否在有效范围内
assert 0 <= weight_light <= 1
assert 0 <= weight_population <= 1
assert 0 <= weight_river_density <= 1

步骤三：处理NoData值

将NoData值转换为0，避免影响叠加结果：

no_data_value = raster_light_projected.no_data_value
weighted_light = Con(IsNull(weighted_light), no_data_value, weighted_light)
weighted_population = Con(IsNull(weighted_population), no_data_value, weighted_population)
weighted_river_density = Con(IsNull(weighted_river_density), no_data_value, weighted_river_density)

步骤四：防止数据溢出

在叠加之前，先将每个栅格的值限制在数据类型允许的最大值之内，这里假设我们使用的是8位灰度图像（最大值为255）：

max_value = 255

weighted_light_clipped = LimitValue(weighted_light, min_value=0, max_value=max_value)
weighted_population_clipped = LimitValue(weighted_population, min_value=0, max_value=max_value)
weighted_river_density_clipped = LimitValue(weighted_river_density, min_value=0, max_value=max_value)

步骤五：叠加

最后，执行加权叠加并输出结果：

final_result = weighted_light_clipped + weighted_population_clipped + weighted_river_density_clipped

arcpy.RasterToPolygon_conversion(final_result, "output_polygon.shp", "VALUE")

以上代码示例使用了Python ArcPy库，你需要在ArcGIS环境中运行这些代码。确保已经安装了ArcGIS Desktop或Pro，并配置了正确的环境变量。如果你在运行过程中遇到任何问题，请提供具体错误信息以便进一步分析。