Finding Coordinates of Multiple Contours

5 次查看（过去 30 天）
Brian DeCicco
Brian DeCicco 2021-5-21
评论： Adam Danz 2021-5-28
Hello, I am working with a dataset and trying to find the coordinates (lat/long) of each contour drawn from my loop. When I run the loop, each new iteration replaces the the previous position data of the last contour drawn, and I wish to get the position data for all 40 lines drawn contained within the same matrix. The position data of the line (and each line) should be stored under the variable "C" which is a 2 x ? double matrix. I appreciated any assistance you may be able to provide! Below is my code thus far:
for n = 1:40
[C,h] = contour(long3,lat3,squeeze(yearly_array2(:,:,n))',[287 287],'y');
end
  2 个评论
Adam Danz
Adam Danz 2021-5-21
编辑：Adam Danz 2021-5-21
There are lots of functions on the file exchange that get the contour line coordinates.
Brian DeCicco
Brian DeCicco 2021-5-21
Thank you Adam, your response here was really helpful as well and I will try these in my coding!

请先登录，再进行评论。

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

采纳的回答

KSSV
KSSV 2021-5-21
C = cell(40,1) ;
H = cell(40,1)
for n = 1:40 ;
[c,h] = contour(long3,lat3,squeeze(yearly_array2(:,:,n))',[287 287],'y');
C{n} = c ;
H{n} = h ;
end
celldisp(C)
  2 个评论
Brian DeCicco
Brian DeCicco 2021-5-21
Thank you! I very much appreciate your help. 1 final question based on this...if I wanted to loop through each point associated with each of the 40 contours and plot the max and min value y-values (i.e. lattitude) at each x-value (i.e. longitude) based on your code above, how would you recommend coding that up most efficiently?
Adam Danz
Adam Danz 2021-5-28
@Brian DeCicco There are some files on the file exchange that organize the contour matrix output into a table, a structure, or a cell array. Here's a list:
https://www.mathworks.com/matlabcentral/answers/836078-finding-coordinates-of-multiple-contours#comment_1534428

请先登录，再进行评论。

更多回答（1 个）

Steven Lord
Steven Lord 2021-5-21
Ran in:
Why create the contours one at a time? Why not pass a vector of contour levels to the contour function and either parse the ContourMatrix property of the contour object or the contour matrix returned as the first output of contour? This will also eliminate the need to break the ambiguity of whether a scalar value represents a number of contours to plot or the height at which to plot one single contour line.
[X, Y, Z] = peaks;
surf(X, Y, Z);
hold on
[C, h] = contour3(X, Y, Z, -5:5, 'r'); % I could also have used surfc
shading interp
C(:, 1)
ans = 2×1
-5 25
This tells me the contour at z = -5 has 25 points.
C(:, 1+25+1) % 1 (first element) + 25 (points on first contour) + 1 (move to start of next contour)
ans = 2×1
-4 35
The contour at z = -4 has 35 points, etc.

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by