---
title: "同一份立方体，怎样读成通道图、谱线与矩图"
description: "从保存的有限通道亮度定义矩，追踪beam混合与地图边缘的漏光，再用掩膜与相关噪声反例检验解释。"
canonical_url: "https://sandboxphysics.com/zh/stories/molecular-line-cube/moments-beams-and-a-saved-ppv-cube/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/molecular-line-cube/moments-beams-and-a-saved-ppv-cube/"
last_modified: "2026-11-04"
---

[分子谱线立方体：天空、速度与隐藏的深度](https://sandboxphysics.com/zh/stories/molecular-line-cube/)

分子谱线立方体 · L3

# 同一份立方体，怎样读成通道图、谱线与矩图

Sandbox Physics · 2026-11-04

通道图、单像素谱线和位置—速度图，都可以是同一份数据的不同切法。矩图再把速度轴压成少量数字。每压缩一次，就应问一次：保留了什么，又把哪种不同情况混在了一起？

## 开始前需要什么

需要透明线面积、光学厚度及有限通道平均的概念。下文用加权平均解释谱线矩，不要求统计学背景。

## 先把三个坐标写对

保存的PPV立方体有两个天空坐标和一个远离速度坐标。固定速度得到通道图，固定天空像素得到谱线；沿图像一行取出所有速度，得到位置—速度图。

交互记录包含$128\times128\times96$个亮度值。计算中的物理长度单位为0.1秒差距，完整天空窗口宽1秒差距，速度边界为$\pm4\,\mathrm{km\,s^{-1}}$。固定观察者位于正深度方向，正远离速度对应气体背离观察者。转动显示相机只换看云的角度，不重算这位观察者的数据。

物理场景中的颜色帮助辨认局部速度。沿深度重叠、吸收和beam混合之后，接收机实际得到的是数据区的亮度；两者不应被当成同一张表面贴图。

所有数据图、探针、有效掩膜和导出都读取最近一次完整接受的记录。新输入尚在计算时，旧记录仍是完整的一份；混合新通道图与旧谱线会制造根本不存在的物理不一致。

对应阅读：Sandbox Physics independent molecular-line science producer and reader package · Mangum & Shirley, How to Calculate Molecular Column Density, arXiv 1501.01703v4 (2017)

## 三个矩各自回答一个问题

用$T_i$表示某像素第$i$个通道的平均亮度，$v_i$是通道中心。零阶矩是线面积，一阶矩是以亮度为权的平均速度，二阶中心矩描述相对这个平均值的速度分散。分母接近零时，后两个量对微小误差特别敏感。

对称双峰的平均速度可能落在几乎没有发光的凹处，那里未必有独立的气体分量。宽度也可能来自多个分量、吸收展宽或beam内速度梯度，不能单凭它宣布温度或湍流能量。

一个可核算的离散例子：透明高斯标准差为$250\,\mathrm{m\,s^{-1}}$，通道宽为$83.333\,\mathrm{m\,s^{-1}}$。用区间平均亮度、却把每个通道的权重放到中心，细通道近似下增加的方差为$\Delta v^2/12\simeq578.7\,\mathrm{m^2\,s^{-2}}$，得到约$251.2\,\mathrm{m\,s^{-1}}$的离散宽度。这个近似要求谱线足够平滑、窗口足够宽；它不是所有窄线的通用修正式。

$$
W=\sum_iT_i\Delta v_i,\quad \bar v=\frac{\sum_i v_iT_i\Delta v_i}{W},\quad \sigma_v^2=\frac{\sum_i(v_i-\bar v)^2T_i\Delta v_i}{W}
$$

三者分别具有亮度乘速度、速度、速度平方的单位。这里使用的是保存通道的中心矩。

对应阅读：Sandbox Physics independent molecular-line science producer and reader package

## beam会把邻居带进一条谱线

有限beam，也就是望远镜的角响应，对天空方向作加权混合。若相邻像素的线峰位于不同速度，同一beam中的谱会叠加，单像素线便可能变宽。当地分子的微观运动保持不变，空间平均已经把谱线铺宽。

这里的云模型使用处处非负、对像素积分的高斯核，先把天空边界外的零值补足，完成卷积后再裁回保存窗口。完整卷积平面上的总量守恒；裁掉窗口外的光以后，保存地图中的总量则可以减少。

独立包保留过一个边缘漏光反例：裁下的地图漏掉约0.589%的积分量。把裁后地图重新归一到原总量会掩盖问题，且改变真实的表面亮度。交互域保留更严格的漏光拒绝门；被拒绝时继续显示上一份有效记录。

空间网格和速度网格需要分别细化。原来每轴64点的天空网格在部分注册病例中没有通过参照，128点版本才通过；把同一片云平移四分之一或半个像素，可以检查亮结构恰巧落在像素中心还是缝隙时，数值结果有多敏感。冷、薄、最大剪切与这些平移相位都有实际对照。整个连续控制域的统一误差定理仍待证明；现有证据覆盖的是这些具体计算。

对应阅读：Sandbox Physics independent molecular-line science producer and reader package · Bendo, Avison & Freestone, Measuring the Rotation of TW Hydrae (19 June 2019)

## 掩膜也会选择噪声

掩膜决定哪些像素和通道参与后续计算。假设一幅没有信号的噪声图围绕零对称。如果只保留高于三倍噪声标准差的像素，留下的噪声全部为正。再求和当然可能出现正的“线面积”。阈值选择改变了统计问题。

同样的阈值也会删掉真实弱线的翼部，线面积和宽度随之下降。独立反例同时保留这两种效应。Lab本身没有测量噪声，采用的是相对峰值的显示掩膜；它只负责隐藏极弱、矩不稳定的像素，不叫检测显著性。

真实ALMA教学cube的相邻通道经过Hanning相关处理。在选定空白环带和速度窗中估出的相邻相关系数约为0.668，两通道和的噪声标准差约是独立通道估计的1.29倍。这个数描述特定窗口的经验协方差，不是任意通道、任意掩膜的认证。

固定窗口的线性求和可以传播完整协方差；根据同一批噪声自适应挑选的掩膜还需要额外处理选择效应。实际练习因此同时保留带符号的固定窗积分和正阈值积分，让差异自己显现。

$$
\operatorname{Var}\!\left(\sum_iw_iT_i\right)=\sum_{ij}w_iC_{ij}w_j=\boldsymbol w^{\mathsf T}C\boldsymbol w
$$

固定权重下，非对角协方差也影响积分误差。独立噪声近似会删掉所有非对角项。

对应阅读：Sandbox Physics independent molecular-line science producer and reader package · Bendo, Avison & Freestone, Measuring the Rotation of TW Hydrae (19 June 2019)

## 练习与解释性答案

1. 左右两峰强度相同，速度分别为正负一公里每秒。一阶矩为零，能说明气体静止吗？

   展开解释

   不能。对称权重的平均恰为零，但两峰都不在零速度。应查看原谱，区分平均值与实际分量；光学厚度还可能进一步改变峰的解释。

2. 完整beam卷积守恒，为什么裁后的地图仍少了总量？

   展开解释

   卷积把部分信号送到了保存窗口以外。守恒成立于完整平面；裁剪改变了积分域。扩大窗口或把漏光单独记账可以检验它，不能靠裁后归一化把它抹掉。

## 动手检验预测

### [追踪一份数据的几种投影](https://sandboxphysics.com/zh/labs/molecular-line-cube/)

1. 在通道图中拖动视线，确认谱线标记与位置—速度图使用同一行。把通道移到线翼，比较局部亮度和积分图。
2. 切换平均速度图，寻找接近零的像素，再看原谱。解释这个零来自哪种对称性。
3. 把云移到允许的端点并改变转轴。读到的谱和所有图应一起更新；若状态被拒绝，检查上一份完整记录是否仍可导出。

[上一篇：气体怎样挡住自己的谱线](https://sandboxphysics.com/zh/stories/molecular-line-cube/opacity-and-the-planck-line-source/) [下一篇：用ALMA给TW Hya的盘做速度切片](https://sandboxphysics.com/zh/stories/molecular-line-cube/replay-alma-tw-hya-line-cube/)

## 来源与计算范围

1. [Mangum & Shirley, How to Calculate Molecular Column Density, arXiv 1501.01703v4 (2017)](https://arxiv.org/abs/1501.01703v4)

   固定v4：传输与亮度方程13—32、配分函数第7节、近似第11—14节及附录A—C。高斯宽度从归一化轮廓推导，另说明A4/A5的印刷问题；不据此声称求出了非LTE。

2. [Bendo, Avison & Freestone, Measuring the Rotation of TW Hydrae (19 June 2019)](https://www.alma.ac.uk/nas/gbendo/Web/Sci/Teach/experiment_twhyarotation.pdf)

   完整八页官方教学指南与具名原图成员。本地教程使用原生头中的频率、beam面积和处理历史；原文的历史距离与倾角不作当前精密值。

3. [Sandbox Physics independent molecular-line science producer and reader package](https://sandboxphysics.com/stories/molecular-line-cube/materials/molecular-line-cube-reader-science-v1.zip)

   自有合成输入、独立频率空间和连续射线参照，以及保留的分辨率失败。原生读者流程在本地取得具名源文件；包内没有观测图、原生表或App/Node引擎。