geopandas如何提取shp数据的经纬度到属性

可以使用GeoPandas库中的`centroid`属性来获取多边形的中心点，然后使用`geometry`属性中的`x`和`y`属性获取中心点的经度和纬度。将这些值存储到新的列中，可以使用`GeoDataFrame`的`assign`方法。

以下是一个例子：

import geopandas as gpd

# 读取shp文件
gdf = gpd.read_file('file.shp')

# 计算中心点经纬度并添加到新列中
gdf = gdf.assign(latitude=gdf.centroid.y, longitude=gdf.centroid.x)

# 保存结果到新的shp文件
gdf.to_file('new_file.shp')

在这个例子中，`gdf.centroid.y`和`gdf.centroid.x`分别提取了中心点的纬度和经度，并将它们存储在名为`latitude`和`longitude`的新列中。最后，使用`to_file`方法将结果保存到新的shp文件中。

相关问题

geopandas如何将所有shp点数据的经纬度到属性

geopandas 可以通过以下步骤将所有shp点数据的经纬度加入到属性中：

1. 读取shp文件，使用geopandas.read_file()函数。

2. 使用geopandas中的apply()方法，遍历所有要素，将要素的几何坐标点转换为经纬度坐标。

3. 将经纬度坐标添加到每个要素的属性中，使用geopandas中的set_geometry()方法。

以下是示例代码：

import geopandas as gpd

# 读取shp文件
gdf = gpd.read_file('path/to/shapefile.shp')

# 定义函数，将几何坐标点转换为经纬度坐标
def get_lon_lat(point):
    lon, lat = point.x, point.y
    return lon, lat

# 将经纬度坐标添加到属性中
gdf['longitude'], gdf['latitude'] = zip(*gdf.geometry.apply(get_lon_lat))

# 将经纬度坐标设置为几何坐标
gdf = gdf.set_geometry(gpd.points_from_xy(gdf.longitude, gdf.latitude))

# 保存修改后的shp文件
gdf.to_file('path/to/new_shapefile.shp')

geopandas如何将点要素的geometry列的经纬度提取到属性

可以使用geopandas中的apply方法和shapely.geometry中的xy属性来提取点要素的经纬度并添加到属性中。具体步骤如下：

1. 读取点要素数据，并将geometry列转换为Point类型：

import geopandas as gpd

# 读取点要素数据
gdf = gpd.read_file('points.shp')

# 将geometry列转换为Point类型
gdf['geometry'] = gdf['geometry'].apply(lambda x: x.centroid)

2. 创建经度和纬度列，并使用apply方法和xy属性提取经纬度：

# 创建经度和纬度列
gdf['longitude'] = None
gdf['latitude'] = None

# 提取经纬度
gdf['longitude'], gdf['latitude'] = gdf['geometry'].apply(lambda x: x.xy[0][0]), gdf['geometry'].apply(lambda x: x.xy[1][0])

3. 查看添加经纬度后的属性表：

print(gdf.head())

输出结果：

   id                 geometry  longitude   latitude
0   1  POINT (116.397 39.907)  116.39722  39.90693
1   2  POINT (116.398 39.907)  116.39822  39.90693
2   3  POINT (116.399 39.907)  116.39922  39.90693
3   4  POINT (116.400 39.907)  116.40022  39.90693
4   5  POINT (116.401 39.907)  116.40122  39.90693