Earth Engine下水体提取

本次内容为Landsat-8遥感影像中水体的提取，用到了New Water Index(NWI)指数，效果较好。

文献（丁凤，2009）中还花一定篇幅对大气校正进行阐述，意在体现其重要性。但个人认为如果只是对分类或提取，大气校正其实是鸡肋的（真香警告），所以本次辨证地对TOA（大气层顶反射率）计算结果NWI(TOA)和SR（地表反射率）计算结果NWI(SR)进行了对比：在案例区域（长江黄冈段）中，NWI(TOA)只要大于0就可以很好地提取出水体，而NWI(SR)则要自己调整阈值。但在后面的普适性验证中，发现NWI(TOA)>0这样的条件虽然很简单，实则并不具备普适性，容易将山体阴影也划分为水体，所以最后还是在案例区域找到了适合NWI(SR)的阈值，并在我国大部分区域（华东，华南，华北，东北，西北，西南，台湾北，台湾南）做了验证，发现-0.6这个阈值很适合Landsat-8地表反射率产品计算的NWI指数做水体提取，对于快速提取水体有很大帮助。有兴趣的可以自行验证。（注：Earth Engine中的SR产品需要除10000才是小于1的反射率值）

以东北地区验证为例，结果如下，蓝色即为水体：

不足：提取结果中包含水田或蓄水池，且对于单像元的水体也可以进行提取，所以提取结果存在碎斑化现象。

/** Extrction of Water* Using Landsat-8 Imagery** Author: Deserts Tsung, Nov 20，2018* Chengdu University of Information Technology*
*/// Function of Water Extraction
// Following New Water Index(NWI) was Proposed by DING Feng in 2009
var waterExtraction = function(inputImage) {var B1 = inputImage.select('B1'), B5 = inputImage.select('B5');var B6 = inputImage.select('B6'), B7 = inputImage.select('B7');var infrareds = B5.add(B6).add(B7);var NWI = B1.subtract(infrareds).divide(B1.add(infrareds)).rename('NWI');var waterPart = inputImage.updateMask(NWI.gte(-0.6));return waterPart;
};var sr = ee.Image('LANDSAT/LC08/C01/T1_SR/LC08_122039_20181026').divide(10000);
var wtrImg = waterExtraction(sr);
Map.addLayer(sr,{bands:['B5','B4','B3'],min:0,max:0.4},'sr');
Map.addLayer(wtrImg.select('B5'),{palette:'#00FFFF'},'water');// East P.R.China Verified
// South P.R.China Verified
// North P.R.China Verified
// North-east Area Verified
// South-west Area Verified
// North-west Area Verified
// North Taiwan Area Verified
// South Taiwan Area Verified

Earth Engine下水体提取相关推荐

Earth Engine下地表温度反演
首先感谢ENVI-IDL技术殿堂在ENVI平台下的指导,原文链接,现将其在GEE平台实现.不同的是,该博文称使用原始DN值计算NDVI对反演影响不大,但通过实践发现还是有差异,所以代码采用了经过大气校 ...
Earth Engine下长时间序列监测
巨流河是清代称呼辽河的名字,她是中国七大江河之一,辽宁百姓的母亲河.哑口海位于台湾南端,是鹅銮鼻灯塔下的一泓湾流,据说汹涌海浪冲击到此,声消音灭. 本次将以哑口海周边区域为对象,研究该海岸近30年来的 ...
Google Earth Engine（GEE）农作物种植结构提取
目录写在前面 1.构建物候特征 2.构建光谱特征 3.将所有影像合并为一幅影像 4.构建随机森林算法进行分类 5.算法的存储 6.面积统计写在前面前段时间因为考研的原因一直没能更新,已经完成了农 ...
Google Earth Engine（GEE）——用两种方法计算NDWI水域面积提取（Landsat 8）
使用 Landsat 8 使用指标 NDWI 使用两种方法计算水域面积 ,这次我们随便选择了陕西附近的一个部分作为我们的研究区, 我们所需用到的函数: ee.Image(args) An object ...
Google Earth Engine（GEE）——利用JRC全球水体数据求湖泊面积以2000—2020年青海湖为例
问题: 之前写了一篇关于JRC水体面积分析和下载的文章,但是前两天有同学问到,如何进行青海湖面积水体的分析,逐月还是逐年的我记不清楚了,但其实都一样,只是在时间筛选的过程中有所不同而已.上次的链接: ...
Google Earth Engine（GEE）——Sentinel-2 Level-2A数据进行水域面积提取（宁夏为例）
本次我们采用哨兵2号卫星数据进行水体提取,提取出的的水体结果显示效果是非常精确的,先来看一下数据:这次我们按照塞上江南来看一眼最后提取的影像结果是非常不错的 Sentinel-2 MSI: Multi ...
GEE_API Docs_Tutorials_1.编程基础和Earth Engine API入门
API Docs_Tutorials_1.编程基础和Earth Engine API入门一.Introduction to JavaScript for Earth Engine(JavaScrip ...
论文笔记（五）FWENet：基于SAR图像的洪水水体提取深度卷积神经网络（CVPR）
FWENet: a deep convolutional neural network for flood water body extraction based on SAR images 作者:J ...
Google Earth Engine学习笔记（一）
基于GEE平台的Landsat8.Sentinel2.MODIS"去云"处理及FUI水色指数运算文章目录基于GEE平台的Landsat8.Sentinel2.MODIS&quo ...

Earth Engine下水体提取

Earth Engine下水体提取相关推荐

最新文章

热门文章