0

我有一张从北距 (X)、东距 (y) 和高程 (z) 数据创建的表面图。我想围绕特定半径的特定点（比如 400 m）画一个圆。圆圈必须显示在表面图上。下面是我创建表面贴图的代码。任何帮助将不胜感激。

% load the map data
x = elevgrid(:,4); % northing
y = elevgrid(:,5); % easting
z = elevgrid(:,3); % elevation
% Put data onto a grid
[qx,qy] = meshgrid(linspace(min(x),max(x)),linspace(min(y),max(y)));
F = TriScatteredInterp(x,y,z);
qz = F(qx,qy);
% plot in matrix form
surf(qx,qy,qz)

编辑 这里是示例数据文本文件的链接 - 格式为经度、纬度、海拔。

4

1 回答 1

0

你应该可以很容易地用 plot3 做到这一点。有关您发布的代码的修改版本，请参见下面的代码。

我不确定您的示例数据是如何表示的（为避免将来出现问题，请尝试将您的问题简化为基本要素，包括输入）。

clear;clc
% create data
xy = -2.5 + 5*gallery('uniformdata',[200 2],0);
x = xy(:,1);
y = xy(:,2);
z = x.*exp(-x.^2-y.^2);

% Put data onto a grid
[qx,qy] = meshgrid(linspace(min(x),max(x)),linspace(min(y),max(y)));
F = scatteredInterpolant(x,y,z);
qz = F(qx,qy);
% plot in matrix form
mesh(qx,qy,qz)

% Plot circle
hold on
r = 0.6;
xp = 0.4;
yp = -1.2;
t = linspace(0,2*pi);
Xp = r*sin(t)+xp;
Yp = r*cos(t)+yp;
plot3(Xp,Yp,F(Xp,Yp))
hold off

BR马格努斯

于 2013-10-29T01:52:55.590 回答