基于空间自相关区域植被净初级生产力时空分异探测方法技术

本发明专利技术公开一种基于空间自相关区域净初级生产力时空分异探测方法，包括如下步骤：S1、获取研究区域研究年份的MOD17A3数据；S2、提取研究区域研究时段的植被NPP图；S3、计算研究区域研究时段的植被NPP变化量；S4、计算研究区域研究时段植被NPP的全局莫兰指数和高/低值聚类指数及相关统计量；S5、判断研究区域数据样本的整体聚类特征；S6：计算研究区域研究时段植被NPP的热点分析指数；S7、判断研究区域研究时段植被NPP的热点分布特征；S8、计算研究区域研究时段逐年植被NPP变化量的局部莫兰指数；S9、判断研究区域研究时段逐年植被NPP变化量的空间聚集特征；本发明专利技术能够准确表达植被NPP的整体聚类特征，并探测出局部空间特征。

Spatial and Temporal Differentiation Detection Method of Net Primary Productivity of Vegetation Based on Spatial Autocorrelation

【技术实现步骤摘要】
基于空间自相关区域植被净初级生产力时空分异探测方法
本专利技术涉及植被生产力探测领域，特别是涉及一种空间自相关区域植被净初级生产力时空分异探测，尤其涉及一种利用卫星遥感数据快速、高效、低廉、准确探测大区域尺度植被净初级生产力时空分异特征的方法。
技术介绍
气候和人类活动对植被生产力水平区域分布有很大影响。准确探测植被净初级生产力空间分布对评估陆地生态系统的结构和功能具有重要意义，净初级生产力(NetPrimaryProductivity)简称为NPP，传统的NPP时空分布情况监测主要依赖于单一像元属性值的计算和比较，忽略了各部分之间的空间相关性，研究结果表示的是多数点的变化情况，这导致在不同的研究尺度下对同一地区的变化描述存在差异，甚至相反。在当前全球气候变化频繁与人类活动加剧的背景下，急需专利技术一种快速、高效、低廉、准确的标准方法来探测和表达区域植被净初级生产力时空分异特征。
技术实现思路
本专利技术的目的是提供一种空间自相关区域植被净初级生产力时空分异探测方法，以解决所述现有技术存在的问题。为实现所述目的，本专利技术提供了如下方案：本专利技术提供一种空间自相关区域植被净初级生产力时空分异探测方法，包括以下步骤：S1、获取研究区域研究年份的MOD17A3数据下载研究区域研究时段对应的MODIS陆地产品MOD17A3数据，然后对MOD17A3数据进行格式转换、地图投影转换，即HDF格式转换为TIFF格式，正弦曲线地图投影转换为WGS84/Albers地图投影；S2、提取研究区域研究时段的植被NPP图将步骤S1中的MOD17A3处理后的数据，用于研究区边界的矢量图对处理后的数据进行掩膜提取，空间分辨率1km×1km，数据格式为TIFF格式，投影为WGS84/AlbersEqualAreaConic；S3、计算研究区域研究时段的逐年植被NPP变化量依据步骤S2中提取的研究区域研究时段植被NPP图，利用后一年的植被NPP减去前一年的植被NPP，生成研究区域研究时段的逐年植被NPP差值图；S4、计算研究区域研究时段植被NPP的全局莫兰指数I和高/低值聚类指数G(d)及相关统计量根据步骤S2中提取的研究时段的植被NPP，计算全局莫兰指数I，公式如下：其中，Xi、Xj分别为点i和点j的观测值；是相关观测值的均值；Wij代表空间权重矩阵，若点i和点j空间相邻，矩阵第i行第j列对应的元素为1，不相邻则为0；n为观测值的数量。再计算每个全局莫兰指数值对应的期望E、P值和z得分；计算高/低值聚类指数G(d)，公式如下：其中，Wij(d)为空间权重矩阵，若点i和点j空间相邻，矩阵第i行第j列对应的元素为1，不相邻则为0；Xi、Xj分别为点i和点j的观测值。再计算每个高/低值聚类指数对应的期望E、P值和z得分；S5、判断研究区域数据样本的整体聚类特征根据步骤S4计算的研究区域研究时段的全局莫兰指数I，高/低值聚类指数G(d)及相关统计量，检验数据能否拒绝零假设，并判断是否存在聚类以及聚类类型：若全局莫兰指数I为正则表示数值相近的NPP之间呈聚集趋势，若全局莫兰指数I为负则表示呈离散趋势；若高/低值聚类指数G(d)大于期望值则表示研究区域的聚类类型为高值聚类，若高/低值聚类指数G(d)小于期望值则表示研究区域的聚类类型为低值聚类；S6、计算研究区域研究时段植被NPP的热点分析指数G*(d)根据步骤S2中提取的研究区域研究时段的植被NPP，计算研究区域研究时段的热点分析指数G*(d)，公式如下：其中，Wij(d)为空间权重矩阵，若点i和点j空间相邻，矩阵第i行第j列对应的元素为1，不相邻则为0；Xj为点j处的NPP观测值。再计算热点分析指数G*(d)对应的期望E、P值和z得分；S7、判断研究区域研究时段植被NPP的热点分布特征根据步骤S6得到的研究区域研究时段植被NPP热点分析指数G*(d)估计值，判断研究区域研究时段植被NPP的热点分布特征：如果Z(Gi*)为正且在显著区间，表明点i周围的值属高值空间集聚(热点区)；反之，如果Z(Gi*)为负且在显著区间，则表明点i周围的值属低值空间集聚(冷点区)；S8、计算研究区域研究时段逐年植被NPP变化量的局部莫兰指数Ii根据步骤S3中计算的研究区域研究时段的逐年植被NPP变化量，计算局部莫兰指数Ii，公式如下：其中，xi为区域i的观测值；为相关观测值的均值；wij为空间权重矩阵，若点i和点j空间相邻，矩阵第i行第j列对应的元素为1，不相邻则为0；n为观测值的数量；S9、判断研究区域研究时段逐年植被NPP变化量的空间聚类特征根据步骤S8中得到的研究区域研究时段逐年植被NPP变化量局部莫兰指数Ii计算结果，判断研究区域研究时段逐年植被NPP变化量的空间聚类特征，包括HH—高值围绕高值、HL—高值围绕低值、LH—低值围绕高值、LL—低值围绕低值四种类型，其中HH和LL这两种为正相关模式，HL和LH为负相关模式。本专利技术公开了以下技术效果：本专利技术基于空间自相关指数对空间信息的表达能力，构建了探测植被NPP时空分异的方法，弥补了传统方法用多数点的趋势概括总体来描述植被净初级生产力变化特征，从而导致在不同的研究尺度下对同一地区的变化描述存在差异的不足。本专利技术能够准确表达植被NPP的整体聚类特征，并探测出局部空间特征，以便将大尺度区域的研究讨论的对象转换为具有空间相关性的几部分，实现了基于卫星遥感数据的大区域范围植被净初级生产力时空分异和演变特征的快速、高效、低廉、准确探测与表达。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术方法的流程示意图；图2为2001年内蒙古呼伦贝尔市与乌兰察布市植被NPP热点分析中冷点和热点占比示意图具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。为使本专利技术的所述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本专利技术作进一步详细的说明。参照图1-2本专利技术提供一种基于空间自相关区域植被净初级生产力时空分异探测方法，包括如下步骤：S1、获取内蒙古2000-2014年的MOD17A3数据下载内蒙古2000-2014年MODIS陆地产品MOD17A3数据，使用MRT软件将MOD17A3数据进行格式转换、地图投影转换，即HDF格式转换为TIFF格式，正弦曲线地图投影转换为WGS84/Albers地图投影；S2、提取内蒙古2000-2014年植被NPP图依据所述步骤S1中获取的MOD17A3陆地产品数据，用内蒙古边界矢量图对处理后的MOD17A3数据进行掩膜提取，空间分辨率1km×1km，数据格式为TIFF格式，投影为WGS84/AlbersEqualAreaConic；S3、计算内蒙古地区逐年植被NPP变化量依据所述步骤S2提取的2000-201本文档来自技高网...

【技术保护点】
1.一种基于空间自相关区域净初级生产力时空分异探测方法，其特征在于：包括以下步骤：S1、获取研究区域研究年份的MOD17A3数据下载研究区域研究时段对应的MODIS陆地产品MOD17A3数据，然后对MOD17A3数据进行格式转换、地图投影转换，即HDF格式转换为TIFF格式，正弦曲线地图投影转换为WGS84/Albers地图投影；S2、提取研究区域研究时段的植被NPP图将步骤S1中的MOD17A3处理后的数据，用于研究区边界的矢量图对处理后的数据进行掩膜提取，空间分辨率1km×1km，数据格式为TIFF格式，投影为WGS84/Albers Equal Area Conic；S3、计算研究区域研究时段的逐年植被NPP变化量依据步骤S2中提取的研究区域研究时段植被NPP图，利用后一年的植被NPP减去前一年的植被NPP，生成研究区域研究时段的逐年植被NPP差值图；S4、计算研究区域研究时段植被NPP的全局莫兰指数I和高/低值聚类指数G(d)及相关统计量根据步骤S2中提取的研究时段的植被NPP，计算全局莫兰指数值I，公式如下：

【技术特征摘要】
1.一种基于空间自相关区域净初级生产力时空分异探测方法，其特征在于：包括以下步骤：S1、获取研究区域研究年份的MOD17A3数据下载研究区域研究时段对应的MODIS陆地产品MOD17A3数据，然后对MOD17A3数据进行格式转换、地图投影转换，即HDF格式转换为TIFF格式，正弦曲线地图投影转换为WGS84/Albers地图投影；S2、提取研究区域研究时段的植被NPP图将步骤S1中的MOD17A3处理后的数据，用于研究区边界的矢量图对处理后的数据进行掩膜提取，空间分辨率1km×1km，数据格式为TIFF格式，投影为WGS84/AlbersEqualAreaConic；S3、计算研究区域研究时段的逐年植被NPP变化量依据步骤S2中提取的研究区域研究时段植被NPP图，利用后一年的植被NPP减去前一年的植被NPP，生成研究区域研究时段的逐年植被NPP差值图；S4、计算研究区域研究时段植被NPP的全局莫兰指数I和高/低值聚类指数G(d)及相关统计量根据步骤S2中提取的研究时段的植被NPP，计算全局莫兰指数值I，公式如下：其中，Xi、Xj分别为点i和点j的观测值；是相关观测值的均值；Wij代表空间权重矩阵，若点i和点j空间相邻，矩阵第i行第j列对应的元素为1，不相邻则为0；n为观测值的数量。再计算每个全局莫兰指数值对应的期望E、P值和z得分；计算高/低值聚类指数G(d)，公式如下：其中，Wij(d)为空间权重矩阵，若点i和点j空间相邻，矩阵第i行第j列对应的元素为1，不相邻则为0；Xi、Xj分别为点i和点j的观测值。再计算每个高/低值聚类指数对应的期望E、P值和z得分；S5、判断研究区域数据样本的整体聚类特征根据步骤S4计算的研究区域研究时段的全局莫兰指数值I，高/低值聚类指数G(d)及相关统计量，检验数据能否拒绝零假设，并判断是否存在聚类以及聚类类型：若全...

【专利技术属性】
技术研发人员：任鸿瑞，尚颖洁，
申请(专利权)人：太原理工大学，
类型：发明
国别省市：山西,14