我有一个纬度/经度向量,另一个包含 Z 值,我需要将其显示在 3D 地图上,如下图右侧所示。我尝试了 bar3,但它很麻烦,因为它需要创建多个图表。
这是一些代码:
S4 = shaperead(filename)
plot([S4.X],[S4.Y],'k'); % plots the map from a shapefile I loaded previously
XX = [-50 -51 ...];
YY = [-1 -2 ...];
ZZ = [ 2.2 3.2 ... ];
stem3(XX,YY,ZZ) % this is an option, but doesn't look good!! :(
关于如何做到这一点有什么想法吗?
谢谢!
所以基于搬迁的能力bar3
无论您想要什么,这里都有一个简单的代码可以完成这项工作。
您不感兴趣但可能感兴趣的事情:
1-填充地图。 2.-设置灯光 3.-创建标签
2、3,都可以轻松完成。 1 我不知道。
Result:
CODE:
clear;clc;
S = shaperead('usastatehi.shp')
S(2)=[]; % Delete Alaska
S(10)=[];% Delete Haway
hold on
plot([S.X],[S.Y],'k'); % plots the map from a shapefile I loaded previously
% syntetic data
Y=[40,45,25];
X=[-100,-85,-80];
Z=[0.5 2.3 1.4];
cmap=colormap(strcat('parula(',num2str(length(Z)),')')); % Create a colormap. Change parula for anything you preffer
ar=abs(diff(ylim))/abs(diff(xlim));
for ii=1:length(X)
h=bar3(Z(ii));
Xaux=get(h,'Xdata');
Xaux=Xaux-1+X(ii);
set(h,'Xdata',Xaux)
Yaux=get(h,'Ydata');
Yaux=Yaux-1+Y(ii);
set(h,'Ydata',Yaux)
set(h,'FaceColor',cmap(ii,:));
end
axis off
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)