3d plotting of Lat, Long and depth

Question

P_L 2019-2-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/443982-3d-plotting-of-lat-long-and-depth

编辑： P_L 2019-2-28

Hi,

I am trying to create a 3D plot - please see attached image.

However, i need to flip the scale so that the depth (km) increases downwards and I need to plot it in a 1:1 aspect ratio as all 3 axis are not equal in reality... Does anyone know how to do this?

Eventually i will need to plot elipsoid errors around these data points - if anyone has any experience with this, then I would be extremely grateful if you could share how you did it!

Many thanks!!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Cris LaPierre 2019-2-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/443982-3d-plotting-of-lat-long-and-depth#answer_360221

编辑：Cris LaPierre 2019-2-8

Try this:

axis equal
ax = gca;
ax.ZDir = 'reverse';

6 个评论
显示 4更早的评论隐藏 4更早的评论

Cris LaPierre 2019-2-9

编辑：Cris LaPierre 2019-2-14

Let's be clear: axis equal is working just fine, it's just not creating the result you want. As the name implies, it makes the step size the same on every axis, resulting in a visually scaled plot. It knows nothing about units.

Since you know the relationship between the two units (1deg = 111km), you can plot with everything being the same units (I'd convert km to degree). I would then update the Z tick labels to be in degrees (current value*111).

Just note, manually updating your labels means they won't automatically update if additional data is added or if the figure is resized. Make sure you have your final image before changing them.

%% Data
Lat=[40.758853	40.763073	40.761671	40.760265	40.763082	40.758853	40.758853	40.756036	40.761671	40.761677	40.763085	40.763082	40.760265	40.758853	40.753203];
Long=[30.396377	30.39452	30.400074	30.401927	30.405623	30.396377	30.396377	30.392681	30.400074	30.407476	30.409325	30.405623	30.401927	30.396377	30.374183];
Depth=[10.623437	11.313867	11.136328	10.248633	10.84043	10.623437	12.201562	10.347266	11.373047	10.426172	10.288086	10.011914	10.011914	9.992187	11.807031];
%% Script
scale = 111; % km/deg
scatter3(Long(:), Lat(:), Depth(:)/scale,'MarkerEdgeColor','k', 'MarkerFaceColor',[0 .75 .75]);
axis equal
ax = gca;
ax.ZDir = 'reverse';
zVal = str2double(ax.ZTickLabel)*scale;
ax.ZTickLabel = num2str(zVal);
xlabel('Long (deg)') 
ylabel('Lat (deg)') 
zlabel('Depth (km)')
%grid on