Retrieving rice (Oryza sativa L.) net photosynthetic rate from UAV multispectral images based on machine learning methods

Tianao Wu; Wei Zhang; Shuyu Wu; Minghan Cheng; Lushang Qi; Guangcheng Shao; Xiyun Jiao

doi:10.3389/fpls.2022.1088499

Retrieving rice (Oryza sativa L.) net photosynthetic rate from UAV multispectral images based on machine learning methods

Front Plant Sci. 2023 Jan 25:13:1088499. doi: 10.3389/fpls.2022.1088499. eCollection 2022.

Authors

Tianao Wu^{1

2

3}, Wei Zhang¹, Shuyu Wu^{1

2

3}, Minghan Cheng⁴, Lushang Qi¹, Guangcheng Shao¹, Xiyun Jiao^{1

2

3}

Affiliations

¹ College of Agricultural Science and Engineering, Hohai University, Nanjing, China.
² State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing, China.
³ Cooperative Innovation Center for Water Safety and Hydro Science, Hohai University, Nanjing, China.
⁴ Jiangsu Key Laboratory of Crop Genetics and Physiology/Jiangsu Key Laboratory of Crop Cultivation and Physiology, Agricultural College, Yangzhou University, Yangzhou, China.

Abstract

Photosynthesis is the key physiological activity in the process of crop growth and plays an irreplaceable role in carbon assimilation and yield formation. This study extracted rice (Oryza sativa L.) canopy reflectance based on the UAV multispectral images and analyzed the correlation between 25 vegetation indices (VIs), three textural indices (TIs), and net photosynthetic rate (Pn) at different growth stages. Linear regression (LR), support vector regression (SVR), gradient boosting decision tree (GBDT), random forest (RF), and multilayer perceptron neural network (MLP) models were employed for Pn estimation, and the modeling accuracy was compared under the input condition of VIs, VIs combined with TIs, and fusion of VIs and TIs with plant height (PH) and SPAD. The results showed that VIs and TIs generally had the relatively best correlation with Pn at the jointing-booting stage and the number of VIs with significant correlation (p< 0.05) was the largest. Therefore, the employed models could achieve the highest overall accuracy [coefficient of determination (R ²) of 0.383-0.938]. However, as the growth stage progressed, the correlation gradually weakened and resulted in accuracy decrease (R ² of 0.258-0.928 and 0.125-0.863 at the heading-flowering and ripening stages, respectively). Among the tested models, GBDT and RF models could attain the best performance based on only VIs input (with R ² ranging from 0.863 to 0.938 and from 0.815 to 0.872, respectively). Furthermore, the fusion input of VIs, TIs with PH, and SPAD could more effectively improve the model accuracy (R ² increased by 0.049-0.249, 0.063-0.470, and 0.113-0.471, respectively, for three growth stages) compared with the input combination of VIs and TIs (R ² increased by 0.015-0.090, 0.001-0.139, and 0.023-0.114). Therefore, the GBDT and RF model with fused input could be highly recommended for rice Pn estimation and the methods could also provide reference for Pn monitoring and further yield prediction at field scale.

Keywords: UAV multispectral remote sensing; machine learning; net photosynthetic rate; rice canopy; textural index; vegetation index.

Grants and funding

This study is financially supported by the Jiangsu Provincial Key Research and Development Program (BE2022390) and the Postgraduate Research and Practice Innovation Program of Jiangsu Province (KYCX21_0537).