基于微波散射计数据自动提取台风高风速云系区域方法

摘要

本发明公开了一种基于微波散射计数据自动提取台风高风速云系区域方法。它包括以下步骤：S1：将从微波散射计获取的台风区域的风速数据、风向数据进行网格化处理，形成网格数据点；S2：计算出每个网格数据点的风矢量；S3：用设定大小的窗口遍历台风区域获得若干个窗口区域，计算每个窗口区域的风速均值；S4：风速均值最大的窗口区域为台风高风速云系区域。本发明能够根据微波散射计数据自动提取出台风高风速云系区域，花费时间短，能够满足业务化发展的要求。

权利要求

1.一种基于微波散射计数据自动提取台风高风速云系区域方法，其特征在于，包括以下步骤：S1：将从微波散射计获取的台风区域的风速数据、风向数据按照经纬度进行网格化处理，形成网格数据点；

S2：计算出每个网格数据点的风矢量；

S3：用设定大小的窗口遍历台风区域获得若干个窗口区域，计算每个窗口区域的风速均值；包括以下步骤：D1：计算每个窗口区域内所有有效网格数据点的风矢量之和的模；

其中，如果网格数据点有风矢量数据，则该网格数据点为有效网格数据点，如果网格数据点没有风矢量数据，则该网格数据点为无效网格数据点；

D2：将每个窗口区域的风矢量之和的模除以其内所有有效网格数据点个数得到每个窗口区域的风速均值；

S4：风速均值最大的窗口区域为台风高风速云系区域。

2.根据权利要求1所述的基于微波散射计数据自动提取台风高风速云系区域方法，其特征在于，所述台风区域的风速数据的网格化处理包括以下步骤：N1：获取台风区域的微波散射计数据，经过预处理，得到风速数据；

N2：将风速数据通过反距离加权插值，得到网格化数据，每个网格的经度长度为a°、纬度长度为a°。

3.根据权利要求1所述的基于微波散射计数据自动提取台风高风速云系区域方法，其特征在于，所述台风区域的风向数据的网格化处理包括以下步骤：M1：获取台风区域的微波散射计数据，经过预处理，得到风向数据；

M2：将风向数据先通过反距离加权插值，再通过单方向线性插值，得到网格化数据，每个网格的经度长度为a°、纬度长度为a°。

说明书

基于微波散射计数据自动提取台风高风速云系区域方法

技术领域

[0001] 本发明涉及遥感技术领域，尤其涉及一种基于微波散射计数据自动提取台风高风速云系区域方法。

背景技术

[0002] 现有台风高风速云系区域提取主要依赖目视解译人工判定，自动提取的方法较为缺乏。主要技术方法有：

[0003] 1、反距离加权法；

[0004] 2、逐步订正法；

[0005] 3、克里金插值法；

[0006] 4、优插值法；

[0007] 5、三维时空插值方法。

[0008] 现有的这些风场反演方法具有以下缺陷：(1)基于人工目视解译提取高风速云系区域的方法受限于判定者的经验，不同的判定人的结果不同；

[0009] (2)单一对风向进行反距离加权插值，容易对某些区域的风向造成误判，且容易造成整轨数据北向倒置；

[0010] (3)在对一轨微波散射计进行气旋中心定位时，若不进行高风速云系区域提取，可能会定位出多个气旋中心。

[0011] 目前自动提取台风高风速云系区域的方法很少，利用海洋二号微波散射计数据进行自动提取的算法极少。海洋二号微波散射计数据提取台风高风速云系区域技术还处于起步阶段，其方法较为单一，多以单一插值方法和求算特定范围的风速均值为主，而微波散射计数据存储方式使得单一的插值方法对风向结果造成误判，结果无法满足提取的需要。

发明内容

[0012] 本发明为了解决现有人工目视解译判定微波散射计高风速云系区域技术所存在的花费时间长、标准不同一、无法做到业务化的技术问题，提供了一种基于微波散射计数据自动提取台风高风速云系区域方法，其能够根据微波散射计数据自动提取出台风高风速云系区域，花费时间短，能够满足业务化发展的要求。

[0013] 为了解决上述问题，本发明采用以下技术方案予以实现：

[0014] 本发明的基于微波散射计数据自动提取台风高风速云系区域方法，包括以下步骤：

[0015] S1：将从微波散射计获取的台风区域的风速数据、风向数据进行网格化处理，形成网格数据点；

[0016] S2：计算出每个网格数据点的风矢量；

[0017] S3：用设定大小的窗口遍历台风区域获得若干个窗口区域，计算每个窗口区域的风速均值；

[0018] S4：风速均值最大的窗口区域为台风高风速云系区域。

[0019] 本方法所处理的对象是来自于微波散射计数据，由一个个的数据点构成，微波散射计数据的数据点不是标准的网格，因此先要对台风区域的风速数据、风向数据进行网格化处理。

[0020] 作为优选，所述台风区域的风速数据的网格化处理包括以下步骤：

[0021] N1：获取台风区域的微波散射计数据，经过预处理，得到风速数据；

[0022] N2：将风速数据通过反距离加权插值，得到的网格化数据，每个网格的经度长度为a°、纬度长度为a°。

[0023] 作为优选，所述台风区域的风向数据的网格化处理包括以下步骤：

[0024] M1：获取台风区域的微波散射计数据，经过预处理，得到风向数据；

[0025] M2：将风向数据先通过反距离加权插值，再通过单方向线性插值，得到的网格化数据，每个网格的经度长度为a°、纬度长度为a°。

[0026] 对风向数据单方向线性插值包括以下步骤：比较待插入点东西相邻点的风速大小，以风速大的相邻点的东西走向的风向作为待插入点东西走向的风向；比较待插入点南北相邻点的风速大小，以风速大的相邻点的南北走向的风向作为待插入点南北走向的风向；根据待插入点东西走向的风向和南北走向的风向计算出待插入点的风向。

[0027] 作为优选，所述计算每个窗口区域的风速均值的方法包括以下步骤：

[0028] D1：计算每个窗口区域内所有有效网格数据点的风矢量之和的模；

[0029] D2：将每个窗口区域的风矢量之和的模除以其内所有有效网格数据点个数得到每个窗口区域的风速均值。

[0030] 作为优选，如果网格数据点有风矢量数据，则该网格数据点为有效网格数据点，如果网格数据点没有风矢量数据，则该网格数据点为无效网格数据点。

[0031] 本发明的有益效果是：可根据微波散射计数据自动提取出台风高风速云系区域，省时省力，更适合在业务中稳定运行。

附图说明

[0032] 图1是本发明的流程图。

具体实施方式

[0033] 下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。

[0034] 实施例：本实施例的基于微波散射计数据自动提取台风高风速云系区域方法，如图1所示，包括以下步骤：

[0035] S1：将从微波散射计获取的台风区域的风速数据、风向数据按照经纬度进行网格化处理，形成网格数据点；

[0036] S2：计算出每个网格数据点的风矢量；

[0037] S3：用设定大小的窗口(经度85°×纬度85°)遍历台风区域获得若干个窗口区域，计算每个窗口区域的风速均值；

[0038] S4：风速均值最大的窗口区域为台风高风速云系区域。

[0039] 台风区域的风速数据的网格化处理包括以下步骤：

[0040] N1：获取台风区域的微波散射计数据，经过预处理，得到风速数据；

[0041] N2：将风速数据通过反距离加权插值，得到网格化数据，每个网格为经度0.25°×纬度0.25°。

[0042] 台风区域的风向数据的网格化处理包括以下步骤：

[0043] M1：获取台风区域的微波散射计数据，经过预处理，得到风向数据；

[0044] M2：将风向数据先通过反距离加权插值，再通过单方向线性插值(只考虑同一方向水平和垂直)，得到网格化数据，每个网格为经度0.25°×纬度0.25°。

[0045] 计算每个窗口区域的风速均值的方法包括以下步骤：

[0046] D1：计算每个窗口区域内所有有效网格数据点的风矢量之和的模；

[0047] D2：将每个窗口区域的风矢量之和的模除以其内所有有效网格数据点个数得到每个窗口区域的风速均值。

[0048] 如果网格数据点有风矢量数据，则该网格数据点为有效网格数据点，如果网格数据点没有风矢量数据，则该网格数据点为无效网格数据点。

[0049] 本方法所处理的对象是来自于海洋二号卫星微波散射计数据，由一个个的数据点构成，海洋二号卫星微波散射计数据的数据点不是标准的网格，因此先要对台风区域的风速数据、风向数据进行网格化处理，包括如下步骤：将台风区域网格化，从微波散射计获取台风区域的风速数据、风向数据，并将这些风速数据、风向数据对应到网格数据点上。既有风速数据又有风向数据的网格数据点就是有效网格数据点，否则就是无效网格数据点。

[0050] 每个网格的经度长度为0.25°、纬度长度为0.25°。网格化后会有缺省点因为海洋二号卫星微波散射计测的是海面风速，会受海表面地物(岛屿)之类的影响，因此部分网格数据点是没有数据的无效点。对风向和风速进行插值是为了使图像更平滑，插值得到的有效数据点的数量增多，人为增大的图像的分辨率。规定一定大小的窗口(经度85°×纬度85°)遍历所有数据，将将每个窗口区域的风矢量之和的模除以其内所有有效网格数据点个数得到每个窗口区域的风速均值，风速均值最大的窗口即为高风速云系区域所在窗口。

[0051] 初始海洋二号卫星微波散射计数据为离散的点，风向的反距离加权插值法与单方向线性插值法的结合能避免风向异常情况出现。每个窗口得到的风矢量之和的模除以其内所有有效网格数据点个数，可排除高风速云系区域出现在图像边缘使致出现无数据点导致判定不准确情况

[0052] 对风向数据单方向线性插值包括以下步骤：比较待插入点东西相邻点的风速大小，以风速大的相邻点的东西走向的风向作为待插入点东西走向的风向；比较待插入点南北相邻点的风速大小，以风速大的相邻点的南北走向的风向作为待插入点南北走向的风向；根据待插入点东西走向的风向和南北走向的风向计算出待插入点的风向。

[0053] 本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

[0054] 尽管本文较多地使用了矢量参数、插值等术语，但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质；把它们解释成任何一种附加的限制都是与本发明精神相违背的。