同时反演在轨高光谱传感器辐射及光谱定标参数的方法
【技术领域】
[0001]本发明涉及遥感传感器技术领域，具体而言，特别涉及同时反演在轨高光谱传感器辐射及光谱定标参数的方法。
【背景技术】
[0002]在轨传感器的性能指标评价，主要用于对在轨运行的高光谱传感器的性能指标进行监测，包括辐射定标参数(增益和偏移)和光谱定标参数(中心波长和波段宽度)，数据生产商或用户根据定标的结果进行遥感数据生产和处理，保证遥感数据本身的准确性和可靠性。
[0003]传感器平台(星载或空载)的抖动、仪器老化及周围环境的变化使得在轨传感器辐射及光谱定标参数较发射前的实验室测量的初始值会发生变化，如不及时纠正，则会影响所接收的遥感图像的质量和准确性，只有准确定标后的数据才能被正确使用，用户才能从图像上获得可靠的信息。
[0004]对传感器定标而言，辐射定标和光谱定标二者理论上应该同步进行，因为辐射地表的参数在光谱定标的时候需要，而光谱定标的参数值辐射定标的时候也需要，二者相互影响，需要同步反演。
[0005]目前，已有的技术辐射定标相对比较成熟，通常利用地面靶标进行在轨传感器增益和偏移参数的反演，基于此可以将遥感图像的数字值(DN，Digital Number)转换成具有物理意义的辐射亮度；而光谱定标则更多的研宄只反演中心波长的偏移，不考虑波段宽度。此外，在进行辐射定标和光谱定标时，二者定标过程是相互独立的，独立定标将导致所得的福射定标参数和光谱定标参数都不准确，影响数据质量。
[0006]针对现有技术中定标方法得到的辐射定标参数和光谱定标参数准确性较低的问题，目前尚未提出有效的解决方法。
【发明内容】
[0007]本发明的主要目的在于提供一种同时反演在轨高光谱传感器辐射及光谱定标参数的方法，以解决现有技术中定标方法得到的辐射定标参数和光谱定标参数准确性较低的冋题。
[0008]依据本发明的一个方面，提供了一种同时反演在轨高光谱传感器辐射及光谱定标参数的方法。该方法包括:通过大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度；根据发射前实验室测量的传感器的光谱定标参数计算传感器的入瞳辐射亮度LI ;获取待反演定标参数的高光谱DN值图像；根据发射前实验室测量的传感器的辐射定标参数和高光谱DN值图像计算高光谱DN值图像的辐射亮度L2 ;对LI和L2分别进行导数计算，依次得到第一求导结果和第二求导结果；对第一求导结果分别进行归一化处理、包络线去除处理和计算光谱角处理，依次得到第一归一化结果al、第一包络线去除结果bl和第一计算光谱角结果Cl ;对第二求导结果分别进行归一化处理、包络线去除处理和计算光谱角处理，依次得到第二归一化结果a2、第二包络线去除结果b2和第二计算光谱角结果c2 ;以al、bl、cl、a2、b2、c2为基础，采用优化算法迭代对比LI和L2的差别，迭代过程中优化算法自动调整待反演参数，直到LI和L2的差别满足预设条件，此时的待反演参数即为反演结果，其中，待反演参数包括增益、偏移、中心波长和波段宽度。
[0009]进一步地，通过大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度具体为:以预设的大气温度、湿度、气溶胶能见度以及在轨传感器的角度和太阳的角度为驱动数据，驱动大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度。
[0010]进一步地，在驱动大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度时，地表的反射率设置为0.1-0.3之间的常数。
[0011]进一步地，大气辐射传输模型为MODTRAN5大气辐射传输模型。
[0012]进一步地，根据发射前实验室测量的传感器的光谱定标参数计算传感器的入瞳辐射亮度LI具体为:以发射前实验室测量的传感器的光谱定标参数为初值，对模拟得到的I纳米分辨率的辐射亮度采用高斯函数进行卷积计算，得到传感器的入瞳辐射亮度LI。
[0013]进一步地，根据发射前辐射定标参数和高光谱DN值图像计算对应的辐射亮度L2具体为:利用发射前实验室测量的传感器的辐射定标参数对高光谱DN值图像进行辐射定标，得到高光谱DN值图像的辐射亮度L2。
[0014]进一步地，预设条件为由al与a2的均方根误差、bl与b2的均方根误差、cl与c2的均方根误差构成的三个均方根误差之和最小。
[0015]通过本发明，提供了一种同时反演在轨高光谱传感器辐射及光谱定标参数的方法，在该方法中，通过大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度，然后根据发射前实验室测量的传感器的光谱定标参数、辐射定标参数分别计算传感器的入瞳辐射亮度LI和待反演定标参数的高光谱图像的辐射亮度L2，再对LI和L2分别进行导数计算后归一化处理、包络线去除处理和计算光谱角处理，最后对处理结果利用优化算法迭代计算，当LI和L2的差别达到最小时即为反演结果，也即增益、偏移、中心波长和波段宽度，从而实现在轨传感器辐射定标参数(增益和偏移)和光谱定标参数(中心波长和波段宽度)的同步反演，弥补已有技术只能确定部分参数的缺陷，为生产高质量遥感数据提供了技术保障。
[0016]上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的【具体实施方式】。
【附图说明】
[0017]通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。在附图中:
[0018]图1是根据本发明第一实施例的同时反演在轨高光谱传感器辐射及光谱定标参数的方法的流程图；
[0019]图2是根据本发明第二实施例的同时反演在轨高光谱传感器辐射及光谱定标参数的方法的流程图。
【具体实施方式】
[0020]下面结合附图和【具体实施方式】对本发明做进一步说明。需要指出的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
[0021]实施例一
[0022]首先，本发明实施例一提供了一种同时反演在轨高光谱传感器辐射及光谱定标参数的方法，参见图1，该方法可以包括以下步骤:
[0023]步骤S102:通过大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度。
[0024]其中，可以预设的大气温度、湿度、气溶胶能见度以及在轨传感器的角度和太阳的角度为驱动数据，驱动MODTRAN5大气辐射传输模型模拟大气层顶部I纳米分辨率的辐射亮度，模拟过程中，地表的反射率设置为0.1-0.3之间的常数。
[0025]步骤S104:根据发射前实验室测量的传感器的光谱定标参数计算传感器的入瞳辐射亮度LI。
[0026]具体地，以发射前实验室测量的传感器的光谱定标参数为初值，对模拟得到的辐射亮度采用高斯函数进行卷积计算，得到传感器的入瞳辐射亮度LI。
[0027]步骤S106:获取待反演定标参数的高光谱DN值图像。
[0028]步骤S108:根据发射前实验室测量的传感器的辐射定标参数和高光谱DN值图像计算高光谱DN值图像的辐射亮度L2。
[0029]具体地，利用发射前实验室测量的传感器的辐射定标参数对高光谱DN值图像进行辐射定标，得到高光谱DN值图像的辐射亮度L2。
[0030]步骤SllO:对LI和L2分别进行导数计算，依次得到第一求导结果和第二求导结果O
[0031]步骤S112:对第一求导结果分别进行归一化处理、包络线去除处理和计算光谱角处理，依次得到第一归一化结果al、第一包络线去除结果bl和第一计算光谱角结果Cl。
[0032]步骤S114:对第二求导结果分别进行归一化处理、包络线去除处理和计算光谱角处