买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：山东省海洋预报减灾中心

摘要：本发明涉及数据处理领域，具体为一种海洋数值模拟产品后处理方法及装置。其包括如下步骤：Step1、读取shapefile格式的数据记录；Step2、读取数值预报产品数据的网格信息，通过网格信息判断网格类型是规则网格还是非规则网格；Step3、初始化预报要素数据有效性数组；Step4、分别针对规则网格和非规则网格计算有效性数组2；Step5、计算有效性数组2和有效性数组1的差异部分，针对差异部分开展数据处理；Step6、经插值处理后，生成与基础地理信息相符的海洋预报产品数据，按照格式要求输出，并开展后续制图输出。本发明能兼容海洋数值预报模式输出网格化产品中的规则网格和非规则网格，提升了普适性。

主权项：1.一种海洋数值模拟产品后处理方法，其特征在于，包括如下步骤Step1-Step6：Step1、读取shapefile格式的数据记录，如果记录类型为multi-polygon格式，将其处理为单个polygon集合，并记录每个polygon的BBox参数；Step2、读取数值预报产品数据的网格信息，通过网格信息判断网格类型是规则网格还是非规则网格，判断依据为分别遍历经度、纬度变量，如果两者为等间距间隔，且经度、纬度间隔值相同，则为规则网格；否则为非规则网格；Step3、初始化预报要素数据有效性数组：有效性数组1初始化：在数值预报模式下，在输出预报结果时将网格区域的非海洋部分设置为固定值NaN，先读取任意预报时刻的预报场，复制其数组，并将NaN值的部分设置为0；有效性数组2初始化：初始化数据组，其维度为n、lat、lon，其中n为Step1中的polygon数量，lat为Step2中的纬度数据、lon为Step2中的经度数据，初始化过程中，全部数据设置为0；Step4、分别针对规则网格和非规则网格计算有效性数组2，针对规则网格计算有效性数组2时，以8倍空间分辨率预生成正方形网格集合，然后逐正方形网格遍历，遍历中的操作为步骤A1-A6：A1、有效性数组2的n为循环控制；A2、在循环内，判断正方形与地理信息数据中的每个多边形的关系，多边形包含正方形，将有效性数组中对应正方形的网格点设置为1；A3、如果多边形与正方形无交集，则跳过该次循环；A4、如果多边形与正方形为相交关系，则将正方形分辨率调整为n2，其中，n不等于1；A5、重复步骤A2-A4，终止条件为：如果n＝1时，判断该多边形是否包括该点，如果是，则设置点所在位置的有效性数组2的值为1；A6、n循环完毕后，将得到的n个lat、lon数组合并，得到有效性数组2；针对非规则网格计算有效性数组2时，进行如下步骤B1-B6：B1、随机抽样3次，每次抽取网格点中5％的网格点，计算网格点之间的平均距离[D1、D2、D3]；B2、确定平均距离Davg所对应的S2空间索引的层级，其中，Davg为3次距离的均值；B3、遍历一遍非规则网格，计算每个网格点的S2空间索引值，并记录每个空间索引值相同的多边形对应的网格点记录集合，遍历完毕后得到空间索引多边形集合；B4、有效性数组2的n为循环控制；B5、在循环内，遍历空间索引多边形，操作如下步骤B51-B53：B51、判断空间索引多边形与地理信息数据中每个多边形的关系，多边形包含空间索引多边形，将有效性数组中对应空间索引多边形的网格点设置为1；B52、如果多边形与空间索引多边形无交集，则跳过该次循环；B53、如果多边形与空间索引多边形为相交关系，则依次判断多边形空间索引多边形中的每个点的关系，如果存在包含关系则设置该网格位置的值为1；B6、n循环完毕后，将得到的n个lat、lon数组合并，得到有效性数组2；Step5、计算有效性数组2和有效性数组1的差异部分，针对差异部分开展数据处理；Step6、经插值处理后，生成与基础地理信息相符的海洋预报产品数据，按照格式要求输出，并开展后续制图输出。

全文数据：

权利要求：

百度查询： 山东省海洋预报减灾中心 一种海洋数值模拟产品后处理方法及装置