一种星星间交叉辐射定标的匹配阈值分析方法

技术领域

本发明属于交叉辐射定标技术领域，具体是一种星星间交叉辐射定标的匹配阈值分析方法。

背景技术

交叉辐射定标,即建立待定标传感器通道与参考传感器通道输出之间的转换关系，是在轨替代定标的一种有效手段，可以充分利用多种数据源实现不同传感器数据的相互转换与同化应用，即利用定标精度较高的遥感仪器作为参考基准，选择相同或相近观测条件下的资料，建立参考遥感仪器入瞳辐亮度与目标遥感仪器测量值之间的关系，实现对目标遥感仪器的定标，交叉定标在我国已经走过了将近二十年的发展历程，从早期的针对不同遥感器的单次交叉定标试验，到目前建立交叉定标业务系统，交叉辐射定标技术已经在我国气象卫星、资源卫星、环境卫星、海洋卫星、军事卫星等领域得到了广泛应用；

然而，由于星星间交叉定标匹配阈值设置不合理，导致匹配样本数量过大、匹配样本代表性不强、匹配样本无法满足交叉定标计算需求，直接影响了星星间交叉定标的精度。

发明内容

本发明的目的在于克服现有技术的缺陷，提供一种星星间交叉辐射定标的匹配阈值分析方法。

为实现上述目的，本发明采用了如下技术方案：

一种星星间交叉辐射定标的匹配阈值分析方法，包括以下步骤：

S1观测文件匹配：基于文件匹配规则，对目标遥感器和参考遥感器的观测数据配对，获得时间匹配的观测文件对；

S2像元时空匹配：基于时间匹配的观测文件对，利用像元匹配规则，实现观测像元的时间、空间和观测角度的匹配，获得像元匹配样本；

S3数据转换：基于光谱响应函数，对像元匹配样本进行数据转换，获得单位统一可直接对比的样本；

S4样本过滤：基于样本过滤阈值，对样本进行过滤，获得可用于交叉定标分析的样本；

S5样本累积：基于卫星自身的稳定性和交叉匹配点的时间分布特征，对多天定标分析样本进行累积，获得若干组样本；

S6质量判断：获取符合质量判断标准的匹配样本。

优选的，所述S1观测文件匹配，包括：

选择目标遥感器星下点经纬度(±DLon，±DLat)以内矩形区域为固定目标区域，获取参考遥感器经过固定目标区域的观测文件，其中，DLat和DLon不超过35°；

结合两卫星扫描规律，将目标遥感器和参考遥感器观测时间最接近的文件进行匹配；

选择目标遥感器与参考遥感器在若干时间阈值、若干空间阈值内的近同时星下点观测文件，作为匹配目标文件，其中，时间阈值小于20分钟，空间阈值小于1000公里。

优选的，所述S2像元匹配规则，包括：像元观测时间匹配规则、空间匹配规则及观测角度匹配规则；

所述观测时间匹配规则，需满足公式：

|t

其中，t

所述像元空间匹配规则，包括空间位置匹配规则及空间距离筛选规则；

所述空间位置匹配规则：选择参考遥感器像元(x

所述空间距离筛选规则：对于空间位置匹配的像元对(x

其中，x为观测像元中心经度，y为观测像元中心纬度，D

所述观测角度匹配规则，需满足公式：

其中，θ

优选的，所述S 3数据转换，包括：

辐射单位转换：将像元匹配样本的辐射单位均转换为毫瓦每平方米每球面度每波数；

光谱匹配：基于参考遥感器观测数据，计算目标遥感器的通道模拟辐射，计算公式为：

其中，L

优选的，所述S 4样本过滤，包括：

空间均匀性过滤：基于对环境场区和等效观测视场均匀性的双重检测，获得均匀情景下的观测结果，环境场区均匀性检测公式为：

SD

其中，SD

等效观测视场均匀性检测公式如下：

|E

其中，E

异常点剔除：基于通道辐亮度的有效物理范围，剔除异常观测样本，获得可用于交叉定标分析的样本。

优选的，所述S 6质量判断，基于质量判断标准，不满足质量判断标准的匹配样本重复步骤S1至步骤S 5，并按照规定步长调整匹配阈值，直至获取符合质量判断标准的匹配样本。

优选的，所述质量判断标准，包括：观测样本累积周期按交叉点出现周期进行统计，样本数大于100，且样本之间线性相关系数大于0.98，目标遥感器搭载于静止轨道卫星平台统计周期为1个月，目标遥感器搭载于极轨卫星平台统计周期不超过7天。

综上所述，由于采用了上述技术方案，本发明的有益效果是：

本发明通过分析星星间交叉辐射定标获取的匹配样本，对关键要素匹配阈值进行迭代和调整，通过匹配样本质量判断来确定合理的阈值，从而根据分析结果，针对不同类型的参考遥感器和目标遥感器的轨道特征、光谱通道、观测模式等载荷平台特性，给出合理的匹配阈值区间和推荐阈值列表；

本发明提出的星星间交叉辐射定标匹配阈值分析方法针对参考遥感器和目标遥感器的特点给出合理阈值范围和阈值列表，可以极大的提高交叉辐射定标匹配样本质量，从而提高定标精度，同时还可以有效节省交叉辐射定标匹配计算的时间，提高定标效率。

附图说明

图1是本发明一种星星间交叉辐射定标的匹配阈值分析方法的流程图。

具体实施方式

以下结合附图1，进一步说明本发明一种星星间交叉辐射定标的匹配阈值分析方法的具体实施方式。本发明一种星星间交叉辐射定标的匹配阈值分析方法不限于以下实施例的描述。

实施例1：

本实施例给出一种星星间交叉辐射定标的匹配阈值分析方法的具体结构，如图1所示，包括以下步骤：

S1观测文件匹配：基于文件匹配规则，对目标遥感器和参考遥感器的观测数据配对，获得时间匹配的观测文件对；

S2像元时空匹配：基于时间匹配的观测文件对，利用像元匹配规则，实现观测像元的时间、空间和观测角度的匹配，获得像元匹配样本；

S3数据转换：基于光谱响应函数，对像元匹配样本进行数据转换，获得单位统一可直接对比的样本；

S4样本过滤：基于样本过滤阈值，对样本进行过滤，获得可用于交叉定标分析的样本；

S5样本累积：基于卫星自身的稳定性和交叉匹配点的时间分布特征，对多天定标分析样本进行累积，获得若干组样本；

S6质量判断：获取符合质量判断标准的匹配样本。

具体地，S1观测文件匹配，包括：

选择目标遥感器星下点经纬度(±DLon，±DLat)以内矩形区域为固定目标区域，获取参考遥感器经过固定目标区域的观测文件，其中，DLat和DLon不超过35°；

结合两卫星扫描规律，将目标遥感器和参考遥感器观测时间最接近的文件进行匹配；

选择目标遥感器与参考遥感器在若干时间阈值、若干空间阈值内的近同时星下点观测文件，作为匹配目标文件，其中，时间阈值小于20分钟，空间阈值小于1000公里。

具体地，像元匹配规则，包括：像元观测时间匹配规则、空间匹配规则及观测角度匹配规则；

观测时间匹配规则，需满足公式：

|t

其中，t

像元空间匹配规则，包括空间位置匹配规则及空间距离筛选规则；

空间位置匹配规则：选择参考遥感器像元(x

空间距离筛选规则：对于空间位置匹配的像元对(x

其中，x为观测像元中心经度，y为观测像元中心纬度，D

观测角度匹配规则，需满足公式：

其中，θ

具体地，S 3数据转换，包括：

辐射单位转换：将像元匹配样本的辐射单位均转换为毫瓦每平方米每球面度每波数；

光谱匹配：基于参考遥感器观测数据，计算目标遥感器的通道模拟辐射，计算公式为：

其中，L

进一步的，S 4样本过滤，包括：

空间均匀性过滤：基于对环境场区和等效观测视场均匀性的双重检测，获得均匀情景下的观测结果，环境场区均匀性检测公式为：

SD

其中，SD

等效观测视场均匀性检测公式如下：

|E

其中，E

异常点剔除：基于通道辐亮度的有效物理范围，剔除异常观测样本，获得可用于交叉定标分析的样本。

进一步的，S 6质量判断，基于质量判断标准，不满足质量判断标准的匹配样本重复步骤S1至步骤S 5，并按照规定步长调整匹配阈值，直至获取符合质量判断标准的匹配样本。

进一步的，质量判断标准，包括：观测样本累积周期按交叉点出现周期进行统计，样本数大于100，且样本之间线性相关系数大于0.98，目标遥感器搭载于静止轨道卫星平台统计周期为1个月，目标遥感器搭载于极轨卫星平台统计周期不超过7天。

工作原理：如图1所示，首先，基于文件匹配规则，对目标遥感器和参考遥感器的观测数据配对，获得时间匹配的观测文件对；

其次，基于时间匹配的观测文件对，利用像元匹配规则，实现观测像元的时间、空间和观测角度的匹配，获得像元匹配样本；

接着，基于光谱响应函数，对像元匹配样本进行数据转换，获得单位统一可直接对比的样本；

然后，基于样本过滤阈值，对样本进行过滤，获得可用于交叉定标分析的样本；

随后，基于卫星自身的稳定性和交叉匹配点的时间分布特征，对多天定标分析样本进行累积，获得若干组样本；

最后，获取符合质量判断标准的匹配样本；

本发明通过分析星星间交叉辐射定标获取的匹配样本，对关键要素匹配阈值进行迭代和调整，通过匹配样本质量判断来确定合理的阈值，从而根据分析结果，针对不同类型的参考遥感器和目标遥感器的轨道特征、光谱通道、观测模式等载荷平台特性，给出合理的匹配阈值区间和推荐阈值列表；

本发明提出的星星间交叉辐射定标匹配阈值分析方法针对参考遥感器和目标遥感器的特点给出合理阈值范围和阈值列表，可以极大的提高交叉辐射定标匹配样本质量，从而提高定标精度，同时还可以有效节省交叉辐射定标匹配计算的时间，提高定标效率。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。