3. 找矿预测因子栅格构建

在完成地质构造线、地层界线、岩体接触带等矢量要素的整理之后，下一步并不是直接进入重分类，而是要把这些“线状、面状、边界状”的矢量信息，转换成栅格空间中能够表达“离控矿要素有多远、有多密、朝哪个方向”的连续数值图层。这一步决定了后续重分类的输入质量：如果距离栅格的空间参考、像元大小或搜索半径设置不当，即使重分类阈值再精细，预测结果也会失真。本节围绕三类预测因子栅格——距离栅格、密度栅格与方向栅格——说明它们的物理含义、构建方法与参数选择，并给出可复现的操作思路。

一、三类因子栅格分别表达什么

距离栅格回答的是“每个像元到最近控矿要素的距离”，常用于表达断裂、岩体接触带、地层界线对成矿的控制具有随距离衰减的特征。密度栅格回答的是“单位面积内要素的发育程度”，例如单位面积内断裂总长度（线密度）或断裂交点数（点密度），用于刻画构造活动的强弱。方向栅格则回答“要素在像元附近的主导走向”，用于区分不同走向断裂与区域应力场、地层走向的配套关系。三者共同构成“远近—强弱—方位”的因子体系，而不是彼此替代。

二、距离栅格：欧氏距离的构建与参数控制

距离栅格最常用的实现是欧氏距离工具。其原理是：对每个输出像元，计算它到最近源要素（线或面边界）的直线距离，并把该距离写入像元值。构建时有三项参数必须明确。第一是像元大小，它应与后续重分类和证据权计算的分辨率一致，例如研究区为 1:5 万比例尺时，常取 25 m 或 30 m；若源数据精度只有 1:20 万，强行用 10 m 像元只会制造虚假精度。第二是最大搜索距离，超过该距离的像元统一赋为最大值或 NoData，避免远景区外围被无意义的大距离值干扰；该值应结合已知矿点到断裂的最大影响距离设定，例如示例中取 3000 m。第三是投影坐标系，距离必须以米为单位计算，若数据为地理坐标系（度），必须先投影转换，否则“度”与“米”混用会使距离因子完全失去物理意义。

实际操作中，可先对构造线做“要素转栅格”得到源栅格，再执行欧氏距离；也可直接以矢量线作为输入。需要注意，线要素的端点、相交节点都会成为距离计算的源，若希望只统计主干断裂，应先按长度或属性筛选，避免短小派生断裂拉低整体距离值。输出后建议用直方图检查距离分布：若大量像元集中在 0–200 m，说明源要素过密或搜索距离过小；若分布过于平缓，则可能搜索距离过大。

三、密度栅格：线密度与点密度的统计口径

密度栅格的核心是“单位面积内的要素量”。线密度常用“线长度/面积”（km/km²），点密度常用“点数/面积”（个/km²）。ArcGIS 中可通过线密度工具直接由折线生成，也可先用“相交”或“近邻”统计每个像元邻域内的要素量再除以邻域面积。关键参数是搜索半径：半径过小，密度图呈斑点状，反映的是局部噪声；半径过大，则不同构造带被平滑成一片，失去区分度。示例中若研究区断裂带宽度约 2–5 km，可取 5 km 搜索半径，使密度值在构造带中心高、两侧递减。

密度因子还需注意统计口径的一致性：若用线密度，所有构造线应使用同一投影和同一最小长度阈值；若用点密度（如断裂交点、岩体接触点），应明确“点”的定义，避免把同一条断裂的多个折点误当交点。密度值通常为连续浮点数，重分类前可先做分位数统计（如 20%、40%、60%、80% 分位）确定分级断点，而不是凭经验硬切。

四、方向栅格：走向统计与角度约定

方向栅格不像距离和密度那样有现成的一键工具，通常需要先计算每条线要素的走向方位角，再在邻域内做方向统计。角度约定必须先在全文统一：本文采用正北为 0°、顺时针增加、走向取 0°–180°（因为走向不分正反）。计算时可用字段计算器对折线起止点求反正切，得到方位角后归一化到 0°–180°。随后以像元为中心开窗口（如 3×3 或 5×5），统计窗口内要素走向的众数或加权平均，作为该像元的方向值。

方向栅格的价值在于区分“与区域主构造平行”和“与主构造斜交”的断裂。例如区域主应力方向为 NE 向时，NE 向断裂往往为控矿断裂，而 NW 向断裂可能为后期破矿构造；若只做距离和密度，二者会被混为一谈。方向栅格输出后，可用“重分类”把 0°–180° 划分为若干方位区间（如 NE 45°–75°、NW 135°–165°），再与距离、密度因子组合。

五、构建流程与质量检查要点

建议按“统一投影 → 筛选要素 → 分别生成三类栅格 → 对齐像元 → 检查分布”的顺序推进。三类栅格必须使用同一像元大小、同一范围、同一投影，否则后续栅格计算器叠加时会出现错位。生成后逐层检查：距离栅格看最小值是否为 0、最大值是否等于设定搜索距离；密度栅格看高值区是否与已知构造带吻合；方向栅格看是否存在大量空值（窗口内无要素时需定义填充规则）。只有三类因子在空间上严格对齐、数值分布合理，才具备进入重分类与证据权计算的条件。