中国耕地熟制变化遥感监测探究.pdfVIP

论文题目：中国耕地熟制变化遥感监测研究 专 业：地图学与地理信息系统 硕 士 生：史 姝 （签名）_____________ 指导教师：汤伏全 （签名）_____________ 李存军 （签名）_____________ 摘 要 多熟种植是我国重要的耕作方式，它对保持和增加粮食的产量以及促进农村经济的 发展都有十分重要的意义。我们国家耕地的熟制格局复杂多样，在过去的三十年里，我 国的耕地资源在数量、质量和利用结构方面因受自然条件以及农村社会状况的影响发生 了复杂的变化，这就需要我们能够很好地掌握全国耕地资源的情况和粮食生产的形势， 以对估计粮食产量的变化及其原因做出科学的决策。 本文以全国范围为研究区域，利用北京师范大学提供的 GLASS LAI 数据，用 TIMESAT 软件对时间序列的GLASS LAI 数据进行重构，得到了能体现不同熟制的时序 GLASS LAI 平滑曲线，运用R 语言编程求平滑曲线极大值个数（即曲线波峰数目）的 方法实现了全国范围内耕地熟制信息的自动提取，将2000 年中国1:10 万的地表覆盖数 据的耕地层数据对结果进行裁剪，最终得到1982 年、1992 年、2002 年、2012 年的全国 耕地熟制格局分布图，对我国近三十年来耕地的熟制变化进行了分析，并计算得出复种 指数，将所得结果分别与统计资料及前人的研究成果进行了对比验证。主要的研究方法 及结果如下： （1）对我国32 个省 （直辖市、自治区）共12 年的时间序列GLASS LAI 数据进行 重投影、拼接、裁剪等预处理，建立时间序列GLASS LAI 曲线，由于影像中存在噪声 的影响，时序GLASS LAI 曲线偶尔会出现锯齿形，利用以Matlab 为运行载体的TIMESAT 3.0 软件中提供的Gaussian 滤波、Logistic 滤波以及 Savitzky-Golay 滤波 （简称S-G 滤 波）对曲线进行重构。结果显示，Gaussian 和Logistic 两种滤波不能很好的拟合二熟、 三熟作物的曲线，这会造成监测结果与实际情况的严重不符，而S-G 滤波能够更好的拟 合不同熟制的时序曲线，这能够较好地展现农作物季节性的变化情况，有利于提取波峰 数量。 （2 ）依据植物的物候特征所形成的时序曲线的变化特点，耕地熟制就相当于时序 曲线上满足条件的极大值点的个数。将上述得到的能体现不同熟制的时序 GLASS LAI 平滑曲线以时序数字的形式存储成文本格式，运用R 语言编程的方法来提取曲线上满足 条件的极大值点的个数。结果显示，此方法可以提取每个像元的时序曲线波峰个数（0 万方数据 代表非耕地，1 代表一年一熟，2 代表一年两熟，3 代表一年三熟）。 （3 ）从最终得到的1982 年、1992 年、2002 年、2012 年的全国耕地熟制格局分布 图中可以明显地看出我国耕地熟制从北到南由单一到复杂，由一熟区到二熟区再到三熟 区的渐变分布，一熟区、二熟区多集中分布，三熟区在零星分布的基础上存在部分小区 域的集中。近三十年来，我国耕地熟制的动态变化情况为：一熟区总体变动不明显，但 区域范围内有变动，部分区域增加的同时部分区域缩减；二熟区明显增加，以河南省为 中心向四周扩展；三熟区呈零星区域的出现并扩展缓慢。 （4 ）从ArcGIS 成图的属性信息中获取我国32 个省 （直辖市、自治区）一熟、二 熟、三熟的像元个数来计算复种指数。结果表明，近三十年来我国的复种指数整体呈持 续上升的趋势，这说明越来越多的耕地实现了一年两熟和一年三熟的耕作方式。将2002 年的复种指数与统计资料以及前人的研究结果进行对比。结果表明，本方法监测的结果 与统计资料和前人监测的结果均呈正相关关系，相关系数分别为0.788 和0.903，验证了 本文结果的可靠性。