干了这么多年环评生态专题,要说哪部分最磨人,制图绝对排前三。特别是土地利用现状图、植被类型图、植被覆盖度图、生态系统类型图、物种适宜生境分布图这五张图,几乎是生态影响评价的标配,也是专家评审时盯得最紧的部分。很多刚入行的同事觉得制图就是“画个图而已”,真上手才发现,坐标系统一、分类体系、解译精度、模型参数、色带设置,每一步都是坑。这篇就把我这些年做环评制图的流程、参数和踩坑记录整理出来,给正在做生态专题的朋友做个参考。
1. 制图前的数据准备与坐标系全局统一
1.1 数据来源与精度选择
环评制图的起点不是软件操作,而是数据。五类图件看似各自独立,实际上共用一套基础数据。我把常用数据源分成三类:
- 遥感影像数据:Landsat系列(30米分辨率)、高分一号/二号(2米~16米)、哨兵二号(10米)、无人机航拍(亚米级)。做区域尺度的环评专题,Landsat和哨兵二号是主力;做建设项目占地范围的小尺度制图,高分影像和无人机才是可靠选择。
- 基础地理数据:行政区划边界、道路、河流水系、居民点、等高线。这类数据通常来自测绘部门或公开地理数据服务,拿到手第一件事是核对边界是否与项目拟占地范围重叠。
- 野外调查数据:样方记录、样线调查、GPS定位点、植被调查表、物种分布记录。这部分数据精度直接决定植被类型图和物种生境图的“可信度”,调查点密度要覆盖所有典型生境类型。
提示:数据精度选择有个简单的判断标准——图件最终出图比例尺决定数据精度下限。1:1万出图至少要保证影像分辨率优于2米,1:5万出图用30米影像还能接受。比例尺定得高,数据跟不上,最后图上全是马赛克,评审一眼就能看出来。
1.2 坐标系与投影:最容易翻车的地方
坐标系是环评制图最基础、也最容易翻车的环节。我在实际工作中见过太多“图上位置看着差不多,叠加分析完全对不上”的情况,十有八九是坐标系问题。
环评图件常用的坐标系有两种:
- 地理坐标系:以经纬度定位,如CGCS2000、WGS84。适合大范围、跨带数据的存储和交换。
- 投影坐标系:把地球曲面展开为平面,如高斯-克吕格投影、UTM投影。环评图件最终出图、面积量算、空间分析都必须在投影坐标系下完成,否则面积会失真。
实际操作时,我一般遵循三条原则:
- 入库统一:所有原始数据不管来自哪里,统一用CGCS2000地理坐标系入库,记录元数据。
- 出图统一:制图开始时统一转成CGCS2000 / 高斯-克吕格投影(3度分带),中央经线按项目区所在带设定。
- 检查统一:每次叠加分析前,先加一层已知坐标的行政区界线验证,能重合再继续后续操作。
注意:投影分带不对是新手最常犯的错。比如项目区在东部某地,结果用了西部的中央经线,整张图会严重倾斜变形。检查方法很简单,把影像和矢量边界叠加,视觉上边界必须贴合地物轮廓,一旦出现“斜切”就要立刻检查投影参数。
2. 土地利用现状图制作:环评图件的“底图骨架”
2.1 分类体系怎么定
土地利用现状图是整个生态专题的底图,后续的植被图、生态系统类型图几乎都要基于它叠加分析。分类体系的选择直接影响解译工作量和后期分析的准确性。
目前环评中应用最普遍的是国家《土地利用现状分类》标准,该标准将土地划分为耕地、园地、林地、草地、商服用地、工矿仓储用地、住宅用地、公共管理与公共服务用地、特殊用地、交通运输用地、水域及水利设施用地、其他土地等一级类,每个一级类下面细分二级类。环评专题制图时,我建议根据项目性质决定分类深度:
- 线性工程(公路、管线):用二级类即可,重点关注耕地、林地、草地的面积占用和切割。
- 区域开发(工业园区、新城建设):建议用二级类甚至三级类,因为需要做占地类型明细统计。
- 生态影响为主的矿山、水电项目:重点关注林地和草地,分类上要把天然林、人工林、灌木林地分开。
2.2 解译标志建立与人工目视解译
分类体系确定后,最核心的工作是解译。我见过有人拿到影像就开画,画到一半才发现同一类地物在不同影像上色调差很远。正确做法是先建立解译标志库。
具体操作是这样:
- 选取样区:在影像覆盖范围内均匀选取若干样区,包括已知的耕地、林地、居民点、水体等地物。
- 记录特征:对每个样区的影像特征进行描述,包括色调、形状、纹理、阴影、位置、布局,并附上影像截图和实地照片。
- 形成标志库:把这些特征整理成一个表格或文档,作为人工目视解译的统一标准。
以我常用的30米Landsat假彩色影像为例(波段组合:近红外、红、绿),水体和湿地呈深蓝黑色,乔木林地呈深红色、纹理粗糙,灌木林地呈浅红褐色、纹理较细,草地呈浅红色到黄色、纹理均匀,耕地呈规则几何块状、颜色随作物季节变化,居民点呈灰白色、边界清晰。
解译时一幅影像我通常要过两遍:第一遍粗略勾勒大类,第二遍细化二级类。这个过程快不了,1:5万出图的小区域项目,人工解译2~3天很正常。
实操心得:解译时不要把每个图斑都抠得死精细。环评图件的核心作用是反映土地利用的空间格局和占地关系,不是测绘地籍图。边界线平滑、拓扑正确、分类准确,比边界位置精确到米更重要。
2.3 制图表达细节
土地利用现状图最后呈现效果好坏,制图表达占一半。
最小图斑面积:航片解译时小于0.15公顷的图斑一般不单独上图(根据不同规范有所调整),并入相邻图斑。影像分辨率越高,最小图斑可以越小。否则图面上全是碎斑,看着像撒了一地芝麻。
颜色配置:土地利用现状图的用色,我建议优先参考国家标准图式规定的颜色,比如耕地常用浅黄色、林地常用绿色系、草地常用浅绿色、水域用蓝色、建设用地用粉红色。同一张图内色相要拉开差距,避免相近颜色相邻导致边界难以辨认。饱和度不要过高,印刷出来会偏色。
边界处理:图斑边界必须闭合、不能有缝隙和重叠,拓扑关系要提前检查。弧段相交、悬挂点这些问题,在GIS软件里一键修复后还要人工抽查。
3. 植被类型图与植被覆盖度图:遥感反演与野外验证
3.1 植被类型图制作流程
植被类型图不是画出来的,是“调查+判读+验证”三遍工序出来的。
分类体系用的是植被生态学的经典方法,从高级到低级依次为植被型组、植被型、群系。环评项目中,一般做到群系层级就够用,比如“暖性针叶林—马尾松林”“落叶阔叶林—栎林”这样的层级。分类过细会导致图面破碎,分类过粗又满足不了生态影响分析的需求。
制作流程分四步:
- 野外调查:沿调查线路布设样方和样线,记录植物种类、群落结构、郁闭度/盖度、GPS坐标。调查点要覆盖所有植被类型,每个类型至少3个以上重复样本。
- 影像判读:结合高分辨率影像的色调、纹理、地形位置,把调查点信息推广到整个研究区。比如阳坡山脊的深红色纹理区域大概率是灌丛,沟谷地带的深绿色高大纹理区域一般是阔叶林。
- 专家知识修正:把判读结果与区域植被分布规律对照。比如某个海拔带理论上应该是针叶林,判读结果却大面积是灌丛,就需要检查是放牧干扰还是判读错误。
- 野外验证:选取判读结果的边界区、可疑区进行野外抽查,用GPS点位验证类型判对率。验证点数量不宜低于调查点数量的30%。
3.2 植被覆盖度:NDVI像元二分模型
植被覆盖度图在环评里主要用于水土流失评价、生态系统质量评估和施工迹地恢复监测。制作方法很多,最常用、最稳定的是NDVI像元二分模型。
原理不复杂,每个像元的植被覆盖度FVC由该像元NDVI值在纯裸土NDVI和纯植被NDVI之间的相对位置决定,公式如下:
$$FVC = \frac{NDVI - NDVI_{soil}}{NDVI_{veg} - NDVI_{soil}}$$
其中NDVI是归一化植被指数,由近红外波段和红波段的反射率计算得出:
$$NDVI = \frac{NIR - Red}{NIR + Red}$$
实际操作中关键在参数选取:
- NDVI_soil:取研究区内裸土、建设用地、无植被覆盖像元的NDVI值,一般取5%分位数。
- NDVI_veg:取研究区内茂密植被像元的NDVI值,一般取95%分位数。
不用全影像的最大值和最小值,而是用分位数,目的是避免云影、水体、噪声像元导致极值失真。
注意:水体在NDVI计算中通常为负值或极低值,会干扰NDVI_soil的统计。计算前建议先用水体掩膜把河流、水库、坑塘剔掉,算完再把水体区域单独赋予0值或保留原始像元。
计算完成后,按环评习惯把覆盖度分为裸地(<10%)、低覆盖(10%~30%)、中覆盖(30%~50%)、高覆盖(50%~70%)、极高覆盖(>70%)五级,再转成专题图。色带从浅黄到深绿渐变,裸地用棕黄色系,视觉上符合“越绿覆盖越好”的直觉。
4. 生态系统类型图制作:从“类型”到“系统”的整合
4.1 生态系统分类的逻辑
很多新手分不清土地利用现状图和生态系统类型图的区别。简单说,土地利用现状图的分类核心是“地表覆盖和利用方式”,生态系统类型图的分类核心是“生物群落及其环境构成的生态功能单元”。
最直观的差别在“水域”上:土地利用分类里河流是“水域及水利设施用地”,而生态系统类型图会把河流区分为河流湿地生态系统,如果河边有滩涂,还会单独分出草本沼泽或灌丛沼泽。农田在土地利用图上是耕地,在生态系统图上则分成农田生态系统(包括水田和旱地两个子类)。
我常用的生态系统分类体系是:
- 森林生态系统(常绿阔叶林、落叶阔叶林、针叶林、针阔混交林、竹林)
- 灌丛生态系统(常绿灌丛、落叶灌丛)
- 草地生态系统(温性草原、热性灌草丛、山地草甸)
- 湿地生态系统(河流湿地、湖泊湿地、沼泽湿地、人工湿地)
- 农田生态系统(水田、旱地)
- 城镇生态系统(居住地、工业用地、交通用地)
- 荒漠生态系统(沙地、裸岩石砾地)
4.2 数据整合策略
生态系统类型图的制作,核心思路是“多源数据整合、按生态学逻辑归并”。
我通常这样操作:以土地利用现状图的图斑为空间框架,叠加植被类型图的群落信息,再结合地形数据(海拔、坡度、坡向)和野外调查数据,逐图斑判定生态系统类型。比如一个图斑在土地利用图中是“有林地”,叠加植被类型后确认是马尾松林,地形上处于低海拔丘陵,就归为森林生态系统中的常绿针叶林亚类。
这个过程看着简单,实际操作有四个容易出问题的点:
- 图斑边界不匹配:土地利用图斑和植被类型图斑边界不完全重合,需要设置一个容差范围,通常不超过一个像元大小。
- 分类冲突:同一边界内土壤和植被信息矛盾时,以野外实际调查记录为准。
- 最小面积控制:生态系统图最小图斑面积一般比土地利用图大,建议不小于1公顷,过小的图斑归并到相邻优势类型。
- 过渡带处理:森林-草地过渡带、湿地-农田过渡带是图斑犬牙交错的重灾区,处理原则是“类型之间保留渐变,不强行画直线边界”。
5. 物种适宜生境分布图:模型与图件的结合
5.1 物种数据与环境变量
物种适宜生境分布图是生态影响评价中预测性最强的图件。它的核心问题不是“现在哪里有什么”,而是“哪些区域具备某物种生存繁衍的条件,工程活动落在哪里会造成影响”。
做图的起点是物种分布数据和环境变量数据。
物种数据来源包括野外调查记录、历史监测资料、文献记录。数据质量要求有三个:一是坐标精度要可靠,至少到秒级;二是要覆盖物种在研究区内的主要分布范围,而不是只在项目区周边采几个点;三是样本量太小时(少于10个有效点),不建议直接建模型,可以用专家知识规则替代。
环境变量一般分四类:
- 气候变量:年均温、最冷季均温、年降水量、降水季节性
- 地形变量:海拔、坡度、坡向、地形湿度指数
- 植被变量:植被类型、植被覆盖度、距森林边缘距离
- 干扰变量:距道路距离、距居民点距离、土地利用类型
5.2 适宜生境模型的常见做法
环评中用的较成熟的模型方法是最大熵模型(Maxent),因为它只需要物种“存在点”数据,对样本量要求相对宽松。另有逻辑回归、随机森林、集成模型等需要“存在+不存在”数据的方法,在数据条件满足时精度也很高。
建模流程大致是:
- 变量筛选:把候选环境变量先做相关性分析,两两相关性高于0.8的变量只保留生态学意义更明确的一个,避免多重共线性。
- 模型训练:用存在点数据和环境变量训练模型,预留25%~30%的点做验证。
- 预测与评价:模型输出每个栅格的适宜性指数(0~1),用AUC值评价模型表现。AUC大于0.8说明模型区分能力良好,可接受。
- 阈值划分:用等敏感度-特异性阈值或10分位训练存在阈值把连续适宜性指数切分为“适宜/不适宜”二值图,再进一步细分为高、中、低适宜区。
提示:阈值选择非常影响最终图面。阈值过低会把大片区域划为适宜生境,图件失去评价意义;阈值过高又会漏掉重要栖息地。我习惯同时输出三个等级(高适宜、中适宜、低适宜),用不同颜色叠加显示,而不是只给二值图。评审专家更愿意看到这种连续变化的表达,因为生境适宜性本身就是渐变的。
5.3 制图表达与不确定性
物种适宜生境分布图最常见的制图问题有两个:一是把模型输出的栅格图直接截图摆放,没有经过任何制图规范处理;二是结果图没有标注模型精度和不确定性,被质疑时拿不出解释。
我的建议是:
- 栅格结果重分类后转矢量,按适宜性等级做面状表达,既方便面积统计,也符合环评制图习惯。
- 图面必须附注:模型方法、样本量、AUC值、使用的环境变量、阈值和分级标准。标注这些不是为了凑字数,是为了让评审能判断图件的可靠性。
- 高适宜区与工程占地范围做叠加分析,计算重叠面积占比,在图上用重点符号标识出“工程影响重点关注区域”。
6. 图件质量控制与常见返工坑
6.1 制图规范自查清单
环评图件的质量检查,我总结了一套“出图前10分钟自查法”,每次都按这个顺序过一遍:
- 图幅是否包含完整图框、图名、图例、比例尺、指北针、坐标系说明?
- 图名是否与报告章节对应,字号是否明显?
- 图例是否覆盖图中所有色块和图斑,图例顺序是否与图面内容对应?
- 比例尺与图幅尺寸是否匹配,1:5万图不能配一个只有1厘米长的比例尺?
- 指北针方向是否与图面实际方向一致?
- 经纬网或公里网格是否清晰,密度是否合适?
- 图面是否存在缺字、压盖、边界断线?
- 导出分辨率是否达到印刷要求(不低于300dpi)?
6.2 高频问题与排查技巧
| 高频问题 | 典型表现 | 排查思路 |
|---|---|---|
| 图层错位 | 植被图斑与影像地物明显偏移 | 检查坐标系是否统一,投影带是否正确 |
| 图例与图斑不对应 | 图例上有林地,图上找不到对应颜色 | 检查分类字段是否存在空值或错误赋值 |
| 面积与报告不一致 | 报告写林地面积500公顷,图上统计700公顷 | 检查是否统计了非项目区的图斑,是否未扣除水域 |
| 边界锯齿严重 | 图斑边界沿像元呈锯齿状 | 需要做边界平滑处理,但注意不能过度平滑导致边界偏移 |
| 出图模糊 | 放大后明显马赛克 | 导出设置分辨率过低,或原始数据精度不足 |
排查时有个实用技巧:把图件导出的矢量数据重新叠加到原始影像上,随机抽10个图斑点核对分类属性。如果错误率超过5%,这张图就得重新解译了。
最后再分享一个我在实际项目中反复体会到的经验:环评制图这份工作,技术门槛真的不算高,但非常磨细致。五张图从数据收集到最终出图,环环相扣,前面一个坐标系的小疏忽,往往要到叠加分析时才会爆炸,那时候返工成本就是按周算了。所以做这套图,我的习惯是“慢前面、快后面”——数据检查和坐标系统一阶段多花三天,后面解译、建模、出图反而顺畅。项目周期再紧,这三天不能省。希望能帮你少踩几个坑。