[发明专利]一种海陆交错带气象等值线自动绘制方法

权利要求书

1.一种海陆交错带气象等值线自动绘制方法，其特征在于包括如下步骤：

1)从多层多变量的海洋气象原始数据文件中读取出待绘制等值线的一层数据，再获取与该数据格式一致的海陆分界数据，利用海陆分界数据对原始数据做叠加分析，分离出海洋区域数据；

2)设定判断噪点的阈值，遍历数据中的每个像元，若当前像元的数值与周围相连的八个像元的方差减去整体的方差大于给定的阈值，则认为该像元为噪点，对该像元进行校正，以此得到海洋区域有效数值范围；

3)根据海洋区域的有效数值范围得到全部区域的一个划分；

4)分别构建海洋区域和其他区域的间隔表，根据间隔表设置各个区域不同疏密的间隔；

5)合并海洋区域和其他区域的间隔；

6)利用MATLAB批量读取数据，并将步骤5)合并后的间隔应用于MATLAB等值线绘制指令，实现自动绘制。

2.根据权利要求1所述的一种海陆交错带气象等值线自动绘制方法，其特征在于所述的步骤1)包括：从一个多层多变量的海洋气象原始数据文件中读取出待绘制等值线的一层数据，该层数据是由n行m列的二维规则格网构成，再获取与该数据格式一致的n行m列的海陆分界数据，海陆分界数据包含与其相对应的二值数据，每个像元都引用相同的地理位置，海陆分界二值数据中0代表其他，1代表海洋，利用海陆分界数据对原始数据做叠加分析，分离出海洋区域数据。

3.根据权利要求1所述的一种海陆交错带气象等值线自动绘制方法，其特征在于所述的步骤2)包括：设定一个判断噪点的阈值，遍历数据中的每个像元，若当前像元的数值f(i，j)与周围相连的八个像元的方差σ²减去整体的方差大于给定的阈值，则认为该像元为噪点，该方法可用以下公式表示：

18Σl=18[g1-f(i,j)]2-σ2‾>ϵ---(1)]]>

式中：g_l为f(i，j)点周围邻点的数值；

当确认某像元为噪点时，用以下公式计算的结果替换原值，

f′(i,j)=(1-P)f(i,j)+P[18Σl=18gl]---(2)]]>

式中：f′(i，j)为像元的新亮度；P为该像元是噪点的概率，

经以上校正后的海洋区域数据可以得到海洋区域有效数值范围。