本实验采用2001年23期时序MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid(简称MODIS13Q1)数据, 每期间隔16天, 属于陆地专题产品, 是MODIS的L3级科学数据集产品之一。MODIS13Q1产品原始文件包含多个数据层, 实验选取第一层归一化植被指数(Normalized Difference Vegetation Index, 简称NDVI)数据文件^[13]。

本实验采用2001年1∶100万中国植被图矢量数据, 可提供实验区植被分布信息; 中国植被区划矢量数据提供实验区范围; 2000年中国土地覆盖1km栅格数据^[14-15](简称WESTDC), 可提供实验区非植被分布信息。其中, WESTDC数据来源于国家自然科学基金委员会“中国西部环境与生态科学数据中心”(http://westdc.westgis.ac.cn)), 该数据借鉴IGBP分类系统, 在保持了中国土地利用数据的总体精度的同时, 补充了中国植被图中对植被类型及植被季相的信息, 更新了中国湿地图, 增加了中国冰川图, 融合了MODIS 2001年土地覆盖产品。

以上矢量、栅格数据均具有空间坐标信息, 且经过影像配准, 精度限制在一个像元内。

本研究选取寒温带针叶林区域作为实验区, 实验区处于127°20′E以西, 49°20 ′N以北的大兴安岭北部及其支脉伊勒呼里山地, 总面积约为211 600 km²。植被覆盖率高, 地带性植被为兴安落叶松林, 有少量灌木植被, 是我国重要的木材产地。由于该区域的寒温带气候, 形成永久冻土层, 水分下渗困难, 形成零星水体, 地表积水使湿地广泛发育, 其中草本沼泽分布在海拔900 m以下地带。该区域平坦谷地区域有农业分布但农业不发达。

首先, 通过ARCGIS软件在寒温带针叶林区域内随机布点, 生成的样点具有中国植被区划文件的空间信息, 每个样点在实验影像上体现为一个像素, 代表实际地物250 m×250 m的样地范围。参考1∶100万中国植被图和中国土地覆盖数据, 分别赋予植被类样点的类型属性, 包括: 乔木、草地、湿地、灌木、农田; 非植被类样点的类型属性, 包括: 建筑物和水体。并对各类样本的初始可分离性进行计算。

然后, 对实验区2001年23期MODIS时序NDVI影像数据产品, 依据不同植被类型的NDVI特征向量不同, 采用非监督分类迭代算法^[15]得到时序NDVI数据波段的集群聚类信息。结合该非监督分类结果, 依据矢、栅数据的空间特征, 将非监督分类的类型信息作为各类样本中的样点属性进行关联, 统计各类样本包含的非监督分类类型和样点个数(j, N_Ri), 并将对应类型的样点比例作为权重值(P_Ri)。权重越大则越能够代表该样本类型, 权重越小则越可能是该类样本中的非纯样点, 按照权重大小降序排列样点(Descend(P_Ri))并计算累积百分比(P_Racc)。其中, j为样点个数、R为植被类型、i为非监督分类类型。

大区域植被样本快速提取总体思路, 如图 1所示:

图  1  大区域植被类型样本快速提取方法流程

Figure 1.  Methodological framework of quick sampling in large area

将初始7类样本对应的非监督分类的结果, 依据每一类样本中包含的非监督分类类型及样点个数, 绘制分布直方图, 如图 2所示:

图  2  各类样本非监督分类类型属性特点直方图

Figure 2.  Histogram of samplings with unsupervised feature

统计结果表明, 研究区域内植被以乔木、草地、湿地类为主, 灌木、农田较少。该区域建筑物在非监督分类中特征明显, 与其他地物很好的区分。水体零星分布, 样点较少。本研究以植被类样本纯化提取为主, 建筑物和水体为非植被类型样本不予讨论。

依据像元类型权重信息和空间信息剔除各类样本中不合要求的样点, 进行样本纯化^[16]。保留权重大的更有代表性的样点, 剔除掉权重比例小于2.0%的样点, 剩余的样点累积百分比为70%±5%, 在此样点比例范围内再结合可分离性指标进行样点的取舍。

对保留样点的非监督分类类型、权重、累积百分比进行统计, 如表 1所示:

最后, 从空间信息上判断, 地类边界过渡区域的像元往往都是混合的非纯像元, 因此剔除相邻类型边界1 000 m范围内过渡区域内的样点。

采用Jeffries-Matusita距离(以下简称J-M距离)和转换分离度(Transformed Divergence, 以下简称T-D)作为测度指标来检测所选取的任意两个类别训练样本之间的分离程度。可分离性指标(Separability)取两者中较小值, 即Separability=Min(J-M, T-D)。可分离性指标值为0~2.0, 大于1.9说明样本之间可分离性好, 1.9~1.8属于合格样本, 小于1.8需要重新选择样本, 小于1考虑将两类样本合成一类样本。

如表 2所示, 对于可分离性指标小于1.8的样本, 认为样本中混有其他类型样点, 因此, 删除此类样本中权重比例小的干扰样点, 再计算可分离性指标。反复操作以上过程, 并比较分析可分离性指标, 决定样点取舍。最终得到具有代表性的植被类型样本。

对比纯化前后的结果可以看出纯化后各类样本的可分离性指标提高到1.9以上, 可分离性良好。只有湿地和草地类型可分离性指标值为1.78, 认为湿地区域覆盖的植被和草地区域覆盖的植被多为草本植物, 受实验中特征影像制约辨识度不高。

对比2001年MODIS全年数据的第33天、第145天、第257天、第353天的影像, 纯化前后的各植被样本的平均值和标准差, 如表 3所示。

通过对比4个时相影像的统计结果可知, 纯化前后平均值的变化无显著规律, 纯化后的样本标准差大多缩小。纯化后的样本数据的分布曲线变窄, 重叠变小, 可分性增大, 更易于类型间的区分。

中国植被图表明实验区植被覆盖以乔木、草地和湿地为主。植被覆盖率为乔木67.40%、草地18.04%和湿地12.48%。灌木0.97%, 农田1.08%两类植被在实验区域较少分布。

为验证纯化后样本的空间分布状况, 统计各植被样本在中国土地覆盖图上分布的混淆矩阵, 如表 3、4所示。

表 4中乔木、草地纯化后样本与中国土地覆盖类型中相应植被类型的空间分布基本一致, 精确程度达到87.53%和90.46%。湿地与草地均以草本植被覆盖为主, 两种植被类型在实验影像上不易区分, 因此样本代表性相对较差。乔木、草地、湿地三者空间分布的总体精确程度为84.82%。

表 5中农田样本与中国土地覆盖类型中相应植被类型的空间分布基本一致, 精确程度达到90.32%, 由于实验区灌木分布极少, 覆盖率仅为0.97%, 同时中国土地覆盖图与中国植被图中灌木类空间分布不一致, 导致灌木样本以WESTDC为参考的空间分布状况精确程度较低。但灌木样本可分离性指标高于1.9, 不影响与其他类型的区分, 故认为灌木样本基本不会影响到各类植被覆盖信息的表达。

实验样本中的植被空间分布与中国土地覆盖图中植被的空间分布相似度超过80%, 因此认为纯化后的样本在该区域具有代表性。

将纯化前后的样本, 分别输入最大似然分类器, 得到实验区域的分类结果, 对比中国植被图图 3(a)中的植被分布, 如图 3所示。

图  3  样本纯化前后分类结果与中国植被图中植被分布的比较分析

Figure 3.  Comparisons of classification results with (a) (b) (c)

验证样点在中国植被图和中国土地覆盖图中植被类型相同的区域随机选取。纯化后样本的最大似然分类结果图 3(c), 总体精度为76.78%;相比未经过纯化的初始样本最大似然分类结果图 3(b), 总体精度44.26%, 提高了32.52%。

由于特征影像、分类器选择会对分类精度造成影响, 验证数据本身也与真实地物存在一定的精度局限性等, 以上因素都可能降低实验影像最终分类结果的精度。因此, 在三者相同的情况下, 纯化后的样本分类结果精度较大提高, 纯化后样本的分类结果与中国植被图中植被分布的总体趋势能够达到基本一致, 认为采样方法是可行的。

(1) 以遥感监督分类方法快速精准的提取大区域植被覆盖信息为研究背景, 本研究结合专家经验和实验数据同质特征向量的聚类信息, 采用统计分析的方法, 实现了针对宏观大区域遥感影像快速提取植被样本的方法。

(2) 实验方法增强了样本的可分离性、缩小了样本标准差, 经验证纯化后样本对实验区域的植被覆盖信息具有代表性, 采样方法可行。所得样本具有实验影像时效特征, 用于遥感监督分类可以缩短大区域植被变化监测时间间隔。

(3) 实验节省了大区域植被类型调查耗费的人力物力资源和时间, 避免了引入多源数据造成的数据冗余和误差影响, 提高了采样效率。