I am using isosurface to plot 3d density maps. I would like to plot multiple maps in the same figure, each one with a different color. How can I do so?

nBins=20;

D=cell(size(A,1),size(A,2));

for i=1:size(A,1)

for j=1:size(A,2)

if ~isempty(A{i,j})

x=A{i,j}(:,1); y=A{i,j}(:,2); z=A{i,j}(:,3);

xBins=linspace(min(x),max(x),nBins);

yBins=linspace(min(y),max(y),nBins);

zBins=linspace(min(z),max(z),nBins);

D{i,j}=zeros(nBins,nBins,nBins);

for ii=1:numel(x)

xi=find((x(ii)>xBins),1,'Last');

yi=find((y(ii)>yBins),1,'Last');

zi=find((z(ii)>zBins),1,'Last');

D{i,j}(xi,yi,zi)=D{i,j}(xi,yi,zi)+1;

end % for ii=1:numel(x)

D{i,j}=smooth3(D{i,j});

%isosurface(D{i,j})

%hold on

end % if ~isempty(A{i,j})

end % for j=1:size(A,2)

end % for i=1:size(A,1)

Answer 1

[Mil Shastri]

You can try this:

[x,y,z,v] = flow;
figure;
p = patch(isosurface(x,y,z,v,-3));
p.FaceColor = 'yellow';
hold on;

p = patch(isosurface(x*3,y*1,z*3,v,-3));
p.FaceColor = 'red';
p.EdgeColor = 'none';

p = patch(isosurface(x*1,y*3,z*3,v,-3));
p.FaceColor = 'green';
p.EdgeColor = 'none';

hold off;

%optional
camlight 
lighting gouraud

See isosurface documentation: https://www.mathworks.com/help/matlab/ref/isosurface.html

[Guy Nir]

To make it clear, here are the changes I made to the script, following your suggestion - thank you!

%isosurface(D{i,j})

p=patch(isosurface(yBins,xBins,zBins,D{i,j},2)); % MATLAB'sstaff suggestion (Mil Shastri)

p.FaceColor = C(j,:); % I added a colormap (not shown)

p.EdgeColor = 'none';

hold on;

Answer 2

[Mil Shastri]

Hi Guy, Please provide a sample for A matrix too. The code doesn't run without it.

Answer 3

[Guy Nir]

I have attached now a sample of the cell array 'A'. I called the sample cell array 'A1', so you would have to rename it as 'A'.

Thank you,

Guy

Answer 4

[Guy Nir]

% Sorry, please use this updated script, otherwise it will plot them all on top of each other.

GridRes=71;

minX=1.2206e+04;

maxX=4.1832e+04;

minY=1.9815e+03;

maxY=39500;

minZ=-913.6920;

maxZ=8.6982e+03;

nBinsX=round((maxX-minX)/GridRes);

nBinsY=round((maxY-minY)/GridRes);

nBinsZ=round((maxZ-minZ)/GridRes);

D=cell(size(A,1),size(A,2));

for i=1:size(A,1)

for j=1:size(A,2)

if ~isempty(A{i,j})

x=A{i,j}(:,1); y=A{i,j}(:,2); z=A{i,j}(:,3);

xBins=linspace(minX,maxX,nBinsX);

yBins=linspace(minY,maxY,nBinsY);

zBins=linspace(minZ,maxZ,nBinsZ);

D{i,j}=zeros(nBinsX,nBinsY,nBinsZ);

for ii=1:numel(x)

xi=find((x(ii)>xBins),1,'Last');

yi=find((y(ii)>yBins),1,'Last');

zi=find((z(ii)>zBins),1,'Last');

D{i,j}(xi,yi,zi)=D{i,j}(xi,yi,zi)+1;

end % for ii=1:numel(x)

D{i,j}=smooth3(D{i,j});

isosurface(D{i,j})

hold on

end % if ~isempty(A{i,j})

end % for j=1:size(A,2)

end % for i=1:size(A,1)

axis equal