使用google earth engine(GEE)提取亚马逊每年森林火灾区域
2019年亚马逊大火还历历在目,作为学遥感的学生,我们可以使用遥感影像来获得着火区。但是使用传统遥感太过于麻烦,我们可以使用GEE,1分钟获得亚马逊10年的着火区域,并统计面积与导出着火区。
代码主要步骤:
1.数据选择与预处理;
数据源选择modis的火点产品数据,该数据是modis的500m月合成数据。
研究区是通过roi工具在GEE上面绘制的,包括了亚马逊森林。
var temp=MCD64.filterDate(start_date,end_date).mosaic().clip(SA)var burn=temp.gte(0)Map.addLayer(burn)
2.下载着火区域;
由于下载区域比较大,需要把maxpixel设置大一点。
Export.image.toDrive({image:burn,maxPixels:1e13,scale: 500,region :SA,fileNamePrefix:NAME,folder:NAME,skipEmptyTiles:true,})
3.统计着火区域面积;
在GEE里面没有直接统计面积的函数,一般有两个办法:
第一种方法是计算每个像素的面积,再统计所有像素的面积;
第二种方法是统计所有像素的个数,再乘以单个像素的面积;
我才用的是第一种方法:
//计算单个像素面积var AreaImage = burn.multiply(ee.Image.pixelArea());//统计面积var AreaNum=AreaImage.reduceRegion({reducer: ee.Reducer.sum(),geometry: SA,scale: 500,maxPixels: 1e14})var BurnedArea=AreaNum.get('BurnDate')print(NAME+"这一年的面积"+"(单位:平方米)")print(BurnedArea)
4.嵌套一个时间循环函数。
因为我们想嵌套一个时间循环,得到2010到2020年的亚马逊森林着火情况.get_monthly_burn是把前3个步骤封装好,并名为get_monthly_burn函数,方便直接调用。
//日期循环
for(var i=0;i<1;i++){var start_date='201'+i+'-1-01';var end_date= '201'+i+'-12-31';get_monthly_water(start_date,end_date) }
根据这个脚本,我们就能得到亚马逊这10年的火灾情况,下表就是统计的着火面积:
时间 | 火灾面积(万平方千米) | 时间 | 火灾面积(万平方千米) | |
2000 | 0.30 | 2019.1 | 0.74 | |
2001 | 19.43 | 2019.2 | 0.72 | |
2002 | 30.69 | 2019.3 | 0.73 | |
2003 | 22.41 | 2019.4 | 0.24 | |
2004 | 28.54 | 2019.5 | 2.71 | |
2005 | 31.20 | 2019.6 | 0.68 | |
2006 | 19.63 | 2019.7 | 1.60 | |
2007 | 44.87 | 2019.8 | 7.28 | |
2008 | 20.46 | 2019.9 | 10.28 | |
2009 | 12.53 | 2019.1 | 3.26 | |
2010 | 47.07 | 2019.11 | 1.40 | |
2011 | 19.27 | 2019.12 | 0.34 | |
2012 | 28.26 | |||
2013 | 13.76 | |||
2014 | 18.18 | |||
2015 | 23.52 | |||
2016 | 21.13 | |||
2017 | 26.30 | |||
2018 | 12.07 | |||
2019 | 27.17 |
可以看出其实亚马逊森林火灾并没有想象中严重,历史上的2007年的火灾更严重。如此引人关注,可能是因为2019年9月份的着火面积太大了。
完整的代码可进入GEE里面查看:https://code.earthengine.google.com/191ce5cf316191bd1878ecb0e9af7fef
可以前往“地信遥感数据汇”获取更多数据。
https://www.gisrsdata.com/
VX:kitmyfaceplease2;欢迎关注公众号:锐多宝的地理空间;
使用google earth engine(GEE)提取亚马逊每年森林火灾区域相关推荐
- 使用google earth engine(GEE)提取2000年到2019年长江下游水体(河流、湖泊)数据
我最近想通过GEE直接统计长时间序列下长江下游流域的水体面积变化情况,如果通过传统做法很复杂,于是想到了使用GEE.下面是提取水体的效果图,时间是2000年的,你也可以设置显示2000-2019年中任 ...
- Google Earth Engine(GEE)——User memory limit exceeded(2)
上一次我们已经知道如何去进行避免这种错误的发生,有关详细内容,如果单单只是解决这个问题我们用到的是limit 和 first,上一次的博客在这里: (207条消息) Google Earth Engi ...
- Google Earth Engine(GEE) 01-中输入提示快捷键Ctrl+space无法使用的问题
Google Earth Engine(GEE) 01-中输入提示快捷键Ctrl+space无法使用的问题 GEE中 Ctrl+space组合键用于代码输入快捷提示,能够提高编码的准确度和速度,但是, ...
- Google Earth Engine(GEE)批量下载代码(以 NDVI数据为例)
下载数据先准备工作(具体细节都能查到):科学上网 谷歌邮箱,谷歌邮箱注册GEE账号. 一:导入需要下载边界shp文件. 标题 找到自己的shp文件,导入除了sbx文件的所有文件. 导入成功 命名ass ...
- 使用Google Earth Engine (GEE)实现MODIS数据批量下载
使用Google Earth Engine GEE实现MODIS数据批量下载 前言 下载数据代码 批量执行run任务 关注公众号,分享GIS知识.ArcGIS教程.SCI论文与科研日常等 前言 上图是 ...
- 基于google earth engine(GEE)下载研究区域影像
基于google earth engine(GEE)下载研究区域影像 当研究需要Landsat数据时,我们可以通过USGS官网或者地理空间数据云平台下载.由于地理空间数据云目前无法下载到较新的数据,可 ...
- 对口型应用WOMBO入驻华为AppGallery;亚马逊云科技中国区域上线新数据传输设备 | 全球TMT...
大公司动向 对口型应用WOMBO上个月入驻华为AppGallery,为华为全球用户带来AI对口型应用.WOMBO可以使用AI快速将一张图片变成搞笑和有趣的配音视频.只需简单三步,华为设备用户即可创建亲 ...
- 使用 Amazon Pricing Calculator 在亚马逊云科技中国区域估算云使用成本
亚马逊云科技为上百种云服务提供了按实际使用量付费的基本定价模式,并在此基础上支持多样的成本节省付费方式. 您可以从亚马逊云科技中国区域定价页面,了解亚马逊云科技各项产品和服务的定价策略.对于不管是从未 ...
- Google Earth Engine(GEE)——可视化动态图
代码: var geometry = /* color: #d63000 *//* shown: false *//* displayProperties: [{"type": & ...
- Google Earth Engine (GEE) ——卫星影像的监督分类(svm)
问题 GEE 提供哪些机器学习技术? 如何对卫星图像进行监督分类? 如何评估分类器的准确性? 如何手动创建自己的几何图形? 目标 练习查找无云图像和使用手绘几何导入 学习训练和应用分类算法所需的基本功 ...
最新文章
- BGP中IBGP和EBGP的区别和联系
- Linux vim打开文件的四种方式
- 输出最大值MXNet实现
- c语言大数相加oj,郑州轻工业大学oj题解(c语言)论如何正确的提高正确率:水题合集(四)...
- 已知法向量 求投影_MIT—线性代数笔记15 子空间投影
- vc如何打开plt图像_opencv_python从zero到hero————之图像基本操作01
- 路由器接口为啥不能配置为任意播地址_计算机三级网络技术(7):路由器配置及使用(下)...
- python全栈薪资这么高,你知道大厂面试都问什么吗?
- 从0开发《工程测绘大师》小程序之6°带中央子午线计算篇(八)
- 汤姆猫代码python_用树莓派实现会说话的汤姆猫
- 如何用ChatGPT做一门课?(包含大纲、脚本、PPT文本)
- 计算机常用程序在DOS中的英文名
- socket缓冲区/sk_buffer/滑动窗口关系
- STM32蓝牙控制循迹避障小车源代码——3.舵机、超声波测距模块
- 安科瑞智能照明控制系统在商业办公楼的应用
- 机械设计基础B【3】凸轮
- Oracle数据库第二课——使用PL/SQL进行初步编程,了解PL/SQL的相关知识
- 百度竞价账户选对关键词究竟有多重要?如何进行关键词筛选
- 【JESD204系列】六、加解扰模块的设计原理
- 我是如何从写不出来,到完成二十万字书稿的?
热门文章
- dataworks 生成表血缘依赖
- python编辑函数求n阶乘fact_python中fact函数的用法
- mac磁盘清理的方法大全
- html5数学公式编辑器,Daum Equation Editor:数学公式编辑器
- 用计算机程序方式打开,怎么使用运行命令符打开电脑程序
- 安装dhcp服务器虚拟2012,windows server 2012 dhcp服务器安装
- java中的坦克大战游戏设计报告论文_毕业论文基于Java语言的坦克大战游戏设计...
- 【光斑PSF处理】激光光斑数据处理含matlab源码-含文档【0102期】
- 微博和微信公众号营销区别
- 纯web端实现二维码识别