直方图算法中用++实现倍频程快速递增的方法

Question

我正在做一个二维直方图算法与仓储。在我的代码中，我分析了以下几行代码，它们占用了很多时间。有没有办法提高性能，通过矢量化或任何其他方法？

for i = 1 : length(r)
    hist_r(r(i))++;
    hist_g(g(i))++;
    hist_b(b(i))++;
endfor

Answer 1

最后，我使用了：

accumarray(r, 1);

谢谢你的所有建议。

Answer 2

我能想到的最简单的方法是

MyHist += 1;

其中MyHist是NBins x 3数组，其中列表示r、g和b。