地质资料处理：MapGIS 图件、钻孔柱状图与等值线
处理地质、矿产、水工环行业常见的三类工作：把 MapGIS 6.x 图件转成通用格式继续分析，按钻孔资料生成单孔柱状图，用钻孔或测点数据插值出等值线图。
你会得到什么
- MapGIS 图件转出的 Shapefile，或由 Shapefile 写回的 MapGIS 文件；
- 每个钻孔的柱状图 PNG；
- 插值格网（.asc、.grd 等）和等值线 Shapefile。
准备
| 需要 | 说明 |
|---|---|
| MapGIS6x 插件 | 内置，读写 WT / WL / WP，读取 MPJ 工程；不需要安装 MapGIS 软件 |
| ZealMap 地质钻孔柱状图插件 | 扩展插件，在插件市场安装；需要安装包含 MFC 的 Microsoft Visual C++ x64 运行库，见制图工具：钻孔柱状图 |
| Surfit 插件 | 扩展插件，在插件市场安装；散点、等值线、断层等多种约束一起插值，见数据处理与空间分析：Surfit |
MapGIS 图件转换
这样说：
看一下
D:\mapgis\地质图.mpj里有哪些图层、各是什么类型、有多少条记录；然后把引用的图层全部转成 Shapefile，中文属性不能乱码；转完统计一下地层面图层里每种岩性的面积。
智能体会怎么做：
- 只读探查 MPJ，列出引用的图层、类型、记录数、范围和文件内记录的投影；
- 把各图层批量导出为 Shapefile，同名不同类型的图层自动加后缀，转换需要审批，默认不覆盖已有文件；
- 需要分析字段内容时，再用 GDAL 查看字段、用 DuckDB 汇总统计。
注意：
- 文件里记录的投影只作线索，没有核实时标为未知，不会据此直接重投影或叠加底图。知道真实坐标系时请告诉它；
- 只转换二维几何，Z / M 值会丢弃；MapGIS 样式只能近似保留；
- 写回 MapGIS 时，精确的投影参数需要提供一个已有的 MapGIS 文件作为模板，例如“投影参数照
old.wp”； - 属性编码可以指定 auto、UTF-8、GBK，中文乱码时可以说“DBF 按 GBK 读”。
钻孔柱状图
这样说：
告诉我钻孔柱状图需要准备哪些表、各有哪些字段；我把 Excel 资料放在
D:\boreholes\原始资料.xlsx，按模板整理成四张 CSV，检查能不能出图；然后生成 ZK01 号孔的柱状图，纵向比例尺 1:500，制图人填“张工”。
智能体会怎么做：
- 从柱状图生成器读出数据模板：钻孔基本信息、原始分层、解释分层、物探四张表的文件名、字段、类型和中文名；
- 把你的资料按模板整理成四张 CSV。缺失的地层、深度、伽马测值不会推测补全，会列出来请你提供；
- 逐表核对文件、字段与列序，列出识别到的孔号；
- 按孔号生成整孔柱状图 PNG，默认 200 DPI，出图需要审批。
注意：
- 纵向比例尺分母 0 表示使用模板值；深孔太细长时取 500 或 1000 压缩图幅，多孔对比时统一取值；
- 深孔整孔 300 DPI 可能超出画布尺寸限制；
- 路径和中文文本需要能用 GBK 表示。
插值等值线
这样说：
用
D:\data\wells.csv（x、y 两列是坐标，z 列是地下水位）插值出 10 米格网，生成 1 米间距的等值线，格网输出 GeoTIFF，等值线输出 Shapefile，坐标系是 CGCS2000 3 度带第 38 带。
智能体会怎么做：
- 读取点数据，Shapefile、GeoJSON 中的点会先用 GDAL 导出成 X、Y、Z 三列文本；
- 用 Surfit 解算格网，可以同时加入已知等值线、断层线等约束；
- 追踪等值线，导出格网和等值线，执行需要审批；
- 用 GDAL 给成果补上坐标系、转成你要的格式，并裁掉数据范围以外的外推区域。
注意：
- Surfit 不处理坐标系，格网步长的单位就是数据坐标的单位，导出的文件需要另外补
.prj； - 格网会在数据外包框外多扩一圈，这部分是外推结果，误差明显更大，交付前应裁掉；
- 插件会拦下读到 0 个点、结果为空面等错误，避免交付整片无效值的结果；
- 百万级格网节点可能超出 Surfit 的内存上限，大范围时适当放大步长。
接着可以说
- “把转出来的地质图在 QGIS 里按岩性分类设色”（见在 QGIS / ArcGIS 里出图）
- “把所有钻孔的柱状图都生成出来”
- “结合 faults.shp 断层线一起插值，断层两侧不连续”
- “把等值线和钻孔点叠在一起做一张水位等值线图”（见制作专题地图）