Intentando extraer xyz de un shapefile multipoligonal

Question

Quiero extraer las x, y y z de los vértices de los polígonos del shapefile donde la geometría es multipoligonal. Una vez extraídos me gustaría tenerlos en un marco de datos geopandas. He visto un par de enfoques sin embargo, el problema parece ser que cuando el shapefile se importa no puedo conseguirlo a una característica que tiene un atributo Z. Por ejemplo, lo siguiente funciona para obtener la x y la y solamente, usando un shapefile que fue generado desde ArcGIS Pro convirtiendo el formato multipolígono a sus partes individuales.

import geopandas as gpd

gdf = gpd.read_file("C:/FG3D_2019/code/Collada_Geospatial/Boston_Mesh_Evaluation/Boston_Bldg1_SingleParts.shp")
print('gdf type', type(gdf))
geom_uu = gdf.unary_union
# Extract the point values that define the perimeter of the polygon
x, y = geom_uu.exterior.coords.xy

Lo que me gustaría es obtener las coordenadas z del campo de la geometría así que esencialmente; x, y, z = geom_uu.exterior.coords.xyz .

¿Alguien ha tenido suerte con esto?

Answer 1

from shapely.geometry import Polygon
polygon = Polygon([[0.0, 0.0, 2.0], [1.0, 0.0, 3.0], [1.0, 1.0, 4.O], [0.0, 1.0, 5.0]])
x,y,z = zip(*polygon.exterior.coords)
print(x)
(0.0, 1.0, 1.0, 0.0, 0.0)
print(y)
(0.0, 0.0, 1.0, 1.0, 0.0)
print(z)
(2.0, 3.0, 4.0, 5.0, 2.0)