一种荒漠草原非生长季枯草重量遥感估算方法与流程

文档序号:11230947阅读:1096来源:国知局
一种荒漠草原非生长季枯草重量遥感估算方法与流程

本发明涉及一种遥感估算荒漠草原枯草重量的方法,尤其是一种利用modis卫星遥感数据快速准确估算荒漠草原非生长季节枯草重量格局与动态的方法。



背景技术:

枯草是荒漠草原非生长季植被重要的组成部分,枯草重量的准确估算对于荒漠草原生态系统碳储量的研究具有重要意义。荒漠草原非生长季枯草重量动态及模拟研究一直是陆地生态学的热点问题。枯草重量被广泛应用于开展荒漠草原植被生产力模型、干旱监测、火灾预警以及畜牧业规划等方面的实践与研究工作。

荒漠草原枯草重量获取的传统方法为实地测量,包括目测法、样点法、样方法以及样线法等,但这些方法一般只适用于小区域尺度,且费时费力、难以大面积展开及具有一定的破坏性。近年来,遥感技术迅猛发展,并以其宏观性、时效性、信息量大等特点在荒漠草原资源调查等领域展示了日渐广阔的应用前景。遥感技术的发展为荒漠草原非生长季枯草重量的获取提供了新的技术手段,尤其是为大范围荒漠草原非生长季枯草重量快速、准确估算提供了可能。

目前,基于卫星遥感技术估算大区域尺度草地非生长季枯草重量的方法主要是利用modis波段6与波段7反射率组合的ndti指数(波段6反射率与波段7反射率的归一化指数)与sti指数(波段6反射率与波段7反射率的比值)。不过,由于荒漠草原植被稀疏,土壤背景对modis波段6与波段7反射率的影响较大,ndti指数估算荒漠草原非生长季枯草重量的精度并不理想。为科学合理利用荒漠草原非生长季枯草资源,亟需提供一种准确估算荒漠草原非生长季枯草重量遥感估算方法。



技术实现要素:

本发明要解决的具体技术问题是如何在荒漠草原植被稀疏,土壤背景对modis波段6与波段7反射率的影响较大,ndti指数估算荒漠草原非生长季枯草重量不准确的问题,并构建估算荒漠草原非生长季枯草重量的遥感光谱指数,提供一种荒漠草原非生长季枯草重量遥感估算方法。

一种荒漠草原非生长季枯草重量遥感估算方法,所述枯草重量遥感估算方法是按照下列步骤进行的:

1)研究区域土地利用数据获取

获取研究区域的土地利用现状遥感监测数据,提取生成研究区域的荒漠草原空间分布图;

2)枯草样地选取

依据上述步骤1)描述方法得到的研究区域的荒漠草原空间分布图,在非生长季选取分布均匀的n个1km×1km的枯草样地,采用手持式全球定位系统记录枯草样地的经纬度坐标;

3)样地枯草重量测定

在每个样地内随机选取m个1m×1m的样方获取其枯草重量,并求取m个样方枯草重量的平均值作为该样地的枯草重量;

4)样地植被指数计算

下载与枯草重量测定日期同步的覆盖样地的mod09ga反射率影像,进行数据格式与地图投影转换;提取位于样地内的mod09ga像元的短波红外第6波段与短波红外第7波段反射率,计算如下光谱指数:

s-ndti-0.1=(1+0.1)×(b6-b7)/(b6+b7+0.1)

s-ndti-0.2=(1+0.2)×(b6-b7)/(b6+b7+0.2)

s-ndti-0.3=(1+0.3)×(b6-b7)/(b6+b7+0.3)

s-ndti-0.4=(1+0.4)×(b6-b7)/(b6+b7+0.4)

s-ndti-0.5=(1+0.5)×(b6-b7)/(b6+b7+0.5)

s-ndti-0.6=(1+0.6)×(b6-b7)/(b6+b7+0.6)

s-ndti-0.7=(1+0.7)×(b6-b7)/(b6+b7+0.7)

s-ndti-0.8=(1+0.8)×(b6-b7)/(b6+b7+0.8)

s-ndti-0.9=(1+0.9)×(b6-b7)/(b6+b7+0.9)

s-ndti-1.0=(1+1.0)×(b6-b7)/(b6+b7+1.0)

式中,b6为mod09ga的短波红外波段第6波段反射率,b7为mod09ga的短波红外波段第7波段反射率;

5)确定估算枯草重量的最佳光谱指数与遥感模型

采用一元线性分析对样地枯草重量与其对应的上述步骤4)中计算的10种光谱指数进行回归拟合,通过比较回归模型的决定系数r2与留一交叉验证均方根误差rmsecv,确定估算枯草重量的最佳植被指数及其对应的遥感模型;

6)计算研究区域待估算日期的步骤5)中确定的最佳光谱指数

获取研究区域待估算日期的mod09ga反射率遥感影像,进行数据格式与地图投影转换,依据上述步骤1)获得的研究区域荒漠草原空间分布图,逐像元计算荒漠草原分布区待估算日期的步骤5)中确定的短波红外第6波段与第7波段反射率组合的最佳光谱指数;

7)研究区域待估算日期枯草重量遥感估算

依据上述步骤6)计算的研究区域荒漠草原分布区待估算日期的最佳光谱指数,并通过上述步骤5)确定的最佳光谱指数对应枯草重量遥感估算模型,逐像元计算研究区域荒漠草原分布区待估算日期的枯草重量;

8)研究区域非生长季枯草重量遥感动态估算研究

利用上述步骤1)-7)的荒漠草原非生长季枯草重量遥感估算方法,对研究区域非生长季枯草重量进行一天、一周或一月的变化过程监测,研究荒漠草原非生长季枯草重量的动态变化规律,对荒漠草原生态系统研究以及可持续利用发展提供科学支撑。

上述技术方案的实施,消除了土壤背景对荒漠草原非生长季枯草重量遥感估算的影响,构建了监测荒漠草原非生长季枯草重量的光谱指数及估算模型,并应用此模型对研究区域的荒漠草原非生长季枯草重量进行了遥感反演,克服了传统的实地调查方法的局限性以及现有的遥感光谱指数估算精度低的问题。本方法弥补了基于modis遥感数据估算荒漠草原非生长季枯草重量的空白,为及时准确估算大区域尺度范围荒漠草原非生长季枯草重量提供了一种新方法,实现了基于modis遥感数据的荒漠草原非生长季大范围区域枯草重量的快速准确监测,为荒漠草原生态研究以及可持续利用提供科学依据。

附图说明

图1是本发明实施案例所述的一种荒漠草原非生长季枯草重量遥感估算方法中不同光谱指数估算枯草重量的精度比较图。

图2是本发明实施案例所述的一种荒漠草原非生长季枯草重量遥感估算方法中最佳光谱指数与实测枯草重量相关关系图。

图3是本发明实施案例所述的一种荒漠草原非生长季枯草重量遥感估算方法中苏尼特荒漠草原枯草重量空间分布图。

图4是本发明实施案例所述的一种荒漠草原非生长季枯草重量遥感估算方法中对比例。

具体实施方式

下面通过一种苏尼特荒漠草原非生长季枯草重量遥感估算方法的应用案例,对本发明的具体实施方式作进一步说明。

实施一种荒漠草原非生长季枯草重量遥感估算方法是按下列步骤进行的:

步骤一、苏尼特荒漠草原土地利用数据获取

获取苏尼特(苏尼特左旗、苏尼特右旗)荒漠草原遥感影像产品mcd12q1,使用mrt软件将mcd12q1数据进行格式与地图投影转换,hdf格式转换为tiff格式,正弦曲线地图投影转换为wgs/albers地图投影,生成苏尼特荒漠草原空间分布图;

步骤二、苏尼特荒漠草原枯草样地选取

依据上述步骤一描述方法得到的苏尼特荒漠草原空间分布图,在非生长季10月中旬选取分布均匀的30个枯草样地(1公里×1公里),采用手持式全球定位系统记录枯草样地的经纬度坐标;

步骤三、苏尼特荒漠草原样地枯草重量测定

在每个样地内随机选取3个1米×1米样方,利用传统收获法获取样方的枯草,65℃恒温烘干至恒重,求取样方枯草重量,计算3个样方枯草重量的平均值作为该样地的枯草重量;

步骤四、苏尼特荒漠草原样地植被指数计算

下载与样地枯草重量测定时间同步的覆盖样地的mod09ga反射率遥感影像,使用mrt软件对mod09ga数据进行格式与地图投影转换,hdf格式转换为tiff格式,正弦曲线地图投影转换为wgs84/albers地图投影;

依据上述步骤二中记录的样地经纬度坐标,提取位于样地内的像元的第6波段(短波红外)与第7波段(短波红外)反射率,计算如下光谱指数:

s-ndti-0.1=(1+0.1)×(b6-b7)/(b6+b7+0.1)

s-ndti-0.2=(1+0.2)×(b6-b7)/(b6+b7+0.2)

s-ndti-0.3=(1+0.3)×(b6-b7)/(b6+b7+0.3)

s-ndti-0.4=(1+0.4)×(b6-b7)/(b6+b7+0.4)

s-ndti-0.5=(1+0.5)×(b6-b7)/(b6+b7+0.5)

s-ndti-0.6=(1+0.6)×(b6-b7)/(b6+b7+0.6)

s-ndti-0.7=(1+0.7)×(b6-b7)/(b6+b7+0.7)

s-ndti-0.8=(1+0.8)×(b6-b7)/(b6+b7+0.8)

s-ndti-0.9=(1+0.9)×(b6-b7)/(b6+b7+0.9)

s-ndti-1.0=(1+1.0)×(b6-b7)/(b6+b7+1.0)

式中,b6为mod09ga第6波段(短波红外波段)反射率,b7为mod09ga第7波段(短波红外波段)反射率;

步骤五、确定苏尼特荒漠草原枯草重量的最佳光谱指数与遥感估算模型

采用一元线性分析对样地枯草重量与其对应的上述步骤四中计算的10种光谱指数进行回归拟合,如附图1所示,通过比较回归模型的决定系数(r2)与留一交叉验证均方根误差(rmsecv),确定估算枯草重量的最佳光谱指数为s-ndti-0.5,如附图2所示,确定的最佳光谱指数s-ndti-0.5对应的遥感估算模型如下:

枯草重量=6.4+844.1×s-ndti-0.5

(r2=0.77,rmsecv=13.5g/m2)

步骤六、计算苏尼特荒漠草原短波红外6波段与7波段组合的s-ndti-0.5光谱指数

获取苏尼特荒漠草原10月下旬的mod09ga反射率遥感影像,使用mrt软件将mod09ga数据进行格式与地图投影转换,hdf格式转换为tiff格式,正弦曲线地图投影转换为wgs84/albers地图投影;

依据上述步骤一描述的方法得到的苏尼特荒漠草原空间分布图,计算10月中旬苏尼特荒漠草原短波红外6波段与7波段反射率组合的s-ndti-0.5光谱指数;

步骤七、苏尼特荒漠草原枯草重量遥感估算

依据上述步骤六描述的方法计算得到的10月中旬苏尼特荒漠草原短波红外6波段与7波段反射率组合的s-ndti-0.5光谱指数,并通过估算模型“枯草重量=6.4+844.1×s-ndti-0.5”,逐像元计算10月下旬苏尼特荒漠草原枯草重量,如附图3所示。

在上述实施例中,根据记录的样地经纬度坐标,提取位于样地内对应mod09ga像元的第6波段(短波红外)与第7波段(短波红外)反射率;计算第6波段与第7波段反射率组合的s-ndti-0.1、s-ndti-0.2、s-ndti-0.3、s-ndti-0.4、s-ndti-0.5、s-ndti-0.6、s-ndti-0.7、s-ndti-0.8、s-ndti-0.9、s-ndti-1.0等10种光谱指数;采用一元线性模型对30个样地枯草重量与其所对应的计算出的10个光谱指数之间的关系进行回归拟合;通过比较拟合方程的决定系数(r2)与留一交叉验证均方根误差(rmsecv),决定系数越大留一交叉验证均方根误差越小,估算精度越高,确定估算枯草重量的最佳光谱指数为s-ndti-0.5。

与背景技术中介绍的用于估算枯草重量的ndti与sti指数相比,见附表图4,本发明方法估算枯草重量的决定系数提高了0.11,留一交叉验证均方根误差分别减小了1.4g/m2,具有重要的实用推广意义。

当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1