第一作者简介: 高闪闪(1987-),女,硕士研究生,主要从事遥感图像处理与分析方面的研究。E-mail:gss1127@sina.cn。
针对现有植被指数不适用于ALOS图像植被信息提取的问题,从分析植被的光谱特征入手,提出了一种基于植被样本的植被指数(vegetation sample-based vegetation index,VSVI),并通过数学公式推导证明了VSVI仅与植被的光谱信息有关,与土壤背景无关,具有一定的消除土壤背景影响的能力。利用该植被指数和采用阈值分割方法提取了南京市某区域ALOS图像中的植被信息,并与差值植被指数,比值植被指数、归一化差值植被指数及土壤调节植被指数等其他植被指数的植被信息提取结果进行了比较。研究结果表明,该文提出的VSVI植被指数能够克服其他植被指数的缺点,植被信息提取精度分别提高了21.7%,27.5%,14%和9.5%。
In view of the phenomenon that the existing vegetation indexes are not suitable for vegetation extraction of ALOS image, this paper, starting with an analysis of the spectral characteristics of the vegetation, puts forward a new vegetation index(vegetation sample-based vegetation index,VSVI) based on the analysis of vegetation samples and proves that this vegetation index is only associated with the spectral information of vegetation but not related to the soil background, thus having a certain capability of eliminating the soil background image with mathematical derivation. The vegetation of ALOS image is extracted by the vegetation index with the method of threshold segmentation and compared with the vegetation indexes(DVI, RVI, NDVI and SAVI). The results show that the vegetation index is capable of overcoming the shortcomings of other vegetation indexes, and the vegetation extraction accuracy can be raised by 21.7%, 27.5%, 14% and 9.5% respectively.
ALOS(advanced land observation satellite)卫星是日本国家空间发展局(NASDA)2006年发射的一颗陆地观测卫星[1], 用于快速获取高精度的地面数据[2]。ALOS卫星携带了3种传感器: 全色立体测图传感器、新型可见光和近红外辐射计AVNIR- 2及相阵型L-波段合成孔径雷达。在波段设置和空间分辨率上, ALOS AVNIR-2 与SPOT5 CCD有很大的相似性。何宇华等[3]针对ALOS 数据图像质量分析了影像纹理特征, 评定了平面点位、影像分类和图斑面积精度, 并将其与SPOT5数据对比, 发现ALOS数据在纹理特征及面积精度方面均明显优于SPOT5数据。这也是近年来ALOS图像被广泛应用于遥感领域研究的重要原因之一。
目前, 关于利用ALOS卫星图像提取植被信息的研究主要采用纹理辅助、支持向量机、决策树和最大似然分类等方法[4, 5, 6, 7], 而植被指数简单实用的特点使其一直是植被遥感研究的重要内容之一[8]。随着遥感技术的发展, 为了满足不同的应用需求, 用于提取植被信息的各种植被指数应运而生, 如差值植被指数(difference vegetation index, DVI)[9]、比值植被指数(ratio vegetation index, RVI)[10]、归一化差值植被指数(normalized difference vegetation index, NDVI)[11]、土壤调节植被指数(soil adjusted vegetation index, SAVI)[12]、全球环境监测指数[13]和三波段梯度差植被指数[14]等。虽然这些植被指数具有一定消除土壤背景或大气影响的能力, 但是利用这些植被指数提取ALOS图像植被信息时存在将建筑物错认为是植被信息和漏提植被信息的现象, 或者植被密度较低的(如农田和草地等)不易提取, 而且消除土壤背景影响能力较差。
针对上述情况, 本文通过分析植被样本数据, 发现植被的光谱数据间存在某种较为固定的关系, 根据这种关系提出了适于ALOS图像的基于植被样本的植被指数(vegetation sample-based vegetation index, VSVI), 并通过植被提取实验证明了VSVI不仅具有很好消除土壤背景影响的能力, 而且植被提取效果好于其他植被指数。
研究区位于南京市浦口区, 地处南京市长江西北岸, E 118° 20'~119° 13', N 31° 14'~ 32° 17'之间; 中部有老山山脉, 西部丘陵起伏, 地貌类型多样化, 集低山、丘陵、平原、岗地、大江和大河为一体。全区总面积约913 km2(其中丘陵山区面积632.7 km2, 平原区面积269.3 km2), 属于亚热带季风气候。区内水系发达, 植被类型丰富, 适合植被信息提取研究。
选择2009年10月获取的ALOS AVNIR-2的level1B2级遥感图像。该图像包括红、绿、蓝和近红外4个多光谱波段, 空间分辨率为10 m。
本文技术路线为: ①对图像进行裁剪、线性拉伸及多尺度分割处理, 并计算图像对象的4个波段的光谱特征均值; ②选取100个植被对象为样本, 对其4个波段光谱特征值进行数据统计与分析; ③建立模型函数, 根据对植被样本的计算分析确定函数模型参数l1, l2, l3和l4的值; ④通过对研究区植被信息的提取, 检验VSVI植被指数是否可行, 并通过与其他植被指数的提取结果的比较进行分析和评价。
本文从分析植被光谱的特点入手, 提出一种基于植被样本光谱特征值数理统计分析的植被指数VSVI; 通过数学公式推导证明该植被指数具有一定消除土壤背景影响的能力, 并通过利用VSVI 从ALOS图像中提取植被信息的实验予以验证; 最后采用其他典型植被指数, 即DVI, RVI, NDVI和SAVI(表1)[9, 10, 11, 12], 分别从ALOS图像中提取了植被信息, 并与利用VSVI提取的效果进行了比较。
研究发现, 植被在不同波段的反射率以不同形式组合一个函数时, 该函数与植被指数有着密切关系。为了研究这种关系, 本文统计分析了ALOS图像中植被在蓝光B、绿光G、红光波段R和近红外波段NIR的光谱特征值范围, 分别为ρ B∈ [121, 135], ρ G∈ [92, 115], ρ R∈ [66, 100]和ρ NIR∈ [84, 120]。当2个波段的光谱特征值存在交叉范围时, 该2个波段比值的分布范围表现为比较零散或者差值正负交替变化。表2列出随机选取的100个较为纯净的植被样本的光谱特征值和比值。
由表2可以看出, 可见光蓝、绿和红波段的植被光谱特征值的标准差较小, 内部变化稳定; 而近红外波段的植被光谱特征值的标准差较大, 与其他波段的比值较为零散, 并且与绿光和红光波段的差值运算正负交替, 导致统计参数时不稳定。因此, 不适合选为统计计算VSVI参数的波段。为了避免这种情况并使植被指数比较稳定, 本文只对可见光蓝、绿和红波段的植被光谱特征值进行统计分析。通过对表2中植被样本光谱特征值的计算分析发现: 不同波段的植被反射率以一定形式可以组合成一个比较稳定的参数, 并且这个参数与植被特性有密切关系。
叶绿素在0.55μ m(绿)波长处为反射峰, 在0.45μ m(蓝)和0.67μ m(红)为吸收带, 在近红外0.7μ m, 0.8μ m, 1.1μ m处为反射陡坡, 4个波段的光谱信息都能一定程度上反映植被信息。假设VSVI与不同波段呈线性关系, 建立四元一次方程, 即
VSVI=l1ρ B+l2ρ G+l3ρ R+l4ρ NIR, (1)
式中l1, l2, l3和l4为待确定的参数。
通过计算分析, 发现ρ R/ρ G, (ρ R-ρ B)/ρ R, (ρ G-ρ R)/ρ G和(ρ B-ρ G)/ρ B能够组合成一个比较稳定的参数。因此, 经多次从ALOS图像中提取植被信息的实验, 表明l1, l2, l3, 和l4分别为ρ R/ρ G, (ρ R-ρ B)/ρ R, (ρ G-ρ R)/ρ G和(ρ B-ρ G)/ρ B时, 植被提取效果最佳。因此,
VSVI=0.75ρ B-0.71ρ G+0.26ρ R-0.21ρ NIR 。 (2)
VSVI具有消除土壤背景影响的能力。假设在可见光蓝、绿、红和近红外4个波段, 植被和土壤面积的比值不随波段变化(尽管该假设比较粗糙, 但一些通用的植被指数在一定程度上也隐含了该假设), 像元的反射率为像元内植被和土壤的面积加权和[7]。若忽略误差项, 则4个波段的反射率分别为
ρ B=Aρ VB+(1+A)ρ SB, (3)
ρ G=Aρ VG+(1+A)ρ SG, (4)
ρ R=Aρ VR+(1+A)ρ SR, (5)
ρ NIR=Aρ VNIR+(1+A)ρ SNIR, (6)
式中: ρ B, ρ G, ρ R及ρ NIR分别为蓝光、绿光、红光及近红外波段的反射率; ρ Vi和ρ Si分别为i波段的植被及土壤反射率(i=B, G, R, NIR); A为植被覆盖度。
式(3)— (6)代入式(2), 化简后得到
VSVI=A(k1-k2)+k2; (7)
式中:
k1=0.75ρ VB-0.71ρ VG+0.26ρ VR-0.21ρ VNIR; (8)
k2=0.75ρ SB-0.71ρ SG+0.26ρ SR-0.21ρ SNIR。 (9)
可见在上述假设条件下, VSVI仅与植被的光谱信息有关, 与土壤背景无关。
植被提取实验中涉及到图像裁剪、波段组合、多尺度分割、波段运算和构造植被提取函数模型等技术。为了证明VSVI具有消除土壤背景影响的能力, 本文采用多阈值分割方法提取植被信息。图1为利用不同植被指数提取植被的结果。
实验结果表明, 利用NDVI, DVI, RVI和SAVI等植被指数均能将植被轮廓提取出来, 但存在一定的漏提或误提等缺点; 只有利用VSVI提取植被克服了其他植被指数的上述缺点。
从图1中可以看出: ①与对原始ALOS图像(图1(a))的植被(红色区域)目视解译结果对比, 图1(b)— (e)中都存在植被信息的漏提和误提现象(蓝色椭圆内为植被信息漏提, 红色椭圆内表示建筑物错提为植被信息, 黄色椭圆内的颜色表示植被密度比较低的植被信息); ②植被提取的实验结果与阈值的选取有关, 在构造模型时阈值的选取要经过多次实验才能确定, 在图1(b)— (e)中, 要么植被密度比较低的农田和草地不易被分割出来, 要么植被与道路和裸地的交叉现象比较严重; ③为了验证VSVI指数的精度, 多尺度分割后统计了植被信息的总对象数为18 769, 并计算了NDVI, DVI, RVI, SAVI和VSVI的提取精度(表3)。与其他植被指数相比, VSVI指数的分类精度大幅提高。
通过分析发现: NDVI, DVI, RVI及SAVI虽然一定程度上能够消除土壤背景的影响, 但是在植被密度较低的情况下植被提取效果不太理想, 并且DVI和RVI存在漏提和错提、NDVI和SAVI有错提现象; 从图1(f)来看, VSVI植被指数克服了其他植被指数的上述缺点, 并且能够较好地消除土壤背景和薄云的影响。实验证明: 本文提出的VSVI植被指数与其他植被指数相比, 从ALOS图像中提取植被的效果更为理想。
1)从植被光谱特征出发, 通过分析ALOS卫星AVNIR-2的4个多光谱波段之间植被样本特征值之间的数学关系, 提出基于植被样本的植被指数(VSVI), 并通过公式推导证明其具有消除土壤背景影响的一定能力。
2)与其他典型植被指数相比, 使用VSVI的植被信息提取效果最好。研究发现, 参数l1, l2, l3和l4的值只要在一定范围内变化, 就不存在对植被信息提取效果的影响。
3)VSVI也有一定的缺陷, 如由于植被样本选取的差异会导致参数值的上下波动, 但总体上并不影响植被信息的提取效果; 由于本文提出植被指数仅局限于对ALOS图像的植被提取, 如何把该植被指数应用于其他卫星数据是今后研究的重点。另外, 植被的提取受阈值选择的影响, 如何更准更快地确定阈值也是今后研究的问题之一。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|