contourc - plotting contour matrix

2012 6 1
2 个回答
回答已采纳
更新时间：2021 9 2
9 次查看（30 天）

1 个投票

I would like to find the contours, manipulate them, and then plot them.
This plots contours
x = 0:0.1:1; y = 0:0.1:1;
[X,Y] = meshgrid(x,y);
Z = sin(X).*cos(Y);
[con_mat, h] = contour(x, y, Z);
To get the contour information without plotting, can use contourc.
con_mat = contourc(x, y, Z);
However, there appears to be no built-in way to plot this "contour matrix".
Any thoughts?
John

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

 采纳的回答

Darin
Darin 2012-6-5

1 个投票

Try this: it uses cplot, which you can find on Matlab central.
If you find something more elegant, please let us all know... it MUST be buried in the Matlab library function SOMEWHERE.
function[] = plot_contourmatrix(C);
% plot_countourmatrix - Plots a contour matrix c as returned from contour
%
% plot_contourmatrix(C)
% Rev History:
% 06-04-12 Created
figure(gcf);
holdstate = ishold;
hold on;
i = 1;
while i<=length(C)
lev = C(1,i);
cnt = C(2,i);
cplot(C(1,i+(1:cnt)),C(2,i+(1:cnt)),repmat(lev,1,cnt));
i = i+cnt+1;
end;
if ~holdstate
hold off;
end;
return
test code - execute as cell
figure(1);
[x,y] = meshgrid(-3:.1:3);
z = sinc(hypot(x,y));
[C,h] = contour(x,y,z);
axis equal
figure(2)
plot_contourmatrix(C);
axis equal

2 个评论
显示 无 隐藏 无

John Mahoney
John Mahoney 2014-12-4
Thanks Darin. I was using something similar, but was checking for integer values in C(2,:). This fails (almost never) when the y value is an integer. Yours does not have this problem.
It would be ideal to have a Matlab function not only for plotting this matrix, but for extracting the points in these contours (would use the same code as above.) as folks often want to process these contours before plotting etc.
Johannes Korsawe
Johannes Korsawe 2015-1-26
Where to find cplot? it is not longer on the FEx ?

请先登录，再进行评论。

更多回答（1 个）

Liviu Ivanescu
Liviu Ivanescu 2018-8-18
编辑：Liviu Ivanescu 2018-8-18

1 个投票

Here is a way to plot contourc data containing several contours of the same value.
cnt = contourc(matrix,[-2 -2]);
szc = size(cnt);
idz = 1;
while idz<szc(2)
izi = cnt(2,idz);
cnt(2,idz) = nan;
idz = idz+izi+1;
end
plot(cnt(1,:),cnt(2,:))

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Nikola Stan
Nikola Stan 2021-1-26
编辑：Nikola Stan 2021-1-26
your example was useful to me. I added some code that allows it to keep different isolines separate, which is what I needed:
cnt = contourc(matrix,[-2 -2]);
szc = size(cntr);
idz = 1;
contourNo = 1;
while idz<szc(2)
izi = cntr(2,idz);
cntr(2,idz) = nan;
contourXY{contourNo} = cntr(:,idz+1:idz+izi);
idz = idz+izi+1;
contourNo = contourNo+1;
end
figure()
hold on
for k = 1:contourNo-1
plot(contourXY{k}(1,:), contourXY{k}(2,:));
end
in the first line did you mean "cntr = contourc(matrix,[-2 -2]);"?
i think you forget the "r" in cntr
thanks

请先登录，再进行评论。

类别

在 帮助中心 和 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