Plotting 3D surface contours of a variable with respect to x, y and z

Question

shravani banerjee 2023-6-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1981619-plotting-3d-surface-contours-of-a-variable-with-respect-to-x-y-and-z

编辑： shravani banerjee 2023-6-22

MIXING_RATIOe.zip

I am trying to plot a 3D surface contour plot a varaible (dust mixing ratio) with respect to lat (x-axis), long (y-axis) and levels (z axis) (attached csv file). But I am unable to plot that. While making the meshgrid I am getting an error Requested 248040x248040x2 (916.8GB) array exceeds maximum array size preference. How to plot this?

2 个评论
显示无隐藏无

Rik 2023-6-12

What code did you try? You can put it right in your question and run the code to show us exactly how to reproduce the error.

And what do you want exactly? A surface plot or a contour plot? Both rely on a single value for the height given an x-y coordinate pair. Your error suggests you have a list of coordinates that should be reshaped instead of using meshgrid.

shravani banerjee 2023-6-12

在 MATLAB Online 中打开

I want a 3D contour plot with x as lat, y as lon, and z as levels and the legend colour will represent my variable. I tried to created a mesh to make a 3D contour plot.

T=readtable('MIXING_RATIOe.csv');
lat = T.lat;
lev = T.levels;
lon = T.lon;
values = T.dust_mixing_ratio;
y = lat;
m = lon;
z = lev;
[Y,M,Z] = meshgrid(y,m,z);

请先登录，再进行评论。

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

Answer 1

Rik 2023-6-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1981619-plotting-3d-surface-contours-of-a-variable-with-respect-to-x-y-and-z#answer_1254134

在 MATLAB Online 中打开

MIXING_RATIOe.zip

It sounds like you want a 3D pointcloud where the color represents a 4th variable.

As you can see below, this might not actually be a very informative plot.

t = tic;%keep track of time

unzip('MIXING_RATIOe.zip')

T=readtable('MIXING_RATIOe.csv');

lat = T.lat;

lev = T.levels;

lon = T.lon;

values = T.dust_mixing_ratio;

cmap = colormap('parula');

N_colors = size(cmap,1);

unique_values = unique(values);

C_ax = [min(values) , max(values)];

N = numel(unique_values);

%for n=1:N

for n=randperm(N) % this is to make sure stopping halfway doesn't distort the plotted colors

L = values == unique_values(n);

col_idx = interp1(...

linspace(C_ax(1),C_ax(2),N_colors),... % data range

1:N_colors,... % colormap indices available

unique_values(n)); % query current value

plot3(lat(L),lon(L),lev(L),...

'.','color',ind2rgb_private(col_idx,cmap))

if n==1,hold('on'),elseif n==N,hold('off'),end

% for the online environment, exit after 20 seconds of work

if toc(t)>20,break,end,hold('on')

end

colorbar,caxis(C_ax) % add colorbar and adjust limits to match the plotted results

function col = ind2rgb_private(col_idx,cmap)

% extend ind2rgb to allow decimal indices

N_colors = size(cmap,1);

col = [...

interp1(1:N_colors,cmap(:,1),col_idx),...

interp1(1:N_colors,cmap(:,2),col_idx),...

interp1(1:N_colors,cmap(:,3),col_idx)];

end

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

shravani banerjee 2023-6-12

Here is an example of the plot I wanted to make.

Rik 2023-6-12

I suspect you could get something like this if you set the Alpha property of a patch object below a certain threshold (and scale that with the value as well).

I don't know how fine your slices need to be before you don't notice them anymore. I suggest you give it a try. I'm not aware of any way to create this using solid voxels in Matlab.

请先登录，再进行评论。