GPU 加速成像流水线能否让艾伦望远镜阵列实现实时 SETI 搜索？

成像以多条视线取代波束合成

传统的技术信号与暂现源搜索通常处理波束合成后的时频数据产品。本文则为艾伦望远镜阵列（ATA）的相关后 UVH5 可见度数据构建了一条干涉成像路径，使时域搜索可以同时在大量图像像素上进行，而不必先提取某个已知位置再做时间分析。论文给出的动机是：成像提供了空间信息，有助于定位候选信号，并区分天体辐射与地面干扰。

该流水线从相关可见度生成 Stokes-I 脏图，将四分位距抑制（IQRM）标记纳入成像权重，对跨越 287.5 MHz 的 48 个频率子带成像，从每个子带和像素中减去整段观测的时间均值，施加基于 IQRM 的子带掩蔽，并使用短历史图像环形缓冲区执行盲试 DM 搜索。色散搜索直接在频率分辨图像上进行，因此每个图像像素都被当作一条独立视线来搜索。

作者将这一路径定位为对波束合成后端的补充而非替代：波束合成的宽带脉冲搜索和实时深度学习方法可以提供第一轮探测，随后用干涉成像检查有希望的候选信号。该成像流水线旨在通过保留空间信息并同时搜索大量像素来支持这类后续分析。

一次脉冲星观测中人为注入脉冲

ATA 在 2026 年 1 月 21 日指向明亮脉冲星 PSR J0332+5434（B0329+54）进行了约 10 分钟的观测。该观测包含七个相邻谱窗（SPW），以 UVH5 格式记录；每个 SPW 有 192 个通道，通道宽度 0.5 MHz，即每个 SPW 为 96 MHz，总记录带宽 672 MHz。数据集包含 6104 次积分，周期约 98.29 ms，以及 406 个基线产品，其中 378 个互相关用于成像，28 个自相关被排除。每次积分独立生成 128x128 像素网格、像元 45 角秒的脏图。

早期的连续谱脏图测试使用单个 SPW（SPW #3），其射频干扰相对较少，脉冲星在该图像时间序列中周期性出现，与其约 0.7145 s 的脉冲周期一致。但该脉冲星的色散量较低，单个 96 MHz SPW 内的带内色散扫过仅占一次积分的很小一部分，无法有力检验图像域 DM 搜索。因此多频段结果改用前三个 SPW，并注入一个明亮、强色散的脉冲信号。

这三个 SPW 提供 576 个连续频率通道，跨越 287.5 MHz，从 1.164 到 1.452 GHz，时间与基线网格与上述相同；其余四个 SPW 留待未来更宽带宽的测试。注入源是一个 DM = 150 pc cm-3 的周期性色散脉冲信号，直接在可见度域加入：在每个选定的（时间、频率、极化）样本上，向每个基线产品加入一个实值复常数，幅度为原生未定标可见度单位下的 10^6，虚部为零。未指定目标输入信噪比。由于加入的相量在各基线上相同，对应的脏图响应在相位中心达到峰值，即 128x128 网格上的像素 (64,64)。每个脉冲的固有宽度为一次积分，约 98.29 ms，脉冲每 100 次积分注入一次，周期约 9.829 s。

图像面信噪比与时间信噪比回答不同问题

论文报告了两个读者不应视为同一测量的信噪比数值。在图像域去色散之后，注入脉冲的中位空间图像面信噪比约为 52，是相对于对应 DM = 150 pc cm-3 去色散图像中无源像素的空间 RMS 测得的。它比较的是源与单张去色散图像中的空间噪声。

相比之下，盲搜索使用从每个像素和 DM 的完整去色散时间序列计算出的时间统计量对候选排序，信噪比定义为最大值减中位数，再除以 1.4826 倍的中位绝对偏差。在三 SPW 测试集上，盲搜索在相位中心找回注入信号，最佳拟合 DM = 150 pc cm-3，时间 MAD 信噪比约为 286；相邻的 140 和 160 pc cm-3 试值约为 272。该时间值比较的是最强脉冲与同一像素上的时间变化，排序选出的是注入周期脉冲串中最强的单个脉冲。

这一区分很重要，因为两个数值描述的是不同的分母和不同的问题：一个是单张图像内的空间对比度，另一个是整段观测中单个像素的时间对比度。读者在将该流水线与波束合成搜索比较时，应先确认引用的是哪个统计量，再对相对灵敏度下结论。

延迟按通道而非按子带计算

频率相关的到达时间遵循冷等离子体色散定律。延迟对 576 个原生频率通道逐一独立计算，而非按子带计算，每个通道频率都与同一个全局参考比较：三个 SPW 中最高的通道中心 1.4515 GHz。因此较低频率通道相对于该参考被延迟。延迟在拼接后的 576 通道频率轴上计算一次，然后在每个 SPW 文件内施加，因此三个频段共享同一个色散参考，而不是各自独立的文件内参考。

连续延迟在注入前四舍五入到最近的整数次积分。对于 DM = 150 pc cm-3，最大延迟来自最低通道中心 1.164 GHz，约为 0.164 s，即约 2 次积分。

在盲搜索中，试色散量为 DM = 10, 20, ..., 200 pc cm-3，共 20 个试值。子带延迟使用相同的冷等离子体关系，但在每个子带中心频率上计算，而非注入时使用的原生通道中心，连续延迟同样四舍五入到最近的整数次积分。施加这些整数时间平移后，子带图像在每个像素上对可用子带取平均，忽略 NaN 平面，形成去色散立方体。

同一套射频干扰标记有两种用法

射频干扰标记遵循 IQRM，它识别可见度统计相对于邻近通道异常的频率通道。掩模生成在拼接后的三 SPW 可见度集上离线执行一次：生成 Stokes I；通过非相干地对跨互相关基线的可见度幅度求和构建时频动态谱；在频率上拼接三个 SPW；在 1000 次积分的块内将每块归约为每通道谱统计量，此处为时间标准差，并施加 IQRM 邻域比较，阈值 3.0，半径约为通道数的 10%（576 通道时为半径 58）；然后生成布尔掩模。在三 SPW 测试集上，联合掩模占用约为时频样本的 29.1%，576 个通道中有 246 个至少被标记一次。

在成像期间，被标记的单个通道在该次积分中对所有基线赋予零权重，而同一子带内未标记的通道仍对脏图有贡献。可见度幅度不被改写，UVH5 文件也不被修改。对于 DM 搜索，作者采用更保守的准则：如果某个子带内的 12 个通道中任意一个在某次积分被标记，整个 128x128 子带平面在该次积分被排除出试 DM 组合。在该数据集中，48 个子带中有 20 个在每次积分都被排除。

作者将该保守规则与两种较宽松的替代方案比较：只要 12 个通道中至少一个未被标记就保留该子带，以及按未标记通道比例对每个子带加权。两种替代方案都保留了多得多的数据，约为子带-时间样本的 89.7%，而保守规则约为 50.3%，但在该数据集上它们让射频干扰主导盲搜索并掩盖了注入脉冲。因此所有报告结果都采用保守的任一通道拒绝规则，作者也指出更温和或空间相关的射频干扰处理仍有待探索。

短图像历史避免存储完整立方体

存储完整的子带立方体 48 x 6104 x 128 x 128 float32 约需 17.4 GB。但对于当前的 DM 网格，最大整数延迟仅为 2 次积分，因此只需要近期子带图像的短历史。流水线使用深度 R = 8 的环形缓冲区，满足 R 大于最大延迟加一，足以形成所有试 DM 对齐，而无需在内存中保留完整的频率×时间×图像立方体。这将常驻图像历史降至约 30 MB。

在环形缓冲区上运行时，试 DM 图像组合比扫描完整存储立方体快约 5.2 倍。作者指出，该缓冲策略受 STRIDE 中流式图像集和环形缓冲区概念的启发，但这里的环形缓冲区仅作为形成色散对齐图像组合的短历史缓冲，并非 STRIDE 的完整实现。

核心科学处理满足周期，输出写入不满足

集成到 CyberEther 流式框架后，核心科学阶段每次积分约需 63.3 ms，低于约 98.29 ms 的原生周期。论文给出的分阶段耗时大致为：回放与 I/O 9.66 ms，IQRM 权重施加 0.09 ms，GPU 子带成像 33.05 ms，均值扣除 0.18 ms，子带掩蔽 0.18 ms，环形缓冲区管理 0.10 ms，CPU 试 DM 组合 20.08 ms。只有脏图成像在 GPU 上运行；回放、IQRM 权重施加、均值扣除、基于 IQRM 的子带掩蔽、环形缓冲区管理和试 DM 组合在当前流水线中仍是 CPU 阶段。

写出完整的去色散数据产品会增加每次积分约 50 至 55 ms，使总处理时间升至每次积分约 124 ms。因此作者将科学处理时间与输出 I/O 区分开来，并明确不把包含 I/O 的运行视为实时。这是论文中的关键阈值：科学路径落在一次积分之内，但当前配置下的完整流水线没有。

论文还报告，在单个 96 MHz SPW 上，端到端连续谱成像（每次积分一张脏图，含显示映射）每次积分约需 2.9 ms，而原生周期为 98 ms。该数值适用于单 SPW 连续谱成像路径，而非完整的三 SPW 图像域 DM 搜索，不应替代 63.3 ms 的科学总耗时。

记录回放尚不是实时相关器后端

当前实现将已记录的相关后可见度数据流经流水线，而不是作为实时 ATA 相关器后端运行。它依赖在整段观测上预先计算的时间均值和 MAD 噪声估计。对于实时运行，这些量需要仅用截至当前积分的数据来估计。候选排序目前在回放后离线执行，而非在流式处理中执行，且当前候选排序阶段未纳入周期性，尽管作者已开发出可浏览去色散时间序列并运行 FFT 和折叠后续分析的工具。

作者还指出，GPU 子带成像被分成 48 个相对较小的操作，伴随重复的数据传输，因此 GPU 利用率尚不理想。迄今的评估侧重于受控注入信号，而非真实数据上的盲搜索。保守的任一通道子带掩模在该数据集上效果良好，但更温和或空间相关的射频干扰处理仍有待探索。

计划中的工作包括：批处理 GPU 子带成像以提升占用率，超越当前 48 个串行子带操作；为实时处理运行背景与噪声估计器；在流式处理中完成候选探测而无需大型去色散数据产品；以及对去色散的逐像素时间序列进行周期性评分。作者还计划对已知、色散更高的脉冲星（包括 DM 约 158 pc cm-3 的 B1933+16）的单个脉冲测试完整的图像域 DM 搜索，并开展覆盖一系列 DM 和信噪比的注入测试，以刻画找回性能。

一个经验证的构件，而非完成的巡天系统

论文的具体主张范围狭窄且可检验：在受控的三 SPW 注入数据集上，流水线在相位中心找回色散脉冲源，DM = 150 pc cm-3，时间信噪比约 286，核心处理在约 63 ms 内完成，低于约 98 ms 的周期。作者明确指出，写出大型去色散数据产品不符合该预算。

对于将本工作与波束合成搜索后端权衡的读者，决策相关的比较不是原始灵敏度，而是该流水线在巡天链条中的位置。论文设想将波束合成第一轮搜索与该成像流水线配对，使波束合成产品中识别出的候选事件可以触发提取相应的可见度或图像数据段，用于图像域定位与确认。这里展示的价值在于，成像与图像域 DM 搜索阶段可以在已记录数据上于 ATA 积分时间内运行，剩余缺口是实时背景与噪声估计、流式候选排序以及输出 I/O。

因此，解读该结果最有用的方式是将其视为一个带有实测时间预算的经验证构件，而非已完成的实时 SETI 系统。63.3 ms 的科学总耗时与约 124 ms 的含 I/O 总耗时这两个数字界定了当前实现向实时运行推进的程度，而保守子带掩蔽结果——仅保留约 50.3% 的子带-时间样本，却是防止射频干扰掩盖注入脉冲所必需的——是任何后续实现在真实数据上都必须重新审视的、本文特有的权衡。

本文引用的文献