[发明专利]一种日盲紫外遥感相机绝对辐射定标系数在轨修正方法

权利要求书

1.一种日盲紫外相机绝对辐射定标系数在轨修正方法，其特征在于步骤如下：

(1)实验室辐射定标，获得绝对辐射定标方程和非均匀校正系数；

111)相机绝对辐射定标

调节相机光源辐亮度由小变大依次为L₁,L₂……L_N,其中L₁<L₂<……<L_N，对应的相机输出平均响应依次为其中得到一组实验室定标点序列，采用最小二乘法线性拟合得到：

L＝KV+C,其中

其中K和C为拟合系数，将相机输出平均响应依次代入上述方程计算辐亮度L值，并计算该相机平均响应下定标残差：

ϵi=L-LiLi=KV‾i+C-LiLi;]]>

依次检查每一点定标残差是否满足|ε_i|<5％，若不满足，则从定标点中剔除该定标点，否则保留该定标点，保留的定标点组成新的定标点序列，如果本次检查中出现剔除定标点，则对新的定标点序列重复最小二乘线性拟合与检查定标残差，直至本次检查中未出现剔除定标点，则本次最小二乘法线性拟合得到的方程为相机绝对辐射定标方程：

L＝KV+C,其中

拟合系数K和C为相机绝对辐射定标系数，分别表示增益和偏置，为定标点序列最小平均响应，为定标点序列最大平均响应，为相机线性响应范围；

112)像元绝对辐射定标

根据相机绝对辐射定标方法，对相机进行逐像元绝对辐射定标，得到相机各像元绝对辐射定标方程为：

L＝k_i,jv_i,j+c_i,j；

其中v_i,j表示像元(i,j)的响应，k_i,j和c_i,j表示像元(i,j)的增益和偏置；

113)非均匀性校正

非均匀校正系数计算公式为：

Gi,j=ki,jK]]>

Qi,j=-CK+Gi,j·ci,jki,j;]]>

其中G_i,j和Q_i,j为非均匀校正系数，即非均匀校正增益和非均匀校正偏移；

(2)获取相机在轨成像参数，包括成像模式、成像时间T₁、卫星轨道根数(a,e,i,Ω,ω,M₀)和在轨图像；

(3)判断相机在轨成像模式是否为星下点成像，若满足，则进入步骤(4)，

否则跳回步骤(2)重新获取新的相机在轨成像参数；

(4)根据成像时间T₁和卫星轨道根数(a,e,i,Ω,ω,M₀)计算获得卫星星下点地理纬度

(5)判断T₁时刻卫星是否位于中纬度地区，判据为：

若满足判据条件，则进入步骤(6)，否则跳回步骤(2)重新获取新的相机在轨成像参数；

(6)计算获得T₁时刻卫星星下点太阳高度角h_s

其中，δ₀表示太阳赤纬，φ表示为星下点纬度，为太阳时角；

(7)判断T₁时刻卫星是否满足光照条件，判据为：

h_s>20°；

若满足判据条件，则进入步骤(8)，否则跳回步骤(2)重新获取新的相机在轨成像参数；

(8)根据T₁时刻卫星星下点纬度星下点太阳高度角h_s，计算获得相机入瞳辐亮度

(9)对在轨图像进行相对非均匀校正，校正公式为：

V_i,j＝G_i,j×v_i,j+Q_i,j；

v_i,j为在轨图像校正前的响应值，V_i,j为在轨图像校正后的响应值；

(10)根据相机地面分辨率a×b，计算获得可用于背景辐亮度反演的最大定标区域行数N_a和列数N_b：

Na=[50a]]]>

Nb=[50b];]]>

其中[·]表示取整；以非均匀校正后的在轨图像中心像素为中心，截取大小为N_row×N_col的矩形区域，N_row、N_col为正整数，取值范围分别是N_row≤N_a和N_col≤N_b，计算获得截取定标区域内图像平均灰度值

(11)判断图像平均灰度值是否位于响应线性范围，判据为：

V‾min<VorbitT<V‾max,]]>

若满足判据条件，则根据相机绝对辐射定标系数K和C，反演入瞳辐亮度

LcalibrationT=(VorbitT-C)/K;]]>

否则跳回步骤(2)重新获取新的相机在轨成像参数；

(12)改变观测时刻依次为T₁,T₂,…T_m，m为大于2的正整数，按照步骤(3)～步骤(8)计算相机入瞳辐亮度序列根据相机对应的绝对辐射定标系数和图像平均响应其中按照步骤(9)～步骤(11)反演入瞳辐亮度序列采用最小二乘法线性拟合得到一元线性回归方程系数K_orbit与C_orbit：

L_orbit＝K_orbitL_calibration+C_orbit；

(13)利用步骤(12)求得的一元线性回归方程系数对步骤(1)得到的相机绝对辐射定标方程进行修正，得到修正后的绝对辐射定标方程为：

L＝K_orbitKV+K_orbitC+C_orbit；

其中K_orbitK为增益，K_orbitC+C_orbit为偏置。