← 返回 2026-04-13

一次一步地构想未来:从单张图像做稀疏轨迹自回归扩散预测 Envisioning the Future, One Step at a Time

Stefan Andreas Baumann, Jannik Wiese, Tommaso Martorella, Mahdi M. Kalayeh, Björn Ommer 📅 2026-04-10 👍 12 2026-07-13 08:36
世界模型 扩散模型 稀疏表示 自回归Transformer 规划 轨迹预测 运动预测

用稀疏点轨迹的逐步自回归扩散,从单张图像生成上千种物理一致的开集未来预测,比视频生成快几个数量级。

前置知识

扩散模型与流匹配 (Flow Matching)

扩散模型通过学习一个由噪声逐步去噪到数据的反向过程来生成样本;流匹配则是它的简化版本,直接用一个常微分方程 (ODE) 的速度场 $v_\phi(\Delta x_{t,\tau}, z)$,把样本沿路径参数 $\tau \in [0,1]$ 从高斯先验平滑过渡到目标数据分布。本文用流匹配作为稀疏轨迹增量 $\Delta x_t$ 的后验参数化形式,相比 GMM(高斯混合)头表示能力更强、训练更稳定,同时其 ODE 采样器允许在每一步独立 cache 条件特征。

论文的后验参数化核心建立在流匹配之上,并用 ODE 求解器(如 Euler)逐步去噪生成下一步轨迹增量;如果不熟悉流匹配,无法理解为什么 FM head 能在保持高维分布的同时支持 KV caching 并实现快速采样,也就无法理解 Myriad 相对 GMM/FPT 的核心进步。

自回归 Transformer 与 KV 缓存

Transformer 通过 self-attention 在序列每个位置捕捉长程依赖;自回归模型则按时间顺序逐个生成 token,每步生成依赖之前所有已生成 token 的键/值 (KV) 缓存来避免重算。本文将时间和轨迹两个维度都做因果因式分解:$p(x_{1:T}\mid x_0, I_0)=\prod_{t=1}^{T}\prod_{i=1}^{K} p(x_t^{(i)} \mid x_t^{<i}, x_{<t}, I_0)$,其中轨迹维度的 '时间过去' KV cache 可复用、'当前时间' 维度只需 cross-attention 到图像与新 token,这使得 32 步 rollout 在 H200 上达到约 2200 samples/min 的吞吐量。

这是 Myriad 能快速并行采样上千条轨迹的关键:一旦理解 KV cache 如何复用历史、当前步骤只需增量新 token,就能体会为何论文强调 '推理时利用时间维度' 的 transformer 块比传统 sequential block 快 2×(batch=4)甚至 3.7×(batch=1)。

稀疏点轨迹 (Sparse Point Trajectories) 与 Point Tracking

区别于稠密光流 (dense optical flow) 或视频像素预测,稀疏点轨迹只关注用户指定或自动检测的少量关键点(如物体上的 16 个点)在 2D 像素平面上的未来位置序列。获取方式常用 CoTracker3、TAPNext 等点跟踪器,它们在每帧输出 '查询点-像素' 的匹配关系。本文用 1024 个随机种子点跑 TAPNext 得到伪真值进行训练,推理时仅对 K 个感兴趣的查询点做预测,输出空间是 $\mathbb{R}^{K \times T \times 2}$,而非视频模型的 $\mathbb{R}^{T \times H \times W \times 3}$。

这是论文 '避免视觉税' 的核心:不预测 $512\times 512$ 像素,而只预测 16 个 2D 点;如果不熟悉 point tracking,无法理解 Myriad 输出空间的轻量性,无法理解其训练数据从何而来,也无法看清它在台球规划里 '枚举 2200 条 16 维轨迹' 的计算可行性。

静态相机假设与重投影

假设摄像机不动时,像素坐标变化完全来自物体在场景中的运动,从而把 2D 轨迹学习与相机运动解耦。论文额外训练了一个变体:先用 V-DPM 跟踪得到 3D 轨迹 + 相机位姿,再把轨迹反投影回第一帧视角,从而在动态相机视频上也能学到 '纯场景运动'。这一假设极大简化评测与可视化(图像背景不变、轨迹变化仅由物理交互产生),但把方法限定在静态或可补偿视角的场景。

这是 OWM benchmark 评测设计与训练数据选择的基础;理解它才能解释论文的限制与未来工作——即如何把方法推广到自运动与场景运动联合预测,例如在自动驾驶或机器人第一视角场景中。

研究动机

现有世界模型(无论是视频生成如 Wan2.2 14B、CogVideo-X 1.5 5B、Stable Video Diffusion 1.1 1.5B,还是特征空间模型如 DINO-WM)几乎都 '画' 一个未来:它们要么在像素空间、要么在 latent 空间生成密集张量,再从密集张量里反解运动。这种视觉税 (visual tax) 在开放世界场景里会带来三个具体问题:第一,模型容量被分配到无关外观细节,导致运动建模能力受限;第二,预测长程多模态未来时,单步跳跃到 t=T 的失败使误差累积、模型崩塌;第三,对每个反事实展开都要渲染整段视频,做台球击球这种'枚举上千次假设'的规划任务时成本不可承受。例如作者在 Table 2 中实测 Wan2.2 与 I2V Diffusion 在 5 分钟内只能产生约 20 次动作评估,准确率仅 16%;而 CogVideo-X 在 OWM 评测中 Best-of-5min 全部 DNF(无法完成)。

本文的目标是本文的具体目标是构建一个能从单张图像出发、在静态相机假设下,以千倍速预测 K 个稀疏查询点未来 2D 轨迹分布的自回归扩散模型——即 Myriad;它要求训练数据来自开放世界视频(而非受控仿真),并在测试时既能被动输出多模态未来,也能通过把初始位移 (poke) 作为条件直接对反事实动作做规划筛选;最终用同一个 665M 参数的模型同时支撑 "预测未来" 与 "为达到目标选择动作" 两类任务,目标是把当前视频生成模型在 OWM 上 DNF 的耗时吞吐量提升约 4–5 个数量级,并在台球规划基准上接近物理引擎 (Oracle) 的 84% 准确率。

与已有工作不同的是,现有的 '图像到轨迹' 类工作(如 Generative Image Dynamics、Track2Act、Flow Poke Transformer)虽然省去外观生成,但它们要么要求模型在不同轨迹上做大规模并行推断(效率差),要么要求模型一步到位地预测整条轨迹(无法显式建模长程多接触交互)。本文把视觉税推到极限——只预测用户定义的 K 个 2D 点而非任何整张图像,并在时间和轨迹两个维度上都做 autoregressive 因式分解,配合高效 fused attention 层、flow matching head 与尺度级联,实现 "像素外 / 物理内" 的细粒度时空推理:从一张台球桌图片出发可在数秒内枚举约 2200 条候选轨迹,比视频模型(0.14 samples/min)快 3–4 个数量级,并把这一稀疏反事实枚举直接转化为决策信号,得出 78% 的击球入袋准确率,逼近物理引擎上界 (Oracle) 的 84%。

核心方法

直觉层面,Myriad 不预测 '未来长什么样',而预测 '未来怎么动':输入一张 I_0 和 K 个查询点 x_0,每步输出所有点的位移 $\Delta x_t$,在线累加得到未来轨迹 $x_t$。技术上是一个 '图像条件自回归扩散':用 DINOv3-L/16 作为图像编码器得到空间特征 $E_{\text{img}}$,每个 (t, i) 对生成一个 motion token——它由 (1) 当前位移的 Fourier 嵌入、(2) 唯一随机单位球 ID、(3) 在原始 $x_0$ 与当前位置 $x_t$ 处的图像特征拼成,再过一个小 MLP;motion token 沿轨迹维度和时间维度都做严格因果自回归,融入共享时空 RoPE 位置编码;最后由一个轻量 Flow Matching head 在 50 步 ODE 求解下输出下一时刻位移分布;可选地,输入一个 'poke'(即前两帧的真实位移)作为初始运动提示。整个架构 665M 参数,KV cache 跨时间步、轨迹条件缓存,使得 32 步 rollout 平均吞吐量约 2200 samples/min,对比 Wan2.2 的 0.141 samples/min 高出四到五个数量级。

核心创新是把 '开放集未来预测' 重写为 '在稀疏轨迹上的逐步自回归分布推断',并配合三项工程优化:(a) **完全避免视觉税**——不生成像素,只预测用户指定的 K 个 2D 点,跳过 VAE 解码和稠密张量运算;(b) **时间和轨迹双维度因果因式分解**——$p(x_{1:T}\mid x_0, I_0)=\prod_t \prod_i p(x_t^{(i)} \mid x_t^{<i}, x_{<t}, I_0)$,使得时间维 KV cache 可复用、轨迹维并行解码;(c) **专为轨迹设计的 FM head + 尺度级联**——用 flow matching 取代 GMM 头以处理多模态重尾分布,并在输入端用 $\tanh(s\cdot \Delta x)$ 把对数尺度的多个 s 拼起来同时表示微小抖动和大跳跃,这在视觉税为零的情况下对稳定训练至关重要。本质区别:FPT (Flow Poke Transformer) 单步跳跃到 t=T,需要无限大模型容量来隐式推算所有中间碰撞;而 Myriad 把这些碰撞显式展开成短时局部可预测的微步,把 '难' 分摊到 16/50 个简单步骤。

方法步骤详情

完整流程分训练与推理两段。训练时:第一步把原始互联网视频居中裁剪到 512×512,用 TAPNext 在 1024 个随机种子点上跟踪得到 16 个连续帧的 2D 轨迹作为伪真值;额外还有 ~1.5M 视频用 V-DPM 做 3D 跟踪再反投影到第一帧得到 '纯场景运动' 标签。第二步用 DINOv3-L/16 编码 I_0 得到图像特征图,对每个 (t, i) 在 $E_{\text{img}}$ 上双线性采样 'what'($x_0$ 处)与 'where'($x_t$ 处)特征。第三步构造 motion token = MLP(Fourier($\Delta x_t^{(i)}$) ‖ shared_ID ‖ what ‖ where)。第四步把所有 motion token 与图像 token 拼成一个序列 $h = [h_{\text{image}} | h_{\text{motion}}]$,构造 mask(图像只关注自身;motion 关注自身 + 图像 + 因果过去 motion),用 24 层 width=1024 的 fused parallel transformer block 做并行 self/cross-attention(每个块内部把 QKV+FFN-up 共享一次 Norm+Proj,再以一次 Down 合并出 +FFN,避免 6 个独立 kernel launch)。第五步 FM head(3 层 FFN,width=1024)以 $z_t^{(i)}$ 和 flow time $\tau$ 为条件,通过 adaptive norm 输出 ODE 速度场 $v_\phi(\Delta x_{t,\tau}^{(i)} \mid z_t^{(i)})$;输入端用对数尺度级联 $[\tanh(s_1 \Delta x_{\tau, 1}) | \ldots | \tanh(s_L \Delta x_{\tau,L})]$ 把跨数个数量级的位移同时稳定表示。第六步用 ELBO $\mathcal{L}_{FM}=\mathbb{E}_{\tau,\Delta x_0,\Delta x_1}\|v_\phi-\frac{\partial \Delta x_\tau}{\partial \tau}\|_2^2$ 端到端训练 400k 步(16×H200,20 小时),bfloat16 混合精度,全局 batch 128,K=16,T=16,AdamW 峰值 lr=3e-5 加 5k 步 warm-up。推理时:用 Euler 解 ODE 50 步解码 $\Delta x_t$,在线更新 $x_{t+1}$,时间维度 KV cache 持续累积,轨迹维度并行;可选先用前 2 帧真值作为 'warm-up poke' 提示。

技术新颖性

技术新颖性主要有四点。第一,把开放集未来预测问题从 '生成未来帧' 转换为 '在稀疏轨迹上的逐步自回归扩散',这是表征层面的范式转移(visual tax → zero),论文通过 FPT 与 GAMED 等 baseline 的对比实验和理论动机('人脑不画未来而追踪关键变化' [13,52])论证其必要性。第二,提出 fused parallel transformer block,把六次 kernel launch(QKV-Proj、Attention、O-Proj、FFN-up、FFN-down 各加 Norm)合并为两次 Launch + 一个 Activation,配套精心设计的 mask,使得 32 步 rollout、batch 4/16 时获得 2× / 3.7× 加速。第三,提出输入尺度级联,用对数级 $\exp(\text{linspace}(\log 0.1, \log 1e5, L))$ 的 $s$ 对输入做 $\tanh(s\cdot \Delta x)$ 饱和,既能表示亚像素细微抖动,也能表示数十像素的快速跳跃,弥补了 FLUID/TiT 等既有 FM head 设计仅适用于图像类轻尾分布的缺陷。第四,引入随机单位球 ID $\text{id}_{traj} \sim U(S^{d-1})$ 替代固定 codebook,论文证明零样本可以从训练时的 16 条轨迹推广到任意数量的查询点(Traj. Emb ablation, Tab. C)。

Motion Token Construction
Figure 2: Motion Token Construction
Positional Encoding Scheme
Figure 3: Positional Encoding Scheme
Fast Reasoning Blocks — (a) Normal Layer vs (b) Our Fused Layer
Figure 4: Fast Reasoning Blocks — (a) Normal Layer vs (b) Our Fused Layer
Our attention mask
Figure 5: Our attention mask
Posterior FM Head
Figure 6: Posterior FM Head
Value distribution
Figure 7: Value distribution

实验结果

论文在四类基准上呈现了 Myriad 的一致优势,且用具体数字量化:**(1) OWM (开放世界)**:Myriad 665M 在 Best-of-5 的总体 MSE 为 0.029,优于 MAGI-1 (4.5B, 0.037)、优于 Wan2.2 (14B, 0.141)、优于 SkyReels V2 (1.3B, 0.304)、优于 SVD (1.5B, 0.714);在 5 分钟时间预算下的 Best-5min 进一步压低到 0.013,比 Wan2.2/CogVideo/SkyReels 高出一个数量级(CogVideoX 直接 DNF),吞吐量 2200 samples/min 对比视频模型 0.051–0.714 提升 ≈ 3000×。**(2) PhysicsIQ**('solid mechanics' 子集 97 样本):Best-of-5 的 MSE=0.115、Best-5min=0.045,与最强视频模型持平;这里视频模型常预测过于简化的轨迹,而 Myriad 能抓真实碰撞细节。**(3) Physion**(134 样本):Best-of-5 MSE=0.048、Best-5min=0.020,反过来,视频模型在此简单刚体场景上幻觉出过度复杂运动,而 Myriad 保留简洁物理。**(4) OWM 子集分解 (Tab. F)**:Myriad 在 Rigidity(rigid 0.007 / non-rigid 0.016)、Multi-Agent (0.017 vs SkyReels 0.075)、Free-Will (0.014 vs SVD 0.069) 等所有切片上均排名第 1 且数量级快于基线;Multi-Agent 场景 Best-5min 仅 0.017,证明其在多体交互下也保持分布准确性。**(5) 台球规划 (Tab. 2)**:在零样本动作选择任务中 Myriad 取得 78% 准确率,逼近物理引擎上界 (Oracle) 的 84%,是 I2V baseline(16%)、Full Trajectory Diffusion(8%)、Flow Poke Transformer(4%)的 5–20×;吞吐量 753 actions/min 是 I2V 模型 20 actions/min 的 38×,是 FPT 13422 actions/min 但准确率仅 4% 的对比性证据。**(6) 校准性 (Fig. 12)**:从像素级误差 (~1/512) 起,预测的不确定度标准差与真实 EPE 在 log-log 空间呈线性,即模型的不确定度是良好校准的,可以作为决策时的置信信号。**(7) 消融**:FM head vs GMM head 同时获得更快收敛(Fig. 6 右)和更低 Best-5(0.029 vs 0.110),尺度级联开启后 Best-5 从 0.033 → 0.029;fused block 提供 2× / 3.7× 加速;ETH-UCY 轨迹预测零样本下与 Trajectron++ 接近,证明方法对经典任务亦具有竞争力。

Open-world & Physical Motion Prediction
Table 1: Open-world & Physical Motion Prediction
Planning Billiard Shots through Future Exploration
Table 2: Planning Billiard Shots through Future Exploration
OWM Composition
Figure 8: OWM Composition
(a) Diverse Actions from a Single Image; (b) Object Coherence
Figure 9: (a) Diverse Actions from a Single Image; (b) Object Coherence
Time-Accuracy Tradeoff on OWM
Figure 10: Time-Accuracy Tradeoff on OWM
Planning a Billiard Shot
Figure 11: Planning a Billiard Shot
Posterior Uncertainty vs. Error
Figure 12: Posterior Uncertainty vs. Error
查看结构化数据
任务指标本文基线提升
OWM 开放世界运动预测 (95 视频混合) Best-of-5 / Best-5min MSE (越低越好) 0.029 / 0.013 (665M 参数,2200 samples/min) MAGI-1 0.037/0.066 (4.5B), Wan2.2 0.141/DNF (14B), CogVideo-X 1.5 0.051/DNF (5B), SkyReels V2 0.304/0.068, SVD 1.1 0.714/0.119 在参数小 7×–20× 的同时 Best-of-5 提升 1.3–24×,Best-5min 提升 5× 且对手大批量 DNF。
PhysicsIQ (solid mechanics, 97 样本) Best-of-5 / Best-5min MSE 0.115 / 0.045 MAGI-1 0.126/0.169, Wan2.2 0.116/DNF, CogVideo-X 0.100/DNF, SkyReels 0.128/0.137, SVD 0.138/0.241 最佳 Best-5min 性能 0.045 vs 第二名 0.169 (≈3.7×), 视频类几乎全部 DNF。
Physion (刚体, 134 样本) Best-of-5 / Best-5min MSE 0.048 / 0.020 MAGI-1 0.061/0.081, Wan2.2 0.069/DNF, CogVideo-X 0.063/DNF, SkyReels 0.069/0.084, SVD 0.070/0.147 Best-5min 0.020 显著优于 SkyReels 0.084 与 SVD 0.147 (4×–7×)。
台球零样本动作规划 (击球入袋) Accuracy (≤1 球径判正确) + 吞吐量 78% accuracy, 753 actions/min (665M) Physics Oracle 84%/55162; I2V Diff (image) 16%/19.8; I2V poke-cond 16%/20.4; AR I2V 12%/22.2; Trajectory diffusion 8%/13.4k; FPT 4%/160.8 在接近 oracle 上界的精度下 (78% vs 84%),吞吐量是 I2V Diffusion 的 38×,且精度比 FPT (4%) 高出 19×。
ETH-UCY 经典轨迹预测 (零样本) Best-of-20 ADE/FDE (越小越好) ETH 0.31/0.80; Hotel 0.17/0.30; Zara01 0.53/1.21; Zara02 0.40/0.91 Trajectron++ ETH 0.39/0.83; Hotel 0.12/0.19; Zara01 0.15/0.33; Zara02 0.11/0.25 在 ETH/Hotel 上与最强专用模型 Trajectron++ 接近甚至更优(Hotel 0.17/0.30 vs 0.12/0.19),但是在 Zara 上弱于专用模型,强调 'in-domain 微调' 仍是差距来源。

局限与改进

作者明确承认两个主要局限。第一,**静态相机假设**:Myriad 默认摄像机固定,避免出现场景运动 vs 自运动的混淆;这一假设简化了评估与可解释性,但限制了它在自动驾驶、机器人第一视角或动态场景下的直接应用。作者在补充材料中报告了一种缓解方案——用 V-DPM 估计相机位姿并把 3D 轨迹重投影回第一帧,但 '联合预测 ego 运动与场景运动' 仍是开放问题。第二,**依赖伪真值跟踪器**:训练轨迹由 TAPNext / V-DPM 等点跟踪器生成,会继承这些模型的偏差、遮挡失败与抖动;同时 OWM benchmark 也必须 '人工核实' 跟踪质量以保证评测可靠。其他隐含观察:(a) 模型对 K(查询点数量)敏感——零样本外推到 24 个查询点时 EPE 从 0.00141 上升到 0.00263(Tab. C),并非完全尺度无关;(b) 在 ETH-UCY 这种 '给定 2D top-down' 的专用 benchmark 上零样本仍逊于任务专用模型(Zara01 0.53 vs 0.15),提示在完全 in-domain 数据上微调仍有空间;(c) FM head 的 50 步 ODE 求解仍是计算主体 (Tab. B:1 NFE=0.00361 vs 50 NFE=0.00138),意味着减少 NFE 仍有提速空间,但收益递减。

独立分析的弱点

独立分析而言,可以归纳出三类弱点和改进方向。第一,**稀疏先验与长时记忆的张力**:当 K 较小(如 1–4 个查询点)或查询点彼此远离(不同物体)时,全序列 attention 难以捕捉跨物体的耦合交互(碰撞、遮挡);可以通过显式加入 pair-wise prior latent 或 scene-level context token 改进。第二,**对动态视角/相机扩展的薄弱**:虽然 V-DPM 重投影可以学到 '去相机后的' 轨迹,但推理时仍假设静态相机;如果要拓展到视频流或自运动预测,需要联合建模 trajectory+ego pose 的双自回归流,例如把相机位姿也当作 token 注入,并把 OWM-equivalent benchmark 扩展到 KITTI/TruckScenes 等动态相机数据集。第三,**训练数据偏差的传递**:训练用的 TAPNext pseudo-GT 在遮挡、快速运动、纹理缺失场景下本身会失效;可以引入不确定性加权的 self-training、在测试时对难以跟踪片段进行自蒸馏,或在 FM head 中加入一个 '置信 dropout' 让推理时自动识别此类样本。

未来方向

作者明确提出的未来方向:将方法推广到动态相机场景进行 ego+scene 联合预测,以及引入更鲁棒的(去跟踪器偏差的)训练信号。基于论文成果可延伸的方向还包括:(a) 与世界模型或 RL agent 闭环集成——Myriad 提供的快速反事实 rollout 可作为 MCTS 类算法的 transition model,例如在机器人操控中将 16 个查询点对应若干抓取坐标;(b) 多智能体扩展:把 K 个查询点扩展为多 agent 的策略模块,每条 trajectory ID 对应一个 agent,从而把 social forecasting 与 OWM 桥接起来;(c) 视频+稀疏轨迹的混合表征:用像素 diffusion 生成关键帧做 '远点 long-term anchor',再用 Myriad 填补中间稀疏轨迹,缓解极端长程(>10s)下的不确定性漂移;(d) 3D 与物理一致性:将训练目标附加接触守恒、能量守恒约束或基于可微物理的 residual 把 FM head 升级为约束流匹配 (constrained flow matching),让分布本身保证物理可能。

复现评估

复现门槛在中等偏高。第一,**代码与数据** — 项目页 (http://compvis.github.io/myriad) 提供代码与 OWM benchmark(论文承诺将 OWM 公开),作者在补充材料 Sec. A 中详细列出了训练数据来源(10M 网络视频 + 1.5M V-DPM-重投影 + billiard simulation)、伪真值来源(TAPNext, V-DPM)以及完整超参 (Tab. A),并给出 PyTorch+Diffusers 风格的实现细节。第二,**算力门槛** — 主模型训练 400k 步 / 16×H200 / 20 小时 ≈ 320 H200-hours;billiard 与 ablation 变体略小。即使只微调,1 张 H200 也需数日,业余研究者难以纯本地复现完整实验。第三,**数据预处理** — 10M 网络爬取视频的去重、静态相机筛选、1024 种子点 seed、TAPNext 推理本身就需要显著算力;论文没有发布这些 raw data,是潜在的可复现性瓶颈。第四,**依赖与精度** — PyTorch + DINOv3-L/16 + AdamW bf16,包含一个自定义 fused parallel transformer 和全可缓存 adaptive-norm FM head;中等规模 GPU(H100 80GB)应足以运行 inference 但训练仍建议多机。第五,**评估** — OWM 评测脚本与 Best-of-N/Best-5min 评测协议在 Appendix A.6 已明确,并附 billiard planning 的 oracle 实现;评测可独立完成。