Python之深入解析一行代码计算每个省面积的神器Geopandas

Python之深入解析一行代码计算每个省面积的神器Geopandas一 前言 GeoPandas 是一个基于 pandas 针对地理数据做了特别支持的第三方模块 它继承 pandas Series 和 pandas Dataframe 实现了 GeoSeries 和 GeoDataFrame 类 使得其操纵和分析平面几何对象非常方便 二 准备 Python 安装在 windows 上安装 Python 下载 Python 的最新版本 访问链接 Python 官网 在 Windows 操作系统上安装 Python3 10 0 大家也可以自行选择最新版

一、前言

  • GeoPandas 是一个基于 pandas,针对地理数据做了特别支持的第三方模块。
  • 它继承 pandas.Series 和 pandas.Dataframe,实现了 GeoSeries 和 GeoDataFrame 类,使得其操纵和分析平面几何对象非常方便。

在这里插入图片描述

二、准备

① Python 安装

  • 在 windows 上安装 Python:
    • 下载 Python 的最新版本:访问链接 Python 官网,在 Windows 操作系统上安装 Python 3.10.0,大家也可以自行选择最新版的下载。

在这里插入图片描述

    • 进入下一个页面后,拉到最下面,选择 Windows embeddable package (64-bit) 下载:

在这里插入图片描述

    • 双击下载完成的应用程序,记得一定要勾选 ADD Python 3.7 to PATH,否则 cmd 无法运行 Python,然后选择 Install Now 进行安装,默认会安装 pip,而且将安装到 C 盘(推荐)。如果希望换一个盘安装,可以点击 Customize installation 在里面选择其他的盘。
  • 在 macOS 上安装 Python:
    • 其实 macOS 上是默认安装了 Python 的,只不过版本是 Python2.7,其相对于现在的 Python 3 缺少了很多特性,而且速度也比 Python 3 慢,因此还是建议大家装 Python 3。
    • 同样地,下载 Python 的最新版本,访问链接:Python 官网,在 macOS 操作系统上安装 Python 3.7.4,大家也可以自行选择最新版的下载。

在这里插入图片描述

    • 进入下一个页面后,拉到最下面,选择 macOS 64-bit universal2 installer 下载:

在这里插入图片描述

    • 下载完成后得到一个 pkg 文件,双击打开,按照指示步骤进行安装即可。

② 安装依赖

  • 请选择以下任一种方式输入命令安装依赖:
    • Windows 环境打开 Cmd (开始 – 运行 – CMD);
    • MacOS 环境打开 Terminal (command + 空格输入 Terminal);
    • 如果使用的是 VSCode 编辑器或 Pycharm,可以直接使用界面下方的 Terminal。
pip install rich 
  • 由于 geopandas 涉及到许多第三方依赖,pip 安装起来非常麻烦,因此这里推荐使用 conda 安装 geopandas,一行语句即可完成安装:
conda install geopandas 

三、基本使用

  • 设定坐标绘制简单的图形:
import geopandas from shapely.geometry import Polygon p1 = Polygon([(0, 0), (1, 0), (1, 1)]) p2 = Polygon([(0, 0), (1, 0), (1, 1), (0, 1)]) p3 = Polygon([(2, 0), (3, 0), (3, 1), (2, 1)]) g = geopandas.GeoSeries([p1, p2, p3]) # g: # result: # 0 POLYGON ((0 0, 1 0, 1 1, 0 0)) # 1 POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)) # 2 POLYGON ((2 0, 3 0, 3 1, 2 1, 2 0)) # dtype: geometry 
  • 这些变量所形成的图形如下:

在这里插入图片描述

  • 这里有一个重要且强大的用法,通过 area 属性,geopandas 能直接返回这些图形的面积:
>>> print(g.area) 0 0.5 1 1.0 2 1.0 dtype: float64 
  • 不仅如此,通过 plot 属性函数,还可以直接生成 matplotlib 图:
>>> g.plot() 
  • 通过 matplot 的 pyplot,可以将图片保存下来:
import matplotlib.pyplot as plt g.plot() plt.savefig("test.png") 
  • 学会上面的基本用法, 就可以进行简单的地图绘制及面积的计算。

四、绘制并算出每个省的面积

  • 此外,它最大的亮点是可以通过 Fiona(底层实现,用户不需要管),读取比如 ESRI shapefile(一种用于存储地理要素的几何位置和属性信息的非拓扑简单格式)。
import geopandas import matplotlib.pyplot as plt from shapely.geometry import Polygon maps = geopandas.read_file('1.shx') # 读取的数据格式类似于 # geometry # 0 POLYGON ((.341 .278, .591... # 1 POLYGON ((-.375 .367, -.3... # 2 POLYGON ((.092 .843, .873 4... # 3 POLYGON ((-34477.046 .963, -41105.128 4... # ... ... maps.plot() plt.savefig("test.png") 
  • 如代码所示,通过 read_file,可以读取 shx、gpkg、geojson 等数据。读取出来的图形如下:

在这里插入图片描述

  • 同样,这个 shapefile 是省级行政区的,每一个省级行政区都被划分为一个区块,因此可以一行语句算出每个省级行政区所占面积:
print(maps.area) # 0 4.e+11 # 1 1.e+12 # 2 1.e+11 # 3 4.e+10 # 4 1.e+12 # 5 1.e+11 # 6 1.e+11 # 7 9.e+10 # 8 1.e+11 # 9 1.e+11 # 10 1.015979e+11 # ... ... 
  • GeoPandas 的更多特性,请参考:GeoPandas。
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/233048.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 支持js代码的博客有哪些?[通俗易懂]

    支持js代码的博客有哪些?[通俗易懂]1、支持js代码的博客谷歌博客,雅虎奇摩,博客大巴,jimdo博客,和讯博客要等级(230000分).2、现在新浪博客怎么才能有效的使用js呢在新浪博客放广告,现在网络上主要流行的两种方法,这里尽量简单化的给您说明。把广告以代码形式放到免费空间里,然后链接到你的博客教程如下:(1).登陆googleadsense后台,生成广告代码。(2).将代码贴到纪事本文件,另存为….

    2022年7月12日
    15
  • Java二叉树前序遍历[通俗易懂]

    Java二叉树前序遍历[通俗易懂]给你二叉树的根节点root,返回它节点值的前序遍历。示例1:输入:root=[1,null,2,3]输出:[1,2,3]示例2:输入:root=[]输出:[]示例3:输入:root=[1]输出:[1]示例4:输入:root=[1,2]输出:[1,2]示例5:输入:root=[1,null,2]输出:[1,2]提示:树中节点数目在范围[0,100]内-100<=Node.val<=100进阶:递归算法很简单

    2025年9月23日
    3
  • navicat15激活碼3月最新在线激活

    navicat15激活碼3月最新在线激活,https://javaforall.net/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

    2022年3月15日
    56
  • 打开虚拟机时出现不能为虚拟电脑打开一个新任务「建议收藏」

    打开虚拟机时出现不能为虚拟电脑打开一个新任务「建议收藏」标题:打开虚拟机时出现不能为虚拟电脑打开一个新任务在用虚拟机打开Ubuntu时出现以下情况解决方法在查找了许多有关资料试用无效后,最终用以下两个步骤解决了该问题1.打开VirtualBox安装文件夹里的\drivers\vboxdrv文件夹2.右键VBoxDrv.inf文件,点击安装;3.安装完成后重启VirtualBox。参考文章在参照该作者成功打开一次后续仍然出现原问题后续发现应该是权限问题每次打开必须用管理员身份,直接双击是不可以的。这样问题就解决啦。…

    2025年6月18日
    5
  • kafka-logs_logstash kafka

    kafka-logs_logstash kafkaKafka的Log存储解析标签(空格分隔):kafka引言Kafka中的Message是以topic为基本单位组织的,不同的topic之间是相互独立的。每个topic又可以分成几个不同的partition(每个topic有几个partition是在创建topic时指定的),每个partition存储一部分Message。借用官方的一张图,可以直观地看到topic和partit

    2022年10月14日
    2
  • async/await 的理解和用法

    async/await 的理解和用法async/await是什么async/await是ES7提出的基于Promise的解决异步的最终方案。asyncasync是一个加在函数前的修饰符,被async定义的函数会默认返回一个Promise对象resolve的值。因此对async函数可以直接then,返回值就是then方法传入的函数。//async基础语法asyncfunctionfun0(){console.log(1);return1;}fun0().then(val=>{conso

    2022年7月12日
    14

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号