专利名称：一种基于绝对辐射定标的多光谱影像辐射校正方法
技术领域：
本发明涉及一种基于绝对辐射定标的多光谱遥感影像辐射校正方法，属于卫星图像处理技术领域，用于光学遥感卫星多光谱相机影像辐射校正处理。
背景技术：
光学遥感影像辐射校正是光学遥感卫星相机影像处理的主要内容之一，其目的是减少或消除遥感影像的各类辐射误差，提升遥感影像辐射精度，增强遥感影像反映实际地物真实物理特性的能力，其内容主要包括两个方面一是消除影像自身像元之间的辐射响应不一致性，称之为相对辐射校正；二是将遥感影像的量化维度由传感器器件的数字响应计数值即DN值转化为实际地物辐射能量值，称之为绝对辐射校正。以上目的与内容适用于采用CCD作为探测器件的可见光全色、近红外、多光谱等各类常用光学遥感影像的辐射校正处理。·对于多光谱遥感影像，其辐射校正处理除具备上述光学遥感影像辐射校正处理的一般特征以外，还应在提升单谱段(波段)遥感影像自身辐射精度的基础上，进一步考虑各谱段(波段)影像之间的辐射相对比例关系与绝对辐射量化程度，减少或消除由于传感器器件光谱响应不一致性、大气光谱透射能力不一致性等光谱响应误差，提升多光谱遥感影像随波长变化的光谱响应精度，一方面在影像的目视观测上提升不同谱段间色彩比例与真实光谱响应程度比例的一致性即相对光谱辐射精度，另一方面在遥感影像量化表达方式上将光谱响应量化程度由传感器器件不同谱段(波段)的数字响应计数值响应曲线转化为实际地物不同谱段(波段)的辐射能量值曲线，增强多光谱影像反映实际地物光谱反射变化真实特性的能力即绝对光谱辐射精度。光学遥感影像绝对辐射校正的主要方法是对遥感卫星相机进行绝对辐射定标，建立相机记录的数字信号计数值(DN值)与对应的辐射能量之间的数学关系。这种数学关系一般情况下采用线性方程表示，该方程包括增益与偏移量两个系数，即绝对定标系数。绝对辐射定标数学公式如下Li = aj XDNj+bj (i = I, 2…η)式中L为辐亮度，DN为响应DN值，i表示多光谱谱段序号。绝对辐射定标的输入条件是遥感卫星相机成像各谱段不同等级的光学辐射能量值序列及其所成影像上对应DN值序列，其中DN值可由直接获取积分球等标准辐射源或靶标在遥感影像对应的DN值得到，因此绝对辐射定标输入参数的重点在于光学辐射能量值的获取。光学辐射能量值可采用辐亮度、辐照度、反射率等参数，其中最常用的参数是辐亮度参数，其单位为W · cm-2 · sr'光学辐射能量值的获取方法主要有三种实验室定标方法、星上在轨定标方法和地面检校场在轨定标方法。实验室定标是指卫星发射前在实验室中应用光谱、辐射等定标系统形成标准辐射源(如积分球定标源、光谱扫描定标等)对遥感卫星相机进行定标；星上在轨定标主要是在卫星发射后在轨飞行期间利用星上定标灯或太阳、月球、恒星等外界辐射光源对遥感卫星相机进行定标；地面检校场在轨定标是通过量测地面铺设的标准反射率祀标的反射福射特性对遥感卫星相机进行定标。常用的绝对辐射校正方法，对于实验室定标数据与地面检校场在轨定标数据的使用是相互孤立的，一般在未开展地面检校场在轨定标之前采用实验室定标数据计算获得的定标系数进行校正，而在开展地面检校场在轨定标之后则直接采用地面检校场在轨定标数据计算获得的定标系数进行校正。直接采用实验室定标数据计算获得的定标系数进行校正的方法局限性在于由于卫星发射后在轨运行阶段与卫星发射前研制阶段相机的辐射响应特性可能存在差异，导致实验室辐射定标结果往往与卫星在轨运行时的相机辐射响应状态不一致，同时实验室辐射定标结果没有考虑地物特性、大气条件、太阳光照条件等外界环境因素的影响，导致实验室辐射定标结果用于卫星在轨运行期间的绝对辐射校正精度难以保证。直接采用地面检校场在轨定标数据计算获得的定标系数进行校正的方法局限性在于由于相机器件特性的限制，影像DN值与辐亮度之间往往是非线性关系，为此需要足够数量的样本点(对应于一定环境条件下不同反射率靶标所对应的标准辐亮度)才能够精 确地描述这种非线性关系；而由于在轨靶标布设可行性、实用性、经济性等限制，地面检校场所布设靶标的反射率等级(对应于标准辐亮度等级)是有限的，往往不能满足精确描述DN值-辐亮度非线性关系的要求，因此将影像DN值直接转化为在轨辐亮度的绝对辐射校正方法往往校正精度不够，制约了定量反演精度的提升。地面检校场在轨定标因靶标反射率等级设置较少往往难以精确地反映相机辐射响应的非线性特性，而实验室定标方法可以较为精确地描述相机辐射响应的非线性特性。为此，应综合利用实验室定标数据与在轨场地定标数据，将实验室定标能够精确描述相机响应非线性特性的优势与在轨场地定标能够真实反应卫星实际成像目标地物物理信息的优势结合起来。
发明内容
本发明的技术解决问题是克服现有技术的不足，提供一种基于绝对辐射定标的多光谱影像辐射校正方法，能够提高遥感卫星多光谱影像的定量反演精度，使多光谱影像的色彩比例真实地反映地物实际光谱特性，减少或消除谱段间的色彩比例偏差。本发明的技术解决方案是一种基于绝对辐射定标的多光谱影像辐射校正方法，实现步骤如下(I)相对辐射校正基于遥感卫星多光谱相机实验室相对定标系数，采用线性变换方法对多光谱原始影像逐谱段逐像元进行相对辐射校正，消除各CCD探测器单元即探元间的辐射响应不一致性，输出相对辐射校正影像；(2)实验室绝对辐射定标获取实验室积分球辐射定标源在不同遥感卫星多光谱相机的工作参数下的多级辐亮度与对应DN值，构建并输出实验室辐亮度-影像DN值查找表；(3)基于实验室绝对定标结果的绝对辐射校正获取步骤(I)得到的相对辐射校正影像成像时刻的遥感卫星多光谱相机的工作参数，基于该工作参数得到对应的实验室辐亮度-影像DN值查找表，通过分段线性插值方法计算各区间对应的实验室绝对辐射定标系数，同时获取相对辐射校正影像中各谱段各像元所在的实验室辐亮度-影像DN值区间，采用所在区间的实验室绝对定标系数对相对辐射校正影像逐谱段逐像元进行实验室绝对辐射校正，输出实验室绝对辐射校正影像；(4)地面检校场在轨绝对辐射定标获取地面检校场不同反射率靶标的在轨辐射定标野外测量数据并以此计算在轨辐亮度，同时获取靶标所在的采用步骤(3)方法处理得到的实验室绝对辐射校正影像，测量不同反射率靶标在该影像上对应的实验室辐亮度，构建并输出不同反射率靶标的在轨辐亮度-实验室辐亮度查找表；(5)基于在轨绝对定标结果的绝对辐射校正获取步骤(3)得到的实验室绝对辐射校正影像成像时刻的遥感卫星多光谱相机的工作参数，基于该工作参数得到对应的在轨辐亮度-实验室辐亮度查找表，通过分段线性插值方法计算各区间对应的在轨绝对辐射定标系数，同时获取实验室绝对辐射校正影像 中各谱段各像元所在的在轨辐亮度-实验室辐亮度区间，采用所在区间的在轨绝对定标系数对实验室绝对辐射校正影像逐谱段逐像元进行在轨绝对辐射校正，输出在轨绝对辐射校正影像。所述步骤(3)中具体实现如下(3. I)获取步骤(I)得到的相对辐射校正影像成像时的遥感卫星多光谱相机工作参数，所述工作参数包含积分级数、增益、积分时间、数字增益和箝位，并根据所述工作参数获取对应的步骤(2)得到的实验室辐亮度-影像DN值查找表；(3. 2)基于步骤(3. I)得到的实验室辐亮度-影像DN值查找表，获取各实验室辐亮度-影像DN值区间的实验室辐亮度上限与下限和对应的影像DN值上限与下限；(3. 3)根据步骤(3. 2)得到的各区间实验室辐亮度上限与下限及对应的DN值上限与下限，以线性方程为数学模型，解算各区间辐亮度-DN值线性方程的系数，所述系数包括增益系数与偏置量系数，所述系数作为各区间实验室绝对辐射定标系数；(3. 4)获取步骤⑴得到的相对辐射校正影像中各谱段各像元的DN值，判断其在实验室辐亮度-影像DN值查找表中对应的实验室辐亮度-影像DN值区间，并进而获取所在区间的实验室绝对福射定标系数；(3. 5)根据步骤(3. 4)得到的包括增益系数与偏置量系数在内的实验室绝对辐射定标系数，采用线性变换方法计算该像元DN值对应的实验室辐亮度。所述步骤(4)中具体实现如下(4. I)获取地面检校场不同反射率靶标的各谱段在轨辐射定标野外测量数据；(4. 2)根据在轨辐射定标野外测量数据计算地面检校场在轨定标各谱段各级反射率靶标的在轨辐亮度；(4. 3)获取靶标所在的采用步骤(3)方法处理得到的各谱段实验室绝对辐射校正
影像；(4. 4)在靶标所在的实验室绝对辐射校正影像上，测量不同反射率靶标对应的实
验室福売度；(4.5)根据不同反射率靶标的地面检校场在轨辐亮度与实验室辐亮度的映射关系，构建在轨辐亮度-实验室辐亮度查找表。
所述步骤(5)中具体实现如下(5. I)获取步骤(3)得到的实验室绝对辐射校正影像成像时的遥感卫星多光谱相机工作参数，所述工作参数包括积分级数、增益、积分时间、数字增益和箝位，并根据所述工作参数获取对应的步骤(4)得到的在轨辐亮度-实验室辐亮度查找表；(5. 2)基于步骤(5. I)在轨辐亮度-实验室辐亮度查找表，获取各在轨辐亮度-实验室辐亮度区间的在轨辐亮度上限与下限和对应的实验室辐亮度上限与下限；(5. 3)根据步骤(5. 2)得到的各区间在轨辐亮度上限与下限和对应的实验室辐亮度上限与下限，以线性方程为数学模型，解算各区间在轨辐亮度-实验室辐亮度线性方程的系数，所述系数包括增益系数与偏置量系数，所述系数作为各区间在轨绝对辐射定标系数；(5. 4)获取步骤(3)得到的实验室绝对辐射校正影像中各谱段各像元的实验室辐亮度，判断实验室辐亮度在在轨辐亮度-实验室辐亮度查找表中对应的在轨辐亮度-实验 室辐亮度区间，并进而获取所在区间的在轨绝对辐射定标系数，所述系数包括增益系数与偏置量系数；(5. 5)根据步骤(5. 4)得到的包括增益系数与偏置量系数在内的在轨绝对辐射定标系数，采用线性变换方法计算该像元实验室辐亮度对应的在轨辐亮度。本发明与现有技术相比的有益效果是(I)本发明通过实验室绝对辐射校正-在轨绝对辐射校正的逐级辐射校正过程，能够提高遥感卫星多光谱影像上的定量反演精度。本发明所采用的先期进行实验室绝对辐射校正，再进行在轨绝对辐射校正的方法，在实验室绝对辐射校正阶段，由于可设置与量测实验室标准辐亮度等级数量较多、因此能够精确描述影像DN值-辐亮度非线性关系的优势，将影像DN值精确转化为实验室辐亮度，消除非线性特性；在在轨绝对辐射校正阶段，由于影像表达方式转化为实验室辐亮度后以消除了非线性响应特性，因此实验室辐亮度与在轨辐亮度之间呈良好的线性关系，克服了标准在轨辐亮度量测样本数量少的缺陷，能够精确地将实验室辐亮度转化为在轨辐亮度，实现精确定量反演。(2)本发明能够实现遥感卫星多光谱相机谱段间的色彩均衡，通过将多光谱影像的各谱段(色彩)比例关系由DN值间的关系转化为地物实际辐亮度间的关系，从而使多光谱影像的色彩比例真实地反映地物实际光谱特性，减少或消除某一个或某几个谱段间的色彩比例偏差。(3)本发明能够较好地消除遥感卫星多光谱相机的非线性响应特性，传统的绝对辐射定标方法是采用线性方程对所有辐亮度-DN值(或在轨辐亮度-实验室辐亮度)进行线性拟合，解决所有辐亮度-DN值(或在轨辐亮度-实验室辐亮度)统一的绝对定标系数；这种方法对于相机的非线性响应存在很大局限性，在非线性相应部分上述绝对定标系数不能正确反映辐亮度-DN值(或在轨辐亮度-实验室辐亮度)间的关系；通过采用本发明的分段线性插值方法，分段解算辐亮度-DN值(或在轨辐亮度-实验室辐亮度)绝对定标系数，能够使系数解算结果与相机非线性响应部分变化趋势相吻合，消除非线性响应特性。(4)本发明能够扩展遥感卫星多光谱相机的成像动态范围，相机辐亮度-DN值响应的非线性特性主要体现在高辐亮度(高DN值)部分，一般表现为辐亮度增大而DN值随之增大的速率显著降低，从而限制了相机的成像动态范围。通过绝对辐射校正，精确地消除相机响应的非线性特性，能够影像在DN值表达方式下原有的非线性部分转化为线性部分，从而扩展了相机的成像动态范围。
图I为本发明的整体流程图；图2为本发明的详细流程图；图3为本发明的相对辐射校正流程图,对应于步骤⑴；图4为本发明的实验室绝对辐射定标流程图，对应于步骤(2);图5为本发明的基于实验室绝对定标结果的绝对辐射校正流程图，对应于步骤 (3)；图6为本发明的地面检校场在轨绝对辐射定标流程图，对应于步骤(4)；图7为本发明的基于在轨绝对定标结果的绝对辐射校正流程图，对应于步骤(5)；图8为基于本发明步骤(2)所述的实验室绝对辐射定标方法构建的遥感卫星相机不同积分级数、同一积分时间工作参数条件下的实验室辐亮度-影像DN值查找表关系示意图，其中横轴表示影像DN值，纵轴表示实验室辐亮度，不同曲线分别对应不同的相机积分级数，图8表明本发明采用的实验室绝对辐射定标方法可以精确地反映相机辐射响应的非线性特性；图9为基于本发明步骤(4)所述的在轨绝对辐射定标方法构建的遥感卫星相机在轨辐亮度-实验室辐亮度查找表关系示意图，其中横轴表示在轨辐亮度，纵轴表示实验室辐亮度，各个子图分别对应相机的不同工作谱段，图9表明基于本发明所述方法计算得到的实验室辐亮度与在轨辐亮度之间呈良好的线性关系，本发明采用的在轨绝对辐射定标方法可以有效地消除相机辐射响应非线性特性，提高绝对辐射校正精度。
具体实施例方式下面结合图I-图9对本发明的具体实施方式
进行进一步的详细描述，其中图I为本发明的整体流程图，图2为包含每一阶段详细实施步骤的详细流程图，图3-图7分别为各个阶段详细实施步骤的流程图，图8为基于本发明步骤(2)所述的实验室绝对辐射定标方法构建的遥感卫星相机实验室辐亮度-影像DN值查找表关系示意图，图9为基于本发明步骤(4)所述的在轨绝对辐射定标方法构建的遥感卫星相机在轨辐亮度-实验室辐亮度查找表关系不意图(I)相对辐射校正如图3所示基于遥感卫星多光谱相机实验室相对定标系数，采用线性变换方法对多光谱原始影像逐谱段逐像元进行相对辐射校正，消除各CCD探测器单元即探元间的辐射响应不一致性，输出相对辐射校正影像；(I. I)获取相对定标系数获取多光谱相机各谱段各探元的相对定标系数gainrelative(kji)与offsetr—m,其中k表示相机谱段序号，i表示相机每一谱段CCD各探兀序号；(I. 2)相对辐射校正按照公式DNrelativeikj 0 = gainrelative(k, 0 · DNorinikj i)+offsetrelative(kj 0 对多光谱相机各谱段各像元的DN值进行校正，其中DNOTin(k, 0表示第k谱段的第i探元的原始DN值，DNrelativeikj 0表示第k谱段的第i探元的相对辐射校正后DN值。(2)实验室绝对辐射定标如图4所示获取实验室积分球辐射定标源在不同遥感卫星多光谱相机的工作参数下的多级辐亮度与对应DN值，构建并输出实验室辐亮度-影像DN值查找表(不同积分级数的相机实验室辐亮度-影像DN值查找表示意图如图8所示)；(2. I)获取多光谱相机各谱段实验室绝对辐射定标数据，主要包括L1(A)—一定标准辐亮度等级随波长变化的实验室积分球在相机入瞳处的标准辐亮度(单位W · m_2 · Sr-1 · μ πΓ1)，I为积分球标准辐亮度等级序号；Sk(A) 一一随波长变化的各谱段相机光谱响应系数(多光谱相机所有谱段统一进行归一化处理)，k为谱段序号；λ i、λ 2——相机光谱响应范围的起、止波长(单位μ m);
(2. 2)根据定标数据计算实验室积分球各级辐亮度，具体公式如下L{kl) = 矣⑷ ⑷以其中Lfcl)为第k谱段第I辐亮度等级的相机实验室辐射定标标准辐亮度；(2. 3)根据定标数据计算实验室积分球第k谱段第I辐亮度等级的相机响应DN值，通常采用计算积分球某一辐亮度等级对应影像所有有效探元的平均DN值表示，记为DNnreana^1),具体计算公式如下
η
_] DN - S 谓"㈨)
mean(kj)^其中DNi(k, υ为第k谱段第I辐亮度等级第i探元的相机实验室定标成像响应DN值，η为某一谱段CXD探元总数；(2. 4)根据实验室积分球各级辐亮度与响应DN值的关系，构建实验室辐亮度-DN值查找表Llab(k, D-DNmeanik, υ，即形成实验室积分球标准辐亮度向量Llab(k,
D = [L(k,i)，L(k,2)...
L(k,m)]与 DN 值向里 DNmean(k, D — [DNmeann(k, D，DNmean(k,2)…DNmean(k,m)]之间的映射关系，其中m为积分球标准辐亮度等级个数。(3)基于实验室绝对定标结果的绝对辐射校正如图5所示获取步骤(I)得到的相对辐射校正影像成像时刻的遥感卫星多光谱相机的工作参数，基于该工作参数得到对应的实验室辐亮度-影像DN值查找表，通过分段线性插值方法计算各区间对应的实验室绝对辐射定标系数，同时获取相对辐射校正影像中各谱段各像元所在的实验室辐亮度-影像DN值区间，采用所在区间的实验室绝对定标系数对相对辐射校正影像逐谱段逐像元进行实验室绝对辐射校正，输出实验室绝对辐射校正影像；(3. I)获取步骤(I)得到的相对辐射校正影像成像时的遥感卫星多光谱相机工作参数，所述工作参数包含积分级数、增益、积分时间、数字增益和箝位，并根据所述工作参数获取对应的步骤(2)得到的实验室辐亮度-影像DN值查找表Llab(k, D-DNmeanikj D，其中k为谱段序号，I为积分球标准辐亮度等级序号；(3. 2)基于步骤(3. I)得到的实验室辐亮度-影像DN值查找表Llabfcl) DNmean(k，I)，获取各实验室辐亮度-影像DN值区间的实验室辐亮度上限Llab(k,t+1)与下限Llab(k,t)和对应的影像DN值上限DNmeanik, t+1)与下限DNmeanik, t)，其中
Llab (k，I) 一 [Llab (k，I)，Llab (k，2) ··· Llab (k，m)]，DNmean(k，i) [DNmean(k，i)，DNmean(k，2) DNmean(k，m) ]，t为某一辐亮度-DN值区间下限所对应的辐亮度或DN值等级，m为积分球标准辐亮度等级个数；(3. 3)根据步骤(3. 2)得到的各区间实验室辐亮度上限
^lab(k, t+1)与下限 Llab(k，t)
及对应的DN值上限 DNmean(k，t+i) 与下限 DNmean(k，t)， 以线性方程为数学模型，解算各区间辐亮度-DN值线性方程的系数，所述系数包括增益系数gainlab(k, t)与偏置量系数0fTSetlab(k, t)，所述系数作为各区间实验室绝对辐射定标系数，具体公式如下
权利要求
1.一种基于绝对辐射定标的多光谱影像辐射校正方法，其特征在于实现步骤如下 (1)相对辐射校正 基于遥感卫星多光谱相机实验室相对定标系数，采用线性变换方法对多光谱原始影像逐谱段逐像元进行相对辐射校正，消除各CCD探测器单元即探元间的辐射响应不一致性，输出相对辐射校正影像； (2)实验室绝对辐射定标 获取实验室积分球辐射定标源在不同遥感卫星多光谱相机的工作参数下的多级辐亮度与对应DN值，构建并输出实验室辐亮度-影像DN值查找表； (3)基于实验室绝对定标结果的绝对辐射校正 获取步骤(I)得到的相对辐射校正影像成像时刻的遥感卫星多光谱相机的工作参数，基于该工作参数得到对应的实验室辐亮度-影像DN值查找表，通过分段线性插值方法计算各区间对应的实验室绝对辐射定标系数，同时获取相对辐射校正影像中各谱段各像元所在的实验室辐亮度-影像DN值区间，采用所在区间的实验室绝对定标系数对相对辐射校正影像逐谱段逐像元进行实验室绝对辐射校正，输出实验室绝对辐射校正影像； (4)地面检校场在轨绝对辐射定标 获取地面检校场不同反射率靶标的在轨辐射定标野外测量数据并以此计算在轨辐亮度，同时获取靶标所在的采用步骤(3)方法处理得到的实验室绝对辐射校正影像，测量不同反射率靶标在该影像上对应的实验室辐亮度，构建并输出不同反射率靶标的在轨辐亮度-实验室辐亮度查找表； (5)基于在轨绝对定标结果的绝对辐射校正 获取步骤(3)得到的实验室绝对辐射校正影像成像时刻的遥感卫星多光谱相机的工作参数，基于该工作参数得到对应的在轨辐亮度-实验室辐亮度查找表，通过分段线性插值方法计算各区间对应的在轨绝对辐射定标系数，同时获取实验室绝对辐射校正影像中各谱段各像元所在的在轨辐亮度-实验室辐亮度区间，采用所在区间的在轨绝对定标系数对实验室绝对辐射校正影像逐谱段逐像元进行在轨绝对辐射校正，输出在轨绝对辐射校正影像。
2.根据权利要求I所述的基于绝对辐射定标的多光谱影像辐射校正方法，其特征在于所述步骤(3)中具体实现如下 (3. I)获取步骤(I)得到的相对辐射校正影像成像时的遥感卫星多光谱相机工作参数，所述工作参数包含积分级数、增益、积分时间、数字增益和箝位，并根据所述工作参数获取对应的步骤(2)得到的实验室辐亮度-影像DN值查找表； (3. 2)基于步骤(3. I)得到的实验室辐亮度-影像DN值查找表，获取各实验室辐亮度-影像DN值区间的实验室辐亮度上限与下限和对应的影像DN值上限与下限； (3. 3)根据步骤(3. 2)得到的各区间实验室辐亮度上限与下限及对应的DN值上限与下限，以线性方程为数学模型，解算各区间辐亮度-DN值线性方程的系数，所述系数包括增益系数与偏置量系数，所述系数作为各区间实验室绝对辐射定标系数； (3.4)获取步骤(I)得到的相对辐射校正影像中各谱段各像元的DN值，判断其在实验室辐亮度-影像DN值查找表中对应的实验室辐亮度-影像DN值区间，并进而获取所在区间的实验室绝对福射定标系数；(3. 5)根据步骤(3. 4)得到的包括增益系数与偏置量系数在内的实验室绝对辐射定标系数，采用线性变换方法计算该像元DN值对应的实验室辐亮度。
3.根据权利要求I所述的基于绝对辐射定标的多光谱影像辐射校正方法，其特征在于所述步骤(4)中具体实现如下 (4. I)获取地面检校场不同反射率靶标的各谱段在轨辐射定标野外测量数据； (4. 2)根据在轨辐射定标野外测量数据计算地面检校场在轨定标各谱段各级反射率靶标的在轨辐亮度； (4. 3)获取靶标所在的采用步骤(3)方法处理得到的各谱段实验室绝对辐射校正影像； (4. 4)在靶标所在的实验室绝对辐射校正影像上，测量不同反射率靶标对应的实验室辐亮度； (4. 5)根据不同反射率靶标的地面检校场在轨辐亮度与实验室辐亮度的映射关系，构建在轨辐亮度-实验室辐亮度查找表。
4.根据权利要求I所述的基于绝对辐射定标的多光谱影像辐射校正方法，其特征在于所述步骤(5)中具体实现如下 (5. I)获取步骤(3)得到的实验室绝对辐射校正影像成像时的遥感卫星多光谱相机工作参数，所述工作参数包括积分级数、增益、积分时间、数字增益和箝位，并根据所述工作参数获取对应的步骤(4)得到的在轨辐亮度-实验室辐亮度查找表； (5. 2)基于步骤(5. I)在轨辐亮度-实验室辐亮度查找表，获取各在轨辐亮度-实验室辐亮度区间的在轨辐亮度上限与下限和对应的实验室辐亮度上限与下限； (5. 3)根据步骤(5. 2)得到的各区间在轨辐亮度上限与下限和对应的实验室辐亮度上限与下限，以线性方程为数学模型，解算各区间在轨辐亮度-实验室辐亮度线性方程的系数，所述系数包括增益系数与偏置量系数，所述系数作为各区间在轨绝对辐射定标系数； (5.4)获取步骤(3)得到的实验室绝对辐射校正影像中各谱段各像元的实验室辐亮度，判断实验室辐亮度在在轨辐亮度-实验室辐亮度查找表中对应的在轨辐亮度-实验室辐亮度区间，并进而获取所在区间的在轨绝对辐射定标系数，所述系数包括增益系数与偏置量系数； (5. 5)根据步骤(5. 4)得到的包括增益系数与偏置量系数在内的在轨绝对辐射定标系数，采用线性变换方法计算该像元实验室辐亮度对应的在轨辐亮度。
全文摘要
一种基于绝对辐射定标的多光谱遥感影像辐射校正方法，步骤为(1)相对辐射校正；(2)实验室绝对辐射定标构建积分球的实验室辐亮度-影像DN值查找表；(3)基于实验室绝对定标结果的绝对辐射校正基于实验室辐亮度-影像DN值查找表，对相对辐射校正影像进行实验室绝对辐射校正；(4)地面检校场在轨绝对辐射定标构建地面检校场靶标在轨辐亮度-实验室辐亮度查找表；(5)基于在轨绝对定标结果的绝对辐射校正基于在轨辐亮度-实验室辐亮度值查找表，对实验室绝对辐射校正影像进行在轨绝对辐射校正。本发明将多光谱影像由DN值精确地转化为地物实际辐亮度，使多光谱影像的色彩真实地反映地物实际光谱特性，为遥感定量化应用奠定基础。
文档编号G01C25/00GK102901516SQ201210378070
公开日2013年1月30日 申请日期2012年9月29日 优先权日2012年9月29日
发明者张晓 , 李海超, 郝胜勇, 付兴科, 马灵霞, 胡沅 申请人:航天恒星科技有限公司
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
