一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法技术

本发明专利技术公开了一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法，包括以下步骤：观测目标区域选择；筛选晴空无耀斑影响的目标区域；均匀目标筛选，输入信号和噪声信号分离；传感器输入信号和噪声的确定；计算传感器信噪比SNR，本发明专利技术适用于卫星观测可见光谱段信噪比在轨估算技术领域，对基于实际卫星遥感器对地观测数据计算信噪比进行了统一的量化方法确定，基于卫星遥感器无云低亮度场景下的在轨观测数据实时估算信噪比；可以直接利用卫星遥感器对地观测数据有效分离观测信号中的噪声，得到对应观测辐射下的仪器噪声和信噪比，用于卫星遥感器在轨信噪比实时估算，获得仪器SNR在轨变化监测结果，用以监视仪器在轨性能的变化。

【技术实现步骤摘要】
一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法
本专利技术属于卫星观测信噪比在轨估算
，具体是一种用于低亮度区域气象卫星仪器可见光谱段的信噪比在轨量化方法。
技术介绍
卫星传感器设计需具有足够的信噪比SNR，以解决对观测区域的细微变化测量。一般使用以下两类方法估算卫星仪器噪声水平：一类是指定仪器等效噪声辐射而不是SNR，此方法带来的误差较大，难以满足高精度的观测要求；二类是在特定输入辐射条件下，通过实验室测量计算出SNR。信噪比计算需确定卫星传感器观测时的能量信号与仪器的噪声。然而,卫星传感器的噪声与观测时的辐射信号有关；现有的在轨计算SNR的技术方案通常采用星上定标源观测数据进行信噪比的测算，其使用了积分球光源、太阳和漫反/透射板光源作为输入信号，不同仪器所设定的输入辐射信号不同，导致在轨测算信噪比所参照的输入辐射信号不同，所得到的SNR不具备可比性。
技术实现思路
本专利技术的目的在于克服现有技术的缺陷，提供一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法。为实现上述目的，本专利技术采用了如下技术方案：一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法，包括以下步骤：观测目标区域选择；筛选晴空无耀斑影响的目标区域；均匀目标筛选，输入信号和噪声信号分离；传感器输入信号和噪声的确定；计算传感器信噪比SNR。优选的，所述观测目标区域选择，选择原则包括：区域观测能量相对集中稳定、潜在的均匀目标多、晴空率高及目标易于获取。优选的，所述筛选晴空无耀斑影响的目标区域，包括：设定太阳天顶角筛选阈值、耀斑角筛选阈值及晴空筛选阈值；基于卫星遥感器定标参数计算得到大气顶反射率，若目标区域太阳天顶角＞太阳天顶角筛选阈值1且太阳天顶角＜太阳天顶角筛选阈值2、耀斑角＞耀斑角筛选阈值且大气顶反射率＜晴空筛选阈值，则保留，反之剔除。优选的，所述大气顶反射率计算公式为：Ref＝(EV-SV)*Slope*d2/cos(SZA)，其中，EV与SV为对地观测和冷空观测数码值，Slope为反射率定标斜率，d为日地距离，SZA为太阳天顶角。优选的，所述均匀目标筛选，输入信号与噪声信号的分离，包括：设定空间窗口及筛选阈值，基于窗口内的阈值方案，若窗口的高于此阈值，剔除该窗口数据，反之则判定该空间窗口为均一目标。优选的，所述均匀目标筛选，输入信号与噪声信号的分离，还包括：对判定为均一目标的空间窗口，计算均值MEAN和标准差STD，进行数据累积。优选的，所述传感器输入信号和噪声的确定，包括：基于累积的标准差STD数据，计算概率分布直方图，进行高斯拟合确定该噪声数据是否符合正态分布，若满足检验的显著性水平，则判定此噪声符合正态分布，以符合正态分布标准差STD数据的均值作为传感器噪声Rnoise，以对应的均值MEAN数据的均值作为传感器输入信号Rsignal。优选的，所述计算传感器信噪比SNR，计算公式为：综上所述，由于采用了上述技术方案，本专利技术的有益效果是：本专利技术对基于实际卫星遥感器对地观测数据计算信噪比进行了统一的量化方法确定，基于卫星遥感器无云低亮度场景下的在轨观测辐射数据实时估算信噪比；本专利技术可以直接利用卫星遥感器对地观测数据有效分离观测信号中的噪声，得到低亮度观测辐射下的仪器噪声和信噪比，可以将此方法用于卫星遥感器在轨信噪比实时估算，获得仪器SNR在轨变化监测结果，用以监视仪器在轨性能的变化。附图说明图1是本专利技术一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法的流程图；图2是本专利技术优选实施例中均匀目标阈值选择为1.0075时的FY-3B/MERSI波段8(412nm)噪声分布直方图；图3是本专利技术优选实施例中FY-3B/MERSI波段8(412nm)、9(443nm)和10(490nm)的SNR和对应信号辐亮度。具体实施方式以下结合附图1-3，进一步说明本专利技术一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法的具体实施方式。本专利技术一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法不限于以下实施例的描述。实施例1：本实施例给出一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法的具体结构，如图1所示，包括以下步骤：观测目标区域选择；筛选晴空无耀斑影响的目标区域；均匀目标筛选，输入信号和噪声信号分离；传感器输入信号和噪声的确定；计算传感器信噪比SNR。具体地，观测目标区域选择，选择原则包括：区域观测能量相对集中稳定、潜在的均匀目标多、晴空率高及目标易于获取。具体地，筛选晴空无耀斑影响的目标区域，包括：设定太阳天顶角筛选阈值、耀斑角筛选阈值及晴空筛选阈值；基于卫星遥感器定标参数计算得到大气顶反射率，若目标区域太阳天顶角＞太阳天顶角筛选阈值1且太阳天顶角＜太阳天顶角筛选阈值2、耀斑角＞耀斑角筛选阈值且大气顶反射率＜晴空筛选阈值，则保留，反之剔除。具体地，大气顶反射率计算公式为：Ref＝(EV-SV)*Slope*d2/cos(SZA)，其中，EV与SV为对地观测和冷空观测数码值，Slope为反射率定标斜率，d为日地距离，SZA为太阳天顶角。具体地，均匀目标筛选，输入信号与噪声信号的分离，包括：设定空间窗口及筛选阈值，基于窗口内的阈值方案，若窗口的高于此阈值，剔除该窗口数据，反之则判定该空间窗口为均一目标。进一步的，均匀目标筛选，输入信号与噪声信号的分离，还包括：对判定为均一目标的空间窗口，计算均值MEAN和标准差STD，进行数据累积。进一步的，传感器输入信号和噪声的确定，包括：基于累积的标准差STD数据，计算概率分布直方图，进行高斯拟合确定该噪声数据是否符合正态分布，若满足检验的显著性水平，则判定此噪声符合正态分布，以符合正态分布标准差STD数据的均值作为传感器噪声Rnoise，以对应的均值MEAN数据的均值作为传感器输入信号Rsignal；进一步的，计算传感器信噪比SNR，计算公式为：实施例2：本实施例给出了风云三号气象卫星B星(FY-3B)的中分辨率光谱成像仪(MERSI)的具体结果，如图1-3所示，包括以下步骤：1.基于设定的贫瘠大洋区域，提取遥感器观测数据：设定观测目标区域范围为-32.3°E～-11°E,-25°N～0°；由FY-3B/MERSIL1B文件读取设定经纬度范围内的数据，提取对地观测数码值EV、冷空观测数码值SV、太阳天顶角SZA，利用定标参数计算大气顶反射率：Ref＝(EV-SV)*Slope*d2/cos(SZA)其中，EV为对地观测，SV为冷空观测数码值，Slope为反射率定标斜率，d为日地距离，SZA为太阳天顶角。2.晴空无耀斑影响的数据筛选：设定太阳天顶角SZA筛选阈值2为70°，本文档来自技高网...

【技术保护点】
1.一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法，其特征在于，包括以下步骤：/n观测目标区域选择；/n筛选晴空无耀斑影响的目标区域；/n均匀目标筛选，输入信号和噪声信号分离；/n传感器输入信号和噪声的确定；/n计算传感器信噪比SNR。/n

【技术特征摘要】
1.一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法，其特征在于，包括以下步骤：
观测目标区域选择；
筛选晴空无耀斑影响的目标区域；
均匀目标筛选，输入信号和噪声信号分离；
传感器输入信号和噪声的确定；
计算传感器信噪比SNR。


2.如权利要求1所述的一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法，其特征在于，所述观测目标区域选择，选择原则包括：区域观测能量相对集中稳定、潜在的均匀目标多、晴空率高及目标易于获取。


3.如权利要求2所述的一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法，其特征在于，所述筛选晴空无耀斑影响的目标区域，包括：
设定太阳天顶角筛选阈值、耀斑角筛选阈值及晴空筛选阈值；
基于卫星遥感器定标参数计算得到大气顶反射率，若目标区域太阳天顶角＞太阳天顶角筛选阈值1且太阳天顶角＜太阳天顶角筛选阈值2、耀斑角＞耀斑角筛选阈值且大气顶反射率＜晴空筛选阈值，则保留，反之剔除。


4.如权利要求3所述的一种用于低亮度区域气象卫星仪器的信噪比在轨量化方法，其特征在于，所述大气顶反射率计算公式为：
Ref＝(EV-SV)*Slope*d2/cos(SZA)，
其中，EV与SV为对地观测和冷空观测数码值...

【专利技术属性】
技术研发人员：孙凌，张鹏，干佳俪，薛文文，袁鸣鸽，
申请(专利权)人：国家卫星气象中心国家空间天气监测预警中心，
类型：发明
国别省市：北京;11