It's simpler in numarray:
> i,j = where(z<0.5)
> For numeric it's more work:
> ind = where(z.flat < 0.5) i = ind//z.shape[1] j = ind %
> z.shape[1]
> One could turn this into a function for Numeric (to handle
> all dimensionalities, etc)
I think you meant to use "nonzero" rather than "where" for Numeric...
from matplotlib.numerix import where, nonzero
from matplotlib.numerix.mlab import rand
z = rand(4,5)
# It's simpler in numarray:
i,j = where(z<0.5)
print i,j
# For numeric it's more work:
ind = nonzero(z.flat < 0.5)
i = ind//z.shape[1]
j = ind % z.shape[1]
print i,j
I also think wrapping this functionality into the find function in mpl
is a good idea... I frequently miss it.
JDH