python读取一个shp面文件，计算面上每个角的夹角度数，单位为角度

可以使用Python中的geopandas库来读取shp面文件，并使用shapely库来计算每个角的夹角度数。具体实现方法如下：

import geopandas as gpd
from shapely.geometry import Polygon

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

# 计算每个面的角度数
for i, row in gdf.iterrows():
    poly = row['geometry']
    angles = []
    for j in range(len(poly.exterior.coords) - 2):
        p1, p2, p3 = poly.exterior.coords[j:j+3]
        angle = Polygon([p1, p2, p3]).exterior.angles[-1]
        angles.append(angle)
    print(f"面{i+1}的角度数为：{angles}")

其中，`gdf.iterrows()`遍历每个面，`poly.exterior.coords`获取面的外环坐标，`Polygon([p1, p2, p3]).exterior.angles[-1]`计算三个点构成的角的夹角度数。