Re: [Matplotlib-users] display a filled lat/lon basemap rectangle.

[Jeff W. <js...@fa...>]

Yeates, Mathew C (388D) wrote:
>
>  
>
> Hi
>
> What is the simplest way to fill  in a 1 degree by 1 degree rectangle 
> on a basemap projection?
>
>  
>
> Mathew
>

Mathew:  Try this (for a 10x10 rectangle, but you get the idea)

from matplotlib.patches import Polygon
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
map = Basemap(projection='moll',lon_0=0)
x1,y1 = map(-10,-10)
x2,y2 = map(-10,10)
x3,y3 = map(10,10)
x4,y4 = map(10,-10)
p = Polygon([(x1,y1),(x2,y2),(x3,y3),(x4,y4)],\
    facecolor='red',edgecolor='blue',linewidth=2)
plt.gca().add_patch(p)
map.drawcoastlines()
map.drawmapboundary()
plt.show()

-Jeff
>