Comment dessiner des rectangles sur un fond de carte

Je suis à la recherche d'un moyen de tracer des rectangles pleins sur un Fond de carte. Je pourrais facilement tracer le rectangle à bords à l'aide de la drawgreatcircle méthode, mais je ne peux pas trouver un moyen de remplir ces rectangles (en spécifiant la couleur et alpha).

source d'informationauteur andreas-h

  1. 27

    Vous pouvez ajouter un matplotlib.les patchs.Polygone() directement à votre axes. La question est de savoir si vous voulez que votre rectangles défini le tracé de coordonnées (lignes droites sur la parcelle) ou en coordonnées cartographiques (grands cercles sur la parcelle). De toute façon, vous spécifiez les sommets de coordonnées de carte et ensuite les transformer pour le tracé des coordonnées en appelant le Fond de carte instance (m() dans l'exemple ci-dessous), construire un Polygone vous-même, et l'ajouter manuellement les axes d'être rendu.

    Pour les rectangles définis dans la parcelle coordonnées, voici un exemple:

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

def draw_screen_poly( lats, lons, m):
    x, y = m( lons, lats )
    xy = zip(x,y)
    poly = Polygon( xy, facecolor='red', alpha=0.4 )
    plt.gca().add_patch(poly)

lats = [ -30, 30, 30, -30 ]
lons = [ -50, -50, 50, 50 ]

m = Basemap(projection='sinu',lon_0=0)
m.drawcoastlines()
m.drawmapboundary()
draw_screen_poly( lats, lons, m )

plt.show()

    Pour les rectangles définis dans les coordonnées de carte, utiliser la même approche, mais interpoler votre ligne de la carte, à l'espace avant de se transformer pour le tracé des coordonnées. Pour chaque segment de ligne, vous aurez à faire:

    lats = np.linspace( lat0, lat1, resolution )
lons = np.linspace( lon0, lon1, resolution )

    Puis de les transformer ces coordonnées pour le tracé des coordonnées (comme ci-dessus, avec m()) et à nouveau de créer un Polygone avec l'intrigue coordonnées.

  2. 3

    Similaire réponse ci-dessus, mais plus de base code: 

    from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon

map = Basemap(projection='cyl')

map.drawmapboundary(fill_color='aqua')
map.fillcontinents(color='coral',lake_color='aqua')
map.drawcoastlines()

x1,y1 = map(-25,-25)
x2,y2 = map(-25,25)
x3,y3 = map(25,25)
x4,y4 = map(25,-25)
poly = Polygon([(x1,y1),(x2,y2),(x3,y3),(x4,y4)],facecolor='red',edgecolor='green',linewidth=3)
plt.gca().add_patch(poly)

plt.show()
  3. 0

    L'aide d'Andrew réponse, je reçois le message d'erreur 

    TypeError: len() of unsized object.

    Cependant, coulée dans le zip pour une liste des correctifs.

    Code complet:

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

def draw_screen_poly( lats, lons, m):
    x, y = m( lons, lats )
    xy = zip(x,y)
    poly = Polygon( list(xy), facecolor='red', alpha=0.4 )
    plt.gca().add_patch(poly)

lats = [ -30, 30, 30, -30 ]
lons = [ -50, -50, 50, 50 ]

m = Basemap(projection='sinu',lon_0=0)
m.drawcoastlines()
m.drawmapboundary()
draw_screen_poly( lats, lons, m )

plt.show()