傅里叶变换光谱仪的复数谱多点辐射定标方法
【技术领域】
[0001] 本发明涉及光电标定技术领域，特别是指一种适用于傅里叶变换光谱仪的复数谱 多点辐射定标方法。
【背景技术】
[0002] 傅里叶变换光谱仪是以双光束干涉测量为基础，利用傅里叶变换原理实现光谱测 量的仪器设备。傅里叶变换光谱仪由于具有高信噪比、高分辨率、快速测量等优点，在化工、 地质、医药、环境、气象、军事、空间技术等领域取得了广泛的应用。现有的傅里叶变换光谱 仪实际测量得到的是入射光的干涉图。在理想状态下，干涉图是对称的，且通过干涉图可以 经过傅里叶变换可以得到实数谱。但是实际中由于仪器误差、环境干扰等因素的影响，干涉 图并不对称，导致光谱仪测量值为复数。
[0003]目前传统的光谱数据辐射定标方法多基于实数光谱，通常是对复数数据进行实数 化处理（如取实部或者幅度值），或是先相位校正然后实数化处理。这样的处理方法会产生 一定的误差，特别是对于相位偏差较大的数据误差更为显著。
【发明内容】
[0004] 针对现有技术中的傅里叶变换光谱仪的存在的问题，本发明要解决的技术问题是 提供一种适用于傅里叶变换光谱仪的复数谱多点辐射定标方法，以获取更优的定标精度。
[0005] 为了达到上述目的，本发明实施例提出了一种傅里叶变换光谱仪复数谱多点辐射 定标方法，包括：
[0006] 步骤1 :分别获取傅里叶变换光谱仪的η个温度点对应的测量电压值，其中 η彡2 ;
[0007] 步骤2 :通过以下公式确定每一温度点的光谱福亮度Li (V):
[0009] 其中ε为所述每一温度点的黑体发射率，第一辐射常数C1 = 1. 1910 X 10 12(W · cm2)，第二辐射常数C2= 1. 4387 (Κ · cm)，T为所述每一温度点的黑体温 度；
[0010] 步骤3、通过最小二乘祛 求解，以得到傅里叶 变换光谱仪在各波数V下的实际线性增益巧〇〇和实际偏移量
[0014] 步骤4、通过以下公式获取傅里叶变换光谱仪输入的目标光谱辐亮度
[0016] 作为上述技术方案的优选，对于η = 2的两点辐射定标情况，所述傅里叶变换光谱 仪测量黑体的一个高温温度点和一个低温温度点，且所述傅里叶变换光谱仪得到对应的两 个测量电压值：高温测量电压值5?Η、低温测量电压值瓦:；
[0017] 其中所述步骤2具体为：通过公式
[0019] 分别计算且所述高温温度点的光谱辐亮度为LH，低温温度点的光谱辐亮度为L c。
[0020] 作为上述技术方案的优选，对于η = 2的两点辐射定标情况，步骤3具体为：
[0021] 通过以下公式计算所述傅里叶变换光谱仪的在各波数V下的实际线性增益G(V) 和实际偏移量J(V):
[0024] 本发明的上述技术方案的有益效果如下：
[0025] 上述方法直接以干涉图傅里叶变换后得到的原始复数谱为基础进行辐射定标处 理，无需实数化或相位校正等步骤，提高了处理效率与精度；且获得的虚数谱可以作为误差 分析与精度评价的重要依据。本发明实施例的方法为设备的误差来源与精度分析研究提供 新的技术途径；且能够支持两点、多点方式的辐射定标，对于黑体文件的选取具有更好的灵 活性。本发明实施例支持所有基于干涉与傅里叶变换原理的光谱设备定标，不仅对光谱仪， 光谱成像仪也同样适用。
【附图说明】
[0026] 图1为本发明实施例的流程示意图；
[0027] 图2和图3为利用本发明实施例的方法进行对30°(：、40°(：、60°(：、70°(：黑体文件计 算出仪器长波探测器的响应函数巧(V)和〇 ^ 和〇 ^ / 的幅度谱和虚部谱；
[0028] 图4和图5为利用本发明实施例的方法中通过和对50°C黑体文件计算 得到L(V)的幅度谱和虚部谱。
【具体实施方式】
[0029] 为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具 体实例进行详细描述。
[0030] 本发明的原理如图1所示。
[0031] 在傅里叶变换光谱仪正常状态时，傅里叶变换光谱仪可以近似为一个线性响应 系统。即：当系统的输入为目标光谱辐亮度L(V)时（目标光谱辐亮度L(V)的单位为W/ cm2 · sr · cm》，则系统的输出为傅里叶变换光谱仪得到的各波数V (单位cm 3下的测量 电压值M(V)(单位为V/cm)，而输入与输出之间的关系即为该傅里叶变换光谱仪的响应函 数。因此辐射定标的过程也就是一个求解仪器响应函数的过程。
[0032] 由上述原理可知，对于傅里叶变换光谱仪这种线性响应系统，其响应函数可表示 为：M(v) = G(v) · L (V)+0(v) (1)。
[0033] 其中两个未知系数G(V)和O(V)分别代表傅里叶变换光谱仪在各波数V下的线性 增益（Gain)和偏移量（Offset)，线性增益和偏移量共同表征了仪器响应函数。L(V)为傅 里叶变换光谱仪输入的目标光谱辐亮度。
[0034] 实际中，傅里叶变换光谱仪的测量值为复数。同时导致响应函数中的线性增益和 偏移量也为复数。因此傅里叶变换光谱仪的实际响应函数可表示为：
[0035] M(v) = G (v) · L(v) + O (v) (2)。
[0036] 其中泣〇〇为实际测量的复数的测量电压值，巧和^?分别代表傅里叶变换光 谱仪在各波数V下的实际线性增益（Gain)和实际偏移量（Offset)，其中实际线性增益和实 际偏移量也都是复数。
[0037] 由公式（2)可知，傅里叶变换光谱仪输入的目标光谱辐亮度
[0039] 因此在求得实际线性增益巧(V)和实际偏移量G(V)的情况下，便可通过公式（3)将 实际测量的复数的测量电压值定标处理为目标光谱辐亮度L(V)，辐射定标的过程便 转化为了求解仪器响应函数。而想要求解实际线性增益巧〇〇和实际偏移量就需要获 得黑体标定文件。即：
[0040] mXv) = G{v)-LXv) + 〇{v) (i=l,2,...,η) (4);
[0041] 其中η为黑体测量的温度点（η彡2);
[0042] 且其中每一个温度点的黑体光谱辐亮度L1 (V)可由公式（5)求得：
[0044] 其中ε为黑体发射率，第一辐射常数C1= I. 1910Χ 10 12(W · cm2)，第二辐射常数 C2= 1. 4387 (K · cm)，T为黑体温度（单位为K)。
[0045] 由此可见定标过程转化为根据黑体参数求得的L1 (V)和光谱仪测得的泣^v)解求 巧(V)和J(V)。这个基本的一次线性拟合问题可以通过最小二乘法求解，BP :
[0047] 求解结果为：
[0051] 其中，下标A表示平均。
[0052] 下面通过一个具体的例子来进行说明。在实际使用中常见的是通过测量黑体的两 个温度点（一个高温温度点和一个低温温度点）来进行标定，即η = 2的情况。这样，傅 里叶变换光谱仪就会得到对应的两个测量电压值：高温测量电压值&、低温测量电压值 &。通过以上的公式就可以分别计算出两个温度点的黑体辐亮度LH、LC，计算公式为：
[0055] 需要注意的是目标测量的相位状态需要与各黑体测量的相位状态一致，即目标文 件与标定文件扫描方向应保持一致，如果扫描方向不一致的数据间用于定标处理，结果可 能导致错误。
[0056] 在求得可v)和可V)的情况下，根据（3)式便可计算得到标定后的光谱辐亮度 L(V)。只有在仪器响应完全线性及完全理想测量条件下，L(V)才为实数，实际中通常取其 复数值的幅度作为最终处理结果，同时虚部谱是设备本身相位误差及测量误差共同积累的 结果，可以作为误差分析与精度评价的重要依据。
[0057] 根据目标温度范围和实际需要，利用光谱仪对不同温度点条件下的黑体进行测 量，以上一章中的实验为例进行说明：
[0058] (1)利用本方法对30°C、40°C、60°C、70°C黑体文件计算出仪器长波探测器的响应 函数巧〇0和歹(V)，巧(V；)和^?/巧(V；)(偏移函数/增益函数反映仪器自身的辐亮度）各自幅 度与虚部谱如图2、图3所示。
[0059] (2)利用句V)和O以对50°C黑体文件计算得到L(V)的幅度与虚部谱如图4、图5 所示。从图4中可以看出本方法处理结果与理论计算结果匹配度很好，仅2360cm 1附近CO2吸收带处误差相对较大。从图5中可以看出L(V)的虚部不完全为0,整体比幅度低约2个 量级，但吸收带处明显增大。
[0060] 结果说明，本方法计算得到的辐亮度结果（幅度谱）准确有效，同时虚部谱又可以 作为误差分析与精度评价的依据。
[0061] 以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员 来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也 应视为本发明的保护范围。
【主权项】
1. 一种适用于傅里叶变换光谱仪的复数谱多点辐射定标方法，包括步骤： 步骤1 :分别获取傅里叶变换光谱仪的n个温度点对应的测量电压值.只^， 其中n彡2 ; 步骤2 :通过以下公式确定每一温度点的光谱辐亮度L1 (V):其中e为所述每一温度点的黑体发射率，第一辐射常数C1= 1.1910X1012(W^cm2), 第二辐射常数C2= 1. 4387(K?cm)，T为所述每一温度点的黑体温度；2. 根据权利要求1所述的适用于傅里叶变换光谱仪的复数谱多点辐射定标方法，其特 征在于，对于n= 2的两点辐射定标情况，所述傅里叶变换光谱仪测量黑体的一个高温温度 点和一个低温温度点，且所述傅里叶变换光谱仪得到对应的两个测量电压值：高温测量电 压值Mff、低温测量电压值Mc7 ; 其中所述步骤2具体为：通过公式分别确定所述高温温度点的光谱辐亮度为Lh，低温温度点的光谱辐亮度为Lc。3. 根据权利要求2所述的适用于傅里叶变换光谱仪的复数谱多点辐射定标方法，步骤 3具体为： 通过以下公式确定所述傅里叶变换光谱仪的在各波数V下的实际线性增益巧〇〇和实 际偏移量0(v):
【专利摘要】本发明公开了一种傅里叶变换光谱仪复数谱多点辐射定标方法，包括：步骤1、分别获取傅里叶变换光谱仪的n个温度点对应的测量电压值其中n≥2；计算每一温度点的光谱辐亮度计算傅里叶变换光谱仪在各波数v下的实际线性增益和实际偏移量<maths num="0001"></maths>通过以下公式获取傅里叶变换光谱仪输入的目标光谱辐亮度
【IPC分类】G01J3/45
【公开号】CN105043548
【申请号】CN201510130410
【发明人】何茜, 王广平, 武敬力, 张亚洲, 雷浩, 刘畅
【申请人】北京环境特性研究所
【公开日】2015年11月11日
【申请日】2015年3月24日