[发明专利]一种光学卫星相对辐射校正系数计算方法、系统及介质

说明书

一种光学卫星相对辐射校正系数计算方法，包括如下步骤：步骤一、对数据源进行第一次挑选，获得第一次挑选后的数据的文件属性索引；步骤二、对第一次挑选后的数据进行分类，获得分类后的数据；步骤三、将步分类后的数据按投入计算节点的数量平均分配，将数据列表和各计算节点必要的执行信息写入到标准化任务清单中；步骤四、初始计算节点将标准化任务清单依次发送到其他各子计算节点上；各子计算节点生成标准化统计数据文件，然后发送回初始计算节点；步骤五、利用直方图匹配方法，获得相对辐射校正系数。本发明方法由在轨卫星的多片CCD综合得到标准的灰度概率密度直方图，能够兼顾单片CCD自身以及多片CCD片间的相对辐射一致性。

技术领域

本发明涉及一种光学卫星相对辐射校正系数计算方法、系统及介质，属于遥感影像处理领域。

背景技术

遥感影像相对辐射校正是一项基础的数据处理工作，是提高遥感卫星图像质量的一个极其重要的途径。相对辐射校正方法主要分为定标法和统计法，相较于定标法，统计法不需要特殊的训练数据，通过分析和统计直接从图像中提取校正系数，从而对影像进行校正，因此当前光学卫星影像地面数据处理中心所分发的0级影像数据适合使用统计法进行相对辐射校正。

采用直方图修整可使图像的灰度间距拉大或使灰度分布均匀，从而使图像的细节更加清晰，达到增强的目的。直方图匹配(直方图规定化)是通过建立直方图查找表，将原始图像的直方图匹配到期望直方图上。一般来说，直方图匹配法要求所选取的图像尽可能多的覆盖各类地物，样本数量充足，使得各个探元的灰度分布比较一致。但是这种在实际应用中也存在一定的局限性。首先，在轨卫星一般有多个成像用的CCD，各CCD的响应系数很难保持一致，使用单片CCD统计得到的标准直方图不能很好的顾及整个影像；在轨卫星影像数据每日都在生产，要覆盖足够多的典型区域会使得统计的数据量极大，一般的单机统计处理耗时长，难以保证相对辐射校正系数的及时更新，无法满足实际生产的需要；此外，依靠人工选取典型影像耗时长、可靠性差，且无法实现整个相对辐射校正系数生产的自动化。

发明内容

本发明要解决的技术问题是：克服现有技术的不足，提供了一种光学卫星相对辐射校正系数计算方法及系统，综合多个成像CCD的成像单元进行统计并获取标准直方图，能够根据特定条件对原始数据进行筛选、分类，采用多节点分布式计算获取相对辐射校正系数。

本发明目的通过以下技术方案予以实现：

一种光学卫星相对辐射校正系数计算方法，包括如下步骤：

步骤一、根据卫星型号、传感器类型、成像模式、成像时间段，对数据源进行第一次挑选，获得第一次挑选后的数据的文件属性索引；

步骤二、基于步骤一中所述数据的文件属性索引中的增益、级数信息，对第一次挑选后的数据进行分类，获得分类后的数据；

步骤三、将步骤二中所述分类后的数据按投入计算节点的数量平均分配，将数据列表和各计算节点必要的执行信息写入到标准化任务清单中；

步骤四、初始计算节点将步骤三中所述的标准化任务清单依次发送到其他各子计算节点上；各子计算节点执行相应的任务清单，生成标准化统计数据文件，然后发送回初始计算节点；

步骤五、初始计算节点接收各子计算节点的标准化统计数据文件后，进行加和统计，获得标准灰度概率密度直方图；然后利用直方图匹配方法，获得相对辐射校正系数。

上述光学卫星相对辐射校正系数计算方法，步骤二中所述分类后的数据包括目标增益和级数组合数据、其他增益和级数组合数据。

上述光学卫星相对辐射校正系数计算方法，步骤三中所述各计算节点必要的执行信息包括卫星型号、传感器类型、成像模式、成像时间段、单景输出路径，综合输出路径。