Question

我对地那达斯说了新鲜事,正在对下面的roid名单进行混淆。我从一个规模更大/更为复杂的CSV档案中删除了这些内容。现在我只是试图将这些要点作为自己的地理数据框架加以划分。

对于有兴趣的人来说,这些要点是密歇根州的温室气体排放源。密执安登基图画的地块只是罚款。

The first five lines of the data:

     Centroid
0   { Geometry : [-83.23064676125937, 42.33149724073557],  SRID : 4326}
1   { Geometry : [-82.4951, 42.7741],  SRID : 4326}
2   { Geometry : [-83.27446123221068, 42.54121682538758],  SRID : 4326}
3   { Geometry : [-86.2007, 42.9103],  SRID : 4326}
4   { Geometry : [-82.4719, 42.7642],  SRID : 4326}

我试图这样做:

import pandas as pd
import geopandas as gpd

mi_sources_centroids = gpd.GeoDataFrame(mi_sources_gdf["Centroid"])
fig, ax = plt.subplots(figsize=(13,13))
mi_sources_centroids.plot(ax=ax, color="red")

但总是存在以下错误:“TypeError:没有数字数据可到地段”。

任何建议都会受到高度赞赏。

Answer 1

请尝试:

import pandas as pd
import geopandas as gpd
from shapely.geometry import Point
import matplotlib.pyplot as plt

# Sample data
data = {
     Centroid : [
        { Geometry : [-83.23064676125937, 42.33149724073557],  SRID : 4326},
        { Geometry : [-82.4951, 42.7741],  SRID : 4326},
        { Geometry : [-83.27446123221068, 42.54121682538758],  SRID : 4326},
        { Geometry : [-86.2007, 42.9103],  SRID : 4326},
        { Geometry : [-82.4719, 42.7642],  SRID : 4326}
    ]
}

# Convert to DataFrame
df = pd.DataFrame(data)

# Extract the coordinates and create Point objects
points = [Point(xy) for xy in df[ Centroid ].apply(lambda x: x[ Geometry ])]

# Create a GeoDataFrame
gdf = gpd.GeoDataFrame(geometry=points, crs="EPSG:4326")

# Plot
fig, ax = plt.subplots(figsize=(13,13))
gdf.plot(ax=ax, color="red")

plt.show()

友情链接