使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]我最近想通过GEE直接统计长时间序列下长江下游流域的水体面积变化情况,如果通过传统做法很复杂,于是想到了使用GEE。下面是提取水体的效果图,时间是2000年的,你也可以设置显示2000-2019年中任意一年的水体数据。代码链接:https://code.earthengine.google.com/2440b9511ba0c1cefaf926c7c47e5ea2读者可以先通过代码看看效果,下面我说一下最主要的思路:1.数据源的选择;2.数据的时间、地点筛选;3.水体的标准设置;4

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

我最近想通过GEE直接统计长时间序列下长江下游流域的水体面积变化情况,如果通过传统做法很复杂,于是想到了使用GEE。下面是提取水体的效果图,时间是2000年的,你也可以设置显示2000-2019年中任意一年的水体数据

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

代码链接:

https://code.earthengine.google.com/de4c42e80c432a909a4167b2b6333e6e使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]https://code.earthengine.google.com/de4c42e80c432a909a4167b2b6333e6e

读者可以先通过代码看看效果,下面我说一下最主要的思路:

1.数据源的选择;

2.数据的时间、地点筛选;

3.水体的标准设置;

4.影像数据导出;

5.统计长江水体面积;

6.时间迭代函数。

下面是每个步骤的做法:

1.数据源的选择:

水体数据源是JRC Monthly Water History, v1.1,这个数据库拥有1987年到现在的全球水体分布图,而这个数据库是根据LANDSAT数据提取的。

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

这个数据只有一个波段,3个值:0,1,2.。其中0代表没有数据;1代表有数据,不是水体;2代表有数据,是水体。

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

我们直接通过GEE的import功能,引进这个数据库:

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

2.数据的时间、地点筛选:

year是研究的年份,yantze_down是我上传的长江中下游的shp图;

//设置需要提取的区域,由于是上传的shp文件,需要转为geometry的格式
    var yantze_down_region = yantze_down.geometry();

    //设置需要提取的年份
    var startDate = ee.Date.fromYMD(year, 1, 1);
    var endDate = ee.Date.fromYMD(year, 12, 31);
    
    //筛选JRC水体数据
    var myjrc = jrc.filterBounds(yantze_down_region).filterDate(startDate, endDate);

3.水体的标准设置;

在这一步时一定要注意,由于JRC数据一年有12景,我们需要根据12景影像来统计一年的合成影像,而不是单单一景。所以合成方法是什么?我考虑到长江下游某些地方绘出翔一些低洼地暂时性积水,但这并不能算作湖泊。而我的目的是稳定性的湖泊、河流。所以我的合成方法是:

某个像素有7个月是水,则该像素为水体。


    //在每个月份的影像中添加一个obs属性的波段,一个像素如果有数据,则为1,没有数据则为0
    myjrc = myjrc.map(function(img){
      var obs = img.gt(0);
      return img.addBands(obs.rename('obs').set('system:time_start', img.get('system:time_start')));
    });

    //在每个月份的影像中添加一个onlywater属性的波段,一个像素如果有水则为1,没有水则为0
    myjrc = myjrc.map(function(img){
      var water = img.select('water').eq(2);
      return img.addBands(water.rename('onlywater').set('system:time_start', img.get('system:time_start')));
    });
    
    //计算每个像素点在一年12景影像中, 有数据的次数
    var totalObs = ee.ImageCollection(myjrc.select('obs')).sum().toFloat();
    
    //计算每个像素点在一年12景影像中, 有水的次数
    var totalWater = ee.ImageCollection(myjrc.select('onlywater')).sum().toFloat();
    
    //统计每个像素点在一年中有水的比例
    var floodfreq = totalWater.divide(totalObs).multiply(100);
    
    //删除没有值的像素
    var myMask = floodfreq.eq(0).not();
    floodfreq = floodfreq.updateMask(myMask);
    
    var viz = {min:0, max:50, palette: ['blue', 'white', 'green']};
    var floodfreq1=floodfreq.clip(yantze_down_region);
    var year_folder=year+"folder_gte";
    
    //如果某个像素一年有7个月有水,则为水体
    var gte60=floodfreq1.gte(60)

4.影像数据导出;

这个就太简单了,直接导出,文件名字为年份,不用详细说。

    
    //导出影像
    Export.image.toDrive({
      image: gte60,
      region: yantze_down_region,
      // fileDimensions:2560,
      scale: 30,
      maxPixels : 1e13,
      folder:year_folder,
      description:year_folder});

Task里面就可以有每一年长江下游区域的水体TIFF图,点击下载

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

5.统计长江水体面积;

这个用reduce函数,直接统计水体面积,也很简单,不详细说。

    //计算计算长江下游水体面积
    var stats2 = gte60.reduceRegion({
      reducer: ee.Reducer.sum(),
      geometry: yantze_down_region,
      scale: 30,
      maxPixels: 1E13
    });

    print(year_folder);
    print(stats2);

控制台里面输出每一年的水体像元个数,每一个像元的面积为30m*30m:

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

6.时间迭代函数。

因为我想统计长时间序列的水体数据,所以进行一个迭代。在这之前,先将前5个步骤封装,传入参数为年份,方便我们在时间迭代函数中进行调用:

//获取哪一年的,如果你想获取2000年到2019年,将条件改为i<20
for(var i=0;i<2;i++){
  if (i<10){  var year='200'+i;}
  if (i>10||i==10){  var year='20'+i;}
  var yearn = parseInt(JSON.parse(year));
  get_yearly_water(yearn);
}

提取水体的过程也就完成了,你按照我这个办法,可以提取全球任意地点长时间序列的水体数据,希望对你有帮助。

可以前往“地信遥感数据汇”​​​​​​​获取更多数据。
https://www.gisrsdata.com/

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

 

 VX:kitmyfaceplease2;欢迎关注公众号:锐多宝的地理空间;

使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据[通俗易懂]

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/193097.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • JetBrains PyCharm 2018.3.1 x64 激活码_通用破解码

    JetBrains PyCharm 2018.3.1 x64 激活码_通用破解码,https://javaforall.cn/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

  • 傅里叶变换(一)——认识傅里叶变换

    傅里叶变换(一)——认识傅里叶变换注:本文为博主参考书籍和他人文章并加上自己的理解所编,作为学习笔记使用并将其分享出去供大家学习。若涉及到引用您的文章内容请评论区告知!如有错误欢迎指正! 参考文章:https://zhuanlan.zhihu.com/p/19763358p.s.本文无论是cos还是sin,都统一用“正弦波”(SineWave)一词来代表简谐波。一、什么是傅里叶变换 时域及频域    在讲…

  • 精美的液晶数字字体素材[通俗易懂]

    精美的液晶数字字体素材[通俗易懂]液晶数字应该比较常见,那么液晶数字字体的应用也是相对广泛了,可以运用于一切需要液晶显示屏上的数字字体显示。对于这样一种有着广泛的应用数字字体,选择使用哪款液晶数字字体也是一个很重要的问题啦!为此,特意为大家收集了几款液晶数字字体供大家选择,喜欢的朋友赶紧收藏起来吧!  DS-Digital字体是一款比较常规的液晶数字字体,这款字体的仅支持数字和大写字母输入,字体端正,结构完整,整体视觉呈现效果…

    2022年10月25日
  • Unity零基础到入门 ☀️| Unity从青铜 到 王者!只差这篇让你学会Unity中最重要的部分——脚本组件✨

    Unity零基础到入门 ☀️| Unity从青铜 到 王者!只差这篇让你学会Unity中最重要的部分——脚本组件✨从倔强青铜到荣耀王者!只差这篇让你学会Unity中最重要的部分——脚本组件✨。本篇内容对Unity脚本做了一个超级详细的介绍,从开发工具到具体实例,API等都有介绍。

  • java缓存技术的介绍

    java缓存技术的介绍一、什么是缓存1、Cache是高速缓冲存储器一种特殊的存储器子系统,其中复制了频繁使用的数据以利于快速访问2、凡是位于速度相差较大的两种硬件/软件之间的,用于协调两者数据传输速度差异的结构,均可称之为Cache二、缓存的分类1、基于web应用的系统架构图2、在系统架构的不同层级之间,为了加快访问速度,都可以存在缓存操作系统磁盘缓存->减少磁盘机械操作

  • 本地数据库同步到云主机上

    本地数据库同步到云主机上同步前的准备:首先你本地跟云主机上都要有数据库、可视化的辅助工具(我用的NavicatPremium,其他的也都一个道理),这里靠的就是这个NavicatPremium工具1.首先在云主机上创建一个链接,建一个数据库,最好是与本地数据库同名2.在本地新建一个连接,可以点击下边的链接测试,测试一下看看是否能连接成功3.找到工具栏里的:工具->数据传输4.经过上述三步你最起码有了两个连接

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号