【问题导向】GEE+Arcgis绘制广西不同时期的水资源分布图

本文主要是介绍【问题导向】GEE+Arcgis绘制广西不同时期的水资源分布图,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

背景

最近遇到的一个小作业,一开始是想按以前熟悉的方法来,后面因为嫌麻烦(懒)尝试了GEE快速做完了。开心
在这里插入图片描述

解决方案

作业是很简单的,就是找到某个地方的自然资源数据,然后绘制成图。
我有以下几个思路

1,DEM制作河流数据

在地理空间数据云下载DEM,然后arcgis中栅格系统的空间分析功能提取河流:首先对DEM数据进行填洼和削峰,然后确定每个栅格单元的水流方向,再根据各栅格单元的水流方向,计算出每个栅格单元的汇流能力,根据汇流能力采用阈值法确定河流。
优点:可以利用在线裁剪功能出一个省的DEM,节约处理的时间。
缺点:
可以下载到不同时期的DEM,不过DEM在几十年的尺度应该变化不大,很难看见不同时期降水的差异。
提取河流的过程比较繁琐,做过几次我就不想做了。

2,国家基础矢量数据

全国地理信息资源目录服务系统提供了1:25万的基础矢量数据,包括了河流、道路、建成区等多种数据。
不足:
1,很多地方只有一年的数据,无法做多年的比较。
2,仅广西一省就有28个shp数据,需要导入arcgis合并。
按之前的经验是把27个shp追加到第一个shp中在这里插入图片描述

3,统计年鉴

在各地的统计局找统计年鉴,然后手动整理数据,并输入shp数据中,利用arcgis出图。这种方法也要耗费较多的时间。
在这里插入图片描述

4,遥感数据分类数据

直接自己干:从地理空间数据云下载不同年份的原始遥感图像(Landsat、哨兵等),利用envi或者arcgis进行图像分类,一般有水、植被、耕地、建设用地、道路等类型。就可以进行分析土地变化。
在这里插入图片描述

5,GEE

GEE提供了一个快速处理数据的方案

GEE

Google Earth Engine是Google提供的对大量全球尺度地球科学资料(尤其是卫星数据)进行在线可视化计算和分析处理的云平台。该平台能够存取卫星图像和其他地球观测数据数据库中的资料并提供足够的运算能力对这些数据进行处理。Google Earth Engine上包含的数据集超过 200 个公共的数据集,超过500万张影像,每天的数据量增加大约 4000张影像,容量超过 5PB。
相比于ENVI等传统的处理影像工具,Google Earth Engine可以快速、批量处理数量“巨大”的影像。通过Google Earth Engine可以快速计算比如NDVI等植被指数,可以预测作物相关产量,监测旱情长势变化,监测全球森林变化等。

1,数据源

JRC Monthly Water History, v1.1,数据库拥有1987年到现在的全球水体分布图(根据LANDSAT数据提取)
在这里插入图片描述
数据只包含一个波段,3种值分别代表该像素点的三种情况: 0为无数据 ;1为有数据但不是水 ;2为水资源。
在这里插入图片描述
然后从中科院地理所下载广西的shp数据,上传到GEE。
在这里插入图片描述

在这里插入图片描述

2,提取数据

提取某一年,在广西境内的图像。若要修改研究区域,上传自己的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,提取水资源

对于一年来说,一共有12副图像,某个地方可能有几月降雨在低洼地形成了池塘,过几个月就没有水了,因此需要进一步处理,去除暂时性水体,只保留经常性水体。

添加两个波段,obs:若该像元点有数据,值为1否则为0。onlywater:若该像元点为水资源,值为1,否则为0。然后计算每个像元点,12个月的总obs值和总onlywater值。

统计每个像元点全年都有水百分比(SUM(onlywater)/SUM(obs))。若占比超过60%也就是该点7个月以上的时间都为水体,就认为该点为水体。

最后利用for循环不断调用函数,就能生成不同年份的图像数据。

//在每个月份的影像中添加一个obs属性的波段,一个像素如果有数据,则为1,没有数据则为0myjrc = myjrc.map(function(img){var obs = img.gt(0);//当第一个值大于第二个值时,返回1。return img.addBands(obs.rename('obs').set('system:time_start', img.get('system:time_start')));});//在每个月份的影像中添加一个onlywater属性的波段,一个像素如果有水则为1,没有水则为0myjrc = 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,导出图像

点击RUN,然后等待执行完毕;在自己账号的谷歌网盘内就能找到图片。下载到本地利用arcgis修饰出图。
在这里插入图片描述

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

代码:https://code.earthengine.google.com/d3a2c7d91d7d08d3fb5c18d86454434b.

参考:
https://blog.csdn.net/weixin_36396470/article/details/106899589.

新开通了本人的公众号,欢迎关注:燕南路GISer ,专注GIS干货分享,不定期更新。
主要兴趣:GIS、时空数据挖掘、python、机器学习深度学习
提问、求资源等都可在公众号后台留言
CSDN的部分内容会重写再搬迁到公众号,欢迎关注!
在这里插入图片描述

这篇关于【问题导向】GEE+Arcgis绘制广西不同时期的水资源分布图的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/362243

相关文章

怎样通过分析GC日志来定位Java进程的内存问题

《怎样通过分析GC日志来定位Java进程的内存问题》:本文主要介绍怎样通过分析GC日志来定位Java进程的内存问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、GC 日志基础配置1. 启用详细 GC 日志2. 不同收集器的日志格式二、关键指标与分析维度1.

Java 线程安全与 volatile与单例模式问题及解决方案

《Java线程安全与volatile与单例模式问题及解决方案》文章主要讲解线程安全问题的五个成因(调度随机、变量修改、非原子操作、内存可见性、指令重排序)及解决方案,强调使用volatile关键字... 目录什么是线程安全线程安全问题的产生与解决方案线程的调度是随机的多个线程对同一个变量进行修改线程的修改操

Redis出现中文乱码的问题及解决

《Redis出现中文乱码的问题及解决》:本文主要介绍Redis出现中文乱码的问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1. 问题的产生2China编程. 问题的解决redihttp://www.chinasem.cns数据进制问题的解决中文乱码问题解决总结

使用Python绘制3D堆叠条形图全解析

《使用Python绘制3D堆叠条形图全解析》在数据可视化的工具箱里,3D图表总能带来眼前一亮的效果,本文就来和大家聊聊如何使用Python实现绘制3D堆叠条形图,感兴趣的小伙伴可以了解下... 目录为什么选择 3D 堆叠条形图代码实现:从数据到 3D 世界的搭建核心代码逐行解析细节优化应用场景:3D 堆叠图

全面解析MySQL索引长度限制问题与解决方案

《全面解析MySQL索引长度限制问题与解决方案》MySQL对索引长度设限是为了保持高效的数据检索性能,这个限制不是MySQL的缺陷,而是数据库设计中的权衡结果,下面我们就来看看如何解决这一问题吧... 目录引言:为什么会有索引键长度问题?一、问题根源深度解析mysql索引长度限制原理实际场景示例二、五大解决

Springboot如何正确使用AOP问题

《Springboot如何正确使用AOP问题》:本文主要介绍Springboot如何正确使用AOP问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录​一、AOP概念二、切点表达式​execution表达式案例三、AOP通知四、springboot中使用AOP导出

Python中Tensorflow无法调用GPU问题的解决方法

《Python中Tensorflow无法调用GPU问题的解决方法》文章详解如何解决TensorFlow在Windows无法识别GPU的问题,需降级至2.10版本,安装匹配CUDA11.2和cuDNN... 当用以下代码查看GPU数量时,gpuspython返回的是一个空列表,说明tensorflow没有找到

解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘问题

《解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘问题》:本文主要介绍解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4... 目录未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘打开pom.XM

IDEA Maven提示:未解析的依赖项的问题及解决

《IDEAMaven提示:未解析的依赖项的问题及解决》:本文主要介绍IDEAMaven提示:未解析的依赖项的问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录IDEA Maven提示:未解析的依编程赖项例如总结IDEA Maven提示:未解析的依赖项例如

Redis分片集群、数据读写规则问题小结

《Redis分片集群、数据读写规则问题小结》本文介绍了Redis分片集群的原理,通过数据分片和哈希槽机制解决单机内存限制与写瓶颈问题,实现分布式存储和高并发处理,但存在通信开销大、维护复杂及对事务支持... 目录一、分片集群解android决的问题二、分片集群图解 分片集群特征如何解决的上述问题?(与哨兵模