基于星地同步观测对比的卫星闪电探测评估方法.pdf

《基于星地同步观测对比的卫星闪电探测评估方法.pdf》由会员分享，可在线阅读，更多相关《基于星地同步观测对比的卫星闪电探测评估方法.pdf（15页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201911054631.0 (22)申请日 2019.10.31 (71)申请人 国家卫星气象中心 （国家空间天气 监测预警中心） 地址 100081 北京市海淀区中关村南大街 46号 (72)发明人 惠雯 (74)专利代理机构 北京纽乐康知识产权代理事 务所(普通合伙) 11210 代理人 田磊 (51)Int.Cl. G01R 29/08(2006.01) G06T 7/62(2017.01) (54)发明名称 一种基于星地同步观测对比的卫星闪电探 测评估方法 (57)摘要。

2、 本发明公开了一种基于星地同步观测对比 的卫星闪电探测评估方法， 其特征在于， 包括如 下步骤， S1:对一个典型区域在特定时间段内星 地同步观测的闪电进行对比， 以说明卫星与地基 闪电探测结果的相关性和差异性； S2:为解决卫 星与地基探测闪电信号不同表现形式的问题， 根 据地闪回击数据的时空粒度特征， 从卫星探测的 逐帧像元数据中提取一种新的卫星探测闪电数 据单元spot， 使其在时空粒度上与地基stroke具 有可比性； S3:为解决卫星与地基观测闪电种类 不完全一致的问题， 找到卫星与地基可以同时观 测闪电的特定位置； S4:对特定区域的卫星与地 基闪电探测结果进行比对， 定义二者的匹。

3、配度。 权利要求书2页 说明书5页 附图7页 CN 110837006 A 2020.02.25 CN 110837006 A 1.一种基于星地同步观测对比的卫星闪电探测评估方法， 其特征在于， 包括如下步骤： S1: 对一个典型区域在特定时间段内星地同步观测的闪电进行对比， 以说明卫星与地 基闪电探测结果的相关性和差异性； S2: 为解决卫星与地基探测闪电信号不同表现形式的问题， 根据地闪回击数据的时空 粒度特征， 从卫星探测的逐帧像元数据中提取一种新的卫星探测闪电数据单元spot， 使其 在时空粒度上与地基stroke具有可比性； S3: 为解决卫星与地基观测闪电种类不完全一致的问题， 找。

4、到卫星与地基可以同时观 测闪电的特定位置； S4: 对特定区域的卫星与地基闪电探测结果进行比对， 定义二者的匹配度。 2.根据权利要求1所述的一种基于星地同步观测对比的卫星闪电探测评估方法， 其特 征在于， 所述步骤S1中的所述特定时间段为2008-2013年。 3.根据权利要求1所述的一种基于星地同步观测对比的卫星闪电探测评估方法， 其特 征在于， 所述步骤S1中的卫星数据使用TRMM LIS数据， 观测区域为38 N38 S之间， 在每个 观测点的平均注视时间为90s， 星下点空间分辨率为36km； 地基数据使用我国雷电监测网 ADTD， 观测范围为我国陆地区域， 单站辐射半径为300km。

5、。 4.根据权利要求3所述的一种基于星地同步观测对比的卫星闪电探测评估方法， 其特 征在于， 研究区域选取我国西南地区。 5.根据权利要求1所述的一种基于星地同步观测对比的卫星闪电探测评估方法， 其特 征在于， 对于所述步骤S3， 找到卫星与地基可以同时观测闪电的特定位置可包括如下步骤： S31: 根据卫星轨道信息， 将卫星在每个观测时间点的ECR坐标转换为大地坐标系下的 经度Slon、 纬度Slat和高程hs， 利用下式来分别转换： 其中， X为ECR横坐标， Y为ECR纵坐标， Z为ECR竖坐标， a为地球长半径， b为地球短半径、 e2 为第一偏心率， 因卫星轨道高度不变， 将上述模型进。

6、行简化为： S32: 通过云检测和云高反演， 获取云边界分布和云高hc； S33: 根据云边界点的经度Clon和纬度Clat， 卫星星下点经度Slon和纬度Slat， 计算卫星相 对于云边界位置的球面距离， 利用下式来计算卫星相对于云边界位置的球面距离： 权利要求书 1/2 页 2 CN 110837006 A 2 其中， SOC为卫星星下点和云边界位置相对于地球球心的夹角； S34: 计算卫星探测闪电的有效半径ls， 利用下式来计算卫星探测闪电的有效半径ls： 其中， R为地球半径， 且将卫星高度hs， 云高hc， 地球半径R， 以及SOC代入上式， 得到， S35: 根据卫星探测闪电的有效。

7、半径ls， 以及卫星在各观测点的观测范围， 确定卫星与 地基可以同时观测闪电的区域， 作为特定区域。 6.根据权利要求1或5所述的一种基于星地同步观测对比的卫星闪电探测评估方法， 其 特征在于， 对于所述步骤S4， 将所述特定区域的卫星与地基闪电探测的匹配度定义为： 。 权利要求书 2/2 页 3 CN 110837006 A 3 一种基于星地同步观测对比的卫星闪电探测评估方法 0001 技术领域 0002 本发明涉及卫星闪电探测技术领域， 具体来说， 涉及一种基于星地同步观测对比 的卫星闪电探测评估方法。 背景技术 0003 闪电是大气中一种超长距离的强放电现象， 是对流性天气系统发展到一定。

8、阶段的 产物， 有效监测闪电信号， 对于研究闪电对强对流天气的指示作用， 以及闪电对大气环境和 气候的影响等具有指导意义。 0004 卫星闪电成像仪通过光谱滤波、 空间滤波、 时间滤波， 以及星上实时事件处理器 （Real Time Event Processor， RTEP） 的背景减光技术， 提取闪电信号。 但是， 卫星闪电探测 有可能造成自然界闪电的人为割裂， 使得一个闪电被分割到两帧数据， 或是被分割到若干 相邻探测单元， 因此， 需要通过地面系统， 采用不同层级数据对闪电信息进行重建。 0005 卫星闪电探测的基本数据单元是闪电事件 （event） ， 卫星闪电成像仪在轨运行中， 因。

9、仪器内部原因或轨道外部环境干扰， 往往造成原始闪电数据里包含大量由非闪电信号形 成的发光像元， 称为虚假闪电event， 根据以往观测经验， 虚假闪电event数量最多会占到全 部闪电event总数的一半以上， 因此在闪电数据的地面处理中需对其进行针对性滤除。 0006 卫星闪电探测过程中， 一方面， RTEP星上处理的背景估计和背景滤除过程可能引 入误差， 另一方面， 对虚假闪电信号的滤除也可能产生误判， 诸如此类因素， 不可避免地会 影响卫星闪电探测结果， 因此， 需要有效方法对卫星闪电探测的正确性进行评估， 可以利用 地基闪电定位网观测的闪电信息验证卫星探测的闪电信号， 通过二者的匹配度。

10、来检验卫星 闪电探测结果， 但卫星与地基观测的目标信号特征存在差异， 无法直接将两种观测数据进 行对比。 发明内容 0007 针对相关技术中的上述技术问题， 本发明提出一种基于星地同步观测对比的卫星 闪电探测评估方法， 能够克服现有技术的上述不足。 0008 为实现上述技术目的， 本发明的技术方案是这样实现的： 一种基于星地同步观测对比的卫星闪电探测评估方法， 其特征在于， 包括如下步骤： S1: 对一个典型区域在特定时间段内星地同步观测的闪电进行对比， 以说明卫星与地 基闪电探测结果的相关性和差异性； S2: 为解决卫星与地基探测闪电信号不同表现形式的问题， 根据地闪回击数据的时空 粒度特征。

11、， 从卫星探测的逐帧像元数据中提取一种新的卫星探测闪电数据单元spot， 使其 在时空粒度上与地基stroke具有可比性； S3: 为解决卫星与地基观测闪电种类不完全一致的问题， 找到卫星与地基可以同时观 说明书 1/5 页 4 CN 110837006 A 4 测闪电的特定位置； S4: 对特定区域的卫星与地基闪电探测结果进行比对， 定义二者的匹配度。 0009 进一步的， 所述步骤S1中的所述特定时间段为2008-2013年。 0010 进一步的， 所述步骤S1中的卫星数据使用TRMM LIS数据， 观测区域为38 N38 S 之间， 在每个观测点的平均注视时间为90s， 星下点空间分辨率。

12、为36km； 地基数据使用我 国雷电监测网ADTD， 观测范围为我国陆地区域， 单站辐射半径为300km。 0011 其中， 研究区域选取我国西南地区。 0012 进一步的， 对于所述步骤S3， 找到卫星与地基可以同时观测闪电的特定位置可包 括如下步骤： S31: 根据卫星轨道信息， 将卫星在每个观测时间点的ECR坐标转换为大地坐标系下的 经度Slon、 纬度Slat和高程hs， 利用下式来分别转换： 其中， X为ECR横坐标， Y为ECR纵坐标， Z为ECR竖坐标， a为地球长半径， b为地球短半径、 e2为第一偏心率， 因卫星轨道高度不变， 将上述模型进行简化为： S32: 通过云检测和云。

13、高反演， 获取云边界分布和云高hc； S33: 根据云边界点的经度Clon和纬度Clat， 卫星星下点经度Slon和纬度Slat， 计算卫星相 对于云边界位置的球面距离， 利用下式来计算卫星相对于云边界位置的球面距离： 其中， SOC为卫星星下点和云边界位置相对于地球球心的夹角； S34: 计算卫星探测闪电的有效半径ls， 利用下式来计算卫星探测闪电的有效半径ls： 其中， R为地球半径， 且将卫星高度hs， 云高hc， 地球半径R， 以及SOC代入上式， 得到， S35: 根据卫星探测闪电的有效半径ls， 以及卫星在各观测点的观测范围， 确定卫星与 说明书 2/5 页 5 CN 110837。

14、006 A 5 地基可以同时观测闪电的区域， 作为特定区域。 0013 进一步的， 对于所述步骤S4， 将所述特定区域的卫星与地基闪电探测的匹配度定 义为： 。 0014 本发明的有益效果： 本发明利用地基闪电定位网观测的闪电信息验证卫星探测的 闪电信号， 通过二者的匹配度来检验卫星闪电探测结果， 对卫星闪电探测结果进行定量评 估， 为更有效地在卫星平台上监测闪电信号， 实现对强对流天气系统的追踪和预警创造有 利条件。 附图说明 0015 为了更清楚地说明本发明实施例或现有技术中的技术方案， 下面将对实施例中所 需要使用的附图作简单地介绍， 显而易见地， 下面描述中的附图仅仅是本发明的一些实施。

15、 例， 对于本领域普通技术人员来讲， 在不付出创造性劳动的前提下， 还可以根据这些附图获 得其他的附图。 0016 图1是根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电探测评估 方法的步骤流程图； 图2是根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电探测评估方法 的我国西南地区卫星flash和地基flash的相对变化量的月变化折线图； 图3是根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电探测评估方法 的我国西南地区卫星flash和地基flash的相对变化量的日变化折线图； 图4是根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电探测评估方法 的我国西南地区卫。

16、星flash和地基flash的绝对变化量的月变化折线图； 图5是根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电探测评估方法 的我国西南地区卫星flash和地基flash的绝对变化量的日变化折线图； 图6是根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电探测评估方法 的我国西南地区2008-2013年星地匹配度月变化折线图； 图7是根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电探测评估方法 的我国西南地区2008-2013年星地匹配度日变化折线图。 具体实施方式 0017 下面将结合本发明实施例中的附图， 对本发明实施例中的技术方案进行清楚、 完 整地描述， 显然， 。

17、所描述的实施例仅仅是本发明一部分实施例， 而不是全部的实施例。 基于 本发明中的实施例， 本领域普通技术人员所获得的所有其他实施例， 都属于本发明保护的 范围。 0018 如图1-7所示， 根据本发明实施例所述的一种基于星地同步观测对比的卫星闪电 探测评估方法， 其特征在于， 包括如下步骤： S1: 对一个典型区域在特定时间段内星地同步观测的闪电进行对比， 以说明卫星与地 基闪电探测结果的相关性和差异性； 说明书 3/5 页 6 CN 110837006 A 6 S2: 为解决卫星与地基探测闪电信号不同表现形式的问题， 根据地闪回击数据的时空 粒度特征， 从卫星探测的逐帧像元数据中提取一种新的。

18、卫星探测闪电数据单元spot， 使其 在时空粒度上与地基stroke具有可比性； S3: 为解决卫星与地基观测闪电种类不完全一致的问题， 找到卫星与地基可以同时观 测闪电的特定位置； S4: 对特定区域的卫星与地基闪电探测结果进行比对， 定义二者的匹配度。 0019 在一具体实施例中， 步骤S1中的所述特定时间段为2008-2013年。 0020 在一具体实施例中， 步骤S1中的卫星数据使用TRMM LIS数据， 观测区域为38 N 38 S之间， 在每个观测点的平均注视时间为90s， 星下点空间分辨率为36km； 地基数据使 用我国雷电监测网ADTD， 观测范围为我国陆地区域， 单站辐射半径。

19、为300km。 0021 优选的， 研究区域选取我国西南地区。 0022 在一具体实施例中， 对于步骤S2， 可从表1中看出我国西南地区2008-2013年LIS spot与ADTD stroke的时空粒度具有相近的数值。 0023 表1 LIS spot与ADTD stroke的时空粒度数值表 时间粒度 （s）空间粒度 （km） LIS spot0.1112.143 ADTD stoke0.1462.271 在一具体实施例中， 对于步骤S3， 为解决卫星与地基观测闪电种类不完全一致的问题， 需考虑云的影响， 找到卫星与地基可以同时观测闪电的特定位置， 包括如下步骤： S31: 根据卫星轨道信。

20、息， 将卫星在每个观测时间点的ECR坐标转换为大地坐标系下的 经度Slon、 纬度Slat和高程hs， 利用下式来分别转换： 其中， X为ECR横坐标， Y为ECR纵坐标， Z为ECR竖坐标， a为地球长半径， b为地球短半径、 e2为第一偏心率， 因卫星轨道高度不变， 将上述模型进行简化为： S32: 通过云检测和云高反演， 获取云边界分布和云高hc； S33: 根据云边界点的经度Clon和纬度Clat， 卫星星下点经度Slon和纬度Slat， 计算卫星相 对于云边界位置的球面距离， 利用下式来计算卫星相对于云边界位置的球面距离： 其中， SOC为卫星星下点和云边界位置相对于地球球心的夹角；。

21、 说明书 4/5 页 7 CN 110837006 A 7 S34: 计算卫星探测闪电的有效半径ls， 利用下式来计算卫星探测闪电的有效半径ls： 其中， R为地球半径， 且将卫星高度hs， 云高hc， 地球半径R， 以及SOC代入上式， 得到， S35: 根据卫星探测闪电的有效半径ls， 以及卫星在各观测点的观测范围， 确定卫星与 地基可以同时观测闪电的区域， 作为特定区域。 0024 在一具体实施例中， 步骤S31中的地球长半径a、 短半径b、 第一偏心率e2为国际大 地测量与地球物理联合会1975年推荐的参考椭球参数。 0025 在一具体实施例中， 对于步骤S4， 将所述特定区域的卫星与。

22、地基闪电探测的匹配 度定义为： 。 0026 为了方便理解本发明的上述技术方案， 以下通过具体使用方式上对本发明的上述 技术方案进行详细说明。 0027 在具体使用时， 根据本发明所述的一种基于星地同步观测对比的卫星闪电探测评 估方法， 首先， 对一个典型区域在特定时间段内星地同步观测的闪电进行对比， 以说明卫星 与地基闪电探测结果的相关性和差异性， 在此步骤， 研究区域选取我国西南地区； 然后， 根 据地闪回击数据的时空粒度特征， 从卫星探测的逐帧像元数据中提取一种新的卫星探测闪 电数据单元spot， 使其在时空粒度上与地基stroke具有可比性， 解决了卫星与地基探测闪 电信号不同表现形式。

23、的问题； 再找到卫星与地基可以同时观测闪电的特定位置； 最后， 对特 定区域的卫星与地基闪电探测结果进行比对， 定义二者的匹配度， 若匹配度接近于1， 说明 利用本算法可明显提高卫星与地基闪电探测数据的可比性。 0028 综上所述， 利用本发明中的算法， 定义出我国西南地区2008-2013年星地匹配度 后， 可得到， 月变化的星地匹配度为0.82， 日变化的星地匹配度为0.81， 均接近于1， 说明利 用本算法可明显提高卫星与地基闪电探测数据的可比性。 因此， 在评估卫星闪电探测结果 时， 将星上实时事件处理器输出的闪电event转化为闪电spot， 再从中提取位于特定位置的 闪电， 将这些。

24、闪电与地基观测结果进行对比， 可以获得更准确的检验结果。 0029 以上所述仅为本发明的较佳实施例而已， 并不用以限制本发明， 凡在本发明的精 神和原则之内， 所作的任何修改、 等同替换、 改进等， 均应包含在本发明的保护范围之内。 说明书 5/5 页 8 CN 110837006 A 8 图1 说明书附图 1/7 页 9 CN 110837006 A 9 图2 说明书附图 2/7 页 10 CN 110837006 A 10 图3 说明书附图 3/7 页 11 CN 110837006 A 11 图4 说明书附图 4/7 页 12 CN 110837006 A 12 图5 说明书附图 5/7 页 13 CN 110837006 A 13 图6 说明书附图 6/7 页 14 CN 110837006 A 14 图7 说明书附图 7/7 页 15 CN 110837006 A 15 。