问题标签 [fiona]

尝试使用时，geopandas.read_file('xyz.shp')我会经常（但不总是）得到一个错误IndexError: list index out of range。它也出现在 Fiona 中，在这两种情况下都是在尝试迭代功能时。（见这里，还有这个 Github 问题）。

有趣的是，这并不总是发生。geopandas 成功加载了一些 shapefile，但通常情况下，它会因 IndexError 而失败。

先前的讨论集中在库依赖上，但它有时在这里起作用的事实使我相信它与 shapefile 本身有关，尽管我不确定是什么。我正在使用新版本的 Geopandas (2.0)。

我正在尝试使用来自该站点的文件：

国家发展计划网站

文件是.e00格式的。当我尝试使用 GeoPandas 打开它们时，我收到一条消息，它们似乎已被压缩。

如果我尝试使用 e00conv 或 AVCE00 解压缩文件，然后用 GeoPandas 打开它们，我会得到一个FionaValueError，表明没有找到数据集。

关于如何打开这些文件以便我可以将它们放入我可以使用的格式的任何建议？

我可以使用加载解压缩的文件，np.fromfile但我所拥有的只是一个向量。

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

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

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

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

fiona在 Python 中使用and处理坐标系时osgeo，似乎有很多方法可以通过导入/导出不同的 crs 格式来定义坐标系，例如：

菲奥娜：

奥斯吉奥：

但是，我注意到，这两个库都允许通过其 EPSG 代码轻松定义 crs，但它们都缺少反函数（将 crs 导出为 ESPG 代码）。

我得到EPSG代码的最接近的是：

但它似乎并不那么可靠，而且其他提议的解决方案似乎也包括大量的调整或至少对Web 服务的依赖。

问题：

1. 有人可以向我展示一种将 CRS 导出为 python 中的 EPSG 代码的简单、直接的方法吗？类似to_epsg()inFiona或ExportToEPSG()in 的东西osgeo？

2. 有人可以解释整个互联网上 EPSG 出口可能性如此短缺的理论背景，尤其是与 EPSG 代码导入的便利性相比。EPSG 代码的全部意义不在于使没有高级地理空间专业知识的人易于识别和使用坐标系吗？它不应该像 CRS 的 ID 一样服务，因此可以轻松检索吗？

我在 ubuntu 14.04 上安装 fiona 有一些问题

我总是得到错误

我已经尝试过：sudo -H pip install fiona以及sudo pip install -U setuptools, python -m pip install fiona, sudo apt-get install python-software-properties,sudo add-apt-repository ppa:ubuntugis/ppa

任何其他想法可以尝试什么？

提前谢谢了。

我正在尝试通过以下方式处理底图中的 shapefile：

虽然代码与另一个 shapefile 配合得很好，但使用这个 shapefile，我得到了错误：

我认为问题在于shapefile可能不干净，但我不知道如何解决这个问题。任何想法？提前谢谢了。

我正在尝试安装geopandas。

我是否做其中之一：

当我尝试使用import geopandasorimport fiona时，我收到以下错误：

关于如何解决这个问题的任何想法？

（顺便说一句，我在 Windows 7 64 位上，以防万一）

编辑 ：

这种行为似乎只发生在 Spyder (v3.1.4) 的 IPython 控制台中。在终端中使用 python3 命令行或 Spyder 的 python 控制台时不会发生这种情况。我仍然希望修复 Spyder 的 IPython 控制台，Spyder 中的 python 控制台真的很慢。

有没有办法使用任何 GDAL/OGR API 或命令行工具来溶解（合并）重叠多边形，同时保持生成的非重叠区域不同？我进行了很多搜索，但找不到任何类似于需要的东西。但是，我认为这个问题还没有解决的可能性很小。

这是我需要的更详细的描述：

• 我的输入由一个单层的形状文件 (ESRI Shapefile) 组成。
• 该层包含无法通过属性区分的多边形。（都具有相同的属性）。
• 他们中的许多人是重叠的，我想得到那些重叠的人的联合。
• 未连接的区域应生成单独的多边形。

这是造成麻烦的最后一点。除了最后一点，我基本上得到了我需要的东西。如果我运行用于溶解形状文件的典型解决方案

我最终得到一个包含所有内容的多边形，即使这些区域没有连接。

更新： 我通过完全放弃 GDAL 解决了这个问题。正如许多消息来源所指出的，通常使用 fiona 和 shapely 处理 shapefile 是一种更好的方法。我在下面发布了我的解决方案。

我正在尝试为输入文件属于某个市区的交叉口创建一个函数，并且使用查询框创建一个输出文件，该输出文件的交叉口仅包含在该查询框中找到的建筑物。

不知道我的代码到目前为止是否正确，但我有两个直接的问题：首先，我不知道如何使用 Fiona 打开输入 shapefile 和我想在输出中生成的新（剪辑）shapefile。我想循环索引列表，选择所需的建筑物，然后将它们写入新文件“out_file”。其次，我收到一个错误：

我有一个.gdb需要使用 Python 过滤的大型地理数据库 ( ) 文件。我的工作机器上有 Python 3.4。

我尝试使用 fiona 来读取 gdb 及其层。我不断收到同样的错误：

gdb = fiona.open('GDB 路径')

错误消息：嵌入的 NUL 字符。

我怎样才能最好地解决这个问题？