在啃《Introduction to Infrared and Electro-Optical Systems》第二章后半部分之前,我一直把“辐射度学”当成一堆单位换算表来背,直到做完这一章的课后习题、看了几组实测数据,才意识到自己之前对红外系统的理解完全是“只知其表”。这一篇笔记想把我踩过的坑、理顺的思路、以及那些让概念真正落地的手算例子都摊开来讲。内容主要集中在第二章(下):从黑体辐射定律到大气传输再到探测器响应,这些是一条完整的链路,任何一个环节的模糊都会在后面的系统设计里变成成倍误差。
1. 辐射度学那堆物理量,别再靠死记硬背了,先搞懂“在测什么”
学这一章之前,我老是被辐射通量、辐照度、辐射强度、辐射亮度这四个名词绕晕。事实上它们之间的关系特别简单,关键在于你问的是“总量、源、落点,还是方向上能看到多少”。
1.1 从总能量到“某根视线上的能量”
第一个物理量是辐射通量,单位是瓦特(W),表示单位时间内通过某个面的辐射能量,是一个纯功率概念。就好比说一个加热器每秒钟发出多少焦耳的热量,不问方向也不问距离。
接着是辐照度,单位是W/m²,等于辐射通量除以接收面的面积。它回答的问题是:落在单位面积上的功率有多少。对于红外系统来说,目标在探测器焦平面上形成的辐照度分布,直接决定了信号的大小,这是整个链路中最常被用来计算的量。
辐射强度则完全换了个视角,它描述点源在单位立体角内发出的辐射通量,单位是W/sr。它不受距离影响,因为立体角已经被归一化。但注意,只有在满足“点源近似”时才能直接使用辐射强度,比如极远处的无人机尾焰,或者恒星。等目标大到可以分辨时,我们就要引入辐射亮度。
辐射亮度是这四个量里最“全”的一个:单位是W/(m²·sr),描述单位投影面积、单位立体角内的辐射通量。它包含了观测方向和表面法线的夹角关系,也正因如此,黑体辐射公式里最常计算的就是这个量。我在笔记里给自己做了一张对照表,每次搞混就翻一下:
| 物理量 | 英文缩写 | 单位 | 关心什么 | 典型例子 |
|---|---|---|---|---|
| 辐射通量 | Φ | W | 总功率 | 加热器总输出 |
| 辐照度 | E | W/m² | 单位接收面积上的功率 | 焦平面接收到的信号 |
| 辐射强度 | I | W/sr | 单位立体角功率 | 点源目标亮度 |
| 辐射亮度 | L | W/(m²·sr) | 投影面积+立体角的综合辐射 | 面源目标成像 |
1.2 为什么很多红外计算要用“光子数”而不是“能量数”
第二章(下)花了不少篇幅把能量和光子速率区分开,我当时有点不解:既然都是表征辐射的物理量,为什么红外系统这么强调光子数?直到后来看红外探测器响应率定义才明白,光电探测器本质上是在“数电子”,光子被吸收后激发出载流子,探测器输出的信号电流与入射光子速率直接相关,而不是与入射能量直接相关。
比如同样功率的2μm短波红外辐射和10μm长波红外辐射,前者的单光子能量约为后者的5倍,这意味着在相同功率下,10μm波段的光子数要多得多,能激发更多的电子。如果误用能量值去估算探测器输出,会极大地低估长波方向的信号。正确的做法是优先把辐照度换算成光子辐照度(单位光子/s·m²),用普朗克关系 (E = h c / \lambda) 逐波段积分。我在做课后题时就遇到一个典型例子:计算某温度黑体在3~5μm与8~12μm两个波段的辐射通量比值,如果直接用能量积分,可能得到两者大致相当;但若考虑到光子速率,8~12μm在同等能量下的光子数是3~5μm的约两倍多,对探测器输出而言,长波信号远没有能量积分的那个印象里“那么小”。
这一下就解释了很多实战中的困惑:为什么制冷型长波探测器明明量子效率不高,却能拿到不错的NETD(噪声等效温差)?秘密就在光子速率上。所以,在学习这一章时,我非常建议把“能量表”和“光子数表”各算一遍,把两者的差异刻在脑子里。这种“双轨计算”的习惯,也让我后面看透了不少系统标称参数。
2. 黑体辐射定律,不要停留在公式背诵上,试着把曲线“读”成故事
黑体辐射这一节是整章的重头戏。教材先给出普朗克公式,再用它推出维恩位移定律和斯特藩-玻尔兹曼定律。但我认为最有用的不是公式本身,而是怎么从中提炼出“让系统工作在哪个波段”的直觉。
2.1 普朗克曲线的形状变化里藏着重要的工程判断
普朗克公式描述了黑体的光谱辐射出射度 (M_{\lambda}(T)) 随波长变化的关系。工程上最常见的两个结论是:
- 温度越高,辐射出射度整体越大,而且峰值波长越短;
- 曲线在高频段急剧上升后迅速衰减,随着温度升高,最大值的幅度增长极快。
我特意把300K、500K、1000K三条曲线画在同一张对数坐标下,才真正意识到斯特藩-玻尔兹曼定律里的四次方关系有多恐怖。300K时,300K物体在10μm附近的辐射远高于峰值波长9.66μm左右两侧的长波区域;而500K时,峰值波长约5.8μm,总辐射能是300K的7.7倍左右。到了1000K,峰值波长跳到2.9μm附近,能量又比500K增长了16倍。也就是说,一个目标温度只要从室温升到几百摄氏度,红外信号不是线性增长,而是四次方爆炸。
这解释了为什么红外制导系统最喜欢锁定的目标是飞机尾喷管和发动机舱壁——它们是整架飞机温度最高的部分,即使被遮挡了一部分,泄漏的热辐射仍然远大于蒙皮背景。相反,隐身设计的一个重要方向就是把排气温度降下来,哪怕只降10K,辐射总功率也能有百分之几十的变化。
维恩位移定律 ( \lambda_{\text{peak}} T = 2898\ \mu\mathrm{m\cdot K} ) 的工程含义非常直接:室温下目标的辐射峰值在10μm附近,所以被动红外系统通常会选择8~14μm波段;而高温目标(如尾焰)可能峰值在3~5μm甚至更短,于是有了中波红外探测的需求。
2.2 不要只背斯特藩-玻尔兹曼公式,要会做“波段积分”
斯特藩-玻尔兹曼定律 (M = \sigma T^4) 给出了全波段的积分辐射出射度,但实际系统都有工作波段限制,不可能接收全谱。于是我们需要学会计算“某波段内的辐射占全波段的百分比”。
这一步最常用的方法是引入 (x = \lambda T) 这个变量,把普朗克公式变成只与 (x) 有关的无量纲形式。我整理了教材里给出的几个重要分界点,方便随手估算:
| 物体温度 | 目标波长 | (x=\lambda T)(μm·K) | 该波长以下能量占比(约) |
|---|---|---|---|
| 300K | 3 μm | 900 | 0.0% |
| 300K | 5 μm | 1500 | 0.1% |
| 300K | 8 μm | 2400 | 约24% |
| 300K | 12 μm | 3600 | 约72% |
| 500K | 3 μm | 1500 | 0.1% |
| 500K | 5 μm | 2500 | 约29% |
上面这个表是我结合教材的辐射函数表整理的。看到这里就能明白一个重要的设计原则:对300K的目标,只看SWIR(1~3μm)几乎什么都看不到,因为能量占比趋近于零;而LWIR(8~14μm)则可以覆盖大部分能量。这也是为什么常温目标成像通常使用长波红外的主要原因。
坑也在这里:很多人以为中波(MWIR, 3~5μm)比长波(LWIR, 8~12μm)波长更短,所以分辨率天然更高,就一定更好。但其实对于300K的背景,MWIR本身的辐射能量很低,如果探测器噪声不压足够低,图像对比度会非常差。这里必须结合后面讲的探测器D*和大气透过率来做系统级的权衡。
3. 真实表面不等于黑体,发射率这个数,稍不留神就让测温结果偏得离谱
教材在引出黑体概念后,马上就转入了真实物体的发射率。这一节内容不长,却是几乎所有工程误判的温床。
3.1 发射率与基尔霍夫定律的实用版本
真实物体的辐射出射度可以写成 (\varepsilon(\lambda, T) M_{\lambda}(T)),其中 (\varepsilon) 是发射率,介于0到1之间。一个很重要的推论是基尔霍夫定律:在热平衡条件下,物体的光谱发射率等于光谱吸收率,而吸收率又等于 (1 - \rho - \tau)((\rho)为反射率,(\tau)为透射率)。对于大多数不透明固体,(\tau=0),所以 (\varepsilon = 1 - \rho)。
这个式子看起来简单,但它在实际测试中的含义非常深刻。我们测一个物体温度时,红外热像仪接收到的是来自物体自身的辐射,同时还有环境中其他物体反射到探测器里的辐射。如果目标表面反射率太高,那探测器看到的就不仅仅是目标本身的温度了,而是把旁边温热物体都“映”进去。这也是为什么很多红外测温场景需要“涂黑”目标表面,或者选用发射率较高的涂料来增强自发热辐射。
一个我印象很深的例子是:用热像仪测抛光铝板温度。铝板的发射率只有0.05左右,它在10μm附近的自身辐射极弱,但反射率接近0.95,因此热像仪测到的几乎是周围环境的反射辐射。如果误以为测到的是铝板本身,误差可能高达十几摄氏度,一不小心就会让人得出“铝板比室温高很多”的错误结论。
3.2 发射率随波长和方向变化,不是随便取个0.9就行
教材反复强调,工程中的发射率不是一个固定常数。许多介质在某个波段发射率高,在另一个波段发射率可能很低,比如白色油漆在可见光波段反射率高,但在长波红外波段其发射率可能接近0.9。而金属表面的发射率往往随波长增加而降低,随表面粗糙度增加而升高,还与观测方向有关。
所以在做系统分析时,至少不能把“发射率0.9”当成万能默认值。我给自己定了一条铁律:除非目标表面已经明确做了高发射率处理,否则在链路预算中一定要把发射率的不确定性留成一块“余量”,并且在项目报告中同时给出“假设发射率为0.9”和“假设发射率为0.7”的两种计算结果,看看最后结论是否还能站得住脚。这个习惯后来帮我避免了不少和需求方之间的扯皮。
另外,方向性与发射率的关系也是探测器型号选择时的潜在变量。教材提到朗伯表面(亮度不随方向变化的漫射面)的发射率与方向无关,但大量真实表面并非朗伯面。对重点目标进行多角度观测时,同一表面在垂直方向和掠射方向上的表观辐射亮度可能完全不同。如果在目标识别算法里把这个因素忽略不计,很容易造成目标特征分布不稳定,从而影响匹配率。
4. 大气窗口不是简单的“透光率”,而是三段命运的岔路口
黑体辐射告诉我们“目标发出什么”,大气传输决定了“我们能收到什么”。第二章(下)关于大气吸收和散射的内容,值得反复读三遍以上。因为实际系统的工作波段选择基本不是由探测器决定,而是被大气窗口卡死的。
4.1 吸收与散射,两种不同性质的“大气损耗”
大气对红外辐射的影响主要分两类:吸收和散射。吸收是分子将光子能量转化为内能的过程,具有明显的波长选择性。水蒸气、二氧化碳、臭氧都有各自的吸收带,比如水蒸气在2.7μm和6.3μm附近有强吸收,二氧化碳在4.3μm附近有很强吸收,这就导致某些波段在大气中几乎“全透不过去”。散射则是辐射被大气中的分子和气溶胶粒子改变传播方向,本身不损耗总能量,但是会让原来沿着视线传播到探测器的辐射被“打散”。
其中瑞利散射的强度与波长的四次方成反比,所以短波方向受影响更大,但这在近红外和红外波段相对可见光已小得多;气溶胶产生的米散射则对波长不太敏感,对红外成像的影响更为复杂。总的来说,大气对红外辐射的损耗是吸收和散射共同作用的结果,透过率随波长变化的曲线就像一条被啃噬的锯齿。
教材中刻画的典型大气窗口大致有三段:
- 短波红外(SWIR)窗口:约1~2.5μm;
- 中波红外(MWIR)窗口:约3~5μm;
- 长波红外(LWIR)窗口:约8~14μm。
这三个窗口之所以被系统设计反复利用,是因为它们恰好处在主要吸收分子的透明度相对较好的区域。但要注意,“窗口”不等于“全透”,窗口内依然存在许多窄带吸收吸收线,尤其在潮湿天气下,中波窗口会被水蒸气吸收迅速压缩至3.8~4.8μm左右。
4.2 大气透过率如何影响“目标看得清还是看不见”
我做一个粗略场景计算时,最有体感的是:即使目标辐射很强,如果大气路径上经过一个强吸收带,那么信号也会断崖式下降。比如在4.3μm附近二氧化碳吸收极强,如果系统没做滤光处理,信号噪声比会被严重拉低;这也是为什么许多气体探测红外相机专门把通道做在3.9~4.1μm之间,去避开吸收带。
大气路径还与观测距离和海拔有关。地面水平路径空气密度高且含有大量水汽,透过率随距离的衰减非常快;而斜向上看的空对地观测路径则能有效缩短等效大气路程。教材里给出了基于MODTRAN类计算程序得到的透过率曲线,但很多学生容易忽略“等效大气厚度”(也称大气路程)这个概念。同样的仰角,同一目标,地平高度差10度,透过率可能相差一倍以上。
我在做这一章习题时,动手算过一个简化例子:假设目标温度500K,波段选择MWIR 3.7~4.8μm,大气透过率0.6;如果波段变为LWIR 8~12μm,大气透过率0.7,但目标辐射按照普朗克积分来看,MWIR波段落在峰值附近,其波段内辐射贡献反而比长波大。可见“选择哪个窗口”没有一个固定答案,必须结合目标温度、探测距离、探测器灵敏度、以及天气统计特性整体仿真,不能单凭透过率比较下结论。
5. 探测器与光学系统的匹配,决定了链路预算最终能不能兑现
第二章(下)越往后越贴近系统设计,开始介绍探测器的主要性能参数,以及光学系统对辐射收集的影响。我以前看参数表只关心分辨率和帧率,这一章彻底改变了我的视角。
5.1 D*、响应率和NETD,这些指标各有各的适用范围
探测器性能的核心参数包括响应率 (\mathcal{R})、探测率 (D^*),以及由这些参数推导出的噪声等效温差(NETD)。
响应率描述的是输出信号与入射功率的比值,单位可以是V/W或A/W。对这个参数要特别留意它的“光谱响应曲线”:很多探测器在峰值波长的响应率很漂亮,但实际波段内的平均响应率可能只有峰值的50%。我们在做链路预算时,不能拿峰值响应率直接替代波段平均响应率,这一点吃过亏的人不在少数。
D*(比探测率)归一化了探测器面积和测量带宽,用来比较不同探测器之间的灵敏度。单位是cm·Hz^1/2·W^-1(琼斯)。D越高,探测器越灵敏。但要注意,D的测量条件(截止频率、背景温度、调制频率)不同,不能直接跨表比较。比如制冷型探测器的D*在长波方向可能高达 (10^{11}) Jones,非制冷探测器大多在 (10^8)~(10^9) Jones,两者相差数百倍,这背后是热噪声和暗电流的巨大差异。
NETD则是最直观的系统级指标,它定义为一个均匀目标区域中,能够被系统分辨的最小温差。换句话说,NETD越小,热像仪看到细微温度差别时越清晰。它把探测器噪声、光学透过率、像元面积、帧积分时间等全部打包成一个数,适合用于终端系统的横向比较。
5.2 光学系统的F数和衍射极限,是在给探测器“喂饭”
光学系统的作用是把目标的辐射会聚到探测器焦平面上。教材用“F数”来描述光学系统的集光能力:(f/# = f/D),其中 f 是焦距,D 是入瞳直径。F数越小,单位面积焦平面上的辐照度越高,信号越强。但也不能无限地减小F数,因为光学系统的衍射极限会给可分辨的最小光斑限定一个半径,典型公式是 (r_{\text{airy}} \approx 1.22 \lambda f/#)。
如果探测器像元尺寸远小于艾里斑半径,光学系统已经达到了衍射极限,再把像元做小并不能提高分辨率,只会增加噪声带宽;如果像元尺寸远大于艾里斑半径,则光学系统没有充分发挥能力,图像会显得模糊。工程上通常将艾里斑直径设计为像元尺寸的2~3倍,这样既兼顾采样率又不会因过采样牺牲灵敏度。教材里的这个匹配分析,让我第一次意识到,光学、探测器和信号处理需要作为一个整体来设计,而不是各自为政。
我在课后计算中实际演练了一个例子:一个长波红外系统,(\lambda=10\mu m),F=2,则艾里斑直径约为 (2.44\lambda F \approx 48.8\mu m)。如果选用的探测器像元尺寸是45μm,那么这就是一个接近衍射极限的“临界匹配”设计。但如果像元尺寸只有12μm,光学系统已经弥散到了4个像元,光能分不在多个像素上,单个像元信噪比反而下降,此时从系统设计角度应该重新评估光学孔径和探测器像元的组合,而不是一味追求小像元。
5.3 链路预算里最容易亏光的三个隐藏缺口
第一是光学系统的透过率不是100%,多片透镜每组表面都有反射和吸收,加上滤光片窄带设计的损耗,实际透过率可能只有70%~80%。有些人只在纸面上填一个85%的平均值,但没考虑冷光阑效率和滤光片的平均透过率,导致最终算出的信号比实测高了一半。
第二是探测器量子效率的带宽平均问题。探测器数据手册给出的峰值量子效率可能很高,但在所关心的波段内平均量子效率通常要打七折甚至五折。正确做法是逐波长相乘再积分。
第三是背景辐射造成的“直流项”。即使在成像时没有明显目标,探测器也会收到来自光学系统和环境的背景辐射。这部分光子会转化为光生载流子,增大散粒噪声。在NETD的计算公式里,如果忽略背景光电流,会低估噪声,最后测出来的NETD总比理论值要差。想弥补这一点,就必须用背景抑制结构(比如制冷型探测器配冷光阑)把背景辐照度压到最低。
这一段的复杂度和实用性都远超我的预期,做完两道完整的链路计算题之后,我对“为什么同一个红外机芯,装在两个不同口径的镜头上性能差异巨大”有了很直接的理解——不是探测器变强了,而是光学系统给探测器“喂”的饭更多了,同时背景噪声也更难压抑了。
6. 我的第二章(下)笔记复盘:概念地图+算例+错题集,三者缺一不可
既然这标题挂的是“学习笔记”,最后我还是想分享一点整理方法上的经验,也许比单纯背公式更能帮到正在啃这本书的朋友。
6.1 别按教材的顺序画脑图,按“辐射流”的顺序串联知识点
很多人喜欢按章节小节来记笔记,从定义1到定义2,再到定律1、定律2,这样的笔记复习时完全串不起来。我第二遍学习时换了个思路:以“辐射从目标出发,穿过大气,经过光学,打到探测器”这条链路为主线,把每个定义、每个定律挂在这个链路对应的位置。这样复习时,只要脑中复现这条链路,所有公式和参数就自然归位了。
举个例子,目标表面的辐射特性用发射率和辐射亮度描述;离开目标后经过大气传输,用了大气透过率;到达光学系统后,被F数和入瞳面积收集;最后聚焦到焦平面,形成辐照度分布;探测器再把辐照度转化为响应电压,除以总噪声得到信噪比。这张流程图上,每一个节点都挂上公式、量纲要义和常考易错点,复习效率比线性笔记高得多。
6.2 亲手做三道完整计算题,胜过抄十遍公式
只算例题是不够的,我特别建议自己变换条件重算。比如教材里有一个例题是求300K黑体在4.5μm处的光谱辐射出射度,大家直接套公式结果都一样。但把温度改成350K后,计算结果和300K时相差了好几倍,这时候你才会真正体会到普朗克公式里的指数项在温度变化时的敏感度。另一个推荐的练习是:给定探测器D*、光学透过率、目标温度和发射率,自己推导并计算出信噪比和NETD。这个综合练习把整章所有知识点全部用上,做完一遍以后再回头听课或读其他资料,很多疑惑都能得到解释。
6.3 错题本里最值得记录的“思维陷阱”
我在错题本里整理了三个最典型的坑,这里直接分享出来:
- 把辐射强度 (I) 和辐射亮度 (L) 混为一谈,不知道辐射强度要乘以观测方向的投影面积才能变成辐射亮度;
- 用全波段辐射出射度替代波段内辐射出射度,导致计算探测器接收功率时大幅度偏大或偏小;
- 忽略了温度单位必须用开尔文,所有符号代入前先换算,否则结果差得离谱。
有一次我用25℃代替298K去算斯特藩-玻尔兹曼,结果还算是有理有据的,但温度一重复计算就完全乱套了。从那以后我习惯把温度大写字母K写在公式旁边,避免自己掉进“摄氏度惯性”。
6.4 下一步学习的衔接思路
第二章(下)学完,算是把辐射学的基础铺好了。接下来第三章讲典型红外系统结构时,就能很自然地对应上链路中的各个环节。我也给自己列了一个延伸方向:结合某种热成像系统,用小目标信噪比公式重新推导一遍探测距离。这个过程虽然只有几步加减乘除,但对建立“从物理量到系统指标”的直觉帮助极大。如果你也在读这本经典教材,不妨沿着这个思路试试,绝对能比单纯做题收获更多。