how to efficiently find the indices?
显示 更早的评论
Hi,
I have a matrix of N-by-M with integers. I need to efficiently find the indices for all of the unique elements in the matrix. The solution I have is via a "for" loop:
uM = unique (M(:));
for i = 1 : length(uM)
I(i) = find(M == uM(i));
end
This works fine, but with a large matrix, this is slow. I wonder if there are better solutions. thanks very much!
1 个评论
Geoff
2012-5-3
One way to make this code faster without changing anything fundamental would be to preallocate the cell array before your loop:
I = cell(length(uM),1);
采纳的回答
Richard Brown
2012-5-2
sort is your friend:
First we sort the entries of A so that like entries are adjacent
[S, iSorted] = sort(A(:));
Then we find the markers - where the value changes, plus the endpoints
markers = [0; find(diff(S)); numel(A)];
All we need to do is find the entries of iSorted corresponding to each block.
for i = 1:numel(markers)-1
uM(i) = S(markers(i)+1);
I{i} = iSorted(markers(i)+1:markers(i+1));
end
2 个评论
Pinpress
2012-5-3
Richard Brown
2012-5-3
for loops get a bit of an unnecessarily bad rap sometimes
更多回答（6 个）
Leah
2012-5-2
Matlab has a nice function built in for this
[B,I,J] = UNIQUE(...) also returns index vectors I and J such
that B = A(I) and A = B(J) (or B = A(I,:) and A = B(J,:)).
Andrei Bobrov
2012-5-3
[uM,n,n] = unique (M(:));
I = accumarray(n,1:numel(n),[],@(x){sort(x)});
2 个评论
Richard Brown
2012-5-3
I think even @(x) {x} would meet the brief for your function
Andrei Bobrov
2012-5-3
Hi Richard! I agree with you.
Richard Brown
2012-5-2
Depending on whether you want the first or the last occurence
[uM, I] = unique(M, 'first')
[uM, I] = unique(M)
Walter Roberson
2012-5-2
0 个投票
Your suggested code will not work if there are any duplicates, as the find() would return multiple values in that case and multiple values cannot be stored into a single numeric array element.
Have you considered using the second or third return value from unique() ?
Pinpress
2012-5-2
0 个投票
2 个评论
Oleg Komarov
2012-5-2
A and J will give you what you want.
Richard Brown
2012-5-2
Not really - J is just A(:), but with the unique elements replaced with 1:nUnique. So it's no better
Geoff
2012-5-2
This comes straight out of some of my own code... I guess it's the reverse of what you want though.
uM = unique(M);
I = arrayfun(@(x) find(uM==x,1), M);
For every element in M, it gives an index into uM. I use this to reduce columns of data in a matrix that are common to multiple targets.
So you seem to want: for every element in uM an array of indices into M. The result of course would be a cell array.
This would be:
I = arrayfun(@(x) find(M==x), uM, 'UniformOutput', false);
2 个评论
Richard Brown
2012-5-2
I think that's what he was trying before, but found the repeated calls to find to be too slow
Pinpress
2012-5-3
类别
在 帮助中心 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
