Abstract:SPAD is a key indicator for evaluating the growth potential and nitrogen nutrition status of millet. To achieve high-precision and wide-coverage dynamic monitoring of the relative chlorophyll content of millet during key growth stages (jointing stage, heading stage, filling stage, and maturity stage), taking the field scale as the research area and integrated multi-source remote sensing data from unmanned aerial vehicles (UAV) and satellites, an SPAD inversion model for optimizing long short-term memory neural networks (LSTM) based on the improved grey wolf optimization algorithm (IGWO) was proposed. During the research process, UAV remote sensing data, satellite remote sensing data and ground point-like SPAD measured values of each key growth period were obtained simultaneously. Through the longitude and latitude coordinates of the ground measured sample points, the point-like SPAD values were spatially matched and associated with the UAV image pixels and satellite image pixels at the corresponding positions to construct an image-measured corresponding data set. The point spread function (PSF) method was adopted for scale upward inference of UAV images. The mean-variance method was combined to preliminarily correct the satellite data. Then the support vector regression (SVR) was used to establish a collaborative correction model for multi-source remote sensing data of UAV and satellite. Based on the corrected high-precision satellite data, the spectral characteristic parameters sensitive to SPAD were screened through Pearson correlation analysis and XG-Boost feature importance ranking. The nonlinear convergence factor was introduced to enhance the hyperparameter optimization ability of the grey wolf optimization algorithm, and finally the IGWO-LSTM SPAD inversion model was constructed. The results showed that compared with other resampling methods, the point spread function method had the least information loss during the scale-up process. The average value and standard deviation of the processed image pixels were 0.103 and 0.056, respectively. The satellite remote sensing data corrected by the histogram matching method effectively retained the original spectral shape, and the spectral angle mapping value (SAM) was as low as 0.062°. The SVR algorithm had the highest model accuracy in the B/G/R/NIR bands during the critical growth period, and the determination coefficients of the four bands were 0.920, 0.961, 0.963 and 0.900, respectively. The coefficient of determination (R2) of the IGWO-LSTM model in the inversion of SPAD during the critical growth period reached 0.985, and the root mean square error (RMSE) was 0.111, which was significantly superior to that of traditional models such as partial least squares regression (PLSR), BP neural network (BPNN), and random forest regression (RF). The research achieved the precise dynamic inversion of SPAD during the key growth period of millet, which was of great significance for the intelligent monitoring of crop growth and the precise application of nitrogen fertilizer.