plotting lat, lon, altitude, variable

2015 3 9
4 个回答
回答已采纳
更新时间：2015 3 16
13 次查看（30 天）

0 个投票

I have 3 vectors of latitude, longitude, and elevation respectively. I also have a corresponding 3d matrix of size latitude*longitude*elevation which contains some data (Temperature). I want to display the Temperature data in 3d space. What is the best way to do this? Note: I don't want to use trisurf because it requires reshaping the matrix into 3 separate vectors, with which I found difficulty in tracking the data points.
Many thanks.
-Sagar

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

 采纳的回答

Chad Greene
Chad Greene 2015-3-11

0 个投票

You can try pcolor3, although it does not work for Matlab 2014b.
% Some data:
[lon,lat,z] = meshgrid(-180:3:180,-90:2:90,0:50:1000);
T = 20 - .1*abs(lat) - .01*z + sind(lon);
% Plot dataset:
pcolor3(lon,lat,z,T);
% Add labels:
xlabel('longitude')
ylabel('latitude')
zlabel('elevation (m)')
title('temperature or something')
axis tight

更多回答（3 个）

Sagar
Sagar 2015-3-12
编辑：Sagar 2015-3-12

0 个投票

Hi Chad, thanks a lot. I tried pcolor3 that you suggested, but with no success.
Please let me know if I did something wrong. I did:
pcolor3(X, Y, Z, permuted_ext_mean); % X, Y, Z and permuted_ext_mean all have 72*85*208 size.
I got following error in Matlab 2013b. Also, it would be great if you could put a colorbar in the output image. Error using interp3/stripnanwrapper (line 204) Insufficient finite values to interpolate.
The closest one I found for my purpose is isosurface, although I am not completely satisfied: surf2 = isosurface(X, Y, Z, permuted_ext_mean, isovalue);

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

Chad Greene
Chad Greene 2015-3-12
编辑：Chad Greene 2015-3-12
Interesting. I'm guessing your permuted_ext_mean variable has a lot of NaN values? Around line 74 of pcolor3 there's a slice line with the 'cubic' interpolation option chosen. There's a small chance that changing 'cubic' to 'linear' or 'nearest' will make it work. In the meantime I may tinker with an interpolationless solution.
Chad Greene
Chad Greene 2015-3-12
I've just updated the function with a 'direct' option, which will plot your data directly instead of interpolating between points.
Sagar
Sagar 2015-3-13
Yes, a lot of NaN values. Just tried the updated one, works well now. It is a bit slow though compared to isosurface, may be it is normal for pcolor.
-Sagar
Sagar
Sagar 2015-3-14
编辑：Sagar 2015-3-14
Here is what I got with pcolor3, I am fully satisfied now. Thanks Chad.

请先登录，再进行评论。

Chad Greene
Chad Greene 2015-3-14

0 个投票

Great! You may also want to add a coast line for context as shown here.

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

Sagar
Sagar 2015-3-15
编辑：Sagar 2015-3-15
Hi Chad, that is a great suggestion and I am trying to do that as below. Is there a way to specify the range of latitude and longitude, say, in my case, 0 to 40N, -20 to 70E (so the length of the two vectors would be unequal). c = load('coast.mat'); plot(c.long,c.lat) view (3)

请先登录，再进行评论。

Sagar
Sagar 2015-3-15

0 个投票

Never mind! I was able to plot as I wanted with a small tweak. This is what I did:
load ('coast.mat');
long (long < -20 | long > 70 = NaN;
lat (lat < 0.00 | lat > 40.00) = NaN;
plot (lat, long);

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

Chad Greene
Chad Greene 2015-3-16
Or similarly,
in = ingeoquad(lat,long,[0 40],[-20 70]);
long(~in) = NaN;
lat(~in) = NaN;

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Climate Science and Analysis 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by