论文部分内容阅读
中国水稻种植面积居世界第二,总产量居世界第一。水稻种植面积的变化、产量丰欠等重要农情信息历来受到各级政府和社会公众的高度重视,是国家制定粮食政策和经济计划的重要依据。及时了解、准确掌握水稻种植面积、产量等对各级政府制定农业生产和农村政策。确保我国乃至世界粮食安全具有重要的意义。国家统计局承担着全国水稻产量的调查职能。经过几十年的努力,形成了较完善的统计调查方法体系。然而,随着我国社会经济的快速发展,政府决策部门、社会公众等统计服务对象对统计的需求提高,以及统计部门面临的问题和压力,如何吸收利用先进技术。逐步改善统计调查手段、降低调查成本、提高调查精度是当前面临和亟待解决的首要问题之一。以遥感为主的3S高新技术具有全覆盖、及时、客观等特点。为改进和提升传统统计手段提供了良好的技术支撑。水稻遥感估产经过多年的发展研究,随着遥感数据的更新以及技术的发展,结合统计局统计抽样调查地块实割实测水稻标准亩产数据,研究统计与遥感相结合的水稻遥感估产,实现统计部门在水稻面积产量监测及预报等业务化运行等方面具有重要意义。湖南省水稻种植面积及总产都居全国首位,且其地形地貌和水稻种植制度在中国水稻种植区具有典型性,因此,本研究以湖南省为研究区,利用MODIS不同数据产品(MOD09A1、MOD09GA、MYD09GA、MOD13Q1、MYD13Q1、MOD15A2、MOD17A2、MOD17A3),在完成水稻遥感估产分区及水稻面积提取的基础上,结合县级、乡镇级统计数据及统计局统计抽样调查地块实割实测数据,以及基于像元水平的MODIS GPP/NPP,研究不同数据完备条件下,水稻遥感估产方法。主要方法与结论如下:(1)基于空间分析与两维图论树算法相结合的湖南省水稻遥感估产分区在分析研究区(湖南省)基本特征的基础上,以水稻种植制度、稻土此、水稻单产水平、地形及县域空间位置作为分区指标,采用空间分析与两维图论树算法相结合的方法,将湖南省分为两个一级区,五个二级区。(2)基于水稻典型发育期光谱特征与多时相MODIS数据湖南省水稻面积遥感信息提取在对MOD09A1数据受云污染的像元进行基于MODIS数据的质量评价(Qualityassessment:QA)信息的相邻时相最大值法插补的基础上,利用水稻移栽期、抽穗期的典型光谱特征、以及相应的EVI与LSWI的关系,排除非水稻植被及其它因素的影响。根据研究区水稻种植制度特点,分别对2000-2008年早稻、晚稻及一季稻面积进行提取。并利用县级统计数据、基于高分辨率遥感影像数字化的水稻田数据、1:1万的土地利用数据,对提取结果进行验证评估。结果显示:湖南省水稻主要分布在洞庭湖流域及湖南省中部地区,相对于以前的研究,本研究的提取精度有较大的提高。(3)基于水稻遥感估产分区、水稻面积提取、县级统计数据的省级水稻总产遥感估算在水稻遥感估产分区、水稻面积提取的基础上,利用1:25万县级行政区划图,提取2000-2008年每年的早稻、晚稻及一季稻对应的各县各时相MOD09A1数据,选取EVI为估产参数。利用水稻生长发育期数据,选择移栽期、孕穗期、抽穗期、乳熟期及成熟期为水稻遥感估产的主要生长期,以未基于分区与基于分区两种思路,建立2000-2007年各县EVI平均值乘以第三章提取的水稻面积的结果(AEVI)与统计总产的一次线性、二次非线性及各生育期逐步回归模型。通过误差分析,选择水稻总产最优遥感拟合模型。在此基础上,利用基于2000-2007年数据得到的水稻总产最优遥感拟合模型预测2008年的省级水稻总产。结果显示:基于分区的估产模型要比未基于分区的水稻遥感估产模型要好,且生育期主要集中于孕穗期到乳熟期,说明这段时间是水稻产量形成的最关键时期。另外。二次非线性模型及逐步回归模型的结果要好。误差分析结果显示,水稻遥感估产模型的拟合结果比预测结果要好,两者省级相对误差都小于5%。(4)基于统计局统计抽样调查地块实割实测数据、MOD13Q1与MYD13Q1EVI的省级水稻单产遥感估算基于国家抽样县抽样地块的实测数据是统计部门计算水稻单产的基础数据,本研究根据地块的空间位置,提取各地块所在地的3×3像元网格对应的MOD13Q1与MYD13Q1相结合的EVI数据,选择移栽期、孕穗期、抽穗期、乳熟期及成熟期的EVI,建立地块实测标准亩产与EVI的一次线性、二次非线性及各生育期逐步回归模型。通过误差分析,选择水稻单产最优遥感拟合模型。在此基础上,利用上一年的水稻总产最优遥感拟合模型预测下一年的省级水稻单产。根据水稻面积提取结果,得到省级水稻总产。与统计值进行比较,结果显示,在省级水平上,基于地块的水稻单产与总产遥感拟合与预测结果相对误差都接近5%。(5)基于乡镇级统计数据、MOD13Q1与MYD13Q1 EVI的的县级水稻总产遥感估算选择醴陵市为研究区,分析基于乡镇级统计数据的水稻遥感估产。选择分别由Terra与Aqua卫星获取得到的MOD13Q1与MYD13Q1 EVI作为水稻估产遥感指标。其空间分辨率约为250 m。时间分辨为8 d。选择移栽期、孕穗期、抽穗期、乳熟期及成熟期的遥感数据,建立EVI乘以水稻面积的结果(AEVI)与乡镇级水稻统计总产的一次线性、二次非线性及各生育期逐步回归模型。通过误差分析,选择水稻总产最优遥感拟合模型。在此基础上,利用上一年的水稻总产最优遥感拟合模型预测下一年的醴陵市水稻总产。估产结果显示,各拟合模型相对误差都小于0.1%。预测结果误差相对于拟合结果较大,但相对误差仍然小于5%。(6)基于像元水平的MODIS GPP/NPP数据的水稻遥感估产、及像元内水稻面积比对估产精度的影响在前人研究的基础上,总结得到基于MODIS GPP/NPP水稻估产方法。利用MODIS 8d的、年的GPP/NPP产品及产品使用说明,估算8 d的NPP,并利用此结果与8d的GPP产品分别进行基于GPP/NPP的水稻遥感估产。借助MODIS GPP/NPP产品生成1×1km网格。研究像元内不同水稻面积比对估产结果的影响。研究结果表明:利用估算年NPP的方法、8 d的LAI和PSNnet来估算每8 d的NPP的结果比MODIS年NPP值较大。基于MODIS NPP的水稻遥感估产结果相对于统计值较小.而基于MODI GPP水稻遥感估产与统计值接近。其估产精度随着像元内水稻面积的增加而减小。像元内水稻面积超过80%,其相对误差小于20%或10%。综上所述,本研究的主要结论如下:基于县级、乡镇级、实测地块数据、MODISGPP数据的水稻遥感估产的结果基本上都能达到95%的精度。从提高精度方面考虑,在乡镇级统计数据具备的前提下,利用基于乡镇统计数据及MOD13Q1、MYD13Q1EVI的方法是最优选择。从水稻估产业务化运行角度来考虑,在地块实测产量及空间位置数据都具备的前提下,利用基于统计局统计抽样调查地块实割实测数据及MOD13Q1、MYD13Q1 EVI的方法较好。在统计局统计抽样调查地块实割实测数据不具备的情况下,利用分区与水稻面积提取结果、以及县级统计数据的水稻遥感估产方法。而对于大面积水稻种植区,利用基于像元水平的MODIS GPP/NPP估产方法较好。