Patches de polígonos de Geopandas em matplotlib Conversão de polígonos

8

Infelizmente, a plotagem de geopandas é extremamente lenta e requer muitos recursos; portanto, eu gostaria de usar o matplotlib para plotagem.

Quando uso o Fiona puro para abrir e ler o shapefile, não tenho problemas para extrair os polígonos como patches do matplotlib, mas agora gostaria de usar como ponto de partida o quadro de dados das geopandas para obter meus polígonos do matplotlib.

Atualmente, estou usando algo como:

with FI.open(df_map_elements, 'r') as layer:
    for element in layer:
        key = int(element['id'])
        if key not in dict_mapindex_mpl_polygon.keys():
            dict_mapindex_mpl_polygon[key]=[]
        for tp in element['geometry']['coordinates']:
            q = np.array(tp)
            polygon = Polygon(q) # matplotlib Polygon NOT Shapely

Para plotar polígonos com matplotlib:

from matplotlib import pyplot as plt
from matplotlib.patches import Polygon
from matplotlib.collections import PatchCollection
Philipp Schwarz
fonte

Respostas:

7

Existe um módulo Python para isso: Descartes (veja Plot shapefile com matplotlib, por exemplo)

from geopandas import GeoDataFrame
test = GeoDataFrame.from_file('poly1.shp')
test.set_index('id', inplace=True)
test.sort()
test['geometry']
testid
0    POLYGON ((1105874.411110075 -6125459.381061088...
1    POLYGON ((1106076.359169902 -6125875.557806003...
2    POLYGON ((1106260.568548799 -6125410.258560049...
3    POLYGON ((1105747.511315724 -6125864.64169466,...
Name: geometry, dtype: object

O tipo da geometria é um polígono bem torneado:

 type(test['geometry'][2])
 shapely.geometry.polygon.Polygon

Agora você pode usar o Descartes para plotar diretamente um polígono bem torneado

import matplotlib.pyplot as plt 
from descartes import PolygonPatch
BLUE = '#6699cc'
poly= test['geometry'][2]
fig = plt.figure() 
ax = fig.gca() 
ax.add_patch(PolygonPatch(poly, fc=BLUE, ec=BLUE, alpha=0.5, zorder=2 ))
ax.axis('scaled')
plt.show()

insira a descrição da imagem aqui

gene
fonte
3

Após a resposta simples e compreensível, criei uma maneira direta de plotar um shp inteiro com o matplotlib. Eu acho que as geopandas devem apenas atualizar sua função de plotagem, porque essa é simples, mas muito mais rápida, incluindo a flexibilidade total do matplotlib - adicionando legenda, título etc.

from descartes import PolygonPatch
import geopandas as gp
import pysal as ps
import numpy as np

# Import libraries for visualization
from matplotlib import pyplot as plt
from matplotlib.patches import Polygon as mpl_Polygon
from matplotlib.collections import PatchCollection

shapefile = 'raw_data/shapefile/yourshapefile.shp'
df_map_elements = gp.GeoDataFrame.from_file(shapefile)

df_map_elements["mpl_polygon"] = np.nan
df_map_elements['mpl_polygon'] = df_map_elements['mpl_polygon'].astype(object)
for self_index, self_row_df in df_map_elements.iterrows():
    m_polygon = self_row_df['geometry']
    poly=[]
    if m_polygon.geom_type == 'MultiPolygon':
        for pol in m_polygon:
            poly.append(PolygonPatch(pol))
    else:
        poly.append(PolygonPatch(m_polygon))
    df_map_elements.set_value(self_index, 'mpl_polygon', poly)

dict_mapindex_mpl_polygon = df_map_elements['mpl_polygon'].to_dict()

E para plotagem:

fig, ax = plt.subplots()
for c_l ,patches in dict_mapindex_mpl_polygon.items():
    p = PatchCollection(patches,color='white',lw=.3,edgecolor='k')
    ax.add_collection(p)
ax.autoscale_view()

plt.show()

insira a descrição da imagem aqui

Philipp Schwarz
fonte