Borrar filtros
Borrar filtros

Plotting portion of a color coded scatter plot using a for loop

3 visualizaciones (últimos 30 días)
AN
AN el 11 de Abr. de 2023
Comentada: AN el 11 de Abr. de 2023
Hi all,
I have created a value-based color coded 4-dimensional scatter plot of a heart (color coded based on activation times corresponding to each xyz point). The xyz coordinates are designated by vectors (x_points, etc...) and the corresponding colors to each xyz point are defined in an nx3 array called rgb_time. I have divided the scatter plot into cubic voxels using discretize, and am attempting to create scatter plots of individual voxels displaying the same color code of the original scatterplot. I am using the following code to loop through the points and determine the ones within the voxel of interest before plotting them:
i_max = 6; %arbitrary value
%graph only the cube containing the most colors (area of interest)
figure()
hold on
for i = 1:size(rgb_time)
if voxelIDs(i) == i_max
scatter3(x_points(i),y_points(i),z_points(i),20,rgb_time(i),'filled')
xlabel('x');ylabel('y');zlabel('z');title('voxel containing isochronal crowding');
str = 'early - red - orange - yellow - green - cyan - blue - purple - pink - late';
end
end
For some reason, however, the resultant voxel plot has distorted colors that are not the same as those displayed on the original scatter plot. Could anyone have an idea as to why this is? I have tried moving the hold on statement to various locations to no avail. Thank you!
  1 comentario
AN
AN el 11 de Abr. de 2023
I have attached examples of the original scatter plot with the desired colored regions, as well as a subset/voxel scatter plot with distorted colors.

Iniciar sesión para comentar.

Respuesta aceptada

Ran Yang
Ran Yang el 11 de Abr. de 2023
Editada: Ran Yang el 11 de Abr. de 2023
rgb_time is a Nx3 array, so you need to call all 3 columns using rgb_time(i, :) to get the color for each voxel. (Right now you're only getting column 1).

Más respuestas (0)

Categorías

Más información sobre Scatter Plots en Help Center y File Exchange.

Productos


Versión

R2017a

Community Treasure Hunt

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

Start Hunting!

Translated by