<div dir="ltr">An extra note here:<div><br></div><div style>One of the great things about numpy (as opposed, to say, MATLAB), is array broadcasting  Thus you generally don't need meshgrid -- why carry all that extra repetitive data around. So here's a version that uses broadcasting instead</div>

<div style><br></div><div style><div><br></div><div style>radius = 5<br></div><div><br></div><div># Center at 0, 0</div><div><br></div><div>x = np.linspace(-5,5,11)</div><div style><br></div><div style># y is re-shaped to be a column vector...</div>

<div>y = np.linspace(-8,8,17).reshape((-1,1))</div><div><br></div><div style># when you add x and y you get the 2-d array for the results...</div><div><br></div><div>out_circle = (x * x + y * y < radius**2)</div><div>
<br>
</div><div>print out_circle</div><div><br></div><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Feb 12, 2014 at 2:43 AM, Wolfgang Draxinger <span dir="ltr"><<a href="mailto:Wolfgang.Draxinger@physik.uni-muenchen.de" target="_blank">Wolfgang.Draxinger@physik.uni-muenchen.de</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On Tue, 11 Feb 2014 22:16:46 +0100<br>
Daπid <<a href="mailto:davidmenhur@gmail.com">davidmenhur@gmail.com</a>> wrote:<br>
<br>
> Here it is an example:<br>
><br>
> import numpy as np<br>
> import pylab as plt<br>
><br>
> N = 100<br>
> M = 200<br>
> x, y = np.meshgrid(np.arange(N), np.arange(M))<br>
><br>
> # Center at 40, 70, radius 20<br>
> x -= 40<br>
> y -= 70<br>
> out_circle = (x * x + y * y < 20**2)<br>
<br>
Neat.<br>
<br>
> Note that I have avoided taking the costly square root of each<br>
> element by just taking the square of the radius. It can also be<br>
> generalised to ellipses or rectangles, if you need it.<br>
<br>
Oops, why did I take the "long" tour anyway?<br>
<br>
> Also, don't use 0 as a False value, and don't force it to be a 0.<br>
> Instead, use "if not ring:"<br>
<br>
'ring' is not a boolean, but a numeric index and the test is if the<br>
index is nonzero. I could have also written 'if ring > 0:' and swapped<br>
the clauses. Yes I know that for a value 0 zero it will test as not,<br>
but in this case I wanted to have it written down, that this tests for<br>
a certain numerical value.<br>
<br>
<br>
Thanks,<br>
<br>
Wolfgang<br>
_______________________________________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org">NumPy-Discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><br>Christopher Barker, Ph.D.<br>Oceanographer<br><br>Emergency Response Division<br>NOAA/NOS/OR&R            (206) 526-6959   voice<br>7600 Sand Point Way NE   (206) 526-6329   fax<br>

Seattle, WA  98115       (206) 526-6317   main reception<br><br><a href="mailto:Chris.Barker@noaa.gov" target="_blank">Chris.Barker@noaa.gov</a>
</div>