一种中波红外焦平面阵列探测器的相对辐射校正方法
【技术领域】
[0001] 本发明涉及中波红外焦平面阵列探测器的相对辐射校正方法，特别涉及一种具备 星上黑体定标装置的中波红外焦平面阵列探测器相对辐射校正方法。
【背景技术】
[0002] 地表的温度一般为300K，地表辐射能量基本上处在3. 0ym以上的波段，中波红外 焦平面阵列探测器利用空间光学系统收集3. 0ym~5. 0ym的中波红外辐射能量，其具有 穿透烟、尘、雾、雪以及识别伪装的能力，可以进行白天和黑夜的准全天候观测，这些特点使 得它获得的高精度面阵热辐射信息在实现导弹预警、红外侦察等军事方面以及森林火灾、 地表温度、洋面温度、昼夜云、海岸线监测等民用方面具有的广阔应用。
[0003] 但是由于中波红外焦平面阵列探测器受材料和工艺水平的限制，以及光电接收阵 列受空间和探元耦合影响很大，使得中波红外焦平面阵列探测器各探测单元响应特性之间 存在普遍的非均匀性，在图像上表现为空间噪声，它将导致对相同地表成像时探测器输出 的数字DN值间存在很大差别，以致使其难以满足中红外成像的应用需求。因此必须将中波 红外焦平面阵列探测器的非均匀性进行归一化相对辐射校正，使各探元输出值调整到同一 基准使得各探元对相同地表具有相同的数字DN输出值。
[0004] 目前，中波红外焦平面阵列探测器相对辐射校正常采用内部温度参考源的星上定 标法（内定标法），即在星上配置可变高温黑体辐射源和低温黑体辐射源，依据定标指令依 次完成对黑体的成像，定标成像完成后利用两点法或多点法实现阵列探测器各探测元响应 相对辐射校正。虽然基于两点法或多点法的相对辐射校正可以有效地去除空间噪声中由热 噪声和散粒噪声等白噪声造成的中波红外焦平面阵列探测器非均匀性，但是无法去除空间 噪声由g-r噪声和布朗噪声等广义1/f噪声造成的中波红外焦平面阵列探测器非均匀性， 故其校正精度不高。
【发明内容】
[0005] 本发明所要解决的技术问题是：克服现有技术的不足，提供一种中波红外焦平面 阵列探测器的相对辐射校正方法，解决了【背景技术】里中波红外焦平面阵列探测器相对辐射 校精度不高的技术问题。
[0006] 本发明的技术方案是：一种中波红外焦平面阵列探测器的相对辐射校正方法，步 骤如下：
[0007] 1)建立探元辐射响应模型：依据中波红外焦平面阵列探测器的成像原理和成像 链路的噪声产生机理，建立中波红外焦平面阵列探测器的辐射响应模型；
[0008] 2)星上黑体定标成像：依次完成对可变高温黑体辐射源和低温黑体辐射源的成 像，获得相应的阵列探测器数字DN值；
[0009] 3)噪声分离：对步骤2)得到的阵列探测器数字DN值进行小波变换，依据广义1/ f噪声的统计自相似特性，估计广义1/f?噪声的小波系数，然后利用小波软阈值方法对得到 的小波系数进行收缩，最后对收缩后的小波系数进行小波逆变换得到分离广义1/f噪声后 的中波红外焦平面阵列探测器数字DN值；
[0010] 4)噪声估计：利用步骤2)得到的可变高温黑体和低温黑体定标数据的数字DN值 和步骤3)得到的分离广义1/f噪声后的数字DN值，计算噪声大小；
[0011] 5)噪声去除：利用步骤4)得到的噪声，在步骤2)得到的可变高温黑体和低温黑 体数据的数字DN值中进行噪声去除，即得到去噪后的可变高温黑体和低温黑体定标数据；
[0012] 6)计算获得相对辐射校正系数：利用步骤5)的去噪后的可变高温黑体和低温黑 体数据，计算获得相对辐射定标系数；
[0013] 7)相对辐射校正：对阵列探测器获得的图像数据，利用步骤4)得到的噪声和步骤 6)得到的相对辐射校正系数，获得阵列探测器的相对辐射校正图像。
[0014] 本发明与现有技术相比的有益效果是：
[0015] (1)本发明利用中波红外焦平面阵列探测器成像时的噪声特点，提出了采用小波 分析法去除广义1/f?噪声，实现了中波红外焦平面阵列探测器高精度的相对辐射校正技术 问题。
[0016] ⑵本发明为中波红外焦平面阵列探测器提出的方法，兼顾了中波红外焦平面阵 列探测器的成像特点与可实行性，对太阳同步轨道卫星或地球同步轨道卫星上具备星上黑 体定标装置的中波红外焦平面阵列探测器遥感图像校正具有很强的工程应用意义。
【附图说明】
[0017] 图1为本发明方法流程图。
【具体实施方式】
[0018] (1)探元辐射响应模型确立：依据中波红外焦平面阵列探测器的成像原理，充分 考虑成像链路的噪声来源，建立中波红外焦平面阵列探测器的辐射响应模型；
[0019] DN(m,n) =f[L(m,n) ] +Nf (m,n) +NW (m,n) ①
[0020] =f[L(m，n) ]+N(m，n)
[0021] 式中：m，n为中波红外焦平面阵列探测器的探元位置mG[1，M]，
[0022] nG[1，N];
[0023] DN(m，n)为探元（m，n)的原始数字DN值；
[0024] f[L(m，n)]为探元（m，n)接收的福射亮度L(m，n)函数；
[0025] Nf (m，n)为探元（m，n)的广义1/f噪声；
[0026] Nw(m，n)为探元（m，n)的白噪声；
[0027] N(m，n) =Nf (m，n)+Nw(m，n)为探元（m，n)的总噪声。
[0028] (2)星上黑体定标成像：依据定标指令依次完成对可变高温黑体辐射源和低温黑 体辐射源的成像，获得相应的中波红外焦平面阵列探测器的原始数字DN值；
[0029] (3)噪声分离：对步骤⑵的数字DN值进行小波变换，依据广义1/f?噪声的统计 自相似特性，估计广义1/f?噪声的小波系数，然后利用小波软阈值方法对得的小波系数进 行收缩，最后对收缩后的小波系数进行小波逆变换得到分离广义1/f噪声后的中波红外焦 平面阵列探测器数字DN值；
[0030] (3. 1)小波变换
[0031] 对步骤⑵的数字DN值进行小波变换，可获得不同小波尺度下的总噪声小波系数 模型为：
[0032] N]{k) =N({k) +Nl{k) ②
[0033] 式中：j为小波尺度，jG[1，J];
[0034] k为第k个小波变换区域，以探元（m，n)为中心的2外像元，kG[1，K]。
[0035] (3. 2)参数估计
[0036] 依据广义1/f噪声的统计自相似特性，其不同小波尺度下