(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111557158.5 (22)申请日 2021.12.18 (71)申请人 国家卫星气象中心(国家空间天气 监测预警中心) 地址 100044 北京市海淀区中关村南大街 46号 (72)发明人 孙凌　张阆　王春晓　胡秀清　 张鹏　 (74)专利代理 机构 北京知企鸿蒙专利代理事务 所(普通合伙) 11692 代理人 张瑞玲 (51)Int.Cl. G01S 17/95(2006.01) G01N 21/59(2006.01) G01N 21/55(2014.01)G06F 30/20(2020.01) (54)发明名称 一种用于反射太阳波段的多场地定标跟踪 方法 (57)摘要 本发明公开了一种用于反射太阳波段的多 场地定标跟踪方法， 包括如下步骤： S1： 数据预处 理； S2： 辐射传输模拟； S3： 定标校验计算； S4： 响 应衰减计算， 此用于反射太阳波段的多场地定标 跟踪方法， 区别于现有技术， 提高了定标校验频 次， 可以满足有效监测仪器 响应衰减和在轨定标 系数更新的需求， 且该方法以稳定目标物的辐射 传输模拟值作为基准， 通过比对模拟值与观测值 计算定标校验系数， 并基于定标校验系数的长期 趋势分析建立了 定标系数 更新模型。 权利要求书2页 说明书5页 附图3页 CN 114265085 A 2022.04.01 CN 114265085 A 1.一种用于反射太阳波段的多场地定标跟踪方法， 其特 征在于， 包括如下步骤： S1： 数据预处 理； S2： 辐射传输模拟； S3： 定标校验计算； S4： 响应衰减计算。 2.根据权利要求1所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 在S1数据预处理时， 读取所选场地周边n ×n个像素点范围的DN值进行计算， 并以选定通道 DN值标准差 /平均值为基准， 剔除值大于 0.02的数据。 3.根据权利要求1所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 在S1数据预处 理时， 对于反射波段， 计算耀斑角并剔除耀斑角＜40 °的观测数据。 4.根据权利要求1所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 在S1数据预处 理时， 将太阳天顶角＞70 °和观测天顶角＞6 0°的观测数据进行剔除处 理。 5.根据权利要求1所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 在S2辐射传输模拟时， 使用矢量辐射传输模型6SV计算晴空目标区的光谱辐射。 6.根据权利要求5所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 在6SV计算时， 采用MODTRAN计算的大气透过率对6SV的结果进行大气层顶反射率的吸收订 正。 7.根据权利要求6所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 对6SV计算模拟结果进行日地距离校正和太阳天顶角余弦校正。 8.根据权利要求7所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 在S3定标校验计算时， 利用6SV联合MODTRAN模 型计算得到的大气层顶表观反射率与卫星观 测经定标处 理得到的反射 率计算定标 校验系数： 其中 为模拟的大气层顶表观反射率， 为观测的经定标处理得到的等效到 日地平均距离和太阳位于天顶时的表观反射率， i表示第i通道， SolZ为太阳天顶角， d为日 地天文单位距离， Sl opei为第i通道的定标 校验系数； 并利用所有场地数据拟合得到给定 定标周期内的定标 校验系数， 并累积其时间序列。 基于定标校验系数时间序列， 采用多项式模型， 如二 次多项式建立定标校正 时变模型， 描述定标 校验系数的长期变化趋势： Slopei＝Fun1(D SL)＝a1iDSL+a2iDSL2+a3i 其中， DSL为自发星日起的日计数， a为定标 校正时变模型系数。 9.根据权利要求8所述的一种用于反射太阳波段的多场 地定标跟踪方法， 其特征在于： 在S4响应衰减计算时， 若 为第i个通道若为观测的经发射前定标处理得到的反射率， 则可定义第i个通道的归一化辐射响应为NResponsei＝a3i/Slopei,描述遥感器辐射响应的 衰变表示 为： NResponsei＝Fun2(D SL)＝c1iDSL+c2iDSL2+c3i权　利　要　求　书 1/2 页 2 CN 114265085 A 2其中， c为辐射响应时变模型系数。 计算得 出遥感器辐射响应的平均年衰变率可为： AnnualRatei＝(Fun2(0) ‑Fun2(DSL))*365/(DSL)。权　利　要　求　书 2/2 页 3 CN 114265085 A 3