地理坐标下网格化数据等值线的处理及编程.pdf

维普资讯 第 2期 (总 第 109期 ) 山 西 地 震 No．2 2O02年 4月 EARTHQ AKE RESEARCH IN SHANXI Apr． · 地震预报 · 文章编号 2002)020020G4 地理坐标下网格化数据等值线的处理及编程 李圣强，姜立新．李闽峰，杨满栋，米宏亮 (中国地震局分析预报 中心 ，北京 100036) 摘要 ：分析 了均匀网格化数据等值线处l理 的具体步骤 ，介绍 r非均匀 网格分布数据 的均 匀同格化 处理 中的 3种方法 按距离加权的最小二乘法 (N P法)，按方位取点加权法和加权 最小二乘法拟 合法 (M s法)，提供 rN P法和 M S珐的 C¨源程序。 关键词 ：地震数据处理 ；网格化数据 ；等值线处理 ；N—P法 ；M～s法 中田分类号：P31j75 文献标识码 ：A b)求交点 判定等值线 与网格棱边是否相交 ，线 引言 性求交 在处理地震前兆观测数据和研究其空间前兆场变 c)等值线在网格内的游动 如图 1所示； 化特征时，往往需要绘制一些观测数据 的等值线 图或 d)终点确定 等值线节点光滑 曲线处理、3次样 色块充填 图。实现这些功能的工具很多，如 SURFER 条曲线拟告和抛物线拟合。 软件、G0IDEN软件等，但是，倘若 自己编写一些较为 专业 的程序模块而且又需要嵌人能绘制空间数据 的功 能时，借助现有的商业软件来实现就显得不够经济和 方便，而且有时还会引起侵权 或非法使用 等纠纷 ．其 实 实现这些功能的接 口模块并不难。本文分析 lr网格 化数据等值线处理的原理和方法 ，并提供 了相应 的用 MicrosoftC 5．0写的源程序函数。 1 网格化数据等值线处理 固 l 等值姥在周秸 内的游动示意围 F曜．．1 Dimgrmmof~ontOllrlineswllkingwithinagrid 1．1 均 匀网格化数据等值线处理 均匀 网格化数据等值线处理的方法是对分布等距 1．2 非均 匀网格分布数据 的均匀用格化处理 的矩形 网格域 ，计算等值线 ， ，Y)一c与矩形网格每 实际工作中，许多采样点的分布是不规则的，所采 一 条边的交点坐标 ，再适当地将这些交点连接起来。一 数据很 大一部分是不均匀的，而且实际采样点的密度 般是先找到等值线 的一个起点 ，由此顺着等值线 一点 往往不能满足绘图要求．这就要求在绘制等值线之前 一 点地求 出交点，同时得到连接顺序。具体步骤如下。 进行数据 的均匀网格化处理和加密数据点．所 以我们 a)起点的寻找 首先找域 的外边框 ，看这些边上 采用曲面处理方法来插值、拟台数据。 有无起点 ，若没有则可 由底 向上或由左 向右逐步检索 1．2．1 按距离加权 的最小二乘法 (N P法) 起点 ； 按距离加权的方法，假定离网格 点越近的数据点 基盒项 目：中国地震局 九五”重点明 日(95 O4一O4—04) 收稿 日期 ：1001—03O5} 回 日期 2001～0808 作者筒彳r：李圣强 (1969- 、．男 ．福建松攫^ 。1992年毕业于福建省福，lI大学地矿系地质