基于遥感与气象数据的冬小麦主产区籽粒蛋白质含量预报

王琳 梁健 孟范玉 孟炀 张永涛 李振海

摘要：开展小麦籽粒蛋白质含量的监测预报研究对于指导农户调优栽培、企业分类收储、期货小麦价格、进口政策调整等具有重要意义。本研究以冬小麦主产区（河南省、山东省、河北省、安徽省和江苏省）为研究区域，构建了冬小麦籽粒蛋白质含量多层线性预测模型，并实现了2019年冬小麦蛋白质含量预报。为了解决预测模型在年际扩展和空间扩展存在偏差的问题，在蛋白质含量估算模型中考虑了气象因素（温度、降水、辐射量）、冬小麦筋型、抽穗一开花期增强型植被指数（EVI）等因素。结果表明，融合3个气象因素的蛋白质含量估算模型建模集精度（R=0.39，RMSE=1.04%）与验证集精度（R=0.43、RMSE=0.94%）均高于融合2个气象因子的估算模型和单个气象因子的估算模型。将蛋白质含量估算模型应用冬小麦主产区的蛋白质含量遥感估算，得到了2019年冬小麦主产区品质预报图，并形成黄淮海地区冬小麦品质分布专题图。本研究结果可同时为后续小麦种植区划和实现绿色、高产、优质、高效粮食生产提供数据支撑。

关键词：冬小麦;籽粒蛋白质含量;遥感;多层线性模型;气象数据

中图分类号：S127;S512.1+1文献标志码：A文章编号：202103-SA007

引用格式：王琳，梁健，孟范玉，孟炀，张永涛，李振海.基于遥感与气象数据的冬小麦主产区籽粒蛋白质含量预报[J].智慧农业（中英文），2021， 3（2）： 15-22.

WANG Lin， LIANG Jian， MENG Fanyu， MENG Yang， ZHANG Yongtao， LI Zhenhai. Estimating grain protein content of winter wheat in producing areas based on remote sensing and meteorological. data[J]. Smart Agriculture，2021， 3（2）： 15-22. （in Chinese with English abstract）

1引言

随着社会经济的日益发展和生活水平的不断提高，当前人们对农产品的需求已经从“有没有”“够不够”转向“好不好”“优不优”。但是，中国农业大而不强、多而不优的问题依然存在。2019年，农业农村部、国家发展和改革委员会、科学技术部等联合发布的《国家质量兴农战略规划（2018—2022）》中提出“推进农业由增产导向转向提质导向”和“走质量兴农之路”等发展目标。籽粒蛋白质含量是小麦品质评价的关键指标，而品质的高低决定小麦收购价格、加工用途和使用价值等[1]。因此，开展小麦籽粒蛋白质含量的监测预报研究，对于指导农户调优栽培、企业分类收储、期货小麦价格、进口政策调整等具有重要意义[2]。

基于作物籽粒蛋白质含量（Grain Protein Content，GPC）遥感监测预报的研究目前已得到探索和初步应用，可以归纳为4类：（1）基于“遥感信息-籽粒蛋白质含量”模式的经验模型，该模型通过分析作物关键生育时期的遥感信息（敏感波段、植被指数、红边参数等光谱特征）直接构建作物GPC统计经验模型[3-7];（2）基于“遥感信息-农学参数-籽粒蛋白质含量”模式的定量模型，该模型根据遥感信息与关键生育期农学参数之间的定量关系及农学参数与GPC之间的定量关系，构建GPC预测模型[8-10];（3）基于遥感数据和生态因子的GPC半机理模型，该模型考虑了作物氮素运转机理和生态因子对GPC的影响，改善了模型的年际扩展性和空间转移性[11-14];（4）基于遥感信息和作物模型结合的机理解释模型，该模型综合考虑籽粒蛋白质形成过程中各种生态因子的影响，通过遥感信息和作物生长模型耦合的同化方法，调整模型模拟变量与遥感观测值的误差达到最小，以调整作物模型的初始参数和状态变量，进而实现GPC的预测[14，15]。综合分析以上各类方法与模型的选择，前两类模型研究较多，操作简单且易实现，但机理解释性不强，在区域间和年際间扩展应用过程中监测结果偏差较大;同化模型方法考虑过多的输入变量，并且复杂的同化算法、品质生长模型本地化精度和高耗时运算时间问题，限制了该类方法的区域大面积应用[16]。Li等[12]和Xu等[13]通过综合考虑影响品质遥感预报模型中环境变异因素及品种筋型因子，引入分层线性模型（Hierarchical. Linear Models，HLM）解释籽粒蛋白质含量-遥感-环境的嵌套问题，构建基于开花期的冬小麦品质预测模型，很好地解决了模型年际扩展和空间扩展存在较大偏差的问题。因此，发展以遥感信息（反映作物养分状况）与环境因子（体现时空变异）结合的GPC半机理预测方法，构建普适性及适用于时空扩展的籽粒蛋白质遥感预测模型具有可行性。

本研究拟在前期构建的冬小麦品质分层遥感预测的基础上[12，13]，初步尝试在全国冬小麦主产区（河南省、山东省、河北省、安徽省和江苏省）开展GPC模型的构建，以期实现冬小麦品质收获前预测。

2数据与方法

2.1冬小麦研究区

试验分别于2008年、2009年和2019年在全国冬小麦5个主产省（河南省、山东省、河北省、安徽省和江苏省）进行数据采集，包括样点取样及测试、小麦品种及筋型确定、气象数据收集和遥感数据获取，共获取样品898个，去除异常点后得到样品864个，其中2008年200个，2009年283个，2019年381个，测样点分布如图1 所示。其中，强筋小麦86个，强筋和中筋混合区小麦249个，中筋小麦380个，中筋和弱筋混合区小麦2个，弱筋小麦147个。

2.2数据获取

2.2.1冬小麦籽粒蛋白质测定

在冬小麦成熟期，获取农户冬小麦样品和样点经纬度，每个农户收割3个样点，每个样点1m，经过自然晾晒、风干、称重和求均值。最后利用近红外谷物分析仪FOSS Infratec 1241（Tecator，赫加奈斯，瑞典）测定冬小麦GPC（14%干基）并记录。

2.2.2冬小麦品种及筋型确定

在冬小麦样品测试过程中，同时记录农户种植冬小麦的品种信息，通过查询品种的信息确定品种的强/中/弱筋特性。区域品质预报过程中，考虑品种的多样性及复杂性，参考2001年农业部发布的《中国小麦品质区划方案》所指定的依据和原则，利用全国农业技术推广服务中心收集汇总的品种推广数据，以及以县为单位发布的当年主推品种，确定每个县种植小麦筋型特性。最终形成以县为单位的冬小麦品质区划图。

2.2.3EOS/MODIS遥感影像

以抽穗—开花期为时间节点进行冬小麦品质预报，本研究选择基于中分辨率成像光谱仪（Moderate-resolution Imaging Spectroradiometer，MODIS）遥感影像构建的增强型植被指数（Enhanced Vegetation Index，EVI）进行建模[13]。如公式（1）。

EVI=2.5×（R-R）/（R+6×R-7.5×R）（1）

其中，R、R和R分别代表近红外波段、红光波段和蓝光波段的反射率。为消除区域物候差异与有云覆盖影响，分别获取2019年和2020年4月—5月MODIS数据，分析获取对应像元最大值作为抽穗—开花期的影像结果。

2.2.4气象数据

冬小麦研究区的气象栅格数据从欧洲中期预报中心（European Center for Medium-Range Weather Forecasts，ECMWF）（http：//www.ecmwf.int/）获取。本研究气象数据主要包括4月份的月积温（℃·日）、月太阳辐射量（MJ/m）和月总降雨量（mm），选择这3个气象因子进行冬小麦品质预报是因为研究区内冬小麦从起身期到开花期约为一个月时间（不分年份），而该生育阶段包括的拔节期、挑旗期、抽穗期及开花期是冬小麦生长最重要的生育时期，对于群体的建成、穗粒数以及养分积累具有决定作用。利用Python 3.7（Python Software Foundation，波特兰，美国）及Matlab 2014（Mathworks Inc，纳蒂克，美国）分别进行数据下载及读取。

2.3冬小麦品质预报模型

分层线性模型（Hierarchical. Linear Model，HLM）是一种考虑了数据嵌套结构（例如学生镶嵌于班级）的最小二乘回归分析，然而相比一般最小二乘回归分析，HLM在针对数据之间的不独立性时，将数据集分层、综合分析层内数据（第一层模型）以及层内与层外数据（第二层模型）之间的关系[17]。

在考虑籽粒品质的区域及年际模型构建时，地域及年际之间的环境因素（第二层模型）对品质的影响会造成遥感信息与品质模型（第一层模型）的差异，因此考虑HLM模型进行品质预测建模具有可行性[12，13]。本研究中，第一层模型是基于GPC与植被指数和籽粒筋型构建，公式如下。

GPC=β+β·EVI+β·Glu+r（2）

其中，GPC为籽粒蛋白质含量;Glu为冬小麦品种筋型值。本研究依据Li等[12]的研究，强筋、中强筋和弱筋小麦品质分别取值1、2和3，后期在区域应用过程中有强筋和中筋种植区与中筋和弱筋种植区，分别取值1.5和2.5。EVI为增强型植被指数，Xu等[13]研究表明EVI与冬小麦品质较其他植被指数具有更好的相关性，因此本研究直接选择了EVI指数作为模型输入参数。β、β和β分别为第一层模型的截距、植被指数系数和品种筋型系数。r表示随机误差。第二层模型基于第一层模型中的模型系数与外部气象数据（温度、降水、辐射量）的归一化值构建，如公式（3）。

β=γ+γ·nRad+γ·nTem+γ·nPre+μ（3）

其中，nRad、nTem和nPre分别为每个栅格4月份月太阳辐射量归一化值、月积温归一化值和月总降雨量归一化值;β分别对应第一层模型中截距、增强型植被指数系数和品種筋型系数;γ为第二层级的截距，γ、γ和γ分别为第二层级总辐射量、总降水量和积温的模型系数，μ则表示为随机误差。气象数据的归一化见公式（4）（以Rad为例）。

nRad=（Rad-Rad）/（Rad-Rad）（4）

其中，Rad、Rad和Rad为分别每个栅格4月份月太阳辐射量、最大太阳辐射量和最小太阳辐射量，MJ/m。为进一步测试模型，考虑不同的气象因子组合并对比模型的精度，以确定最优的气象因子。模型的构建与验证通过Matlab 2014实现。

2.4统计分析

经过异常点剔除后获得2008年、2009年和2019年样点数据共计864组，将三年数据混合后随机选择其中80%数据用于建模（共计691组），剩下20%数据用于模型精度验证（共计173组）。利用Matlab 2014实现对冬小麦GPC预测值与实测值之间的统计分析。采用决定系数R和均方根误差验证预测模型的精确度与可靠性。R值越接近于1，说明模型拟合程度越高。RMSE值越小，说明模型预测能力越强，稳定性好且可靠性高。另外，为进一步分析冬小麦品质预测模型中各个自变量因子对最终GPC的贡献度，采用相对重要性分析法（Dominance Analysis）[18]对模型中各个自变量因子对GPC的相对重要性进行计算和排序。

3结果与分析

3.1冬小麦品质区划分析

参考《中国小麦品质区划方案》所指定的依据和原则，结合5个冬小麦主产省各自的品质区划，以及当前主推品种的面积及市县分布，最终形成冬小麦品质区划分布图（图2）。

总体来看，强筋小麦主要以北部地区种植为主，如河北中南部、河南省北部、山东省临沂市和滨州市等。中筋小麦以山东省种植面积最广、河南南部地区、江苏省北部和安徽省北部等也有分布。弱筋小麦则主要分布在江苏南部及沿海地区、安徽中南部及河南南部地区等。通过调查发现，冬小麦的品质区划除了与生态因子（气象及土壤等）相关以外，也与品种推广和政府决策相关，例如2019年山东省济麦22、鲁原502和山农28号等中筋品种推广占总面积的95.3%，虽然山东省适合强筋和中强筋小麦种植，但综合来看强筋小麦品种占比不大。另外，在政府决策方面，2019年河北省发布的《河北省强筋小麦产业提质增效推进方案（2019—2022年）》确定以藁城、赵县、栾城、元氏、柏乡、隆尧、任县、宁晋、南和、沙河、大名、肥乡等为主要推广区县。本研究所确定的小麦品质区划专题图是在小麦品质区划和2019年各市县的主推品种调查分析后形成的。在确定过程中，受县级单位信息发布不全、部分市县品种复杂等影响，在最终确定品种筋型过程中，以市为单位确定较多，后续的研究及调查过程中需要更精确的品种推广范围，以确定准确的冬小麦筋型分布信息，为后续的品质预测模型提供可靠的筋型空间分布信息。

3.2GPC预报模型的气象因子选择

以不同的气象因子参与分层线性模型进行冬小麦GPC预报模型构建，并对比模型精度进行气象因子的评价（表1）。结果表明，单独一个气象因子进行品质预报时，以nRad为输入变量构建的冬小麦GPC预报模型精度（R=0.31，RMSE=1.12%）均高于分别以nTem和nPre为输入变量构建的冬小麦GPC预报模型。两个气象因子参与冬小麥GPC预报模型的精度有一定的提高，以nRad和nPre组合构建的冬小麦GPC预报模型较好（R=0.38，RMSE=1.05%）。三个气象因子共同参与冬小麦GPC预报精度最佳，模型的R和RMSE分别为0.39和1.04%。考虑到冬小麦籽粒蛋白质形成过程中受辐射量、温度和降水的综合影响，在区域大面积模型构建过程中，气象因子存在较大的空间变异性。有关气象因子对GPC预报模型在不同尺度范围上的影响，有必要在品质关联的气象要素的指标筛选及关联度分析等方面进行进一步研究，特别是气象因素在区域的空间变异。

进一步通过相对重要性分析GPC预报模型中各个自变量因子对GPC预测的贡献。从表2结果中可以看出，冬小麦品种对GPC的影响是最大的，Glu相对重要性达到75.31%，品种的遗传特性直接决定品质的高低，考虑到较大的研究区域，品种的多样性对模型结果具有较高的相关性;气象因子对GPC的重要性依次为辐射量（nRad）、温度（nTem）和降水量（nPre），相对重要性分别为10.05%、7.61%和6.54%;遥感信息重要性最低。但是，当以局部区域（以河北省冬小麦研究区为例）为研究区时，局部区域的品种多样性减少，Glu的相对重要性降低到48.22%;同时由于局部区域的气象条件变异减小，除nRad比重相对增加外，nTem和nPre依次减少;反之，体现作物长势及养分差异的EVI的相对重要性增加（22.65%）。

通过以上分析初步得到以下结论。模型中各个自变量因子对GPC的贡献度在不同的研究尺度有所差异：局部区域品种与气象因子差异较小，植株长势及养分的变异对GPC的结果影响更大;反之，大范围区域品种与气象因子的较大差异，导致植株长势及养分的变异对GPC的结果影响降低。当然，本研究在样点收集过程中通过示范区（2008年和2009年）和农户协助收集（2019年）获得，样点空间代表性不强，对结论的解释还有待进一步验证。另外，本研究气象因子选择4月份数据作为模型输入，主要考虑到该生育阶段包括的拔节期、挑旗期、抽穗期及开花期是冬小麦生长最重要的生育时期，对于群体的建成、穗粒数以及养分积累具有决定作用。对于区域应用而言，由于区域间该生育时期时长差异，后期在进一步品质预测的研究中，需结合区域物候提取，进而精确获取冬小麦关键生育时期的气象数据。

3.3冬小麦GPC预报模型

以3个气象因子为输入变量构建冬小麦GPC预报模型作为初步的区域预报模型。图3是GPC 估算模型建模集和验证集的实测值与估算值关系图。建模集的R为0.39，RMSE为1.04%，可以看出样点数据中高估的样点数据多于被低估的样点数据，尤其是部分2009年北京样点和江苏样点，但大部分样点是集中在1：1线的周围;验证集的R为0.43， RMSE为0.94%;所有样本数据（建模集+验证集）实测值与估算值的R为0.40，RMSE为1.03%。结果表明，结合遥感与气象数据的分层线性预测模型在区域GPC预测上具有可行性，并且模型的预测稳定性较高。

3.4冬小麦GPC区域预报

将GPC模型由样点数据推及到面域范围，得到2019年冬小麦GPC估算图（图4）。越红的地区代表GPC越高，绿色代表GPC偏低，基本呈现从北到南冬小麦品质逐渐降低的趋势。2019 年黄淮海地区冬小麦GPC最高的地区是河北省、山东省北部以及河南省中部和北部，基本在13.6%以上，其中河北省大部分及河南省北部GPC达到14.6%以上;山东省西部、河南省东部与西南部地区、安徽省北部与江苏省北部地区GPC基本在12.3%～13.6%之间。GPC最低的地区是江苏省东南部沿海地区与河南省东南部地区，GPC基本低于12.3%。

4结论

本研究在864个冬小麦实测品质样点（2008年、2009年和2019年）基础上，利用分层线性模型构建了冬小麦GPC估算模型对2019年冬小麦主产区（河南省、山东省、河北省、安徽省和江苏省）范围内的冬小麦GPC进行估算。估算模型分为两层：第一层模型是基于GPC、EVI和冬小麦籽粒筋型构建;第二层模型是基于第一层模型中的模型系数与外部气象数据（温度nTem、降水nPre、辐射量nRad）的归一化值构建。GPC估算模型建模集精度R为0.39，RMSE为1.04%;验证集精度R为0.43，RMSE为0.94%。将GPC估算模型应用于冬小麦主产区的GPC遥感估算，最终得到2019年冬小麦主产区品质预报图。结果证明，构建的籽粒蛋白质遥感预测模型具有较好的普适性和时空扩展性。

本研究以MODIS影像（空间分辨率为250m）为例初步在大尺度黄淮海地区冬小麦种植区开展应用，粗分辨率的遥感影像数据在冬小麦面积提取方面存在混合像元的问题在本研究中未进行深入探讨。因此，在接下来的研究中有必要开展高分辨率影像数据（Landsat-TM. Sentinel-2或者高分系列卫星数据等）在冬小麦品质预报的应用。

参考文献：

[1]王纪华，李存军，刘良云，等.作物品质遥感监测预报研究进展[J].中国农业科学，2008， 41（9）： 2633-2640. WANG J， LI C， LIU L， et al. Estimating grain protein content in winter wheat with multi-temporal. hyperspectral. measurements[J]. Scientia Agricultura Sinica， 2008， 41（9）：2633-2640.

[2]李少昆，谭海珍，王克如，等.小麦籽粒蛋白质含量遥感监测研究进展[J].农业工程学报，2009， 25（2）： 310-315.

LI S， TAN H， WANG K， et al. Research progress in wheat grain protein content monitoring using remote sensing[J]. Transactions of the CSAE， 2009， 25（2）： 310-315.

[3] JIN X， XU X， FENG H， et al. Estimation of grain protein content in winter wheat by using three methods with hyperspectral. data[J]. International. Journal. of Agriculture & Biology， 2014， 16（3）： 498-504.

[4] ONOYAMA H， RYU C， SUGURI M， et al. Estimation of rice protein content before harvest using ground- based hyperspectral. imaging and region of interest analysis[J]. Precision Agriculture， 2017， 19： 721-734.

[5] RODRIGUES F A， BLASCH G， DEFOURNY P， et al. Multi-temporal. and spectral. analysis of high-resolution hyperspectral. airborne imagery for precision agriculture： Assessment of wheat grain yield and grain protein content[J]. Remote Sensing， 2018， 10： ID 930.

[6] XU X， LI C， DONG Y， et al. Estimating grain protein content in winter wheat with multi-temporal. hyperspectral. measurements[J]. Sensor Letters， 2014， 12： 1-5.

[7] BARMEIER G， HOFER K， SCHMIDHALTER U. Mid-season prediction of grain yield and protein content of spring barley cultivars using high-throughput spectral. sensing[J]. European Journal. of Agronomy， 2017， 90： 108-116.

[8]陳鹏飞，王吉顺，潘鹏，等.基于氮素营养指数的冬小麦籽粒蛋白质含量遥感反演[J].农业工程学报，2011，27（9）：75-80.

CHEN P， WANG J， PAN P， et al. Remote detection of wheat grain protein content using nitrogen nutrition index[J]. Transactions of the CSAE， 2011， 27（9）： 75-80

[9]屈莎，李振海，邱春霞，等.基于开花期氮素营养指标的冬小麦籽粒蛋白质含量遥感预测[J].农业工程学报，2017， 33（12）：186-193.

QU S， LI Z， QIU C， et al. Remote sensing prediction of winter wheat grain protein content based on nitrogen nutrition index at anthesis stage[J]. Transactions of the CSAE， 2017，33（12）： 186-193.

[10] ZHAO H， SONG X， YANG G， et al. Monitoring of nitrogen and grain protein content in winter wheat based on sentinel-2A data[J]. Remote Sensing， 2019， 11（14）： ID 1724.

[11]李振海，徐新剛，金秀良，等.基于氮素运转原理和GRA-PLS算法的冬小麦籽粒蛋白质含量遥感预测[J].中国农业科学，2014， 47（19）： 3780-3790.

LI Z， XU X， JIN X， et al. Remote sensing prediction of winter wheat protein content based on nitrogen translocation and GRA-PLS method[J]. Scientia Agricultura Sinica， 2014， 47（19）： 3780-3790.

[12] LI Z， TAYLOR J， YANG H， et al. A hierarchical. interannual. wheat yield and grain protein prediction model using spectral. vegetative indices and meteorological. data[J]. Field Crops Research， 2020， 248： ID 107711.

[13] XU X， TENG C， ZHAO Y， et al. Prediction of wheat grain protein by coupling multisource remote sensing imagery and ECMWF data[J]. Remote Sensing， 2020，12（8）： ID 1349.

[14] LI Z， JIN X， ZHAO C， et al. Estimating wheat yield and quality by coupling the DSSAT-CERES model andproximal. remote sensing[J]. European Journal. of Agronomy， 2015， 71： 53-62.

[15] LI Z， WANG J， XU X， et al. Assimilation of two variables derived from hyperspectral. data into the DSSAT- CERES model for grain yield and quality estima- tion[J]. Remote Sensing， 2015， 7（9）： 12400-12418.

[16]李振海，杨贵军，王纪华，等.作物籽粒蛋白质含量遥感监测预报研究进展[J].中国农业信息，2018， 30（1）： 46-54.

LI Z， YANG G， WANG J， et al. Remote sensing of grain protein content in cereal： A review[J]. China Agricultural. Informatics， 2018， 30（1）： 46-54.

[17] FELDSTAIN A， WOLTMAN H， MACKAY J C， et al. An introduction to hierarchical. linear modeling[J]. Tutorials in Quantitative Methods for Psychology， 2012， 8 （1）： 62-69.

[18] YE D， NG Y K， LIAN Y. Culture and happiness[J]. Social. Indicators Research， 2015， 123（2）： 519-547.

Estimating Grain Protein Content of Winter Wheat in Producing Areas Based on Remote Sensing and Meteorological. Data

WANG Lin1，2， LIANG Jian3， MENG Fanyu4， MENG Yang1，2， ZHANG Yongtao5， LI Zhenhai1，2*

（1. Key Laboratory of Quantitative Remote Sensing in Ministry of Agriculture and Rural. Affairs/Beijing Research Center for Information Technology in Agriculture， Beijing 100097， China; 2. National. Engineering Research Center for Information Technology in Agriculture， Beijing 100097， China; 3. National. Agro-tech Extension and Service Center， Beijing 100125， China; 4. Beijing Agriculture Technology Extension Station， Beijing 100029， China; 5. Jiangsu Nonidt Agricultural. Science and Technology Co. Ltd， Nanjing 210001， China）

Abstract： With the rapid development of economy and people's living standards， people's demands for crops have changed from quantity to quality. The rise and rapid development of remote sensing technology provides an effective method for crop monitoring. Accurately predicting wheat quality before harvest is highly desirable to optimize management for farmers， grading harvest and categorized storage for the enterprise， future trading price， and policy planning. In this research， the main producing areas of winter wheat （Henan， Shandong， Hebei， Anhui and Jiangsu provinces） were chosed as the research areas， with collected 898 samples of winter wheat over growing seasons of 2008， 2009 and 2019. A Hierarchical. Linear model （HLM） for estimating grain protein content （GPC） of winter wheat at heading-flowering stage was constructed to estimate the GPC of winter wheat in 2019 by using meteorological. factors， remote sensing imagery and gluten type of winter wheat， where remote sensing data and gluten type were input variables at the first level of HLM and the meteorological. data was used as the second level of HLM. To solve the problem of deviation in interannual. and spatial. expansion of GPC estimation model， maximum values of Enhanced Vegetation Index （EVI） from April to May calculated by moderate-resolution-imaging spectroradiometer were computed to represent the crop growth status and used in the GPC estimation model. Critical. meteorological. factors （temperature， precipitation，radiation） and their combinations for GPS estimation were compared and the best estimation model was used in this study. The results showed that the accuracy of GPC considering three meteorological. factors performed higher accuracy （Calibrated set： R=0.39， RMSE=1.04%; Verification set： R=0.43， RMSE=0.94%） than the others GPC model with two meteorological. factors or single meteorological. factor. Therefore， three meteorological. factors were used as input variables to build a winter wheat GPC forecast model for the regional. winter wheat GPC forecast in this research. The GPC estimation model was applied to the GPC remote sensing estimation of the main winter wheat-producing areas， and the GPC prediction map of the main winter wheat producing areas in 2019 was obtained， which could obtain the distribution of winter wheat quality in the Huang-Huai-Hai region. The results of this study could provide data support for subsequent wheat planting regionalization to achieve green， high- yield， high-quality and efficient grain production.

Key words： winter wheat; grain protein content （GPC）; remote sensing; hierarchical. linear model （HLM）; meteorological. data

（登陸www.smartag.net.cn免费获取电子版全文）

作者简介：王琳（1995—），女，硕士研究生，研究方向为遥感信息处理与分析。E-mail：17852320332@163.com。

*通讯作者：李振海（1989—），男，博士，高级工程师，研究方向为农业遥感与信息技术。电话：010-51503215。E-mail：lizh323@126.com。