农业机械学报

作者：马战林，文枫，周颖杰，鲁春阳，薛华柱，李长春

单位：河南城建学院,河南城建学院,河南城建学院,河南城建学院,河南理工大学,河南理工大学,

关键词：冬小麦;估产;Sentinel-2卫星;集合卡尔曼滤波;PyWOFOT作物生长模型

基金：国家自然科学基金面上项目（41721333）、河南省哲学社会科学规划项目(2022BJJ026、2021BJJ062）和河南城建学院大学生创新创业训练计划专项（202211765012、202211765018）

为精准、高效、实时地实现区域冬小麦产量估算，以河南省鹤壁市淇县桥盟乡石桥村为研究区，基于分辨率10m的Sentinel-2多时相光学遥感影像，利用集合卡尔曼滤波（Ensemble Kalman filter，EnKF）算法同化PROSAIL辐射传输模型反演的多期叶面积指数（Leaf area index，LAI）到PyWOFOST作物生长模型中实现一定数量不同长势单点产量的估测，最后利用建立的机器学习模型和面域数据反演区域冬小麦产量，实现作物生长模型与机器学习算法的应用耦合及一种新的区域冬小麦估产模式。研究基于Sobol参数敏感性分析法量化对贮藏器官总干重质量（Total dry weight of storage organs，TWSO）与LAImax的敏感性参数，并基于反演的多期LAI和粒子群优化（Particle swarm optimization，PSO）算法优化与LAImax相关的TDWI、TBASE、CVS、CVL敏感性参数，将其输入到PyWOFOST模型中，利用EnKF算法和时序LAI数据调整对TWSO相关的AMAXTB1、TDWI、TSUMEM、CVO敏感性参数，实现单点产量的估算；与实测单点产量相比，该方法估算的R2、RMSE、MAE、Bias分别为0.8665、468.64kg/hm2、385.70kg/hm2和103.08，为建立随机森林回归（Random forest regression，RFR）区域估产算法提供准确的单点产量训练数据。针对研究区（309.32hm2），基于不同长势人工样点产量数据建立的RFR区域估产算法，区域估产精度为99.44%，每公顷算法运行用时1.55s；应用EnKF算法同化多时期面域LAI到PyWOFOST作物生长模型中的区域估产精度为89.01%，每公顷算法运行用时约0.47h；耦合PyWOFOST作物生长模型与RFR机器学习算法的区域估产精度达到95.58%，每公顷算法运行用时8.85s（训练数据的单点产量计算占总时长约81.35%），显著降低机器学习算法所需的人工成本和同化变量过程计算的时间及算力成本。研究结果为准确、快速的大区域作物估产提供理论支持和技术参考。

来源：2023年第6期

《农业机械学报》期刊编辑部

查看农业机械学报杂志2023年第6期