应用粒子群算法的遥感信息与水稻生长模型同化技术

在研究遥感信息和水稻生长模型的同化过程中, 最小化遥感反演与生长模型(RiceGrow)输出的水稻生长 信息差值绝对值时引入了一种新的优化算法-粒子群算法(PSO), 并对比了其与模拟退火算法(SA)的优缺点; 探讨 了叶面积指数(LAI)和叶片氮积累量(LNA)分别作为同化参数时的同化效果。结果表明, PSO 无论是从同化效率还是 反演精度上都要好于SA, 粒子群优化算法是一种可靠的遥感与模型同化算法; LAI 和LNA 作为外部同化参数时各 有优势, LAI 作为同化参数可获得较准确的播期及播种量, 而LNA 作为同化参数可获得更为准确的施氮量信息。但 是LAI 作为外部同化参数时的反演结果总体要优于利用LNA 作为同化参数时的反演结果。利用试验资料对该技术 进行了测试和检验, 结果显示反演的模型初始参数的平均值与真实值的相对误差(RE)均小于2.5%, 均方根误差 (RMSE)为0.7—2.2, 产量模拟值与实测值之间的相对误差为5%左右, 模拟与实测相关指标值吻合度较高, 该同化 技术具有较好的适用性。从而为生长模型从单点扩展到区域尺度应用奠定了基础。     

ASAR数据与水稻作物模型同化制作水稻产量分布图

提出了利用雷达数据进行水稻估产的技术方法,并以ASAR数据为例,探讨了雷达数据在水稻估产中的可行性.首先利用ASAR数据进行水稻制图,从各时相ASAR数据中提取水稻后向散射系数.随后,基于像元尺度,采用同化方法,以LAI为结合点,将水稻作物模型ORYZA2000与半经验水稻后向散射模型结合,建立嵌套模型模拟水稻后向散射系数.选择水稻出苗期和播种密度为参数优化对象,利用全局优化算法SCE-UA对0RYZA2000模型重新初始化,使模拟的水稻后向散射系数值与实测值误差最小,并由优化后的ORYZA2000模型计算每个像元的水稻产量,生成水稻产量分布图.结果表明,水稻产量分布图能够描绘研究区水稻实际产量的分布趋势,但由于采用潜在生长条件模拟,模拟的水稻平均产量比实测平均值高约13%,验证点的水稻产量模拟值与实测值相对误差为11.2%.由于半经验水稻后向散射模型存在对LAI变化不够敏感和对水层的简化处理,增加了水稻估产的误差.但从总体上看,利用该方法进行区域水稻估产是可行的,并为多云多雨地区的水稻遥感监测提供了重要参考.     

基于WOFOST模型与UAV数据的玉米生长后期地上生物量估算

地上生物量能够有效反映作物的生长状态,其信息的实时估算对产量预测和农田生产管理都有重要意义。作物生长模型因其详尽的生理生化基础和对生长过程数字化描述能力,成为生物量估算的理想模型。近年来,研究人员利用数据同化算法将时间序列遥感数据同化到作物生长模型中,实现了作物模型由基于气象站的点模拟到区域尺度面模拟的外推,使生物量模拟结果同时具备大范围和机理性两个方面的特点。这一模式下,时间序列的遥感数据质量将对生物量模拟精度产生直接影响,作物生长后期受到光谱饱和的影响,遥感数据的作物冠层信息获取能力会出现明显下降,因此有必要对该阶段遥感数据和作物模型的结合方式进行优化,提升生物量模拟精度。本文针对东北地区春玉米生物量遥感估算存在的问题,提出了利用WOFOST作物模型结合无人机（UAV）遥感数据实现作物生长后期生物量准确估算的新思路。新思路首先利用多光谱遥感数据获取WOFOST模型具备较高空间异质性的土壤速效养分参数以提升模型的空间信息模拟能力,使其能在一定程度上摆脱点尺度模拟的限制。同时,结合集合卡尔曼滤波算法将生长前期无人机（UAV）遥感数据同化到模型中,以缩短模型单独运行时间,减少模型运行过程中的参数误差累积,实现无遥感数据参与下的短期作物生长模拟,并输出生长后期相应的生物量模拟结果。最后,本文利用地面实测数据对新方法的生物量模拟精度进行了评价。结果表明,与全生育期数据同化相比,新方法的生物量估算精度有了明显的提升（全生育期同化：R² = 0.45,RMSE = 4254.30 kg/ha;新方法：R²= 0.86,RMSE = 2216.79 kg/ha）。     

GF-6影像下的南方路域针叶植被LAI反演研究

为了更好应用国产高分辨率遥感影像监测评价南方路域植被环境,研究南方路域针叶植被叶面积指数遥感反演.该文以长益高速研究区域的高分六号影像(GF-6)为基础,提出了可适用于针叶叶片的LIBERTY+ SAIL耦合模型并结合多元线性回归、局部加权回归反演路域植被针叶LAI的方法.研究中以耦合模型模拟的冠层光谱反射率、GF-6影像和野外实测生化参数为数据源,通过相关性分析,将与LAI相关性较高的SAVI、RVI和EVI 3种植被指数作为反演因子,结合组合模型反演LAI并评定模型的反演精度.结果 表明,耦合模型对南方路域针叶植被LAI的估算精度整体较高,对比分析两种叶面积指数的组合预测模型,耦合模型结合局部加权回归组合反演LAI具有优越性,可更好地反演路域植被针叶LAI.     

基于遥感数据与作物生长模型同化的冬小麦长势监测与估产方法研究

本文以LAI作为结合点，讨论了利用复合型混合演化（SCE—UA）算法实现CERES—Wheat模型与遥感数据同化的可行性。CERES—Wheat模型同化后主要生育期和产量的模拟值分别与真实条件下模型相应模拟值以及实测值进行比较。结果表明，同化后CERES—Wheat模型的模拟精度对LAI外部同化数据的误差并不十分敏感。并且在LAI同化数据较少时，也可获得较好的同化结果。这一特点体现了SCE—UA算法应用于同化过程的优越性，为同化策略在区域冬小麦长势监测及估产中的应用提供了基础。     

利用双变量同化与PCA-Copula法的冬小麦单产估测

为了进一步提高冬小麦产量估测的精度,基于集合卡尔曼滤波算法和粒子滤波（particle filter, PF）算法,对CERES–Wheat模型模拟的冬小麦主要生育期条件植被温度指数（vegetation temperature condition index,VTCI）、叶面积指数（leaf area index, LAI）和中分辨率成像光谱仪（moderate-resolution imaging spectroradiometer, MODIS）数据反演的VTCI、LAI进行同化,利用主成分分析与Copula函数结合的方法构建单变量和双变量的综合长势监测指标,建立冬小麦单产估测模型,并通过对比分析选择最优模型,对2017—2020年关中平原的冬小麦单产进行估测。结果表明,单点尺度的同化VTCI、同化LAI均能综合反映MODIS观测值和模型模拟值的变化特征,且PF算法具有更好的同化效果;区域尺度下利用PF算法得到的同化VTCI和LAI所构建的双变量估产模型精度最高,与未同化VTCI和LAI构建的估产模型精度相比,研究区各县（区）的冬小麦估测单产与实际单产的均方根误差降低了56.25 kg/hm2,平均相对误差降低了1.51%,表明该模型能有效提高产量估测的精度,应用该模型进行大范围的冬小麦产量估测具有较好的适用性。     

GEE环境下的玉米低温冷害损失快速评估

大范围、及时、准确的灾害损失评估与制图对防灾减灾、农业保险和粮食安全等至关重要。针对传统灾害损失评估方法空间尺度单一、泛化能力差、时效性低,可操作性弱等问题,本文建立了一种遥感产品耦合作物模型的多尺度的灾害损失评估方法MDLA （a Multiscale Disaster Loss Assessment）。该方法利用作物模型的多情景模拟产生大量的灾害样本,结合对应日期的遥感指标构建灾害脆弱性模型,依托Google Earth Engine（GEE）平台将其应用到高分辨率遥感影像和格点灾害指标进行逐象元评估。以鄂伦春自治旗玉米为例,基于精细校准的CERES-Maize模型的模拟,利用两个生长季窗口的LAI和冷积温（CDD）建立统计模型来刻画低温对最终产量的影响,结合Sentinel-2数据逐格点计算完成高精度损失制图。结果显示,校准后的CERES-Maize模拟物候和产量的NRMSE 分别为3.3%和8.9%。冷害情景模拟结果表明不同类型和生育期的低温冷害对玉米产量的影响不尽相同,其中生长峰值期（出苗—吐丝和吐丝—灌浆）最为敏感。回代检验显示,MDLA方法估算精度为11.4%,与历史冷害年份的实际损失相吻合。经评估,鄂伦春2018-08-09的冷害导致玉米减产23.7%,受灾面积1.86×10⁴ ha,其中高海拔地区损失较重（减产率>25%）,低温冷害对该区玉米生产构成了严重的威胁。与现有的统计回归、作物模型模拟以及同化等技术相比,其优势在于：（1）结合遥感观测和作物模型模拟技术能更好地刻画了灾害对产量的影响过程;（2）利用GEE平台快速处理海量遥感数据,提高了灾害损失评估的时效性;（3）不受地面实测数据的限制,易操作,可实现动态、多尺度（象元、田块、村,县等）的损失评估,这为防灾减损、维持粮食丰产稳产提供了保障,也为农业保险的业务化运行提供了思路。     

ENVISAT ASAR 数据用于水稻监测和参数反演

用雷达后向散射模型模拟了水稻生长周期内入射角对雷达后向散射的影响关系。用模拟结果归一化雷达数据的后向散射系数，得到同一入射角下水稻周期内后向散射系数时间序列值。分析了归一化ASAR数据与水稻生物参数的关系，实验结果表明，ASAR数据可以用来估测水稻参数。     

基于数据同化的元胞自动机

提出基于集合卡尔曼滤波(EnKF)的元胞自动机(CA)模型。在CA模型中,由于不同的样本会训练出不同参数值 的转换规则,且获取的转换规则在整个模拟过程中不能改变等原因,误差在模拟过程中会不断累积。本文在CA模型中 引入集合卡尔曼滤波的数据同化方法,建立了基于集合卡尔曼滤波的数据同化CA模型,同化遥感观测数据,根据得出 的同化值修正模拟结果使之向真实情况逼近。利用该模型模拟了广东省东莞市的发展情景(1995年—2005年),实验表 明,与传统CA模型相比,基于集合卡尔曼滤波的CA模型能够融合遥感观测数据,并能更有效地模拟城市扩张过程,达 到良好的模拟效果。     

农作物长势综合遥感监测方法

作物收获之前进行大范围作物生长状况评价 ,可以尽早的获得有关作物产量信息。介绍了中国农情遥感监测系统的综合作物长势监测方法。以遥感数据标准化处理、云标识、云污染去除和非耕地去除为基础 ,生成质量一致的遥感数据产品集 ,提取区域作物生长过程。作物长势监测分为实时作物长势监测和作物生长趋势分析。实时的作物长势监测可以定性和定量地在空间上分析作物生长状况 ,分级显示作物生长状况 ,分区域统计水田和旱地中不同长势占的比重。作物生长趋势分析可以进行年际间的生长过程对比 ,从时间轴上反映作物持续生长的差异性 ,统计全国、主产区、省和区划单元 4个尺度的耕地、水田、旱地作物生长过程曲线年际间差异 ,从而为早期的产量预测提供信息。通过处理流程的系统化 ,建设了运行化的作物长势遥感监测分析系统 ,为用户构建了综合的作物实时生长状况 ,苗情的生长趋势分析环境。同时可以依据野外地面实测信息对遥感监测结果进行标定和检验。 1998年以来 ,系统在满足日常运行的前提下 ,技术方法逐渐改进和完善 ,监测结果的精度和可靠性不断得到提高。     

Crop growth modelling and crop yield forecasting using satellite-derived meteorological inputs

Distributed crop simulation models are typically confronted with considerable uncertainty in weather variables. In this paper the use of MeteoSat-derived meteorological products to replace weather variables interpolated from weather stations (temperature, reference evapotranspiration and radiation) is explored. Simulations for winter-wheat were carried for Spain, Poland and Belgium using both interpolated and MeteoSat-derived weather variables. The results were spatially aggregated to national and regional level and were evaluated by comparing the simulation results of both approaches and by assessing the relationships with crop yield statistics over the periods 1995–2003 from EUROSTAT. The results indicate that potential crop yield can be simulated well using MeteoSat-derived meteorological variables, but that water-stress hardly occurs in the water-limited simulations. This is caused by a difference in reference evapotranspiration which was 20–30% smaller in case of MeteoSat. As a result, the simulations using MeteoSat-derived meteorological variables performed considerably poorer in a regression analyses with crop yield statistics on national and regional level. Our results indicate that a recalibration of the model parameters is necessary before the MeteoSat-derived meteorological variables can be used operationally in the system.     

Remote sensing-based global crop monitoring: experiences with China's CropWatch system

Monitoring the production of main agricultural crops is important to predict and prepare for disruptions in food supply and fluctuations in global crop market prices. China's global crop-monitoring system (CropWatch) uses remote sensing data combined with selected field data to determine key crop production indicators: crop acreage, yield and production, crop condition, cropping intensity, crop-planting proportion, total food availability, and the status and severity of droughts. Results are combined to analyze the balance between supply and demand for various food crops and if needed provide early warning about possible food shortages. CropWatch data processing is highly automated and the resulting products provide new kinds of inputs for food security assessments. This paper presents a comprehensive overview of CropWatch as a remote sensing-based system, describing its structure, components, and monitoring approaches. The paper also presents examples of monitoring results and discusses the strengths and limitations of the CropWatch approach, as well as a comparison with other global crop-monitoring systems.     

Monitoring winter wheat growth in North China by combining a crop model and remote sensing data

Both of crop growth simulation models and remote sensing method have a high potential in crop growth monitoring and yield prediction. However, crop models have limitations in regional application and remote sensing in describing the growth process. Therefore, many researchers try to combine those two approaches for estimating the regional crop yields. In this paper, the WOFOST model was adjusted and regionalized for winter wheat in North China and coupled through the LAI to the SAIL–PROSPECT model in order to simulate soil adjusted vegetation index (SAVI). Using the optimization software (FSEOPT), the crop model was then re-initialized by minimizing the differences between simulated and synthesized SAVI from remote sensing data to monitor winter wheat growth at the potential production level. Initial conditions, which strongly impact phenological development and growth, and which are hardly known at the regional scale (such as emergence date or biomass at turn-green stage), were chosen to be re-initialized. It was shown that re-initializing emergence date by using remote sensing data brought simulated anthesis and maturity date closer to measured values than without remote sensing data. Also the re-initialization of regional biomass weight at turn-green stage led that the spatial distribution of simulated weight of storage organ was more consistent to official yields. This approach has some potential to aid in scaling local simulation of crop phenological development and growth to the regional scale but requires further validation.     

Integration of maximum crop response with machine learning regression model to timely estimate crop yield

Timely and reliable estimation of regional crop yield is a vital component of food security assessment, especially in developing regions. The traditional crop forecasting methods need ample time and labor to collect and process field data to release official yield reports. Satellite remote sensing data is considered a cost-effective and accurate way of predicting crop yield at pixel-level. In this study, maximum Enhanced Vegetation Index (EVI) during the crop-growing season was integrated with Machine Learning Regression (MLR) models to estimate wheat and rice yields in Pakistan's Punjab province. Five MLR models were compared using a fivefold cross-validation method for their predictive accuracy. The study results revealed that the regression model based on the Gaussian process outperformed over other models. The best performing model attained coefficient of determination (R2), Root Mean Square Error (RMSE, t/ ha), and Mean Absolute Error (MAE, t/ha) of 0.75, 0.281, and 0.236 for wheat; 0.68, 0.112, and 0.091 for rice, respectively. The proposed method made it feasible to predict wheat and rice 6– 8 weeks before the harvest. The early prediction of crop yield and its spatial distribution in the region can help formulate efficient agricultural policies for sustainable social, environmental, and economic progress.     

基于作物缺水指数的土壤含水量估算方法

为研究江苏省徐州市的土壤水分时空分布及动态变化,基于MODIS数据和站点气象数据,利用蒸散发双层模型和考虑土壤水分可供率的改进双层模型分别计算实际蒸散发量,利用Penman-Monteith模型计算区域潜在蒸散发量,计算获得作物缺水指数(crop water stress index,CWSI),并与2010年7月和11月的土壤相对含水量实测数据分别进行回归分析建模,得到了土壤含水量分布图。结果表明:基于蒸散发双层模型的土壤含水量估算结果与实测值的决定系数分别为0.53和0.72,平均相对误差分别为5.89%和9.6%;对双层模型进行改进后,土壤含水量估算结果与实测值的决定系数都为0.84,平均相对误差分别为3.47%和6.03%,利用改进后的双层模型对土壤相对含水量进行估算效果更好。     

机载Lidar数据的农作物覆盖度及LAI反演

虽然Lidar点云数据已被广泛应用于获取森林各项结构参数,但这些方法并不适合于低矮的灌丛、林地和农作物。本文以玉米为研究对象,提出利用机载Lidar点云数据的强度信息和全波形数据中的距离与扫描天顶角信息,反演农作物覆盖度和LAI的方法。在黑河进行的飞行实验和地面验证表明,该方法具有较高精度,也表明Lidar在低矮自然植被监测和农业应用上有较大潜力。     

Assessing the benefit of satellite-based Solar-Induced Chlorophyll Fluorescence in crop yield prediction

Large-scale crop yield prediction is critical for early warning of food insecurity, agricultural supply chain management, and economic market. Satellite-based Solar-Induced Chlorophyll Fluorescence (SIF) products have revealed hot spots of photosynthesis over global croplands, such as in the U.S. Midwest. However, to what extent these satellite-based SIF products can enhance the performance of crop yield prediction when benchmarking against other existing satellite data remains unclear. Here we assessed the benefits of using three satellite-based SIF products in yield prediction for maize and soybean in the U.S. Midwest: gap-filled SIF from Orbiting Carbon Observatory 2 (OCO-2), new SIF retrievals from the TROPOspheric Monitoring Instrument (TROPOMI), and the coarse-resolution SIF retrievals from the Global Ozone Monitoring Experiment-2 (GOME-2). The yield prediction performances of using SIF data were benchmarked with those using satellite-based vegetation indices (VIs), including normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near-infrared reflectance of vegetation (NIRv), and land surface temperature (LST). Five machine-learning algorithms were used to build yield prediction models with both remote-sensing-only and climate-remote-sensing-combined variables. We found that high-resolution SIF products from OCO-2 and TROPOMI outperformed coarse-resolution GOME-2 SIF product in crop yield prediction. Using high-resolution SIF products gave the best forward predictions for both maize and soybean yields in 2018, indicating the great potential of using satellite-based high-resolution SIF products for crop yield prediction. However, using currently available high-resolution SIF products did not guarantee consistently better yield prediction performances than using other satellite-based remote sensing variables in all the evaluated cases. The relative performances of using different remote sensing variables in yield prediction depended on crop types (maize or soybean), out-of-sample testing methods (five-fold-cross-validation or forward), and record length of training data. We also found that using NIRv could generally lead to better yield prediction performance than using NDVI, EVI, or LST, and using NIRv could achieve similar or even better yield prediction performance than using OCO-2 or TROPOMI SIF products. We concluded that satellite-based SIF products could be beneficial in crop yield prediction with more high-resolution and good-quality SIF products accumulated in the future.     

Images difference of ASAR data for rice crop mapping in Fuzhou,China

This paper introduces ENVISAT ASAR data application on rice field mapping in the Fuzhou area, using multi-temporal ASAR dual polarization data acquired in 2005. The procedure for ASAR data processing here includes data calibration, image registration, speckle reduction and conversion of data format from amplitude to dB for backscatter. The backscatter of rice increases with the rice growing stages, which was much different from other land covers. Based on image difference techniques, 6 schemes were designed with ASAR different temporal and polarization data for rice field mapping. Difference images between images in the early period of rice crop and growing or ripening period, are more suitable for rice extraction than those difference images between different polarizations in the same date. The most accurate result of late rice extraction was achieved based on the difference of HH polarization data acquired in October and August. Therefore, for rice field mapping, the temporal information is more important than polarization information. The data during the early growing season of rice is very important for high accuracy rice mapping.     

农作物长势综合监测——以印度为例

提出农作物长势综合监测方法,利用卫星遥感得到的NDVI时间序列数据,综合采用实时监测、过程监测和时间序列聚类监测方法,明确不同方法适用的监测尺度及监测目的,对不同范围农作物长势进行监测。改进了Crop Watch全球农情遥感速报系统运行化作物长势监测方法,克服了原有作物长势监测中实时监测方法无法反映相同区域苗情在整个生长过程中的连续变化情况的缺点。实现对相同区域作物长势连续变化的定量描述,可对作物长势进行更准确的判断。利用官方发布的作物单产变幅数据,对单产变幅较大的12个作物主产省区作物长势监测结果的准确性进行判断,结果表明:6个邦的实时监测和聚类监测方法所得结果一致,都符合作物单产变化的实际状况;4个邦的聚类监测方法所得结果对作物长势监测更为准确,更符合该区域作物单产的实际变化;1个邦实时监测结果对作物长势监测比聚类监测方法更为准确;只有1个邦采用两种方法对作物长势的监测存在误差,聚类监测方法在对农作物生长过程的连续监测及空间分布的定量化表述方面,比实时监测更为准确。3种方法可以综合使用,实现业务化运行的农作物长势监测。