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

可以使用geopandas库中的"centroid"函数来获取多边形的中心点经纬度，然后将其添加到属性表中。以下是示例代码：

import geopandas as gpd

# 读取shp文件
data = gpd.read_file('path/to/your/file.shp')

# 获取中心点经纬度
data['lon'] = data.centroid.x
data['lat'] = data.centroid.y

# 保存修改后的文件
data.to_file('path/to/your/new_file.shp')

这将在原始属性表中添加"lon"和"lat"列，并将其填充为每个多边形的中心点经纬度。最后，您可以将修改后的数据保存为新的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