Maskenbereich außerhalb der importierten Shapefile (Grundkarte/Matplotlib)

Question

Ich zeichne Daten auf einer Grundkarte der Ostküste der USA und Kanada durch Matplotlib. Zusätzlich zu der Basisschicht (einem gefüllten Konturdiagramm) überlagerte ich ein Shapefile dieser Fokusregion auf den Daten unter Verwendung von Matplotlibs readshapefile-Werkzeug.

Ich frage mich, wie ich alle gerasterten Daten außerhalb der Shape-Datei maskieren kann. Natürlich kann ich durch Matplotlib einen Maskocean-Befehl ausführen, aber ich würde immer noch die gefüllten Konturen westlich des St. Lawrence haben. Kann jemand jetzt wie das? Ich hatte nicht viel Glück, online zu suchen.

def make_map(lon,lat,param): 
    fig, ax = plt.subplots() 
    ax.axis('off') 
    x1 = -83. 
    x2 = -57. 
    y1 = 37. 
    y2 = 50. 
    projection='merc' 
    resolution='h' 
    m = Basemap(projection=projection, llcrnrlat=y1, urcrnrlat=y2, llcrnrlon=x1, 
       urcrnrlon=x2, resolution=resolution) 
    x,y = m((lon-360.),lat) 
    m.ax = ax 
    my_cmap = cm.get_cmap('coolwarm') 
    pp = m.contourf(x, y, param, 30, cmap=my_cmap, extend='both') 
    m.drawmapscale(-67, 39.5, -70, 43.5, 500, fontsize=8, barstyle='fancy') 
    return fig, m, x, y 

def drawstates(ax, shapefile='../StateProv_UTMrp'): 
     shp = m.readshapefile(shapefile, 'states',zorder = 1, drawbounds=True) 
     for nshape, seg in enumerate(m.states): 
      poly = Polygon(seg, facecolor='w',alpha=0.0, edgecolor='k') 
      ax.add_patch(poly) 

fig, m, x, y = make_map(lon, lat, param) 
drawstates(m.ax) 

Answer 1

Ich denke, dieser Artikel, den ich gerade gefunden habe, könnte Ihnen etwas Hilfe geben. Aber ich bin mir nicht sicher, dass es eine vollständige Antwort ist.

http://basemaptutorial.readthedocs.org/en/latest/clip.html