---
title: "读懂 EHT 成像与动态重建的偏好"
description: "逐图阅读固定版本论文，区分真实观测与模拟检验，追踪空间和时间偏好怎样进入重建。"
canonical_url: "https://sandboxphysics.com/zh/stories/interferometric-imaging/read-eht-imaging-and-dynamic-priors/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/interferometric-imaging/read-eht-imaging-and-dynamic-priors/"
last_modified: "2026-10-28"
---

[从电波延迟走到 EHT 干涉成像](https://sandboxphysics.com/zh/stories/interferometric-imaging/)

AC05 · Paper guide

# 读懂 EHT 成像与动态重建的偏好

Sandbox Physics · 2026-10-28

把《M87 事件视界望远镜首批结果 IV》的 arXiv:1906.11241v1 放在旁边。从图 1、3、5、7，顺着采样基线走到成像选择。再打开 Johnson 等人的《干涉测量动态成像》arXiv:1711.01286v1，看看光源在阵列转完以前就开始变化时，问题增加了哪一层。

## 开始前需要什么

先读 L1–L3，熟悉可见度振幅、有向闭合相位和脏波束。阅读固定在这里列出的论文版本，不把后来其他 EHT 结果混进同一组图。

## 图 1 先让我们看见缺了哪些测量

先看 Paper IV 图 1 的空间频率坐标轴，再去看后面的图像。一个位置表示用波长计的投影基线。离中心越远，对应越细的角度；方向决定它对光源哪一部分形状敏感。已有轨迹之间的空隙，也要认真看。

如果真有一面填满地球尺度的口径，就能测到连续得多的站间距离。EHT 的站点很稀疏，最长基线带来极高分辨能力，也留下大片空缺。波束核心的宽度给出一种有用尺度，却不表示该尺度上的每种结构都同样分得清。

接着用同样的覆盖去比较观测日和候选源。重建形状改变，可能是天空变了，也可能是采样方向变了，或方法对缺失频率的反应变了。这张图提示后面必须区分哪些可能。

对应阅读：EHT Collaboration, First M87 Results IV, arXiv:1906.11241v1 · NRAO Essential Radio Astronomy, Chapter 3, aperture synthesis

## 图 3 读的是站点不能随意改变的相位

图 3 为若干三角形画出闭合相位随格林尼治平恒星时的变化。选一个三角形，沿着标签读各面板。每个点组合三条同时测得的有向基线，不是某台天线自己的相位。可见度接近零时，相位对噪声特别敏感，角度很大也要结合振幅和误差一起判断。

回到式 7 与第 2.1.2 节，找到站点增益模型。三角形乘积能消去站点相位，振幅里的站点增益仍然存在。闭合振幅另用多条基线的比值消去振幅增益。读数据项时，要把这两类标定问题分开。

教程 CSV 的时间列是十进制 UTC 小时，图 3 用的是 GMST。要叠图，必须正确转换时间、匹配观测日、三角形方向和平均方式。教程只在原生时间戳下重算一个真实三角形，没有把它冒充成整条论文时间序列。

对应阅读：EHT Collaboration, First M87 Results IV, arXiv:1906.11241v1 · EHT 2019-D01-01 M87 data release, pinned commit 27ae76a

## 图 5 给算法一片真值已知的天空

图 5 用合成源模型，以及这些模型会产生的可见度。先看已知输入，再顺着采样走到重建。这里真值可知，因为源由作者构造。这项检验回答：经过同样类型的不完整观测，算法能把这些形状恢复到什么程度？

比较时先匹配角分辨率。重建特别锐利，可能只因为含有观测支持不了的细节；把参考模糊到共同尺度，更容易看清大结构是否恢复。也要比较预测信号。在一个训练源上恢复得好很有用，但还要换其他可能的源继续测。

真正的 M87 观测随后进入另一阶段。将模拟检验单独放清楚，作者就能调试、评价成像方法，而不必先认定某一张 M87 重建就是真值。

对应阅读：EHT Collaboration, First M87 Results IV, arXiv:1906.11241v1

## 式 14 和图 7 把成像偏好摆出来

式 14 把数据符合程度与带权的图像正则项放在一起。权重常叫超参数，决定程序有多偏好平滑等性质。附录 A 给出数学定义。比较权重前，先读符号和归一化；离开这些规定，一个数字没有通用含义。

图 7 改变 MEM 与 TV 权重，展示一组重建。MEM 是基于熵的偏好，TV 用总变差惩罚局部图像变化。先找到能合理符合观测的设置，再看哪些图像特征在改变设置后仍然留下。

这些面板来自不同成像选择，不是从后验分布随机抽出的图。它们之间的视觉分散程度，不能直接当作已经标定的误差带。论文还结合模拟检验与独立流程比较，第 6.2.2 节有助于理解这些选择怎样接受评价。

对应阅读：EHT Collaboration, First M87 Results IV, arXiv:1906.11241v1

## 地球转动也会混合时间与空间

翻到动态成像论文第 2.3–2.4 节。一条基线上的信号变了，可能是它走到了新的 $(u,v)$ 位置，也可能是光源恰好在这时改变。静态成像把两种变化都交给同一张固定图像。源变得足够快时，这个假设就可能把时间变化塞成空间结构。

论文讨论怎样用接近的空间频率测量识别变化，再去重建一段序列。看一部漂亮动画之前，先问每帧有多少直接观测。某帧很少有数据，甚至根本没有时，连接邻帧的假设便格外重要。

对应阅读：Johnson et al., Dynamical Imaging with Interferometry, arXiv:1711.01286v1

## 式 3 说明相邻画面怎样相连

动态目标函数在数据项、普通图像正则项之外，再加时间正则项。下面用教学记号把三者分开；论文自己的式 3 给出具体符号和归一化，第 3.1–3.3 节分别定义时间选择。

一种偏好要求相邻帧连续，一种偏好让每帧只小幅偏离时间平均图像，还有一种用流动模型偏好近似稳定的亮度输运。缺少数据时，它们会给出不同画面。时间权重接近零，各帧更独立；强烈偏好静止的权重若取太大，也会压掉真实变化。

$$
J(\{I_t\})=D(\{I_t\})+\alpha R_{\mathrm{image}}(\{I_t\})+\beta R_{\mathrm{time}}(\{I_t\})
$$

这里用统一教学记号表示数据误差、图像偏好和时间偏好。不同时间正则项仍然对应不同的源变化假设。

对应阅读：Johnson et al., Dynamical Imaging with Interferometry, arXiv:1711.01286v1

## 图 7–12 要逐一看清观测身份

动态论文图 7–10 使用模拟黑洞系统观测，其中也有只用振幅和闭合信息的重建。已知模拟画面允许直接检查恢复程度。画面中的运动属于合成实验，不是 2017 年 EHT 实测黑洞电影。

图 12 则使用真实的 7 毫米 VLBA M87 喷流跨历元观测。看清图注里的波长、阵列和时间间隔。它检验的是另一种源尺度和采样条件，不能与教程里的 1.3 毫米 EHT 摘录混在一起。

读下一篇成像论文时，先记五件事：实际测量什么、采到哪些空间频率、空间和时间上加入什么偏好、做了哪些已知真值检验、改变假设后哪个特征还在。写清这几项，就能沿着图像追到一个可以检查的科学结论。

对应阅读：Johnson et al., Dynamical Imaging with Interferometry, arXiv:1711.01286v1 · EHT Collaboration, First M87 Results IV, arXiv:1906.11241v1

## [重现一次空间偏好比较](https://sandboxphysics.com/zh/labs/stellar-interferometer/)

1. 选“基线挤在一个方向”，比较同一组信号的两种正则化重建。
2. 加入地球自转覆盖，找到改变观测条件后仍保留的源特征。
3. 用这个静态例子说明，动态论文为什么还得交代不同时间帧怎样相连。

## 练习与解释性答案

1. 某个特征在同一组正则化扫描里始终出现，下一步问什么？

   展开解释

   看看使用不同正则项的独立方法是否也保留它，以及实际采样下的合成检验能否恢复类似特征。同一组扫描共用代码、数据选择和许多假设，一致性有意义，也有共同来源。

2. 动画在两次观测之间多出一帧，它的信息来自哪里？

   展开解释

   来自其他帧的测量和连接各帧的时间正则项。应把这一帧理解为指定假设下的推断插值，再检查换一种时间偏好后，想解释的运动是否仍然存在。

[上一篇：拟合一个真实 EHT 闭合三角形](https://sandboxphysics.com/zh/stories/interferometric-imaging/fit-a-real-eht-closure-triangle/)

## 来源与计算范围

1. [NRAO Essential Radio Astronomy, Chapter 3, aperture synthesis](https://www.cv.nrao.edu/~sransom/web/Ch3.html)

   小视场傅里叶关系与地球自转综合，重点为式 3.180–3.200。浏览器采用球形地球、单一波长和有限像素天空。

2. [EHT 2019-D01-01 M87 data release, pinned commit 27ae76a](https://github.com/eventhorizontelescope/2019-D01-01/tree/27ae76aacb9067a83388395c3826012b97aac371)

   实测数据为 2017 年 4 月 10 日 EHT-HOPS SR1 网络标定后的 Stokes I，低频段中心频率 227.0707 GHz。教程保留发布者 CSV，并选取同时刻 AA–AZ–LM 三角形。数据采用 PDDL 1.0。

3. [EHT Collaboration, First M87 Results IV, arXiv:1906.11241v1](https://arxiv.org/abs/1906.11241v1)

   固定为 2019 成像论文的 v1，DOI 10.3847/2041-8213/ab0e85。阅读图 1、3、5、7，站点增益式 7、正则成像式 14 及附录 A。本页链接原图，不另行复制论文图片。

4. [Johnson et al., Dynamical Imaging with Interferometry, arXiv:1711.01286v1](https://arxiv.org/abs/1711.01286v1)

   固定为 2017 年 11 月 3 日版本，DOI 10.3847/1538-4357/aa97dd。第 2.3–3 节与式 3 说明时间正则化；图 7–10 使用模拟黑洞观测，图 12 使用跨历元的真实 7 毫米 VLBA M87 数据。