Computing findpeaks along 3-Dimensional matrix
显示 更早的评论
I have a 3-dimensional matrix, after some code manipulation on a cell array. The matrix is of size 101x151x476.
The first number represents the y space, the second number is the x space, and the third number is the number of time steps. So essentially, the matrix contains magnitudes along a y-x plane, for varying timesteps.
I am trying to compute the peaks at each (y,x) across the timesteps. I am doing this using findpeaks. This is what I have so far:
for yIndex= 1:1:101
for xIndex= 1:1:151
[peaks{y,x}] = findpeaks(a3Dmatrix(yIndex,xIndex,:));
end
end
However, the cell "peaks" is incorrect, as it is not of the size 101x151 (y by x). Any help with this?
采纳的回答
Star Strider
2020-10-20
I am guessing here, since I do not have your data. Since findpeaks is going to return vectors for the peaks, locations (and othher outputs if you want them), I would do something like this:
for yIndex= 1:1:101
for xIndex= 1:1:151
[peaks{yIndex,xIndex,:},locs{yIndex,xIndex,:}] = findpeaks(a3Dmatrix(yIndex,xIndex,:));
end
end
I cannot test that because I do not have your data, so I am posting it as UNTESTED CODE.
6 个评论
Vivek
2020-10-20
Star Strider
2020-10-20
I do not understand what line threw that error. Cell arrays can certainly contain empty vectors (denoted as {[]}), and it is certainly possible that there are no identifiable peaks in your signal (constant, or that the maximum is at one or the other end of the vector, so not a ‘peak’ as findpeaks defines it). That would result in findpeaks returning an empty output for that particular set of coordinates.
Perhaps a better option is:
for yIndex= 1:1:101
for xIndex= 1:1:151
[pks,lcs] = findpeaks(a3Dmatrix(yIndex,xIndex,:));
peaks{yIndex,xIndex} = pks;
locs{yIndex,xIndex} = lcs;
end
end
That might make it easier to understand what the problem is, since you could look at the size function output of size(pks).
If that indicates an empty vector, then:
[max,idx] = maxk(a3Dmatrix(yIndex,xIndex,:), 2);
could suggest what the problem is.
You could put that in an if block to trap the empty vector, for example:
for yIndex= 1:1:101
for xIndex= 1:1:151
[pks,lcs] = findpeaks(a3Dmatrix(yIndex,xIndex,:));
if (numel(pks) < 1))
[pks,lcs] = maxk(a3Dmatrix(yIndex,xIndex,:), 2);
end
peaks{yIndex,xIndex} = pks;
locs{yIndex,xIndex} = lcs;
end
end
Then later look at ‘locs’. If they are the first or last 2 points in the vector (or the first and last points), there are no peaks in that segment.
Vivek
2020-10-20
Star Strider
2020-10-20
The problem is that I did not test my code.
Try this version:
for yIndex= 1:size(a3Dmatrix,1)
for xIndex= 1:size(a3Dmatrix,2)
vect = squeeze(a3Dmatrix(yIndex,xIndex,:));
[pks,lcs] = findpeaks(vect);
if (numel(pks) < 1)
[pks,lcs] = maxk(vect, 2);
end
peaks{yIndex,xIndex} = pks;
locs{yIndex,xIndex} = lcs;
end
end
It ran without error with a random 3D matrix, and appeared to produce reasonable results.
Vivek
2020-10-20
Star Strider
2020-10-20
My pleasure!
If my Answer helped you solve your problem, please Accept it!
.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
