Quantile/Median/2D binning in Python

знаете ли вы быстрое/элегантное решение Python/Scipy/Numpy для следующей проблемы :У вас есть набор координат x, y со связанными значениями w (all 1D arrays ). Теперь бины x и y помещаются в двумерную сетку (размера BINSxBINS )и вычисляются квантили (, такие как медиана )значений w для каждого бина, что в конце должно привести к двумерному массиву BINSxBINS с необходимые квантили.

Это легко сделать с помощью вложенного цикла, но я уверен, что есть более элегантное решение.

Спасибо, Марка