理，依次得到第二归一化结果a2、第二包络线去除结果b2和第二计算光谱角结果c2。
[0033]步骤S116:以al、bl、cl、a2、b2、c2为基础，采用优化算法迭代对比LI和L2的差另IJ，迭代过程中优化算法自动调整待反演参数，直到LI和L2的差别满足预设条件，此时的待反演参数即为反演结果。
[0034]其中，待反演参数包括增益、偏移、中心波长和波段宽度，预设条件为由al与a2的均方根误差、bl与b2的均方根误差、Cl与c2的均方根误差构成的三个均方根误差之和最小。
[0035]采用该实施例，能够实现在轨高光谱传感器辐射定标参数(增益和偏移)和光谱定标参数(中心波长和波段宽度)的同步反演，弥补已有技术只能确定部分参数的缺陷，为生产高质量遥感数据提供了技术保障。
[0036]实施例二
[0037]该实施例为在实施例一的基础上进一步优选的同时反演在轨高光谱传感器辐射及光谱定标参数的方法，参见图2，该方法可以包括以下步骤:
[0038]高光谱传感器由于波段宽度较窄对中心波长与波段宽度的变化比较敏感，尤其在大气吸收波段周围光谱定标参数的变化会在光谱辐射亮度中明显体现出来，这为利用吸收波段反演高光谱辐射和光谱定标参数提供了可能。本技术要同时反演的参数共4个，S卩，辐射定标参数:增益和偏移，光谱定标参数:中心波长和波段宽度。
[0039](I)设定典型大气温度、湿度、气溶胶能见度以及传感器和太阳的角度驱动MODTRAN5大气辐射传输模型模拟大气层顶部高分辨(I纳米)辐射亮度。具体模拟过程中地表的反射率设置为0.1-0.3之间的常数；
[0040](2)以先前实验室测量的传感器光谱定标参数为初值，采用高斯函数模拟得到传感器入瞳辐射亮度LI ；
[0041](3)利用先前实验室的辐射定标参数对待反演定标参数的高光谱DN值图像进行定标，得到辐射亮度L2 ;
[0042](4)将LI和L2分别进行导数计算并归一化、包络线去除以及计算光谱角三种处理，以处理后的数据为基础，采用优化算法迭代对比LI和L2的差别，迭代过程中优化算法自动调整待反演参数(增益、偏移、中心波长和波段宽度)，直到LI和L2之间的均方根误差最小(上述三种处理方式的均方根误差之和最小)，此时所得到增益、偏移、中心波长和波段宽度4个参数即为反演结果。
[0043]采用该实施例，首次实现在轨高光谱传感器增益、偏移、中心波长以及波段宽度共4个性能参数的同步反演技术；同时，该技术操作简单，不受地表类型限制，反演过程也无需实时的地面靶标测量数据，仅依靠遥感图像即可实现在轨传感器定标参数监测。
[0044]从以上各实施例的描述中，可以看出，本发明实施例实现了如下技术效果:同时实现在轨高光谱传感器辐射定标参数(增益和偏移)和光谱定标参数(中心波长和波段宽度)的反演，监测传感器的性能变化，以便根据性能参数的变化调整遥感数据处理方法，为生产质量可靠的遥感数据提供技术保障。
[0045]需要说明的是，上述装置或系统实施例属于优选实施例，所涉及的单元和模块并不一定是本申请所必须的。
[0046]本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。对于本申请的装置实施例而言，由于其与方法实施例基本相似，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
[0047]以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此，任何熟悉该技术的人在本发明所揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
【主权项】
1.一种同时反演在轨高光谱传感器辐射及光谱定标参数的方法，其特征在于，包括: 通过大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度； 根据发射前实验室测量的所述传感器的光谱定标参数计算所述传感器的入瞳辐射亮度LI ; 获取待反演定标参数的高光谱DN值图像； 根据发射前实验室测量的所述传感器的辐射定标参数和所述高光谱DN值图像计算所述高光谱DN值图像的辐射亮度L2 ； 对LI和L2分别进行导数计算，依次得到第一求导结果和第二求导结果； 对所述第一求导结果分别进行归一化处理、包络线去除处理和计算光谱角处理，依次得到第一归一化结果al、第一包络线去除结果bl和第一计算光谱角结果Cl ; 对所述第二求导结果分别进行归一化处理、包络线去除处理和计算光谱角处理，依次得到第二归一化结果a2、第二包络线去除结果b2和第二计算光谱角结果c2 ; 以al、bl、cl、a2、b2、c2为基础，采用优化算法迭代对比LI和L2的差别，迭代过程中优化算法自动调整待反演参数，直到LI和L2的差别满足预设条件，此时的待反演参数即为反演结果，其中，所述待反演参数包括增益、偏移、中心波长和波段宽度。2.根据权利要求1所述的同时反演在轨高光谱传感器辐射及光谱定标参数的方法，其特征在于，通过大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度具体为: 以预设的大气温度、湿度、气溶胶能见度以及在轨传感器的角度和太阳的角度为驱动数据，驱动大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度。3.根据权利要求2所述的同时反演在高光谱轨传感器辐射及光谱定标参数的方法，其特征在于，在驱动大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度时，地表的反射率设置为0.1-0.3之间的常数。4.根据权利要求3所述的同时反演在轨高光谱传感器辐射及光谱定标参数的方法，其特征在于，所述大气辐射传输模型为MODTRAN5大气辐射传输模型。5.根据权利要求1所述的同时反演在轨高光谱传感器辐射及光谱定标参数的方法，其特征在于，根据发射前实验室测量的所述传感器的光谱定标参数计算所述传感器的入瞳辐射亮度LI具体为: 以发射前实验室测量的所述传感器的光谱定标参数为初值，对模拟得到的I纳米分辨率的辐射亮度采用高斯函数进行卷积计算，得到所述传感器的入瞳辐射亮度LI。6.根据权利要求1所述的同时反演在轨高光谱传感器辐射及光谱定标参数的方法，其特征在于，根据发射前实验室测量的所述传感器的辐射定标参数和所述高光谱DN值图像计算所述高光谱DN值图像的辐射亮度L2具体为: 利用发射前实验室测量的所述传感器的辐射定标参数对所述高光谱DN值图像进行辐射定标，得到所述高光谱DN值图像的辐射亮度L2。7.根据权利要求1所述的同时反演在轨高光谱传感器辐射及光谱定标参数的方法，其特征在于，所述预设条件为由al与a2的均方根误差、bl与b2的均方根误差、cl与c2的均方根误差构成的三个均方根误差之和最小。
【专利摘要】本发明公开了一种同时反演在轨高光谱传感器辐射及光谱定标参数的方法。该方法包括：通过大气辐射传输模型模拟大气层顶部1纳米分辨率的辐射亮度；以发射前光谱定标参数为初值,通过优化算法不断调整光谱定标参数,并计算对应的传感器入瞳辐射亮度L1；获取待反演定标参数的高光谱DN值图像；根据发射前辐射定标参数和高光谱DN值图像计算对应的辐射亮度L2；对L1和L2分别进行导数计算，并对两个求导结果分别进行归一化处理、包络线去除处理和计算光谱角处理，以处理结果为基础，采用优化算法迭代对比L1和L2的差别，直到优化满足预设条件，此时的待反演参数即为反演结果。通过本发明，能够实现在轨高光谱传感器辐射定标参数和光谱定标参数的同步反演。
【IPC分类】G01S7/497
【公开号】CN104914424
【申请号】CN201510239765
【发明人】王天星, 阎广建, 穆西晗, 任华忠, 刘耀开
【申请人】中国科学院遥感与数字地球研究所
【公开日】2015年9月16日
【申请日】2015年5月12日