Estimating chlorophyll with thermal and broadband multispectral high resolution imagery from an unmanned aerial system using relevance vector machines for precision agriculture

Precision agriculture requires high-resolution information to enable greater precision in the management of inputs to production. Actionable information about crop and field status must be acquired at high spatial resolution and at a temporal frequency appropriate for timely responses. In this study, high spatial resolution imagery was obtained through the use of a small, unmanned aerial system called AggieAir^TM. Simultaneously with the AggieAir flights, intensive ground sampling for plant chlorophyll was conducted at precisely determined locations. This study reports the application of a relevance vector machine coupled with cross validation and backward elimination to a dataset composed of reflectance from high-resolution multi-spectral imagery (VIS–NIR), thermal infrared imagery, and vegetative indices, in conjunction with in situ SPAD measurements from which chlorophyll concentrations were derived, to estimate chlorophyll concentration from remotely sensed data at 15-cm resolution. The results indicate that a relevance vector machine with a thin plate spline kernel type and kernel width of 5.4, having LAI, NDVI, thermal and red bands as the selected set of inputs, can be used to spatially estimate chlorophyll concentration with a root-mean-squared-error of 5.31 μg cm⁻², efficiency of 0.76, and 9 relevance vectors.     

An easy and low-cost method for preprocessing and matching small-scale amateur aerial photography for assessing agricultural land use in Burkina Faso

In recent decades, the Kou watershed in south-western Burkina Faso has suffered from poor water management. Despite the abundance of water, most water users regularly face water shortages because of the increase in the amount of land under irrigation. To help them achieve a more equitable allocation of irrigated land, local stakeholders need an easily managed low-cost tool for monitoring and mapping these irrigated zones. The aim of this study was to develop a fast and low-cost procedure for mosaicing and geo referencing amateur small-scale aerial photographs for land-use surveys. Sets of tens (2009) and hundreds (2007) of low-altitude aerial photographs, with a resolution of 0.4 m and 0.8 m, respectively, were used to create a detailed land-cover map of typical African small-scale irrigated agriculture. A commercially available stitching tool and GIS allowed geo referenced 'mono-images’ to be constructed; both mosaics were warped on a high-resolution SPOT image with a horizontal root mean square error (RMSE) of about 11 m. The RMSE between the two image datasets was 2 m. This approach is less sensitive to atmospheric conditions that are non-predictable in programming satellite imagery.     

Detecting shadows in multi-temporal aerial imagery to support near-real-time change detection

Multi-temporal aerial imagery captured via an approach called repeat station imaging (RSI) facilitates post-hazard assessment of damage to infrastructure. Spectral-radiometric (SR) variations caused by differences in shadowing may inhibit successful change detection based on image differencing. This study evaluates a novel approach to shadow classification based on bi-temporal imagery, which exploits SR change signatures associated with transient shadows. Changes in intensity (brightness from red–green–blue images) and intensity-normalized blue waveband values provide a basis for classifying transient shadows across a range of material types with unique reflectance properties, using thresholds that proved versatile for very different scenes. We derive classification thresholds for persistent shadows based on hue to intensity ratio (H/I) images, by exploiting statistics obtained from transient shadow areas. We assess shadow classification accuracy based on this procedure, and compare it to the more conventional approach of thresholding individual H/I images based on frequency distributions. Our efficient and semi-automated shadow classification procedure shows improved mean accuracy (93.3%) and versatility with different image sets over the conventional approach (84.7%). For proof-of-concept, we demonstrate that overlaying bi-temporal imagery also facilitates normalization of intensity values in transient shadow areas, as part of an integrated procedure to support near-real-time change detection.     

Downscaling Landsat-8 land surface temperature maps in diverse urban landscapes using multivariate adaptive regression splines and very high resolution auxiliary data

ABSTRACT

We propose a method for spatial downscaling of Landsat 8-derived LST maps from 100(30?m) resolution down to 2–4?m with the use of the Multiple Adaptive Regression Splines (MARS) models coupled with very high resolution auxiliary data derived from hyperspectral aerial imagery and large-scale topographic maps. We applied the method to four Landsat 8 scenes, two collected in summer and two in winter, for three British towns collectively representing a variety of urban form. We used several spectral indices as well as fractional coverage of water and paved surfaces as LST predictors, and applied a novel method for the correction of temporal mismatch between spectral indices derived from aerial and satellite imagery captured at different dates, allowing for the application of the downscaling method for multiple dates without the need for repeating the aerial survey. Our results suggest that the method performed well for the summer dates, achieving RMSE of 1.40–1.83?K prior to and 0.76–1.21?K after correction for residuals. We conclude that the MARS models, by addressing the non-linear relationship of LST at coarse and fine spatial resolutions, can be successfully applied to produce high resolution LST maps suitable for studies of urban thermal environment at local scales.     

Visualizing the ill-posedness of the inversion of a canopy radiative transfer model: A case study for Sentinel-2

Monitoring biophysical and biochemical vegetation variables in space and time is key to understand the earth system. Operational approaches using remote sensing imagery rely on the inversion of radiative transfer models, which describe the interactions between light and vegetation canopies. The inversion required to estimate vegetation variables is, however, an ill-posed problem because of variable compensation effects that can cause different combinations of soil and canopy variables to yield extremely similar spectral responses. In this contribution, we present a novel approach to visualise the ill-posed problem using self-organizing maps (SOM), which are a type of unsupervised neural network. The approach is demonstrated with simulations for Sentinel-2 data (13 bands) made with the Soil-Leaf-Canopy (SLC) radiative transfer model. A look-up table of 100,000 entries was built by randomly sampling 14 SLC model input variables between their minimum and maximum allowed values while using both a dark and a bright soil. The Sentinel-2 spectral simulations were used to train a SOM of 200 × 125 neurons. The training projected similar spectral signatures onto either the same, or contiguous, neuron(s). Tracing back the inputs that generated each spectral signature, we created a 200 × 125 map for each of the SLC variables. The lack of spatial patterns and the variability in these maps indicate ill-posed situations, where similar spectral signatures correspond to different canopy variables. For Sentinel-2, our results showed that leaf area index, crown cover and leaf chlorophyll, water and brown pigment content are less confused in the inversion than variables with noisier maps like fraction of brown canopy area, leaf dry matter content and the PROSPECT mesophyll parameter. This study supports both educational and on-going research activities on inversion algorithms and might be useful to evaluate the uncertainties of retrieved canopy biophysical and biochemical state variables.     

Combining UAV-based plant height from crop surface models,visible, and near infrared vegetation indices for biomass monitoring in barley

In this study we combined selected vegetation indices (VIs) and plant height information to estimate biomass in a summer barley experiment. The VIs were calculated from ground-based hyperspectral data and unmanned aerial vehicle (UAV)-based red green blue (RGB) imaging. In addition, the plant height information was obtained from UAV-based multi-temporal crop surface models (CSMs). The test site is a summer barley experiment comprising 18 cultivars and two nitrogen treatments located in Western Germany. We calculated five VIs from hyperspectral data. The normalised ratio index (NRI)-based index GnyLi (Gnyp et al., 2014) showed the highest correlation (R² = 0.83) with dry biomass. In addition, we calculated three visible band VIs: the green red vegetation index (GRVI), the modified GRVI (MGRVI) and the red green blue VI (RGBVI), where the MGRVI and the RGBVI are newly developed VI. We found that the visible band VIs have potential for biomass prediction prior to heading stage. A robust estimate for biomass was obtained from the plant height models (R² = 0.80–0.82). In a cross validation test, we compared plant height, selected VIs and their combination with plant height information. Combining VIs and plant height information by using multiple linear regression or multiple non-linear regression models performed better than the VIs alone. The visible band GRVI and the newly developed RGBVI are promising but need further investigation. However, the relationship between plant height and biomass produced the most robust results. In summary, the results indicate that plant height is competitive with VIs for biomass estimation in summer barley. Moreover, visible band VIs might be a useful addition to biomass estimation. The main limitation is that the visible band VIs work for early growing stages only.     

Optimising three-band spectral indices to assess aerial N concentration,N uptake and aboveground biomass of winter wheat remotely in China and Germany

Remotely and accurately quantifying the canopy nitrogen status in crops is essential for regional studies of N budgets and N balances. In this study, we optimised three-band spectral algorithms to estimate the N status of winter wheat. This study extends previous work to optimise the band combinations further and identifies the optimised central bands and suitable bandwidths of the three-band nitrogen planar domain index (NPDI) for estimating the aerial N uptake, N concentration and aboveground biomass. Analysis of the influence of bandwidth change on the accuracy of estimating the canopy N status and aboveground biomass indicated that the suitable bandwidths for optimised central bands were 37 nm at 846 nm, 13 nm at 738 nm and 57 nm at 560 nm for assessing the aerial N uptake and were 37 nm at 958 nm, 21 nm at 696 nm and 73 nm at 578 nm for the assessment of the aerial N concentration and were 49 nm at 806 nm, 17 nm at 738 nm and 57 nm at 560 nm for the estimation of aboveground biomass. The optimised three-band NPDI could consistently and stably estimate the aerial N uptake and aboveground biomass of winter wheat in the vegetative stage and the aerial N concentration in the reproductive stage compared to the fixed band combinations. With suitable bandwidths, the broadband NPDI demonstrated excellent performance in estimating the aerial N concentration, N uptake and biomass. We conclude that the band-optimised algorithm represents a promising tool to measure the improved performance of the NPDI in estimating the aerial N uptake and biomass in the vegetative stage and the aerial N concentration in the reproductive stage, which will be useful for designing improved nitrogen diagnosis systems and for enhancing the applications of ground- and satellite-based sensors.     

CONSIDERATIONS OF SPECIFICATIONS FOR SMALL FORMAT AERIAL PHOTOGRAPHY

Following the success of the symposium held by the Photogrammetric Society on 31st March, 1987, a second one-day symposium was convened to give further support to the rapidly growing interest in the development of small format photography for aerial survey. This was held at Burlington House, London on 1st November, 1988. Over 70 people attended. Six of the papers are reproduced here, dealing with a range of topics including specifications for and navigation of small format aerial photography, the creation of digital terrain models from 35 mm photography, wildlife surveys in the Arctic and survey photography from model aircraft and tethered balloons.     

Spectral Characteristics of Domestic and Wild Mammals

Few studies have recorded the spectral signatures of domesticated live animals and in particular few have examined wild species. Using in situ radiometry we acquired visual and near infrared spectral signatures of wild elk (Cervus elaphus) and domesticated cattle (Bos taurus) and horses (Equus caballus). Signatures were significantly different among species across all bands with the exception of cattle and horses in the red band. Further research is needed to determine if the shallower slopes in the red-shift region of the animal signatures would allow for distinction from vegetation using various remote sensors. Application of in situ spectral signatures to remotely sensed imagery could provide an efficient method for counting wildlife.     

Ocean Internal Wave Observations Using Space Shuttle and Satellite Imagery

Abstract

Space shuttle photographs and satellite radar (SAR) images provide an excellent view of high‐contrast ocean features such as internal waves, fronts, eddies, oil slicks, and cloud patterns which contain the signatures of atmospheric processes. Since ocean internal waves generate local currents which modulate surface wavelets and slicks, we have been able to detect packets of internal waves in space shuttle photographs and radar imagery of the Atlantic, Pacific, and Indian Oceans. A global database on internal waves has been developed at our center with support from ONR and NASA, and is accessible on the Internet. The database includes visible and radar imagery. To test the database, digitally orthorectified images were used for dynamic and statistical analysis of internal waves. In the deep ocean we found the wavelength distribution to be Gaussian while in the coastal ocean it is Rayleigh. Results have also been applied to non‐linear evolution studies of ocean internal waves, atmospheric solitary waves and to estimate ocean currents.     

Estimation of green grass/herb biomass from airborne hyperspectral imagery using spectral indices and partial least squares regression

The main objective was to determine whether partial least squares (PLS) regression improves grass/herb biomass estimation when compared with hyperspectral indices, that is normalised difference vegetation index (NDVI) and red-edge position (REP). To achieve this objective, fresh green grass/herb biomass and airborne images (HyMap) were collected in the Majella National Park, Italy in the summer of 2005. The predictive performances of hyperspectral indices and PLS regression models were then determined and compared using calibration (n = 30) and test (n = 12) data sets. The regression model derived from NDVI computed from bands at 740 and 771 nm produced a lower standard error of prediction (SEP = 264 g m⁻²) on the test data compared with the standard NDVI involving bands at 665 and 801 nm (SEP = 331 g m⁻²), but comparable results with REPs determined by various methods (SEP = 261 to 295 g m⁻²). PLS regression models based on original, derivative and continuum-removed spectra produced lower prediction errors (SEP = 149 to 256 g m⁻²) compared with NDVI and REP models. The lowest prediction error (SEP = 149 g m⁻², 19% of mean) was obtained with PLS regression involving continuum-removed bands. In conclusion, PLS regression based on airborne hyperspectral imagery provides a better alternative to univariate regression involving hyperspectral indices for grass/herb biomass estimation in the Majella National Park.     

Use of Kendall's coefficient of concordance to assess agreement among observers of very high resolution imagery

Ground-based vegetation monitoring methods are expensive, time-consuming and limited in sample size. Aerial imagery is appealing to managers because of the reduced time and expense and the increase in sample size. One challenge of aerial imagery is detecting differences among observers of the same imagery. Six observers analysed a set of 1-mm ground sample distance aerial imagery for graminoid species composition and important ground-cover characteristics. Kendall's coefficient of concordance (W) was used to measure agreement among observers. The group of six observers was concordant when assessed as a group. When each of the observers was assessed independently against the other five, lack of agreement was found for those graminoid species that were difficult to identify in the aerial images.     

Assessing the transferability of airborne laser scanning and digital aerial photogrammetry derived growing stock volume models

Three-dimensional (3D) data from airborne laser scanning (ALS) and, more recently, digital aerial photogrammetry (DAP) have been successfully used to model forest attributes. While multi-temporal, wall-to-wall ALS data is not usually available, aerial imagery is regularly acquired in many regions. Thus, the combination of ALS and DAP data provide a sufficient temporal resolution to properly monitor forests. However, field data is needed to fit new forest attribute models for each 3D data acquisition, which is not always affordable. In this study, we examined whether transferability of growing stock volume (GSV) models may provide an improvement in the efficiency of forest inventories updating. We used two available ALS datasets acquired with different characteristics in 2009 and 2010, respectively, generated two DAP point clouds from imagery collected in 2010 and 2017, and utilized field data from two ground surveys conducted in 2009 and 2016-2017. We first analyzed the stability of point cloud derived metrics. Then, Support Vector Regression models based on the most stable metrics were fitted to assess model transferability by applying them to other datasets in four different cases: (1) ALS-ALS, (2) DAP-DAP temporal, (3) ALS-DAP and (4) ALS-DAP temporal. Some metrics were found to be enough stable in each case, so they could be used interchangeably between datasets. The application of models to other datasets resulted in unbiased predictions with relative root mean square error differences ranging from -8.27% to 14.59%. Results demonstrated that 3D-based GSV models may be transferable between point clouds of the same type as well as point clouds acquired using different technologies such as ALS and DAP, suggesting that DAP data may be used as a cost-efficient source of information for updating ALS-assisted forest inventories.     

Role of aerial photographs and landsat imagery in forest resources surveys

Many of the data needs for efficient management of forest resources can be met by aerial photographs. Commercially important tree species can be distinguished from other less important miscellaneous species with the help of aerial photographs. Forests can be classified according to their height and density classes. Aerial photographs have become indispensable for mapping of forests and preparation of forest inventories. A comparison of interpretation results obtained from landsat imagery and aerial photographs (1 ∶ 10,000 Black and White panchromatic photography) with respect to forestry interpretation is given. It is pointed out that the imagery obtained from satellities can be used for reconnaissance of a region and for deciding the priorities for carrying out more detailed surveys of forest resources with the help of air photointerpretation techniques     

VEGETATION AS A CRITERION FOR INTERPRETING THE MOISTURE REGIME OF A LANDSCAPE

Vegetation signatures on aerial photographs and space imagery are used as indicators of soil moisture differences in a Siberian taiga landscape. The authors use remote sensing products to identify all major habitat types resulting from variable soil moisture regimes that were observed on the ground. These types are described, as are their interpretation keys and the effects of scale change on habitat discrimination. Translated from: Kosmicheskiye metody izucheniya prirodnoy sredy Sibiri i Dal'nego Vostoka, Novosibirsk, 1983, pp. 63-74.     

热带森林植被的动态变化遥感监测

本文讨论了以热带森林植被为主体的再生资源的面积动态变化监测。研究中包括两个部分。首先,我们利用多时相遥感图像对大面积的西双版纳州进行地类判读,系统地分析了森林植被的动态变化。其次,利用Landsat MSS和TM数据对自然保护区的动态变化进行了包含无监督分类和归一化差值植被指数分析的数字图像处理,变化分类也相当符合实际。总的实验结果表明,这种监测方法是很有效的,可在再生资源监测中特别是在森林植被监测中加以推广应用。     

Resolution vs. image quality in pre-tsunami imagery used for tsunami impact models in Aceh,Indonesia

Land cover roughness coefficients (LCRs) have been used in multivariate spatial models to test the mitigation potential of coastal vegetation to reduce impacts of the 2004 tsunami in Aceh, Indonesia. Previously, a Landsat 2002 satellite imagery was employed to derive land cover maps, which were then combined with vegetation characteristics, i.e., stand height, stem diameter and planting density to obtain LCRs. The present study tested LCRs extracted from 2003 and 2004 Landsat (30 m) images as well as a combination of 2003 and 2004 higher spatial resolution SPOT (10 m) imagery, while keeping the previous vegetation characteristics. Transects along the coast were used to extract land cover, whenever availability and visibility allowed. These new LCRs applied in previously developed tsunami impact models on wave outreach, casualties and damages confirmed previous findings regarding distance to the shoreline as a main factor reducing tsunami impacts. Nevertheless, the models using the new LCRs did not perform better than the original one. Particularly casualties models using 2002 LCRs performed better (δAIC > 2) than the more recent Landsat and SPOT counterparts. Cloud cover at image acquisition for Landsat and low area coverage for SPOT images decreased statistical predictive power (fewer observations). Due to the large spatial heterogeneity of tsunami characteristics as well as topographic and land-use features, it was more important to cover a larger area. Nevertheless, if more land cover classes would be referenced and high resolution imagery with low cloud cover would be available, the full benefits of higher spatial resolution imagery used to extract more precise land use roughness coefficients could be exploited.     

Interpretation of Synthetic Aperture Radar (SAR) imagery for geological appraisal: A comparative study in Anantapur district of Andhra Pradesh

The paper describes the details of a comparative study of geological interpretations carried out from Synthetic Aperture Radar (SAR) imagery, Landsat MSS (B & W) imagery and Aerial Photographs, covering 2100 sq km of area in Anantapur district of Andhra Pradesh. The area comprises Peninsular—Gneissic Complex and rocks of Dharwar and Cuddapah Super Groups beside the Quaternary alluvial deposits along the Penneru river and its tributaries. Geomorphologically the areas is represented by denudational, fluvial and structural landforms. The study indicates that the details of the geological and geomorphological maps prepared from SAR imagery and aerial photographs are comparable despite the smaller scale of SAR imagery while the same are not exhibited in Landsat imagery mainly due to its low resolution. Although broad lithological units are possible to be discriminated on SAR as well as aerial photographs, some of the finer rock types viz. gabbroic dykes could be discriminated from the delerite dykes in the SAR imagery due to their different surface roughness. Stereoscopic coverage and enhanced micro-relief of SAR imagery gives better geomorphological details in comparison to aerial photographs. A detailed study of lineaments has also been carried out which shows that in SAR imagery there is over-representation of short lineaments due to enhanced micro-relief and radarshadow effects across the look direction and under-representation of lineaments along the look direction. Landsat imagery is perhaps the best for demarcating lineaments of regional magnitude while aerial photographs are good for depicting shorter lineaments. However, certain lineaments seen in SAR imagery are often not continuously seen on aerial photographs.     

The brightness temperature adjusted dust index: An improved approach to detect dust storms using MODIS imagery

Moderate Resolution Imaging Spectroradiometer (MODIS) imagery provides a good data source for timely and accurate monitoring of dust storms. However, effective MODIS-based dust indices are inadequate. In this study, we proposed an improved brightness temperature adjusted dust index (BADI) by integrating the brightness temperatures of three thermal infrared MODIS bands: band₂₀ (3.66–3.84 μm), band₃₁ (10.78–11.28 μm) and band₃₂ (11.77–12.27 μm). We used the BADI to monitor several representative dust storms over the Northeast Asia between 2000 and 2011. When compared to commonly used MODIS-based dust indices, such as the brightness temperature difference index in band₃₂ and band₃₁ (BTD_32-31) and the normalized difference dust index (NDDI), the BADI captured the spatial extent and density of dust storms more accurately. The BADI detected dust storm extent with an overall accuracy >90%, which was 7% and 29% higher than the results derived from BTD_32-31 and NDDI, respectively. The BADI also demonstrated good agreement with the density indicator of MODIS Deep Blue Aerosol Optical Depth (R² = 0.59, P < 0.01). We suggest that the BADI is an effective tool to monitor large-scale dust storms.