GIS研习社 GIS研习社
GIS基础理论 Band Math ENVI

ENVI怎么计算NDWI指数?Band Math怎么用?

GIS基础理论
Dr.GIS
wowwwai GIS研习社 · 工具流程与项目排障

引言：很多同学在做水体提取时会问：ENVI怎么计算NDWI指数?Band Math怎么用? 这其实是一个典型的遥感指数计算问题。只要影像波段对应正确，并且在 ENVI 的 Band Math 中写对表达式，就可以快速生成 NDWI 水体指数栅格，用于河流、湖泊、水库、湿地等水体识别。

本文以 ENVI 软件为例，讲清楚 NDWI 指数的计算公式、Band Math 的使用方法、不同卫星影像的波段选择、结果验证方式，以及常见报错和误区。适合刚接触遥感影像处理、GIS 制图和水体提取的读者。

ENVI计算NDWI指数 Band Math使用流程
ENVI 使用 Band Math 计算 NDWI 指数的基本流程：确认波段、输入公式、输出指数栅格。

背景：为什么要在 ENVI 中计算 NDWI 指数

背景：NDWI，全称 Normalized Difference Water Index，通常翻译为归一化水体指数。它利用水体在绿光波段反射较强、在近红外波段反射较弱的特点，把水体和非水体拉开差异。

在实际 GIS 和遥感项目中，NDWI 常用于：

  • 从 Landsat、Sentinel-2、高分、高分辨率无人机多光谱影像中提取水体。
  • 辅助制作河流、湖泊、水库、湿地分布图。
  • 进行洪水淹没范围、岸线变化、水域面积变化分析。
  • 作为分类前的特征波段，提高水体识别效果。

ENVI 的 Band Math 是一个非常常用的栅格计算工具。它类似于 ArcGIS 栅格计算器或 QGIS 栅格计算器，可以对不同波段进行加、减、乘、除、逻辑判断和条件运算。计算 NDWI 时，Band Math 是最直接的方法之一。

原理：NDWI 公式和波段应该怎么选

原理：常用的 McFeeters NDWI 公式如下：

NDWI = (Green - NIR) / (Green + NIR)

其中：

  • Green 表示绿光波段。
  • NIR 表示近红外波段。
  • NDWI 结果通常在 -1 到 1 之间。
  • 水体一般表现为较高的 NDWI 值，植被、建筑物、裸地通常较低。

不同遥感影像的绿光波段和近红外波段编号不一样，计算前一定要先确认波段。下面是常见数据源的参考关系：

数据源 绿光波段 Green 近红外波段 NIR 常用 NDWI 表达式
Landsat 8/9 OLI Band 3 Band 5 (B3 – B5) / (B3 + B5)
Landsat 5/7 TM/ETM+ Band 2 Band 4 (B2 – B4) / (B2 + B4)
Sentinel-2 MSI Band 3 Band 8 (B3 – B8) / (B3 + B8)
国产高分多光谱数据 需查看元数据 需查看元数据 按传感器波段表确定

需要注意，NDWI 有多个变体。例如有些文章会使用改进归一化差异水体指数 MNDWI，它通常使用绿光波段和短波红外波段：

MNDWI = (Green - SWIR) / (Green + SWIR)

本文重点讲 ENVI 怎么计算 NDWI 指数，也就是使用绿光波段和近红外波段的常规 NDWI。不要把 NDWI 和 MNDWI 的波段混用。

步骤：ENVI Band Math 计算 NDWI 指数完整流程

步骤：下面以多光谱影像已经完成辐射定标、大气校正或至少已经具备可比波段值为前提，介绍 ENVI 中使用 Band Math 计算 NDWI 的操作方法。

第 1 步：打开影像并检查波段信息

  1. 启动 ENVI。
  2. 选择 File，打开需要计算 NDWI 的多光谱影像。
  3. 在 Layer Manager 或 Data Manager 中查看影像波段列表。
  4. 确认哪个波段是 Green，哪个波段是 NIR。

如果是 Landsat 8/9 OLI，一般 Green 是 Band 3，NIR 是 Band 5。如果是 Sentinel-2，Green 一般是 Band 3，NIR 一般是 Band 8。但在实际项目中，不建议只凭经验操作，最好查看影像元数据、波段名称或传感器说明。

第 2 步：打开 Band Math 工具

  1. 在 ENVI 主界面中，找到 Toolbox。
  2. 依次进入 Raster Management 或使用工具搜索框。
  3. 搜索并打开 Band Math。

不同版本 ENVI 的菜单位置可能略有差异。如果找不到，可以直接在 Toolbox 搜索框输入 Band Math。

第 3 步：输入 NDWI 表达式

在 Band Math 的表达式窗口中输入公式。为了避免整数除法或数据类型问题，建议在表达式中加入浮点写法。

通用写法如下：

float(b1 - b2) / float(b1 + b2)

这里的 b1 和 b2 只是变量名，并不是固定波段号。通常可以设置：

  • b1 对应 Green 绿光波段。
  • b2 对应 NIR 近红外波段。

如果你使用 Landsat 8/9，公式逻辑相当于：

float(Band3 - Band5) / float(Band3 + Band5)

如果你使用 Sentinel-2，公式逻辑相当于：

float(Band3 - Band8) / float(Band3 + Band8)

但在 ENVI Band Math 中，一般不是直接输入 Band3、Band5，而是输入变量 b1、b2，然后在后续界面把变量绑定到对应波段。

第 4 步：为变量选择对应波段

  1. 点击确认表达式后，ENVI 会要求为变量 b1、b2 指定输入数据。
  2. 选择 b1 时，指定绿光波段 Green。
  3. 选择 b2 时，指定近红外波段 NIR。
  4. 检查变量和波段是否对应正确。

这是 ENVI 计算 NDWI 指数最容易出错的地方。公式写对了，但 b1、b2 绑定反了，结果就会变成相反的指数，水体可能显示为低值甚至负值。

第 5 步：设置输出文件

  1. 选择输出路径。
  2. 建议输出为 ENVI 标准栅格或 GeoTIFF，方便后续在 ArcGIS Pro、QGIS 中使用。
  3. 给结果文件命名，例如 ndwi_landsat8_202405.tif。
  4. 点击运行。

输出完成后，ENVI 会生成一幅单波段 NDWI 栅格。该栅格可以用于显示、阈值分割、面积统计或进一步分类。

第 6 步：显示和拉伸 NDWI 结果

NDWI 是单波段指数影像，默认显示效果可能比较灰。可以进行拉伸显示：

  • 使用灰度或蓝色渐变色带显示 NDWI。
  • 查看水体区域是否明显偏亮或偏高值。
  • 使用 Cursor Value 或像元查询工具检查水体、植被、裸地的 NDWI 值。

一般情况下，水体 NDWI 值相对较高，但具体阈值不是固定的。常见经验阈值可以从 0 附近开始试验，例如 NDWI 大于 0 初步判定为水体，但必须结合影像类型、季节、阴影、浑浊水体和研究区环境进行调整。

常见坑：ENVI 计算 NDWI 指数容易出错的地方

常见坑：NDWI 公式本身很简单，但在 ENVI Band Math 里经常出现结果异常、全黑、值域不对、水体提取不准等问题。下面这些问题建议逐项检查。

1. 绿光波段和近红外波段选错

这是最常见的问题。不同传感器波段编号不一样，不能把 Landsat 的波段编号直接套到 Sentinel-2 或其他国产影像上。

  • Landsat 8/9：Green 是 Band 3，NIR 是 Band 5。
  • Sentinel-2：Green 是 Band 3，NIR 是 Band 8。
  • Landsat 5/7：Green 是 Band 2，NIR 是 Band 4。

如果波段选错，NDWI 结果可能仍然能算出来，但数值没有物理意义。

2. b1 和 b2 变量绑定反了

公式写成 float(b1 - b2) / float(b1 + b2) 时，b1 必须是 Green，b2 必须是 NIR。如果反过来，结果相当于取了相反数，水体高值和低值会颠倒。

3. 没有使用浮点计算

如果输入影像是整型数据，直接相除可能导致结果精度异常。建议使用：

float(b1 - b2) / float(b1 + b2)

这样可以让 ENVI 按浮点方式计算，输出结果更符合 NDWI 的连续指数特征。

4. 分母为 0 或异常值没有处理

理论公式中分母是 Green + NIR。如果某些像元两个波段值都为 0，可能出现无效值。通常在遥感影像的背景区、无数据区、裁剪外区域比较常见。

如果影像有明显背景值，建议先设置 NoData、做掩膜，或者在后续阈值提取时排除背景区域。

5. 直接用未校正影像做严格分析

如果只是课堂练习或快速制图，使用原始 DN 值有时也能看到水体差异。但如果要做时间序列对比、水域面积变化监测或论文分析，建议使用经过辐射定标和大气校正后的地表反射率数据。

6. 阈值固定为 0 导致误分

NDWI 大于 0 常被用作水体初始阈值，但并不是所有地区都适用。城市阴影、山体阴影、湿地植被、浑浊水体、浅水区都可能影响结果。

建议结合直方图、目视解译、样本点和高分辨率底图调整阈值。

方法比较：Band Math、Classification 和阈值分割怎么选

方法比较：在 ENVI 中提取水体，不只有 Band Math 一种方式。Band Math 主要负责计算 NDWI 指数，后续是否要生成水体矢量或水体掩膜，还需要根据任务选择其他方法。

方法 适用场景 优点 局限
Band Math 计算 NDWI 生成水体指数、观察水体分布 公式清晰，操作简单，适合教学和快速分析 只生成指数，不自动得到最终水体边界
NDWI 阈值分割 把 NDWI 转成水体/非水体二值图 流程直观，便于面积统计 阈值需要按区域调整
监督分类 复杂地物分类、水体与阴影混杂区域 可以结合多个波段和样本 需要训练样本，结果依赖样本质量
面向对象分类 高分辨率影像水体提取 能利用形状、纹理和邻域信息 参数较多，学习成本更高
MNDWI 城市水体、建筑物干扰明显区域 对建筑物背景有时更稳健 需要短波红外波段，不适合所有传感器

如果你的目标是学习 ENVI Band Math 怎么用，建议先完成 NDWI 计算，再尝试阈值分割。如果你的目标是正式水体面积统计，应进一步做掩膜、阈值优化、斑块清理和精度验证。

检查清单：计算结果是否正确怎么判断

检查清单：ENVI 计算 NDWI 指数后，不要只看软件是否成功输出文件，还要检查结果是否合理。

  • 是否确认了影像数据源，例如 Landsat 8、Sentinel-2 或其他传感器。
  • 是否确认 Green 和 NIR 的具体波段编号。
  • Band Math 表达式是否为 float(b1 - b2) / float(b1 + b2)。
  • b1 是否绑定到绿光波段。
  • b2 是否绑定到近红外波段。
  • 输出结果是否为单波段浮点栅格。
  • 水体区域的 NDWI 值是否明显高于植被、裸地和建筑物。
  • 无数据区、黑边、裁剪外区域是否被排除。
  • 如果做面积统计，是否先确定合理阈值。
  • 如果做多期对比，是否使用同一数据源、同一预处理流程和同一空间范围。

一个简单的验证方法是：在影像上找几个明显水体点、植被点、裸地点和建筑物点，分别查询 NDWI 值。如果水体没有表现出相对高值，就要回头检查波段、公式和预处理。

FAQ：ENVI 计算 NDWI 和 Band Math 常见问题

FAQ：下面整理几个与 ENVI 怎么计算 NDWI 指数、Band Math 怎么用直接相关的问题。

1. ENVI Band Math 里可以直接写 Band3 和 Band5 吗？

通常不建议这样写。ENVI Band Math 更常见的做法是先写变量表达式，例如 float(b1 - b2) / float(b1 + b2)，然后在变量选择界面把 b1 指定为 Green，把 b2 指定为 NIR。这样更清晰，也不容易因为文件波段顺序不同而混乱。

2. NDWI 结果为什么有负值？

NDWI 本来就可能出现负值。NDWI 的理论范围通常在 -1 到 1 之间。非水体、植被、裸地和建筑物经常表现为较低值或负值。只要水体区域相对高于周边地物，结果通常是正常的。

3. NDWI 大于 0 就一定是水体吗？

不一定。NDWI 大于 0 可以作为初始经验阈值，但不同地区、不同影像、不同季节会有差异。山体阴影、建筑阴影、湿地植被和云阴影都可能干扰判断。正式项目中应结合样本验证和阈值调整。

4. Landsat 8 计算 NDWI 应该用哪两个波段？

Landsat 8 OLI 常规 NDWI 使用 Band 3 作为 Green，Band 5 作为 NIR，公式为 (Band3 - Band5) / (Band3 + Band5)。在 ENVI Band Math 中可以写成 float(b1 - b2) / float(b1 + b2)，然后把 b1 绑定 Band 3，把 b2 绑定 Band 5。

5. Sentinel-2 计算 NDWI 应该用哪两个波段？

Sentinel-2 常规 NDWI 通常使用 Band 3 作为 Green，Band 8 作为 NIR，公式为 (Band3 - Band8) / (Band3 + Band8)。如果不同分辨率波段参与计算，要注意重采样和空间分辨率一致性。

6. NDWI 和 MNDWI 有什么区别？

NDWI 常用 Green 和 NIR，适合一般水体增强。MNDWI 使用 Green 和 SWIR，常用于建筑物背景复杂的城市水体提取。两者公式不同，适用场景也不同。本文讲的是 ENVI 中用 Band Math 计算常规 NDWI。

7. ENVI 计算 NDWI 后怎么提取水体面积？

可以先对 NDWI 结果设置阈值，把水体和非水体分成二值栅格，然后进行栅格统计或转矢量后计算面积。面积统计前要确保影像坐标系单位合理，最好使用投影坐标系而不是经纬度坐标系。

结论：ENVI 计算 NDWI 的关键是波段和公式都要对

结论：ENVI 怎么计算 NDWI 指数并不复杂，核心流程就是：打开影像、确认 Green 和 NIR 波段、进入 Band Math、输入 float(b1 - b2) / float(b1 + b2)、绑定正确波段并输出结果。

真正需要注意的是三个细节：第一，不同传感器的波段编号不同；第二，Band Math 中 b1 和 b2 的变量绑定不能反；第三，NDWI 阈值不能机械套用，必须结合研究区和影像质量验证。

如果你只是学习 ENVI Band Math 怎么用，建议先用 Landsat 8 或 Sentinel-2 做一幅 NDWI 练习图；如果你要做正式水体提取，还应继续完成阈值分割、噪声清理、面积统计和精度验证。