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

2010-04-07 Thread Yeates, Mathew C (388D)
Hi What is the simplest way to fill in a 1 degree by 1 degree rectangle on a basemap projection? Mathew -- Download Intel#174; Parallel Studio Eval Try the new software tools for yourself. Speed compiling, find bugs

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

2010-04-07 Thread Jeff Whitaker
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

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

2010-04-07 Thread Mathew Yeates
I think this will only work with some projections but not all. I looked at the code for tissot. It's pretty hairy but it almost does what I want. (It draws projected circles instead of projected rectangles. On Wed, Apr 7, 2010 at 1:54 PM, Jeff Whitaker jsw...@fastmail.fm wrote: Yeates, Mathew

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

2010-04-07 Thread Jeff Whitaker
Mathew Yeates wrote: I think this will only work with some projections but not all. I looked at the code for tissot. It's pretty hairy but it almost does what I want. (It draws projected circles instead of projected rectangles. Mathew: You said you wanted a NxN degree polygon - that's what

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

2010-04-07 Thread Mathew Yeates
lets say I want to shade the area with lat/lon corners 34.-117 and 35,-116 but my map was created with projection='aeqd' The shade area will not be a rectangle. In fact the edges will be curved. See the basemap code for tissot. I think every point on the boundary of the lat/lon box has to

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

2010-04-07 Thread Jeff Whitaker
Mathew Yeates wrote: lets say I want to shade the area with lat/lon corners 34.-117 and 35,-116 but my map was created with projection='aeqd' The shade area will not be a rectangle. In fact the edges will be curved. See the basemap code for tissot. I think every point on the boundary