基于无人机数据的滨海湿地互花米草%28Spartina alterniflora%29地上生物量反演研究.pdf

生态与农村环境学报，（）：收稿日期：基金项目：国家自然科学基金（，）；江苏省科技计划项目社会发展重大科技示范项目（）；江苏省自然资源发展专项基金（）通信作者：苏孟园，吴涵，何柯欣，等基于无人机数据的滨海湿地互花米草（）地上生物量反演研究生态与农村环境学报，（）：，（）：基于无人机数据的滨海湿地互花米草（）地上生物量反演研究苏孟园，吴 涵，何柯欣，李昕阳，李玉凤（南京师范大学海洋科学与工程学院，江苏 南京）摘要：互花米草（）快速入侵的关键在于其生长繁殖以及对环境的适应能力，互花米草生物量的大小正是它入侵能力的重要体现。以无人机遥感影像为数据源，通过建立关于植被覆盖度和植株高度的模型对互花米草地上生物量进行估算研究。研究结果表明，在条子泥研究区内互花米草地上生物量在 之间，地上生物量平均值为 ，的互花米草地上生物量小于 ，地上生物量为 以及 的互花米草分别占 和 ，地上生物量较高的互花米草主要分布在靠近海堤区域与近海一侧潮滩，中间区域地上生物量较低。遥感影像景观分类结果 系数为 ，互花米草地上生物量估算的均方根误差为 ，反演结果能够在一定程度上反映研究区互花米草地上生物量分布情况。互花米草分布多呈斑块状，部分呈丛分布，处于互花米草快速扩张阶段。研究所采用的无人机遥感能够同步获取采样时间点研究区的影像数据，可以很好地应用于时效性较强的生物量遥感反演。该研究为互花米草地上生物量的定量化反演提供了新的方法。关键词：无人机；地上生物量；植被覆盖度；植株高度；互花米草中图分类号：文献标志码：文章编号：（），（，）：，（），；，第 期 苏孟园等：基于无人机数据的滨海湿地互花米草（）地上生物量反演研究 ：（）；作为典型的盐生植物，互花米草（）原产于北美大西洋沿岸，因具有保滩护岸、加快淤积的作用，于 年 月引入中国，在我国沿海迅速繁衍扩张、竞争，取代土著植物，占领光滩，影响入侵区地表沉积速率与水动力过程，对原生海滨湿地生态系统及当地居民经济生产活动产生巨大影响，成为近年来我国有关生物入侵问题争论的焦点。互花米草快速入侵的关键在于其生长繁殖能力以及对环境的适应能力，互花米草生物量的大小正是它入侵能力的重要体现。同时，生物量反演对于互花米草的科学治理及其生物固碳研究具有重要意义。然而江苏淤泥质滨海湿地的难以进入性限制了对互花米草的调查和研究。目前遥感技术已经被广泛应用于植被监测等研究领域，而受限于滨海区域复杂的气候条件，难以获取目标时相高质量的卫星遥感影像，且野外采样与卫星遥感数据的获取不能完全同步，导致数据获取存在时间差，信息无法完全匹配。相比于卫星遥感，无人机遥感的飞行高度较低，空间分辨率高，能够实时获取研究区采样时间点的影像数据，可对难以进入的研究区进行拍摄，在时效性较强的生物量遥感反演方面有很好的应用前景。在目前的研究中，基于光学遥感数据的生物量反演多是根据植物光谱特性，形成植被指数，通过分析植被指数与生物量之间的统计特征关系建立数学估算模型，而湿地植被生物量的决定因素不仅包括叶片的平均密度，植被茎秆也是生物量的重要组成部分，因此仅根据植被光谱特征和纹理特征的生物量估算存在难以避免的误差。目前多是根据单一的植株高度等形态指数或植被指数建立反演模型，对互花米草地上生物量进行反演，而以一个综合的模型结合植被覆盖度与植株高度这两方面信息对互花米草地上生物量进行反演的研究较少。研究基于无人机遥感获得的研究区可见光影像数据，分别提取互花米草高度信息及覆盖度信息，根据植被地上生物量计算公式对研究区互花米草地上生物量进行反演，为互花米草地上生物量的定量化反演研究提供了新的方法。材料与方法 研究区概况研究区位于江苏省盐城市东台区的条子泥湿地（，），该区域地处我国典型季风气候区，受海洋性和大陆性气候双重影响，年平均气温 ，年平均降水量 ，条子泥湿地是原始和典型的淤泥质滩涂，生长的植被保持自然演替状态，生物多样性丰富，拥有世界上面积最大的潮间带湿地，是全球最重要的滨海湿地生态系统之一。研究区受互花米草入侵，互花米草广泛分布于潮滩，属于区域的优势物种（图）。图 研究区位置 生 态 与 农 村 环 境 学 报第 卷 数据来源及处理 无人机遥感影像获取及处理研究采用的是大疆精灵 无人机，无人机配备 英寸 万像素 影像传感器，可获取红、绿、蓝真彩色 格式影像。于 年 月 日：对研究区进行拍摄，在设定航高下影像空间分辨率为 。月研究区大部分互花米草虽然已经干枯，但是互花米草植株茎秆比较坚挺，且冬季处于低潮期，研究区受潮汐影响有限，互花米草地上生物量的损失比较少，对研究的影响不大。互花米草生物量实测数据于 年 月 日：在研究区所选样点处进行同步采样，布设 个大 小 为 的样方，齐地收割各样方内全部互花米草植株。现场测得互花米草鲜重的平均值、最大、最小值分别为 、，各样方互花米草的平均株高为 ，最大株高为 ，最小株高为 ，植被覆盖度的平均值、最大、最小值分别为 、。研究方法基于无人机遥感影像数据对研究区互花米草地上生物量进行反演的流程主要包括 部分：无人机遥感影像景观分类、互花米草高度提取、互花米草覆盖度计算及计算植被地上生物量。无人机遥感影像景观分类无人机影像的高空间分辨率在加强地物细节的同时也增大了同种地物之间的差异以及不同地物间的相似程度，导致分类精度降低，而面向对象的分类方法与传统的分类方法相比，不再仅仅依靠地物的光谱特征，而是更多地利用地物的几何信息与结构信息，将分类对象分割成为内部具有相同或相似特征的基本单元，后续的分析与处理也都基于基本单元进行，这样的分类方法能够有效平滑影像噪音和不必要的细节，提高分类结果的精度。使用 对研究区无人机影像进行面向对象的分类，选择多尺度分割算法，经过反复小区域试验结果，分割尺度设置为，形状因 子 参 数 设 置 为 ，紧 凑 度 参 数 设置为，其余参数为默认值，得到的对象大小适中、内部光谱变异小，与其他地物的边界清晰准确，分割效果较好。分割得到的结果采用最邻近特征算法，手动选择样本对研究区影像进行分类，导出分类结果。另外，在野外调查的基础上，通过目视解译的方法对错分漏分的地物进行修正，以确保分类精度。高程与淹水深度等变化对互花米草高度与生物量有重要影响。由于区域较为平坦，高程变化不大，总体表现为由海堤向海呈带状缓慢降低的趋势，因此按照距离海堤远近的等距离划分法，由陆向海将互花米草集中分布的区域划分为距离堤坝：、及 共 个带状区域。区域已处于研究区互花米草分布最外围，互花米草分布稀疏零散，因此将此区域划分至同一条带，其他分区面积大致相同（图）。图 研究区景观分类结果 互花米草高度提取研究使用 软件导入无人机遥感影像，添加控制点对无人机影像进行空中三角测量，后进行空三检查、校正错误，提交重建任务，选择三维网络进行原始 建模，在此基础上选择正射影像为建模目的，进行数字地表模型（，）类型的图像建模，即可得到研究区。在面向对象分类结果的基础上，使用 工具提取出 分类结果中的互花米草部分，利用缓冲区工具设置 缓冲区，根据缓冲区对研究区 模型进行提取。借助分区统计工具计算缓冲区高度值作为互花米草生长区域地面高度，提取研究区 模型中互花米草所在区域高度，两者相减即可得到互花米 第 期 苏孟园等：基于无人机数据的滨海湿地互花米草（）地上生物量反演研究草高度。互花米草覆盖度计算无人机携带的 影像传感器所拍摄的影像仅包含可见光波段。研究区位于滨海湿地，背景土壤中吸附水分较多，会造成蓝光波段反射率较高。通常情况下使用的基于无人机可见光波段数据的植被指数，如可见光差异植被指数（，）、过绿指数（，）等在研究区的适用性较差，不能很好地区分不同植被覆盖度的互花米草。该研究在分辨效果较好的归一化绿蓝指数（，）基础上进行改进，构建归一化红蓝指数（，）来反映植被信息，计算公式为 。（）式（）中，为 指数值；为红光波段数值；为蓝光波段数值。在此基础上，使用像元二分法模型计算研究区互花米草植被覆盖度。（）式（）中，为植被覆盖度；为像元 指数值；为裸地或无植被覆盖像元 指数值；为全覆盖纯植被像元 指数值。研究区 指数计算结果的范围为 ，由于无人机影像数据计算过程可能存在图像误差，且影像中人造堤坝等地物可能会对指数计算结果产生影响，不能直接将计算所得到的 指数的最大值与最小值代入为 和，因此统计各指数值的像元占比，对 指数置信度区间进行分析，取累计频率为 的 指数值为，取值为 ，取累计频率为 的 指数值为，取值为 。互花米草地上生物量计算地上生物量指一定时间内地表单位面积存在的有机质总量，以覆盖度来表示单位面积上的植被面积，则互花米草地上生物量（，）可以用下式进行计算：。（）式（）中，为互花米草名义生物量，即覆盖度为 时单位体积互花米草的鲜重，；为互花米草植被覆盖度；为互花米草平均株高，。通过计算互花米草植被覆盖度及互花米草的高度即可估算出互花米草地上生物量。无人机反演结果精度验证方法采用基于混淆矩阵的 系数来对无人机遥感影像景观分类结果的精度进行检验。在研究区随机选取 个样点，根据实测数据与高空间分辨率的遥感影像，通过人工目视解译的方法确定所选取样点的景观类型作为地物的实际景观类别，遥感分类结果作为预测景观类别，构建混淆矩阵计算景观分类结果的 系数。利用照相法和专家判读相结合的方法确定 个实测点位的植被覆盖度，将实测植被覆盖度与利用像元二分法模型计算得到的植被覆盖度进行比较，通过计算相对误差与平均误差对研究区互花米草植被覆盖度的反演结果精度进行检验。，（）。（）式（）（）中，为相对误差；为遥感估算的植被覆盖度；为实测植被覆盖度；为平均误差；为实测样本总数。根据 个实测点位的测量数据对研究区互花米草植株高度和地上生物量反演结果的精度进行分析，以均方根误差（，）来定量化表示实测值与反演结果之间的关系。（）。（）式（）中，为实测数值；为反演数值；为实测样本总数。结果与分析 互花米草分布特征 景观分类结果研究区互花米草斑块及面积占比如图 所示，互花米草在靠近堤坝区域集中连片呈条带状分布，中部区域为互花米草与碱蓬交错带，互花米草零散、破碎化分布，到近海区域又有呈条带状分布的互花米草，向海有互花米草斑块零散分布，呈现出向外扩张的趋势。研究区互花米草总面积为 ，互花米草总斑块数量约 个，平均斑块面积为 ，最大斑块面积为 ，最小斑块面积为 。生 态 与 农 村 环 境 学 报第 卷图 研究区互花米草斑块及面积占比 互花米草空间分布特征研究区互花米草最大斑块位于距堤坝 以内区域。各个划分区域内互花米草面积与平均斑块面积由陆向海总体上逐渐减小，仅在 范围内有所增大随后又逐渐减小。由图 可知，由陆向海互花米草斑块数量先增加后减少，外的斑块数量最少，共有 个；范围内斑块数量最多，有 个；在各区域内分布的互花米草占研究区互花米草总面积的比例中，范围内互花米草面积占比最大，为 ，其次是 范围，占比为 ，以外占比最小，为 。表明研究区互花米草空间分布上由陆向海经历了集中连片到零散破碎再到集中连片的过程。互花米草高度和覆盖度研究区互花米草高度与覆盖度计算结果如图 所示。互花米草高度范围为 ，平均高度为 ，靠岸与近海区域互花米草较高，中间区域互花米草较矮。区域互花米草的平均高度最大，为 ；其次是 区域，为 ，再次是 以及 区域，均为 ；范围内最小，为 。高度在 的互花米草所占面积最大，为 ，占研究区互花米草总面积的 ，和 的互花米草分别占 和 ，的互花米草仅占 （表）。研究区互花米草高度呈现出由陆向海逐渐降低后又逐渐升高的变化趋势。图 研究区互花米草高度与植被覆盖度 研究区互花米草覆盖度平均值为 ，覆盖度为 的互花米草占比最大，占研究区互花米草总面积的 ，其次是 和 范围，分别占 和 ，覆盖度的互花米草所占面积最小，为 。覆盖度平均值最大的区域是 范围，为 ，第 期 苏孟园等：基于无人机数据的滨海湿地互花米草（）地上生物量反演研究范围为 ，与 范围均为 ；范围最小，为 。研究区互花米草生长区域覆盖度总体较高，由陆向海呈现出先减小后增大的变化趋势。表 研究区互花米草植株高度 植株高度 像元数面积 占比 互花米草地上生物量反演结果 基于无人机遥感反演结果精度检验遥感影像景观分类结果经过计算 系数为，说明研究区景观分类结果与实际地物类型已经高度一致，符合做进一步研究的精度要求。研究区互花米草植被覆盖度反演结果的平均误差为，在误差允许的范围内。以研究区实测互花米草植株高度对高度提取结果的精度进行检验，结果显示互花米草植株高度提取结果的均方根误差为 ，其误差产生的原因主要是互花米草自然状态下受到环境扰动而出现倒伏以及无人机影像分辨率的限制。研究区互花米草地上生物量估算结果的均方根误差为 ，误差来源主要为引入的互花米草植株高度与植被覆盖度所造成的累积误差，反演结果对于反映研究区互花米草地上生物量具有一定的参考价值。互花米草地上生物量空间分布特征根据实测数据中互花米草地上生物量数据，获得研究区互花米草名义生物量（）为 ，再根据上述分析所得到的互花米草高度与植被覆盖度信息进行波段计算，即可得到研究区互花米草地上生物量。研究区互花米草地上生物量在 之间，平均值为 ，的互花米草地上生物量 ，地上生物量在 之间的互花米草占 ，地上生物量 的互花米草占 。在划定的各个区域中，互花米草地上生物量的平均值最高，为 ；其次是 范围，平均值为 ，再次是 范围，平均值为 ，范围互花米草地上生物量的平均值最低，为 （图）。在靠近堤坝的带状区域以及近海一侧零散分布的斑块互花米草地上生物量普遍较高，位于中部的互花米草地上生物量则相对较低。图 研究区互花米草地上生物量反演结果 不同剖面互花米草地上生物量变化特征在根据植株高度与植被覆盖度得到的互花米草地上生物量反演结果基础上，垂直于岸线设置 条剖面线（图），使其相对均匀地分布于研究区范围，对互花米草由陆向海地上生物量变化情况进行分析。图 为各剖面互花米草地上生物量变化图。在垂直于岸线方向上，由陆向海 范围内互花米草地上生物量较高，为 左右；随着离堤距离的增大，互花米草地上生物量快速减小；在 范围内互花米草地上生物量主要在 左右，在距离堤坝较远一侧区域互花米草地上生物量在 之间，部分区域互花米草地上生物量在 以上，互花米草地上生物量分布沿剖面线变化趋势为高低高。根据图 可以看出，在平行于岸线方向上，研究区北侧 剖面区域近堤一侧互花米草地上生物量大于远堤一侧。逐渐向南的过程中，近堤一侧互花米草地上生物量逐渐减小，而远堤一侧互花米草地上生物量逐渐增大。研究区互花米草地上生物量最高的区域集中于西北、东南方向，西南、东北方向互花米草地上生物量较高，中部区域互花米草地上生物量普遍较低。生 态 与 农 村 环 境 学 报第 卷图 不同剖面互花米草地上生物量变化 讨论 互花米草地上生物量空间分异原因有研究表明，定居时间的长短会影响互花米草植株高度。一般情况下，多年定居互花米草群落的植株较高。对于互花米草的有性繁殖，幼苗定居成功率与裸露斑块大小呈正相关，幼苗很难在成熟植被下定居。互花米草生长繁殖受所在区域土壤水盐条件影响较大，在一定的范围内，互花米草株高和克隆存活率会随着土壤含水量的提高而升高，而且潮汐淹没频率的提高也会刺激互花米草无性分株的生成。研究区近堤一侧互花米草定居时间较长，植株较高，地上生物量较高。位于中部区域的互花米草定居时间比靠近堤坝区域短，在生存环境方面比近海一侧土壤含水量与潮汐淹没频率低，且该区域碱蓬的分布也会给互花米草带来一定的竞争压力，互花米草地上生物量较低。互花米草在向外扩张的过程中借助繁殖体进行短程“流”式传播的同时利用种子进行近距离跳跃式扩散，远堤一侧互花米草密度较低，裸露滩涂较多且无碱蓬竞争，距离堤坝越远潮汐淹没频率越高，土壤含水量也相应增大，有利于互花米草生长，地上生物量较高。研究区互花米草分布形态不规则，多数呈斑块状分布，部分呈丛分布，处于互花米草快速扩张阶段。沿海潮滩生境极其严酷，互花米草入侵初期作为群落演替的先锋植物占据空生态位，江苏沿海 标准海岸线中 以上属于淤泥质海岸，在互花米草入侵前多为裸露的泥滩，且潮滩滩面广阔、坡度平缓，为互花米草的扩张提供了广阔的空间。研究区互花米草通过克隆繁殖，沿平行于海岸线方向进行短程“流”式扩张，同时利用潮汐将种子和幼苗运送至潮滩前沿，生成零散的新生斑块，进行垂直于海岸线的扩张，这与其他研究的互花米草在侵入地的扩张策略相符合。互花米草地上生物量反演结果反演得到的研究区互花米草地上生物量在 之间，平均值为 ，的互花米草地上生物量 ，地上生物量在 之间的互花米草占 ，的占 。周在明等通过遥感影像估算三沙湾滩涂互花米草地上生物量，发现 的互花米草地上生物量在 之间，总体上高于笔者研究。前者研究影像的获取时间为 年 月，处于一年中互花米草生物量稳定较高的时间，研究区处于半封闭内湾，泥沙条件好且气候属于中亚热带季风湿润气候，更加适宜互花米草生长；而笔者研究影像获取时间为 月，大部分互花米草老植株处于生长末期已经干枯而新植株尚未萌发，导致地上生物量普遍较低。虞海英 第 期 苏孟园等：基于无人机数据的滨海湿地互花米草（）地上生物量反演研究在使用机载雷达和高光谱数据对大丰市沿海滩涂湿地湿生植被地上生物量进行反演的研究中发现，研究区互花米草地上生物量为 ，分布在 的互花米草占比极高，研究区影像获取时间为 年 月，在气候条件相近、数据获取时间相差 个月左右的条件下，互花米草地上生物量反演结果最大值较大，总体上笔者研究与其相近。韩爽等在对盐城自然保护区核心区植被地上生物量进行的遥感反演研究中发现，互花米草地上生物量湿重主要集中在 以上，干重集中在 之间，研究影像获取时间为 年 月，此时气温升高，正是互花米草快速生长时期，地上生物量与该研究结果相比较高。结论基于无人机遥感数据提取互花米草植株高度与覆盖度信息，通过植被地上生物量计算公式对研究区互花米草地上生物量进行反演，结论如下：（）研究采用的反演方法所得到的结果能够在一定程度上体现研究区互花米草地上生物量分布情况。与实测数据相比，互花米草地上生物量估算的均方根误差为 。（）研究区互花米草地上生物量在 之间，平均值为 ，靠近堤坝与近海一侧地上生物量较高，中部区域较低，是互花米草定居时长和分布区域的环境条件不同导致的。（）研究区互花米草分布形态不规则，多数呈斑块状分布，部分呈丛分布，处于互花米草快速扩张阶段。研究引入植株高度与植被覆盖度对互花米草地上生物量进行估算。在今后的研究中将进一步结合遥感影像数据以及研究区环境因素构建更多因素的反演模型，以提高互花米草地上生物量反演精度，同时注意长时间序列的遥感影像与实测数据以及研究区潮汐淹没频率、土壤含水量、盐度等数据的收集，以便对研究区互花米草分布及原理进行进一步分析。参考文献：，（）：毛立，孙志高，陈冰冰，等闽江河口互花米草入侵湿地土壤无机硫赋存形态及其影响因素生态学报，（）：，（）：潘卫华，陈家金，王岩近 年福建红树林和互花米草群落时空变化及景观特征生态与农村环境学报，（）：，（）：陈思明互花米草（）在中国沿海的潜在分布及其对气候变化的响应生态与农村环境学报，（）：，（）：时宇，黄黄，冉珊珊，等氮磷添加对互花米草土壤系统碳分配的影响生态与农村环境学报，（）：，（）：陈怀璞，张天雨，葛振鸣，等崇明东滩盐沼湿地土壤碳氮储量分布特征生态与农村环境学报，（）：，（）：周在明，杨燕明，陈本清基于无人机影像的滩涂入侵种互花米草植被信息提取与覆盖度研究遥感技术与应用，（）：，（）：，（）：崔万新，李锦荣，司前程，等基于无人机可见光数据荒漠灌木覆盖度提取方法研究水土保持研究，（）：，（）：，刘琳，郑兴明，姜涛，等无人机遥感植被覆盖度提取方法研究综述东北师大学报（自然科学版），（）：，（），（）：赵天舸，于瑞宏，张志磊，等湿地植被地上生物量遥感估算方法研究进展生态学杂志，（）：生 态 与 农 村 环 境 学 报第 卷，：，（）：周在明，杨燕明，陈本清三沙湾滩涂互花米草地上生物量估算湿地科学，（）：，（）：虞海英融合机载 和高光谱数据的滨海湿生植被生物量反演方法研究北京：中国测绘科学研究院，：，韩爽，甄艳，谭清梅，等盐城海滨湿地植被地上生物量遥感反演研究林业资源管理，（）：，（）：刘潞，余夏杨，唐洪根，等围垦对条子泥互花米草种群年季扩张特征的影响农业资源与环境学报，（）：，（）：刘舒，朱航基于超高空间分辩率无人机影像的面向对象土地利用分类方法农业工程学报，（）：，（）：王启田，林祥国，王志军，等利用面向对象分类方法提取冬小麦种植面积的研究 测绘科学，（）：，（）：李营，陈云浩，陈辉，等 影像的翅碱蓬植被指数构建武汉大学学报（信息科学版），（）：，（）：李丽纯，林瑞坤，吴振海，等闽江口湿地互花米草地上生物量的月动态特征资源科学，（）：，（）：周在明，杨燕明，陈本清基于 遥感影像的滩涂湿地入侵种互花米草植株高度的反演研究海洋学报，（）：，（）：任丽娟，王国祥，仇乐，等江苏潮滩湿地不同生境互花米草形态与生物量分配特征生态与农村环境学报，（）：，（）：，（），（）：马喜君，陆兆华，刘黎华江苏省盐城市海滨湿地互花米草入侵的生态风险评价江苏农业科学，（）：，（）：，（）：王娟，刘红玉，李玉凤等入侵种互花米草空间扩张模式识别与景观变化模拟生态学报，（）：，（）：程成，雷金睿，何松等海口市互花米草分布现状和扩散动态湿地科学，（）：，（）：作者简介：苏孟园（），女，安徽临泉人，研究方向为海洋资源开发技术。：（责任编辑：陈 昕）