[发明专利]一种基于统计样本定权的光学遥感卫星相对辐射校正方法

说明书

本申请公开了一种基于统计样本定权的光学遥感卫星相对辐射校正方法，该方法包括：根据预设起始时间、预设结束时间以及预设统计周期采集0级数据作为第N月对应的统计样本，根据统计样本进行直方图统计和匹配得到第一灰度查找表，并根据第一灰度查找表对第N月对应的遥感卫星影像进行校正得到初始校正后的影像；判断初始校正后的影像中搭接区是否存在条带噪声；若存在，则调整预设起始时间、结束时间和统计周期；根据调整后的起始时间、结束时间和统计周期对第N月对应的遥感卫星影像进行校正得到校正后的影像。本申请解决了现有技术中对探测器状态发生跳变时所采集的影像的校正效果较差的技术问题。

技术领域

本申请涉及航天光学遥感影像处理技术领域，尤其涉及一种基于统计样本定权的光学遥感卫星相对辐射校正方法。

背景技术

光学遥感卫星相机一般采用多个全色多光谱TDICCD进行反射镜拼接，而CCD反射镜拼接是利用分光反射镜将像平面分割成空间分离的多个像面，可在各个像平面交错布置多个CCD，并使交错相邻两片CCD边界处像元满足重叠要求，形成等效的大视场探测器。但是，采用反射镜拼接方法在搭接区会产生渐晕现象，使得搭接区的辐射能量呈逐渐下降趋势，像元之间的辐射响应一致性显著下降，形成搭接区条带噪声，进而影响遥感卫星影像。因此，为了避免搭接区条带噪声对光学遥感卫星影像的影响，需要对遥感卫星进行辐射校正。

目前，常见了解决搭接区条带噪声的方法包括查表法、逐行扫描法和函数逼近法等方法。参见申请号2017106446289《一种基于大数据统计的遥感卫星相对辐射定标处理方法》的专利介绍了一种基于查表法对0级数据进行相对辐射校正的方法，该方法的具体过程为：采用遥感卫星0及条带数据全统计的方法，对数据进行直方图统计，然后根据直方图统计，进行概率的解算，并对DN值为0和DN值为最大量化值进行边界处理，再对直方图和统计概率进行匹配，形成查找表；根据查找表对0级数据进行相对辐射校正。但是，该专利中所介绍的方法是统计一段时间内每个探元的归一化直方图，将所有非搭接区的探元直方图加和并归一化作为期望直方图。搭接区和非搭接区的每个探元的直方图均与期望直方图进行匹配，生成查找表。统计周期一般为两个月。用第N-1月上半月至第N月下半月两个月的0级数据作为统计样本进行直方图统计和匹配，计算处的灰度查找表用来校正第N月上半月的数据；用第N-1月下半月至第N+1月上半月2个月的0级数据作为统计样本进行直方图统计和匹配，计算出的灰度查找表用来校正第N月下半月的数据；以此类推。对于探测器成像状态缓慢变化，不存在震荡性跳变的情况下，该方法校正效果良好。但当探测器成像状态在统计周期时间内出现一定频率的跳变时，生成的查找表对这些探测器状态发生跳变时采集的影像的校正效果较差。

发明内容

本申请解决的技术问题是：针对现有技术中对探测器状态发生跳变时所采集的影像的校正效果较差。本申请提供了一种基于统计样本定权的光学遥感卫星相对辐射校正方法，本申请实施例所提供的方案中，当根据预设起始时间、预设结束时间和预设统计周期对第N月对应的遥感卫星影像进行辐射校正所得到初始校正后的影像中搭接区存在条带噪声时，通过调整预设起始时间、预设结束时间和预设统计周期，直到调整后的起始时间、结束时间和统计周期满足预设需求，然后再根据调整后的起始时间、结束时间和统计周期对第N月对应的遥感卫星影像进行校正得到校正后的影像，即本申请实施例所提供的方案中，通过校正后的影像中搭接区信息来动态调整统计的预设起始时间、预设结束时间和预设统计周期，避免当探测器成像状态在统计周期时间内出现一定频率的跳变时，生成的查找表对这些探测器状态发生跳变时采集的影像的校正效果较差的问题。

第一方面，本申请实施例提供一种基于统计样本定权的光学遥感卫星相对辐射校正方法，该方法包括：

根据预设起始时间、预设结束时间以及预设统计周期采集0级数据作为第N月对应的统计样本，根据所述统计样本进行直方图统计和匹配得到第一灰度查找表，并根据所述第一灰度查找表对第N月对应的遥感卫星影像进行校正得到初始校正后的影像，其中，N[1,12]；

判断所述初始校正后的影像中搭接区是否存在条带噪声；