基于 GEE Landsat 与 Sentinel 数据的归一化水体指数 NDWI 计算和水体提取

目录

[1 归一化水体指数(NDWI)](#1 归一化水体指数(NDWI))

[2 完整代码](#2 完整代码)

[3 运行结果](#3 运行结果)

[4 运行结果可视化对比](#4 运行结果可视化对比)


1 归一化水体指数(NDWI)

归一化水体指数(NDWI)是遥感影像分析中常用的指数,主要用于识别和提取水体。NDWI通过比较可见光波段(绿波段)与近红外波段(NIR)的反射率差异,突出水体与非水体区域的区别。水体在近红外波段的反射率较低,而在绿波段的反射率较高,因此NDWI能够有效地识别水域。该指数在水体监测、水资源管理以及环境变化研究中有着广泛应用。

2 完整代码

javascript 复制代码
// 定义研究区域
var studyArea = table;

/**
 * 计算NDWI指数,并将其作为新波段添加到输入图像中
 * @param {ee.Image} image - 输入的影像
 * @param {Array<string>} bands - 用于计算NDWI的波段名称数组
 * @returns {ee.Image} - 添加了NDWI波段的影像
 */
function calculateNDWI(image, bands) {
    var ndwi = image.normalizedDifference(bands).rename('NDWI');
    return image.addBands(ndwi);
}

/**
 * 处理Sentinel - 2影像数据,计算并返回指定时间和区域内的平均NDWI
 * @returns {ee.Image} - Sentinel - 2数据的平均NDWI影像
 */
function processSentinel2() {
    return ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED')
       .filterDate('2023-03-01', '2023-05-31')
       .filterBounds(studyArea)
       .map(function (image) {
            // 云掩码处理
            var cloudMask = image.select('QA60')
               .bitwiseAnd(1 << 10).eq(0)
               .and(image.select('QA60').bitwiseAnd(1 << 11).eq(0));
            return image.updateMask(cloudMask);
        })
       .map(function (image) {
            return calculateNDWI(image, ['B3', 'B8']);
        })
       .select('NDWI')
       .mean()
       .clip(studyArea);
}

/**
 * 处理Landsat 8影像数据,计算并返回指定时间和区域内的平均NDWI
 * @returns {ee.Image} - Landsat 8数据的平均NDWI影像
 */
function processLandsat8() {
    return ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
       .filterDate('2023-03-01', '2023-05-31')
       .filterBounds(studyArea)
       .map(function (image) {
            // 光学波段和热红外波段的辐射定标
            var opticalBands = image.select(['SR_B.*'])
               .multiply(0.0000275).add(-0.2);
            var thermalBands = image.select(['ST_B.*'])
               .multiply(0.00341802).add(149.0);
            return image.addBands(opticalBands, null, true)
               .addBands(thermalBands, null, true);
        })
       .map(function (image) {
            return calculateNDWI(image, ['SR_B3', 'SR_B5']);
        })
       .select('NDWI')
       .mean()
       .clip(studyArea);
}

// 计算Sentinel - 2和Landsat 8数据的NDWI
var sentinel2NDWI = processSentinel2();
var landsatNDWI = processLandsat8();

// 设置地图中心并添加图层
Map.centerObject(studyArea, 10);
Map.addLayer(sentinel2NDWI, { min: -1, max: 1, palette: ['green', 'white', 'blue'] }, 'Sentinel-2 NDWI');
Map.addLayer(landsatNDWI, { min: -1, max: 1, palette: ['green', 'white', 'blue'] }, 'Landsat 8 NDWI');

// 导出Sentinel - 2数据中NDWI大于0的区域
Export.image.toDrive({
    image: sentinel2NDWI.gt(0).selfMask(),
    description: "Water_Bodies_Study_Area_2023_Sentinel2",
    scale: 30,
    region: studyArea,
    fileFormat: 'GeoTIFF',
    maxPixels: 1e9
});

// 导出Landsat 8数据中NDWI大于0的区域
Export.image.toDrive({
    image: landsatNDWI.gt(0).selfMask(),
    description: "Water_Bodies_Study_Area_2023_Landsat8",
    scale: 30,
    region: studyArea,
    fileFormat: 'GeoTIFF',
    maxPixels: 1e9
});

3 运行结果

点击RUN即可下载运行结果

4 运行结果可视化对比

点击RUN下载好运行结果栅格数据,然后利用ArcMap软件打开下载好的数据进行可视化对比,效果如下图所示。
Water_Landsat8
Water_Sentinel2

相关推荐
Johny_Zhao3 小时前
Docker + CentOS 部署 Zookeeper 集群 + Kubernetes Operator 自动化运维方案
linux·网络安全·docker·信息安全·zookeeper·kubernetes·云计算·系统运维
鸭鸭鸭进京赶烤4 小时前
大学专业科普 | 云计算、大数据
大数据·云计算
SH11HF8 小时前
小菜狗的云计算之旅,学习了解rsync+sersync实现数据实时同步(详细操作步骤)
学习·云计算
搞笑的秀儿11 小时前
信息新技术
大数据·人工智能·物联网·云计算·区块链
泊浮目18 小时前
未来数据库硬件-网络篇
数据库·架构·云计算
UI设计和前端开发从业者1 天前
UI前端大数据处理策略优化:基于云计算的数据存储与计算
前端·ui·云计算
我是小bā吖2 天前
阿里云服务网格ASM实践
网络·阿里云·云计算·服务发现
长征coder2 天前
AWS MySQL 读写分离配置指南
mysql·云计算·aws
Johny_Zhao2 天前
Docker 一键安装部署 JumpServer 堡垒机
linux·网络安全·信息安全·云计算·shell·jumpserver·ldap·yum源·系统运维
艾伦_耶格宇2 天前
【ACP】阿里云云计算高级运维工程师--ACP
运维·阿里云·云计算