一种更新全球作物总产量栅格数据的方法技术

本发明专利技术公开了一种更新全球作物总产量栅格数据的方法，包括数据准备、统计分析、高通滤过处理、数据导出、土地利用重分类、聚合分析等步骤。利用本发明专利技术的方法能够对全球175种不同作物的总产量栅格数据进行及时更新，数据时效性较强，数据精度较高。

A Method for Updating the raster Data of Global Crop Yield

【技术实现步骤摘要】
一种更新全球作物总产量栅格数据的方法
本专利技术属于资源环境
，尤其涉及一种基于空间分析技术更新作物总产量栅格数据的方法。
技术介绍
地球上的农田面积约占1500万平方公里，为人类提供了大部分的食物，对人类福祉至关重要。尽管大多全球土地覆盖的卫星群涉及农田的数据集只有几个类别，但对于回答从生物多样性保护到粮食安全再到生物地球化学循环等问题等关键问题提供了至关重要的信息。有关农作物选择、产量和肥料使用等农业土地利用措施的信息是更有限。Monfreda等提出了通过结合国家、州的土地使用数据集和县级人口普查统计，更新的全球农田数据集为5*5分(大约10公里乘10公里)经纬度网格，生成了2000年全球175种不同作物的总产量栅格数据。然而，该数据已经10多年没有再进行更新，数据缺乏时效性，加上数据分辨率低，为5*5分的分辨率，因此，提高数据精度和更新作物总产量栅格数据显得尤为重要。
技术实现思路
本专利技术所要解决的技术问题是针对现有技术的不足提供一种更新全球作物总产量栅格数据的方法。本专利技术的技术方案如下：一种更新全球作物总产量栅格数据的方法，包括以下步骤：A1数据准备：选取用作更新的土地利用栅格数据；A2利用ArcGIS的工具CalculateStatistics进行统计分析，便于进行分类，输出文件landuse_stat.gif；A3高通滤过处理；为减少数据的碎化，使用MajorityFilter工具进行高通滤波处理，选择参数为周围邻近的8个像元，替代阈值方法为majortiy(众数)，输出栅格层landuse_major.itf；A4数据导出；将A3的landuse_major.gif数据导出，输出TIFF格式文件，文件名landuse_export.gif，默认输出栅格大小为0.008333333，栅格大小为原产量栅格的栅格大小(0.08333333)的1/10；A5对土地利用重分类；在landuse_export.gif基础上，利用ArcGIS中的Reclassify工具对土地利用进行重分类，对数据进行二值化分类，1表示选取的土地利用，0表示不选取的土地利用，输出文件名为landuse2value.gif；A6聚合分析；在上一步基础上，使用ArcGIS中的Aggregate工具进行聚合分析，根据情况，选择适合的栅格大小因子；输出文件为landuse.gif；A7总产量空间数据以及作物列表、国家表数据准备；准备2000年不同作物总产量栅格数据，要求作物名作为文件名，文件扩展名为：tif，如wheat.gif，统一放在文件夹production目录下，便于调用；准备作物列表数据，表名称为croplist；准备国家表country_table,必须包含国家名称字段Countryname；A8在2000年作物总产量栅格数据的基础上，按照国家提取每种作物的总产量数据；A9产量比值文件准备；准备作物产量统计文件，若用2010年的数据对2000年的产量栅格数据进行更新，则需要2010年的作物产量统计数据，生成2010年与上一步提取的2000年的作物总产量的比值数据Production_ratio；A10国家行政区划数据准备及挂接；准备世界国家行政区划空间数据Countries，通过Join命令，将countries与production_ratio基于国家名称字段进行挂接关联，生成country_production；A11迭代栅格数据；使用iterateraster迭代器，对产量栅格数据进行迭代，以方便进行循环操作，栅格目录选择production文件夹，栅格格式选择TIF，输出文件为每次迭代时的产量栅格数据，迭代输出变量为name；若迭代输出的文件名为wheat.gif，则name为wheat.gif；A12土地利用更新；上一步的输出数据与步骤A6生成的土地利用数据(landuse.gif)进行Times(乘积)运算，输入栅格数据1为：landuse.gif，输入栅格数据2为迭代器输出的产量栅格数据，栅格大小为0.041666665，数据范围选择：unionofinput(输入数据的并集范围)，输出文件夹为landupdate文件夹，输出文件名为land_％name％，％name％为表示引用name变量；若迭代输出的文件名为wheat.gif，则name为wheat.gif，输出的文件名则为：land_wheat.gif；A13进行产量数据校正；由于在进行栅格转换时，原始的栅格大小栅格大小0.08333333，新的栅格大小栅格大小为0.041666665，原始栅格大小是新栅格大小的2倍，新产量栅格数据应除以4；使用divide工具，进行数据校正，输入数据1为：land_％name％，输入数据2为：4，输出数据为：divide_％name％；A14记录提取；基于步骤A7的croplist表，用Tableselect工具提取记录；输入表为：croplist，输出位置为：Scratch.mdb，输出表名为：Table_select，输出表达式为：[CROPNAME]＝'％cropname％'；A15字段提取；使用Getfieldvalue工具，提取Table_select表cropname字段中的记录，并作为字段；输入表为：Table_select，字段为：cropname，数据类型为：Field，输出变量为：value；A16进行栅格转换；使用Featuretoraster工具，输入图层为：country_production，字段选取value变量，输出栅格大小为：0.041666665，输出图层为：Production0_％value％；A17对数据进行后处理；进行Nodata(空值)检测，并将nodata值设置为0；使用Isnull工具判别是否空值，生成临时栅格文件Null_％value％.gif；再使用con(条件查询)工具，表达式为"VALUE"＝1，条件为真时(即为Nodata数据时)，值为0，否则等于Production0_％value％值，输出文件为Production_％value％.gif；A18总产量数据更新；在步骤A13和步骤A17输出结果的基础上，进行乘积运算，输入文件1为：Production_％value％.gif，输入文件2为：divide_％name％；，输出Production_up_％value％.gif；A19循环步骤A11-18，直至读取完所有作物总产量栅格数据；完成后，即完成基于土地利用和2010年统计数据对2000年作物总产进行更新，并且提高了数据精度。所述的更新全球作物总产量栅格数据的方法，所述步骤A8包括以下步骤：1)使用iterateraster迭代器，对总产量栅格数据进行迭代，2)通过parsepath工具对迭代输出的文件名进行路径解析，解析类型选择name；提取的文件名不包含扩展名，获得作物名，输出变量为cropname，变量引用方式为：％cropname％；3)使用Zonalstatisticsastable(分区统计)工具按照国家进行提取，输入分区数据为：Countries.shp，分区统计字段为：FAO，输入栅格数据为：每次迭代的总产量栅格数据本文档来自技高网...

【技术保护点】
1.一种更新全球作物总产量栅格数据的方法，其特征在于，包括以下步骤：A1数据准备：选取用作更新的土地利用栅格数据；A2利用ArcGIS的工具Calculate Statistics进行统计分析，便于进行分类，输出文件landuse_stat.tif；A3高通滤过处理；为减少数据的碎化，使用Majority Filter工具进行高通滤波处理，选择参数为周围邻近的8个像元，替代阈值方法为majortiy(众数)，输出栅格层landuse_major.itf；A4数据导出；将A3的landuse_major.tif数据导出，输出TIFF格式文件，文件名landuse_export.tif，默认输出栅格大小为0.008333333，栅格大小为原产量栅格的栅格大小(0.08333333)的1/10；A5对土地利用重分类；在landuse_export.tif基础上，利用ArcGIS中的Reclassify工具对土地利用进行重分类，对数据进行二值化分类，1表示选取的土地利用，0表示不选取的土地利用，输出文件名为landuse2value.tif；A6聚合分析；在上一步基础上，使用ArcGIS中的Aggregate工具进行聚合分析，根据情况，选择适合的栅格大小因子；输出文件为landuse.tif；A7总产量空间数据以及作物列表、国家表数据准备；准备2000年不同作物总产量栅格数据，要求作物名作为文件名，文件扩展名为：tif，如wheat.tif，统一放在文件夹production目录下，便于调用；准备作物列表数据，表名称为croplist；准备国家表country_table,必须包含国家名称字段Countryname；A8在2000年作物总产量栅格数据的基础上，按照国家提取每种作物的总产量数据；A9产量比值文件准备；准备作物产量统计文件，若用2010年的数据对2000年的产量栅格数据进行更新，则需要2010年的作物产量统计数据，生成2010年与上一步提取的2000年的作物总产量的比值数据Production_ratio；A10国家行政区划数据准备及挂接；准备世界国家行政区划空间数据Countries，通过Join命令，将countries与production_ratio基于国家名称字段进行挂接关联，生成country_production；A11迭代栅格数据；使用iterate raster迭代器，对产量栅格数据进行迭代，以方便进行循环操作，栅格目录选择production文件夹，栅格格式选择TIF，输出文件为每次迭代时的产量栅格数据，迭代输出变量为name；若迭代输出的文件名为wheat.tif，则name为wheat.tif；A12土地利用更新；上一步的输出数据与步骤A6生成的土地利用数据(landuse.tif)进行Times(乘积)运算，输入栅格数据1为：landuse.tif，输入栅格数据2为迭代器输出的产量栅格数据，栅格大小为0.041666665，数据范围选择：union of input(输入数据的并集范围)，输出文件夹为landupdate文件夹，输出文件名为land_％name％，％name％为表示引用name变量；若迭代输出的文件名为wheat.tif，则name为wheat.tif，输出的文件名则为：land_wheat.tif；A13进行产量数据校正；由于在进行栅格转换时，原始的栅格大小栅格大小0.08333333，新的栅格大小栅格大小为0.041666665，原始栅格大小是新栅格大小的2倍，新产量栅格数据应除以4；使用divide工具，进行数据校正，输入数据1为：land_％name％，输入数据2为：4，输出数据为：divide_％name％；A14记录提取；基于步骤A7的croplist表，用Table select工具提取记录；输入表为：croplist，输出位置为：Scratch.mdb，输出表名为：Table_select，输出表达式为：[CROPNAME]＝'％cropname％'；A15字段提取；使用Get field value工具，提取Table_select表cropname字段中的记录，并作为字段；输入表为：Table_select，字段为：cropname，数据类型为：Field，输出变量为：value；A16进行栅格转换；使用Feature to raster工具，输入图层为：country_production，字段选取value变量，输出栅格大小为：0.041666665，输出图层为：Production0_％value％；A17对数据进行后处理；进行Nodata(空值)检测，并将nodata值设置为0；使用Isnull工具判别是否空值，生成临时栅格文件Null_％value％....

【技术特征摘要】
1.一种更新全球作物总产量栅格数据的方法，其特征在于，包括以下步骤：A1数据准备：选取用作更新的土地利用栅格数据；A2利用ArcGIS的工具CalculateStatistics进行统计分析，便于进行分类，输出文件landuse_stat.gif；A3高通滤过处理；为减少数据的碎化，使用MajorityFilter工具进行高通滤波处理，选择参数为周围邻近的8个像元，替代阈值方法为majortiy(众数)，输出栅格层landuse_major.itf；A4数据导出；将A3的landuse_major.gif数据导出，输出TIFF格式文件，文件名landuse_export.gif，默认输出栅格大小为0.008333333，栅格大小为原产量栅格的栅格大小(0.08333333)的1/10；A5对土地利用重分类；在landuse_export.gif基础上，利用ArcGIS中的Reclassify工具对土地利用进行重分类，对数据进行二值化分类，1表示选取的土地利用，0表示不选取的土地利用，输出文件名为landuse2value.gif；A6聚合分析；在上一步基础上，使用ArcGIS中的Aggregate工具进行聚合分析，根据情况，选择适合的栅格大小因子；输出文件为landuse.gif；A7总产量空间数据以及作物列表、国家表数据准备；准备2000年不同作物总产量栅格数据，要求作物名作为文件名，文件扩展名为：tif，如wheat.gif，统一放在文件夹production目录下，便于调用；准备作物列表数据，表名称为croplist；准备国家表country_table,必须包含国家名称字段Countryname；A8在2000年作物总产量栅格数据的基础上，按照国家提取每种作物的总产量数据；A9产量比值文件准备；准备作物产量统计文件，若用2010年的数据对2000年的产量栅格数据进行更新，则需要2010年的作物产量统计数据，生成2010年与上一步提取的2000年的作物总产量的比值数据Production_ratio；A10国家行政区划数据准备及挂接；准备世界国家行政区划空间数据Countries，通过Join命令，将countries与production_ratio基于国家名称字段进行挂接关联，生成country_production；A11迭代栅格数据；使用iterateraster迭代器，对产量栅格数据进行迭代，以方便进行循环操作，栅格目录选择production文件夹，栅格格式选择TIF，输出文件为每次迭代时的产量栅格数据，迭代输出变量为name；若迭代输出的文件名为wheat.gif，则name为wheat.gif；A12土地利用更新；上一步的输出数据与步骤A6生成的土地利用数据(landuse.gif)进行Times(乘积)运算，输入栅格数据1为：landuse.gif，输入栅格数据2为迭代器输出的产量栅格数据，栅格大小为0.041666665，数据范围选择：unionofinput(输入数据的并集范围)，输出文件夹为landupdate文件夹，输出文件名为land_％name％，％name％为表示引用name变量；若迭代输出的文件名为wheat.gif，则name为wheat.gif，输出的文件名则为：land_wheat.gif；A13进行产量数据校正；由于在进行栅格转换时，原始的栅格大小栅格大小0.08333333，新的栅格大小栅格大小为0.041666665...

【专利技术属性】
技术研发人员：周卫，何萍，徐新朋，雷秋良，梁霞，梁国庆，艾超，
申请(专利权)人：中国农业科学院农业资源与农业区划研究所，
类型：发明
国别省市：北京,11