5个用于地理空间数据分析的Python包

大家好,地理空间数据分析在城市规划、环境研究、农业和交通等行业中至关重要。不断增长的需求导致Python包在各种地理数据分析要求上的使用增加,例如分析气候模式、研究城市发展或跟踪疾病传播等,评估和选择具有快速处理、修改和可视化能力的正确工具对于有效分析和可视化地理空间数据至关重要。

首先,了解什么是地理空间数据至关重要。地理空间数据是具有地理成分的数据,表示地球表面上对象、特征或事件的位置和特性。它描述了物理宇宙中各种物体的空间连接、分布和属性。地理空间数据主要有两种类型:

  • 栅格数据(Raster data):适用于没有固定边界的连续信息,表示为具有值的单元格网格,这些值表示观察到的特征。通常在固定的时间间隔上进行监测,并进行插值以创建连续表面。

  • 矢量数据(Vector data):使用点、线和多边形来表示空间属性,包括兴趣点、交通网络、行政边界和土地地块等,通常用于具有精确位置或硬约束的离散数据。

地理空间数据可以以多种格式存储,如 ESRI Shapefile、GeoJSON、Erdas Imagine Image File Format(EIF)、GeoTIFF、Geopackage(GPKG)、GeoJSON、Light Detection、Ranging(LiDAR)等。

地理空间数据涵盖了各种类型,例如卫星图像、高程模型、点云、土地利用分类和基于文本的信息,为跨行业的空间分析和决策提供了有价值的洞察,微软、谷歌、Esri和亚马逊网络服务等大型公司都利用地理空间数据来获得有价值的洞察。接下来一起了解一下地理空间数据分析的五大Python包,这些包能够进行数据读取/写入、操作、可视化、地理编码和地理索引,适用于初学者和有经验的用户,为地理空间数据的有效探索、可视化和洞察提供支持。

1. Geopandas

适用于:矢量数据

Geopandas是一个广泛使用的Python库,用于处理矢量地理空间数据,在Pandas DataFrames中提供直观的地理数据处理,它支持Shapefile和GeoJSON等格式,并提供合并、分组和空间连接等空间操作。Geopandas可以与流行的库(如Pandas、NumPy和Matplotlib等)无缝集成,它可以处理大型数据集,但这可能带来挑战。Geopandas包通常用于空间数据分析任务,包括空间连接、查询以及缓冲区和交叉分析等地理空间操作,需要使用不同的包,如Shapely处理几何操作,Fiona访问文件,以及matplotlib绘图。

例如,可以使用Geopandas来探索房地产数据,以确定城市中最昂贵的街区,或者分析人口数据以可视化不同社区的增长和迁移模式。

可以使用以下pip命令安装该包:

pip install geopandas

使用GeoPandas进行绘图,通过以下代码查看内置地图:

import geopandas 
# 检查可用地图
geopandas.datasets.available

本文将使用GeoPandas加载世界地图数据集,提取美国的形状文件,并使用以下代码将其绘制在图表上:

# 选择特定地图
geopandas.datasets.get_path('naturalearth_lowres')
# 打开所选地图-GeoDataFrame
world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))
# 创建GeoDataFrame的子集
usa = world[world.name == "United States of America"]
# 绘制子集
usa.plot();

以上代码将打印出一个子集数据框的地图:

图片

2. Folium

适用于:点云

Folium是一个用于创建带有标记、弹出式窗口、choropleth和其他地理空间可视化的Python库,它与Leaflet JavaScript库集成,并允许将地图导出为HTML。并且可以与Geopandas和Cartopy结合使用,并使用Map Tiles处理大型数据集。Folium在简洁性、美观性以及与其他地理空间库的集成方面表现出色,然而对于高级地理空间分析和操作,Folium可能存在一些局限性。

例如,Folium可用于供应链和物流领域,用于可视化分销网络、优化路线和监测货物位置及装运地点。可以使用以下命令安装Folium:

pip install folium

现在可以使用以下几行代码在以[0, 0]为中心的示例交互式地图上放置一个标记:

import folium
# 以中心坐标(0, 0)生成Folium地图
map = folium.Map(location=[0, 0], zoom_start=2)
# 定位坐标(0, 0)
folium.Marker([0, 0]).add_to(map)
# 显示地图
map

图片

可以根据特定的地理空间数据进一步自定义此地图,例如添加标记、图层或样式选项。

3. ipyleaflet

适用于:点云,交互式

ipyleaflet包能够在Python中轻松创建交互式地图,特别适用于Jupyter笔记本,在其中允许用户生成和共享具有各种底图、标记和其他地理空间操作的交互式地图。ipyleaflet基于Leaflet JavaScript库构建,支持GeoJSON和WMS图层、CSS和JavaScript样式以及地理空间计算。虽然ipyleaflet在交互式小部件方面表现出色,但由于其依赖于JavaScript,对于纯Python项目可能不是理想选择。

例如,ipyleaflet可以应用于环境监测,以可视化传感器数据、监测空气质量并实时评估环境变化。要安装ipyleaflet,可以使用以下pip命令:

pip install ipyleaflet

使用ipyleaflet绘图,用以下代码在地图上创建一个标记,该标记位于纽约市的坐标(40.7128,-74.0060)处,以表示感兴趣的点:

from ipyleaflet import Map, Marker
# 创建地图
m = Map(center=(40.7128, -74.0060), zoom=12)
# 添加标记
marker = Marker(location=(40.7128, -74.0060))
m.add_layer(marker)

以下是代码的输出结果:

图片

4. Rasterio

适用于:栅格数据

Rasterio是一个功能强大的Python库,用于处理地理空间栅格数据,具有高效的性能和广泛的操作,如裁剪、重投影和重采样。它支持各种栅格格式,并与其他地理空间库良好集成,尽管在处理矢量数据和复杂分析任务方面存在一些限制。然而,Rasterio是Python中高效操作和处理栅格数据的必备工具。

例如,Rasterio可用于读取和写入卫星图像、进行地形分析、从数字高程模型中提取数据以及进行遥感分析等任务。

!pip install rasterio

rasterio.open()函数用于打开文件,read()方法将图像读取为NumPy数组。最后,使用Matplotlib中的plt.imshow()函数显示图像,plt.show()在输出中显示绘图。

import rasterio
from rasterio.plot import show

使用rasterio库打开并可视化来自kaggle数据集“High-resolution GeoTIFF images of climatic data”中的sample.tif文件的栅格图像,将红色通道(图像中的一种颜色通道)作为子图使用Reds颜色映射进行显示,将原始图像(包含多个颜色通道)作为另一个子图使用viridis颜色映射进行显示。也可以使用相同的方法可视化其他颜色通道,如绿色和蓝色。

src = rasterio.open('/content/sample.tif')
plt.figure(figsize=(15,10))
fig, (axr, axg) = plt.subplots(1,2, figsize=(15,7))
show((src, 1), ax=axr, cmap='Reds', title='red channel')
show((src), ax=axg, cmap='viridis', title='original image')
plt.show()

图片

原始GeoTIFF图像(右)来源:Kaggle.com

https://www.kaggle.com/datasets/abireltaief/highresolution-geotiff-images-of-climatic-data

在地理空间分析中,分析特定的颜色通道(如红色、蓝色和绿色)是为了关注或提取与图像中那些颜色分量相关的特定属性、特征或特征所表示的有价值信息,可以用于遥感中的植被健康、植被指数或水体等方面。

5. Geoplot

适用于:矢量数据,交互式

Geoplot是一个用户友好型的Python库,可以快速创建具有吸引力的地理空间可视化效果,包括等值线地图和散点图,它与流行的数据处理库(如Pandas等)无缝集成,并支持多种地图投影。然而,Geoplot在交互式地图支持和绘图类型方面存在一些限制,比专门的地理空间库要少。尽管如此,它仍然对于快速实现地理空间数据可视化和获取空间模式的洞察力非常有价值。

!pip install geoplot

本文将使用Geoplot绘制一个等值图可视化,根据世界形状文件中的“continent”属性选择亚洲国家,根据“pop_est”属性分配颜色强度,并使用带图例的“icefire”彩色地图进行绘制,图例的尺寸大小为10x5。

import geoplot
# 绘制亚洲人口分布图
asia = world.query("continent == 'Asia'")
geoplot.choropleth(asia, hue = "pop_est", cmap = "icefire",legend=True, figsize = (10, 5));

图片

例如,geoplot包可以创建choropleth地图,用于可视化人口密度、绘制犯罪事件的空间模式、显示环境因素的分布,并根据地理数据分析疾病的传播。

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

【移远QuecPython】EC800M物联网开发板的MQTT协议腾讯云数据上报

【移远QuecPython】EC800M物联网开发板的MQTT协议腾讯云数据上报 文章目录 导入库初始化设置MQTT注册回调订阅发布功能开启服务发送消息函数打包调用测试效果附录:列表的赋值类型和py打包列表赋值BUG复现代码改进优化总结 py打包 导入库 from TenCentYun import TX…

【Java】线程的调度、生命周期及状态转换

🌺个人主页:Dawn黎明开始 🎀系列专栏:Java ⭐每日一句:夜色难免黑凉,前行必有曙光 📢欢迎大家:关注🔍点赞👍评论📝收藏⭐️ ​ 文章目录 一.&…

day08_子网划分与子网掩码

什么是子网划分? 1、概念:借主机位给网络位使用,以此来达到把一个大网段划分为n个儿子网段的目的,2. 为何要进行子网划分?3、子网掩码:就是对ip地址打记号4、 网络地址的计算机方式:ip地址与子网掩码都转换成二进制&a…

链表面试题

💓作者简介👏:在校大二迷茫大学生 💖个人主页🎉:小李很执着 💗系列专栏:Leetcode经典题 每日分享:其实要过那条马路并不难,就看谁在对面等你❣️❣️❣️ 目录…

MySQL数据库约束

目录 数据库约束 1.NULL约束 2.UNIQUE:唯一约束 3.DEFAULT:默认值约束 4.PRIMARY KEY:主键约束 5.FOREIGN KEY:外键约束 数据库约束 以下为本篇文章会介绍的约束 (1)NOT NULL - 指示某列不能存储 NULL 值。 (2)UNIQUE - …

Spark SQL编程

1. Spark SQL概述 1.1 什么是Spark SQL Spark SQL是用于结构化数据处理的Spark模块。与基本的Spark RDD API不同,Spark SQL提供的接口为Spark提供了有关数据结构和正在执行的计算的更多信息。在内部,Spark SQL使用这些额外的信息来执行额外的优化。与Spa…

qt+opengl 着色器VAO、VBO、EBO(四)

文章目录 一、顶点着色器和片段着色器代码分析1. 着色器12. 顶点着色器2 二、使用步骤1. 使用着色器12. 使用着色器23. 在着色器2中使用EBO 三、完整代码 一、顶点着色器和片段着色器代码分析 1. 着色器1 用到的坐标矩阵, 四个四边形顶点坐标 float vertices_data[36] {// 所…

PlantUML基础使用教程

环境搭建 IDEA插件下载 打开IEDA系列IDE,从FIle–>Settings–>Plugins–>Marketplace 进入到插件下载界面,搜索PlantUML,安装PlantUML Integration和PlantUML Parser两个插件,并重启IDE 安装和配置Graphviz 进入官网…

C/C++轻量级并发TCP服务器框架Zinx-框架开发001: 读取标准输入,回显到标准输出

文章目录 完整代码实现参考-非项目使用项目使用的代码 - 乱-但是思路与上面的相同创建Kernel类添加删除修改epoll&#xff0c;才能写run方法创建stdin_Channel类在Kernel类中实现run方法 完整代码实现参考-非项目使用 #include <errno.h> #include <signal.h> #in…

蓝桥杯每日一题2023.11.14

题目描述 题目分析 此题目的最终目标是将字母都填上数使等式符合条件&#xff0c;实际我们发现可以使用搜索将所有符合条件的进行判断&#xff08;答案&#xff1a;29&#xff09; 由于小数可能会出现错误故我们将其进行简单变化进行搜索 #include<bits/stdc.h> using…

No207.精选前端面试题,享受每天的挑战和学习

🤍 前端开发工程师(主业)、技术博主(副业)、已过CET6 🍨 阿珊和她的猫_CSDN个人主页 🕠 牛客高级专题作者、在牛客打造高质量专栏《前端面试必备》 🍚 蓝桥云课签约作者、已在蓝桥云课上架的前后端实战课程《Vue.js 和 Egg.js 开发企业级健康管理项目》、《带你从入…

Pytorch自动混合精度的计算:torch.cuda.amp.autocast

1 autocast介绍 1.1 什么是AMP? 默认情况下&#xff0c;大多数深度学习框架都采用32位浮点算法进行训练。2017年&#xff0c;NVIDIA研究了一种用于混合精度训练的方法&#xff0c;该方法在训练网络时将单精度&#xff08;FP32&#xff09;与半精度(FP16)结合在一起&#xff…

Dart利用私有构造函数_()创建单例模式

文章目录 类的构造函数_()函数dart中构造函数定义 类的构造函数 类的构造函数有两种&#xff1a; 1&#xff09;默认构造函数&#xff1a; 当实例化对象的时候&#xff0c;会自动调用的函数&#xff0c;构造函数的名称和类的名称相同&#xff0c;在一个类中默认构造函数只能由…

post 和get参数 请求

json参数 post请求格式 RestController public class HelloController { //json参数 post 请求RequestMapping("/jsonParam")public String jsonParam(RequestBody User user){System.out.println(user);return "OK";} } postman 接口测试工具…

spring cloud alibaba 简介

微服务搭建组件选型 1.服务注册中心 Nacos(spring-cloud-alibaba) 2.服务通信 OpenFeign(spring-cloud) 3.服务熔断、降级、限流 Sentinel(spring-cloud-alibaba) 4.网关 Gateway(spring-cloud) 5.服务配置中心 …

MySQL被攻击后创建数据库报错1044 - Access denied for user ‘root‘@‘%‘ to database ‘xxx‘

MySQL被攻击后创建数据库报错1044 - Access denied for user root% to database xxx 一、问题二、解决过程1、正常过程2、踩坑&#xff08;已经解决问题的可以不看&#xff09; 一、问题 最近数据库被攻击了&#xff0c;业务数据库都没了 还好也不是有重要数据&#xff0c;但再…

【HUST】网安纳米|2023年研究生纳米技术考试参考

目录 1 纳米材料是什么 2 纳米材料的结构特性 3 纳米结构的其他特性 4 纳米结构的检测技术 5 纳米材料的应用 打印建议&#xff1a;PPT彩印&#xff08;这样重点比较突出&#xff09;&#xff0c;每面12张PPT&#xff0c;简单做一下关键词目录&#xff0c;亲测可以看清。如…

Uniapp开发 购物商城源码 在线电商商城源码 适配移动终端项目及各小程序

lilishop电商商城系统 商城移动端&#xff0c;使用Uniapp开发&#xff0c;可编译为所有移动终端项目及各小程序 源码下载&#xff1a;https://download.csdn.net/download/m0_66047725/88487579 源码下载2&#xff1a;关注我留言

02 # 类型基础:强类型与弱类型

宽泛的定义 在强类型语言中&#xff0c;当一个对象从调用函数传递到被调用函数时&#xff0c;其类型必须与被调用函数中声明的类型兼容 – Liskov, Zilles 1974 通俗定义 强类型语言不允许改变变量的数据类型&#xff0c;除非进行强制类型转换 比如下面 Java 里不能将布尔类…

最小二乘法及参数辨识

文章目录 一、最小二乘法1.1 定义1.2 SISO系统运用最小二乘估计进行辨识1.3 几何解释1.4 最小二乘法性质 二、加权最小二乘法三、递推最小二乘法四、增广最小二乘法 一、最小二乘法 1.1 定义 1974年高斯提出的最小二乘法的基本原理是未知量的最可能值是使各项实际观测值和计算…