Блог

У меня есть geodataframe, который выглядит следующим образом:

 ID_0    ISO NAME_0  ID_1    NAME_1  ID_2    NAME_2  TYPE_2  ENGTYPE_2   NL_NAME_2   VARNAME_2   geometry    soyb_a  percent percent_sum
1489    33  BRA Brazil  12  Mato Grosso 1490    Nova Mutum  Município   Municipality    0   0   POLYGON ((-56.61388 -12.87704, -56.57753 -12.8...   1078374.8   2.923144    2.923144
1405    33  BRA Brazil  11  Mato Grosso do Sul  1406    Sapezal Município   Municipality    0   0   POLYGON ((-57.82408 -19.11719, -57.78419 -19.0...   1027233.8   2.784516    5.707660
1529    33  BRA Brazil  12  Mato Grosso 1530    Sapezal Município   Municipality    0   0   POLYGON ((-58.92996 -12.64107, -58.93618 -12.6...   1027233.8   2.784516    8.492176
  

Я вижу список точек в столбце «геометрия», но хотел бы извлечь эти точки и поместить их в список. Например, в pandas вы можете сделать что-то вроде df[‘column’].to_list() . Однако при попытке этого я получаю сообщение об ошибке:

 gdf.iloc[0]['geometry'].to_list()

AttributeError: 'Polygon' object has no attribute 'to_list'
  

Есть идеи, как я могу удалить обозначение «Полигон» и буквально просто получить список точек, составляющих этот полигон? Чтобы было ясно, мне не нужны внешние или границы полигона, мне нужны все точки внутри границы.

Это общая функция, которую я использую для проверки содержимого полигона — не уверен, что это именно то, что вы ищете. Я считаю, что полигоны могут иметь произвольную сложность, поэтому может случиться так, что у вас могут быть части внутри частей:

 def listPoints(someGeometry):
    '''List the points in a Polygon in a geometry entry - some polygons are more complex than others, so accommodating for that'''    
    pointList = []
    try:
        #Note: might miss parts within parts with this
        for part in someGeometry:
            x, y = part.exterior.coords.xy
            pointList.append(list(zip(x,y)))
    except:
        try:
            x,y = someGeometry.exterior.coords.xy
            pointList.append(list(zip(x,y)))
        except:
            #this will return the geometry as is, enabling you to see if special handling is required - then modify the function as need be
            pointList.append(someGeometry)
    return pointList
  

Затем применить как лямбда:

 gdf.geometry.apply(lambda x: listPoints(x)).values.tolist()