做生态影响评价这一行的人都有个共识报告写得再漂亮评审会上专家翻到附图那几页心里基本就有数了。项目区位置图、土地利用现状图、地表水系图、植被类型图、植被覆盖度图、归一化植被指数图、生态系统类型图、土壤侵蚀图、生物量分布图——这九张图不是凑数的装饰品它们是整个生态影响评价的空间证据链。位置图回答在哪土地利用和生态系统类型图回答是什么水系图回答水怎么走植被三件套回答绿得怎么样、长得好不好土壤侵蚀图回答土稳不稳生物量图回答生态家底有多厚。九张图串起来就是一套完整的生态现状空间表达体系。我这些年做过不少这类专题图件从最初用GIS软件一张张硬抠到后来搭模板、写批处理、建符号库踩过的坑能写满一个笔记本。这篇文章就把这套图件的制作逻辑、数据来源、关键参数、实操流程和避坑经验完整拆一遍适合刚接手生态评价图件的新人照着做也适合做了几年但总觉得图不够专业的老手对照查漏。核心关键词就围绕生态影响评价、土地利用现状、植被覆盖度、归一化植被指数、土壤侵蚀、生物量分布这些展开不绕弯子直接上干货。1. 九张图件的整体设计与制图逻辑拆解1.1 先想清楚图给谁看再决定图怎么画很多人一上来就打开软件开始描图这是最容易返工的做法。生态影响评价的附图有两个明确的读者一是评审专家二是审批部门的技术审查人员。专家看的是逻辑是否自洽、分类是否规范、结论是否有数据支撑审查人员看的是图件是否齐全、要素是否完整、图例坐标是否规范。这两类读者的关注点不同决定了制图时科学性和规范性必须同时满足。所以我的习惯是动手前先把九张图分成三组来理解。第一组是本底组项目区位置图、土地利用现状图、生态系统类型图这三张回答项目落在什么环境里。第二组是植被组植被类型图、植被覆盖度图、归一化植被指数图、生物量分布图这四张回答植被质量如何。第三组是过程组地表水系图、土壤侵蚀图这两张回答生态过程是否稳定。分组的意义在于同组图件的数据源、分类体系、配色逻辑要相互呼应跨组之间要有承接关系。比如土地利用现状图的分类要能映射到生态系统类型图植被类型图的边界要和生物量分布图的高值区对得上否则专家一眼就能看出图件是各画各的。提示九张图不是九个独立任务而是一个数据体系的不同切面。制图前先画一张数据关系草图标明哪张图用哪个数据源、哪两张图之间需要保持边界一致能省掉后期大量返工。1.2 坐标系、投影与制图范围的三条硬规矩生态评价图件最容易在基础设置上翻车。第一条规矩所有图件必须使用统一的坐标系和投影。国内项目通常用CGCS2000地理坐标系或对应的3度带高斯-克吕格投影面积量算类图件土地利用、生态系统类型、土壤侵蚀建议用投影坐标系因为地理坐标系下直接算面积会有变形误差。第二条规矩制图范围要一致。九张图的外框范围应当统一通常以项目区为中心外扩一定缓冲距离缓冲距离根据生态影响评价等级确定一级评价范围大、二级三级依次收窄具体以评价技术导则和项目实际影响半径为准。第三条规矩比例尺要匹配精度。位置图可以用小比例尺概览但土地利用、植被类型这类专题图的比例尺要和数据源精度匹配用30米分辨率的遥感数据硬出1:10000的图边界锯齿会非常明显。我见过一个典型问题某项目的土地利用现状图用的是地理坐标系面积统计出来和报告正文对不上差了将近百分之三。排查半天才发现是投影设置的问题。所以我的做法是建一个统一的工程模板坐标系、范围、比例尺、图框、图例位置全部预设好九张图都在这个模板里出从源头上杜绝不一致。1.3 分类体系的选择跟着导则走别自创门派土地利用现状图的分类主流做法是参照国家标准土地利用现状分类体系结合项目区实际归并到评价需要的层级。生态系统类型图则通常参照生态系统分类的相关标准把土地利用类型映射到森林、草地、湿地、农田、城镇等生态系统大类。这里的关键是映射关系要写清楚报告里要有一张对照表说明土地利用的哪些二级类归并到了哪个生态系统类型否则两张图的面积对不上专家必然追问。植被类型图的分类要结合项目区实际的植被区划用建群种或优势种命名比如常绿阔叶林落叶阔叶灌丛山地草甸这类。土壤侵蚀图则按水力侵蚀、风力侵蚀等类型和强度分级来表达分级标准参照土壤侵蚀分类分级标准。归一化植被指数图和植被覆盖度图是定量图前者直接反映NDVI值分布后者通常由NDVI通过像元二分模型转换而来。生物量分布图则是估算图数据来源可以是实测样方外推、经验模型估算或遥感反演。注意分类体系一旦确定九张图的相关图例、配色、统计口径都要统一。最忌讳的是土地利用图用一套分类、生态系统图用另一套分类、报告正文又用第三套口径三处对不上评审基本过不了。2. 数据源准备与核心参数解析2.1 遥感数据、地形数据与专题数据的选型逻辑九张图的数据源可以归为三类。第一类是遥感影像数据主要用于土地利用现状解译、植被类型判读、NDVI计算和植被覆盖度反演。常用的中分辨率数据包括陆地卫星系列和国产高分系列空间分辨率从10米到30米不等。选型时看两个指标空间分辨率和时相。生态评价通常要求使用近期影像时相上优先选植被生长旺盛季这样植被类型判读和NDVI计算才有代表性。如果项目区跨气候带或地形复杂单一时相可能不够需要多时相辅助判读。第二类是地形数据主要是数字高程模型用于水系提取、坡度计算和土壤侵蚀模型。常用的有30米和12.5米分辨率的公开DEM数据精度要求高的项目可以用更高分辨率的测绘成果。DEM的精度直接影响水系提取的准确性和坡度因子的计算进而影响土壤侵蚀图的可靠性。第三类是专题数据包括项目区边界、行政区划、水系矢量、土壤类型图、气象数据等。项目区边界必须用正式的勘界定标成果这是位置图的基准。土壤类型数据用于土壤侵蚀模型中的可蚀性因子气象数据用于降雨侵蚀力因子。这些数据往往来自不同部门格式、坐标系、精度各不相同整合工作量不小。数据类别典型来源主要用途关键精度要求遥感影像陆地卫星、国产高分系列土地利用解译、NDVI、植被覆盖度分辨率与成图比例尺匹配时相为生长季数字高程模型公开DEM或测绘成果水系提取、坡度、土壤侵蚀分辨率不低于30米高精度项目用12.5米项目区边界勘界定标成果位置图基准、范围裁剪必须为正式成果坐标系统一土壤类型数据土壤普查成果土壤可蚀性因子图斑边界与制图范围匹配气象数据气象站点或格网数据降雨侵蚀力因子系列长度满足模型要求2.2 NDVI计算与植被覆盖度反演的参数细节归一化植被指数是近红外波段与红光波段反射率的差值比上和公式本身不复杂但实操中有几个细节决定成败。第一大气校正做不做。原始影像的反射率受大气影响直接算NDVI会有偏差尤其是跨时相比较多时不做大气校正的NDVI不可比。第二云和阴影的掩膜。生长季影像云量往往不低云和云影必须剔除否则NDVI图上会出现大片异常低值。第三水体、建筑等非植被地物的处理这些地物NDVI值本身偏低统计时要单独说明不能混进植被覆盖度的均值里。植被覆盖度反演常用像元二分模型核心是确定纯裸地和纯植被的NDVI阈值。这两个阈值的取法直接影响覆盖度结果常见做法是取NDVI累积频率的5%和95%分位数作为裸地和纯植被端元也有用置信区间截断的做法。我的经验是阈值不能照搬别的项目必须结合本项目区的NDVI直方图来定干旱区和湿润区的端元值差别很大。定完阈值后覆盖度按NDVI减裸地端元值除以植被端元值减裸地端元值计算结果截断在0到1之间。提示NDVI端元阈值确定后建议在报告里写明取值方法和具体数值并附上NDVI直方图。这是评审常问的点写清楚了能省掉现场解释。2.3 土壤侵蚀模型与生物量估算的因子准备土壤侵蚀图通常用通用土壤流失方程或其修正版本估算核心是降雨侵蚀力、土壤可蚀性、坡长坡度、植被覆盖与管理、水土保持措施五个因子相乘。降雨侵蚀力因子由气象数据计算土壤可蚀性因子由土壤理化性质查算坡长坡度因子由DEM计算植被覆盖与管理因子和覆盖度、土地利用类型相关水土保持措施因子按措施类型赋值。五个因子逐个算完再相乘得到侵蚀模数然后按分级标准分强度等级。生物量估算有几种路径。有实测样方数据的可以用样方生物量结合植被类型面积外推没有实测数据的用遥感反演的植被指数结合经验模型估算比如用NDVI或增强植被指数与生物量的统计关系。无论哪种路径都要说明估算方法和不确定性。生物量分布图的价值在于反映生态系统的碳储量和生产力水平是生态现状评价的重要支撑。3. 九张图件的实操制作流程3.1 项目区位置图一张图讲清在哪、周边有什么位置图是九张图的门面制作逻辑是由大到小、层层定位。通常用三级嵌套第一级是区域位置展示项目在较大地理范围内的方位第二级是项目区周边展示项目与周边城镇、交通干线、自然保护地等敏感目标的关系第三级是项目区范围用醒目边界标出。三级嵌套可以用一张图内嵌小图的方式也可以分幅表达。实操上先加载行政区划底图和项目区边界设置好投影和范围然后叠加地形晕渲或影像底图增强可读性。项目区边界用粗线加填充透明色突出周边敏感目标用点符号加标注。比例尺、指北针、图例、经纬网四要素缺一不可。我习惯在位置图上标注项目区到最近敏感目标的直线距离这个信息专家很关注能直观反映区位敏感性。3.2 土地利用现状图与生态系统类型图解译、分类与映射这两张图是本底组的核心制作流程是影像解译、图斑勾绘、分类赋值、拓扑检查、符号化出图。解译方式有目视解译和自动分类两种生态评价图件对精度要求高通常以目视解译为主、自动分类为辅。解译时先建立解译标志比如不同地类在影像上的色调、纹理、形状特征形成解译标志表后再大面积勾绘。图斑勾绘完成后必须做拓扑检查查重叠、查缝隙、查悬挂点。这一步最枯燥但最不能省图斑有缝隙或重叠面积统计必然出错。分类赋值按前面确定的分类体系来赋值完成后做面积汇总和报告正文的统计表核对。生态系统类型图在土地利用图基础上做映射映射关系用对照表固定下来保证两张图面积一致。符号化时土地利用图用标准配色耕地、林地、草地、水域、建设用地、未利用地各用约定俗成的色系生态系统类型图用另一套配色但保持逻辑对应。图例要完整列出所有出现的类型不能只列主要的。3.3 地表水系图从DEM提取到矢量核对的完整链路水系图的制作有两条路径一是用现有水系矢量数据二是从DEM自动提取。实操中往往是两者结合自动提取的结果用来核对和补充现有数据。从DEM提取水系的流程是填洼、流向计算、流量累积、河网提取、河网分级关键参数是最小汇水面积阈值阈值设小了河网过密、设大了支流丢失需要结合项目区实际水系密度反复调试。提取出的河网要和影像、现有水系图叠加核对修正明显错误。水系图要区分常年河、时令河、湖泊、水库、坑塘等类型用不同线型和填充表达。项目区内的水系要标注名称项目区外的干流支流也要适当表达体现水系连通关系。如果项目涉及涉水工程还要标出工程与水域的位置关系。注意水系提取的最小汇水面积阈值没有万能值平原区和山区的合理阈值差别很大。我的做法是先按经验值跑一版和影像比对后调整通常要试三到五版才能定下来。3.4 植被类型图与植被覆盖度图判读边界与分级表达植被类型图在土地利用图的基础上细化重点是把林地、草地细分到植被型或群系组。判读时结合影像纹理、色调和地形部位有条件的辅以实地调查样方。植被类型边界要和土地利用图的林地草地边界保持一致不能出现土地利用图是林地、植被类型图却是灌丛的矛盾。植被覆盖度图是分级图通常按覆盖度高低分几级比如低覆盖、中低覆盖、中覆盖、中高覆盖、高覆盖。分级阈值结合项目区实际和评价需要确定分级后用渐变色系表达从浅绿到深绿。分级图的关键是分级标准要写清楚图例要标注每级的覆盖度区间。3.5 NDVI图与生物量分布图定量数据的可视化处理NDVI图是连续值栅格可视化时用渐变色带表达从低值到高值对应不同颜色。出图前要做拉伸通常用2%到98%的累积频率截断避免极值把色带拉偏。图例要标注NDVI的取值范围和色带对应关系。如果NDVI图用于对比分析色带分级要统一不能一张图一个色带。生物量分布图同样是连续值或分级值制作时先完成生物量估算得到栅格或图斑数据再分级符号化。分级可以用自然断点法或等间距法前者更能反映数据分布特征。图例标注生物量单位和分级区间报告正文说明估算方法和精度。3.6 土壤侵蚀图因子计算、分级与制图表达土壤侵蚀图的制作链路最长五个因子逐个计算再综合。降雨侵蚀力因子用气象站点数据计算后插值成面土壤可蚀性因子用土壤类型图关联理化性质查算坡长坡度因子用DEM计算植被覆盖与管理因子结合土地利用和覆盖度赋值水土保持措施因子按措施类型赋值。五因子相乘得到侵蚀模数栅格再按分级标准分微度、轻度、中度、强烈、极强烈、剧烈等级。制图时按侵蚀强度分级符号化用从浅到深的色系表达图例标注每级对应的侵蚀模数区间。项目区内的侵蚀斑块要清晰可辨和土地利用图、植被覆盖度图叠加分析能看出侵蚀高值区往往对应植被覆盖度低值区这种图间一致性是评审看重的。4. 常见问题排查与实操避坑经验4.1 图件之间打架的典型场景与解决思路九张图出自同一套数据体系但实操中最常见的问题就是图间不一致。典型场景有这么几个土地利用图面积和生态系统类型图面积对不上多半是映射关系没固定或图斑边界有差异植被类型图的林地边界和土地利用图的林地边界错位多半是两次勾绘没有共用底图NDVI高值区和植被覆盖度高值区对不上多半是端元阈值取值不当或云掩膜不彻底土壤侵蚀高值区和植被覆盖度低值区不吻合多半是因子赋值有误。解决思路是建立一致性检查清单出图前逐项核对面积汇总是否一致、边界是否共用、分级是否对应、统计口径是否统一。我习惯把九张图叠在一个工程里用卷帘工具来回对比边界错位一眼就能看出来。问题现象可能原因排查方法解决措施土地利用与生态系统面积不符映射关系未固定、图斑边界差异面积汇总对比、边界叠加检查固定映射对照表、共用底图勾绘植被类型与土地利用林地边界错位两次勾绘未共用底图卷帘对比、边界叠加统一底图、边界继承NDVI与覆盖度高值区不吻合端元阈值不当、云掩膜不彻底检查直方图、检查掩膜重定端元、完善掩膜侵蚀高值区与覆盖度低值区不吻合因子赋值错误逐因子核查修正因子赋值4.2 出图规范与评审关注的细节清单图件规范是评审的硬门槛。必备要素包括图名、图例、比例尺、指北针、经纬网或公里网、制图单位、制图日期。图名字体、字号要统一图例位置要一致九张图放在一起要像一套。坐标网标注要规范投影信息要在图廓外注明。评审常关注的细节还有项目区边界是否醒目、敏感目标是否标注、分类是否完整、图例是否与图面一致、面积统计是否与正文一致。我见过图例漏列一个地类、图面却有该地类图斑的情况这种低级错误最伤专业形象。出图前逐张过一遍要素清单能避免大部分问题。提示九张图建议用同一套图框模板图名、图例、比例尺、指北针的位置全部固定只换图面内容。这样出图效率高成套性也好。4.3 数据精度与成图比例尺的匹配经验数据精度和成图比例尺不匹配是隐性问题评审不一定当场指出但会影响图件可信度。30米分辨率影像解译的土地利用图成图比例尺不宜大于1:50000硬出1:10000的图图斑边界锯齿明显细碎地类无法表达。DEM提取的水系30米分辨率适合中小比例尺大比例尺水系图需要更高精度数据或实测矢量。我的经验是先确定成图比例尺再倒推需要的数据精度。位置图可以用1:100000到1:50000专题图根据评价等级和项目规模定通常1:50000到1:10000。比例尺定了数据源精度、图斑最小上图面积、线状要素最小长度都跟着定整套图件的精度体系就统一了。4.4 批量出图与模板化提效的实操心得九张图如果每张都从头设置工作量翻倍还容易出错。我的做法是建工程模板预设坐标系、制图范围、图框、图例框、比例尺、指北针、图名字体九张图共用。数据准备好后逐张替换图面内容导出时用批量导出功能统一分辨率和格式。符号库也要提前建好土地利用、生态系统、植被类型、土壤侵蚀各建一套配色方案存成样式文件下次项目直接调用。NDVI和覆盖度的色带也存成模板保证同类图件风格一致。这套模板化做法熟练之后九张图从数据准备到出图工作量能压缩不少而且成套性和规范性明显提升。5. 从图件到结论空间证据链的闭环九张图做完不是终点图件要能支撑报告结论。位置图支撑区位敏感性分析土地利用和生态系统类型图支撑生态现状本底评价水系图支撑水生态影响分析植被三件套支撑植被质量评价土壤侵蚀图支撑水土流失影响分析生物量图支撑生态系统生产力评价。每张图在报告正文里都要有对应的文字解读图说数据、数据说结论形成闭环。我在实际项目中的体会是图件质量的分水岭不在软件操作而在制图前的数据体系设计和制图后的一致性核对。软件操作学几天就会但数据关系想不清楚、图间矛盾不核对图做得再花哨也经不起评审推敲。另外一个小技巧九张图全部完成后把缩略图拼成一张总览图自己看一遍哪张图风格突兀、哪张图信息密度异常一眼就能发现比逐张检查效率高得多。