代码之家 › 专栏 › 技术社区 › Kernel

如果指定了crs参数,shapefile的geopandas.read_file将出错

shapefile geopandas python-3.x python

Kernel · 技术社区 · 5 月前

所有, 我使用ESRI世界国家通用形状文件,该文件可用 here 使用GeoPandas

shp_file =gpd.read_file('World_Countries/World_Countries_Generalized.shp')
print(shp_file.crs)

我得到的CRS是 EPSG:3857 ,但是一旦我将CRS添加到gpd.read_file中,如下所示

shp_file1 =gpd.read_file('../../Downloads/World_Countries/World_Countries_Generalized.shp',crs='EPSG:3857')

我得到了以下错误 /opt/anaconda3/envs/geo_env/lib/python3.12/site-packages/pyogrio/raw.py:198: RuntimeWarning: driver ESRI Shapefile does not support open option CRS return ogr_read(

你知道我为什么会遇到这个错误吗?这是否意味着文件读取不正确?

谢谢

1 回复 | 直到 5 月前

Pieter 5 月前

自geopandas 1.0以来,默认情况下使用另一个更快的底层库来读取geopandas中的文件: pyogrio 。更多信息可以在这里找到: fiona vs pyogrio 当使用这个新库时 crs 参数不受支持。

对于这种特定情况,最简单的解决方案是删除 crs='EPSG:3857' 参数原样无论如何都是无用的,因为crs已经被正确读取了?

shp_file1 = gpd.read_file('../../Downloads/World_Countries/World_Countries_Generalized.shp')

如果你想读取一个没有.prj文件或错误的shapefile,你可以使用 set_crs 功能 GeoDataFrame 设定或推翻 crs 读取文件后:

shp_file2 = gpd.read_file('wrong_crs.shp').set_crs(crs='EPSG:3857', allow_override=True)

推荐文章

ÎÎÎ½Î· ÎÎ®Î¹Î½Î¿Ï · Python lxml.html语法错误:使用lxml find时XPATH的谓词无效

3 月前

Cam · Pandas列表日期到日期时间

3 月前

RASEL MAHMUD · 为什么以及如何在is_even()函数内的IF条件中递归X变量在满足0后递增?[副本]

4 月前

jjkennedy · Pandas文本文件导入:当每个文件中存在多个表时,自动选择1个表

4 月前

LMC · Numpy数组布尔索引以获取包含元素

4 月前

vr8ce · 非成对标记中特定字符的正则表达式

5 月前

Kernel · 如果指定了crs参数,shapefile的geopandas.read_file将出错

5 月前

ShaAnder · 为什么sqllachemy返回的是类而不是字符串

5 月前

sixtytrees · detectron2软件包未安装(没有名为“torch”的模块),但我安装了torch

5 月前

Pernoctador · Python映射可以复制吗?我需要参考地图

5 月前