GEE学习笔记 八十二:湖北旱情分析(2019年8月-9月)

news2024/10/5 15:23:43

2019年8月1日至9月15日,湖北平均雨量为近53年来同期最少,仅69.4毫米。尤其是中东部地区降水异常偏少,一个半月里大多数地方的累计降水量均不足70毫米,其中洪湖、英山降水仅4毫米。鄂西东部一直到鄂东地区的降水量均比常年同期偏少6成~9成,黄冈、武汉、黄石、咸宁等地偏少8成以上。8月以来,湖北中东部地区始终没有出现有效降水,加上今年高温天气持续时间长,进一步加剧土壤水分蒸发,导致该省中东部地区的气象干旱持续发展加重。

来源:

百度安全验证

我们接下来就用GEE结合气象数据来做一个简单的数据分析,看看目前湖北干旱程度到底如何?


一、寻找合适的数据源

    这里用到的主要是气象数据,我这里使用的是NOAA的CFSV2这个6小时的气象数据:CFSV2: NCEP Climate Forecast System Version 2, 6-Hourly Products(https://developers.google.com/earth-engine/datasets/catalog/NOAA_CFSV2_FOR6H),主要用到的波段是降雨量波段和气温波段。

二、加载处理数据

    处理降雨数据和气温数据有些不太一样,降雨数据需要计算在时间段内的降雨总量,而气温我们寻找的是最高气温。气温计算比较简单我就不做代码展示,这里我简单展示一下如何计算降雨总量。

  1. //imgCol就是影像集合,其中pre为降雨量波段

  2. var preImgCol = imgCol.select("pre");

  3. //使用join将影像集合合并为每日数据

  4. var join = ee.Join.saveAll("matches");

  5. var filter = ee.Filter.equals({

  6.   leftField: "date",

  7.   rightField: "date"

  8. });

  9. var joinImgs = join.apply(preImgCol.distinct("date"), preImgCol, filter);

  10. var preImgs = joinImgs.map(function(image) {

  11.   var _imgList = ee.List(image.get("matches"));

  12.   var _tempCol = ee.ImageCollection.fromImages(_imgList);

  13.   //由于是6小时,所以时间需要计算6*60*60

  14.   var _dayImg = _tempCol.sum().clip(roi).multiply(6 * 60 * 60);

  15.   var _date = image.get("date");

  16.   _dayImg = _dayImg.set("system:time_start", ee.Date.parse("yyyy-MM-dd", _date).millis());

  17.   return _dayImg;

  18. });

  19. preImgs = ee.ImageCollection(preImgs);

三、分析展示结果

下面展示的是2019年8月1日到2019年9月15日湖北全省每日最高温度和降雨变化情况

每日最高温度,单位摄氏度

每日降雨量

总体状况

 

    通过上面图表可以了解到这一个半月湖北一直是高温少雨的情况。接下来可以看一下空间分布图,可以明显看到在湖北东部,高温少雨干旱比较严重。

(1)最高温度分布图:

(2)降雨总量分布图

进一步我们统计一下各个市的状况:

(1)降雨量

(2)气温

四、对比历史数据

    其中蓝色线是实际的情况,红色线是拟合的趋势线,使用的是GEE中linearFit()来做线下回归计算得到的结果,相关代码如下:

  1. var dataList = yearImgs.reduceColumns(ee.Reducer.toList(2), ["year", "tem"])

  2.                        .get("list");

  3. dataList = ee.List(dataList);

  4. var result = dataList.reduce(ee.Reducer.linearFit());

  5. print("linear fit result", result);

  6. var offset = ee.Number(ee.Dictionary(result).get("offset"));

  7. var scale = ee.Number(ee.Dictionary(result).get("scale"));

  8. //计算预测数据

  9. dataList = dataList.map(function(data) {

  10.   data = ee.List(data);

  11.   var predict = ee.Number(data.get(0)).multiply(scale).add(offset);

  12.   data = data.add(predict);

  13.   return data;

  14. });

(1)对比40多年同时期的降雨数据

 

其中趋势拟合结果如下,非常明显可以看到降雨量总体趋势是在下降,也就是降雨量越来越少。

(2)最高气温

其中趋势拟合结果如下,气温总体趋势是在上升,最高气温是越来越高,天气越来越热。

上面所有的是总体情况对比,下面是各个主要的市对比历史情况。

    这几个主要的市也是目前干旱最严重的地方,从对比数据结果就可以看出降雨量同比下降了8成左右,同历史相比也下降了8成左右。

五、在进一步分析一下

    这个分析有点牵强,不过还是将分析结果直接贴出来,是关于湖北省雨季总的降雨量的一个统计分析,时间切割点在大坝修筑前后,可以看到总体降雨量少了将近100mm,具体为啥会这样我也不知道原因。

来源请引用:地理遥感生态网科学数据注册与出版系统 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/350934.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

IP协议相关特性

日升时奋斗,日落时自省 目录 1、地址管理 1.1、子网掩码 2、路由选择 针对网络层的IP协议理解 TCP/IP协议栈TCPIP 此处详细解析IPv4协议(v4版本) 4位版本:此处的取值只有两个v4 ,v6 4位首部长度 :描…

【WPF】WindowChrome 自定义窗口完美实现

WindowChrome 自定义窗口完美实现简介效果图自定义最小化、最大化、关闭按钮布局实现结语简介 Microsoft官网关于 WindowChome 的介绍 截取Microsoft文章的一段话:   若要在保留其标准功能时自定义窗口,可以使用该 WindowChrome 类。 该 WindowChrome…

多元回归分析 | CNN-BiLSTM卷积双向长短期记忆神经网络多输入单输出预测(Matlab完整程序)

多元回归分析 | CNN-BiLSTM卷积双向长短期记忆神经网络多输入单输出预测(Matlab完整程序) 目录 多元回归分析 | CNN-BiLSTM卷积双向长短期记忆神经网络多输入单输出预测(Matlab完整程序)预测结果评价指标基本介绍程序设计参考资料预测结果 评价指标 训练结束: 已完成最大轮…

新建一个完整的react项目和完善初始项目

一:新建一个完整的react项目 1.环境准备 目前我的环境是 node:16.17.1 npm: 8.15.0 查看环境:1):打开命令提示符工具,利用node -v和npm -v 查看一下自己的环境,如果觉得重新卸载、安装node比较…

用记事本实现“HelloWorld”输出

一、在任意文件夹中创建一个新的文本文档文件并写入以下代码 public class Hello{public static void main (String[] args){System.out.print("Hello,World!");} } 二、修改文件名称及文件类型为 Hello.java 特别注意:文件命名必须与代码中类的名称相同…

为什么伟大的产品只专注做一件事

uber 不允许你预订出租车。亚马逊一开始只是卖书。谷歌只是一个搜索引擎。麦当劳没有餐具。不知为什么,我们仍然相信一个产品要想成功,它必须做很多事情。这通常发生在两种情况下:当新产品试图让市场相信它们是值得的,或者当公司提…

儿童蜡笔出口美国CPC认证CPSIA+ASTM963测试

蜡笔是将颜料掺在蜡里制成的笔,可有数十种颜色,画画用。蜡笔没有渗透性,是靠附着力固定在画面上,不适宜用过于光滑的纸、板,亦不能通过色彩的反复叠加求得复合色。它是儿童学习色彩画的理想工具,一些画家用…

【Spark分布式内存计算框架——Spark SQL】5. DataFrame(下)

3.3 Row DataFrame中每条数据封装在Row中,Row表示每行数据,具体哪些字段位置,获取DataFrame中第一条数据。 如何构建Row对象:要么是传递value,要么传递Seq,官方实例代码: import org.apache.…

百分点科技宣布接入百度文心一言能力

2月16日, 百分点科技宣布成为百度文心一言(英文名:ERNIE Bot)首批生态合作伙伴。后续,百分点科技将全面体验并接入文心一言的能力。百分点科技专注于数据科学理论和技术的创新实践,长期坚持基础技术和行业应…

初始QML

Qt Quick的介绍 : Qt Quick是QML的标准类型和功能库。它包括视觉类型,交互类型,动画,模型和视图,粒子效果和着色器效果。QML 应用程序开发人员可以通过单个导入语句访问所有这些功能,简单来说Qt Quick是一…

使用 Hashnode API、Typescript 和 GraphQL 将博客文章添加到您的 React 站点

在本文中,我们将:使用 Next.js 引导一个 React.js Typescript 项目。设置 Apollo GraphQL 客户端并将其集成到我们的项目中。设置 GraphQL Codegen 以生成我们可以在整个应用程序中使用的类型、类型安全查询和自定义挂钩。创建一个索引页面,其…

Failed at the node-sass@4.14.1 postinstall script

vue项目启动,安装node14.18.0版本,构建时报错: Failed at the node-sass4.14.1 postinstall script 其实在构建过程中,还出现了其他组件的各种报错,最后反思了一下,觉得是nodeJs的版本问题,最…

BIM技巧 | Revit中如何给房间填充颜色?就5步

大家在Revit平面创建好房间后,有没有觉得各房间因为没有着色而区分不明显、视觉效果一般呢? 一、今天就教给大家如何给房间填充上颜色。 01 第一步 首先,将各个房间创建好; 02 第二步 在【建筑】-【房间和面积】单击下拉菜单…

深入理解vue2.x中Object.defineproperty()和vue3.x中Proxy

前言 vue2.x中数据的双向绑定主要通过Object.defineproperty()方法实现,data中的数据改变通过Object.defineProperty()对属性设置set属性,获取通过get属性,Object.defineProperty的作用就是劫持一个对象的属性,通常我们对属性的getter和sett…

魔兽世界私服架设教程——如何搭建魔兽世界私服

TrinityCore是一个魔兽世界服务端模拟器,我们可以通过TrinityCore来学习大型网络游戏服务端的编写,从中汲取营养来编写我们自己的游戏。一、前期准备工作CPU需要支持SSE2指令集Boost版本大于等于1.59.0MySQL数据库版本大于等于5.1.0OpenSSL版本为1.0.xCM…

基于机器学习LSTM的古代汉语切分标注算法及语料库研究 完整代码+数据+论文

完整代码:https://download.csdn.net/download/qq_38735017/87382302摘 要近年来,深度学习的浪潮渗透在科研和生活领域的方方面面,本文主要研究深度学习在自然语言处理,尤其是古汉语自然语言处理方面的应用。本文旨在利用计算机帮…

C#中GDI+的矩形功能扩展

原文出处:https://haigear.blog.csdn.net/article/details/129060020 GDI发展到GDI绘制函数中的参数往往都有矩形这个参数(除绘制直线和路径),所以我们用好了矩形绘图就容易多了。 一、中心定位绘制图形 但当我们绘制一个图形时…

Towards Adversarial Attack on Vision-Language Pre-training Models

摘要虽然视觉-语言预训练模型(VLP)在各种视觉-语言(VL)任务上表现出革命性的改进,但关于其对抗鲁棒性的研究在很大程度上仍未被探索。本文研究了常用VLP模型和VL任务的对抗性攻击。首先,我们分析了不同设置下对抗性攻击的性能。通过研究不同扰动对象和攻…

HHDESK图片管理——批量重命名及递归搜索

HHDESK作为一款国产桌面软件,考虑到国人的操作及阅读习惯。因此我们开发了一些有意义的新功能,比如今天要介绍的图片批量重命名及递归搜索功能 1.图片批量重命名功能 网上下载的图片名称大多杂乱无章,一眼望去毫无头绪。 而windows自带的…

第41天|LeetCode198. 打家劫舍、LeetCode213. 打家劫舍II、LeetCode337. 打家劫舍III

1.题目链接:198. 打家劫舍 题目描述: 你是一个专业的小偷,计划偷窃沿街的房屋。每间房内都藏有一定的现金,影响你偷窃的唯一制约因素就是相邻的房屋装有相互连通的防盗系统,如果两间相邻的房屋在同一晚上被小偷闯入&a…