2017-10-08 15 views
1

Je travaille avec un fichier de formes. Je n'ai aucun problème à le lire, à le tracer et à rendre la carte plutôt jolie. Cependant, lorsque je le trace (après l'avoir reprojecté au bon EPSG en utilisant QGIS), les bords sont tous dentelés (comme indiqué ci-dessous). Y a-t-il un moyen de lisser en utilisant Python? enter image description hereLissage d'une sortie de fichier de formes - Basemap python

import matplotlib.pyplot as plt 
from mpl_toolkits.basemap import Basemap 
from matplotlib.patches import Polygon 
from matplotlib.collections import PatchCollection 
import numpy as np 

#insert code for basemap setup m = Basemap(...) 
m.arcgisimage(service = 'ESRI_StreetMap_World_2D', xpixels = 1000, verbose = True) 
states_info = m.readshapefile('shapefiles/states', 'states') 
spc_info = m.readshapefile('shapefiles/corrected_epsg', 'spc', drawbounds = False) 

patches = [] 
ax = plt.gca() 

for info, shape in zip(m.spc_info, m.spc): 
    x, y = zip(*shape) 
    if info['DN'] == 2: 
     color = '#80c580' 
     zorder = 2 
     patches.append(Polygon(np.array(shape), True)) 
    if info['DN'] == 5: 
     color = '#f7f780' 
     zorder = 3 
     patches.append(Polygon(np.array(shape), True)) 
    ax.add_collection(PatchCollection(patches, facecolor= color, zorder=zorder, alpha = 0.7)) 

Source pour ces shapefiles.

Répondre

1

This question Les réponses expliquent comment le paquet Shapely a une méthode Simplify basée sur l'algorithme de Douglas-Puecker.