2. 磁法数据准备与质量控制

磁法网格化的质量上限，在数据进入 Oasis montaj 的“Gridding”对话框之前就已经被决定了。网格化算法只能忠实地把离散测点插值成规则网格，它无法识别一个坐标被误录成 116.3° 而不是 116.03° 的点，也无法自动剔除同一位置重复采集造成的虚假高值。因此，本节把数据准备视为一道独立的质控工序：先确认坐标与磁异常值本身可信，再处理重复点与缺失值，最后用裁剪和投影把数据限定在正确的空间参照下。

一、测点坐标检查：先看参照系，再看数值范围

坐标检查的第一步不是看数值大小，而是确认坐标系定义。磁法测点常见的坐标形式有两类：地理坐标（经纬度，单位度或度分秒）和投影平面坐标（如高斯-克吕格 3° 带、UTM，单位米）。Oasis montaj 的数据库对 X、Y 字段本身不携带坐标系语义，坐标系信息保存在工程或地图的投影设置中。因此必须保证“字段里的数值”与“工程声明的投影”一致：若 X 字段存的是经度 116.35、Y 字段存的是纬度 39.92，而工程按投影坐标（米）解释，网格范围会瞬间变成几十米的微小区域，或者直接报错。

数值范围检查可以用统计量快速完成。对 X、Y 分别计算最小值、最大值、均值与标准差，并绘制 X-Y 散点图。正常测网应呈现规则的测线—测点结构：散点图上是等间距的平行线束，或规则网格点阵。若出现少量远离主体的孤立点，多半是坐标录入错误（多一位、少一位、小数点错位、经纬度与投影坐标混用）。教学示例中，某测区 X 应在 500000～512000 m、Y 应在 4400000～4412000 m 之间，若发现 X=51200 的点，即可判定为漏写一位，应回到原始记录核对而非直接删除。

磁异常值检查关注三件事：量纲与单位、异常范围、以及异常与坐标的对应关系。总场磁异常 ΔT 通常以 nT 为单位，幅值从几 nT 到数千 nT 不等；若字段中出现 0.00005 这类数值，往往是单位换算成了 T 而未乘 10⁹。用直方图观察分布形态，正常磁异常多为近似正态或右偏分布，若出现极端离群值（如 ±100000 nT），应核查是否为仪器饱和、探头碰撞或数据拼接错误。此外，把磁异常值作为颜色绘制在 X-Y 散点图上，可以直观发现“整条测线数值恒定”“相邻测点跳变数百 nT”等采集问题。

二、重复点与缺失值的处理原则

重复点指坐标完全相同或极其接近（例如距离小于 1 m）的多个测点。磁法野外作业中，重复点可能来自复测、GPS 漂移或数据合并时的重复导入。重复点对网格化的影响是局部的但很典型：同一位置存在两个不同 ΔT 值时，插值算法会在该处产生一个虚假的尖峰或“牛眼”异常，误导后续解释。处理方式取决于重复点的来源：若是同一测点的多次观测，取算术平均值最合理；若其中一次明显受干扰（如探头碰撞导致跳变），则按数据质量优先级保留可信值并记录剔除理由。Oasis montaj 中可用数据库的“Find Duplicate”类工具按 X、Y 组合查找重复记录，再人工确认。

缺失值包括坐标缺失和磁异常值缺失两类。坐标缺失的点无法参与网格化，只能剔除；磁异常值缺失的点，若坐标完整，可以保留坐标而将 ΔT 置为空（Dummy），让网格化算法按邻域插值填补，但必须在报告中说明空值比例。经验上，若某测线空值比例超过 10%～15%，该测线的插值结果可信度明显下降，应考虑整线剔除或补测。处理缺失值时切忌用 0 填充——0 nT 在磁异常中是一个有物理含义的数值，会被当作真实低值参与插值，造成虚假的负异常条带。

三、数据裁剪与投影设置

裁剪的目的是把数据限定在有效测区范围内，同时为网格化留出合理的边缘缓冲。裁剪范围通常取测点包络矩形再向外扩展 1～2 个点距（示例中约 20～40 m），这样网格边缘的插值不会因邻域点不足而失真。裁剪过紧会把边缘测点切掉，裁剪过松则引入大片无数据区，网格化时产生外推假象。在 Oasis montaj 中可通过坐标范围选择或地图窗口交互裁剪，裁剪后应核对剩余点数与预期是否一致。

投影设置是容易被忽视却影响全局的一步。需要明确三件事：投影类型（如横轴墨卡托、高斯-克吕格）、中央经线与假东偏移、以及坐标单位（米）。若测区跨带，应统一到同一带内或改用地理坐标并说明。投影设置完成后，检查地图窗口中测点分布的形状是否与实地测网一致——若测网被明显拉伸或旋转，说明投影参数有误。最后，把 X、Y、ΔT 三个字段整理为网格化所需的干净通道，并保存一份处理前的原始数据备份与处理日志，保证整个质控过程可追溯、可复核。