用Python进行栅格数据的分区统计和批量提取


Posted in Python onMay 27, 2021

有时候我们会有这样的想法,就是针对某个区域的栅格数据,要提取它的平均值或者其他统计指标,比如在一个省内提取多年的降雨数据,最后分区域地计算一些统计值,或者从多个栅格数据中提取某个区域的数值形成一个序列。为了方便,画一个示意图看看,比如就像提取这个区域中的某一个市的区域,然后形成一个序列数据,这就可以使用rasterstats库了,此外的分区统计也可以用这个库

这个实验使用的数据格式分别是栅格(*.tif)和矢量(.shp),之后的分区统计操作和栅格数据的提取都是源于这两类数据。为了能使用上这个rasterstats库,选择了在google colab平台运行脚本,因为安装库实在是太方便了,在win上老是安装不上的,在google notebook立马就搞定了,而且可以把数据存储到谷歌云盘,直接在notebook中就是可以链接使用的

用Python进行栅格数据的分区统计和批量提取

那么现在就开始做测试,使用的数据就是左侧的栅格和矢量数据集
导入相关的模块

import geopandas as gpd
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import rasterio
import rasterstats
from rasterio.plot import show
# show()方法用来展示栅格图形
from rasterio.plot import show_hist
# 用来展示直方图
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter

使用geopandas和rasterio分别读取矢量和栅格数据

# 使用geopandas读取矢量数据
districts = gpd.read_file('/content/drive/MyDrive/Datashpraster/Data/Districts/districts.shp')

# 使用rasterio读取栅格数据,栅格数据和矢量数据的坐标投影需要一致
raster = rasterio.open('/content/drive/MyDrive/Datashpraster/Data/Rainfall Data Rasters/2020-4-1.tif')
# 把矢量数据和栅格数据绘制到一个axis上,这个axis不是坐标轴,而是图形
plt.rcParams['font.family'] = 'Times New Roman'
plt.rcParams['font.size'] = 20

fig, (ax1,ax2) = plt.subplots(1,2,figsize=(15,6))

show(raster, ax=ax1,title='Rainfall')
# 读取进来的矢量数据可以直接调用gpd的plot()方法绘制
districts.plot(ax=ax1, facecolor='None', edgecolor='red')
show_hist(raster,ax=ax2,title='hist')

plt.show()

先绘制一下结果看看

用Python进行栅格数据的分区统计和批量提取

读取栅格数据:

# 提取雨量栅格值到numpy数组
# 遵循GDAL规则从第一波段读取
rainfall_data = raster.read(1)
rainfall_data

用Python进行栅格数据的分区统计和批量提取

开始分区统计:

# 设置坐标变换信息
affine = raster.transform

# 准备开始进行空间分区计算
# 第一个参数是矢量分区,第二个是栅格,第三个是坐标变换信息,第四个是统计均值
avg_rallrain = rasterstats.zonal_stats(districts,rainfall_data,affine=affine,stats=['mean'],geojson_out=True)
# avg_rallrain

# 除了统计平均值之外,还有最大最小值那些

绘制一下,只是一个简单的图形而已

用Python进行栅格数据的分区统计和批量提取

当然第二部分更有意思,就是从多个分散的栅格数据中提取数据形成一个序列

,就是这些tif数据

用Python进行栅格数据的分区统计和批量提取

loop这些栅格数据集:

用Python进行栅格数据的分区统计和批量提取

获得提取到的结果,没错,就是这么一个序列数据,然后就是绘图了

用Python进行栅格数据的分区统计和批量提取

转换数据格式

# 将Date列转为时间型
data['Date'] = pd.to_datetime(data['Date'], infer_datetime_format=True)

# print(data)

data['Date'] = data['Date'].dt.date
print(data)

用Python进行栅格数据的分区统计和批量提取

绘图结果就是简单的图形而已

# 准备绘制图形
fig,(ax1,ax2)= plt.subplots(2,1,figsize=(18,6))
plt.rcParams['font.size'] = 15

data.plot(x='Date', y='Average_RF_Porto', ax=ax1, kind='bar', title='Avg_Rail_Porto')
data.plot(x='Date', y='Average_RF_Faro', ax=ax2, kind='bar', title='Avg_Rail_Faro',color='red')

#自动调整图形的分布
plt.tight_layout()
plt.show()

用Python进行栅格数据的分区统计和批量提取

结果就这样一个序列图,目的就是从栅格提取指定的研究区,然后提取栅格的值,再来绘图

虽然感觉不是那么花里胡哨的图,但这个应该还是比较实用的,特别是大批量提取栅格值的时候。由于在google colab里面操作的步骤比较多,中间可能有省略的地方,但重要的应该都在文中了,当然也可以迁移运用到其他地方,也可以查看一下这个第三方库的教程,比如read(1)是什么意思,官网的docs就写得有,实在是很方便的

以上就是用Python进行栅格数据的分区统计和批量提取的详细内容,更多关于Python 栅格数据的分区统计和批量提取 的资料请关注三水点靠木其它相关文章!

Python 相关文章推荐
python控制台英汉汉英电子词典
Apr 23 Python
python的类方法和静态方法
Dec 13 Python
简单介绍Python中的几种数据类型
Jan 02 Python
Python定时器实例代码
Nov 01 Python
详解Python自建logging模块
Jan 29 Python
Tensorflow使用支持向量机拟合线性回归
Sep 07 Python
Python unittest单元测试框架总结
Sep 08 Python
对python 命令的-u参数详解
Dec 03 Python
Python3enumrate和range对比及示例详解
Jul 13 Python
基于python的itchat库实现微信聊天机器人(推荐)
Oct 29 Python
python变量的作用域是什么
May 26 Python
Python破解极验滑动验证码详细步骤
May 21 Python
手把手教你怎么用Python实现zip文件密码的破解
Python基础学习之奇异的GUI对话框
教你使用pyinstaller打包Python教程
Matplotlib绘制混淆矩阵的实现
只需要100行Python代码就可以实现的贪吃蛇小游戏
PyTorch dropout设置训练和测试模式的实现
May 27 #Python
pytorch Dropout过拟合的操作
You might like
php页面消耗内存过大的处理办法
2013/03/18 PHP
Linux下安装oracle客户端并配置php5.3
2014/10/12 PHP
php使用PDO操作MySQL数据库实例
2014/12/30 PHP
PHP创建word文档的方法(平台无关)
2016/03/29 PHP
PHP入门教程之正则表达式基本用法实例详解(正则匹配,搜索,分割等)
2016/09/11 PHP
PHP数据对象PDO操作技巧小结
2016/09/27 PHP
PHP后台实现微信小程序登录
2018/08/03 PHP
PHP使用POP3读取邮箱接收邮件的示例代码
2020/07/08 PHP
prototype 1.5 & scriptaculous 1.6.1 学习笔记
2006/09/07 Javascript
索趣科技的答案
2007/02/07 Javascript
JQUERY设置IFRAME的SRC值的代码
2010/11/30 Javascript
jQuery学习笔记之jQuery选择器的使用
2010/12/22 Javascript
jQuery 幻灯片插件(带缩略图功能)
2011/01/24 Javascript
jQuery中closest()函数用法实例
2015/01/07 Javascript
简单介绍JavaScript的变量和数据类型
2015/06/03 Javascript
JavaScript中length属性的使用方法
2015/06/05 Javascript
超精准的javascript验证身份证号的具体实现方法
2015/11/18 Javascript
jQuery使用getJSON方法获取json数据完整示例
2016/09/13 Javascript
JS实现快速的导航下拉菜单动画效果附源码下载
2016/11/01 Javascript
AngularJS封装指令方法详解
2016/12/12 Javascript
实例解析angularjs的filter过滤器
2016/12/14 Javascript
详解nodejs微信公众号开发——5.素材管理接口
2017/04/11 NodeJs
Vue 进阶教程之v-model详解
2017/05/06 Javascript
bootstrap日期插件daterangepicker使用详解
2017/10/19 Javascript
JS 遍历 json 和 JQuery 遍历json操作完整示例
2019/11/11 jQuery
javascript实现页面的实时时钟显示示例
2020/08/06 Javascript
Python 命令行非阻塞输入的小例子
2013/09/27 Python
Python探索之创建二叉树
2017/10/25 Python
python3使用pandas获取股票数据的方法
2018/12/22 Python
Python argparse模块使用方法解析
2020/02/20 Python
亚马逊印度站:Amazon.in
2017/10/15 全球购物
越南母婴用品购物网站:Kids Plaza
2020/04/09 全球购物
关于保护环境的建议书
2014/08/26 职场文书
边城读书笔记
2015/06/29 职场文书
2015年电气技术员工作总结
2015/07/24 职场文书
python pandas 解析(读取、写入)CSV 文件的操作方法
2022/12/24 Python