这是 Python 库 GeoPandas 的用武之地。 本文和大家一起学习如何使用 GeoPandas有效地可视化地理空间数据。...CRS/坐标参考系统告诉我们如何（使用投影 或数学方程）将圆形地球上的位置（坐标）转换为扁平的二维坐标系（例如计算机屏幕或纸张）上的相同位置地图）。最常用的 CRS 是“EPSG:4326”。...▲ GeoPandas – GeoDataFrame 和 GeoSeries 在GeoPandas的主要数据结构是GeoDataFrame延伸的PandasDataFrame。...GeoDataFrame包含一个或多个GeoSeries（延伸PandasSeries）每个都包含在一个不同的几何形状的投影（GeoSeries.crs）。...数据准备 在导入 GeoPandas 之前阅读Teams数据集，数据集和代码可以在公众号『数据STUDIO』回复【GeoPandas】获取。
Berlin_transport.head() 图4 对站点的空间分布进行简单的可视化： 图5 接着我们就利用sjoin()将区划面作为左表，站点作为右表，在op='intersects'参数设置下进行空间连接，再衔接groupby..., right_df=Berlin_transport.query("fclass=='bus_stop'"), op='intersects') \ .groupby...keep_geom_type：同叠加分析overlay中的同名参数 基于实际例子进行演示，我们读入数据berlin_footway_WGS84.shp，包含了柏林全部的步道路网线数据，并转换到适合柏林地区的投影...空间数据分析系列第一篇文章，到今天这篇为止，geopandas中全部实用的主线内容（截至0.7.0版本），都在这断断续续撰写完成的9篇文章中介绍完毕，不敢说是geopandas中文资料里最好的，但穿插了众多例子和举一反三的内容...相关内容学习精进的过程记录下来，通过博客与微信公众号与广大的读者朋友共同交流学习，期间认识了很多业内大牛和朋友，收获了很多很多。
图1 而创建缓冲区时也需要遵循一定的参数，从而决定怎样向几何对象外进行缓冲，geopandas中buffer()和shapely中的buffer()方法参数一致，主要参数如下： distance：...用于指定向外缓冲的距离，单位与矢量数据自带单位保持一致，在常见的投影坐标系如Web Mercator（EPSG:3857）下就是以米为单位，因此需要注意一定要先将矢量数据转换为合适的投影坐标系之后，再进行缓冲区分析才是合理有效的...图11 translate() translate()用于实现矢量对象的平移操作，其主要参数有xoff和yoff，分别控制在x维度和y维度上的平移距离（与对应的投影单位保持一致）: ?...2.4 空间融合与拆分 有时候我们希望对矢量数据按照某些字段进行分组，再分别对非矢量列与矢量列进行聚合及合并，类似于pandas中的groupby.agg()；而有些时候我们希望把矢量类型为Multi-xxx...通过geopandass中的dissolve()和explode()方法，我们就可以实现上述功能： dissolve() dissolve()用于对矢量数据进行融合，可以理解为对矢量数据进行groupby
关于其介绍和安装等请参考其github主页，本文不再赘述。 三、子区域数据分类统计 直接进入正题，现有某省的分类统计数据shp文件以及此省的行政区划数据shp文件。...geopandas进行投影转换很简单。...代码如下： gpd_new=gpd.to_crs(crs) 其中gpd表示原始数据，gpd_new为转换后的数据，crs表示需要转换的投影参数，在https://github.com/geopandas...所以我们只需要将上述两种数据转换到同一投影即可，问题是假设我们不知道它们的投影类型，那么也很容易，只需要将其都转换成4326或其他投影即可，这样就能保证二者转换后为同一投影。...当然并一定局限于省和市，比如全球和国家或者国家和省等。只要存在包含关系即可通过此种方式进行处理。这是鸡年的第一篇博客，愿所有人今年都能有个好的结果！
图5 接着我们就利用sjoin()将区划面作为左表，站点作为右表，在op='intersects'参数设置下进行空间连接，再衔接groupby，以统计出各区划面内部的公交站点数量： gpd.sjoin..., right_df=Berlin_transport.query("fclass=='bus_stop'"), op='intersects') \ .groupby...keep_geom_type：同叠加分析overlay中的同名参数 基于实际例子进行演示，我们读入数据berlin_footway_WGS84.shp，包含了柏林全部的步道路网线数据，并转换到适合柏林地区的投影...篇文章中介绍完毕，不敢说是geopandas中文资料里最好的，但穿插了众多例子和举一反三的内容，绝对是帮助大家理解学习geopandas非常实在的参考资料。 ...相关内容学习精进的过程记录下来，通过博客与微信公众号与广大的读者朋友共同交流学习，期间认识了很多业内大牛和朋友，收获了很多很多。
作为基于geopandas的空间数据分析系列文章的第二篇，通过本文你将会学习到geopandas中的坐标参考系管理。...上述的这些用于在不同情况下定义对象位置信息的坐标系统，就称为坐标参考系统（Coordinate Reference System，下文统称CRS）： 图2 CRS可细分为地理坐标系和投影坐标系。...地理坐标系虽然解决了我们在地球球面上定位的问题，但纬度和经度位置没有使用统一的测量单位。...在上一篇文章（数据科学学习手札74）基于geopandas的空间数据分析——数据结构篇中我们介绍了创建GeoSeries和GeoDataFrame的方法。...实际上，现实的空间分析计算任务中，必须要为数据设置合适的CRS，在geopandas.GeoSeries()和geopandas.GeoDataFrame()中就包含参数crs。
作为基于geopandas的空间数据分析系列文章的第二篇，通过本文你将会学习到geopandas中的坐标参考系管理。...就又需要有一套将地球球面展平的方法，上述的这些用于在不同情况下定义对象位置信息的坐标系统，就称为坐标参考系统（Coordinate Reference System，下文统称CRS）： 图2 CRS可细分为地理坐标系和投影坐标系...3 geopandas中的坐标参考系管理 至此，我们已经对CRS有了较为全面的了解，打好了基础，接下来我们来正式学习geopandas中的坐标参考系管理，geopandas调用pyproj作为CRS... 在上一篇文章（数据科学学习手札74）基于geopandas的空间数据分析——数据结构篇中我们介绍了创建GeoSeries和GeoDataFrame的方法，实际上，现实的空间分析计算任务中，必须要为数据设置合适的...CRS，在geopandas.GeoSeries()和geopandas.GeoDataFrame()中就包含参数crs，下面我们举例说明，还是先用到geopandas自带的世界国家地区数据，我们从中选择中国
()和shapely中的buffer()方法参数一致，主要参数如下： distance：用于指定向外缓冲的距离，单位与矢量数据自带单位保持一致，在常见的投影坐标系如Web Mercator（EPSG:3857...这时对矢量数据进行简化就非常有必要，geopandas中沿用shapely中的simplify()方法，帮助我们对过于复杂的线和面进行简化，和QGIS中简化矢量的方法一样，simplify()使用了科学的...，分别控制在x维度和y维度上的平移距离（与对应的投影单位保持一致）: 图12 2.3 叠加分析 geopandas基于shapely中的overlay()，为GeoDataFrame赋予了同样的可以作用到整个矢量列的...2.4 空间融合与拆分 有时候我们希望对矢量数据按照某些字段进行分组，再分别对非矢量列与矢量列进行聚合及合并，类似于pandas中的groupby.agg()； 而有些时候我们希望把矢量类型为Multi-xxx...通过geopandass中的dissolve()和explode()方法，我们就可以实现上述功能： dissolve() dissolve()用于对矢量数据进行融合，可以理解为对矢量数据进行groupby
计划涵盖geopandas的数据结构、投影坐标系管理、文件IO、基础地图制作、集合操作、空间连接与聚合。...geopandas的安装和使用需要若干依赖包，如果不事先妥善安装好这些依赖包而直接使用pip install geopandas或conda install geopandas，可能会引发依赖包相关错误导致安装失败...2 数据结构 geopandas作为pandas向地理分析计算方面的延拓，基础的数据结构延续了Series和DataFrame的特点，创造出GeoSeries与GeoDataFrame两种基础数据结构：...pandas中的Series，GeoSeries在被创建完成之后也拥有很多实用的地理属性，下面对其中较为常用的进行列举： area area属性返回与GeoSeries中每个元素一一对应的面积值（这里的面积单位和下文涉及的长度单位取决于投影坐标系...中的投影坐标系管理，敬请期待。
官方用户手册中文版用户手册投影操作等查询 Pysal-空间计量库 Pysal是一个面向地理空间数据科学的开源跨平台库，重点是用python编写的地理空间矢量数据。...它支持空间分析高级应用程序的开发，例如 空间簇、热点和异常点的检测从空间数据构建图形地理嵌入网络的空间回归与统计建模空间计量经济学探索性时空数据分析 官方手册中文文档空间计量综述 Geopandas...-空间数据分析库 Geopandas可以在Python下更方便的处理地理空间数据。...GeoPandas进一步依赖于 fiona进行文件存取和 descartes ，matplotlib 进行绘图。 ...本机投影支持：地理空间绘图的最基本特性是投影：如何以正确的方式将球体展开到平坦的表面（地图）上取决于要描绘的内容。geoplot提供这些选项。
图1 2.1.1 shapefile 作为非常常见的一种矢量文件格式，geopandas对shapefile提供了很好的读取和写出支持，下面分为不同情况来介绍： 完整的shapefile 如图2...图3 缺少投影的shapefile 当shapefile中缺失.prj文件时，使用geopandas读入后形成的GeoDataFrame会缺失crs属性： ?...图14 2.1.4 过滤 geopandas在0.1.0版本中新增了bbox过滤，在0.7.0版本中新增了蒙版过滤和行过滤功能，可以辅助我们根据自己的需要读入原始数据中的子集，下面一一进行介绍： bbox...，只需要将driver参数设置为GPKG即可，这里需要注意一个bug：在使用geopandas导出GeoPackage文件时，可能会出现图21所示错误： ?...图21 但我观察到即使出现了上述错误，GeoPackage文件也是成功保存到路径下的且整个程序并未被打断，因此可以无视上述错误： ?
本系列文章就将围绕geopandas及其使用过程中涉及到的其他包进行系统性的介绍说明，每一篇将尽可能全面具体地介绍geopandas对应方面的知识，计划涵盖geopandas的数据结构、投影坐标系管理、...geopandas的安装和使用需要若干依赖包，如果不事先妥善安装好这些依赖包而直接使用pip install geopandas或conda install geopandas可能会引发依赖包相关错误导致安装失败...2 数据结构 geopandas作为pandas向地理分析计算方面的延拓，基础的数据结构延续了Series和DataFrame的特点，创造出GeoSeries与GeoDataFrame两种基础数据结构...中的Series，GeoSeries在被创建完成之后也拥有很多实用的地理属性，下面对其中较为常用的进行列举： area area属性返回与GeoSeries中每个元素一一对应的面积值（这里的面积单位和下文涉及的长度单位取决于投影坐标系...图37 以上就是本文的全部内容，如有笔误望指出，系列文章下一篇将详细介绍geopandas中的投影坐标系管理，敬请期待。
： geopandas 绘制中国地图 matplotlib add_axes()添加南海小地图 绘图文件分享 geopandas 读取中国地图文件 geopandas提供了非常方便的read_file...我们进行投影转换(epsg=2343)和进行一些简单的设置，代码如下： fig, ax = plt.subplots(figsize=(12, 8),dpi=80) ax = china_main.geometry.to_crs...2343.geometry[3].y) # 移除子图坐标轴刻度， ax_child.set_xticks([]) ax_child.set_yticks([]) 可以发现，除了显示范围的不同外，其他的和绘制主题部分的代码一致...总结 本期推文使用了Python-geopandas进行了中国地图的绘制,讲解了数据标记，投影转换等内容。...但需指出的是： geopandas 的安装较为麻烦，建议使用 conda install --channel conda-forge geopandas 进行安装。
大多数 GIS 软件和数据库都提供了一种机制来计算和使用数据图层的空间索引。QGIS 和 PostGIS 使用基于 R-Tree 数据结构的空间索引方案 - 它使用几何边界框创建分层树。...在这篇文章中，我将向你展示如何创建使用点密度图geopandas和h3-py库在Python。 国家地理空间情报局的海事安全信息门户以反航运活动消息的形式提供所有海盗事件的形状文件。...import geopandas as gpd from h3 import h3 GeoPandas 允许直接从 zip 文件中读取数据层。...我们groupby在h3列上使用 Panda 的函数，并count在输出中添加一个新列，其中包含每个 H3 id 的行数。...counts = incidents.groupby(['h3']).h3.agg('count').to_frame('count').reset_index() 我们现在知道每个网格单元的盗版事件数量
本期我们试着使用Python-geopandas包绘制空间地图，主要的知识点如下： geopandas 绘制中国地图 matplotlib add_axes()添加南海小地图 绘图文件分享 geopandas...读取中国地图文件 geopandas提供了非常方便的read_file()方法用于读取geojson文件，我们直接进行默认投影(WGS84)的绘制，代码如下： file = r"中国省级地图GS（2019...(epsg=2343)和进行一些简单的设置，代码如下： fig, ax = plt.subplots(figsize=(12, 8),dpi=80) ax = china_main.geometry.to_crs...总结 本期推文使用了Python-geopandas进行了中国地图的绘制,讲解了数据标记，投影转换等内容。...但需指出的是： geopandas 的安装较为麻烦，建议使用 conda install --channel conda-forge geopandas 进行安装。
最近在梳理Python中可以制作数据地图的可视化工具包，分别实践了geopandas、folium、Basemp,通过对比发现，静态地图中最为成熟的最终还得是Basemap工具，它是...Basemap工具在地理信息读写、坐标映射、空间坐标转化与投影等方面做的要比geopandas更加成熟，它可以使用常规的地图素材数据源（shp）作为底图进行叠加绘图，效果与精度控制比较方便，图表质量堪比...:,["NickName","Province","Signature"]] 2、聚合计算好友地区分布： aggResult = mydata.groupby
，而在日常空间数据分析工作中矢量文件的读入和写出，是至关重要的环节。...本文使用到的所有数据都可以在文章开头提及的Github仓库对应本文路径下找到： 图1 2.1.1 shapefile 作为非常常见的一种矢量文件格式，geopandas对shapefile提供了很好的读取和写出支持...geometry/china_provinces/china_provinces.shp') print(data.crs) # 查看数据对应的crs data.head() # 查看前5行 图3 缺少投影的...下面我们来读入图13所示的文件： 图13 图14 2.1.4 过滤 geopandas在0.1.0版本中新增了bbox过滤，在0.7.0版本中新增了蒙版过滤和行过滤功能，可以辅助我们根据自己的需要读入原始数据中的子集...这里需要注意一个bug：在使用geopandas导出GeoPackage文件时，可能会出现图21所示错误： 图21 但我观察到即使出现了上述错误，GeoPackage文件也是成功保存到路径下的且整个程序并未被打断
今天要讲解的主角是R语言中的sf包和Python中的geopandas库。...---- 为什么今天把geopandas和R语言空间数据可视化写在一起，因为他们很巧合的用到了相同的地理信息处理技术，无论是数据源的支持上、还是空间数据的结构存储还是投影设置上都是如此。...geopandas包同时支持导入shp素材和json素材，导入之后得数据结构与R语言中得sf导入之后得结构是一致得，地理多边形边界点信息都被压缩成了一个非常整齐的列表存储，列表内每一个单独的子项目都代表着一个多边形...china_map.plot(column="AREA",figsize=(20,12),cmap="Greens") 如果你想要在此图层上添加另外一个散点图层，则需设定两个具有同样投影信息的GeoDataFrame...###设置相同的投影： china_map_ploygon.crs={'init': 'epsg:3395'} china_map_point.crs ={'init': 'epsg:3395'} #
，我们已经对geopandas的基础知识、基础可视化，以及如何科学绘制分层设色地图展开了深入的学习，而利用geopandas+matplotlib进行地理可视化固然能实现常见的地图可视化，且提供了操纵图像的极高自由度...nyc-collision-factors.geojson') nyc_collision_factors.head() 图3 首先我们使用geoplot中的polyplot来绘制纽约行政区划，这里使用geoplot自带的Albers等面积投影作为投影...当然你也可以自由地通过legend_values和legeng_labels这两个参数自定义图例内容。...值得注意的是，因为常见在线地图如谷歌地图、OpenStreetMap、高德地图等的投影均为EPSG:3857也就是我们常说的Web Mercator，所以一旦要使用webmap，则投影锁死为EPSG:3857...如有疑问和意见，欢迎留言或在我的Github仓库中发起issues与我交流。 -END-
领取专属 10元无门槛券
手把手带您无忧上云
云服务器
即时通信 IM
ICP备案
云直播
实时音视频
