星等这个词我最早是从一张手抄的星图上认识的当时以为它就是个亮度排名数字越小越亮稀里糊涂记下来就完事。真正被它绊住是在郊外做一次目视观测按星图去找一颗标着 8 等的目标漫天星点死活找不到后来才发现自己拿目视极限当成了星图极限手里的设备根本到不了那个亮度。再往后接触 CCD 测光误差一度飘到 0.3 等以上回头刨根问底问题全出在对星等这套体系的理解上——它是怎么定出来的、在不同波段下代表什么、定标链路里哪一步没做对。这篇东西就是把这些年关于星等的经验和踩过的坑整理一遍。从两千多年前靠肉眼排出来的六档亮度到今天用像素阵列把信号量化到毫星等中间那套对数规则其实一直没变变的只是测量手段和误差控制。内容既适合刚入门、想把几等星这个概念彻底搞清楚的朋友也适合已经架起望远镜、开始做简易测光、但结果老是和别人对不上的同好。我会把公式、参数、实操步骤全写出来能直接抄的部分直接抄但也会标清楚哪些地方是我自己实测出来的经验值别照搬。1. 星等这套古老体系为什么今天还在用1.1 从一等星到六等星说起一个反直觉的刻度最早那套分级完全是凭眼睛来的。两千多年前的一批天文学家把天上能看到的星按亮度分成六档最亮的叫一等最暗的叫六等。这个划分没有任何仪器参与纯靠人眼对比排序所以它本质上是一个感知排序不是物理量测量。人眼对亮度的响应本来就接近对数——亮度翻一倍你感觉到的不是翻一倍而是增加了一点点。这也是为什么这套粗糙的分级居然用到了今天。问题在于肉眼排序的结果是离散的一共就六档相邻两档之间实际差多少亮度根本说不清。到了十九世纪中期一位英国天文学家做了件关键的事他拿仪器实际测了一批星的亮度发现传统分级的一等和六等之间亮度差大约是一百倍。于是他提了个假设——既然一等到六等差 5 个星等、亮度差 100 倍那就让星等和亮度的关系遵循一个固定的对数规则。这个假设后来成了整个星等体系的基石也是今天所有测光工作的起点。这里有个必须说清的反直觉点**星等是越小越亮而且是负得越多越亮。**太阳的视星等约 -26.7满月约 -12.7全天最亮的那颗恒星约 -1.46。初学者最容易犯的错就是把它当分数理解觉得 -1.46 比 6.5 小所以暗正好反了。我在带新人的时候一般会让他们先背一句话数字越小越亮负号代表特别亮。1.2 对数刻度的数学内核与那几个必须记住的换算星等和亮度的关系式写出来很简单m₁ - m₂ -2.5 × lg(E₁ / E₂)m 是星等E 是能流密度或者说亮度、光通量具体用哪个量取决于你测的是什么。这个 -2.5 系数不是随便定的它就是5 个星等对应 100 倍亮度这个约定反推出来的——因为 lg(100) 2要让 5 2.5 × 2 成立系数就得是 2.5前面再加个负号保证越亮星等越小。由此可以直接推出几个必须刻在脑子里的换算关系星等差亮度比直观理解0.1 等约 1.096 倍约差 10%测光精度的常见门槛1 等约 2.512 倍相邻两档的亮度关系2.5 等10 倍换算时最常用的跳板5 等100 倍体系定义的源头10 等10000 倍大口径和小口径的差距量级反过来如果知道亮度比要算星等差就用 m₁ - m₂ -2.5 × lg(E₁/E₂)。这个式子在写脚本的时候天天用我一般会封成两个函数放工具箱里import math def mag_diff_to_flux_ratio(dm): 已知星等差求亮度比 E1/E2 return 10 ** (-0.4 * dm) def flux_ratio_to_mag_diff(ratio): 已知亮度比 E1/E2求星等差 m1-m2 return -2.5 * math.log10(ratio) # 举例13 等星比 0 等星暗多少倍 print(mag_diff_to_flux_ratio(13)) # 约 6.31e-06 print(flux_ratio_to_mag_diff(1e6)) # 约 -15.0实测下来新手最容易在这里出的岔子不是公式记错而是符号方向搞反。比如目标比参考星暗 1 等结果算出来写成 -1 等整条光变曲线就翻过来了。我的习惯是先在纸上把谁减谁写清楚再往代码里填宁可多写一行注释。注意lg 是常用对数以 10 为底不少文献里写成 log别和自然对数 ln 混用。写脚本时 Python 的 math.log 默认是自然对数必须用 math.log10。1.3 视星等、绝对星等、表面亮度三个别搞混的量这三个概念混起来用是重灾区我自己也被坑过。视星等m从地球上看到的亮度取决于天体本身的光度和它离我们多远。它是一个观测值跟你在哪、用什么设备、测哪个波段都有关系。所有人说的这颗星几等默认都是视星等。绝对星等M把天体强行挪到 10 秒差距约 32.6 光年处它该有的视星等。这是天体本身光度的直接体现和距离无关。换算用距离模数公式m - M 5 × lg(d) - 5其中 d 的单位是秒差距pc把常见的距离代进去可以做成一张速查表距离距离模数 m-M说明10 pc0定义点100 pc5银河系内近邻1 kpc10银河系内100 kpc20星系晕、近邻星系1 Mpc25本星系群附近10 Mpc30需要大口径反过来如果已知视星等和绝对星等可以反推距离d 10^((m-M5)/5)单位 pc。这个式子在实际工作里用得极多我在处理变量星样本的时候基本每次都要跑一遍。表面亮度这个是最容易被忽略的。点状天体恒星讲星等没问题但星系、星云这类延展天体光被摊在一大片面积上你报一个总星等其实意义有限。这时候要用表面亮度单位是星等每平方角秒mag/arcsec²。举个实测感受一个总星等 9 等的星系如果它的视面积是 100 平方角秒那平均表面亮度大约是 9 2.5 × lg(100) 14 mag/arcsec²但如果摊到 10000 平方角秒就是 19 mag/arcsec²接近城市夜空背景肉眼几乎看不到。心得判断一个延展天体能不能看到看总星等是没用的得看表面亮度和当地夜空背景亮度的对比。我用一台 200mm 的镜子在城郊看一个总星等 9.5 等的星系什么都看不到换到暗夜环境下同样的设备同样的目标就出来了。差别全在背景。2. 现代测光的完整链路从光子进筒到数字出炉2.1 测光系统与滤光片为什么不能只说这颗星几等一个裸的几等是不够用的。因为天体辐射是连续谱不同波段能量分布不同你用不同颜色的滤光片去测得到的星等数值会差很多。一颗偏红的恒星和一颗偏蓝的恒星在蓝波段和红波段的星等差可以有好几等。所以正规的做法是报星等必须带上波段标记比如 V 11.24意思是 V 波段测出来 11.24 等。二十世纪中叶天文界建立了一套经典的三色测光系统用 U、B、V 三个字母命名分别对应紫外、蓝、可见光偏黄绿三个宽波段。它的设计思路其实很朴素把人眼敏感的波段定为 V再往两边补两个形成三个互相独立的通道。后来又扩展出 R、I 等更红的波段。这套系统一直用到现在原因是它的滤光片响应曲线被定义得非常明确全世界的观测数据能互相比较。再往后出现了更多套系统比如某大型数字巡天项目使用的五个波段u、g、r、i、z特点是滤光片更宽、覆盖波长范围更连续好处是能一次拿到天体谱能量分布的粗轮廓代价是波段之间重叠较多颜色项处理更麻烦。系统常用波段特点适合场景经典三色系统U、B、V带宽窄、定义明确、历史数据多变星、标准星定标红端扩展R、I适合偏红目标红矮星、红化严重的星场宽波段五色u、g、r、i、z覆盖连续、流量定标方便巡天、星系测光窄带Hα、OIII 等只挑特定谱线星云、发射线天体实际选滤光片的时候我一般遵循两条一是要和已有数据兼容你测出来的数得能跟别人的对得上否则毫无意义二是要考虑目标本身的颜色测一颗极红的星用 U 波段信号弱得可怜曝光时间翻三倍都不一定够。2.2 定标体系Vega 系统、AB 系统与流量单位星等是个相对量必须有个零点。历史上最自然的做法是拿一颗亮度稳定的星当0 等标准这就是所谓的 Vega 系统——那颗天琴座最亮的 0 等星在 V 波段被人为定义为 0 等左右实际上约 0.03 等。这套系统直观、好记但有个致命问题那颗星本身在红外波段能量偏低导致红外波段的零点非常别扭。于是有了另一套基于物理流量的系统叫 AB 系统。它的核心思想是不管什么波段零点都定成一个固定流量 3631 央斯基Jy。这样一来各波段的星等能直接反映物理流量做谱能量分布拟合特别方便。换算公式m_AB -2.5 × lg(f_ν) - 48.60f_ν 的单位是 erg·s⁻¹·cm⁻²·Hz⁻¹。验证一下f_ν 3631 Jy 3.631 × 10⁻²⁰ erg·s⁻¹·cm⁻²·Hz⁻¹-2.5 × lg(3.631e-20) 48.60所以 m_AB 0对上了。这两套系统的数值在 V 波段附近差别很小零点差零点零几等但越往红外差别越大。我踩过的坑是拿了一个星表的 AB 星等去和自己算的 Vega 系统零点比发现偏了 0.3 等以上还以为是测量错了查了半天才发现是系统不同。注意拿到任何一组星表数据第一件事是查清楚它是 Vega 系统还是 AB 系统。这个信息一般在星表的说明文档里不看清楚后面全白干。2.3 大气消光与颜色项不做这两步校正数据基本废了这两个校正是业余测光和能出学术数据之间最大的分水岭很多人做到测光就停了数据看着漂亮但和别人永远对不上问题就在这。大气消光天体光穿过大气层会被吸收和散射越靠近地平线穿过的空气越厚损失越多。空气厚度用大气质量airmass表示天顶方向约等于 1高度角 30 度时约等于 2越接近地平线涨得越快粗略用 X sec z 估算z 是天顶距但这个近似在低高度角会失效需要更复杂的公式。消光量 消光系数 k × 空气质量 X。V 波段的消光系数在海平面好天气下典型值约 0.15 mag/airmass高海拔可以低到 0.10污染重的地方能到 0.3 以上。B 波段更大约 0.25U 波段能到 0.4 以上。这意味着一颗在天顶测出来 12.00 等的星挪到高度角 30 度X2再测就变成 12.30 等了。差 0.3 等已经超过大多数变星研究的精度要求。颜色项不同滤光片实际响应曲线和标准系统并不完全一致加上大气消光本身也和颜色有关蓝光被散射得更厉害所以需要按目标的颜色指数做二次修正。完整定标式通常写成m_std m_inst ZP k × X C × (颜色指数)m_std 是标准星等m_inst 是仪器星等ZP 是零点k 是消光系数X 是空气质量C 是颜色项系数。C 的典型值在 0.03 到 0.10 之间看起来很小但对于毫星等精度的研究就是主要误差源。实操上求这些系数得观测一组标准星覆盖不同的高度角和不同的颜色然后做多元线性回归。我一般会选 6 到 10 颗标准星高度角从 30 度到接近天顶铺开颜色指数跨度尽量大。一个晚上跑下来能把 ZP、k、C 三个系数同时解出来。3. 用小型望远镜做一次靠谱的测光实战3.1 器材搭配与场地选择不是非得大口径才能做测光我见过用 80mm 折射镜做得非常漂亮的变星光变曲线。关键是搭配合理、流程规范。部件建议配置选型理由主镜150–250mm 反射或折射焦比 f/4–f/6焦比太快像质难保证太慢效率低赤道仪带自动导星的赤道仪承重要留 30% 余量长曝光必须稳风一吹就废相机制冷型单色传感器像元 3–5 μm单色量子效率高制冷压低暗流滤光片至少有 V 波段最好配齐 B、R兼容历史数据能做颜色项导星独立导星镜或离轴导星焦距 1000mm 以上必须导星场地无月夜夜空背景优于 20 mag/arcsec²背景亮一倍信噪比掉一半关于场地我想多说两句。很多人纠结设备其实夜空背景才是决定你能测多暗的第一因素。同一个设备在背景 21.5 mag/arcsec² 的地方和 18 mag/arcsec² 的地方极限星等能差 2 到 3 等相当于口径差了两三倍。我自己的经验是宁可开车两小时去暗一点的地方也别花钱升级设备。3.2 拍摄参数怎么定从信噪比公式倒推曝光时间这是整个流程里最工程师的一步也是最能体现专业度的地方。信噪比的经典公式CCD 方程是SNR S / sqrt(S n_pix × (B D RN²))S 是目标总电子数n_pix 是参与计算的像素数含背景环B 是每像素背景电子数D 是暗流RN 是读噪声。这个公式的意思是噪声有三个来源——光子本身的散粒噪声sqrt(S)、背景和暗流的散粒噪声、读噪声三者平方相加。实际走一遍。假设一台 200mm 口径、焦距 1000mm、焦比 f/5 的反射镜配一台像元 3.76 μm、读噪声 3 e⁻、制冷后暗流可忽略的单色相机。第一步算像素尺度像素尺度 206265 × 像元尺寸(mm) / 焦距(mm) 206265 × 0.00376 / 1000 ≈ 0.78 角秒/像素这个值和典型视宁度2–3 角秒比偏小略微过采样但可以接受。如果想更匹配视宁度可以用 2×2 合并代价是损失分辨率。第二步估算零等星在入瞳处的光子数。V 波段 0 等星在地面每秒每平方米约 10¹⁰ 个光子这个数记住了很好用。口径 0.2 m有效面积约 0.0314 m²光学效率取 0.7量子效率取 0.8综合下来10¹⁰ × 0.0314 × 0.7 × 0.8 ≈ 1.76 × 10⁸ e⁻/s。第三步算目标电子数。一颗 13 等星流量是 0 等星的 10^(-0.4×13) ≈ 6.3 × 10⁻⁶ 倍所以 1.76e8 × 6.3e-6 ≈ 1110 e⁻/s。曝光 120 秒得到约 1.33 × 10⁵ 个电子。第四步算背景。假设夜空背景 20.5 mag/arcsec²像素面积 0.78² ≈ 0.61 平方角秒所以每像素对应 20.5 2.5 × lg(0.61) ≈ 21.0 等。相对 0 等星的流量比是 10^(-0.4×21.0) ≈ 1.0 × 10⁻⁸·⁴换算成电子率约 0.6 e⁻/s/px。120 秒下来 72 e⁻/px。取 n_pix 60孔径加背景环的有效像素数背景噪声 sqrt(60 × 72) ≈ 66 e⁻。第五步算总信噪比。散粒噪声 sqrt(1.33e5) ≈ 365。总噪声 sqrt(365² 66² 60×3²) ≈ sqrt(133225 4356 540) ≈ 372。SNR ≈ 133000 / 372 ≈ 357。远超需求。这说明什么说明 200mm 加 CCD 在 120 秒下测 13 等星绰绰有余。那极限在哪按 SNR 5 反解能算出这个配置在这套假设下的理论极限在 19 到 20 等附近。但实测一般会低 1 到 2 等因为实际光学效率、视宁度、云、月光、平场误差都会咬掉一截。我自己的经验是理论算出来 19.5实际能稳定拿到 17.5 就已经不错了。import math def estimate_snr(mag_target, mag_zero_count_rate, t_exp, sky_mag_per_arcsec2, pixel_scale, n_pix, read_noise3.0, dark_current0.0): mag_zero_count_rate: 0 等星每秒总电子数含光学效率与QE flux_ratio 10 ** (-0.4 * mag_target) S mag_zero_count_rate * flux_ratio * t_exp pix_area pixel_scale ** 2 sky_mag_per_pix sky_mag_per_arcsec2 - 2.5 * math.log10(pix_area) sky_ratio 10 ** (-0.4 * sky_mag_per_pix) B mag_zero_count_rate * sky_ratio * t_exp D dark_current * t_exp noise math.sqrt(S n_pix * (B D read_noise ** 2)) return S / noise z 1.76e8 print(round(estimate_snr(13, z, 120, 20.5, 0.78, 60), 1)) # 约 357提示理论计算只是给个量级别拿它当准绳。真正靠谱的做法是同一目标拍 3 组不同曝光时间看信噪比是不是按 sqrt(t) 增长不增长说明有系统性噪声比如未制冷的暗流、背景梯度先把系统问题解决掉。3.3 图像预处理本底、暗流、平场的正确顺序这一步最容易出错顺序错了后面全废。正确的处理链条是本底帧Bias零秒曝光记录读出电路的固定偏置。需要拍 20 到 50 张中值合成。制冷相机的本底非常稳定拍一次能用很久但换了增益或温度就得重拍。暗流帧Dark与科学帧相同曝光时间、相同温度、快门关闭。同样 20 张以上中值合成。制冷到 -10℃ 以下暗流通常小到可以忽略可以跳过但要先确认。平场帧Flat这一步最容易被敷衍。平场的作用是消除光学系统的渐晕、灰尘阴影、像素响应不均。必须用均匀光源拍而且要在和科学帧相同的对焦、相同的滤镜、相同的相机角度下拍。黄昏天光平场是业余条件下最实用的方案因为天顶方向的天空亮度足够均匀。平场拍不好结果就是同一颗星在视场中心测出 12.00 等挪到边缘测出 12.15 等你还以为自己测错了。处理后合成主平场Master Flat。标准流程是每个科学帧先减主本底再减主暗流然后除以归一化后的主平场。归一化就是把平场除以它自己的中位值让它平均值约等于 1这样除法不会改变整体流量标度。注意平场千万别用饱和的帧。天光平场的典型曝光要控制在中位计数 30% 到 50% 满阱之间太暗噪声大太亮会饱和。我一般会拍 15 到 20 张不同曝光时间的事后挑中间那批。3.4 孔径测光与零点求解预处理完的图丢进测光软件选测光方式。两种主流孔径测光在星像位置画一个圆孔径把圆内像素加起来当信号再在圆外画一个同心环取环内中位数当背景乘以孔径面积后减去。优点是简单快、对暗星友好缺点是当星像密集时会互相污染。点扩散函数拟合测光先用若干亮星拟合出一个二维点扩散函数模型再把模型套到每颗星上。优点是能处理拥挤星场、精度更高缺点是需要足够多的亮星、计算量大、参数调不好容易翻车。对业余条件我建议先用孔径测光把流程跑通再考虑升级。孔径大小的选择有个经验值孔径半径取点扩散函数半高全宽的 1.5 到 2 倍。如果视宁度是 3 角秒、像素尺度 0.78 角秒/像素那半高全宽约 3.8 像素孔径半径取 6 到 8 像素比较合适。背景环的内径取孔径半径加 3 到 5 像素外径取内径的 2 倍左右。仪器星等算出来是m_inst -2.5 × lg(计数率) 常数。这个常数是任意的真正的标准星等要靠零点定标ZP m_std - m_inst经过消光校正后系综零点比单星零点靠谱得多。做法是视场里如果有 10 颗星表里的已知星分别算它们的零点取中位数或加权平均得到这一帧的零点。这样能抵消单颗星表星本身的误差。3.5 结果校验拿已知星和星表对表测完不算完必须校验。我一般会做三件事第一残差检查。把视场里所有星表已知星的测光结果和星表值相减看残差分布。如果残差是随机分布的、标准差在 0.02 到 0.05 等之间说明流程正常。如果残差随星等呈现趋势亮星偏正、暗星偏负说明有非线性或者背景扣除有问题。第二空间分布检查。把残差按视场位置画出来如果边缘系统偏负八成是平场没做好。这一条特别灵我几次出问题都是靠这个发现的。第三重复性检查。同一目标同一晚上连拍三组看结果的一致性。组间标准差明显大于单组内部误差说明有系统性问题比如云飘过、跟踪抖动、焦点漂移。4. 极限星等肉眼、双筒、望远镜、CCD 的实战数据4.1 肉眼极限星等与夜空亮度分级肉眼能看到多暗的星取决于瞳孔直径、视觉适应程度、天气透明度和夜空背景亮度。无月条件下的理想值大约 6.5 等极暗的地方能看到 7 等甚至 7.5 等城市里普遍在 3 到 4 等。业余观测里常用一套从 1 到 9 级的暗夜分级来描述夜空质量1 级是最暗的原始天空9 级是市中心。我做了个对照表分级环境描述典型肉眼极限背景亮度mag/arcsec²1无人区无任何人工光7.522.0 以上3乡村边缘6.821.55城郊6.020.57城市边缘4.5–5.019.09市中心3.0–4.017.0 以下这套分级不是拿来装专业的它特别实用你只要报出自己场地的分级别人就能大致判断你的数据质量上限。我给自己常去的两个场地分别做了测定一个 4 级一个 6 级同一套设备在这两处的 CCD 极限星等差了将近 2 等。4.2 望远镜目视极限星等的经验公式目视望远镜的极限星等有个流传很广的经验公式m_lim ≈ 2.1 5 × lg(D)D 是口径单位毫米。代入一下100mm 给 12.1 等200mm 给 13.6 等300mm 给 14.5 等。这个公式在中口径范围内比较准大口径会低估。另一个更符合物理的算法是光收集增益法望远镜口径 D 相对于人眼瞳孔 d暗适应后约 7mm光收集面积大了 (D/d)² 倍星等增益是 2.5 × lg((D/d)²) 5 × lg(D/d)。200mm 对 7mm 瞳孔增益 5 × lg(28.6) ≈ 7.3 等加上肉眼极限 6.5 等得到 13.8 等。和上面那个公式基本一致。但这里有个巨大的坑这两个公式都假设理想条件。实际使用中光学镀膜损失、副镜遮挡、视宁度、光污染都会打折扣。我在 5 级场地上用 200mm 实际能看到的极限大概 12.5 到 13 等比公式低了 0.5 到 1 等。设备理论极限4 级场地实测7 级场地实测肉眼6.56.04.510×50 双筒10.09.57.5100mm 折射12.111.59.5200mm 反射13.813.011.0300mm 反射14.914.012.04.3 CCD 极限星等的估算与实测偏差CCD 的极限星等和目视完全是两套逻辑。目视是人眼能否感知到差别CCD 是信号能否超过噪声到统计显著的程度通常以 SNR 5 为探测门槛。理论上它比目视深得多但也更容易被各种系统性因素拖累。前面 3.2 节算过200mm CCD 在 20.5 mag/arcsec² 背景下 120 秒的理论 5σ 极限约 19.5 等。但这是单帧理论值实际还有几个折损视宁度星像摊开后信号分散到更多像素同一孔径内有效信号减少。视宁度从 2 角秒恶化到 4 角秒极限可能掉 0.5 等。背景梯度月光、城市光造成大尺度背景亮度不均孔径测光的背景估计会偏。平场残余平场不完美会引入 1% 到 3% 的误差对暗星来说是致命的因为暗星信号本来就弱。叠加的方式单帧极限是一回事多帧叠加后极限能提升。叠加 N 帧理论极限提升 2.5 × lg(sqrt(N)) ≈ 1.25 × lg(N)。叠加 10 帧提升 1.25 等叠加 100 帧提升 2.5 等。但前提是每帧都对得准、背景扣除干净。我自己实测的数据200mm f/5 制冷相机20.5 背景单帧 180 秒SNR5 的极限约 17.8 等叠加 20 帧后能到 19.3 等左右。理论算出来是 19.5 和 21.1实际分别低了 1.7 和 1.8 等。这个差距我反复验证过主要来自视宁度和平场残余。心得看到文献上吹的极限星等别太上头。那些数据通常是在顶级台址、顶级设备、超长叠加、精心处理的条件下得到的。业余条件下老老实实按实测数据规划观测目标能省下大量无用功。5. 常见问题排查与踩坑记录5.1 测光误差来源速查做了这么久我把遇到的误差源整理成了一张表出问题时从上往下捋基本能定位到。误差源典型量级表现特征排查方法光子散粒噪声0.001–0.01 等随机、不可消除增加曝光或叠加帧数天空背景噪声0.005–0.05 等随机与背景亮度正相关换暗夜场地、缩小孔径平场残余0.01–0.05 等随视场位置系统性变化检查平场是否同焦点同滤镜大气消光未校正0.05–0.5 等与高度角强相关加消光系数校正颜色项未校正0.01–0.1 等与目标颜色相关用多颗不同颜色标准星回归零点选错0.05–0.3 等整体偏移核对星表的星等系统饱和0.1 等以上亮星偏暗、平坦检查峰值计数是否接近满阱宇宙线击中随机跳变单帧异常亮点多帧中值叠加自动剔除云0.05–1 等整体或局部变暗无规律看导星曲线、看背景计数用这张表的时候有个技巧先看误差是不是随空间变化。随空间变八成是平场或光学问题随高度角变是消光随颜色变是颜色项随机的是噪声或天气。5.2 几个让人抓狂的真实场景场景一同一颗星前后两晚测出来差 0.15 等。我一开始怀疑是设备问题后来发现那两晚的透明度不同。解决办法是每晚上都重新测零点不要复用前一晚的零点。更严谨的做法是观测标准星场把当晚的消光系数也解出来。这个坑我踩过两次后来养成习惯——一晚上至少拍三组标准星分布在不同的高度角。场景二视场中心的星测得准边缘的星系统偏暗 0.1 等。经典平场问题。具体原因可能是平场拍的时候对焦点变了或者平场光源不均匀或者拍摄过程中相机角度动过。检查办法是把残差按视场坐标画散点图如果呈现明显的径向或者环形结构基本可以确诊。场景三暗星测出来的星等太好比星表还亮。八成是背景扣除过度或者是附近有一颗暗星污染了孔径。前者是因为背景环里混进了别的星中位数偏高减多了后者要靠检查星图。我的做法是对每颗测光星先看它周围 3 倍孔径半径内有没有别的源有就标记出来单独处理。场景四叠加之后极限星等没提升。常见原因是帧与帧之间的背景水平不一致叠加时用了平均值而不是中值或者叠加前没做背景归一化。还有一个更隐蔽的原因跟踪有系统性漂移星像位置每帧都在动叠加时被抹平了。这种情况得先做星点对齐再叠加。5.3 一套可复用的自检清单每次出数据前我会按这个清单走一遍基本上能把 90% 的问题挡在发布之前。本底、暗流、平场是否都是同一温度、同一增益、同一滤镜下拍的。平场的中位计数是否在满阱的 30% 到 50% 之间有没有饱和。科学帧的峰值计数是否超过满阱的 70%超了就缩短曝光重拍。每帧的背景中位值是否稳定波动超过 10% 要查云。星像半高全宽是否稳定突然变大查对焦和视宁度。定标用的标准星是否至少有 5 颗残差标准差是否小于 0.05 等。残差有没有随位置、星等、颜色的系统性趋势。是否记录了观测时刻、空气质量、天气、月相。第 7 条特别关键。残差无结构说明流程干净残差有结构说明还有系统误差没消掉。这句话我贴在显示器边上。6. 星等的延伸应用从变量星到系外行星6.1 光变曲线星等随时间变化的价值星等本身是个静态量但一旦加上时间轴它能讲的故事就多了。变量星的亮度变化周期和振幅能直接反映它的物理性质比如脉动变星的周期和光度之间存在确定关系可以当标准烛光用来测距。这套方法在二十世纪初奠定了宇宙距离尺度的基础到今天还在用。业余条件下做变星测光门槛其实不高。目标选 8 到 13 等的脉动变星用 150mm 以上的镜子配 V 波段滤光片每晚拍几十帧处理完画出光变曲线能看到明显的周期性起伏。我第一次做出完整光变曲线的时候那种成就感比拍一张漂亮的天体照片强多了——因为这是真正的数据不是图像。做光变曲线有个关键要求同一颗变星前后几次观测必须用同一套参照星和同一套流程。因为星等是相对的参照星一换整条曲线的零点就飘了。我一般会在视场里固定选 3 到 5 颗参照星整个观测季都不换并且定期检查它们的星等是否稳定。6.2 凌星测光毫星等精度的挑战当一颗行星从恒星前面经过恒星的亮度会短暂下降下降幅度通常在 0.01 等即 1%甚至更小。要测出这种变化精度必须做到毫星等这是业余测光的天花板级挑战。难度在哪1% 的亮度变化意味着信噪比必须远高于 100而且所有系统误差都要压到 0.001 等以下。具体来说平场必须极其精准任何 0.5% 的残余都会直接翻译成 0.005 等的假信号。不能有饱和目标星必须落在传感器的线性区间内。必须做差分测光也就是同时测目标星和若干颗亮度相近的参照星然后算差值。这样能抵消天气、消光、高度角变化带来的共模误差。这是业余凌星测光能成功的最关键技巧。时间精度要高曝光中点时刻要精确到秒级否则光变曲线的形状会变形。我自己尝试过几次最好的一次测到了一个约 0.02 等的下降曲线形状能看出来但噪声还是偏大。后来复盘主要问题是平场不够精细、没有用制冷相机导致暗流不稳。这件事让我彻底理解了那句话毫星等测光拼的不是望远镜是误差控制。6.3 光污染测量星等也是环境指标这个应用可能很多人没想过。星等还能用来量化光污染在同一地点、同一时间数一数天顶附近某个天区内能看到的星数和最暗的那颗的星等就能反推出当地的夜空亮度。具体做法很朴素选一个固定的天区比如某个星座围成的区域数清里面肉眼可见的星查出它们的星等最暗的那颗的星等就是当天的肉眼极限星等。连续记录几个月能看出明显的变化趋势。我用这个方法在自己常去的场地做了两年的记录肉眼极限从 6.0 等降到了 5.6 等对应的背景亮度上升了约 0.4 mag/arcsec²肉眼能明显感觉星星变少了。如果手头有相机还可以更精确用广角镜头对着天顶拍固定曝光测出天光背景的每平方角秒星等再换算成物理亮度单位。这套方法的好处是不依赖人眼、可重复、可对比比我觉得今天天挺黑这种描述可靠得多。最后分享一个我自己一直在用的小习惯每次观测结束不管多晚多累我都会在记录本上手写三行——今晚的零点、今晚的消光系数、今晚觉得最可疑的一个环节。这三行字看起来没什么但过几个月回头看往往就是这些零碎的记录帮我把当时想不通的问题想通了。星等这套东西说到底是个人造的对数刻度它本身不带任何物理意义但正是因为它足够简单、足够稳定、全世界通用才撑起了整个天文学最基础的测量体系。把它吃透后面所有的测光和数据分析都会顺很多。
