MATLAB Answers

How to plot volumetric concentration points in a slice meshgrid?

6 views (last 30 days)
Timofey Broslav
Timofey Broslav on 4 May 2018
Edited: Timofey Broslav on 4 May 2018
I have a 2-D array A(time, concentration value at a z value). I would like to use this for my mesh grid for each time step. i.e for grid 1 it will be (1, concentration values in z) then grid two would have (2, concentration values in z). I made the assumption that the concentration at each z point will be equal for the y and x plane. So far I have the code shown below, however, the plot i get is a blank white plot. I think i need to make v dependent on xp,yp, and zp, but I do not know how to do that. Anyone have any ideas how to implement this?
Thank you
dx = 6.35/120;
dy = 5.08/120;
dz = (0.01+0.0111)/120;
xp = 0:dx:6.35;
yp = 0:dy:5.08;
zp= 0 :dz:(0.01+0.0111);
[xp,yp,zp] = meshgrid(xp,yp,zp);
for j = 1:10000
v = B(:,j);
xslice = [3.175,6.35]; % location of y-z planes
yslice = [2,5.089]; % location of x-z plane
zslice = [0,.0111]; % location of x-y planes
Where B is a (120 x 80000) double. With 120 z time steps and 80000 time steps.
Timofey Broslav
Timofey Broslav on 4 May 2018
Would you know how I would be able to create a more smooth image while keeping all concentration points (120x120x120)? This just looks ugly. I am also looking into displaying a color scale, if you have any hints on that.

Sign in to comment.

Accepted Answer

Wick on 4 May 2018
Ok, looks like you've got your volumetric plots going. It's easy from here.
First, to display a color scale, the command is
Second, to explain how to get rid of the black lines you've got to understand something called handles. MATLAB represents its graphics in a series of nested structures called graphics handles. When you create a plot object you can usually assign the output to a variable. That will assign that variable the handle of that object. You can then interrogate the properties of that object by using the 'get' command. You can 'set' properties to change them. I recommend looking at the halp for more detail. However, for your case, it's just going to be:
h = slice(X,Y,Z,V,Sx,Sy,Sz);
  1 Comment
Timofey Broslav
Timofey Broslav on 4 May 2018
Thank you Chad, here is what I got. (btw I had to flip around my z limits because I realized it was upside down haha)

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!