2000-2023年全球每日归一化水指数NDWI数据集

1. NDWI

归一化水指数(NDWI,Normalized Difference Water Index)是通过对遥感影像中的特定波段进行归一化差值计算,以突出影像中的水体信息的一种指数。NDWI的通用计算公式为:

NDWI = ((G - IR)/(IR + G))。G - 来自绿波段的像素值;IR - 来自红外波段的像素值。

NDWI对植被冠层中液态水含量的变化非常敏感,通常用于监测植被的水分状况、检测水体以及研究水文过程。NDWI的值范围从 -1.0 到 1.0,正值通常表示水体或湿润的植被区域,而负值则代表干燥的土壤或非水体区域。该指数的计算基于近红外(NIR)波段和一个短波红外(SWIR)波段,通常为≈1.24μm的波段,当该波段不可用时,会选择最近的可用红外波段来替代。MODIS Terra 每日 NDWI 产品(MODIS/006/MOD09GA)是通过表面反射率合成数据生成的。NDWI的主要应用包括监测植被的水分含量、识别水体、以及研究湿地变化和洪水等水文应用。在数据技术上,NDWI波段的像素大小为463.313米,通过近红外和短波红外波段的组合,帮助有效检测植被中的水分变化及水体的分布情况。

2. 数据获取

本次分享的是部分年份的全球NDWI年数据,后台回复MOD09GANDWI获取网盘链接。

"gee"(Google Earth Engine)是一个强大的在线数据分析平台,常用于地理空间数据分析,包括遥感图像处理。要计算体面积并基于Sentinel数据应用大津(Otsu's method)二值化,然后制作20162023的变化折线图,你可以按照以下步骤编写JavaScript代码: ```javascript // 导入必要的库和数据集 var SentinelData = ee.ImageCollection('COPERNICUS/S2') var waterIndex = ee.ImageCollection('JRC/GSW1_3'); // 定义日期范围 var startDate = '2016-01-01'; var endDate = '2023-12-31'; // 获取指定时间范围内的Sentinel数据,并选择合适的波段 var filteredData = SentinelData.filterDate(startDate, endDate) .select(['B4', 'B8']); // 计算NDWI (归一化差分指数),这是一个常用来识别体的指标 var ndwi = filteredData.normalizedDifference(['B8', 'B4']); // 应用大津阈值分割(Otsu's thresholding) var otsuThreshold = ee.Algorithms.OtsuThresholding(nda); var binaryWaterMask = ndwi.gt(otsuThreshold); // 提取每个影像的体面积 var waterAreaSeries = filteredData.map(function(image) { var year = image.date().get('year'); return image.updateMask(binaryWaterMask).select(null, 'water').reduceRegion({ reducer: ee.Reducer.sum(), geometry: image.geometry(), // 使用影像几何进行计算 scale: 100, // 缩放因子 bestEffort: true, maxPixels: 1e13 // 处理大数据量 }).set('year', year); }); // 合并所有份的结果为一个FeatureCollection var waterAreaFC = ee.FeatureCollection(waterAreaSeries.flatten()); // 制作折线图 Map.addLayer(waterAreaFC, {color: 'blue'}, 'Water Area Over Time'); var chart = ui.Chart.image.seriesByYear( waterAreaFC, 'year', 'water', {xTitle: 'Year', yTitle: 'Water Area (sq. km)', seriesTitle: 'Sentinel Data - Water Body Evolution'} ); // 输出图表到UI或导出成文件 chart.draw(); ``` 这只是一个基本示例,实际操作可能会根据具体的需求(如处理不同分辨率的数据、数据预处理等)有所调整。运行这段代码需要在Google Earth Engine平台上。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

GIS遥感数据处理应用

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值