免费文献传递   相关文献

Research on validating MODIS/LAI product in a typical semi-arid steppe of Xilinhot in Northern China

基于野外观测和TM数据的锡林浩特典型草原MODIS/LAI产品验证



全 文 :书基于野外观测和犜犕数据的锡林浩特典型
草原 犕犗犇犐犛/犔犃犐产品验证
王圆圆1,王猛2,3,李贵才1,王军邦2,杨忠东1,戎志国1
(1.中国气象局中国遥感卫星辐射测量和定标重点开放实验室 国家卫星气象中心,北京100081;
2.中国科学院地理科学与资源研究所,北京100101;3.中国科学院研究生院,北京100049)
摘要:通过2008年夏季生长季的3期LAI野外观测和TM 图像的处理分析,对锡林浩特国家观象台附近24km2
典型草原区的 MODIS/LAI产品精度从时空2种维度进行了评估。结果显示,在时间维度上,MODIS/LAI产品可
以较为准确的体现草地植被生长轨迹;在空间维度上,MODIS/LAI产品与TM/LAI反演结果有相似的变化趋势。
以TM/LAI结果为参考真值,MODIS/LAI产品的相对均方根误差约为24%,并存在10.8%的正偏差,表明有一定
程度的高估。误差分析结果显示,MODIS/LAI产品精度主要受到土壤背景光谱以及地表均匀程度的影响,在LAI
较小的区域,以及在草地空间异质性因受人类活动而增强的7月初和8月末,MODIS/LAI产品误差更大。
关键词:MODIS/LAI产品;TM/ETM+数据;典型草原
中图分类号:S812.3  文献标识码:A  文章编号:10045759(2011)04025209
  叶面积影响着植物光合、土壤呼吸、地表蒸散等生物物理过程,是重要的植被结构参数[13]。遥感具有大范
围、准实时、多尺度获取陆表参数的优势[4]。中分辨率成像光谱仪(moderateresolutionimagingspectroradiome
ter,MODIS)叶面积指数(leafareaindex,LAI)产品自发布以来就受到了多学科研究人员的密切关注,并被作为
关键参数输入到遥感、植被、生态、气候等领域的模型或模式中[5]。
MODIS/LAI产品的应用效果受制于产品质量,因此产品精度的独立验证非常重要;同时,产品精度验证也
将积极促进反演算法的改进[6]。MODIS/LAI产品验证的难点在于无法直接从地表获取同等像元尺度(1km×1
km)上的观测,因此高空间分辨率数据常被作为中间桥梁,将地表的单点观测转换至区域结果。验证方法通常包
括3个步骤:1)野外观测获取LAI;2)利用高空间分辨率遥感图像反演LAI;3)将高空间分辨率LAI尺度转换至
低空间分辨率,和MODIS/LAI产品进行对比。截至到目前,已有不少国内外学者都做过MODIS/LAI产品验证
工作,但由于研究区域不同,这些验证工作得出的结论也不尽相同。Cohen等[7,8]曾利用ETM+(enhancedthe
maticmapper)图像和地表实测LAI值,验证了北美4种不同植被覆盖区域(农田、草原、针叶林、混交林)的 MO
DIS/LAI产品,发现与旧版算法(Colection3)相比,新版算法(Colection4)可明显改善农田和草原植被覆盖区
的LAI估算精度,但是仍然高估了森林地区的LAI;Cohen等[7]在文章中还探讨了不同空间尺度下LAI与 MO
DIS反射率(波段1~7)之间的相关性,提出在反演算法中增加短波红外信息,将有助于提高LAI反演精度。
Pisek和Chen[9]也针对北美地区 MODIS/LAI产品开展了验证研究,得出 MODIS/LAI产品的相对绝对误差中
值变化范围为34%~88%,并发现当LAI处于高值区时,产品值的波动较大;他们认为增加亚像元地表覆盖信
息,或森林植被的聚集指数信息,或林下的背景光谱信息可以改善LAI产品的反演精度;Abuelgasim等[10]在加
拿大森林覆盖地区,验证了3种LAI产品[分别基于 MODIS,SPOT4/VGT(satelitepourl’observationdela
terre,VEGETATIONproduct)和POLDER1(POLarizationanddirectionalityoftheearth’sreflectances)遥感
数据生成],结果显示VGT/LAI产品与实测结果更为接近,另2种LAI产品精度较低,他们提出为改善产品质
量,很有必要在LAI反演算法中增加区域调节因素。Fensholt等[11]针对塞内加尔地区的Savanna草原进行了
252-260
2011年8月
   草 业 学 报   
   ACTAPRATACULTURAESINICA   
第20卷 第4期
Vol.20,No.4
 收稿日期:20100611;改回日期:20100731
基金项目:公益性行业科研专项[GYHY(QX)200769,GYHY200906036]和国家自然科学基金(40501048)资助。
作者简介:王圆圆(1981),女,安徽六安人,助研,博士。Email:wangyuany@cma.gov.cn
通讯作者。Email:ligc@cma.gov.cn
MODIS/LAI产品验证,认为产品可以较好体现植被的季相变化,但存在高估,幅度为2%~15%。从目前已经完
成的 MODIS/LAI产品验证工作来看,学者们较为关注浓密植被和生长旺季[79],因为在上述条件下,基于辐射传
输理论反演LAI的主程序经常无法获取有效值,继而求助基于归一化植被指数 NDVI(normalizeddifference
vegetationindex)的后备程序,但NDVI在LAI值较高的情况下存在饱和,使反演结果具有很大不确定性。饱和
现象固然对LAI反演有重要影响,土壤背景也是一个不容忽视的关键因素[11,12],尤其是在植被覆盖较为稀疏的
地区,可是当前针对这类区域的MODIS/LAI产品验证工作还比较少;此外,在已有很多研究中,验证多是基于区
域LAI均值的时间序列展开的,较少有空间格局的对比。
本研究以冠层开放、覆盖稀疏的内蒙古锡林浩特典型草原区为研究对象,通过在2008年夏季开展的3次野
外观测试验,获取大量样方LAI,与TM(thematicmapper)/ETM+光谱数据建立相关,生成高精度的LAI分布
图,从时间和空间2个维度上对MODIS/LAI产品进行验证,并探讨了其误差分布特征。研究结果对于了解该地
区的 MODIS/LAI产品质量具有重要意义,对前人LAI验证研究亦是很好的补充。
1 材料与方法
1.1 研究区和样方布设
研究区位于内蒙古自治区锡林浩特市东北部25km处的锡林浩特国家气候观象台野外站,面积大约24
km2,海拔约1000m(图1)。该区为半干旱气候,冬季受蒙古高压气流控制,寒冷干燥,夏季受季风影响,较为温
暖湿润,年均气温0.5~1.0℃,无霜期约100d,年均降水350mm,集中在6-9月[13]。从植被类型上来看,该区
是亚欧大陆的典型草原区,建群种以克氏针茅(犛狋犻狆犪犽狉狔犾狅狏犻犻)、羊草(犔犲狔犿狌狊犮犺犻狀犲狀狊犻狊)为主。土壤为栗钙土。
图1 犜犕假彩色合成图像(2008年8月25日)上研究区80个样方(黄点)示意图
犉犻犵.1 犐犾狌狊狋狉犪狋犻狅狀狅犳80狆犾狅狋狊(狔犲犾狅狑狆狅犻狀狋狊)狅犳狋犺犲狊狋狌犱狔犪狉犲犪狅狀犜犕犳犪犾狊犲犮狅犿狆狅狊犻狋犲犻犿犪犵犲(犪犮狇狌犻狉犲犱狅狀犃狌犵狌狊狋25,2008)
 白色格网表示等经纬度投影下 MODIS/LAI产品像元Thewhitegrids(1km×1km)arethepixelsofMODIS/LAIproduct
研究区中立有CO2 涡度塔(N44°08′31.2″,E116°18′42.5″),海拔1158m。按照循环采样的模式[14],在涡
度塔周围布设了80个25m×25m样方(如图1中黄点)。循环采样是一种非常高效的采样方式,在降低了采样
频率的条件下,仍然能够在不同采样距离上保持采样密度,避免了在相同采样距离上的多余采样。具体的采样方
案见表1,在每一个循环长度上,最优的采样位置是能够实现以最少的取样数目,提供足够量的相差不同距离的
352第20卷第4期 草业学报2011年
数据对数。本研究以涡度塔为中心,沿东西和南北方
向分别采用5/21和4/13的采样设计方案,采样的基
本步长是25m,循环长度分别是21和13基本步长
处,分别在0,1,4,14,16和0,1,3,9基本步长的倍数
处设置采样点(图2)。在通量贡献区约950m×550
m的区域上,共设置80个采样点(受围栏和道路的影
响,实际有效样方数为78个)。
除了上述80个密集样方,另有20个25m×25m
样方随机分布在涡度塔周围5km×5km范围内,以
准确获取整个研究区LAI的数据变异特征。样方通
过经纬仪和激光测距仪精确定位,在样方的4个角点
设置木桩作为标志。每个样方内分别在4个角和中心
位置设置5个0.5m×0.5m的亚样方。计算LAI
时,将5个亚样方的平均值作为样方的值。
表1 循环采样方案的选点方式[14]
犜犪犫犾犲1 犇犻犳犳犲狉犲狀狋狊犪犿狆犾犻狀犵狊犮犺犲犿犲狊犫犪狊犲犱
狅狀犮狔犮犾犻犮狊犪犿狆犾犻狀犵犱犲狊犻犵狀
循环方式
Cycle
definition
循环长度
Cycle
length
采样点数
Samplingpoint
number
采样位置
Sampling
position
1/1 1 1 0
2/3 3 2 0,1
3/7 7 3 0,1,3
4/13 13 4 0,1,3,9
5/21 21 5 0,1,4,14,16
6/31 31 6 0,1,3,8,12,18
7/37 37 7 0,1,6,10,17,23,35
图2 基于循环采样的80个样方的空间分布格局
犉犻犵.2 犔犪狔狅狌狋狅犳狋犺犲80狆犾狅狋狊犻狀狊狋狌犱狔狉犲犵犻狅狀犫犪狊犲犱狅狀犮狔犮犾犻犮狊犪犿狆犾犻狀犵犱犲狊犻犵狀
  亚样方内的LAI通过称重收割法进行测量。具体做法是,用0.5m×0.5m的样方框,把植物按针茅、羊草、
知母(犃狀犲犿犪狉狉犺犲狀犪犪狊狆犺狅犱犲犾狅犻犱犲狊)、沙葱(犃犾犾犻狌犿犿狅狀犵狅犾犻犮狌犿)、杂类草分别齐根剪下带回实验室称量鲜重,
70℃下烘48h,获取样方内每个物种的干重,物种干重与比叶面积(SLA,specificleafarea)相乘即可得到物种叶
面积指数,样方内所有物种的叶面积指数相加可得样方的总叶面积指数。为测量SLA,每物种取10个叶片样
本,平铺到扫描仪上进行扫描,扫描时使羊草、知母等宽叶植物尽量展开,针茅、葱等针状的或者卷曲的植物,直接
用投影面积进行测量。扫描图片用 Matlab软件进行计算,得到叶片面积,叶片面积除以叶片干重,即可获取每个
物种的SLA。
452 ACTAPRATACULTURAESINICA(2011) Vol.20,No.4
1.2 数据及研究方法
2008年夏季,国家卫星气象中心组织相关工作人员进行了3次野外观测试验,时间分别为7月初(7月2-7
日),7月末(7月20-26日)和8月末(8月25-30日)。7月初,草场处于快速生长早期,7月末,草场属于生长
最旺盛时期,到8月末,针茅开始出现抽穗,部分草场已开始收割。
根据野外观测时间,购买了2幅TM影像(获取时间分别为2008年7月8日和2008年8月25日),下载1
幅ETM+影像(获取时间为2008年8月1日,http://glovis.usgs.gov)。总计3幅影像,图像轨道号为124/29,
在研究区域上空均是无云晴空(ETM+图像上则不存在数据丢失线)。首先对3幅影像进行了几何精纠正
(RMSE<0.5像元);其次,根据头文件信息做了辐射定标,获得表观反射比;最后,是大气纠正工作,由于2008
年7月8日获取了CE318观测数据和探空资料,可知气溶胶光学厚度为0.30,大气水汽含量为2.018g/cm2,因
此利用6s模型(secondsimulationofasatelitesignalinthesolarspectrum)对该天TM影像作了大气校正,通
过向6s模型输入图像获取时间、太阳天顶角、太阳方位角、气溶胶光学厚度、大气水汽含量、表观反射比、地表高
程等参数,即可获得除去大气吸收、散射等影响的真实地表反射率。鉴于在8月1日和8月25日都没有大气探
测资料,同时为了保证多时相图像数据在光谱上的一致性,该2日内的图像利用稳定目标进行了相对的大气校
正[15],选择的稳定目标包括具有低反射特性的水体和高反射特性的锡林浩特机场跑道,假设这2种地物的反射
率在短时间内不发生变化,建立表观反射比和地表反射率的线性回归方程,即可反推其他地物的地表反射率。
图3 犕犗犇犐犛/犔犃犐产品验证流程
犉犻犵.3 犉犾狅狑犮犺犪狉狋狅犳犕犗犇犐犛/犔犃犐狆狉狅犱狌犮狋狏犪犾犻犱犪狋犻狅狀
 a:基于密集采样区的 MODIS/LAI产品验证 MODIS/LAIvalidationthroughcomparisonwithfieldmeasurements;b:基于 TM/LAI反演结果的
MODIS/LAI产品验证 MODIS/LAIvalidationthroughcomparisonwithTM/LAIinversionresults
552第20卷第4期 草业学报2011年
  针对 MODIS/LAI产品(8d合成、1km空间分辨率),下载了10期数据,时间跨度是从第169天至第241
天,包括了整个夏季生长季。
本研究分2种途径验证 MODIS/LAI产品,一个是时间序列的分析,直接对比80个密集采样区的样方LAI
观测值与 MODIS/LAI产品(图3a),这样做可以直观看到2种不同观测尺度对研究区LAI的刻画特征,虽然从
空间上来看有不妥之处,但因为地表采样数量较多,对时间序列比较的影响不大;另一个是空间格局的分析,以1
kmTM/LAI反演结果为参照,和研究区内 MODIS/LAI产品进行逐像元的对比(图3b),虽然有不少学者建议,
考虑到定位的不确定性,应该在图像分割后对比目标之间的LAI[16],由于研究区为均匀草地,不存在其他土地覆
盖类型,对定位错误的敏感性较低,因此采取逐像元对比是可行的。
2 结果与分析
2.1 密集采样区对比分析结果
经过检查,发现 MODIS/LAI产品质量控制参数(QC)在研究区内均为0,说明LAI反演值的可信度很高。
密集采样区在产品图像上跨越了4个相邻像元(图1),其平均值和标准差见图4中蓝线所示,从中可以看出,
MODIS/LAI产品时间序列可较好反映草地植被在生长季由稀疏到茂密的生长过程;从6月下旬-7月初,LAI
迅速增加,由0.4上升到1以上,到8月下旬,LAI回落至0.7左右。
地表实测LAI(图4中红线所示)同样也表现出开始迅速增加而后缓慢下降的趋势,其和 MODIS/LAI均值
的吻合程度在7月初最好,在7月末和8月末,偏离程度较大,而且偏离的方向也不一致,7月末,MODIS/LAI产
品值偏小,8月末,产品值则偏大。地表实测的LAI还具有标准差高、变异系数大的特点,显示出样方尺度的LAI
观测包含更多细节信息(表2),7月初地表实测LAI变异系数明显高于其他2个时期。
表2 密集采样区实测犔犃犐和 犕犗犇犐犛/犔犃犐比较
犜犪犫犾犲2 犆狅犿狆犪狉犻狊狅狀犫犲狋狑犲犲狀狅犫狊犲狉狏犲犱犔犃犐犪狀犱犕犗犇犐犛/犔犃犐狆狉狅犱狌犮狋犪狋犻狀狋犲狀狊犻狏犲狊犪犿狆犾犲犪狉犲犪
生长时期
Growthstage
实测LAI(狀=78)
ObservedLAI(狀=78)
平均值 Mean 标准差Std 变异系数CV
MODIS/LAI(狀=4)
平均值 Mean 标准差Std 变异系数CV
7月初(7月2-7日)EarlyJuly(July2-7) 0.6800 0.2632 0.3871 0.6500 0.0577 0.0888
7月末(7月20-26日)LateJuly(July20-26) 0.9959 0.2466 0.2476 0.8750 0.0957 0.1094
8月末(8月25-30日)LateAugust(August25-30) 0.8849 0.1950 0.2203 0.9750 0.0500 0.0513
 狀为样本数狀isthesamplesize.
2.2 基于TM/LAI反演结果的 MODIS/LAI产品验证
通过对样方LAI和TM各波段光谱以及植被指数进行相关分析后,发现LAI与位于短波红外的第5波段
(中心波长1.6μm)相关性最好,但是相关系数的绝对值仍然不是很高,只有0.5。如果采用最小二乘法则会产
生方差压缩的问题,为了保持LAI观测值的变异程度,采取了RMA(reducedmajoraxis)回归方法[17]
犢=(α+β犡+ε),α=珚犢-(σ犢/σ狓)珡犡,β=σ犢/σ狓
式中,犢 和犡 为因变量和自变量,α、β为系数,珚犢和珡犡 分别是犢 和犡 的均值,σ犢 和σ狓 是犢 和犡 的标准差。
图5显示了TM第5波段和LAI的散点图以及回归式。利用该回归式可以得到3期的TM/LAI分布图,将
其按照 MODIS/LAI产品像元网格(等经纬度投影)进行平均,获取1km空间分辨率的TM/LAI结果,即可与
MODIS/LAI产品进行逐像元的对比。
图6显示了 MODIS/LAI和1km空间分辨率TM/LAI的散点图及线性拟合方程,可看出两者变化趋势一
致,具有较好的相关性(狉=0.6),但是 MODIS/LAI产品数值反差较小,值域为[0.4,1.1],而TM/LAI的数值值
域为[0.2,1.3]。从回归结果来看,LAI较小的样本更偏离拟合线,说明 MODIS/LAI的反演误差在草地覆盖稀
疏区会更大。表3中的误差统计结果也说明了这一点,当LAI>0.85时,相对均方根误差为19.69%,而当
LAI≤0.85时,相对均方根误差达到44.36%。
652 ACTAPRATACULTURAESINICA(2011) Vol.20,No.4
从不同时期的MODIS/LAI产品来看,7月初和8
图4 犕犗犇犐犛/犔犃犐产品时间序列(蓝色)和
地面实测犔犃犐(红色)对比
犉犻犵.4 犆狅犿狆犪狉犻狊狅狀犫犲狋狑犲犲狀犕犗犇犐犛/犔犃犐狋犻犿犲
狊犲狉犻犲狊(犫犾狌犲)狑犻狋犺犳犻犲犾犱犿犲犪狊狌狉犲犿犲狀狋狊(狉犲犱)
月末的LAI产品精度低于7月末(表4),这与草场的
均匀性有关。7月初是草地生长早期,上一年草场利
用活动对这一时期的草地格局仍有重要影响;这一时
期围封草地LAI明显高于刈割区,草地有很强的空间
异质性(见表2中变异系数)。到7月末,上一年被刈
割的草场完全长出新草,自然因素尤其是地形及其对
降水的再分配,对空间格局的形成起到主导作用。由
于研究区地形平缓,草地生长比较均匀,加之草地生长
茂盛,LAI值大,使得7月末的MODIS/LAI产品精度
最高。到了8月末,又值草场收割时间,从 TM 影像
上可明显看出草地分布呈斑块状(图1),因为是拖拉
机收割作业,速度很快,草场的时空格局变化剧烈。而
MODIS/LAI是8d合成产品,因此用一期TM 影像
评估 MODIS/LAI产品精度存在较大局限,也容易得
出产品精度较低的结论。
图5 基于犜犕波段5的犔犃犐估算回归式
犉犻犵.5 犛犮犪狋狋犲狉狆犾狅狋犪狀犱狉犲犵狉犲狊狊犻狅狀犾犻狀犲犫犲狋狑犲犲狀犜犕犫犪狀犱
5(1.6μ犿)犪狀犱犳犻犲犾犱犔犃犐犿犲犪狊狌狉犲犿犲狀狋狊
图6 犕犗犇犐犛/犔犃犐与犜犕/犔犃犐的散点图
犉犻犵.6 犛犮犪狋狋犲狉狆犾狅狋犫犲狋狑犲犲狀犕犗犇犐犛/犔犃犐
犪狀犱犜犕/犔犃犐
 红点、绿点、蓝点分别表示7月初(7月2-7日),7月末(7月20-26日)和8月末(8月25-30日)的观测数据 Red,greenandbluepointsrepre
sentobservationsfromearlyJuly(July2-7),lateJuly(July20-26)andlateAugust(August25-30)respectively
3 讨论
3.1 TM/LAI估算模型
在TM光谱和LAI的相关分析中,发现位于短波红外的第5波段(1.6μm)与LAI的相关系数最大,这在其
他文献中也有提及[7,18]。在植被冠层覆盖率低、地表背景对遥感观测信号有较大影响的情况下,近红外和NDVI
等参数与LAI几乎无关,而短波红外则表现出很好的预测能力[18,19]。之所以可以用于估算LAI是因为短波红
外对植被体内的液态水含量非常敏感,对于正常生长的多数植被而言,植被水分都是和LAI紧密相关的[12,20],这
一结论对于今后利用遥感手段监测典型草地生物物理参数具有重要指导作用。
752第20卷第4期 草业学报2011年
表3 研究区 犕犗犇犐犛/犔犃犐产品针对不同类型像元的误差统计(以犜犕/犔犃犐结果为真值)
犜犪犫犾犲3 犈狉狉狅狉狊狋犪狋犻狊狋犻犮狊狅犳犕犗犇犐犛/犔犃犐狆狉狅犱狌犮狋犳狅狉犱犻犳犳犲狉犲狀狋狋狔狆犲狊狅犳狆犻狓犲犾狊狌狊犻狀犵犜犕/犔犃犐犪狊狋犺犲狉犲犳犲狉犲狀犮犲
误差
Error
所有像元(狀=72)
Althepixels
LAI值≤0.85的像元(狀=48)
OnlypixelswhoseLAIvalues
≤0.85areconsidered
LAI值>0.85的像元(狀=24)
OnlypixelswhoseLAIvalues
>0.85areconsidered
均方根误差Rootmeansquareerror(RMSE) 0.2388 0.2527 0.2082
相对均方根误差RelativeRMSE(%) 32.61 44.36 19.69
 狀为样本数狀isthesamplesize.
表4 研究区不同时期 犕犗犇犐犛/犔犃犐产品误差统计(以犜犕/犔犃犐结果为真值)
犜犪犫犾犲4 犈狉狉狅狉狊狋犪狋犻狊狋犻犮狊狅犳犕犗犇犐犛/犔犃犐狆狉狅犱狌犮狋犪狋犱犻犳犳犲狉犲狀狋犵狉狅狑狋犺狊狋犪犵犲狊狌狊犻狀犵犜犕/犔犃犐犪狊狋犺犲狉犲犳犲狉犲狀犮犲
误差
Error
7月初(7月2-7日)
EarlyJuly(July2-7)
7月末(7月20-26日)
LateJuly(July20-26)
8月末(8月25-30日)
LateAugust(August25-30)
均方根误差Rootmeansquareerror(RMSE) 0.1979 0.2084 0.2975
相对均方根误差RelativeRMSE(%) 35.89 21.18 44.95
3.2 MODIS/LAI产品精度评价
以1km分辨率TM/LAI作为参照(即相对更准确的结果),评价MODIS/LAI产品时发现,两者在空间上的
变化趋势比较一致。但是当LAI较小时,MODIS/LAI产品误差偏大,说明地表背景光谱干扰对产品精度有明
显影响。从均值来看,MODIS/LAI产品的均值为0.81,TM/LAI的均值为0.73,说明产品值存在10.8%的正
偏差。
3.3 与同类研究的比较
Fensholt等[11]得出 MODIS/LAI产品在草原高估2%~15%,本研究得出的10.8%高估水平与此具有一致
性。Cohen等[8]提出,在叶片生物量较小的地方,MODIS/LAI产品和地表实测结果的一致性更好,本研究结果
也印证了这一观点,这一方面可能是由产品算法自身特点决定的,另一方面可能与2008年锡林浩特地区雨水较
多、草地长势良好、土壤裸露度较低有关。今后如果能在干旱年份也开展相关地表观测试验,将有助于更加全面
认识 MODIS/LAI产品质量。
3.4 MODIS/LAI产品误差分析
影响MODIS/LAI产品质量的主要因素包括:1)用于建立查找地表的物理模型的不确定性;2)土地覆盖类型
的不确定性;3)地表反射率的不确定性。区域性的影响因素是后2个,考虑到研究区土地覆盖类型几乎全是草
地,对本研究结果影响最大的是地表反射率的不确定性。在7月初和8月末,草地异质性更高,地表反射率不确
定性增大,故导致了 MODIS/LAI产品精度降低。
3.5 关于研究区范围的选择
本研究野外试验的观测样方80%分布于涡度塔周围1km2 范围内,另有20%分布于涡度塔周围25km2 范
围内,考虑到这种空间分布特征,选定了涡度塔周围20~30km2 作为研究区。虽然面积较小,但在该区域内,样
本观测频率在空间和时间上都很密集,足够的样本量确保了研究结论具有很高的可信度;此外,由于锡林浩特地
区的植被类型绝大部分为典型草原,结论在该区域也具有一定代表性。为了在更大空间范围上准确评估锡林浩
特地区的LAI产品质量,需要在以后的野外试验中扩大采样点分布范围(不仅仅是涡度塔周围)、增加植被观测
类型(如在该地区东部,存在其他类型草地、灌木、森林等)。就目前所获得的试验数据而言,较小的验证区域较为
合理。
4 结论
本研究通过2008年夏季在锡林浩特典型草原区开展的3期LAI观测试验,辅以TM/LAI反演、尺度转换以
852 ACTAPRATACULTURAESINICA(2011) Vol.20,No.4
及空间匹配,实现了从时空2个维度上验证了 MODIS/LAI产品的数据质量。研究结果显示,尽管在典型草原
区,植被冠层开放,LAI值偏小(LAI<2),但 MODIS/LAI产品质量仍较令人满意:在时间维度上,可准确体现植
被的生长轨迹,在空间维度上,一定程度反映了地表变化的细节。误差分析结果显示,MODIS/LAI产品相对均
方根误差约为24%,在LAI值较小或是草地均匀性较差的情况下,产品误差增加,表明土壤背景以及地表均匀程
度是该区域 MODIS/LAI产品精度的主要影响因素,而这2个因素又与当地人类活动紧密相关。
致谢:王英舜、国振杰、师桂花、骆俊、陈洁、郑照军、毕研盟、陆其锋、武胜利、张勇、李元等同志参加了部分野外观
测,在此表示感谢!
参考文献:
[1] BerterretcheM,HudakAT,CohenWB.Comparisonofregressionandgeostatisticalmethodsformappingleafareaindex
(LAI)withLandsatETM+dataoveraborealforest[J].RemoteSensingofEnvironment,2005,96:4961.
[2] 周萍,刘国彬,薛?.草地生态系统土壤呼吸及其影响因素研究进展[J].草业学报,2009,18(2):184193.
[3] 谢贤健,兰代萍,白景文.三种野生岩生草本植物的抗旱性综合评价[J].草业学报,2009,18(4):7580.
[4] 梁天刚,崔霞,冯琦胜,等.2001-2008年甘南牧区草地地上生物量与载畜量遥感动态监测[J].草业学报,2009,18(6):
1222.
[5] ReichP,TurnerD,BolstadP.Anapproachtospatialydistributedmodelingofnetprimaryproduction(NPP)attheland
scapescaleanditsapplicationinvalidationofEOSNPPproducts[J].RemoteSensingofEnvironment,1999,70:6981.
[6] YangWZ,TanB,HuangD,犲狋犪犾.MODISleafareaindexproducts:Fromvalidationtoalgorithmimprovement[J].IEEE
TransactionsonGeoscienceandRemoteSensing,2006,44(7):18851897.
[7] CohenWB,MajerspergerTK,YangZQ,犲狋犪犾.ComparsionoflandcoverandLAIestimatesderivedfromETM+andMO
DISforfoursitesinNorthAmerica:Aqualityassessmentof2000/2001provisionalMODISproducts[J].RemoteSensingof
Environment,2003,88:233255.
[8] CohenWB,MaierspergerTK,TurnerDP,犲狋犪犾.MODISlandcoverandLAIcolection4productqualityacrossninesitesin
thewesternHemisphere[J].IEEETransactionsonGeoscienceandRemoteSensing,2006,44:18431857.
[9] PisekJ,ChenJM.ComparisonandvalidationofMODISandVEGETATIONglobalLAIproductsoverfourBigfootsitesin
NorthAmerica[J].RemoteSensingofEnvironment,2007,109:8194.
[10] AbuelgasimRA,FernandeRA,LeblancSG.EvaluationofnationalandglobalLAIproductsderivedfromopticalremote
sensinginstrumentsoverCanada[J].IEEETransactionsonGeoscienceandRemoteSensing,2006,44:18721884.
[11] FensholtR,SandholtI,SchultzRasmussenM.EvaluationofMODIS/LAI,fAPARandtherelationbetweenfAPARand
NDVIinasemiaridenvironmentusinginsitumeasurements[J].RemoteSensingofEnvironment,2004,91:490507.
[12] HouborgR,SoegaardH,BoeghE.Combiningvegetationindexandmodelinversionmethodsfortheextractionofkeyvegeta
tionbiophysicalparametersusingTerraandAquaMODISreflectancedata[J].RemoteSensingofEnvironment,2007,106:
3958.
[13] 李素英,李晓兵,莺歌,等.基于植被指数的典型草原区生物量模型———以内蒙古锡林浩特市为例[J].植物生态学报,
2007,31:2331.
[14] BurrowsSN,GowerST,ClaytonMK,犲狋犪犾.Applicationofgeostatisticstocharacterizeleafareaindex(LAI)fromflux
towertolandscapescalesusingacyclicsamplingdesign[J].Ecosystems,2005,5:667679.
[15] SongC,WoodcockC,SetoK,犲狋犪犾.ClassificationandchangedetectionusingLandsatTMdata:whenandhowtocorrectat
mosphericeffects[J].RemoteSensingofEnvironment,2001,75:230244.
[16] TianY,WoodcockCE,WangY,犲狋犪犾.MultiscaleanalysisandvalidationoftheMODISLAIproductⅠ.Uncertaintyas
sessment[J].RemoteSensingofEnvironment,2002,83:414430.
[17] CohenWB,MaierspergerTK,GowerST,犲狋犪犾.AnimprovedstrategyforregressionofbiophysicalvariablesandLandsat
ETM+data[J].RemoteSensingofEnvironment,2003,84:561571.
952第20卷第4期 草业学报2011年
[18] BerterretcheM,HudakAT,CohenmWB,犲狋犪犾.ComparisonofregressionandgeostatisticalmethodsformappingLeafAr
eaIndex(LAI)withLandsatETM+dataoveraborealforest[J].RemoteSensingofEnvironment,2005,96:4961.
[19] TunerD,CohenW,KennedyR,犲狋犪犾.RelationshipsbetweenLAIandLandsatTMspectralvegetationindicesacrossthree
temperatezonesites[J].RemoteSensingofEnvironment,1999,70:5268.
[20] CeccatoP,FlasseS,TarantolaS,犲狋犪犾.Detectingvegetationleafwatercontentusingreflectanceintheopticaldomain[J].
RemoteSensingofEnvironment,2001,77:2223.
犚犲狊犲犪狉犮犺狅狀狏犪犾犻犱犪狋犻狀犵犕犗犇犐犛/犔犃犐狆狉狅犱狌犮狋犻狀犪狋狔狆犻犮犪犾狊犲犿犻犪狉犻犱狊狋犲狆狆犲狅犳犡犻犾犻狀犺狅狋犻狀犖狅狉狋犺犲狉狀犆犺犻狀犪
WANGYuanyuan1,WANGMeng2,3,LIGuicai1,WANGJunbang2,
YANGZhongdong1,RONGZhiguo1
(1.KeyLabofRadiometricCalibrationandValidationforEnvironmentalSatelites,ChinaMeteorological
Administration(LRCVES/CMA),NationalSateliteMeteorologicalCenter,Beijing100081,China;
2.InstituteofGeographicSciencesandNaturalResources,CAS,Beijing100101,China;
3.GraduateSchoolofChineseAcademyofSciences,Beijing100049,China)
犃犫狊狋狉犪犮狋:MODIS/LAIproductinatypicalsemiaridsteppeofXilinhotinNorthernChinawasvalidatedbased
onthreecomprehensivefieldsurveyswhichwerecarriedoutduringJuly-August2008inthestudyarea.LAI
wasmeasuredin25m×25mplotsbyadestructiveharvestmethodof80intensiveplotsand20scatteredplots.
ThreeTM/ETM+imagesacquiredduringthefieldsurveyperiodswerecolectedandprocessed.Thevalidation
workcoveredtwoparts:Timeseriesanalysisandspatialpatternanalysis.Intimeseriesanalysis,MODIS/LAI
trackintheintensivesamplingregionextractedfromDOY169toDOY241werecomparedwithgroundmeas
uredvalues.Inspatialpatternanalysis,LAImapwith30mspatialresolutionwerefirstderivedbybuildinga
regressionlinebetweenLAIandTMband5(1.6μm),andthenaggregatedto1kmspatialresolutionbasedon
MODISgrids.ResultsindicatedthatinthetemporaldomainMODIS/LAIproductcancorrectlyreflectgrass
growthtrends;InthespatialdomainitshowedasimilarpatternwithTM/LAI(1km).UsingtheTM/LAIas
reference,therelativeRMSEofMODIS/LAIwasabout24%andapositivebias(10.8%)existed,suggesting
slightoverestimation.ErrordistributionanalysisrevealedthattheMODIS/LAIproductwaslessaccuratein
regionswithlowerLAIvalues,probablyduetosoilcontamination.Moreover,therelativeRMSEoftheMO
DIS/LAIproductwashigheratthebeginningofJulyandtheendofAugustwhentheheterogeneityofgrass
coverincreasedasaresultofhumanactivities.
犓犲狔狑狅狉犱狊:MODIS/LAIproduct;TM/ETM+data;typicalsteppe
062 ACTAPRATACULTURAESINICA(2011) Vol.20,No.4