On 12/8/2010 9:48 AM, [email protected] wrote:
> Just a thought since I'm not doing spatial statistics.
>
> If you can create (integer) labels that assigns each point to a zone,
> then you can treat it essentially as a 1d grouped data, and you could
> use np.bincount to calculate some statistics, or alternatively
> scipy.ndimage.measurements for some additional statistics.
>
> This would avoid any python loop, but require a full label array.
Josef,
The measurements module did the trick; thanks for the pointer. I just
stumbled across a very similar thread on the scipy listserv that you
answered basically the same question with some nice code (sorry for the
redundancy):
http://mail.scipy.org/pipermail/scipy-user/2009-February/019850.html
BTW, the OP on that thread (Jose Gomez-Dans) has a script out there for
doing just this type of operation that I was after:
http://sites.google.com/site/spatialpython/zonal-statistics
thanks, matt
_______________________________________________
NumPy-Discussion mailing list
[email protected]
http://mail.scipy.org/mailman/listinfo/numpy-discussion
