如何通过在 R 或 MATLAB 中一次添加每个点来动画 3D 散点图答案

【问题标题】：How to animate 3D scatter plot by adding each point at a time in R or MATLAB如何通过在 R 或 MATLAB 中一次添加每个点来动画 3D 散点图
【发布时间】：2019-04-24 23:04:44
【问题描述】：

我有一组 3D 坐标 here。数据有 52170 行和 4 列。每行代表一个点。第一列是点索引号，从 1 增加到 52170。第二到第四列分别是 x、y 和 z 轴的坐标。前10行如下：

seq    x               y        z
1   7.126616    -102.927567 19.692112
2   -10.546907  -143.824966 50.77417
3   7.189214    -107.792068 18.758278
4   7.148852    -101.784027 19.905006
5   -14.65788   -146.294952 49.899158
6   -37.315742  -116.941185 12.316169
7   8.023512    -103.477882 19.081482
8   -14.641933  -145.100098 50.182739
9   -14.571636  -141.386322 50.547684
10  -15.691803  -145.66481  49.946281

我想创建一个 3D 散点图，其中每个点都使用 R 或 MATLAB 按顺序添加到该图中。先添加第一行表示的点，再添加第二行表示的点，...，一直到最后一个点。

另外，我希望控制添加点的速度。

对于二维散点图，我可以使用以下代码：

 library(gganimate)
 x <- rnorm(50, 5, 1)
 y <- 7*x +rnorm(50, 4, 4)
 ind <- 1:50
 data <- data.frame(x, y, ind)

ggplot(data, aes(x, y)) + geom_point(aes(group = seq_along(x))) + transition_reveal(ind)

但我找不到有关如何为 3D 散点图执行此操作的信息。谁能告诉我如何做到这一点？谢谢。

【问题讨论】：

您应该简化您的数据并将其与您尝试过的内容一起发布在这里。
这个answer 很好地描述了这个过程。只需适应您的需求...
你试过these solutions吗？
感谢您的建议。我搜索了使用 MATLAB 执行此操作的方法，但似乎没有直接的答案。我是 MATLAB 新手，所以我很难适应代码。我对 R 比较熟悉。欢迎提出任何使用 R 的建议。
如果您不想要 MATLAB 解决方案，请不要标记为 MATLAB。鉴于您已标记为 MATLAB，因此有数百个示例脚本可以随心所欲地执行，所有这些都可以通过 Google 访问。

标签： r matlab animation ggplot2 gganimate

【解决方案1】：

这是MATLAB的答案

在一般情况下，可以按照相同的方法对绘图（或 3d 图、散点图、曲面或其他图形对象）进行动画处理：

执行第一个 plot/plot3/scatter/surf，并检索其句柄。第一个图可以包含第一个“初始”点集，甚至可以为空（使用NaN 值创建一个带有不可见数据点的图）。
设置轴限制和所有其他将被固定的可视化选项（视点、摄像机角度、闪电......）。无需设置动画期间将要进化的选项。
在循环中，更新最小的绘图对象属性集：XData、YData（ZData 如果是 3D 绘图，CData 如果绘图对象有一些并且您想要动画颜色）。

下面的代码是适用于您的情况的上述方法的实现：

%% Read data and place coordinates in named variables
csvfile = '3D scatter plot.csv' ;
data = csvread(csvfile,2) ;
% [optional], just to simplify notations further down
x = data(:,2) ;
y = data(:,3) ;
z = data(:,4) ;

%% Generate empty [plot3] objects
figure
% create an "axes" object, and retrieve the handle "hax"
hax = axes ;
% create 2 empty 3D point plots:
% [hp_new]   will contains only one point (the new point added to the graph)
% [hp_trail] will contains all the points displayed so far
hp_trail  = plot3(NaN,NaN,NaN,'.b','Parent',hax,'MarkerSize',2) ;
hold on
hp_new    = plot3(NaN,NaN,NaN,'or','Parent',hax,'MarkerSize',6,'MarkerEdgeColor','r','MarkerFaceColor','g','LineWidth',2) ;
hold off

%% Set axes limits (to limit "wobbling" during animation)
xl = [min(x) max(x)] ;
yl = [min(y) max(y)] ;
zl = [min(z) max(z)] ;
set(hax, 'XLim',xl,'YLim',yl,'ZLim',zl)

view(145,72)    % set a view perspective (optional)

%% Animate
np = size(data,1) ;

for ip=1:np
    % update the "new point" graphic object
    set( hp_new , 'XData',x(ip), 'YData',y(ip), 'ZData',z(ip) )

    % update the "point history" graphic object
    % we will display points from index 1 up to the current index ip
    % (minus one) because the current index point is already displayed in
    % the other plot object
    indices2display = 1:ip-1 ;
    set(hp_trail ,...
        'XData',x(indices2display), ...
        'YData',y(indices2display), ...
        'ZData',z(indices2display) )

    % force graphic refresh
    drawnow

    % Set the "speed"
    % actually the max speed is given by your harware, so we'll just set a
    % short pause in case you want to slow it down
    pause(0.01) % <= comment this line if you want max speed
end

这将产生：

【讨论】：