Next Article in Journal
Commercial Cultivation of Australian Wild Oryza spp.: A Review and Conceptual Framework for Future Research Needs
Next Article in Special Issue
Estimation of Winter Wheat Residue Coverage Based on GF-1 Imagery and Machine Learning Algorithm
Previous Article in Journal
Utilization of By-Product of Groundnut Oil Processing for Production of Prodigiosin by Microbial Fermentation and Its Novel Potent Anti-Nematodes Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Methodology for the Automated Delineation of Crop Tree Crowns from UAV-Based Aerial Imagery by Means of Morphological Image Analysis

by
Juan Manuel Ponce
1,*,
Arturo Aquino
1,
Diego Tejada
1,
Basil Mohammed Al-Hadithi
2,3 and
José Manuel Andújar
1
1
Research Center in Technology, Energy and Sustainability (CITES), University of Huelva, Huelva-Palos de la Frontera Rd., Palos de la Frontera, 21819 Huelva, Spain
2
Centre for Automation and Robotics, UPM-CSIC. J., Gutierrez Abascal, 2., 28006 Madrid, Spain
3
Department of Electrical, Electronics, Control Engineering and Applied Physics, School of Industrial Design and Engineering, Technical University of Madrid, Ronda de Valencia, 3., 28012 Madrid, Spain
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(1), 43; https://doi.org/10.3390/agronomy12010043
Submission received: 22 November 2021 / Revised: 20 December 2021 / Accepted: 22 December 2021 / Published: 25 December 2021
(This article belongs to the Special Issue Remote Sensing in Smart Agriculture)

Abstract

:
The popularisation of aerial remote sensing using unmanned aerial vehicles (UAV), has boosted the capacities of agronomists and researchers to offer farmers valuable data regarding the status of their crops. This paper describes a methodology for the automated detection and individual delineation of tree crowns in aerial representations of crop fields by means of image processing and analysis techniques, providing accurate information about plant population and canopy coverage in intensive-farming orchards with a row-based plant arrangement. To that end, after pre-processing initial aerial captures by means of photogrammetry and morphological image analysis, a resulting binary representation of the land plot surveyed is treated at connected component-level in order to separate overlapping tree crown projections. Then, those components are morphologically transformed into a set of seeds with which tree crowns are finally delineated, establishing the boundaries between them when they appear overlapped. This solution was tested on images from three different orchards, achieving semantic segmentations in which more than 94% of tree canopy-belonging pixels were correctly classified, and more than 98% of trees were successfully detected when assessing the methodology capacities for estimating the overall plant population. According to these results, the methodology represents a promising tool for automating the inventorying of plants and estimating individual tree-canopy coverage in intensive tree-based orchards.

1. Introduction

Modern agricultural practices developed around the precision agriculture (PA) paradigm, demand data-collecting systems for assembling information regarding the spatial and temporal variability of those factors of influence in agricultural production [1]. This data-driven approach is aimed at developing decision-making frameworks to help farmers with their daily tasks [2], thus supporting the eventual optimisation of the farming-inputs management and favouring the improvement of the overall agricultural activity in terms of crop productivity, sustainability, and profitability [3].
Hence, non-invasive data acquisition regarding crops condition has become of great interest to the agricultural sector [4,5]. Thus, extensive related research has been carried out on the development of solutions aimed at improving farming processes, by developing decision support models by means of remotely sensed data [6,7,8,9,10,11]. In this sense, aerial remote sensing has had a pivotal role in agriculture over time. Indeed, its capacities for systematically assessing large tracts of land soon caught the attention of researchers, so the first studies on the use of aerial images in agriculture date back to the 1950s [12]. Notwithstanding, it is in the 1970s when their use began to stand out when satellite images started to become popular as data sources for agronomic researchers [1]. Initial applications, which focused mainly on the land cover classification [13,14] distinguishing between types of crops and vegetation present in the surveyed fields, soon became more complex as technology advances allowed higher spatial resolution and increased computational power, necessary for processing and analysing the collected data.
However, in recent years, there has been an increasing interest in aerial remote sensing technologies within modern agriculture, mostly due to the popularisation of unmanned aerial vehicles (UAV). Its advantages over other forms of aerial imaging, although debated, have caught the attention of farmers and agronomists, and UAV-based remote sensing has strengthened as a data-acquisition technology for PA [15,16,17,18], with a vast amount of related research published, that has consistently underlined its usefulness in agronomy, and its applicability in diverse related domains. Among them, UAV-based remotely sensed imagery has shown huge potential and applicability in crop field mapping and characterisation [19]. Low-altitude aerial imagery can offer some advantages when compared to satellite-based imaging, in terms of both spatial and temporal resolution, or when performing image acquisition under cloudy conditions [20,21]. By conveniently processing these high-spatial resolution data by means of diverse image analysis techniques, it can be yielded accurate representations of the fields surveyed [22], thus enabling the possibility to discriminate relevant information regarding the orography of the cultivation area, biological characteristics of the vegetation, dendrometric features of the canopies, or phenological plant status, among others [23,24,25,26].
Amongst the different issues that can be approached using UAV-based imagery, individual tree detection and characterization are of special interest. Indeed, accurate knowledge about plant population, and the geometric characteristics of the cultivation areas, regarding the spatial relationship between plants and their individual canopy-related features, can offer useful information for farmers. On the one hand, accurate data about tree populations can lead to meaningful statistics, such as crop density, directly involved in the development of early yield estimation models, or in the enhancement of watering/fertilization strategies [27,28]. Equally, geometric features of plant canopies can offer relevant indicators, also related to the yield capacities of the crop, or the water and nutritional status of the plants [29]. Due to their high-throughput within PA context, and given the difficulties in collecting these data in traditional manners, all mostly based on in-field visual inspections and costly labour-intensive manual measurements [30], the use of UAV-based imagery in conjunction with computer vision methodologies has gained relevance within the research community, seeking solutions for obtaining this type of information in the most autonomous and accurate way possible. Thus, Zortea et al. [31] proposed a method for detecting citrus-trees from aerial images collected with a quadcopter. The authors undertook the detection of problem processing those images by means of deep learning, with a double convolutional neural network (CNN) configuration, first identifying the tree rows, and then classifying the set of candidates of tree-occurrences computed. A similar approach was followed by Csillik et al. [32], using CNNs and simple linear iterative clustering (SLIC) for the purpose of detecting citrus-trees. In the same vein, Ampatzidis and Partel [33] also drew on deep learning techniques for automating citrus-tree detection from aerial imagery and additionally estimating different features of the plants identified, as the area of their canopy projections. Recio et al. [34] focused on citrus crops too, proposing an image analysis procedure which, from an initial k-means clustering-based pixel classification of the aerial captures, allowed individual plant identification and tree crown-canopy coverage estimations. Ok and Ozdarici-Ok [35] presented a method for identifying individual tree canopies, applicable in citrus tree-based crops. They proposed an ad-hoc orientation-based radial symmetry transform for detecting tree candidates from digital surface models (DSM) of crop plots, and then individually delineating them by means of active contour models. Modica et al. [36] proposed a geographic object-based image analysis (GEOBIA) approach for segmenting crop vegetation from UAV-based multispectral images. To that end, they developed and compared different solutions based on four different machine-learning object-classification algorithms: K-nearest neighbour, support vector machines, random forests, and normal bayes, all tested in a citrus-tree orchard and in an onion field. Dong et al. [37] combined digital height models (DHM) with vegetation indices for individually segmenting apple and pear-tree crowns from aerial representations of cultivation areas. Likewise, Marques et al. [38] proposed the generation of a model based on elevation information, obtained from initial aerial captures, and the use of vegetation indices for detecting chestnut trees, and subsequently extracting different morphological characteristics from each of them. Salamí et al. [39] developed a framework for detecting and counting olive trees by processing aerial captures on-the-fly, using and comparing two different segmentation approaches: colour-based and stereo vision-based segmentations. Jiang et al. [40] used UAV imagery with the aim of automating the identification of papaya and lemon trees. To that end, they proposed a detection algorithm that processes images, previously transformed into the Lab colour space, on the basis of scale-space filtering (SSF). Johansen et al. [41] explored the impact of the height at which aerial imagery acquisition is performed when extracting geometrical characteristics of lychee-trees, for assessing the effects of pruning in the corresponding orchard.
Although extensive research has been carried out on this subject, there is still room for improvement, and some critical related issues have only been partially resolved. On one hand, most studies to date tend to focus on specific crops, and the generalisation of the proposed solutions and their usefulness with different types of crop trees remains uncertain. On the other hand, it can be appreciated a lack of research conducted in high density orchards, where plant arrangement is prone to the occurrences of overlapping canopies between adjacent trees, thus hindering the task of automatically detecting the plants at an individual level by means of computer vision techniques. Indeed, this issue could be of great interest given the trend towards intensification of agriculture in recent years, and the increasingly common adoption of super-intensive planting schemes [42,43,44].
Within this context, this paper presents a novel methodology for identifying crop trees in low-altitude UAV-based aerial images by means of image analysis techniques, with the ability to automatically delineate those tree crowns detected, thus enabling the estimation of features related to the size and morphology of each individual canopy. This solution is aimed at being applicable on high-density orchards with a regular row-based planting arrangement. Moreover, it was designed to be applicable to scenarios where plant arrangement is conducive to the occurrence of overlapping between the canopies of adjacent trees. Even in those situations in which the degree of overlapping was so extreme that lead to the loss of the shape of the tree crowns, in their corresponding pixel representations, the methodology showed outstanding performance. In addition, whereas most of the related research to date has tended to focus on specific crops, the developed methodology was designed to work independently of the type of trees under study. To this end, it was conceived dispensing with the use of vegetation indices, or any kind of machine learning-based approach, which would have to be modelled from any fixed training configuration, obtained from each specific crop. Therefore, a comprehensive framework based on morphological image analysis is proposed, which after transforming multispectral aerial images of an orchard according to the procedure proposed by Sarabia et al. [45], is capable of individually segmenting the pixel area of each crown projection, despite the presence of intra-row overlapping tree crowns, yielding estimations of the overall plant population as well as of individual tree canopy coverage.

2. Materials and Methods

2.1. Case Study Site

The developed methodology has been tested with three different types of tree-based crops: olive (Olea europaea L.), lemon (Citrus limon L.), and orange (Citrus sinesis L.) trees. To that end, three different orchards, one per crop, were surveyed.
In the case of the olive trees, the grove under study is in Gibraleón, province of Huelva (Andalusia, Southwest Spain). This land plot, centred in DMS coordinates 7°02′48.44″ W and 37°20′39.80″ N, spreads over an area of approximately 17.5 ha, with a total of 3916 olive trees planted. On the other hand, the two citrus-tree orchards can also be found in this province, in the surroundings of the village of La Redondela, Isla Cristina. The lemon-tree orchard is located centred in the coordinates 7°17′06.93″ W and 37°14′07.24″ N, and it has an approximate extent of 1.37 ha and 552 trees planted. The orange-tree orchard can be found in the coordinates 7°18′06.75″ W and 37°13′49.55″ N, with an area of 1.94 ha and 781 trees. Figure 1 shows aerial images of these three orchards. It should be remarked that none of these aerial images properly belong to the present study, as they were taken from the Spanish National Geographic Institute (“Instituto Geográfico Nacional”) [46]. They were not used at any moment throughout the development, nor the validation of the methodology described in this paper, so they are proposed just for illustrative purposes.
It should be underscored that they all show a regular row-based planting pattern. In this respect, it is worth noting that in the case of the orange grove, there is a very high degree of overlap between intra-row adjacent trees. This is the reason it was selected for the study, with the aim of carrying out a sort of stress test with which to evaluate the developed methodology.

2.2. Image Acquisition Equipment

Aerial images of the land plots under study were collected with a multispectral sensor MicaSense RedEdge-MTM (MicaSense, Inc., Seattle, WA, USA). This device is able to simultaneously capture in five different discrete spectral bands, and its main features are summarised in Table 1.
The sensing equipment was completed with a GPS module, directly connected to the camera for georeferencing the images as they were captured, and a 5-band downwelling light sensor (DLS). From the data collected by this latter sensor, it is possible to make corrections in flight regarding global lighting changes. In addition, a MicaSense Reflectance Panel was included in the setup. Captures of this panel provide accurate information regarding the amount of light reaching the ground when performing image acquisition, so these data are used for transforming raw pixels values to absolute reflectance when later processing the collected images.
As UAV, a DJITM Matrice 100 (SZ DJI Technology Co., Ltd., Shenzhen, Guangdong, China) quadcopter was used for the olive case study. On the other hand, for the two remaining cases corresponding to the citrus-tree orchards, image acquisition was performed with a DJITM Phantom 3 (SZ DJI Technology Co., Ltd., Shenzhen, Guangdong, China). Because of the dimensions of both orchards, this task could be carried out with a smaller UAV with lesser capabilities. The main characteristics of both platforms are listed in Table 2. The UAVs, along with the sensing equipment specified in Table 1, are also shown in Figure 2.

2.3. Flights Planning and Execution

All those flights carried out for image acquisition were planned with the DJIFlightPlannerTM (AeroScientific-Spatial Scientific Pty. Ltd., Adelaide, Australia; Photometrix, Kew, Australia; Dechant Consulting Services Inc., Beaux Arts Village, WA, USA) software. They all shared the same configuration, aiming at capturing aerial images with a forward overlap of 85%, a lateral overlap of 65%, and a ground sample distance (GSD) of 4.8 cm. Flight missions were performed with a height of 70 m with regard to the take-off point, a cruising speed of 4.2 m/s, and time between captures of 1.5 s. They were executed autonomously in GPS mode and monitored with the ground control software Litchi (VC Technology, Ltd., London, UK).
The olive grove was flown on 13 June 2019, approximately between 11 a.m. and 1 p.m. 8865 images were taken per spectral band, yielding a total of 44,325 captures. In the case of citrus trees, two different flights were made, one for each type of crop. Both flights took place on 1 March 2020, in the same time slot as in the olive grove case.
On the other hand, a total of 1520 images, 304 per spectral band, were taken in the orange grove. 1580 (316 per spectral band) images were acquired in the case of the lemon-tree orchard. In Figure 3, some examples of images captured during those flights are shown. Specifically, they correspond to the orange-tree orchard.

2.4. Developed Methodology

2.4.1. Image Preprocessing

In order to accomplish the main goal of the developed methodology, which is the individual delineation of the tree crowns boundaries in aerial images, it is necessary to transform initial captures into a data structure suitable for that purpose. The first pre-processing stage was carried out by the application of the procedure proposed by Sarabia et al. [45]. It puts forward a comprehensive method to automatically process multispectral aerial images of a land plot, with the objective of obtaining a unique binary image representative of it, where crop tree-belonging pixel regions are segmented from the ground, and location points of each individual plant are estimated. To that end, drawing from the set of aerial captures, an initial photogrammetric imaging procedure generates a 3D-point representation of the plot to survey. In this case, PIX4DmapperTM (Pix4D, Prilly, Switzerland) was used as photogrammetry software for the purpose of yielding a high-density point cloud for each case study addressed. It should be noted that this software was configured such that each resulting 3D point could be correctly re-projected in a minimum of 3 different images, and no image scaling was applied. Once the 3D-point representations were obtained, corresponding digital height models (DHM) were computed by applying an inverse distance weighing (IDW) interpolation [47]. This task was conducted using ArcGisTM 10.3 (Esri, Inc., Redlands, CA, USA) along with its Geostatistical Analyst Tools extension. For this procedure, the software was configured for a fixed neighbourhood search, with a number of input neighbours up to 4, a search radius of 10 m, and a weighting exponent of 2. Then, each DHM was approached as a greyscale image. Height information at each point was treated as a pixel-intensity value, thus the models were processed by means of morphological analysis and statistical thresholding with the aim of segmenting those pixel-areas belonging to crop-tree canopies from the rest of that image. To that end, a background estimation was computed for each DHM-image, aiming at filtering each tree-crown projection by replacing the grey values of their corresponding pixels with the minimum intensity value found in the set of ground pixels present in their nearest surroundings. For that purpose, the process is carried out by iteratively applying a set of morphological opening operations [48], using a disk-shaped structuring element whose radius, β, is increased at each step, in an attempt to offer the optimum filtering for each tree canopy. The stop condition occurs when this kernel is large enough to completely contain the biggest tree-crown projection in the image. Thus, the maximum value to be reached by the radius of the structuring element corresponds to one of the two parameters that need to be set in order to apply this methodology, so it should be configured to guarantee that the corresponding structuring element is able to cover the largest individual canopy. The second required parameter relates to the height below which regional maxima of the image are discarded, assuming that tree crowns are represented at higher elevations, thus favouring their segmentation. Hence, it should be configured in accordance with the height expected for the tree crops. Table 3 shows the values of both configuration parameters for each case study addressed.
After subtracting the background estimations from the corresponding grayscale representations of the DHMs, Otsu’s method [49] was applied for eventually binarising those images, thus yielding the eventual canopy mass-ground segmentation. Finally, the estimation of the tree location points was undertaken by analysing the size of each of the connected components in the resulting segmentations, with regard to the dominant tree crown size in the binary image. These connected components can be defined as those groupings of pixels with the same intensity value and a neighbourhood relationship (8-connectivity was considered in this study) [50], which, in this case, potentially represent the projections of tree canopies. In Figure 4, the workflow followed by this methodology is illustrated.

2.4.2. Image Processing

The starting point of the procedure here presented, which comprises the core of the methodology, is a binary image representative of the orchard understudy, and where tree canopies are segmented from the ground, along with a set of potential location points for the individual crop plants in that image. Figure 5 shows, as an example, the binary image resulting from pre-processing the aerial captures collected in one of the study sites considered in this work. Additionally, the estimated tree location points are represented on it.
Hence, given both pieces of information, the general idea behind the developed methodology goes through the individual processing of each of the connected components. It focuses on transforming those components into small disjoint groupings of pixels, somehow representative of each individual tree contained within them. These sorts of pixel-seeds are aimed at backing the eventual splitting of the connected component according to the plants contained. Hence, they are used to support the final individual tree coverage segmentation, by a Watershed transform-based operation. This process of transforming and segmenting the connected components is mainly based on the application of different morphological operators, which are formally defined in Appendix A, so the reader is encouraged to consult it for further information. The whole developed methodology that is presented in detail hereafter is briefly described in the flowchart shown in Figure 6.
Given the binary image IBIN yielded by the preprocessing, each of its connected components is individually processed, as mentioned above. First, it is checked the number of trees in the component, estimated during preprocessing. Given a connected component of IBIN, cci, if the estimated number of trees contained in it, nTi, is just one, it is not necessary to conduct any processing at all. When nTi is greater than one, its corresponding pixel region is initially isolated in an alternative intensity matrix, embedding it in the smallest all-black-pixel box that can contain it. For this auxiliary image, its Euclidean distance transform [51] is calculated. Thus, given this binary image ICCi, which encloses cci in its proper bounding box, its distance transform, IDTi, is calculated in accordance with the following definitions:
I D T i = D T ( I C C i ) ,
[ D T ( I C C i ) ] ( p ) = min { d ε ( p , q ) | f ( q ) = 0 } ,
d ε ( p , q ) = ( p x q x ) 2 + ( p y q y ) 2
where, given the binary image f and the pixel p, [DT(f)(p)] refers to the Euclidean distance from p to its closest background pixel. The resulting matrix, IDTi, can be seen and subsequently treated as a greyscale image, as illustrated in Figure 7.
Subsequently, the regional maxima of the resulting intensity matrix IDTi are explored as tree representative. The goal here is to isolate those more representative regional maxima, as they can be expected to take place where the centres of mass of the trees within the component are potentially located. This is carried out on the basis of the h-maxima transform, which allows us to remove from IDTi all maxima with depth lower or equal than scalar h. It can be computed as the morphological reconstruction by dilation ( R δ ) of IDTi from the marker image resulting from subtracting h from the intensity value of each pixel in IDTi. Thus, irrelevant maxima are firstly removed from IDTi by computing:
I D T i , h = H M A X h ( I D T i ) = R I D T i δ ( I D T i h )
The surviving regional maxima in IDTi,h are those considered representatives and, therefore, are extracted from the image. It is achieved by subtracting from IDTi,h the result of suppressing from it all the remaining regional maxima:
I R M A X i , h = I D T i , h R I D T i , h δ ( I D T i , h 1 )
Note that the rightmost term of the subtraction is the h-maxima transform setting the value of h to 1. This value enables the removal of all remaining regional maxima from IDTi,h, as all of them will have, at least, that height. Then, by thresholding IRMAXi,h, it can be achieved a binary mask containing those maxima representing the centre of mass of the trees:
[ I R M A X B I N i , h ] ( p ) = { 255   i f   I R M A X i , h ( p ) > 0 0               i n   a n y   o t h e r   c a s e
The result of this process depends on the definition of the h value in Equation (4). The irregularity of the shape of trees, which can be corroborated in Figure 7, may force the appearance of peripheral regional maxima of reduced height. To avoid the selection of any of these non-representative maxima, the procedure described in Equations (4)–(6) is iteratively performed for increasing h values starting from 1; being the value selected that fulfilling to be the maximum h value, such that the number of segmented maxima contained in IRMAXBINi,h is equal to the number of estimated trees nTi:
I M A X i = I R M A X B I N i , h   | # C C ( I R M A X B I N i , h ) = n T i # C C ( I R M A X B I N i , h + 1 ) < n T i
where for a given binary image f, #CC(f) refers to the number of connected components existing in it. Figure 8 shows the maxima extracted at different h values for the connected component previously proposed in Figure 7.
From this binary representation of the maxima, IMAXi, an ad-hoc image is automatically generated, and eventually used as a marker image to perform the individual segmentation of each tree crown within the initial component cci. This new image, named hereafter ISEEDSi, results from transforming IMAXi such that each of its components is replaced by an all-white pixels-circle with the same area. This is, for each connected component belonging to IMAXi, a corresponding circle-shaped connected component can be found in ISEEDSi, with the same pixel-area and centre of mass. Figure 9 shows an example of this ad-hoc image, obtained in this specific case from the initially connected component presented in Figure 7.
As stated before, ISEEDSi contains the set of markers with which the initially connected component cci is finally segmented at the individual tree crown level. ISEEDSi has been generated according to the hypothesis that the isolated maxima, from which it is obtained, are somehow representative of the size and positions of the overlapped treetops in the component. Therefore, the goal at this point is the achievement of an artifact of separation between the area of influence of each circle-shaped pixel-region, assuming that it will subsequently provide a proper segmentation for the tree crowns. This approach will be exploited in the application of the Watershed transform. Indeed, given a greyscale image, this segmentation tool approaches the intensity levels as values of altitude, so the darkest regions can be seen as catchment basins. By “flooding” them, the lines of water convergence can be detected. These ridgelines are determined by the morphological features of those catchment basins, such as their shape or depth, so they can adequately split the zones of influence of the basins. Hence, the idea is to approach the circled areas as those catchment basins by employing the distance transform formulated in Equation (2). Notwithstanding, in this case, the pixel-regions corresponding to the circles are expected to be the darkest (the deepest) areas in the transformed image, so the distance transform must be redefined to compute here the distance of each pixel to the closest white one. Thus, given ISEEDSi, the grayscale image IDTSEEDSi, over which the Watershed segmentation is subsequently performed, can be mathematically defined as follows:
I D T S E E D S i = D T ( I S E E D S i ) ,
[ D T ( f ) ] ( p ) = min { d ε ( p , q ) | f ( q ) = 1 } ,
Next, the Watershed transform is applied on IDTSEEDSi. The mathematical formulation of this morphological transform can be quite extensive, so the reader is encouraged to consult related literature already published for further study [52,53]. In this sense, it should be underscored that the Watershed transform implementation used in this study is based on an algorithm proposed in [54]. Thus, after its application on IDTSEEDSi, those ridges establishing the boundaries of the basins can be extracted as follows:
[ I R I D G E S i ] ( p ) = { 0   i f   [ W S ( I D T S E E D S i ) ] ( p ) = 0 1   i n   a n y   o t h e r   c a s e
where WS denotes the Watershed transform. Finally, the binary image IRIDGESi is used to separate the tree crown projections contained in the initially connected component cci:
I F I N A L i = min { I C C i , I R I D G E S i }
Figure 10 illustrates this final stage of segmentation based on the Watershed transform.
Each connected component cci containing multiple trees in the original binary image IBIN, is replaced by those in its corresponding processed sub-image IFINALi, resulting in a final binary image, ISEG, where every tree crown appears individually segmented.

2.4.3. Image Postprocessing

When undertaking the last stage of segmentation, whereby initially connected components are properly split according to the ridgelines previously yielded, over segmentation phenomena can occur [55,56,57]. Indeed, given the irregularity in terms of the silhouette of the tree crowns, it is possible for non-significant parts of a plant to be outside the boundaries established by the Watershed transform for its corresponding projection. This would potentially provoke the generation of small connected components wrongly disconnected (according to the 8-connectivity considered) from the main ones representing tree crowns (see Figure 1a). These anomalous pixel-regions are removed by applying a morphological opening followed by a morphological reconstruction to exactly recover the surviving greater connected components corresponding to tree crowns:
I F I N A L P O S T i = R I F I N A L i δ ( γ β ( I F I N A L i ) )
where Rδ refers to the reconstruction by dilation, and γβ is the morphological opening by a disc-shaped structuring element β. The size of β must be large enough to remove the anomalous small connected components, but also smaller than tree crowns to preserve them and enable their reconstruction. The considerable difference in terms of size between the anomalous and tree-crown components meant this decision was not critical, due to a wide range of values being valid. In this context, the radius of β was set to 5.
The explained postprocessing aimed at removing irrelevant components is illustrated in Figure 11.

2.5. Evaluation of the Developed Methodology

As mentioned in Section 2.1, the developed methodology was tested on images acquired in three different crop fields. The approach chosen for assessing its performance was based on the comparison of the binary images resulting from each case study with handmade binary segmentation of the tree crowns. These latter were obtained from ad-hoc images resulting from orthomosaics representatives of each land plot under study, which were in turn generated from the aerial captures collected during the image acquisition stage. Indeed, given the set of multispectral aerial captures collected for each case study, they were processed with photogrammetry software (Pix4DTM Mapper), in order to achieve high-resolution aerial representations of the corresponding land plots. Resulting in five different orthomosaics after processing, one per spectral band of capture, three of them (blue, red-edge, and near-IR) were conveniently combined to finally yield a representative aerial image of the crop, wherein it became rather easy to differentiate the individual tree crowns for a human observer. Then, by using an image editor, a manual delineation of each tree crown was carried out, producing a binary image which, as a gold-standard, was eventually compared at pixel level with the corresponding segmentation provided by the developed methodology. Figure 12 shows a proposal of the colour images generated for each case study.
It should be noted that any other combination of the available orthomosaics may be valid for this purpose. In this sense, the sole objective was to obtain an image comfortable to work with when editing it for manually segmenting the tree crowns using an image editor.
Hence, for each case study, given the binary image resulting from applying the developed methodology, and its corresponding ground-truth image, every pixel in the first one was labelled in accordance with the next four categories:
  • True Positive (TP): that foreground pixel (white pixel) for which its analogous one in the ground-truth image was categorised as tree crown-belonging pixel.
  • False Positive (FP): that foreground pixel which was labelled as a non-tree crown-belonging pixel in the ground-truth image.
  • True Negative (TN): that background pixel (black pixel) such that the corresponding one in the ground-truth image was categorised as a non-tree crown-belonging pixel.
  • False Negative (FN): that background pixel for which its analogous one was labelled as tree crown-belonging pixel in the ground-truth image.
Based on these pixel categories, the following metrics are proposed for the purpose of assessing the accuracy of the segmentation resulting from the application of the developed methodology:
  • Precision (PR): this metric refers to the probability with which a given foreground pixel was correctly categorised. It can be formulated as follows:
P R = t p t p + f p
where tp and fp are, respectively, the number of pixels categorised as TP and the number of them labelled as FP.
  • Recall (RC): it represents the ratio between the number of foreground pixels correctly classified and the whole set of instances of actual tree-belonging pixels in the image. Mathematically:
R C = t p t p + f n
where fn is the number of pixels labelled as FN.
  • F-score: as the harmonic mean of these two metrics just proposed:
F s c o r e = 2 × P R × R C P R + R C
  • Overall Accuracy (OA): it proposes the percentage of pixels correctly classified.
O A = t p + t n t p + t n + f p + f n
  • Intersection-over-Union (IoU): also known as the Jaccard index, this metric signifies the similarity between the resulting segmentation and its corresponding ground-truth image. It can be defined with the following expression:
I o U = | I R E S U L I G T | | I R E S U L I G T |
where IRESUL and IGT are the final segmentation and its binary ground-truth image, respectively. They should be approached as the sets of their corresponding foreground instances, so this metric is computed as the ratio between the number of elements in their intersection and the cardinality of their union.
In addition, a secondary analysis of the results was undertaken, focused on the capacities of the proposed solution for individual tree crown identification, by comparing the connected components yielded in the resulting binary images and the actual tree crowns, represented throughout the connected components in the corresponding ground-truth image. The categories defined below can be instantiated on the basis of this comparison:
  • True Positive, at tree-level (TPt): that connected component in the final segmentation corresponding to an actual tree crown in the binary ground-truth image.
  • False Positive, at tree-level (FPt): that connected component in the final segmentation for which it cannot be found an analogous component, representative of a tree crown, in the ground-truth image.
  • False Negative, at tree-level (FNt): that actual tree crown in the ground-truth segmentation not represented in the final segmentation. In other words, those tree crowns are not detected by the algorithm.
Thus, metrics equivalent to those proposed in Equations (13)–(15) were computed for the defined cases.
On the other hand, it should be noted that methodology was not only tested according to the parametrisation proposed in Section 2.4.1, but also with different configurations regarding the maximum size of the filtering kernel used during preprocessing, in order to evaluate the influence of this value on its performance.

3. Results

As mentioned in Section 2.5, the assessment of the developed methodology was approached in two ways: first, analysing at pixel level the semantic segmentation provided by the final binary images; next, by quantifying its performance in terms of individual tree crown detection. All this, for each of the three different cases considered for testing.
Regarding the capacities of the developed methodology for providing accurate segmentation of individual tree crowns in a binary representation of the field plot, Table 4 summarises the results calculated for the metrics proposed for that purpose.
As it can be observed in Table 4, the results obtained are rather similar in the three cases tested, achieving F-score values over 90% for all of them. Assuming that this latter measure along with IoU may be the most significant metrics to quantify the performance of the developed methodology in terms of pixel classification, the overall lower result corresponded to the lemon-tree orchard case study. Curiously, because of its characteristics regarding plant spatial arrangement and the average tree crown size, the occurrences of overlapping tree crowns were significantly lower in this case. This is in contrast to the orange grove case, where a high level of overlapping between adjacent canopies implied the processing of a considerable number of large connected components, in order to split the individual plant pixel-regions which appeared fused. Indeed, this case study was considered for stress testing because of the extreme conditions of its planting arrangement. However, and despite that, a subtle better accuracy was achieved in the segmentation when compared with the test performed for the lemon-tree orchard, yielding promising results when delineating individual tree crowns in very large, connected components, as can be appreciated in Figure 13. Notwithstanding, the differences are so insignificant for the three crops that they are hardly attributable to the developed methodology, and they may be related to the ground-truth images. The generation of those ground-truth images was performed by hand, involving an arduous task in which it was often difficult to visually differentiate the tree crowns from other kinds of vegetation when manually delineated by a human observer. In this sense, given the nature of the images under study and the complexity of this process, it is not possible to guarantee that the different case studies were represented with the same degree of accuracy by their corresponding ground-truth images.
On the other hand, Table 5 presents the performance of the developed methodology when detecting individual tree crowns projections, and subsequently providing an estimation of the plant population. As can be corroborated, the performance is consistent for the three crops studied, showing only slight deviations.
At this point, it should be remembered that the initial estimation of the trees embedded in the connected components is realised on the basis of a dominant tree crown size within the whole population of tree projections in the binary representation to process. This reference measure, obtained during preprocessing, is used as a basis for partitioning the major diameter of each component, yielding such an estimation. Specifically, this value is obtained from the initial binary image representative of the orchard, computing for each connected component the ellipse with its same normalised second central moment [58]; it is taken from the whole set of minor axis lengths of all the ellipses obtained, as the maximum one. This reference value is increased by 20% before being used for considering certain tolerance. Finally, the number of trees in a given connected component is determined by dividing its major axis length by this value. There are some related phenomena that may occur during this subprocess, that can condition the result of the individual tree crown delineation. Thus, small tree crowns that are actually overlapped may be detected as unique plants. It should be emphasised that these small tree-aggregations happen very occasionally, and mostly at the end of the planting rows, where for some reason there is a certain tendency to plant the trees closer to each other. Likewise, regular adjacent tree crowns with a high degree of canopy overlapping may exhibit this problem. This may provoke an underestimation, which is subsequently translated to the eventual individual segmentation, as the number of predicted trees is used as a threshold when retrieving relevant maxima during the procedure. Consequently, individual tree crown-pixel regions may not be properly split when underestimating the number of crop trees in a given connected component. This will cause an increase in false negatives at the tree level (FNt). Furthermore, the opposite phenomenon can occur, i.e., the overestimation of the number of potential tree projections embedded in each component. This may happen mostly with large connected components, with occurrences of rather big individual tree canopies. This overestimation involves the generation of a larger number of segmentation seeds than should be expected, and a subsequent oversegmentation of the corresponding connected components, which in turn results in the appearance of new false positives (FPt). Both issues are illustrated in Figure 14.
Nevertheless, prior estimations of the trees contained in each connected component were highly and comparably accurate in all the three case studies, thus limiting the specific occurrences of the above-mentioned segmentation problems, and the overall performance of the methodology may be concluded as highly positive.
Finally, Table 6 collects results of applying the whole procedure for different sizes of the filtering-kernel radius, all in a range of plus/minus 30 pixels with respect to the initial parameter-value fixed for each case study (Table 3), for the purpose of assessing its sensibility and impact on the overall performance.
As commented in Section 2.4.1, the maximum size to be reached by the structuring element used for image filtering should ideally be the smallest possible to cover the largest tree crown in the image. However, this parameter is manually selected, and estimating an appropriate value for it can be complex, especially when the number of trees under study is large. As expected, the poorest results were obtained when underestimating the maximum size of the kernel. Indeed, if the largest kernel is not able to cover some of the individual canopies, these may not be properly filtered, so resulting residual maxima, being then part of the background estimation, are wrongly subtracted from the initial DHM intensity matrix, thus affecting the subsequent canopy-ground segmentation, and consequently the overall performance of the methodology. The accuracy of the segmentation at pixel level is compromised by wrongly discarding those regions of the canopies that could not be properly filtered out, increasing the number of pixels classified as FN. On the other hand, tree crowns can hardly be uniquely represented by a single connected component due to oversegmentation, impoverishing the capabilities of the proposed solution to automate plant detection as the number of false positives (FPt) explodes. Indeed, some additional tests were carried out, using extremely small values as the limit for the kernel radius, with the aim of confirming this issue. Thus, for the olive case study, the methodology was applied with a maximum radius of 10 pixels, resulting in a final segmentation with a total number of 8795 connected components. Given an actual population of 3916 trees in the olive groove, this segmentation would reach 4879 FPt in the best-case scenario, implying a very poor performance in terms of individual tree identification. For the purpose of illustrating this issue, Figure 15 shows a canopy-ground segmentation of the orange-tree orchard, carried out with a maximum size for the filtering kernel too small to satisfy the requirements of the methodology (15 pixels), as an extreme example of the problems arising from underestimating this parameter, where it can be visualised some of the issues just commented.
Contrary, when overestimating the maximum size of the kernel it was not observed significant impoverishment of the results. In fact, even when performing stress tests with extremely large radius values, in the order of several hundred pixels, the results remained within similar ranges. However, caution must be exercised in drawing conclusions in this respect, as it may be possible that other case studies with important differences in terms of ground elevation or significant variability of trees’ height would yield inadequate results if too large structuring elements are used. In general, based on the results reported, there is a sweet spot for this parameter, which it may be found around the average radius of the population of tree-crowns projections under study. The criterion suggested by the authors of the procedure, based on reaching a size large enough to cover the largest of the trees, seems appropriate, given that according to the results, the performance of the procedure is mainly compromised when the estimation of the background is carried out with filtering structuring elements that are too small, and not necessarily when this parameter is overestimated.

4. Discussion

Reported results regarding the assessment of the developed methodology, suggest its viability as a tool for the automated counting and delineation of individual tree canopy projections in aerial representations. Furthermore, its performance was comparable in the three cultivars tested, reinforcing its potential applicability in cultivation areas regardless of the type of crop. Notwithstanding, despite these compelling results, there are some points that are noteworthy.
In accordance with results reported, the number of undetected trees after individual canopy delineation is quite low, even in the orange orchard-case study, in which it can be observed large connected components where the pixel subregions to be split appear deeply fused, comprising a scenario prone to the underestimation of the tree crowns enclosed in those components during preprocessing; and by extension, to the occurrence of false negatives (FNt). Despite this, it must not be overlooked that any inaccuracy of the initial estimation of the number of trees contained in each connected component, which is performed at preprocessing stage, may affect the quality of the final segmentation. In this vein and as already mentioned, the number of trees in the aerial binary representation is estimated on the basis of a reference value during this preprocessing stage. This value, after automatically computed, is used to rather split the major diameter of each connected component predicting the number of plants contained. The obtention of this reference value might be approached in different ways for the purpose of slightly improving the overall performance. Indeed, different possibilities were considered during the experimentation process, as individual reference values were calculated from each connected component instead of a global one, but no tangible improvements were observed with any of the alternatives proposed. That said, despite there being little room for improvement in this sense, further research may focus on this topic, exploring new ways to predict the potential number of tree crowns embedded in the components to split.
With regard to the generation of ad-hoc segmentation seeds representatives of each of the pixel subregions to be split, there are also some aspects to be discussed. Indeed, in a near-ideal scenario, with a high level of homogeneity in terms of form and size of the tree crowns and a regular spacing pattern between plants, the location points provided by the methodology used during preprocessing, approached as single white pixels, could properly work as representative of the centre of mass of the corresponding tree projections. They could subsequently be used as proper segmentation seeds from which to determine the boundaries of each subcomponent, thus making it unnecessary to calculate ad-hoc segmentation seeds, as representative artifacts of the spatial influence of each plant within each component. However, this can hardly be expected in a real setting, and occurrences of adjacent overlapping tree crowns with unbalanced sizes are usual. In those cases, the corresponding estimated location points may appear displaced from the hypothetical centres of mass of the tree projections. This is because of the way they are calculated from the aforementioned global reference value, somehow assuming every crop tree has the same size within a given plot. Hence, they cannot properly define the pixel-region influenced by each tree, and subsequently, they cannot ensure an accurate individual tree crown segmentation.
Another issue to be addressed is related to the fact that the solution here proposed does not use vegetation indices, nor any kind of machine learning technology. As already mentioned, one of the main goals pursued during the research has been the applicability of the methodology regardless of the type of tree-based crop to be studied. In this sense, machine learning approaches are based on the development of classification models from training sets, which in turn are fed with representative data of the visual features of the elements to be identified; in this case, tree canopies. So, there is a dependency between the classification models achieved and the kinds of plants with which the corresponding training sets have been generated, thus restricting the scope of application, and hence compromising the generality of the solution. Notwithstanding, this does not invalidate these kinds of approaches. Indeed, proposals based on this form of artificial intelligence for detecting and delineating tree crowns from aerial representations are quite common and effective, always under the assumption that their replicability with different types of crops may be limited. Thus, Csillik et al. [32] detected citrus-trees from multispectral imagery with precision and recall of 94.59% and 97.94%, respectively. Another example of the use of CNNs for individual tree identification is proposed by Ampatzidis and Partel [33]. They undertook the detection of citrus-trees by means of CNN, achieving a precision of 99.90% and a recall of 99.70%, and overall accuracy of 85.5% when segmenting the canopy cover in the grove tested. In this latter work, after tree detection, an alternative approach for estimating the individual canopy projections is proposed based on the normalized difference vegetation index (NDVI). In fact, vegetation indices are a common tool for extracting plant features, included dendrometric characteristics, as it is the case. Thus, Marques et al. [38] used a CIR (Colour-InfraRed) imagery-based vegetation index along with elevation data for individually detecting chestnut-trees, reporting a precision of 99.44%, a recall of 97.80%, and an F-score of 98.61%, and segmenting the canopy coverage with accuracies from 93% to 99% in all case studies addressed. Despite being proven effective in certain scenarios, this kind of approach, based on vegetation indices, can be penalised in terms of performance when a high presence of weed occurs, as their spectral reflectance signature is hardly distinguishable from that of the plants to identify [59,60,61,62]. Indeed, for the purpose of illustrating this issue, vegetation index-based segmentations for two of the case studies addressed in this work are proposed in Figure 16. For each citrus-tree orchard, its corresponding NDVI information, computed from the initial multispectral images, is presented as a grey-scale image, along with the binary representation resulting from the attempt to segment the tree-belonging canopies by exploiting this information. It should be noted that this segmentation has been approached by using statistical thresholding [49] with the corresponding NDVI data.
As can be seen, the result obtained in the case of the lemon-tree orchard is more than acceptable. In fact, this segmentation could potentially be used as a basis for applying the methodology proposed here to delineate each tree crown individually. However, this is not the case for the orange grove, where the presence of a considerable amount of weeds, with a spectral response apparently similar to that offered by the orange trees, makes it difficult to discriminate between tree-canopy projections and soil. This does not disqualify the use vegetative indices within this area of research, since in this case it may be achieved another combination of spectral reflectance measurements that could improve segmentation results for the orange-tree orchard; but it does exemplify the problems that could be involved in using this type of indicators. On this basis, and because of the initial premise of yielding a methodology applicable without any adjustment, in as many diverse scenarios as possible, the use of any kind of vegetation indices was discarded in this work.
In fact, despite the use of vegetation indices and machine learning techniques may be considered as a trend in this field, several studies have addressed the problem of discriminating individual tree-canopies by dispensing this kind of approach, and still reporting significant results. Thus, as previously introduced, Salamí et al. [39] proposed colour and stereo-vision -based segmentations to solve the automated detection of olive trees from aerial captures. They reported an F-score of 99.9% when detecting trees from a dataset with a total of 332 plants. Notwithstanding, it should be noted that the method was tested in olive orchards where plant arrangement guaranteed some spacing between trees, so the problem of overlapping adjacent canopies was not addressed. On the other hand, the study proposed by Ok and Ozdarici-Ok [35] may be of particular interest, given that their proposal was aimed at individually segmenting citrus-tree crowns, by means of morphological image analysis, and it also faced the overlapping canopies issue. After evaluating their methodology, they reported a precision of 91.1%, and a recall of 91.3% when assessing its performance at the pixel level. Positive results were also reported when evaluating the quality of the solution in terms of object detection, achieving an F-score value of 91.2%. In order to establish comparisons with this study, as well as with other works previously discussed, it can be underlined that the methodology here proposed achieved overall precision and recall values of 91.3% and 94.1%, respectively, when assessing the semantic segmentations in terms of pixel classification. Likewise, tests conducted for evaluating individual tree detection outperformed them yielding an F-score of 99.27%. In addition, it should be remarked that, contrary to most of the related published literature, the proposed methodology was conceived to be applicable as a generic solution regardless of the type of crop. Hence, it was tested with three different types of crop trees, obtaining promising results in all the threes case studies addressed.

5. Conclusions

In this paper, an image analysis methodology for the processing of aerial multispectral sensed data collected in crop tree-based intensive cultivation areas has been proposed, generating binary representations of them, and automatically delineating each tree crown that appears, thus providing individual plant canopy cover information.
Different case studies were set out to evaluate the effectiveness of the methodology, and to assess its viability as a solution independent from the type of tree-based crop cultivated, facing diverse scenarios in terms of the variability in the size of the tree population, or the characteristics of the plant’s spatial arrangement. In all cases, including the borderline test realised in an orange-tree grove which presented challenging large intra-row canopy aggregations, results provided accurate segmentations at individual tree crown levels. Based on these results, along with the performance shown in terms of individual tree detection, the framework proposed here can be of broad use to farmers and agronomists for inventorying and monitoring crop trees in cultivated land plots.
Further research could be undertaken to investigate different approaches aimed at detecting and treating those exceptional cases reported, responsible for specific situations in which the under/overestimation of trees embedded in those connected components to process occurs. In addition, despite the fact that the methodology was tested with three different types of tree-based crops, it is assumable that undertaking new case studies focused on other additional crops would reinforce the applicability of the proposed solution, as independent of the kinds of trees farmed.

Author Contributions

Conceptualization, J.M.P. and A.A.; methodology, J.M.P. and A.A.; software, J.M.P., A.A. and B.M.A.-H.; validation, J.M.P., A.A., D.T. and B.M.A.-H.; formal analysis, J.M.P. and A.A.; investigation, J.M.P. and A.A.; resources, J.M.A.; data curation, J.M.P., D.T. and B.M.A.-H.; writing—original draft preparation, J.M.P.; writing—review and editing, A.A. and J.M.A.; visualization, J.M.P. and A.A.; supervision, J.M.A.; project administration, J.M.A.; funding acquisition, J.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

Grant P18-RTJ-4539 funded by the Regional Government of Andalusia and by “PAIDI, Plan Andaluz de Investigación, Desarrollo e Innovación”.

Acknowledgments

Authors would like to thank “Virgen de la Oliva” olive-oil cooperative for generously offering their orchards to conduct this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Mathematical morphology (MM) is a well-established technique for the nonlinear analysis of spatial structures, with remarkable applicability in image processing. Thus, morphological image analysis is rather based on the basic idea of transforming the image to process by its interaction with a simpler structure, called structuring element (SE) and for which its shape and size are known a-priori, resulting in more expressive images according to certain aimed characterization. The morphological operators used in the methodology proposed here are formally described below. In addition, literature related to morphological image analysis can be consulted for deeper study [48,63].
Let f be a greyscale image representing a mapping from a subset D f of 2 , which defines its domain, to a bounded subset of nonnegative integers N 0 :
f : D f 2 N 0 = { 0 ,   , t m a x } .
where t m a x is the maximum reachable value according to the type of data used (e.g., 1 for binary images, 256 for 8-bit images, etc.). Hence, f maps the correspondence by element between two sets, the first being composed of spatially ordered elements ρ (pixels), ρ D f and denoted by a pair of coordinates ( x , y ) , while the second is built with an ordered set of possible values.
With the previous definitions, the intersection of two greyscale images, f and g , is defined as:
f ( ρ ) g ( ρ ) = m i n [ f ( ρ ) , g ( ρ ) ] .
being m i n the minimum operation. The union of those two images can be conversely defined as follows:
f ( ρ ) g ( ρ ) = m a x [ f ( ρ ) , g ( ρ ) ] ,
where m a x denotes the maximum operation.
As already mentioned, morphological operations basically consist of probing the image to process with a certain SE. Mathematically, a SE element can be seen as a binary image β , defining a mapping from a subset D β of 2 to the subset of integer binary values B 0 :
β : D β 2 B 0 = { 0 ,   1 } .
Therefore, β maps the correspondence between the spatially ordered pixels ρ , ρ D β and referenced by a pair of coordinates ( x , y ) , and their values. This mapping must be designed so as to morphologically describe the object to be analysed, with # ( D β ) < # ( D f ) being necessary for its application. Common shapes implemented with SEs include circles, lines, diamonds, etc. In practice, SE is used as a kernel, with its origin in its central pixel. Hence, an image is probed pixel by pixel with this kernel, modifying the pixel in the image at every step, matching with the central pixel of the kernel, according to a given operation.
The morphological erosion of image f by an SE β centred in pixel ρ , is given by the following expression:
[ ε β ( f ) ] ( ρ ) = m i n { f ( ρ + b ) | b D β } .
As a result, pixel ρ in image f is modified with the minimum value of its neighbourhood according to the filter implemented by SE β . Therefore, the effect of erosion is the expansion of darker regions, conditioned by the shape defined in SE.
The dual operator of erosion is dilation. The morphological dilation of image f by a SE β , centred in pixel ρ , is formulated as:
[ δ β ( f ) ] ( ρ ) = m a x { f ( ρ + b ) | b D β } .
Conversely, dilation expands brighter regions in f according to the morphology of SE.
Combining erosion and dilation, two new operators called opening ( γ ) and closing ( φ ) can be defined:
γ β ( f ) = δ β ( ε β ( f ) ) ,
φ β ( f ) = ε β ( δ β ( f ) ) .
Opening removes those brighter objects in the image that can be completely covered by β . Inversely, closing operation removes those darker objects in the image which are completely covered by this SE.
The operators described are complemented by geodesic transformations. The geodesic dilation is the iterative dilation of an image f , called marker, with respect to the mask image g , using a unitary SE. Marker f must be contained within mask g . Mathematically, the operator is defined as:
δ g ( n ) ( f ) = δ g ( 1 ) [ δ g ( n 1 ) ( f ) ] ,   being   δ g ( 1 ) ( f ) = δ ( f ) g , where : # ( D f ) = # ( D g ) ,   and   f ( ρ ) g ( ρ ) ,   ρ D f , D g .
Conversely, the geodesic erosion of image marker f constrained by mask g is:
ε g ( n ) ( f ) = ε g ( 1 ) [ ε g ( n 1 ) ( f ) ] ,   being   ε g ( 1 ) ( f ) = ε ( f ) g , where : # ( D f ) = # ( D g ) ,   and   f ( ρ ) g ( ρ ) ,   ρ D f , D g .
Geodesic dilation and erosion are the basis for building morphological reconstructions. Indeed, the morphological reconstruction by dilation of mask g by marker f , is the geodesic dilation of f constrained by g until idempotence, and it is denoted by:
R g δ ( f ) = δ g ( i ) ( f ) , where i is such that : δ g ( i ) ( f ) = δ g ( i + 1 ) ( f ) .
Dually, the morphological reconstruction by erosion of mask g by marker f , is the geodesic erosion of f constrained by g until idempotence:
R g δ ( f ) = δ g ( i ) ( f ) , where i is such that : ε g ( i ) ( f ) = ε g ( i + 1 ) ( f ) .

References

  1. Mulla, D.J. Twenty five years of remote sensing in precision agriculture: Key advances and remaining knowledge gaps. Biosyst. Eng. 2013, 114, 358–371. [Google Scholar] [CrossRef]
  2. Zhai, Z.; Martínez, J.F.; Beltran, V.; Martínez, N.L. Decision support systems for agriculture 4.0: Survey and challenges. Comput. Electron. Agric. 2020, 170, 105256. [Google Scholar] [CrossRef]
  3. Srinivasan, A. Handbook of Precision Agriculture: Principles and Applications; CRC: New York, NY, USA, 2006; ISBN 9781482277968. [Google Scholar]
  4. Liaghat, S.; Balasundram, S.K. A review: The role of remote sensing in precision agriculture. Am. J. Agric. Biol. Sci. 2010, 5, 50–55. [Google Scholar] [CrossRef] [Green Version]
  5. Wójtowicz, M.; Wójtowicz, A.; Piekarczyk, J. Application of remote sensing methods in agriculture. Commun. Biometry Crop Sci. 2016, 11, 31–50. [Google Scholar]
  6. Aquino, A.; Millan, B.; Diago, M.-P.; Tardaguila, J. Automated early yield prediction in vineyards from on-the-go image acquisition. Comput. Electron. Agric. 2018, 144, 26–36. [Google Scholar] [CrossRef]
  7. Aquino, A.; Ponce, J.M.; Andújar, J.M. Identification of olive fruit, in intensive olive orchards, by means of its morphological structure using convolutional neural networks. Comput. Electron. Agric. 2020, 176, 105616. [Google Scholar] [CrossRef]
  8. Noguera, M.; Millán, B.; Pérez-Paredes, J.J.; Ponce, J.M.; Aquino, A.; Andújar, J.M. A New Low-Cost Device Based on Thermal Infrared Sensors for Olive Tree Canopy Temperature Measurement and Water Status Monitoring. Remote Sens. 2020, 12, 723. [Google Scholar] [CrossRef] [Green Version]
  9. Qureshi, W.S.; Payne, A.; Walsh, K.B.; Linker, R.; Cohen, O.; Dailey, M.N. Machine vision for counting fruit on mango tree canopies. Precis. Agric. 2017, 18, 224–244. [Google Scholar] [CrossRef]
  10. Zhang, J.; Tian, H.; Wang, D.; Li, H.; Mouazen, A.M. A Novel Approach for Estimation of Above-Ground Biomass of Sugar Beet Based on Wavelength Selection and Optimized Support Vector Machine. Remote Sens. 2020, 12, 620. [Google Scholar] [CrossRef] [Green Version]
  11. Zhao, Y.; Gong, L.; Zhou, B.; Huang, Y.; Liu, C. Detecting tomatoes in greenhouse scenes by combining AdaBoost classifier and colour analysis. Biosyst. Eng. 2016, 148, 127–137. [Google Scholar] [CrossRef]
  12. Colwell, R.N. Determining the prevalence of certain cereal crop diseases by means of aerial photography. Hilgardia 1956, 26, 223–286. [Google Scholar] [CrossRef] [Green Version]
  13. Bauer, M.; Cipra, J. Identification of Agricultural Crops by Computer Processing of ERTS MSS Data. LARS Tech. Rep. 1973, 20, 205–212. [Google Scholar]
  14. Hoffman, R.O.; Edwards, D.M.; Eucker, C.C. Identifying and Measuring Crop Type Using Satellite Imagery. Trans. ASAE 1976, 19, 1066–1070. [Google Scholar] [CrossRef]
  15. Colomina, I.; Molina, P. Unmanned aerial systems for photogrammetry and remote sensing: A review. ISPRS J. Photogramm. Remote Sens. 2014, 92, 79–97. [Google Scholar] [CrossRef] [Green Version]
  16. Hunt, E.R.; Daughtry, C.S.T. What good are unmanned aircraft systems for agricultural remote sensing and precision agriculture? Int. J. Remote Sens. 2018, 39, 5345–5376. [Google Scholar] [CrossRef] [Green Version]
  17. Yang, G.; Liu, J.; Zhao, C.; Li, Z.; Huang, Y.; Yu, H.; Xu, B.; Yang, X.; Zhu, D.; Zhang, X.; et al. Unmanned Aerial Vehicle Remote Sensing for Field-Based Crop Phenotyping: Current Status and Perspectives. Front. Plant Sci. 2017, 8, 1111. [Google Scholar] [CrossRef]
  18. Zhang, C.; Kovacs, J.M. The application of small unmanned aerial systems for precision agriculture: A review. Precis. Agric. 2012, 13, 693–712. [Google Scholar] [CrossRef]
  19. Hassler, S.C.; Baysal-Gurel, F. Unmanned Aircraft System (UAS) Technology and Applications in Agriculture. Agronomy 2019, 9, 618. [Google Scholar] [CrossRef] [Green Version]
  20. Gevaert, C.M.; Tang, J.; García-Haro, F.J.; Suomalainen, J.; Kooistra, L. Combining hyperspectral UAV and multispectral Formosat-2 imagery for precision agriculture applications. In Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing (WHISPERS); IEEE: Piscataway, NJ, USA, 2014. [Google Scholar] [CrossRef]
  21. Matese, A.; Toscano, P.; Di Gennaro, S.F.; Genesio, L.; Vaccari, F.P.; Primicerio, J.; Belli, C.; Zaldei, A.; Bianconi, R.; Gioli, B. Intercomparison of UAV, Aircraft and Satellite Remote Sensing Platforms for Precision Viticulture. Remote Sens. 2015, 7, 2971–2990. [Google Scholar] [CrossRef] [Green Version]
  22. He, F.; Xiong, W.; Habib, A. A Structure-from-Motion Approach Using UAV-based Imagery for Precision Agriculture Applications. In Proceedings of the 10th International Conference on Mobile Mapping Technology, Cairo, Egypt, 6–8 May 2017. [Google Scholar]
  23. Chapman, S.C.; Merz, T.; Chan, A.; Jackway, P.; Hrabar, S.; Dreccer, M.F.; Holland, E.; Zheng, B.; Ling, T.J.; Jimenez-Berni, J. Pheno-Copter: A Low-Altitude, Autonomous Remote-Sensing Robotic Helicopter for High-Throughput Field-Based Phenotyping. Agronomy 2014, 4, 279–301. [Google Scholar] [CrossRef] [Green Version]
  24. Marino, S.; Alvino, A. Detection of Spatial and Temporal Variability of Wheat Cultivars by High-Resolution Vegetation Indices. Agronomy 2019, 9, 226. [Google Scholar] [CrossRef] [Green Version]
  25. Zaman-Allah, M.; Vergara, O.; Araus, J.L.; Tarekegne, A.; Magorokosho, C.; Zarco-Tejada, P.J.; Hornero, A.; Albà, A.H.; Das, B.; Craufurd, P.; et al. Unmanned aerial platform-based multi-spectral imaging for field phenotyping of maize. Plant Methods 2015, 11, 1–10. [Google Scholar] [CrossRef] [Green Version]
  26. Roth, L.; Streit, B. Predicting cover crop biomass by lightweight UAS-based RGB and NIR photography: An applied photogrammetric approach. Precis. Agric. 2018, 19, 93–114. [Google Scholar] [CrossRef] [Green Version]
  27. Jin, X.; Liu, S.; Baret, F.; Hemerlé, M.; Comar, A. Estimates of plant density of wheat crops at emergence from very low altitude UAV imagery. Remote Sens. Environ. 2017, 198, 105–114. [Google Scholar] [CrossRef] [Green Version]
  28. Aparna, P.; Ramachandra, H.; Sounder, H.; Harshita, M.P.; Nandkishore, K.; Vinod, P.V. CNN Based Technique for Automatic Tree Counting Using Very High Resolution Data. In Proceedings of the 2018 International Conference on Design Innovations for 3Cs Compute Communicate Control (ICDI3C), Bangalore, India, 25–26 April 2018; IEEE: Piscataway, NJ, USA, 2018; pp. 127–129. [Google Scholar]
  29. Rosell, J.R.; Sanz, R. A review of methods and applications of the geometric characterization of tree crops in agricultural activities. Comput. Electron. Agric. 2012, 81, 124–141. [Google Scholar] [CrossRef] [Green Version]
  30. Ma, Q.; Su, Y.; Guo, Q. Comparison of Canopy Cover Estimations from Airborne LiDAR, Aerial Imagery, and Satellite Imagery. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 4225–4236. [Google Scholar] [CrossRef]
  31. Zortea, M.; Macedo, M.M.G.; Mattos, A.B.; Ruga, B.C.; Gemignani, B.H. Automatic Citrus Tree Detection from UAV Images based on Convolutional Neural Networks. In Proceedings of the 31th Sibgrap/WIA—Conference on Graphics, Patterns and Images, Foz do Iguaçu, Brazil, 29 October–1 November 2018. [Google Scholar]
  32. Csillik, O.; Cherbini, J.; Johnson, R.; Lyons, A.; Kelly, M. Identification of citrus trees from unmanned aerial vehicle imagery using convolutional neural networks. Drones 2018, 2, 39. [Google Scholar] [CrossRef] [Green Version]
  33. Ampatzidis, Y.; Partel, V. UAV-Based High Throughput Phenotyping in Citrus Utilizing Multispectral Imaging and Artificial Intelligence. Remote Sens. 2019, 11, 410. [Google Scholar] [CrossRef] [Green Version]
  34. Recio, J.A.; Hermosilla, T.; Ruiz, L.A.; Palomar, J. Automated extraction of tree and plot-based parameters in citrus orchards from aerial images. Comput. Electron. Agric. 2013, 90, 24–34. [Google Scholar] [CrossRef]
  35. Ok, A.O.; Ozdarici-Ok, A. 2-D delineation of individual citrus trees from UAV-based dense photogrammetric surface models. Int. J. Digit. Earth 2017, 11, 583–608. [Google Scholar] [CrossRef]
  36. Modica, G.; De Luca, G.; Messina, G.; Praticò, S. Comparison and assessment of different object-based classifications using machine learning algorithms and UAVs multispectral imagery: A case study in a citrus orchard and an onion crop. Eur. J. Remote Sens. 2021, 54, 431–460. [Google Scholar] [CrossRef]
  37. Dong, X.; Zhang, Z.; Yu, R.; Tian, Q.; Zhu, X. Extraction of Information about Individual Trees from High-Spatial-Resolution UAV-Acquired Images of an Orchard. Remote Sens. 2020, 12, 133. [Google Scholar] [CrossRef] [Green Version]
  38. Marques, P.; Pádua, L.; Adão, T.; Hruška, J.; Peres, E.; Sousa, A.; Sousa, J.J. UAV-Based Automatic Detection and Monitoring of Chestnut Trees. Remote Sens. 2019, 11, 855. [Google Scholar] [CrossRef] [Green Version]
  39. Salamí, E.; Gallardo, A.; Skorobogatov, G.; Barrado, C. On-the-Fly Olive Trees Counting Using a UAS and Cloud Services. Remote Sens. 2019, 11, 316. [Google Scholar] [CrossRef] [Green Version]
  40. Jiang, H.; Chen, S.; Li, D.; Wang, C.; Yang, J. Papaya Tree Detection with UAV Images Using a GPU-Accelerated Scale-Space Filtering Method. Remote Sens. 2017, 9, 721. [Google Scholar] [CrossRef] [Green Version]
  41. Johansen, K.; Raharjo, T.; McCabe, M. Using Multi-Spectral UAV Imagery to Extract Tree Crop Structural Properties and Assess Pruning Effects. Remote Sens. 2018, 10, 854. [Google Scholar] [CrossRef] [Green Version]
  42. Diez, C.M.; Moral, J.; Cabello, D.; Morello, P.; Rallo, L.; Barranco, D. Cultivar and tree density as key factors in the long-term performance of super high-density olive orchards. Front. Plant Sci. 2016, 7, 1226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Menzel, C.M.; Le Lagadec, M.D. Increasing the productivity of avocado orchards using high-density plantings: A review. Sci. Hortic. 2014, 177, 21–36. [Google Scholar] [CrossRef]
  44. Majid, I.; Khalil, A.K.; Nazir, N. Economic Analysis of High Density Orchards. Int. J. Adv. Res. Sci. Eng. 2018, 7, 821–829. [Google Scholar] [CrossRef]
  45. Sarabia, R.; Aquino, A.; Ponce, J.M.; López, G.; Andújar, J.M. Automated Identification of Crop Tree Crowns from UAV Multispectral Imagery by Means of Morphological Image Analysis. Remote Sens. 2020, 12, 748. [Google Scholar] [CrossRef] [Green Version]
  46. Instituto Geográfico Nacional. Available online: http://www.ign.es/web/ign/portal (accessed on 19 November 2021).
  47. Shepard, D. A two-dimensional interpolation function for irregularly-spaced data. In Proceedings of the 1968 23rd ACM National Conference, New York, NY, USA, 27–29 August 1968; ACM Press: New York, NY, USA, 1968; pp. 517–524. [Google Scholar]
  48. Soille, P. Morphological Image Analysis: Principles and Applications; Springer: Berlin/Heidelberg, Germany, 2004; ISBN 9783662050880. [Google Scholar]
  49. Otsu, N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans. Syst. Man. Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef] [Green Version]
  50. Rosenfeld, A. Connectivity in Digital Pictures. J. ACM 1970, 17, 146–160. [Google Scholar] [CrossRef]
  51. Maurer, C.R.; Qi, R.; Raghavan, V. A linear time algorithm for computing exact Euclidean distance transforms of binary images in arbitrary dimensions. IEEE Trans. Pattern Anal. Mach. Intell. 2003, 25, 265–270. [Google Scholar] [CrossRef]
  52. Beucher, S.; Lantuejoul, C. Use of watersheds in contour detection. In Proceedings of the International Workshop on Image Processing, Real-Time Edge and Motion Detection, Rennes, France, 17–21 September 1979; pp. 17–21. [Google Scholar]
  53. Beucher, S.; Meyer, F. The Morphological Approach to Segmentation: The Watershed Transformation. In Mathematical Morphology in Image Processing, 1st ed.; Dougherty, E., Ed.; CRC: New York, NY, USA, 1992; Volume I, pp. 433–481. ISBN 9780824787240. [Google Scholar]
  54. Meyer, F. Topographic distance and watershed lines. Signal Process. 1994, 38, 113–125. [Google Scholar] [CrossRef]
  55. Gies, V.; Bernard, T.M. Statistical solution to watershed over-segmentation. Proc. Int. Conf. Image Process. ICIP 2004, 3, 1863–1866. [Google Scholar] [CrossRef]
  56. Frucci, M. From Segmentation to Binarization of Gray-level Images. J. Pattern Recognit. Res. 2008, 1, 1–13. [Google Scholar] [CrossRef]
  57. Najman, L.; Schmitt, M. Geodesic saliency of watershed contours and hierarchical segmentation. IEEE Trans. Pattern Anal. Mach. Intell. 1996, 18, 1163–1173. [Google Scholar] [CrossRef] [Green Version]
  58. Jain, A.K. Fundamentals of Digital Image Processing; Prentice Hall: Upper Saddle River, NJ, USA, 1989; ISBN 0133361659. [Google Scholar]
  59. Pérez, A.J.; López, F.; Benlloch, J.V.; Christensen, S. Colour and shape analysis techniques for weed detection in cereal fields. Comput. Electron. Agric. 2000, 25, 197–212. [Google Scholar] [CrossRef]
  60. Thorp, K.R.; Tian, L.F. A Review on Remote Sensing of Weeds in Agriculture. Precis. Agric. 2004, 5, 477–508. [Google Scholar] [CrossRef]
  61. López-Granados, F. Weed detection for site-specific weed management: Mapping and real-time approaches. Weed Res. 2011, 51, 1–11. [Google Scholar] [CrossRef] [Green Version]
  62. Torres-Sánchez, J.; Peña-Barragán, J.M.; Gómez-Candón, D.; De Castro, A.I.; López-Granados, F. Imagery from unmanned aerial vehicles for early site specific weed management. Precis. Agric. 2013, 13, 193–199. [Google Scholar] [CrossRef]
  63. Serra, J. Image Analysis and Mathematical Morphology; Academic Press Inc.: Cambridge, MA, USA, 1982; Volume I, ISBN 9780126372427. [Google Scholar]
Figure 1. Aerial images of the case study sites: (a) lemon-tree orchard; (b) orange-tree orchard; (c) olive grove. Note in each case the corresponding plot under study enclosed in yellow.
Figure 1. Aerial images of the case study sites: (a) lemon-tree orchard; (b) orange-tree orchard; (c) olive grove. Note in each case the corresponding plot under study enclosed in yellow.
Agronomy 12 00043 g001
Figure 2. Aerial image acquisition equipment: (a) DJI Matrice 100-based setup; (b) DJI Phantom 3 Advance-based setup.
Figure 2. Aerial image acquisition equipment: (a) DJI Matrice 100-based setup; (b) DJI Phantom 3 Advance-based setup.
Agronomy 12 00043 g002
Figure 3. Examples of aerial images acquired for the case study corresponding to the orange-tree orchard, captured in different spectral bands: (a) Blue band; (b) Green band; (c) Red band; (d) NIR band; (e) Red-Edge band.
Figure 3. Examples of aerial images acquired for the case study corresponding to the orange-tree orchard, captured in different spectral bands: (a) Blue band; (b) Green band; (c) Red band; (d) NIR band; (e) Red-Edge band.
Agronomy 12 00043 g003
Figure 4. Flowchart of the methodology proposed by Sarabia et al. [45], aimed at generating a binary image from aerial multispectral captures of a land plot, where canopy coverage is shown segmented from the ground, and the set of location points of the crop trees in the plot is generated.
Figure 4. Flowchart of the methodology proposed by Sarabia et al. [45], aimed at generating a binary image from aerial multispectral captures of a land plot, where canopy coverage is shown segmented from the ground, and the set of location points of the crop trees in the plot is generated.
Agronomy 12 00043 g004
Figure 5. Binary image resulting from pre-processing aerial captures of the orange grove-case study, according to the methodology proposed by Sarabia et al. [45], represented together with the estimated tree location points (red circles). Note in the zoomed area, highlighted in red, the occurrences of overlapping tree crowns.
Figure 5. Binary image resulting from pre-processing aerial captures of the orange grove-case study, according to the methodology proposed by Sarabia et al. [45], represented together with the estimated tree location points (red circles). Note in the zoomed area, highlighted in red, the occurrences of overlapping tree crowns.
Agronomy 12 00043 g005
Figure 6. Flowchart of the morphological image analysis algorithm for individually segmenting the crowns of the trees in the image aimed at processing.
Figure 6. Flowchart of the morphological image analysis algorithm for individually segmenting the crowns of the trees in the image aimed at processing.
Agronomy 12 00043 g006
Figure 7. (a) Sub-image of the initial binary image, IBIN; (b) the resulting image of enclosing the connected component highlighted in red in (a) in its corresponding bounding box; (c) intensity matrix, IDTi, resulting from the computing of the Euclidean distance transform of the binary image represented in (b).
Figure 7. (a) Sub-image of the initial binary image, IBIN; (b) the resulting image of enclosing the connected component highlighted in red in (a) in its corresponding bounding box; (c) intensity matrix, IDTi, resulting from the computing of the Euclidean distance transform of the binary image represented in (b).
Agronomy 12 00043 g007
Figure 8. Maxima extracted from the distance-transform image corresponding to the connected component proposed in Figure 7, computed at different heights h: (a) h = 1; (b) h = 5; (c) h = 10; (d) h = 15; (e) h = 16. Note how the proper maxima are obtained for a height of 15, since the number of trees previously estimated for the component is 2, and this value for h is the highest with which it is achieved a maxima segmentation with 2 connected components.
Figure 8. Maxima extracted from the distance-transform image corresponding to the connected component proposed in Figure 7, computed at different heights h: (a) h = 1; (b) h = 5; (c) h = 10; (d) h = 15; (e) h = 16. Note how the proper maxima are obtained for a height of 15, since the number of trees previously estimated for the component is 2, and this value for h is the highest with which it is achieved a maxima segmentation with 2 connected components.
Agronomy 12 00043 g008
Figure 9. Segmentation markers for finally splitting the connected component proposed in Figure 7: (a) binary image corresponding to the representative maxima of those trees enclosed in the component, IMAXi; (b) binary image generated from (a), where each connected component is replaced by a disc-shaped pixel-region with the same area and location, and which comprises the final segmentation markers.
Figure 9. Segmentation markers for finally splitting the connected component proposed in Figure 7: (a) binary image corresponding to the representative maxima of those trees enclosed in the component, IMAXi; (b) binary image generated from (a), where each connected component is replaced by a disc-shaped pixel-region with the same area and location, and which comprises the final segmentation markers.
Agronomy 12 00043 g009
Figure 10. Final segmentation at individual tree crown level of the connected component cci: (a) distance transform of ISEEDSi, based on the Euclidean distance of each pixel to the closest white one; (b) Watershed transform of (a); (c) ridge line obtained from the Watershed transform-image (b); (d) the image resulting from splitting the initial component cci with the ridge line in (c). It should be noted that, for the sake of facilitating its visualization, the image display range has been modified for subimages (a,b), conveniently increasing their contrast.
Figure 10. Final segmentation at individual tree crown level of the connected component cci: (a) distance transform of ISEEDSi, based on the Euclidean distance of each pixel to the closest white one; (b) Watershed transform of (a); (c) ridge line obtained from the Watershed transform-image (b); (d) the image resulting from splitting the initial component cci with the ridge line in (c). It should be noted that, for the sake of facilitating its visualization, the image display range has been modified for subimages (a,b), conveniently increasing their contrast.
Agronomy 12 00043 g010
Figure 11. Postprocessing aimed at filtering anomalous small connected components generated after segmentation: (a) initial subimage ICCi of the connected component to be segmented; (b) ridge line-image IRIDGESi obtained with the Watershed transform for segmenting (a); (c) resulting segmentation after splitting the component in (a) using (b) -note in the zoomed area the generation of anomalous disconnected artifacts; (d) segmentation-image IFINALPOSTi obtained after postprocessing-note in the zoomed area that all elements are connected, considering 8-connectivity, and how the shape of the crowns are not affected as a consequence of the postprocessing.
Figure 11. Postprocessing aimed at filtering anomalous small connected components generated after segmentation: (a) initial subimage ICCi of the connected component to be segmented; (b) ridge line-image IRIDGESi obtained with the Watershed transform for segmenting (a); (c) resulting segmentation after splitting the component in (a) using (b) -note in the zoomed area the generation of anomalous disconnected artifacts; (d) segmentation-image IFINALPOSTi obtained after postprocessing-note in the zoomed area that all elements are connected, considering 8-connectivity, and how the shape of the crowns are not affected as a consequence of the postprocessing.
Agronomy 12 00043 g011
Figure 12. Ad-hoc colour representations of each case study, obtained from the combination of the Blue, Red-Edge, and NIR spectral bands, and used to generate the binary ground-truth images. (a) Lemon-tree orchard; (b) orange-tree orchard; (c) olive grove.
Figure 12. Ad-hoc colour representations of each case study, obtained from the combination of the Blue, Red-Edge, and NIR spectral bands, and used to generate the binary ground-truth images. (a) Lemon-tree orchard; (b) orange-tree orchard; (c) olive grove.
Agronomy 12 00043 g012
Figure 13. Sub-images of the final segmentations obtained for each considered case study: (a) lemon-tree orchard; (b) orange-tree orchard; (c) olive grove.
Figure 13. Sub-images of the final segmentations obtained for each considered case study: (a) lemon-tree orchard; (b) orange-tree orchard; (c) olive grove.
Agronomy 12 00043 g013
Figure 14. Problems with the individual tree crown delineation: (a) non-segmentation of adjacent tree-crown projections due to underestimation of the number of trees contained in the connected component; (b) oversegmentation of adjacent tree-crown projections due to overestimation of the number of trees contained in the connected component. Note highlighted in red those subcomponents involved with the issues reported.
Figure 14. Problems with the individual tree crown delineation: (a) non-segmentation of adjacent tree-crown projections due to underestimation of the number of trees contained in the connected component; (b) oversegmentation of adjacent tree-crown projections due to overestimation of the number of trees contained in the connected component. Note highlighted in red those subcomponents involved with the issues reported.
Agronomy 12 00043 g014
Figure 15. Resulting segmentation for the orange-tree case study, executed with a maximum size of 15 pixels for the radius of the filtering kernel used during preprocessing for estimating the background of the corresponding DHM-greyscale image.
Figure 15. Resulting segmentation for the orange-tree case study, executed with a maximum size of 15 pixels for the radius of the filtering kernel used during preprocessing for estimating the background of the corresponding DHM-greyscale image.
Agronomy 12 00043 g015
Figure 16. NDVI-based segmentation of the canopy coverage in the citrus-tree orchards: (a) lemon-tree orchard NDVI data, represented as a greyscale image; (b) binary representation of the lemon-tree orchard resulting from NDVI-based segmentation; (c) orange-tree orchard NDVI data, represented as a greyscale image; (d) binary representation of the orange-tree orchard resulting from NDVI-based segmentation.
Figure 16. NDVI-based segmentation of the canopy coverage in the citrus-tree orchards: (a) lemon-tree orchard NDVI data, represented as a greyscale image; (b) binary representation of the lemon-tree orchard resulting from NDVI-based segmentation; (c) orange-tree orchard NDVI data, represented as a greyscale image; (d) binary representation of the orange-tree orchard resulting from NDVI-based segmentation.
Agronomy 12 00043 g016
Table 1. Main features of the MicaSense RedEdge-MTM multispectral camera.
Table 1. Main features of the MicaSense RedEdge-MTM multispectral camera.
Spectral BandsCentre Wavelength (nm)Bandwidth (nm)
Blue47520
Green56020
Red66810
Red-Edge71710
NIR84040
Ground Sample Distance (GSD)8 cm per pixel (per band) at 120 m AGL 1
Max Capture Rate1 capture per second (all bands), 12-bit RAW
Field of View47.2° HFOV 2
Imager Resolution1280 × 960 pixels
1 Above Ground Level. 2 Horizontal Field of View.
Table 2. Main features of the UAVs used for image acquisition.
Table 2. Main features of the UAVs used for image acquisition.
DJITM Matrice 100DJITM Phantom 3
Diagonal Wheelbase650 mm350 mm
Max Take-off Weight3600 g1280 g 1
Max Speed17 m/s (GPS mode, no payload, no wind)
22 m/s (ATTI mode, no payload, no wind)
16 m/s (ATTI mode, no payload, no wind)
Max Wind Resistance10 m/s10 m/s
Operating Temperature10° to 40°0° to 40°
1 Value not provided by the manufacturer. The indicated one corresponds to the weight of the aircraft, batteries and propellers included.
Table 3. Configuration of the parameter values of the preprocessing procedure, taken for each case study.
Table 3. Configuration of the parameter values of the preprocessing procedure, taken for each case study.
Case StudyFiltering-Kernel Maximum Radius-Size, rβ (pixels/cm 1)Relevant Maxima-Threshold, h (m)
Lemon-tree42.5/2041.5
Orange-tree55/2641.1
Olive70/3361
1 Based on a GSD of 4.8 cm.
Table 4. Accuracy of the semantic segmentation provided by the developed methodology.
Table 4. Accuracy of the semantic segmentation provided by the developed methodology.
Case StudyPR 1RC 2F-score 3OA 4IoU 5
Lemon-tree0.919930.897500.908580.989860.83247
Orange-tree0.931690.942800.937210.982910.88185
Olive0.889180.982490.933510.981850.87530
1 Precision; PR = tp/(tp + fp). 2 Recall; RC = tp/(tp + fn). 3 F-score =2 × (PR × RC)/(PR + RC). 4 Overall Accuracy; OA = (tp + tn)/(tp + tn + fp + fn). 5 Intersection-over-Union; IoU = |IRESULIGT|/|IRESULIGT|.
Table 5. Results for the assessment of the developed methodology in terms of individual crop tree identification.
Table 5. Results for the assessment of the developed methodology in terms of individual crop tree identification.
Case Studytpt 1fpt 2fnt 3PRt 4RCt 5F-scoret 6
Lemon-tree5480410.992750.99636
Orange-tree7582230.997370.970550.98377
Olive39065100.998720.997440.99808
1 Number of CCs categorised as TPt. 2 Number of CCs categorised as FPt. 3 Number of CCs categorised as FNt. 4 Precision; PRt = tpt/(tpt + fpt). 5 Recall; RCt = tpt/(tpt + fnt). 6 F-scoret = 2 × (PR × RC)/(PR + RC).
Table 6. Results of the developed methodology in terms of accuracy of semantic segmentation and individual crop-tree identification, regarding the maximum radius size of the preprocessing filtering kernel, rβ.
Table 6. Results of the developed methodology in terms of accuracy of semantic segmentation and individual crop-tree identification, regarding the maximum radius size of the preprocessing filtering kernel, rβ.
Semantic Segmentation AccuracyIndividual Tree Detection
Case Studyrβ (pixels)PR 1RC 2F-score 3OA 4IoU 5PRt 6RCt 7F-scoret 8
Lemon-tree12.50.903570.304660.455670.959150.295060.644980.862320.73798
22.50.904980.569470.699050.972480.537340.752110.967390.84627
32.50.920300.896360.908170.989830.8317910.992750.99636
42.50.919930.89750.908580.989860.8324710.992750.99636
52.50.919930.897510.908580.989860.8324810.992750.99636
62.50.919800.897670.908600.989860.8325010.992750.99636
72.50.919650.897750.908570.989860.8324510.992750.99636
Orange-tree250.890270.39770.549800.911870.379120.609450.941100.73981
350.913340.563250.696790.933670.534670.967410.950060.95866
450.931700.927550.929620.981000.868500.992170.973110.98255
550.931690.942800.937210.982910.881850.997370.970550.98378
650.931670.942840.937220.982910.881860.997370.970550.98378
750.931520.943350.937400.982950.882170.997370.970550.98378
850.932450.941310.936860.982830.881210.996070.973110.98446
Olive400.893260.925400.909050.975990.833260.968500.997190.98264
500.890900.980050.933350.981850.875030.998720.997190.99796
600.889180.982490.933510.981850.875300.998470.997450.99796
700.889180.982490.933510.981850.875300.998720.997450.99808
800.888680.982740.933350.981800.875020.998720.997700.99821
900.888100.983040.933160.981740.874700.999230.997960.99859
1000.888230.982960.933200.981750.874760.999230.997700.99847
1PR = tp/(tp + fp). 2 RC = tp/(tp + fn). 3 F1 = 2 × (PR × RC)/(PR + RC). 4 OA = (tp + tn)/(tp + tn + fp + fn). 5 IIoU = |IRESULIGT|/|IRESULIGT|. 6 PRt = tpt/(tpt + fpt). 7 RCt = tpt/(tpt + fnt). 8 F-scoret = 2 × (PRt × RCt)/(PRt + RCt).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ponce, J.M.; Aquino, A.; Tejada, D.; Al-Hadithi, B.M.; Andújar, J.M. A Methodology for the Automated Delineation of Crop Tree Crowns from UAV-Based Aerial Imagery by Means of Morphological Image Analysis. Agronomy 2022, 12, 43. https://doi.org/10.3390/agronomy12010043

AMA Style

Ponce JM, Aquino A, Tejada D, Al-Hadithi BM, Andújar JM. A Methodology for the Automated Delineation of Crop Tree Crowns from UAV-Based Aerial Imagery by Means of Morphological Image Analysis. Agronomy. 2022; 12(1):43. https://doi.org/10.3390/agronomy12010043

Chicago/Turabian Style

Ponce, Juan Manuel, Arturo Aquino, Diego Tejada, Basil Mohammed Al-Hadithi, and José Manuel Andújar. 2022. "A Methodology for the Automated Delineation of Crop Tree Crowns from UAV-Based Aerial Imagery by Means of Morphological Image Analysis" Agronomy 12, no. 1: 43. https://doi.org/10.3390/agronomy12010043

APA Style

Ponce, J. M., Aquino, A., Tejada, D., Al-Hadithi, B. M., & Andújar, J. M. (2022). A Methodology for the Automated Delineation of Crop Tree Crowns from UAV-Based Aerial Imagery by Means of Morphological Image Analysis. Agronomy, 12(1), 43. https://doi.org/10.3390/agronomy12010043

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop