问题标签 [geopandas]

我是python的初学者。在提问之前，对不起我的英语不好。

我有两个 shapefile（建筑物 - 多面体，城市 - 多面体）我试图将这两个相交，因为我想获取 city1 的建筑物位置信息。

我将向您展示我使用的两个不同的代码。

另一个是

我认为这两者的不同之处在于上面的类型相交sg_g(shapely.geometry.polygon.Polygon)并且sb(geopandas.geodataframe.GeoDataFrame)

另一个是sg_g(geopandas.geodataframe.GeoDataFrame)和sb(geopandas.geodataframe.GeoDataFrame)

第一个的输出都是假的（错误的），但第二个是正确的

1. 你能告诉我两者的区别吗？

2. 有什么办法可以一次性转换(geopandas.geodataframe.GeoDataFrame)吗(shapely.geometry.polygon.Polygon)？

在 centos 7.3 上，尝试安装 shapely 或 geopandas 需要访问https://github.com/libgeos/libgeos libgeos。

试图通过安装这个sudo yum install libgeos-dev告诉我这个包不可用。看来我缺少存储库。到目前为止，我一直无法找到一个有效的，因为http://trac.osgeo.org/geos指向https://yum.postgresql.org/repopackages.php#pg96的 RPM 但仍然在 rpm -Uvh 之后https://download.postgresql.org/pub/repos/yum/9.6/redhat/rhel-7-x86_64/pgdg-centos96-9.6-3.noarch.rpm找不到安装候选。

我有一个包含 GeoJSON 对象的列表，如下所示：

geometry关键，如何将 转换为coordinates传统的经纬度，最好是在 Python 中？

非常感谢。

我使用geopandas'sto_file()方法将 shapefile 读入 geopandas 对象。shapefile 有一个.prj带有 ESRI WKT 样式投影信息的有效文件：

这样，创建的地理数据框将 crs 属性设置为字典，与 proj4 字符串或 epsg 代码相比，我发现它很难使用：

Geopandas投影文档清楚地表明该.crs方法接受许多不同形式的 crs 信息（epsg 代码、字典、proj4 字符串……），但在将 geopandas 写入 shapefile 时似乎无法控制所需的格式。

问题： 有没有办法指定所需的 crs 格式或任何内置方法来在 crs 属性的不同格式之间切换？

我已将 shapefile 读入 aGeoDataFrame并对其进行了一些修改：

但是，我还想osgeo.ogr在其上应用模块的一些功能：

问题：有什么方法可以直接使用或转换目前以 GeoDataFrame 形式存在的内存中的 shapefile osgeo.ogr.DataSource？

到目前为止，我这样做的方法是将 GeoDataFrame 保存到文件中to_file()，然后osgeo.ogr.Open()再次保存，但这对我来说似乎有点多余。

我正在使用 jupyter notebook 和 Anaconda 2 环境来运行脚本。我最近安装了 geopandas 包。geopandas 包在 conda 列表下列出：geopandas 0.2.1 py27_3 conda-forge

当我尝试导入 geopandas 时，会出现以下错误：

谢谢你

当我将数据框保存为带有几何形状的 shapefile 时，出现以下错误。

我无法找出这个错误的含义。

我试图找到两组多边形的交点，其中一组是从人口普查数据产品（poly1）中选择的一些人口普查区域，另一组由我自己定义的多边形（poly2）。我有兴趣查看哪些人口普查区域落入我的每个多边形（交叉点）中，并从那里计算相交的人口普查区域的面积等等。

我使用了以下代码：

这适用于涉及的人口普查区域不多的情况。我精心挑选了一些案例（5000~40,000 个人口普查区）并在 ArcMap 中进行了交集计算并交叉检查了结果。似乎 geopandas.overlay 给出的结果不太准确：有时给定的“相交”人口普查区不应该算作相交；有时缺少应该相交的人口普查区域。

此外，如果我只从我的自定义多边形 (poly2) 中选择一个多边形并选择它应该与之相交的正确人口普查区域，geopandas.overlay 可以正常工作。

我不知道为什么会这样。似乎只有关于 geopandas 覆盖速度的讨论（https://github.com/geopandas/geopandas/issues/400），但我还没有发现任何关于准确性的信息。我错过了什么吗？

我现在正在尝试在我的 Mac OS X 上安装 geopandas。

我做了什么：

•在此处下载并安装 Anaconda for Python 2.7 。图形安装程序。

•python --version返回以下结果：Python 2.7.13 :: Anaconda 4.3.0 (x86_64)

• 使用这条线安装geopandas conda install -c conda-forge geopandas：。据我了解，此命令还安装了所有其他必需的软件包。这是安装的所有命令的列表：

• geopandas 在我跑步时出现conda list

• 我在test.py 中写了一些代码：

问题：

• 当我在终端中运行执行python test.py时，我收到以下错误：

• 我尝试的第二种方法是通过 PyCharm。我已将项目解释器设置为 python 2.7。geopandas 没有出现在包列表中，因此 PyCharm 甚至在此之前阻塞：

希望任何人都知道如何解决这两个问题。

更新 1

安装 Anaconda 后brew doctor，现在返回：

似乎与这个问题有关。我该如何解决这些问题以使我的东西正常工作？

更新 2

我们可以忘记 PyCharm 问题，结果不知何故我有另一个 python 2.7，这就是他建议的那个。我已将其更改为 Anaconda python 2.7，现在它给出了与终端输出相同的错误。

我现在正在处理shp文件，但遇到了投影问题。让我在下面给你我的代码。

将 pandas 作为 pd 导入 geopandas 作为 gpd 从 geopandas 导入 GeoSeries、GeoDataFrame 导入 o​​s

这是我的数据构建

当我执行时print(base_map.crs), print(Aelly_to_crs.crs)，我得到如下结果。

上图是执行Aelly.plot(ax=base_map.plot()) 的结果，可以看到两张图不匹配。我该如何解决这个问题？

- - - - - -编辑

我想要的输出如下图。

这是我的数据链接：http: //blog.naver.com/khm2963/220929301892

下面的图片是下载我的苍蝇的过程