Re: [Numpy-discussion] [enhancement] sum_angle() and sum_polar()
On Tue, May 29, 2012 at 11:03 AM, Stéfan van der Walt ste...@sun.ac.za wrote: On Mon, May 28, 2012 at 11:53 AM, Travis Oliphant tra...@continuum.io wrote: I could see these functions going into scipy.ndimage but again because they are not necessarily just image processing functions, and the fact that they are so simple, perhaps they are best put into NumPy itself. I'm wondering about the general applicability of these functions. Can anyone suggest some use cases? An example from solid state physics: If you have a spin chain with some long-range interaction and you have the known, dense, coupling matrix J, sum_angle(J, pi/4) gives you a view at the distance dependence of the interaction. -- Robert Jordens. ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
Re: [Numpy-discussion] [enhancement] sum_angle() and sum_polar()
On Tue, May 29, 2012 at 12:06 PM, Charles R Harris charlesr.har...@gmail.com wrote: I'd like to see these functions is scipy somewhere. The function names aren't very descriptive and the one line summaries don't give a very good idea of what they do, so I think those bits could use improvement. Mention of the Hough/Radon transform would help, I had to pull out that connection by reading the code... I'll fix the descriptions. What more descriptive names did you have in mind? -- Robert Jordens. ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
[Numpy-discussion] [enhancement] sum_angle() and sum_polar()
Hi everyone, I am proposing to add the the two following functions to numpy/lib/twodim_base.py: sum_angle() computes the sum of a 2-d array along an angled axis sum_polar() computes the sum of a 2-d array along radial lines or along azimuthal circles https://github.com/numpy/numpy/pull/230 Comments? When I was looking for a solution to these problems of calculating special sums of 2-d arrays I could not find anything and it took me a while to figure out a (hopefully) useful and consistent algorithm. I can see how one would extend these to higher dimensions but that would preclude using bincount() to do the heavy lifting. Looking at some other functions, the doctests might need to be split into real examples and unittests. Best, -- Robert Jordens. ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
