固定具体产品，不只记住档案首页
来源为 UVIS_HSP_2008_261_ALPARA_E_TAU01KM.TAB，是 2008 年 9 月 17 日出掩星观测，径向采样一千米，波段一百一十至一百九十纳米。产品标签和整卷说明都标为 CO-SR-UVIS-HSP-2/4-OCC-V3.0，数据集发表于 2020 年，产品创建于 11 月 27 日。UVIS 通用网页仍在介绍第二版，这里以实际下载标签为准。
固定表格保留十二个原始列，从十二万五千到十二万六千千米，共 1001 条逐字节来源记录，每条一百二十一字节。SHA-256 为 6baadf1d91701ad8e868745f7aad9e49d853cec169ce45fea51b9a4f30fdd94c。原标签、原说明和出处可一同下载。下载时只提取这段半径所需的 HTTP 字节范围，整文件哈希尚未验证，输入本身也已是档案分箱产品。
保留的质量标志全部为零。试选的其他观测带有双星或未建模背景标志，因此没有用于本例。换一个产品时，不要因为曲线好看就删掉非零标志，应保留原值，并明确分析规则。
对应阅读：Cassini UVIS alpha Arae2008-261 egress, original PDS3 product label · NASA PDS Cassini UVIS archive README, version3.0 · Colwell et al., Cassini UVIS stellar occultation observations of Saturn’s rings
取对数前先认清单位
档案的平均信号，是每次原始积分的光子计数在径向箱内取均值；这一列尚未扣背景。恒星模型给出未遮挡恒星信号，背景模型另列。从平均信号扣一次背景，除以恒星模型，再取对数，乘上绝对环面仰角 54.4 度的正弦。
已发布法向光学厚度列已经经过档案处理，这里用来独立对照。在它上面再扣背景没有物理含义。重算值与打印值的最大差约 0.00171；我们使用的是舍入后的信号、模型和仰角，没有原始未舍入几何，这种细小差异因此应保留下来。
半径是环面内到土星中心的距离；经度与环面视线角以度计。两列时间为自 J2000 起算的 TDB 秒，光子先离开环面，后到 Cassini。以后若要接上卫星图样相位，这个差别不能丢。光的波长固定在紫外波段，本次搜索的波长则是沿半径量出的空间间距。
对应阅读：Cassini UVIS alpha Arae2008-261 egress, original PDS3 product label · Colwell et al., Cassini UVIS stellar occultation observations of Saturn’s rings
带上计数误差，别虚构标定精度
样本数列说明每个径向箱平均了多少次原始积分。若各次为独立泊松计数，均值方差就是平均计数除以样本数。再按下式，把它传播过背景扣除与对数。
这个条件估计将档案背景与未遮挡恒星模型固定，未包括它们的误差、恒星有限大小、衍射、残余仪器效应及径向相关。标签里的最大可探测厚度是指定计数统计灵敏度，不是逐行标准差，不能拿来替换误差条。
没有样本、有坏标志或净信号非正的箱，需要明确删失或采用似然处理。脚本核对保留形状与标志，推导函数保留无效状态，不会把删失厚度换成一个方便的高值。
对应阅读：Cassini UVIS alpha Arae2008-261 egress, original PDS3 product label · Colwell et al., Cassini UVIS stellar occultation observations of Saturn’s rings
完整跑一次固定分析
将 Python 脚本、精确表格、标签、说明与出处放在一起。命令核对表格哈希，重算光学柱，保存全部实测数组与标志，再完成声明的局部搜索和敏感性比较。输出为 cassini-replay-rerun.json，可与保留的 cassini-replay.json 对照。
公开的小型科学重跑与受保护的交互引擎分开。可选的另一条命令，独立重分析保存的合成粒子位置与计数。两者都不下载新观测，不导入实验引擎，也不把模拟冰粒变成 Cassini 实测粒子。
python3 -W error replay-cassini.py
python3 -W error replay-synthetic.py
| 窗口起点 / 千米 | 一千米分箱尺度 / 千米 | 两千米分箱尺度 / 千米 | 约化卡方 |
|---|---|---|---|
对应阅读：Cassini UVIS alpha Arae2008-261 egress, original PDS3 product label · NASA PDS Cassini UVIS archive README, version3.0
量一个局部尺度，再看残差
显示的一千千米片段包含 Janus 四比三共振区域附近的强结构。预先规定三个八十千米窗口，分别从 125310、125410、125510 千米开始。每次用仅含计数误差的权重，同时拟合正弦、余弦与二次径向基线。空间周期范围为三至六十千米，频率均匀网格有 1200 点，实际数组全部保存。
主导分量分别为 16.53、5.943、5.172 千米。另一个三箱平滑波峰识别器得到不整齐的相邻间距，包括二十六、二十、十六和九千米；哪些峰通过阈值与最小间距规则，会影响结果。这些测量支持径向结构总体变细，却没有证明每个峰都属于同一列线性波。
关键诊断是约化卡方，三个窗口约为一百零四、一百零七和五十四，计数噪声解释不了残差。非线性波峰和重叠结构仍在。加入 256 组独立计数噪声扰动后，前两个窗口的敏感性范围很窄，但它们不是物理波长置信区间，最后一个窗口已经会在竞争分量之间切换。
对应阅读：Colwell et al., Cassini UVIS stellar occultation observations of Saturn’s rings · Salo & Mondino-Llermanos, Photometric modelling of self-gravity wakes and overstable oscillations · Goldreich & Tremaine, The excitation and evolution of density waves
用实际计数改变分辨率
将相邻径向箱按原始积分次数合并，用同样权重平均源与背景模型，累加积分次数，再从合并强度重算光学厚度。直接平均已经取过对数的厚度，会得到另一种统计量。
两千米分箱的搜索，将最短波长限制为不短于实际径向间距中位数的两倍。前两个主导分量仍约 16.49 与 5.916 千米，最后一个却变为约 4.006 千米。如此依赖分辨率与局部模型的结果，不能只附一个很窄的计数噪声范围就报告。
重跑得到的是实测光学剖面与指定的局部空间尺度。研究级分析还需分离合适的波，建模土星引力与共轨卫星随时间变化的强迫，并先验证非线性结构处理。
对应阅读：Cassini UVIS alpha Arae2008-261 egress, original PDS3 product label · Colwell et al., Cassini UVIS stellar occultation observations of Saturn’s rings · Goldreich & Tremaine, The excitation and evolution of density waves
让结果带着自己的来源历史
引用 Colwell、Jerousek、Becker 与 Esposito 的 2020 年 PDS 第三版数据集，以及 2010 年 UVIS 分析论文。原始说明保留了引用请求。这是可公开访问的 NASA PDS 科学数据，教程没有替作者声明另一份 Creative Commons 授权。
换一次掩星时，先固定精确产品和标签，让计数归一、几何与标志跟着每一步走。再换一个窗口或分箱宽度比较空间尺度。一个有用的结果，既包括恢复出来的结构，也包括提醒你停止解释的那次拟合。
对应阅读：Cassini UVIS alpha Arae2008-261 egress, original PDS3 product label · NASA PDS Cassini UVIS archive README, version3.0 · Colwell et al., Cassini UVIS stellar occultation observations of Saturn’s rings