从数字矩阵创建直方图的矩阵

Question

我有一个矩阵的值从1到32，我需要创建一个直方图的矩阵，使得（i，j）直方图将由零组成，除了它将具有1在（i，j）上。我可以用循环来做到这一点，但是有没有一种方法可以将此操作矢量化？

这里是循环的变种：

hist = zeros(size(a,1),size(a,2), 32); 
for i = 1 : size(a, 1) 
    for j = 1 : size(a, 2) 
     hist(i,j,a(i,j)) = 1 
    end 
end 

Answer 1

你可以做到这一点与accumarray。

[y x]=ndgrid(1:size(a,1),1:size(a,2)); 
hist=accumarray([y(:) x(:) a(:)],1);