how to use output indexes from the MAX function
6 views (last 30 days)
I want to find the indexes of the maximum and the second maximum in 3D array for a certain dimension. My strategy is to use the MAX function for finding the indexes of the first maximum. Use these indexes to make the value of the 3D array –Inf and call again the MAX function in order to find the indexes of the second maximum of the 3D array
My question is how to make optimal use of the returned indexes from the MAX function (without for loops).
cube1 = 10*rand([5 6 7]); %create 3D array
cube2 = cube1; %make copy for com parison
[dummy i1st] = max(cube1,,3); %find first max across 3th dimension
cube1(i1st)=-inf; % this wont work
% for loop does work but is not elegant, require recoding when dimension in changed
cube1(r,c,i1st(r,c)) = -inf;
[dummy i2nd] = max(cube1,,3); %find second max across 3th dimension
Does anyone has some ideas how to get rid of the for loop?
Sven on 11 Dec 2011
Hi Patrick, I think this is what you want to do:
arrSz = [5 6 7]; % Size of matrix
cube1 = 10*rand(arrSz); % Create the matrix
[vals1, idxs1] = max(cube1,,3); % find first max across 3th dimension
% Make a 3D mask corresponding to indices in idxs1
idxsMask = bsxfun(@eq,reshape(1:arrSz(3), 1,1,), idxs1);
cube1(idxsMask) = -inf; % Set first max vals to -inf
[vals2, idxs2] = max(cube1,,3); % find second max across 3th dimension
The trick here is the bsxfun() call. I make a set of indices along the 3rd dimension, then bsxfun does its magic asking if the 2D matrix idxs1 matches that array. bsxfun() is a strange function if you haven't seen it before, but once you understand what it does it can be very useful. And very efficient - much faster than the loop.
Walter's question is a valid one... what if the max value isn't unique along the 3rd dimension? Should your second call to max() return the next index that had the same value as the first call? Or should the first call turn all elements that had the max value to -inf?
If you want to do this second way, simply replace the idxsMask line above with the following:
idxsMask = bsxfun(@eq,vals1,cube1);