小编M.O*_*.O.的帖子

有人可以解释一下这不起作用吗？

有问题的形状文件在这里，在代码中它的读作是

shp=gpd.read_file("Microdatos_Censo_2017_Manzana/Microdatos_Censo_2017_Manzana.shp")\nshp.crs="epsg:4326"\nbreakpoint()\nshp=shp.to_crs(epsg=3857)##\xc2\xa0Error here\n

我只是不明白发生了什么。我有Python 3.8.5，geopandas 0.8.1，pyproj 2.6.1.post1。不确定哪些其他软件包对于了解其版本很重要。

谢谢！

编辑：

1.- 修复我弄错的 shapefile 链接。

2.- 它与作为重复项发布的问题不同，因为正如您在图像上看到的，打印语句返回shp.crs正确的 crs 信息，而不是“无”。我定义了一个 crs 但to_crs不起作用。

好的,如果您使用它,您将获得列表的可能组合

stuff = [1, 2, 3]
for L in range(0, len(stuff)+1):
    for subset in itertools.combinations(stuff, L):
        print(subset)

结果是

()
(1)
(2)
(3)
(1, 2)
(1, 3)
(2, 3)
(1, 2, 3)

有没有办法改变这个,以便结果

(0, 0, 0)
(1, 0, 0)
(0, 2, 0)
(0, 0, 3)
(1, 2, 0)
(1, 0, 3)
(0, 2, 3)
(1, 2, 3)

或其他代码会这样做？