[发明专利]一种空间异常极值无偏等值图绘制方法

摘要

本发明提供了一种空间异常极值无偏等值图绘制方法，对原始采样点集A提取极值点，得到极值点点集A’；所述极值点为局部邻域分析范围内，比周围所有点的空间属性值高或低的点；对原始采样点集A进行规则矩形网格化，得到规则矩形网格节点点集B；合并点集A’和点集B，得到点集C；删除点集C中导致病态三角形生成的畸形网格节点，得到点集D；利用点集D建立不规则三角网，绘制等值线图G。利用本发明绘制的等值图空间异常中心与极值点位置吻合，曲线光滑美观，在保证可视化效果的同时，提高了等值线精度。

主权项

一种空间异常极值无偏等值图绘制方法，利用离散的原始采样点集绘制，其特征在于具体包括以下步骤：(1)对原始采样点集A提取极值点，得到极值点点集A’；所述极值点为局部邻域分析范围内，比周围所有点的空间属性值高或低的点；(2)对原始采样点集A进行规则矩形网格化，得到规则矩形网格节点点集B及规则矩形网；(3)合并点集A’和点集B，得到点集C；(4)删除点集C中导致病态三角形生成的畸形网格节点，得到点集D；具体通过以下步骤删除畸形网格节点：(401)已知点集C中，来自规则矩形网格节点点集B的点的数量为N，来自极值点点集A’的点的数量为M；(402)对于来自极值点点集A’的一点a_m，且m∈[1,M]，确定点a_m在规则矩形网中所在的单元格R，确定单元格R中距离点a_m最近的点b_n，且点b_n来自规则矩形网格节点点集B，n∈[1,N]，并通过以下公式计算点b_n与点a_m的距离d：$d=\sqrt{{(x_{bn}-x_{am})}^2+{(y_{bn}-y_{am})}^2}\mathrm{...}\left(1\right)$其中，x_bn和y_bn分别为点b_n的x坐标和y坐标，x_am，y_am分别为点a_m的x坐标和y坐标；(403)通过以下公式计算网格R的对角线长度d’：$d^{\′}=\sqrt{X^2+Y^2}\mathrm{...}\left(2\right)$其中，X和Y分别为规则矩形网格节点点B中任意2个相邻点在x方向上的点距以及在y方向上的点距；(404)判断是否删除点b_n：给定阈值β，若d＜β×d′，则删除点b_n；否则保留点b_n；(405)重复步骤402到步骤405，直到对对极值点点集A’中的每个点判断一次是否删除距离其最近的网格顶点，最终剩下的所有点构成点集D；(5)利用点集D建立不规则三角网，绘制等值线图G，绘制结束。