2. 矿区多源地质数据准备与QGIS工程初始化
剖面线布设与剖面数据提取的精度,首先取决于进入QGIS的各类资料是否"同源同框"。矿区资料通常来自测绘、地质、物探、钻探等不同部门,格式、比例尺、坐标系、高程基准往往互不一致:地形图可能是CGCS2000的DWG,地质图是北京54的扫描栅格,钻孔是Excel坐标表,物探剖面则是带测点桩号的文本或图片。若直接把这些数据拖入工程,会出现图层错位、剖面线穿错地层、高程值偏差数十米等问题。因此本节的任务不是"打开软件",而是建立一套坐标统一、属性可追溯、可复用的工程底座。
一、资料清点与格式归类
动手前先做一次资料台账,把每份数据的关键元信息记录下来,这是后续判断坐标是否可信的依据。建议至少记录:数据名称、来源部门、格式、比例尺/精度、坐标系与高程基准、覆盖范围、时效。以某铅锌矿区教学示例为例,可整理为下表(示例数据,仅用于说明字段组织方式):
| 资料 | 格式 | 坐标系/基准 | 用途 |
|---|---|---|---|
| 1:2000地形图 | DWG | CGCS2000 / 1985国家高程 | 底图、等高线、剖面高程 |
| 1:10000地质图 | 扫描TIFF+配准点 | 北京54 / 1956黄海高程 | 地层界线、断层 |
| 钻孔台账 | XLSX | 含X、Y、孔口高程、孔深 | 地层控制、剖面约束 |
| 物探剖面 | CSV(桩号、测点坐标、视电阻率) | 独立桩号系统 | 推断构造、异常体 |
归类后要明确两类问题:一是坐标系不一致(如北京54与CGCS2000),二是高程基准不一致(1956黄海与1985国家高程,两者在局部地区可差数十厘米至一米量级,需按矿区实测转换参数处理,不可随意忽略)。
二、建立统一坐标系的QGIS工程
启动QGIS后,先通过"工程 → 属性 → CRS"设定工程坐标系。矿区大比例尺制图一般选用以矿区中央经线为准的高斯-克吕格投影带,例如CGCS2000 3度带第34带(中央经线102°E),EPSG代码需按实际带号选取。工程CRS的作用是"显示与计算基准",它不会自动改变源数据的坐标值,因此源数据本身仍需转换或正确声明。
对矢量数据,优先使用"矢量 → 数据管理工具 → 重投影图层"生成新文件,而不是仅修改图层CRS声明——后者只改标签不改坐标,是错位的最常见原因。对扫描地质图,使用"栅格 → 配准"工具,选取至少4个分布均匀的控制点(建议图廓角点+中部交叉点),变换类型选一阶多项式,残差控制在1~2个像元内再输出GeoTIFF。钻孔表通过"图层 → 添加图层 → 添加分隔文本图层"导入,X字段填东坐标、Y字段填北坐标,CRS指定为台账所记坐标系,导入后用"导出 → 要素另存为"落到工程CRS。
三、工程属性与图层组织
在"工程 → 属性"中填写标题、作者、摘要,并把坐标系、高程基准、资料来源写入元数据,便于多人协作时追溯。工程设置里建议统一长度单位为米、面积单位为平方米,坐标显示小数位设为3位。图层按"底图—地质—勘探—物探—剖面设计"分组,命名采用"类型_比例尺_来源"规则,如"地质_1万_扫描配准"。同时为钻孔图层设置标注(孔号)与符号(空心圆),为后续剖面线选取提供直观参照。
四、初始化后的自检与固化
工程搭好后做三项自检:其一,用"测量"工具量取地形图上已知边长,与图上注记比对,验证配准精度;其二,打开钻孔属性表,检查孔口高程与地形等高线是否吻合,若系统性偏高或偏低,说明高程基准未统一;其三,缩放到全图范围,确认所有图层覆盖范围一致、无整体偏移。自检通过后,将工程另存为.qgz,并把转换后的中间数据统一放入工程目录的data子文件夹,使用相对路径保存,避免换机后图层丢失。
完成上述工作后,工程中已具备统一坐标、统一高程基准、分层可管理的底图与勘探数据,剖面线布设时即可直接在地形与地质底图上拾取端点、读取里程与高程,并为后续剖面数据提取提供可靠的坐标与属性基础。