重新投影GeoDataframe未正确转换几何体

Question

我正在使用Geopandas通过读取一个已投影到EPSG:3735（俄亥俄州平面南部，美国英尺）的形状文件（.shp）来创建一个geodataframe。现在我想把另一个geodataframe连接到这个。第二个geodataframe是通过读取包含lon，lat格式的x，y坐标的.csv文件创建的。因此，我将第二个地理数据框的crs定义为EPSG:4326（WGS 1984）。要使用空间连接函数，两个geodataframe的crs必须相同，因此我尝试使用To_crs将第一个geodatafame重新投影到EPSG:4326。在重新投影之前，Geodataframe gdf如下所示：

如果我写：

gdf_proj = gdf.to_crs("EPSG:4326")

结果是这样的：

我还尝试将第二个gdf从EPSG:4326转换为EPSG:3735，但第二个gdf中的点几何体也将成为点（inf，inf）。我不确定这里出了什么问题。有人能推荐一种在geodataframe中重新投影多边形的好方法吗？

我刚试过其他的crs。其中一些会起作用（不会变成无穷大），比如37853035。我也试过其他一些国家飞机的电脑，但其中一些工作，而另一些没有。俄亥俄州南部（3735）和北部（3734）不工作。印第安纳州东部（2965）不起作用，但西部（2966）起作用。此外，NAD83地面军事系统（4269）的项目不起作用。

Answer 1

我想我已经明白了。不应使用WGS 1984（4326）初始化.csv文件中的第二个gdf，而应将其设置为NAD 1983（4269）。两个地面军事系统都使用lat和long，所以我想从一个转换到另一个会带来一些麻烦。然后，我可以使用to_crs将gdf从EPSG:4269转换为EPSG:3735，因为3735是使用NAD83的PC。因此，我想当使用Geopandas中的其他地面军事系统将数据从一个地面军事系统直接投影到另一台计算机时，可能会引起一些麻烦。我没有检查源代码，但为了避免类似的问题，我想说，如果目标crs位于另一个地面军事系统中，在尝试使用投射到\u crs时要小心。

Answer 2

您可以检查gdf是否设置了正确的坐标系。您可以使用gdf.crs实现这一点。如果不是EPSG:3735，而是其他东西，这可能会导致你看到的无穷大。

如果这是错误的，调查原因可能是有意义的。但在任何情况下，您都可以使用set\u crs覆盖错误的值：

gdf = gdf.set_crs("EPSG:3735")
gdf_proj = gdf.to_crs("EPSG:4326")