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地图,用于可视化人口密度、绘制犯罪事件的空间模式、显示环境因素的分布,并根据地理数据分析疾病的传播。

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

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

相关文章

【Liunx】部署Ansible自动化运维工具

Ansible自动化运维工具 概述安装部署1.通过yum下载Ansible2.对自己做免密配置3.修改ansiable host配置对服务器进行分组4.测试:对所有服务器进行ping命令5.写playbook6.执行我们写的playbook脚本7.验证 概述 ansible是新出现的自动化运维工具,基于Pytho…

Python math 标准库学习总结

math --- 数学函数 该模块提供了对C标准定义的数学函数的访问。 这些函数不适用于复数;如果你需要计算复数,请使用 cmath 模块中的同名函数。将支持计算复数的函数区分开的目的,来自于大多数开发者并不愿意像数学家一样需要学习复数的概念。…

IDEA安装配置SceneBuilder

1、下载 SceneBuilder 地址: Scene Builder - Gluon ​​​​​​​ Scene Builder | JavaFX中文官方网站 选择符合自己系统版本,JAVA版本的下载 2、安装SceneBuilder 执行下载的 SceneBuilder 安装程序,并按照安装向导的指示进行…

Linux内核分析(十四)--内存管理之malloc、free 实现原理

目录 一、引言 二、malloc实现方式 ------>2.1、动态内存分配的系统调用:brk / sbrk ------>2.2、malloc实现思路 ------------>2.2.1、最佳适应法 ------------>2.2.2、最差适应法 ------------>2.2.3、首次适应法 ------------>2.2.4、下一个适应…

701. 二叉搜索树中的插入操作

原题链接:701. 二叉搜索树中的插入操作 思路: 因为是二叉搜索树,我们可以利用这个特性来确定val到底是在那一边的子树,而不是全部遍历二叉搜索树 直接使用递归,如果查询到root为NULL,则直接创建新结点即可…

[linux] 由创建用户开始

用户创建 【添加用户】 方案一&#xff08;自动化添加&#xff09;&#xff1a;adduser <用户名> 后按提示设计密码、全名等 方案二&#xff08;个性化添加&#xff09;&#xff1a;useradd -m <用户名>&#xff0c; 设置密码 passwd <用户名>。当出现切换用…

Feign调用返回值统一处理

因为接口一般有固定的返回格式&#xff0c;有数据、返回码和异常时错误信息。 不做任何处理情况下&#xff0c;需要在调用Feign的业务代码处通过判断返回码来获取数据。这种重复的代码可以抽出来统一处理。 解决方法1 使用自定义Decoder来统一处理&#xff0c;重写Object de…

高防IP可以抵御哪些恶意攻击

高防IP协议可以隐藏用户的站点&#xff0c;使得攻击者无法发现恶意攻击的目标网络资源&#xff0c;从而提高了源站的安全性。能够有效抵御常见的恶意攻击类型ICMPFlood、UDPFlood、 TCPFlood、SYNFlood、ACKFlood等&#xff0c;帮助游戏、金 融、电子商务、互联网、政企等行业抵…

实际使用Elasticdump工具对Elasticsearch集群进行数据备份和数据还原

文/朱季谦 目录一、Elasticdump工具介绍二、Elasticdump工具安装三、Elasticdump工具使用 最近在开发当中做了一些涉及到Elasticsearch映射结构及数据导出导入的工作&#xff0c;怕以后会把这过程忘记&#xff0c;可谓好记性不如烂笔头&#xff0c;故而记录成一篇博文。 玩El…

CCF编程能力等级认证GESP—C++4级—样题1

CCF编程能力等级认证GESP—C4级—样题1 单选题&#xff08;每题 2 分&#xff0c;共 30 分&#xff09;判断题&#xff08;每题 2 分&#xff0c;共 20 分&#xff09;编程题 (每题 25 分&#xff0c;共 50 分)第一题 绝对素数第二题 填幻方 参考答案单选题判断题编程题1编程题…

上海某公司医疗实施工程师笔试题

文章目录 一、Oracle数据库基础二、作为一个实施人员&#xff0c;用户如果提出我们软件功能所不具备的需求&#xff0c;你该如何处理&#xff1f;简单列出你可能采取的几种措施三、简答题三、答案1.Oracle数据库基础2.作为一个实施人员&#xff0c;用户如果提出我们软件功能所不…

基于mindspore的大模型llama2-7b---微调/推理

环境信息&#xff1a; 1.mindformers指导链接&#xff1a; https://gitee.com/mindspore/mindformers/blob/dev/docs/model_cards/llama2.md#%E6%A8%A1%E5%9E%8B%E6%9D%83%E9%87%8D%E4%B8%8B%E8%BD%BD%E4%B8%8E%E8%BD%AC%E6%8D%A2 2.镜像&#xff1a;mindformers0.8.0_ms2.2.0…

钉钉统计部门个人请假次数go

前言 最近小组需要统计部门各种请假次数&#xff0c;写了一个方法&#xff0c;第一次实战中用到递归函数&#xff0c;简单记录一下。 效果展示 这些数据不需要返回json&#xff0c;这里这样是为了方便测试。可以通过这些数据完成其它的操作。 功能实现 钉钉服务端调试工具A…

短剧公众号系统:连接创作者与观众的桥梁

随着数字媒体的迅速发展&#xff0c;短剧公众号系统作为一种新型的内容传播方式&#xff0c;逐渐在社交媒体平台上崭露头角。它通过将创作者与观众紧密连接&#xff0c;为双方提供了丰富的内容与互动体验。本文将探讨短剧公众号系统的价值&#xff0c;以及如何通过这种系统连接…

Java主流分布式解决方案多场景设计与实战

Java的主流分布式解决方案的设计和实战涉及到多个场景&#xff0c;包括但不限于以下几点&#xff1a; 分布式缓存&#xff1a;在Java的分布式系统中&#xff0c;缓存是非常重要的一部分。常用的分布式缓存技术包括Redis、EhCache等。这些缓存技术可以用来提高系统的性能和响应…

多个Obj模型合并

MergeObj&#xff08;合并Obj模型&#xff09; 1 概述 由于项目原因&#xff0c;需要下载谷歌地图上的模型&#xff0c;关于谷歌模型下载的&#xff0c;见我的CSDN博客. 由于下载谷歌地图上的数据&#xff0c;会分多个模块下载。下载完成后&#xff0c;怎么合并&#xff0c;在…

【数据结构初阶】链表OJ

链表OJ 题目一&#xff1a;移除链表元素题目二&#xff1a;反转链表题目三&#xff1a;链表的中间节点题目四&#xff1a;链表中倒数第k个结点题目五&#xff1a;合并两个有序链表题目六&#xff1a;链表分割题目七&#xff1a;链表的回文结构题目八&#xff1a;相交链表题目九…

B2B企业如何打造独立站:从策略到实施的全面指南

随着数字化转型的加速&#xff0c;B2B企业越来越认识到独立站的重要性。然而&#xff0c;如何建设一个优秀的独立站&#xff0c;以及如何将独立站与企业的整体战略相结合&#xff0c;是许多企业面临的挑战。本文将详细探讨B2B企业如何从策略到实施打造一个成功的独立站。 一、…

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

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

2656. K 个元素的最大和 --力扣 --JAVA

题目 给你一个下标从 0 开始的整数数组 nums 和一个整数 k 。你需要执行以下操作 恰好 k 次&#xff0c;最大化你的得分&#xff1a; 从 nums 中选择一个元素 m 。 将选中的元素 m 从数组中删除。 将新元素 m 1 添加到数组中。 你的得分增加 m 。 请你返回执行以上操作恰好 k …