Evaluation of the Influence of Field Conditions on Aerial Multispectral Images and Vegetation Indices

Remote sensing is a method used for monitoring and measuring agricultural crop fields. Unmanned aerial vehicles (UAV) are used to effectively monitor crops via different camera technologies. Even though aerial imaging can be considered a rather straightforward process, more focus should be given to data quality and processing. This research focuses on evaluating the influences of field conditions on raw data quality and commonly used vegetation indices. The aerial images were taken with a custom-built UAV by using a multispectral camera at four different times of the day and during multiple times of the season. Measurements were carried out in the summer seasons of 2019 and 2020. The imaging data were processed with different software to calculate vegetation indices for 10 reference areas inside the fields. The results clearly show that NDVI (normalized difference vegetation index) was the least affected vegetation index by the field conditions. The coefficient of variation (CV) was determined to evaluate the variations in vegetation index values within a day. Vegetation index TVI (transformed vegetation index) and NDVI had coefficient of variation values under 5%, whereas with GNDVI (green normalized difference vegetation index), the value was under 10%. Overall, the vegetation indices that include near-infrared (NIR) bands are less affected by field condition changes

Evaluation of the Influence of Field Conditions on Aerial Multispectral Images and Vegetation Indices

Remote sensing is a method used for monitoring and measuring agricultural crop fields. Unmanned aerial vehicles (UAV) are used to effectively monitor crops via different camera technologies. Even though aerial imaging can be considered a rather straightforward process, more focus should be given to data quality and processing. This research focuses on evaluating the influences of field conditions on raw data quality and commonly used vegetation indices. The aerial images were taken with a custom-built UAV by using a multispectral camera at four different times of the day and during multiple times of the season. Measurements were carried out in the summer seasons of 2019 and 2020. The imaging data were processed with different software to calculate vegetation indices for 10 reference areas inside the fields. The results clearly show that NDVI (normalized difference vegetation index) was the least affected vegetation index by the field conditions. The coefficient of variation (CV) was determined to evaluate the variations in vegetation index values within a day. Vegetation index TVI (transformed vegetation index) and NDVI had coefficient of variation values under 5%, whereas with GNDVI (green normalized difference vegetation index), the value was under 10%. Overall, the vegetation indices that include near-infrared (NIR) bands are less affected by field condition changes

Ground Truth Validation of Sentinel-2 Data Using Mobile Wireless Ad Hoc Sensor Networks (MWSN) in Vegetation Stands

Satellite-based remote sensing (RS) data are increasingly used to map and monitor local, regional, and global environmental phenomena and processes. Although the availability of RS data has improved significantly, especially in recent years, operational applications to derive value-added information products are still limited by close-range validation and verification deficits. This is mainly due to the gap between standardized and sufficiently available close-range and RS data in type, quality, and quantity. However, to ensure the best possible linkage of close-range and RS data, it makes sense to simultaneously record close-range data in addition to the availability of environmental models. This critical gap is filled by the presented mobile wireless ad hoc sensor network (MWSN) concept, which records sufficient close-range data automatically and in a standardized way, even at local and regional levels. This paper presents a field study conducted as part of the Durable Environmental Multidisciplinary Monitoring Information Network (DEMMIN), focusing on the information gained with respect to estimating the vegetation state with the help of multispectral data by simultaneous observation of an MWSN during a Sentinel-2A (S2A) overflight. Based on a cross-calibration of the two systems, a comparable spectral characteristic of the data sets could be achieved. Building upon this, an analysis of the data regarding the influence of solar altitude, test side topography and land cover, and sub-pixel heterogeneity was accomplished. In particular, variations due to spatial heterogeneity and dynamics in the diurnal cycle show to what extent such complementary measurement systems can improve the data from RS products concerning the vegetation type and atmospheric conditions

Hyperspectral reflectance measurements from UAS under intermittent clouds: Correcting irradiance measurements for sensor tilt

One great advantage of optical hyperspectral remote sensing from unmanned aerial systems (UAS) compared to satellite missions is the possibility to fly and collect data below clouds. The most typical scenario is flying below intermittent clouds and under turbulent conditions, which causes tilting of the platform. This study aims to advance hyperspectral imaging from UAS in most weather conditions by addressing two challenges: (i) the radiometric and spectral calibrations of miniaturized hyperspectral sensors; and (ii) tilting effects on measured downwelling irradiance. We developed a novel method to correct the downwelling irradiance data for tilting effects. It uses a hybrid approach of minimizing measured irradiance variations for constant irradiance periods and spectral unmixing, to calculate the spectral diffuse irradiance fraction for all irradiance measurements within a flight. It only requires the platform's attitude data and a standard incoming light sensor. We demonstrated the method at the Palo Verde National Park wetlands in Costa Rica, a highly biodiverse area. Our results showed that the downwelling irradiance correction method reduced systematic shifts caused by a change in flight direction of the UAS, by 87% and achieving a deviation of 2.78% relative to a on ground reference in terms of broadband irradiance. High frequency (< 3 s) irradiance variations caused by high-frequency tilting movements of the UAS were reduced by up to 71%. Our complete spectral and radiometric calibration and irradiance correction can significantly remove typical striped illumination artifacts in the surface reflectance-factor map product. The possibility of collecting precise hyperspectral reflectance-factor data from UAS under varying cloud cover makes it more operational for environmental monitoring or precision agriculture applications, being an important step in advancing hyperspectral imaging from UAS.Innovation Fund Denmark/[7048-00001B]/IFD/DinamarcaAgricultural Water Innovations in the Tropics/[]/AgWIT/CanadáUniversidad de Costa Rica/[805-C0-603]/UCR/Costa RicaUCR::Vicerrectoría de Docencia::Ciencias Básicas::Facultad de Ciencias::Escuela de Físic

Terrestrial plant productivity and soil moisture constraints

Dolman, A.J. [Promotor]Jeu, R.M.H. de [Copromotor]Werf-, G.R. van der [Copromotor

Illumination Geometry and Flying Height Influence Surface Reflectance and NDVI Derived from Multispectral UAS Imagery

Small unmanned aerial systems (UAS) have allowed the mapping of vegetation at very high spatial resolution, but a lack of standardisation has led to uncertainties regarding data quality. For reflectance measurements and vegetation indices (Vis) to be comparable between sites and over time, careful flight planning and robust radiometric calibration procedures are required. Two sources of uncertainty that have received little attention until recently are illumination geometry and the effect of flying height. This study developed methods to quantify and visualise these effects in imagery from the Parrot Sequoia, a UAV-mounted multispectral sensor. Change in illumination geometry over one day (14 May 2018) had visible effects on both individual images and orthomosaics. Average near-infrared (NIR) reflectance and NDVI in regions of interest were slightly lower around solar noon, and the contrast between shadowed and well-illuminated areas increased over the day in all multispectral bands. Per-pixel differences in NDVI maps were spatially variable, and much larger than average differences in some areas. Results relating to flying height were inconclusive, though small increases in NIR reflectance with height were observed over a black sailcloth tarp. These results underline the need to consider illumination geometry when carrying out UAS vegetation survey

Landsat 5 Thematic Mapper Reflectance and NDVI 27-year Time Series Inconsistencies Due to Satellite Orbit Change

The Landsat 5 Thematic Mapper (TM) sensor provided the longest single mission terrestrial remote sensing data record but temporally sparse station keeping maneuvers meant that the Landsat 5 orbit changed over the 27 year mission life. Long-term Landsat 5 TM reflectance inconsistencies may be introduced by orbit change induced solar zenith variations combined with surface reflectance anisotropy, commonly described by the Bi-directional Reflectance Distribution Function (BRDF). This study quantifies the local overpass time and observed solar zenith angle changes for all the Landsat 5 TM images available at two latitudinally separated locations along the same north-south Landsat path (27) in Minnesota (row 26) and Texas (row 42). Over the 27 years the Landsat 5 orbit changed by nearly 1 h and resulted in changes in the Landsat 5 observed solar zenith angle of N10°. The Landsat 5 orbit was relatively stable from 1984 to 1994 and from 2007 to 2011, but changed rapidly from 1995 to 2000, and from 2003 to 2007. Rather than directly examine Landsat 5 TM reflectance time series a modelling approach was used. This was necessary because unambiguous separation of orbit change induced Landsat reflectance variations from other temporal variations is non-trivial. The impact of Landsat 5 orbit induced observed solar zenith angle variations on the red and near-infrared reflectance and derived normalized difference vegetation index (NDVI) values were modelled with respect to different Moderate-Resolution Imaging Spectroradiometer (MODIS) BRDF land cover types. Synthetic nadir BRDF-adjusted reflectance (NBAR) for the Landsat 5 TM observed and a modelled reference year 2011 solar zenith were compared over the 27 years of acquisitions. Ordinary least squares linear regression fits of the NBAR difference values as a function of the acquisition date indicated an increasing trend in red and near-infrared NBAR and a decreasing trend in NDVI NBAR due to orbit changes. The trends are statistically significant but small, no more than 0.0006 NDVI/year. Comparison of certain years of Landsat 5 data may result in significant reflectance and NDVI differences due only to Landsat 5 orbit changes and cause spurious detection of “browning” vegetation events and underestimation of greening trends. The greatest differences will occur when 1995 Landsat 5 TM data are compared with 2007 to 2011 data; NDVI values could be up to 0.11 greater in 1995 than in 2011 for anisotropic land cover types and up to 0.05 greater for average CONUS land cover types. A smaller number of Landsat 5 TM images were also examined and provide support for the modelled based findings. The paper concludes with a discussion of the implications of the research findings for Landsat 5 TM time series analyses

Assimilation de données satellitaires pour le suivi des ressources en eau dans la zone Euro-Méditerranée

Une estimation plus précise de l'état des variables des surfaces terrestres est requise afin d'améliorer notre capacité à comprendre, suivre et prévoir le cycle hydrologique terrestre dans diverses régions du monde. En particulier, les zones méditerranéennes sont souvent caractérisées par un déficit en eau du sol affectant la croissance de la végétation. Les dernières simulations du GIEC (Groupe d'Experts Intergouvernemental sur l'Evolution du Climat) indiquent qu'une augmentation de la fréquence des sécheresses et des vagues de chaleur dans la région Euro-Méditerranée est probable. Il est donc crucial d'améliorer les outils et l'utilisation des observations permettant de caractériser la dynamique des processus des surfaces terrestres de cette région. Les modèles des surfaces terrestres ou LSMs (Land Surface Models) ont été développés dans le but de représenter ces processus à diverses échelles spatiales. Ils sont habituellement forçés par des données horaires de variables atmosphériques en point de grille, telles que la température et l'humidité de l'air, le rayonnement solaire et les précipitations. Alors que les LSMs sont des outils efficaces pour suivre de façon continue les conditions de surface, ils présentent encore des défauts provoqués par les erreurs dans les données de forçages, dans les valeurs des paramètres du modèle, par l'absence de représentation de certains processus, et par la mauvaise représentation des processus dans certaines régions et certaines saisons. Il est aussi possible de suivre les conditions de surface depuis l'espace et la modélisation des variables des surfaces terrestres peut être améliorée grâce à l'intégration dynamique de ces observations dans les LSMs. La télédétection spatiale micro-ondes à basse fréquence est particulièrement utile dans le contexte du suivi de ces variables à l'échelle globale ou continentale. Elle a l'avantage de pouvoir fournir des observations par tout-temps, de jour comme de nuit. Plusieurs produits utiles pour le suivi de la végétation et du cycle hydrologique sont déjà disponibles. Ils sont issus de radars en bande C tels que ASCAT (Advanced Scatterometer) ou Sentinel-1. L'assimilation de ces données dans un LSM permet leur intégration de façon cohérente avec la représentation des processus. Les résultats obtenus à partir de l'intégration de données satellitaires fournissent une estimation de l'état des variables des surfaces terrestres qui sont généralement de meilleure qualité que les simulations sans assimilation de données et que les données satellitaires elles-mêmes. L'objectif principal de ce travail de thèse a été d'améliorer la représentation des variables des surfaces terrestres reliées aux cycles de l'eau et du carbone dans le modèle ISBA grâce à l'assimilation d'observations de rétrodiffusion radar (sigma°) provenant de l'instrument ASCAT. Un opérateur d'observation capable de représenter les sigma° ASCAT à partir de variables simulées par le modèle ISBA a été développé. Une version du WCM (water cloud model) a été mise en œuvre avec succès sur la zone Euro-Méditerranée. Les valeurs simulées ont été comparées avec les observations satellitaires. Une quantification plus détaillée de l'impact de divers facteurs sur le signal a été faite sur le sud-ouest de la France. L'étude de l'impact de la tempête Klaus sur la forêt des Landes a montré que le WCM est capable de représenter un changement brutal de biomasse de la végétation. Le WCM est peu efficace sur les zones karstiques et sur les surfaces agricoles produisant du blé. Dans ce dernier cas, le problème semble provenir d'un décalage temporel entre l'épaisseur optique micro-ondes de la végétation et l'indice de surface foliaire de la végétation. Enfin, l'assimilation directe des sigma° ASCAT a été évaluée sur le sud-ouest de la France.More accurate estimates of land surface conditions are important for enhancing our ability to understand, monitor, and predict key variables of the terrestrial water cycle in various parts of the globe. In particular, the Mediterranean area is frequently characterized by a marked impact of the soil water deficit on vegetation growth. The latest IPCC (Intergovernmental Panel on Climate Change) simulations indicate that occurrence of droughts and warm spells in the Euro-Mediterranean region are likely to increase. It is therefore crucial to improve the ways of understanding, observing and simulating the dynamics of the land surface processes in the Euro-Mediterranean region. Land surface models (LSMs) have been developed for the purpose of representing the land surface processes at various spatial scales. They are usually forced by hourly gridded atmospheric variables such as air temperature, air humidity, solar radiation, precipitation, and are used to simulate land surface states and fluxes. While LSMs can provide a continuous monitoring of land surface conditions, they still show discrepancies due to forcing and parameter errors, missing processes and inadequate model physics for particular areas or seasons. It is also possible to observe the land surface conditions from space. The modelling of land surface variables can be improved through the dynamical integration of these observations into LSMs. Remote sensing observations are particularly useful in this context because they are able to address global and continental scales. Low frequency microwave remote sensing has advantages because it can provide regular observations in all-weather conditions and at either daytime or night-time. A number of satellite-derived products relevant to the hydrological and vegetation cycles are already available from C-band radars such as the Advanced Scatterometer (ASCAT) or Sentinel-1. Assimilating these data into LSMs permits their integration in the process representation in a consistent way. The results obtained from assimilating satellites products provide land surface variables estimates that are generally superior to the model estimates or satellite observations alone. The main objective of this thesis was to improve the representation of land surface variables linked to the terrestrial water and carbon cycles in the ISBA LSM through the assimilation of ASCAT backscatter (sigma°) observations. An observation operator capable of representing the ASCAT sigma° from the ISBA simulated variables was developed. A version of the water cloud model (WCM) was successfully implemented over the Euro-Mediterranean area. The simulated values were compared with those observed from space. A more detailed quantification of the influence of various factors on the signal was made over southwestern France. Focusing on the Klaus storm event in the Landes forest, it was shown that the WCM was able to represent abrupt changes in vegetation biomass. It was also found that the WCM had shortcomings over karstic areas and over wheat croplands. It was shown that the latter was related to a discrepancy between the seasonal cycle of microwave vegetation optical depth (VOD) and leaf area index (LAI). Finally, the direct assimilation of ASCAT sigma° observations was assessed over southwestern France

What is cost-efficient phenotyping? Optimizing costs for different scenarios

Progress in remote sensing and robotic technologies decreases the hardware costs of phenotyping. Here, we first review cost-effective imaging devices and environmental sensors, and present a trade-off between investment and manpower costs. We then discuss the structure of costs in various real-world scenarios. Hand-held low-cost sensors are suitable for quick and infrequent plant diagnostic measurements. In experiments for genetic or agronomic analyses, (i) major costs arise from plant handling and manpower; (ii) the total costs per plant/microplot are similar in robotized platform or field experiments with drones, hand-held or robotized ground vehicles; (iii) the cost of vehicles carrying sensors represents only 5–26% of the total costs. These conclusions depend on the context, in particular for labor cost, the quantitative demand of phenotyping and the number of days available for phenotypic measurements due to climatic constraints. Data analysis represents 10–20% of total cost if pipelines have already been developed. A trade-off exists between the initial high cost of pipeline development and labor cost of manual operations. Overall, depending on the context and objsectives, “cost-effective” phenotyping may involve either low investment (“affordable phenotyping”), or initial high investments in sensors, vehicles and pipelines that result in higher quality and lower operational costs