Mathematica 快速 2D 分箱算法

我在 Mathematica 中开发适当快速的分箱算法时遇到一些麻烦。我有一个很大的（~100k 元素）数据集，其形式为 T={{x1,y1,z1},{x2,y2,z2},....} 我想将其分入大约 100x100 个 bin 的 2D 数组中，bin 值由每个 bin 中的 Z 值之和给出。

目前，我正在迭代表的每个元素，使用 Select 根据 bin 边界列表挑选出它应该位于哪个 bin，并将 z 值添加到占用该 bin 的值列表中。最后，我将 Total 映射到 bin 列表上，对它们的内容求和（我这样做是因为我有时想做其他事情，比如最大化）。

我曾尝试使用 Gather 和其他此类函数来执行此操作，但上述方法速度快得离谱，尽管也许我使用 Gather 的效果不佳。无论如何，按照我的方法进行排序仍然需要几分钟，我觉得 Mathematica 可以做得更好。有人有一个方便的高效算法吗？

这是一种基于 Szabolcs 帖子的方法，速度大约快一个数量级。

data = RandomReal[5, {500000, 3}];
(*500k values*)
zvalues = data[[All, 3]];

epsilon = 1*^-10;(*prevent 101 index*)
(*rescale and round (x,y) coordinates to index pairs in the 1..100 range*)
indexes = 1 + Floor[(1 - epsilon) 100 Rescale[data[[All, {1, 2}]]]];

res2 = Module[{gb = GatherBy[Transpose[{indexes, zvalues}], First]}, 
    SparseArray[
     gb[[All, 1, 1]] -> 
      Total[gb[[All, All, 2]], {2}]]]; // AbsoluteTiming

给出大约 {2.012217, Null}

AbsoluteTiming[
 System`SetSystemOptions[ 
  "SparseArrayOptions" -> {"TreatRepeatedEntries" -> 1}];
 res3 = SparseArray[indexes -> zvalues];
 System`SetSystemOptions[ 
  "SparseArrayOptions" -> {"TreatRepeatedEntries" -> 0}];
 ]

给出大约 {0.195228, Null}

res3 == res2
True

"TreatRepeatedEntries" -> 1 添加重复位置。