全 文 :书基于航空高光谱影像的青海省玛多县
高寒草原景观格局特征分析
焦全军1,2,张兵1,赵晶晶1,2,刘良云1,2,胡勇1,2,3
(1.中国科学院对地观测与数字地球科学中心,北京100094;2.中国科学院数字地球重点实验室,北京100094;
3.西安科技大学测绘科学与技术学院,陕西 西安710054)
摘要:景观生态学在草地资源管理、生态安全评价方面已得到广泛应用,开展不同尺度的景观格局特征分析对景观
生态过程研究十分重要。获取不同覆盖度草地的景观格局特征是青藏高原草地退化等生态系统演变研究的基础
工作之一。从高分辨率遥感的角度提取中小尺度上草地景观格局特征,是开展景观格局尺度分析的有益尝试。本
研究通过高空间分辨率的PIS航空高光谱影像,分析不同草地类型的光谱特征差异,开展玛多实验区草地类型精
细分类,并从景观异质性、景观优势度和景观破碎度3个景观层次,分析米级分辨率下不同覆盖度高寒草原样地的
景观格局特征差异,发现高寒草原样地的景观优势度、景观异质性与草地覆盖度分别呈正、负相关,而破碎度指标
则随着高寒草原草地覆盖度的降低出现先升高后降低的现象。
关键词:黄河源;景观格局;高寒草原;航空高光谱
中图分类号:S81205 文献标识码:A 文章编号:10045759(2012)02004308
景观生态学以景观结构、功能和动态特征为主要研究对象,开展不同时空尺度下格局及过程、人类作用为主
导的景观演化研究[1,2]。针对草地退化、荒漠化方面的问题,景观生态学在宏观尺度下的草地资源管理、生态评
价中具有重要作用[3,4]。作为青藏高原的主要生态类型,草地生态系统对青藏高原的生态安全、经济效益十分重
要。数十年来,在气候变化和人为干扰的双重作用下,青藏高原腹地的黄河源区出现草地覆盖度下降、土地荒漠
化加重的问题。以植被覆盖度分级的黄河源草地退化监测及其景观分析研究已得到广泛重视[5,6]。
景观生态学强调空间异质性、等级结构和尺度在研究生态学格局与过程中的重要性,不同尺度的景观格局差
异以及转换关系能够与植被生态过程联系起来,景观格局与生态过程的多尺度研究为景观生态学的核心方向之
一[79]。遥感技术能够为景观分析提供地物覆盖类型的空间分布数据,已被广泛用于景观格局分析研究。当前景
观生态研究存在的基础数据积累不足的问题,阻碍了景观格局与生态功能的关联研究,原因之一为遥感数据的利
用程度不足,特别是多尺度遥感数据的使用尤其缺乏[10]。景观格局的尺度研究需要不同尺度空间数据的支
持[11],当前的草地景观分析研究以利用LandsatTM 图像[12,13]和 MODIS影像[14]为主,基于高空间分辨率遥感
影像的中小尺度景观分析则鲜有开展。借助高分辨率遥感数据开展植被景观分析,能够提供中小尺度下景观格
局数据,有助于从景观分析的角度更深入理解不同景观格局特性对应的生态过程[1517]。
遥感分类误差对景观分析影响极大,高精度的植被类型遥感分类结果对景观分析至关重要[18]。运用多光谱
遥感手段开展青藏高原草地精细分类尚存在一定困难,高光谱遥感则能够探测更精细的植被理化属性[19],并有
助于提高草地类型的分类精度[2022]。在景观分析中使用高光谱遥感的地物类型空间填图,能够提高景观分析结
果的准确性。本研究将基于青藏高原黄河源星-机-地同步实验获取的航空高光谱影像,在草地类型精细识别
的基础上,运用景观生态学的方法开展不同覆盖度草地的景观格局特征分析,为开展青藏高原高寒草原退化监测
工作提供可靠的数据支持和科学依据。
第21卷 第2期
Vol.21,No.2
草 业 学 报
ACTAPRATACULTURAESINICA
43-50
2012年4月
收稿日期:20101229;改回日期:20110308
基金项目:国家重点基础研究发展计划 (2009CB723902)和国家自然科学基金项目 (41001263)资助。
作者简介:焦全军(1981),男,河南巩义人,助理研究员,博士。Email:qjiao@@ceode.ac.cn
通讯作者。Email:zhangbing@ceode.ac.cn
1 材料与方法
1.1 实验区概况
玛多县属青海省果洛藏族自治州,位于青海省南部,地处青藏高原腹地的黄河源区,是黄河流经的第一县,地
理位置为33°50′~35°40′N,96°55′~99°20′E。气候属青藏高原气候系统,为典型的高原大陆性气候。实验区河
流密布,湖泊、沼泽众多,主要分布着高山草甸和高寒草原,沿河及湖泊周围低洼区域分布有藏北嵩草沼泽。高寒
草原是青藏高原分布面积最大、地带性最为典型的草地类型。草地覆盖程度是研究草地退化程度的重要指
标[23],本研究选择玛多县高寒草原不同覆盖度草地为研究对象。结合实地考察和青藏高原高寒草原调查的相关
研究成果[5,2426],将高寒草原区细分为高覆盖草地(覆盖度大于50%)、中覆盖草地(覆盖度为30%~50%)、低覆
盖草地(覆盖度低于30%),同时实验区存在沙化区域,定义为沙化稀疏杂草地[6],将实验区地物分为高寒草原高
覆盖草地、中覆盖草地、低覆盖草地以及沙化稀疏杂草地、沙地、水体六类。
1.2 高光谱数据获取与处理
于2009年8-9月实施了青藏高原黄河源玛多星-机-地同步观测实验,地面开展了植被类型调查、覆盖度
调查及自动气象站实况数据采集,以鄂陵湖北岸的紫花针茅(犛狋犻狆犪狆狌狉狆狌狉犲犪)高寒草原区域为航空飞行样区,利
用航空高光谱仪器和地面光谱仪,获得航空高光谱影像和地面光谱数据。地面光谱仪采用 ASD公司生产的
Fieldspec3光谱仪,光谱范围为350~2500nm,光谱波段采样间隔为:1.4nm(350~1050nm)、2nm(1000~
2500nm)。航空高光谱飞行样区分为28个航带,共飞行了10个架次。航空高光谱仪器为面阵推扫式成像光谱
仪(pushbroomimagingspectrometer,PIS)[27],仪器参数见表1。
表1 面阵推扫式成像光谱仪(犘犐犛)的基本参数
犜犪犫犾犲1 犘狉犻犿犪狉狔犮犺犪狉犪犮狋犲狉犻狊狋犻犮狊狊狌犿犿犪狉狔狅犳犘狌狊犺犫狉狅狅犿犐犿犪犵犻狀犵犛狆犲犮狋狉狅犿犲狋犲狉(犘犐犛)
仪器Instruments CCD像元数NumberofCCDpixels 空间分辨率Spatialresolution(m)波段范围 Wavelength(nm)波段数Numberofband
PIS 1400 2 450~900 800
航空高光谱影像经几何校正和大气校正2个处理步骤,获得地表反射率数据。几何校正充分利用航空POS
系统,通过分析相机坐标系、惯性测量组件(inertialmeasurementsunits,IMU)坐标系、局地地理坐标系、WGS84
(worldgeodeticsystem)大地坐标系、投影坐标系之间的关系,分别建立直接法和间接法几何纠正算法[28],获得
鄂陵湖北岸航带的高寒草原实验区的2009年8月30日航空高光谱影像校正结果(图1)。基于大气辐射传输模
型的大气校正方法已比较成熟,本研究使用ENVI4.5的FLAASH(fastlineofsightatmosphericanalysisof
spectralhypercubes)大气校正模块,基于当天的实测大气参数,开展了航空高光谱影像的大气校正并获得反射率
数据。
本次研究利用ENVI4.5的光谱角度填图(spectralanglemapping,SAM)方法进行实验区高光谱遥感草地
类型填图。SAM是一个基于自身的波谱分类方法,通过计算地物基准光谱与影像光谱之间的广义夹角来表征其
匹配程度,匹配的夹角越小,表明地物与基准光谱越相似[29]。光谱角度匹配法能够反映出波形的相似度,适合开
展高光谱遥感的地物类型分类[30,31]。在地物光谱先验知识的支持下,角度匹配法具有较高的地物分类精度。
1.3 景观格局指数
本研究根据草地覆盖度划分差异化的样区,开展景观格局特征分析。根据实地考察和草地覆盖度航空高光
谱反演结果,选择不同植被覆盖度状况的草地样区,并对不同植被覆盖度的样区进行景观格局指数的计算。首
先,在经过几何校正和大气校正后的航空高光谱影像基础上,将利用归一化植被指数(normalizeddifferenceveg
etationindex,NDVI)的二分像元模型估算植被覆盖度(fractionalvegetationcover,FVC)[32]:
犉犞犆= 犖犇犞犐-犖犇犞犐狊狅犻犾犖犇犞犐狏犲犵-犖犇犞犐狊狅犻犾
(1)
式中,NDVIveg表示完全覆盖植被的NDVI值,NDVIsoil表示无植被覆盖裸地的NDVI值。
44 ACTAPRATACULTURAESINICA(2012) Vol.21,No.2
图1 玛多实验区的航空高光谱影像(假彩色显示)
犉犻犵.1 犜犺犲犳犪犾狊犲犮狅犾狅狉犱犻狊狆犾犪狔狅犳犪犻狉犫狅狉狀犲犺狔狆犲狉狊狆犲犮狋狉犪犾犻犿犪犵犲狉狔犻狀犕犪犱狌狅狊狋狌犱狔犪狉犲犪(犳犪犾狊犲犮狅犾狅狉犱犻狊狆犾犪狔)
以景观几何特征为基础的景观格局指数可以有效地反映景观生态过程[2]。通过划分不同植被覆盖度的研究
样区,并在景观格局指数选取上着重于从景观层次的景观异质性、景观优势度和景观破碎度3个角度,来表征不
同覆盖度高寒草原样地的景观格局特征。
(1)景观异质性。从景观内部包含类型的多样性及各类型分布结构的复杂性来分析景观异质性,利用多样性
指数和均匀度指数来表征。
1)Shannon多样性指数(Shannon’sdiversityindex,SHDI):
犛犎犇犐=-∑
狀
犽=1
犘犽ln(犘犽) (2)
式中,Pk为斑块类型,犽为在景观中出现的概率(通常以该类型占有的栅格细胞或像元数占景观栅格细胞总数的
比例来估算),狀为景观中斑块类型的总数。多样性指数建立在信息论基础之上,往往用来度量系统结构组成的
复杂程度。
2)Shannon均匀性指数(Shannonevennessindex,SHEI)
均匀度指数反映景观中各斑块在面积上分布的不均匀程度,通常以多样性指数和其最大值的比来表示。以
Shannon均匀度指数为例,均匀度可表达为:
犛犎犈犐= 犛犎犇犐犛犎犇犐犿犪狓=
-∑
狀
犽=1
犘犽ln(犘犽)
ln(狀)
(3)
式中,犛犎犇犐是Shannon多样性指数,犛犎犇犐max是其最大值,即每种生态类型所占的景观面积比例都是1/狀。显
然,当犛犎犈犐趋于1时,景观斑块分布的均匀程度也趋于最大。犛犎犈犐值越大,代表景观要素类型愈丰富,景观
多样性越大;犛犎犈犐值越小,代表景观中各斑块在面积上分布越不均匀。在样区内类别狀确定后,犛犎犈犐与犛犎
犇犐存在定量的转换关系。
(2)景观优势度。景观类别的优势度通过最大斑块指数百分比(largestpatchindex,LPI)来表示,分别在类
层和景观层上计算,类层计算结果与景观层计算结果相同的类为优势度最大的类别。在类层,LPI指景观内每类
别的最大斑块面积与景观总面积的百分比,即第1类都有1个LPI;在景观层,LPI指所有类别中最大的那个斑
块面积与景观总面积的百分比,即景观层的LPI只有1个。景观优势度决定着景观中的优势种类及内部种类的
54第21卷第2期 草业学报2012年
丰度等生态特征。
(3)景观破碎度。景观的破碎程度通过景观平均斑块面积(areamean,AM)和周长面积分形维数指数(per
imeterareafractaldimension,PAFD)来表示。
1)平均斑块面积(AM)的定义为:景观层所有类别的斑块面积总和与斑块数量总和之比,反映景观区域内平
均斑块大小,从一定程度上来指示景观破碎度的大小。景观平均斑块面积越小,景观破碎程度越大。
2)周长面积分形维数指数(perimeterareafractaldimension,PAFD):是景观空间格局复杂程度的度量,通过
计算某一斑块形状与相同面积的圆或正方形之间的偏离程度来测量其形状的复杂程度。PAFD值越大,代表景
观内部结构越复杂,景观破碎程度越高。
2 结果与分析
2.1 光谱分析及分类结果
基于2009年8-9月地面采集的光谱数据,同时辅以PIS影像提取草地和其他地物的反射率光谱(图2),比
较不同草地类型和其他地物类型的反射率光谱,发现高寒草原高覆盖草地的红光波段吸收强,近红外波段反射
高;高寒草原中覆盖草地的红光吸收强,近红外反射则不高;高寒草原低覆盖草地在可见光波段的反射率较高覆
盖草地高,光谱的“红谷”不明显,而近红外波段的反射率较高覆盖草地低,原因是矮火绒草(犔犲狅狀狋狅狆狅犱犻狌犿狀犪
狀狌犿)叶片呈灰白色,其在高寒草原低覆盖草地中比例相对较大,因此矮火绒草光谱对高寒草原低覆盖草地光谱
的影响较大;土壤沙质化的稀疏杂草草地受到沙地反射率的影响较大,同时杂草类叶片的叶绿素含量较高,红外
波段的反射率高于高寒草原低覆盖草地的反射率;沙地的反射率从可见光至近红外波段呈上升趋势,其近红外波
段的反射率低于高寒草原高覆盖草地和沙化稀疏杂草草地;水体在近红外波段的反射率趋近于0。基于实验区
的地物光谱分析,发现6种地物间存在一定的光谱特征差异,从而为在黄河源开展基于高光谱数据的草地类型精
细填图提供了可能。
在光谱提取和分析的基础上,利用光谱角度匹配方法,通过ENVI4.5进行实验区PIS高光谱影像的草地类
型精细填图(图3a)。地面调查和影像人工解译共获得各类测试样本集,通过计算混淆矩阵的方法获取相对标准
类别底图的总体精度(84.3%),说明光谱角度匹配法适合开展高光谱遥感影像的地物类型精细分类。结果显示,
高寒草原的高覆盖草地主要分布于地势较缓的山地,并且临近河流;高寒草原中覆盖草地分布于大面积的平缓区
域,与河流和湖泊距离较远,并存在大量的低覆盖草地斑块。经过统计得出,航区飞行区内中、低覆盖草地占
14%左右,稀疏杂草地占5%左右。稀疏杂草草地主要出现在高寒草原的严重退化区域和河漫滩区域。
2.2 景观格局指数分析
基于航空遥感影像提取裸地NDVIsoil和完全覆盖草地NDVIveg,代入植被覆盖度二分模型中(公式1),获得
实验区草地覆盖度填图结果。根据高寒草原草地覆盖度的差异,以中覆盖草地为基质层,在鄂陵湖北岸的高寒草
原实验区选取草地覆盖度依次降低的3个样区(图3),分别为样区1(平均覆盖度:0.37)、样区2(平均覆盖度:
0.31)和样区3(平均覆盖度:0.28)。根据实验方法中涉及的3类景观格局指数,使用fragstats软件运行得到每
个样区的景观格局指数计算结果,按类层和景观层统计结果如表2。
2.3 景观格局特征分析
(1)景观异质性。3个样区的景观层SHDI及SHEI的结果显示,植被覆盖度指标最低的样区3的景观多样
性丰富度值最高(0.7860),且景观中各斑块在面积上分布相对其他2个样区更均匀(0.7154),异质性最强;植被
覆盖度指标最高的样区1的景观多样性丰富度最低(0.5501),且景观中各斑块在面积上分布非常不均匀
(0.5007),异质性最弱;样区2的景观多样性丰富度值(0.7306)稍弱于样区3,其景观中各斑块在面积上分布均
匀性值(0.6651)介于样区1和样区3之间,说明样区2的景观异质性居中。
(2)景观优势度。景观优势度有助于确定样区的优势类型和描述优势类比例状况。通过景观层的优势度
LPI值与基质层的LPI值对比,高寒草原实验区3个样区的中覆盖草地基质层均保持着优势。样区1的基质层
占居绝对优势,其LPI值高达0.7745;植被覆盖度居中的样区2基质层优势度有所下降,其LPI值仅为0.6611;
植被覆盖度最低的样区3内虽基质层仍占居优势,但基质层的优势度有了大幅的降低,仅为0.4445,较植被覆盖
64 ACTAPRATACULTURAESINICA(2012) Vol.21,No.2
度最高的样区1基质层的LPI值下降了42.61%。景观优势度指标和覆盖度指标关系密切,由于3个样区的基
质层都为中覆盖草地类,相对低覆盖草地和稀疏杂草草地类型,基质层的植被覆盖度较高,因此,本实验区内样区
景观优势度越高,代表了样区的植被覆盖度值越大。
图2 实验区典型地物的反射率光谱曲线
犉犻犵.2 犛狆犲犮狋狉犪狅犳犱犻犳犳犲狉犲狀狋犾犪狀犱犮狅狏犲狉狋狔狆犲狊犻狀狊狋狌犱狔犪狉犲犪
图3 基于犘犐犛数据的草地类型填图(犪)与不同覆盖度的实验样区(犫)
犉犻犵.3 犌狉犪狊狊犾犪狀犱狋狔狆犲犿犪狆狆犻狀犵犫犪狊犲犱狅狀犘犐犛犻犿犪犵犲狉狔(犪)犪狀犱狊犪犿狆犾犲犪狉犲犪狊狑犻狋犺犱犻犳犳犲狉犲狀狋犳狉犪犮狋犻狅狀犪犾狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉(犫)
(3)景观破碎度。从景观层和景观类型2个层次分析不同样区的景观破碎度。从样区景观层的破碎度指标
来看,3个样区基质层的景观平均斑块面积(AM)受到非基质层的影响,数值都较小,其中覆盖度居中的样区2的
AM值最低,说明样区2的斑块数目最多,而样区2的PAFD值最小,总体破碎度较高。样区中基质层的景观破
碎度是重点关注的对象,样区2基质层的AM值最大,为1.0335,这是由于在基质层退化的同时保存了较大的
斑块,另外其PAFD值相对样区3和样区1为最小,显示基质层的形状异质结构最为复杂;植被覆盖度最低的样
区3基质层的AM值最小,仅为0.2147,基质层斑块最为破碎,而PAFD值高于样区2,与样区1的PAFD值接
近,显示其基质层的形状结构变得相对简单。总的来说,样区1的基质层景观破碎度居中,样区2的基质层景观
破碎度最小,样区3的基质层景观的破碎度最大。
74第21卷第2期 草业学报2012年
表2 不同样区景观层及主要草地类型的景观格局指数统计结果
犜犪犫犾犲2 犜犺犲狊狋犪狋犻狊狋犻犮狊狅犳犾犪狀犱狊犮犪狆犲狆犪狋狋犲狉狀犻狀犱犻犮犲狊犻狀犱犻犳犳犲狉犲狀狋狊犪犿狆犾犲犪狉犲犪狊
景观层次与类别
Landscapelevelandtypes
景观格局指数
Landscapepatternindex
样区1
SampleareaNo.1
样区2
SampleareaNo.2
样区3
SampleareaNo.3
景观层Landscapelevel SHDI 0.5501 0.7306 0.7860
SHEI 0.5007 0.6651 0.7154
LPI 0.7745 0.6611 0.4445
AM 0.1264 0.1183 0.1300
PAFD 1.2677 1.2313 1.2720
高寒草原中覆盖草地类
Moderatecoveralpinesteppe
LPI 0.7745 0.6611 0.4445
AM 0.6706 1.0335 0.2147
PAFD 1.3135 1.3876 1.3257
高寒草原低覆盖草地类
Lowcoveralpinesteppe
LPI 0.0001 0.0617 0.1168
AM 0.0074 0.0615 0.1478
PAFD 1.2101 1.2248 1.2806
沙化稀疏杂草类
Alpinedesertifiedgrassland
LPI 0.0670 0.0046 0.0039
AM 0.0369 0.0095 0.0107
PAFD 1.2631 1.2037 1.1689
SHDI:Shannon多样性指数Shannon’sdiversityindex;SHEI:Shannon均匀性指数Shannonevennessindex;LPI:最大斑块指数百分比 Largest
patchindex;AM:平均斑块面积 Areamean;PAFD:周长面积分形维数指数Perimeterareafractaldimension.
综上所述,从样区景观层的整体而言,植被覆盖度越低,以中覆盖草地为基质层的高寒草甸景观优势度越低,
而景观异质性则越高,覆盖度居中的样地2破碎度最高;从中覆盖草地基质层来看,植被覆盖度指标最低的样区
3中基质层景观内部结构多样化程度最高,基质层的优势度处于消失的边缘,生态环境承载力最弱;样区2的基
质层景观的内部形状结构多样化程度最高,破碎度最大,基质层的优势度仍然明显,生态环境承载力较强;样区1
的基质层优势度最高,多样性最弱,而其破碎度较样区2较好,生态环境承载力最好。从3个样区在整个实验区
的分布状态来看,该区的草地生态环境承载力从由东北至西南逐渐下降,这与样区与鄂陵湖及河流的距离有较强
的相关性,土壤水分条件空间分布是高寒草原不同区域植被生长状况差异的重要因素[33]。
3 结论与讨论
利用航空高分辨率遥感数据,尝试用于青藏高原典型的高寒草原景观格局特征分析,相对于LandsatTM 影
像的大尺度景观格局分析来讲,是在中小尺度上草原生态系统的景观分析的重要补充。通过高寒草原同一区域
内不同覆盖度草地样区的景观格局指数分析,发现草地生长状况与景观格局指数之间存在一定的规律,表现在:
样区基质层的优势度随着高寒草原覆盖度的降低而逐渐减低,覆盖度越低其异质性越高;另外,3个样区的破碎
度统计显示,随着样区覆盖度的降低,破碎度出现先升高后降低的现象,分析其原因,一方面是中覆盖草地基质层
中出现大量低覆盖草地斑块,另一方面覆盖度持续降低,作为非基质层的低覆盖草地和更为稀疏的沙化杂草草地
小斑块出现合并,即草地总体覆盖度、景观格局特征的变化是和草地内生态过程发展紧密联系的。草地退化过程
中的覆盖度变化往往伴随着草地内不同斑块形态变异[34],开展草地盖度和景观格局指数的综合分析,有利于发
挥空间观测技术监测大范围的草地退化状况。综合不同覆盖度草地的景观格局特征,可能为评价草地退化程度
提供新的测算依据,当然,如何选择和构建景观格局指数诊断草地生态过程,还需要很长的路要走。
高光谱遥感有助于提高地物类型遥感分类精度,从而保障景观分析具有坚实的数据基础。实验区内不同生
长状态的高寒草原草地、稀疏杂草草地等类型存在一定的光谱差异,其光谱差异是基于高光谱遥感影像的草地类
型填图可行的基础条件。本实验时间接近8月底,所以本组光谱分析的结果也具有明显的季相特征,在推广应用
中需要考虑草地光谱的时相变异[20]。
84 ACTAPRATACULTURAESINICA(2012) Vol.21,No.2
分析了中小尺度上的玛多高寒草原不同退化程度样区的景观格局特征,而不同分辨率遥感影像的景观格局
特征间差异与尺度转换关系以及生态学过程含义,还需要深入研究。在地面层次开展不同空间尺度下的样地调
查,获取空间尺度连续变化的景观格局指数曲线[9],则能够为草地景观格局特征遥感监测提供先验知识和生态学
过程解释。另外本研究关注单一时相的草地样区,开展中小尺度上高寒草原景观变化监测对于研究草地演替与
气候变化、人类活动的关系十分重要,需要寻找合适的数据开展进一步研究。
参考文献:
[1] TurnerM.Spatialandtemporalanalysisoflandscapepatterns[J].LandscapeEcology,1990,4(1):2130.
[2] 邬建国.景观生态学———格局、过程、尺度与等级[M].北京:高等教育出版社,2001.
[3] 田育红,刘鸿雁.草地景观生态研究的几个热点问题及其进展[J].应用生态学报,2003,14(3):427433.
[4] 陈全功.江河源区草地退化与生态环境的综合治理[J].草业学报,2007,16(1):1015.
[5] 王根绪,丁永建,王建,等.近15年来长江黄河源区的土地覆被变化[J].地理学报,2004,59(2):163173.
[6] 王辉,袁宏波,徐向宏,等.黄河源区沙化草地区域景观特征[J].应用生态学报,2006,17(9):16651670.
[7] GaoQ,YangX.Arelationshipbetweenspatialprocessesandapartialpatchinessindexinagrasslandlandscape[J].Land
scapeEcology,1997,12(5):321330.
[8] 傅伯杰,吕一河,陈利顶,等.国际景观生态学研究新进展[J].生态学报,2008,28(2):798804.
[9] 鲁为华,朱进忠,王东江,等.天山北坡两种退化梯度下伊犁绢蒿种群空间分布的点格局分析[J].草业学报,2009,18(5):
142149.
[10] 何东进,洪伟,胡海清.景观生态学的基本理论及中国景观生态学的研究进展[J].江西农业大学学报,2003,25(2):276
282.
[11] WuJ.Effectsofchangingscaleonlandscapepatternanalysis:scalingrelations[J].LandscapeEcology,2004,19(2):125
138.
[12] 潘竟虎,刘普幸,赵军.黑河下游土地利用与景观格局时空特征分析[J].土壤,2008,40(2):306311.
[13] 张宏斌,杨桂,霞黄青,等.呼伦贝尔草甸草原景观格局时空演变分析———以海拉尔及周边地区为例[J].草业学报,2009,
18(1):134143.
[14] 毕晓丽,周睿,刘丽娟,等.泾河沿岸景观格局梯度变化及驱动力分析[J].生态学报,2009,25(5):10411047.
[15] ChoppingM,RangoA,HavstadK,犲狋犪犾.Canopyattributesofdesertgrasslandandtransitioncommunitiesderivedfrom
multiangularairborneimagery[J].RemoteSensingofEnvironment,2003,85(3):339354.
[16] 熊小刚,韩兴国.运用状态与过渡模式讨论锡林河流域典型草原的灌丛化[J].草业学报,2006,15(2):913.
[17] SadlerR,HazeltonM,BoerM,犲狋犪犾.Derivingstateandtransitionmodelsfromanimageseriesofgrasslandpatterndynam
ics[J].EcologicalModeling,2010,221(3):433444.
[18] LangfordW,GergelS,DietterichT,犲狋犪犾.Mapmisclassificationcancauselargeerrorsinlandscapepatternindices:Exam
plesfromhabitatfragmentation[J].Ecosystems,2006,9(3):474488.
[19] 钱育蓉,杨峰,李建龙,等.利用高光谱数据快速估算高羊茅牧草光合色素的研究[J].草业学报,2009,18(4):94102.
[20] 张风丽,尹球,匡定波,等.草地光谱分类最佳时相选择分析[J].遥感学报,2006,10(4):482488.
[21] 周磊,辛晓平,张宏斌,等.高光谱遥感在草原监测中的应用[J].草业科学,2007,26(4):2027.
[22] IrisarriJ,OesterheldM,VeronS,犲狋犪犾.Grassspeciesdifferentiationthroughcanopyhyperspectralreflectance[J].Interna
tionalJournalofRemoteSensing,2009,30(22):59595975.
[23] 李文龙,薛中正,郭述茂,等.基于3S技术的玛曲县草地植被覆盖度变化及其驱动力[J].兰州大学学报(自然科学版),
2010,46(1):8591.
[24] 周华坤,周立,刘伟,等.青海省玛多县草地退化原因及畜牧业可持续发展[J].中国草地,2003,25(6):6267.
[25] 张镱锂,刘林山,摆万奇,等.黄河源地区草地退化空间特征[J].地理学报,2006,61(1):314.
[26] 薛娴,郭坚,张芳,等.高寒草甸地区沙漠化发展过程及成因分析———以黄河源区玛多县为例[J].中国沙漠,2007,27(5):
725732.
[27] 马智宏,王北洪,朱大洲,等.自主研究发明的成像光谱仪近地应用及建模初探[J].现代科学仪器,2010,(2):101104.
94第21卷第2期 草业学报2012年
[28] 张浩,张兵,卫征.基于POS的高光谱影像快速几何纠正[J].测绘通报,2009,1:1417.
[29] KruseF,LefkoffJ,BoardmanK,犲狋犪犾.Thespectralimageprocessingsystem(SIPS)interactivevisualizationandanalysisof
imagingspectrometerdata[J].RemoteSensingofEnvironment,1993,44:145163.
[30] SilvestriS,MaraniM,MaraniA.Hyperspectralremotesensingofsaltmarshvegetation,morphologyandsoiltopography[J].Phys
icsandChemistryoftheEarth,2003,28:1525.
[31] 童庆禧,张兵,郑兰芬.高光谱遥感———原理技术与应用[M].北京:高等教育出版社,2006.
[32] QiJ,MarsettR,MoranM,犲狋犪犾.SpatialandtemporaldynamicsofvegetationintheSanPedroRiverbasinarea[J].Agricultural
andForestMeteorology,2000,105(13):5568.
[33] 王军德,王根,陈玲.高寒草甸土壤水分的影响因子及其空间变异研究[J].冰川冻土,2006,28(3):428433.
[34] 李以康,林丽,张法伟,等.嵩草草甸退化和恢复过程中主要牧草演替和地表特征变化[J].草业学报,2010,19(5):179
185.
犔犪狀犱狊犮犪狆犲狆犪狋狋犲狉狀犪狀犪犾狔狊犻狊狅犳犪犾狆犻狀犲狊狋犲狆狆犲犫犪狊犲犱狅狀犪犻狉犫狅狉狀犲犺狔狆犲狉狊狆犲犮狋狉犪犾犻犿犪犵犲狉狔
犻狀犕犪犱狌狅犮狅狌狀狋狔,犙犻狀犵犺犪犻狆狉狅狏犻狀犮犲
JIAOQuanjun1,2,ZHANGBing1,ZHAOJingjing1,2,LIULiangyun1,2,HUYong1,2,3
(1.CenterforEarthObservationandDigitalEarth,ChineseAcademyofSciences,Beijing100094,China;
2.KeyLaboratoryofDigitalEarth,ChineseAcademyofSciences,Beijing100094,China;
3.ColegeofGeomatics,Xi’anUniversityofScience&
Technology,Xi’an710054,China)
犃犫狊狋狉犪犮狋:Landscapeecologyhasbeenwidelyusedingrasslandresourcemanagementandinecologicalsecurity
evaluation.Amultiscaleanalysisoflandscapestatisticsisveryimportantforcorrelationwithecologicalproces
ses.Acquiringfeaturesoflandscapepatternsfromdifferentgrasslandtypeswithvariousfractionalvegetation
coverisonetypeofbasicworkinlandscapescaleanalysisofgrasslanddegradation.Extractionofsmalscale
landscapepatternsfromhighresolutionremotesensingisausefulattempttocarryoutlandscapescaleanalysis.
Thepresentstudyanalyzeddifferentspectralcharacteristicsofdifferentgrasslandtypesandtookprecision
grasslandtypesmappingbasedonairbornehighspatialresolutionhyperspectralPISimageryinMaduocounty,
Qinghaiprovince.Metreresolutionlandscapepatterndifferencesinalpinesteppestudyareaswithdifferent
fractionalvegetationcoverwereexploredthroughthreelandscapeindices(landscapeheterogeneity,landscape
fragmentationandlandscapedominancelevel).Theresultsshowedthatlandscapeheterogeneityindexand
landscapedominanceindexarecloselyrelatedtotheleveloffractionalvegetationcoverinalpinegrassland.
However,withreducedfractionalvegetationcoverofalpinesteppesampleareas,fragmentationindexinitialy
decreasedbutthenincreased.
犓犲狔狑狅狉犱狊:sourceregionofYelowRiver;landscapepattern;alpinesteppe;airbornehyperspectral
05 ACTAPRATACULTURAESINICA(2012) Vol.21,No.2