请问怎么Qgis河流提取出来图长这样

Question

按照网上的教程在Qgis里面对地块依据dem进行河流一键提取，为什么生成的图长这样不是那种连贯的？同时生成流向图之后会这样糊成一片黑又是为什么？ಠ_ಠ

Answer 1

下午好🌅🌅🌅
本答案参考通义千问

你提到在QGIS中使用DEM数据进行河流提取后，生成的河流图看起来不连贯，流向图呈现一片黑色，这通常是由以下原因导致的。以下是详细的分析和解决方案：

---

一、问题分析

1. 河流提取结果不连贯

• 可能原因：
  • DEM数据质量差（如噪声多、缺失值未处理）。
  • 河流提取参数设置不合理（如阈值太小或太大）。
  • 河流网络生成时未正确连接断开的河道。
  • 没有进行流域划分或河网提取的后续处理。

2. 流向图糊成一片黑

• 可能原因：
  • 流向图的渲染方式不正确（如未正确设置颜色映射）。
  • 数据范围过大，导致图像无法清晰显示。
  • 未正确使用流向插件（如SAGA、GRASS等）生成流向数据。
  • 输出文件格式不支持或未正确导出。

---

二、解决方法（详细步骤）

1. 检查并预处理DEM数据

重点：确保DEM数据干净且无空值

• 打开QGIS，加载你的DEM数据。
• 使用 “Raster Calculator” 工具，将空值替换为0或其他合理数值：

  ("dem@1" > 0) * "dem@1"
  

  或者使用 “Fill NoData” 插件填充空值。

加粗重点： 请务必对DEM数据进行预处理，否则后续所有操作都会受到影响。

---

2. 使用正确的工具进行河流提取

推荐工具：

• SAGA GIS 的“Stream Network”模块
• GRASS GIS 的 “r.watershed” 模块

步骤如下：

1. 计算流向（Flow Direction）

   • 在QGIS中，使用 SAGA 或 GRASS 的“Flow Direction”工具。
   • 输入你的DEM数据，输出一个流向栅格（Flow Direction）。

2. 计算流量（Flow Accumulation）

   • 使用 SAGA 的“Flow Accumulation”工具，输入流向栅格。
   • 这一步会得到一个表示单位面积内汇水的栅格。

3. 提取河流（Stream Network）

   • 使用 SAGA 的“Stream Network”工具，设置合适的阈值（Threshold）。
   • 阈值越大，提取的河流越少；阈值越小，提取的河流越多。
   • 建议根据实际地形调整，比如从500开始尝试。

4. 将栅格转为矢量（可选）

   • 使用 “Raster to Vector” 工具，将河流栅格转为线状矢量。

---

3. 修正河流不连贯的问题

重点：确保河流网络是连续的

• 如果河流网络出现断裂，可以尝试以下方法：
  • 调整阈值，使其更合适。
  • 使用 “Vector Geometry” 工具中的 “Snap Geometries to Grid” 或 “Merge Lines” 工具，手动连接断裂的河流段。
  • 使用 “v.clean” 工具（GRASS）清理矢量数据。

---

4. 处理流向图显示异常

重点：确保流向图正确渲染

• 如果流向图显示为一片黑，可能是由于：
  • 颜色映射错误：在图层属性中选择正确的渲染类型（如“Categorized”或“Singleband Pseudocolor”）。
  • 数据范围过大：尝试将流向图的最大值和最小值手动设置为0-8（因为流向一般为0~7的整数）。
  • 输出格式问题：确保输出的是GeoTIFF，而不是其他格式。

示例代码（Python + GDAL）：

from osgeo import gdal
import numpy as np

# 读取流向栅格
ds = gdal.Open('flow_direction.tif')
band = ds.GetRasterBand(1)
arr = band.ReadAsArray()

# 设置颜色映射
color_table = gdal.ColorTable()
for i in range(8):
    color_table.SetColorEntry(i, (255, 255, 255, 255))  # 白色

# 保存为新的GeoTIFF
driver = gdal.GetDriverByName('GTiff')
out_ds = driver.Create('flow_direction_colored.tif', ds.RasterXSize, ds.RasterYSize, 1, band.DataType)
out_band = out_ds.GetRasterBand(1)
out_band.WriteArray(arr)
out_band.SetColorTable(color_table)
out_ds.FlushCache()

---

三、推荐流程总结（有序列表）

1. 预处理DEM数据

   • 填充空值、去除噪声、确保数据完整。

2. 计算流向（Flow Direction）

   • 使用SAGA或GRASS的“Flow Direction”工具。

3. 计算流量（Flow Accumulation）

   • 用于识别集水区和河流网络。

4. 提取河流（Stream Network）

   • 使用SAGA的“Stream Network”工具，设置合适的阈值。

5. 将河流转为矢量（可选）

   • 方便后续编辑和可视化。

6. 检查流向图渲染

   • 确保颜色映射正确，数据范围合理。

7. 修复河流断开问题

   • 使用矢量工具合并或连接断开的河流段。

---

四、常见错误排查表

| 问题 | 原因 | 解决方案 | |------|------|----------| | 河流不连贯 | DEM数据质量差或阈值设置不当 | 预处理DEM，调整阈值 | | 流向图一片黑 | 渲染方式错误或数据范围问题 | 设置颜色映射，限制数值范围 | | 河流提取失败 | 工具未正确配置 | 检查SAGA/GRASS插件是否启用 |

---

如果你能提供具体的DEM数据截图或部分操作日志，我可以进一步帮你定位具体问题。希望这些信息对你有帮助！