Re: [Numpy-discussion] moving window in numpy
Matthew Perry wrote: Hi all, I'm not sure if my terminology is familiar but I'm trying to do a moving window analysis (ie a spatial filter or kernel) on a 2-D array representing elevation. For example, a 3x3 window centered on each cell is used to calculate the derivate slope of that cell. Can this easily be implemented using numpy? If I understand you correctly, scipy.ndimage does exactly what you ask for. See http://www.scipy.org/SciPyPackages/Ndimage. Cheers, VS. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion
Re: [Numpy-discussion] moving window in numpy
On Montag 26 November 2007, Matthew Perry wrote: Hi all, I'm not sure if my terminology is familiar but I'm trying to do a moving window analysis (ie a spatial filter or kernel) on a 2-D array representing elevation. For example, a 3x3 window centered on each cell is used to calculate the derivate slope of that cell. Did you have a look at ndimage? http://www.scipy.org/SciPyPackages/Ndimage Currently I have tried implementing in pure python loops (too slow) You could probably have a faster, pure python version with some indexing tricks. E.g. if your gradient filter is static (i.e. the weights do not adapt to the local image content), you can apply it to whole images that contain shifted image contents. This way, you're using the fast C loops of NumPy. -- Ciao, / /.o. /--/ ..o / / ANS ooo signature.asc Description: This is a digitally signed message part. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion
Re: [Numpy-discussion] moving window in numpy
Hi Matthew On Sun, Nov 25, 2007 at 05:00:48PM -0800, Matthew Perry wrote: I'm not sure if my terminology is familiar but I'm trying to do a moving window analysis (ie a spatial filter or kernel) on a 2-D array representing elevation. For example, a 3x3 window centered on each cell is used to calculate the derivate slope of that cell. Take a look at Anne Archibald's post: http://projects.scipy.org/pipermail/numpy-discussion/2006-November/024760.html and its attachment here: http://projects.scipy.org/pipermail/numpy-discussion/attachments/20061120/d057cb0b/attachment.py Regards Stéfan ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion