Finding the highest mountain peak

6 次查看（过去 30 天）
Brantosaurus
Brantosaurus 2022-9-29
评论： Brantosaurus 2022-9-30
How might i find the highest mountain peak, and its name, from a defined set of named mountains?
I have data for a set of 10 mountains as 4 matrices.
[X] and [Y] are single row and column matrices defining distances along orthogonal horizontal coordinates.
[Z] is a 3D matrix (:,:,i) containing corresponding x,y height values and an index i=1:10 to define the mountain set.
[Name] is a 10 element 1D text matrix of of the mountain names.
Using a for loop I can plot the surfaces, find the peaks of each mountain and their names:
surf(X,Y,Z(:,:,i))
maximum=max(max(Z(:,:,i)))
Name(i)
How might i find the overall highest peak, and its name in one hit (without the loop)?

请先登录，再回答此问题。

采纳的回答

Image Analyst
Image Analyst 2022-9-29
Why can't you just do this to find the row and column of the max for each of the 10 slices of data:
for k = 1 : size(Z, 3)
% Extract this mountain set.
thisMountainSet = Z(:, :, k);
% Find the max value of it
maxValue(k) = max(thisMountainSet(:));
% Find row(s) and column(s) where it occurs.
[rowsOfMax, colsOfMax] = find(thisMountainSet == maxValue(k));
% Now do something with this information....
end
  3 个评论
显示 1更早的评论
Image Analyst
Image Analyst 2022-9-30
Personally I prefer using find() to get the rows and columns directly and immediately rather than using the confusing, intermediate function ind2sub(). Note also about the solution you accepted, max only gives you one max location (if you ask for the index), not all of them like find() does, so it your max occurs in more than one location, you may want to know that, and max won't tell you that.
Brantosaurus
Brantosaurus 2022-9-30
Had a play with your script. I now get your point! Working in the loop you suggest on a step-by-step basis looks the safest option to identify potential multi peaks and their locations.
I managed to produce a summary after the loop using:
[M] = max(maxValue(:));
[I] = find(maxValue == M);
It appears to find the (overall) max and case(s) where they occur to refer back as necessary.
Many thanks for the advice. Much appreciated.

请先登录，再进行评论。

更多回答（2 个）

KSSV
KSSV 2022-9-29
编辑：KSSV 2022-9-29
[val,idx] = max(Z,[],3) ;
Get the maximum along third dimension, idx gives you the thirs index where max occurs. Now find the max of val, get it's index.
Or you ccan also use.
A = rand(3,3,3);
[M,I] = max(A(:));
[I_row, I_col, I_page] = ind2sub(size(A),I)
  1 个评论
Brantosaurus
Brantosaurus 2022-9-30
It was the linear (column) index idx that was causing my confusion. Your second example using ind2sub sorted it out. Just what i needed. Much appreciated. :)

请先登录，再进行评论。

Brantosaurus
Brantosaurus 2022-9-29
Thanks. The below worked to get the overall highest peak:
[val,idx]=max(Z,[],3);
[value,index]=max(val);
[gvalue,gindex]=max(value); #gvalue is correct
Struggling with the name index! Tried to find correct 'position' with:
[row,col]=find(Z==gvalue)
but it doesn't work. Included is val and idx matrix if it helps.

社区

更多回答在  ThingSpeak Community

类别

在 Help Center 和 File Exchange 中查找有关 Pulsed Waveforms 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by