4. 环境空间数据获取与格式转换

环境影响专题图的质量，首先取决于“底料”是否干净、来源是否可靠、坐标是否一致。本节围绕三条常见数据来源——公开环境数据库、地面监测站点、遥感影像——说明如何获取数据，并在 QGIS 中完成矢量与栅格格式转换、坐标系统一，为后续的缓冲区分析、叠加分析与制图出图打好基础。

一、三类数据源的特征与获取方式

公开环境数据库通常以矢量为主，例如行政区划、水系、自然保护区边界、排污口点位、土地利用分类等，常见格式为 Shapefile、GeoJSON、GeoPackage 或 CSV（含经纬度字段）。获取时优先选择带明确坐标系说明、带属性字段字典、带更新时间的来源；若只有 CSV，需要确认经纬度字段的列名、单位（度还是度分秒）与坐标顺序（经度在前还是纬度在前），这是后续最容易出错的地方。

地面监测站点数据多为“点位 + 时间序列”的组合：站点位置是一张点表，污染物浓度（如 PM₂.₅、SO₂、COD、氨氮）是另一张按小时或按日记录的表，两者通过站点编号关联。获取后应先在表格软件中检查站点编号是否唯一、是否存在重复站点、坐标是否落在合理范围内（例如经度 73°～135°、纬度 3°～54° 大致对应中国范围），再导入 QGIS。

遥感影像常见为 GeoTIFF，可能附带 .tfw 世界文件或已内嵌地理参考。获取时要关注空间分辨率、成像时间、云量、投影方式（如 UTM 分带或 WGS 84 地理坐标）以及波段顺序。若影像没有坐标信息，需要借助控制点或已知范围进行地理配准，否则后续无法与其他图层对齐。

二、在 QGIS 中把 CSV 站点表变成点图层

以“某市 2023 年国控站点 PM₂.₅ 年均浓度”为例（示例数据，仅用于演示流程）：CSV 含 station_id, station_name, lon, lat, pm25 五列。在 QGIS 中通过“图层 → 添加图层 → 添加分隔文本图层”打开文件，分隔符选逗号，几何图形选“点坐标”，X 字段选 lon，Y 字段选 lat，坐标系先选 EPSG:4326（WGS 84 地理坐标）。导入后右键图层“属性 → 源”确认坐标系，再用“识别要素”工具点几个站点，核对坐标是否落在正确城市范围内。若点位整体偏移到海上或邻省，多半是经纬度顺序颠倒或单位错误，应回到 CSV 修正，而不是在 QGIS 里硬调。

导入成功后建议立即“另存为”GeoPackage（.gpkg），把点图层固化下来。GeoPackage 单文件、支持中文、字段长度不受 Shapefile 的 10 字符限制，适合作为中间成果的容器。

三、矢量与栅格的相互转换

矢量转栅格最典型的场景是把土地利用面转成栅格，用于与遥感影像做叠加或参与栅格计算。使用“栅格 → 转换 → 栅格化（矢量转栅格）”，关键参数有三个：字段选择（分类面选类别编码字段，如 lu_code）、输出栅格分辨率（应与目标遥感影像一致，例如 10 m）、输出范围（选“使用图层范围”并指定参考影像）。分类数据必须保证编码为整数，否则栅格化后类别会混乱。

栅格转矢量常见于把遥感反演的浓度栅格提取为等值面。使用“栅格 → 转换 → 多边形化（栅格转矢量）”，字段选波段值，得到按像元值分区的面。若栅格是连续型（如 PM₂.₅ 浓度），多边形化前应先“重分类”为若干等级（如 0–35、35–75、75–115 μg/m³），否则会生成大量碎小多边形，既慢又难读。

四、坐标统一的实操要点

坐标统一的核心是“先定目标，再逐个重投影，最后验证”。环境影响专题图若面向国内管理与制图，常用 CGCS2000 地理坐标（EPSG:4490）或 CGCS2000 3° 带高斯-克吕格投影（如 EPSG:4547，中央经线 114°E）。地理坐标适合全国尺度，投影坐标适合市域尺度且能直接量算面积与距离。

矢量重投影用“矢量 → 数据管理工具 → 重投影图层”，指定目标 CRS 后另存为新文件；栅格重投影用“栅格 → 投影 → 变形（重投影）”，重采样方法按数据类型选择：连续型（高程、浓度、温度）用双线性或三次卷积，分类型（土地利用、覆盖类型）必须用最近邻，避免产生不存在的类别值。项目本身也应在“项目 → 属性 → 坐标系”中设为同一目标 CRS，这样所有图层在同一画布下才能正确叠加。

五、结果验证与常见问题

验证坐标是否真正统一，可用三种方法：一是叠加比对，把重投影后的站点与行政区划叠加，看点位是否落在正确区县；二是量测校验，用测量工具量一段已知长度的道路或已知面积的行政区，与官方公布值比对，偏差应在合理范围内；三是范围检查，查看图层范围是否与目标区域吻合，若出现跨半球或超大范围，说明 CRS 设置有误。

常见问题包括：CSV 经纬度顺序颠倒导致点位偏移；Shapefile 中文乱码（应设置编码为 UTF-8 或 GBK）；栅格重采样方法选错导致分类值失真；不同来源数据基准面不一致（如 WGS 84 与 CGCS2000 差异虽小，但在高精度制图中仍需统一）。处理原则是：先查源数据说明，再在 QGIS 中做最小必要转换，每一步另存新文件并记录参数，保证过程可追溯、结果可复现。