# Using a Structural Root System Model to Evaluate and Improve the Accuracy of Root Image Analysis Pipelines

Guillaume Lobet, Iko T. Koevoets, Manuel Noll, Patrick E. Meyer, Pierre

Tocquin, Loic Pagès, Claire Perilleux

► To cite this version:

Guillaume Lobet, Iko T. Koevoets, Manuel Noll, Patrick E. Meyer, Pierre Tocquin, et al. Using a Structural Root System Model to Evaluate and Improve the Accuracy of Root Image Analysis Pipelines. Frontiers in Plant Science, 2017, 8, 10.3389/fpls.2017.00447. [https://doi.org/10.3389/fpls.2017.00447](https://doi.org/10.3389/fpls.2017.00447). [arXiv:1705.03047](https://arxiv.org/abs/1705.03047)

HAL Id: hal-01517235

https://hal.science/hal-01517235

Submitted on 2 May 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Using a Structural Root System Model to Evaluate and Improve the Accuracy of Root Image Analysis Pipelines

Guillaume Lobet¹,²*†, Iko T. Koevoets³†, Manuel Noll¹†, Patrick E. Meyer¹, Pierre Tocquin¹, Loïc Pagès⁴ and Claire Périlleux¹

¹InBioS-PhytoSYSTEMS, University of Liège, Liège, Belgium, ²Institut für Bio-und Geowissenschaften: Agrosphere, Forschungszentrum Jülich, Jülich, Germany, ³Plant Cell Biology, Swammerdam Institute for Life Sciences, University of Amsterdam, Amsterdam, Netherlands, ⁴INRA, Centre d'Avignon, UR 1115 PSH, Avignon, France

Root system analysis is a complex task, often performed with fully automated image analysis pipelines. However, the outcome is rarely verified by ground-truth data, which might lead to underestimated biases. We have used a root model, ArchiSimple, to create a large and diverse library of ground-truth root system images (10,000). For each image, three levels of noise were created. This library was used to evaluate the accuracy and usefulness of several image descriptors classically used in root image analysis softwares. Our analysis highlighted that the accuracy of the different traits is strongly dependent on the quality of the images and the type, size, and complexity of the root systems analyzed. Our study also demonstrated that machine learning algorithms can be trained on a synthetic library to improve the estimation of several root system traits. Overall, our analysis is a call to caution when using automatic root image analysis tools. If a thorough calibration is not performed on the dataset of interest, unexpected errors might arise, especially for large and complex root images. To facilitate such calibration, both the image library and the different codes used in the study have been made available to the community.

**Keywords:** image analysis, root structural model, benchmarking, image library, machine learning

## INTRODUCTION

Roots are of utmost importance in the life of plants and hence selection on root systems represents great promise for improving crop tolerance to biotic and abiotic stresses (as reviewed in Koevoets et al., 2016). As such, their quantification is a challenge in many research projects. This quantification is usually two-fold. The first step consists in acquiring images of the root system, either using classic imaging techniques (CCD cameras) or more specialized ones (microCT, X-Ray, fluorescence,...). The next step is to analyse the pictures to extract meaningful descriptors of the root system.

To paraphrase the famous Belgian surrealist painter, René Magritte: “**Figure 1A** is not a root system.” **Figure 1A** is an image of a root system and that distinction is important. An image is indeed a two-dimensional representation of an object, which is usually three-dimensional.

## OPEN ACCESS

**Edited by:**
Alexander Bucksch,
University of Georgia, USA

**Reviewed by:**
Christopher N. Topp,
Donald Danforth Plant Science Center,
USA
Tony Pridmore,
University of Nottingham, UK

*Correspondence:
Guillaume Lobet
g.lobet@fz-juelich.de

† These authors have contributed
equally to this work.

**Specialty section:**
This article was submitted to
Plant Biophysics and Modeling,
a section of the journal
Frontiers in Plant Science

**Received:** 14 September 2016

**Accepted:** 15 March 2017

**Published:** 03 April 2017

**Citation:**
Lobet G, Koevoets IT, Noll M, Meyer PE, Tocquin P, Pagès L and Périlleux C (2017) Using a Structural Root System Model to Evaluate and Improve the Accuracy of Root Image Analysis Pipelines. *Front. Plant Sci.* 8:447. doi: 10.3389/fpls.2017.00447
<figure><img src="image_2.png" /><figcaption>**FIGURE 1** | (A) Image of a 2-week old maize root system grown in rhizotron. (B) Close-up showing overlapping roots. (C) Close-up showing crossing roots.</figcaption></figure>

Nowadays, measurements are generally not performed on the root systems themselves, but on the images, and this raises some issues.

Image analysis is the acquisition of traits (or descriptors) describing the objects contained in a particular image. In a perfect situation, these descriptors would accurately represent the biological object of the image with negligible deviation from the biological truth (or data). However, in many cases, artifacts might be present in the images so that the representation of the biological object is not accurate anymore. These artifacts might be due to the conditions under which the images were taken or to the object itself. Mature root systems, for instance, are complex branched structures, composed of thousands of overlapping (Figure 1B), and crossing segments (Figure 1C). These features are likely to impede image analysis and create a gap between the descriptors and the data.

Root image descriptors can be separated into two main categories: morphological and geometrical descriptors. Morphological descriptors refer to the shape of the different root segments forming the root system (Table 1). They include, among others, the length and diameter of the different roots. For complex root system images, morphological descriptors are difficult to obtain and are prone to error as mentioned above. Geometrical descriptors give the position of the different root segments in space. They summarize the shape of the root system as a whole. The simplest geometrical descriptors are the width and depth of the root system. Since these descriptors are mostly defined by the external envelope of the root system, crossing and overlapping segments have little impact on their estimation and hence they can be considered as relatively errorless. Geometrical descriptors are expected to be loosely linked to the actual root system topology, since identical shapes could be obtained from

<table><caption>TABLE 1 | Root system parameters used as ground-truth data.</caption><thead><tr><th>Name</th><th>Description</th><th>Unit</th></tr></thead><tbody><tr><td>tot_root_length</td><td>The cumulative length of all roots</td><td>mm</td></tr><tr><td>tot_1_order_length</td><td>The cumulative length of all root axes</td><td>mm</td></tr><tr><td>tot_2+_order_length</td><td>The cumulative length of all lateral roots</td><td>mm</td></tr><tr><td>mean_1_order_length</td><td>The mean first-order roots length</td><td>mm</td></tr><tr><td>mean_2+_order_length</td><td>The mean lateral root length</td><td>mm</td></tr><tr><td>n_1_orders</td><td>The total number of first order roots</td><td>–</td></tr><tr><td>n_2+_orders</td><td>The total number of lateral roots</td><td>–</td></tr><tr><td>mean_2+_order_density</td><td>The mean lateral root density: for each first-order root, the number of lateral roots divided by the axis length (total length).</td><td>mm-1</td></tr><tr><td>mean_1_order_diam</td><td>The mean diameter of the first-order roots</td><td>mm</td></tr><tr><td>mean_2+_order_diam</td><td>The mean diameter of the lateral roots</td><td>mm</td></tr><tr><td>mean_2+_order_angle</td><td>The mean insertion angle of the lateral roots</td><td>°</td></tr></tbody></table>

different root systems (the opposite is true as well). They are usually used in genetic studies, to identify genetic bases of root system shape and soil exploration.

Several automated analysis tools were designed in the last few years to extract both types of descriptors from root images (Armengaud et al., 2009; Galkovskyi et al., 2012; Pierret et al., 2013; Bucksch et al., 2014). However, the validation of such tools is often incomplete and/or error prone. For technical reasons, the validation is usually performed on a small number of ground-truth images of young root systems. In agreement, most analysis tools are specifically designed for this kind of root systems. In the few cases where validation is performed on large and complex root systems, it is usually not on ground-truth images, but in comparison with previously published tools (measurement of X with tool A compared with the same measurement with tool B). This might seem a reasonable approach, regarding the scarcity of ground-truth images of large root systems. However, the inherent limitations of these tools, such as scale or root system type (fibrous- vs. tap-roots) are often not known. Users might not even be aware that such limitations exist and apply the provided algorithm without further validation on their own images. This can lead to unexpected errors in the final measurements.

One strategy to address the lack of in-depth validation of image analysis pipelines would be to use synthetic images generated by structural root models (models designed to recreate the physical structure and shape of root systems). Many structural root models have been developed, either to model specific plant species (Pagès et al., 1989), or to be generic (Pagès et al., 2004, 2013). These models have been repeatedly shown to faithfully represent the root system structure (Pagès and Pellerin, 1996). In addition, they can provide the ground-truth data for each synthetic root system generated, independently of its complexity. However, they have not been used for validation of image analysis tools (Rellán-Álvarez et al., 2015), with one exception performed on young seedling unbranched roots (Benoit et al., 2014).

Here we (i) illustrate the use of a structural root model, Archisimple, to systematically analyse and evaluate an image analysis pipeline and (ii) use the model-generated images to improve the estimation of root traits.
# MATERIALS AND METHODS

## Nomenclature Used in the Paper

**Ground-truth data:** The real (geometrical and morphometrical) properties of the root system as a biological object. They are determined by either manual tracking of roots or by using the output of simulated root systems.

**(Image) Descriptor:** Property of the root image. It does not necessarily have a biological meaning.

**Root axes:** First order roots, directly attached to the shoot.

**Lateral roots:** Second- (or lower) order roots, attached to another root.

## Creation of a Root System Library

We used the model ArchiSimple, which was shown to allow the generation of a large diversity of root systems with a minimal amount of parameters (Pagès et al., 2013). To produce a large library of root systems, we ran the model 10,000 times, each time with a random set of parameters (Figure 2A). For each simulation, the growth and development of the root system were constrained in two dimensions.

The simulations were divided into two main groups: fibrous and tap-rooted. For the fibrous simulations, the model generated a random number of root axes and secondary (radial) growth was disabled. For tap-root simulations, only one root axis was produced and secondary growth was enabled (the extent of which was determined by a random parameter).

The root system created in each simulation was stored in a Root System Markup Language (RSML) file. Each RSML file was then read by the RSML Reader plugin from ImageJ to extract ground-truth data for the library (Lobet et al., 2015). These ground-truth data included geometrical and morphological parameters (Table 1). For each RSML data file, the RSML Reader plugin also created three JPEG images (at a resolution of 300 DPI) for each root system. To simulate one type of image degradation, we added different levels of noise to the images (using the Salt and Pepper Filter in ImageJ) (Figure 2D). For each root system, we computed overlapping index as the number of root segments having an overlap with other root segments over the total number of root segments.

## Root Image Analysis

Each generated image was analyzed using a custom-made ImageJ plugin, Root Image Analysis-J (or RIA-J). For each image, we extracted a set of classical root image descriptors, such as the total root length, the projected area, and the number of visible root tips (Figure 2E). In addition, we included shape descriptors such as the convex-hull area or the exploration ratio (see Supplemental File 1 for details of RIA-J). The list of traits and algorithms used by our pipeline is listed in Table 2. Distribution of the different descriptors is given in the Supplemental Figure 2.

## Data Analysis

Data analysis was performed in R (R Core Team)<fnref n="1" />. Plots were created using ggplot2 (Wickham, 2009) and lattice (Sarkar, 2008).

The Mean Relative Errors (MRE) were estimated using the equation:

$$MRE = \frac{\sum_{1}^{n} \frac{|\bar{y}_{i} - y_{i}|}{y_{i}}}{n}$$

where n is the number of observations, $\bar{y}_i$ is the ground-truth and $y_i$ is the estimated ground-truth.

## Random Forest Framework

A random forest is a state-of-the-art machine learning algorithm typically used for making new predictions (in both classification and regression tasks). Random Forests can perform non-linear predictions and, thus, often outperform linear models. Since its introduction by Breiman (2001), Random Forests have been widely used in many fields from gene regulatory network inference to generic image classification (Huynh-Thu et al., 2013; Marée et al., 2016). Random Forest relies on growing a multitude of decision trees, a prediction algorithm that has shown good performances by itself but, when combined with other decision trees (hence the name forest), returns predictions that are much more robust to outliers and noisy data (see bootstrap aggregating, Breiman, 1996).

In a machine learning setting one is given a set $D = \{(x_1, y_1), (x_2, y_2), ..., (x_n, y_n)\}$,

where $x_i = (x_i^1, x_i^2, ..., x_i^t)$ is an element of a s-dimensional feature space X,

and $y_i = (y_i^1, y_i^2, ..., y_i^t)$ an element of a t-dimensional response space Y.

The learning task is to find a model

$$M: X \rightarrow Y$$

that predicts the data in a good way, where goodness is measured with regard to an error function L.

A decision tree $T_D$ is a machine learning method that, for a dataset D, constructs a binary tree with each node representing a binary question and each leaf a value of the response space. In other words, a prediction can be made from an input value by looking at the set of binary questions that leads to a leaf (e.g., is the first-order root bigger than q1 and if yes is the number of second-order roots smaller than q2 and if no, ...).

Each decision is based upon exactly one feature and is used for deciding which branch of the tree a given input value must take. Hence a decision tree splits successively the set D into smaller subsets and assigns them a value $y_i = T_D(x_i)$ of the response space.

The choice of the feature used for splitting depends on a relevance criterion. In our setting, the default relevance criterion from the randomForest R package (CRAN randomForest, 2015), namely the Gini index, has been used.

A Random Forest

$$F_D = (T_{D,k})_{k \in I} \text{ where } I = \{1, 2, ..., l\} \quad (1)$$

<footnote n="1">¹R Core Team R: A Language and Environment for Statistical Computing.</footnote>
<figure><img src="image_3.png" /><figcaption>TABLE 2 | Root image descriptors extracted by RIA-J.</figcaption></figure>

<table><thead><tr><th>Name</th><th>Description</th><th>Unit</th><th>Reference</th></tr></thead><tbody><tr><td colspan="4"><strong>MORPHOLOGY</strong></td></tr><tr><td>area</td><td>Projected area of the root system</td><td>mm<sup>2</sup></td><td>Galkovskyi et al., 2012</td></tr><tr><td>length</td><td>Length of the skeleton of the root system image</td><td>mm</td><td>Galkovskyi et al., 2012</td></tr><tr><td>tip_count</td><td>Number of end branches in the root system skeleton</td><td>-</td><td></td></tr><tr><td>diam_mean</td><td>Mean diameter of the root object in the image</td><td>mm</td><td></td></tr><tr><td colspan="4"><strong>GEOMETRY</strong></td></tr><tr><td>width</td><td>The maximal width of the root system</td><td>mm</td><td>-</td></tr><tr><td>depth</td><td>The maximal depth of the root system</td><td>mm</td><td>-</td></tr><tr><td>width_depth_ratio</td><td>Ratio between the width and the depth of the root system</td><td>-</td><td>Galkovskyi et al., 2012</td></tr><tr><td>com_x-com_y</td><td>Relative coordinates of the center of mass of the root system</td><td>-</td><td>Galkovskyi et al., 2012</td></tr><tr><td>convexhull</td><td>Area of the smallest convex shape encapsulating the root system</td><td>mm<sup>2</sup></td><td>Galkovskyi et al., 2012</td></tr><tr><td>exploration</td><td>Ratio between the convex hull area and the projected area</td><td>-</td><td>Galkovskyi et al., 2012</td></tr></tbody></table>
consists of *l* decision trees $T_{D,k}$, where several key parameters such as the feature space, are chosen randomly (hence the word Random in the algorithm name). While using a random subspace strongly accelerates the growth of a single tree, it can also decrease its accuracy. However, the use of large number of trees counterbalance advantageously those two effects. The final prediction for each input value $x_i$ corresponds to the majority vote of all the decision trees of the forest $T_{D,k}(x_i)$ in a classification setting while an average of all predicted values is used in a regression task.

## Framework Description

Our method consists of three typical steps:

1. A preprocessing step, where we replace missing values of the training set.

2. A model generation step where, for each response variable, we generate different models according to two Random Forest parameters (number of trees and number of splits).

3. A model selection step, where we choose the best performing pair of parameters of the previous step for each one of the response variable.

## Preprocessing

Missing values in our dataset might arise due to highly noisy images, where the measurement of certain descriptors has been infeasible. To deal with this issue, we first replaced missing values.

This is done using the imputation function of the randomForest R package. It replaces all missing values of a response variable by the median and then a Random Forest is applied on the completed data to predict a more accurate value. We favored 10 trees for computing the new value over the default value of 300 as we found that it offered sufficiently accurate results for our application while being much faster.

## Model Generation

In the model generation step, for each of the response variables, several forests with different number of trees and different number of splits ($t_i, m_j$) are tested. In practice, the training set $D_{train}$ is divided into $m_j$ disjunct subsets $D_{train}^{m_j}$ and on each of those, a Random Forest $F_{D_{train}^{m_j}}$ is trained on a growing number of $t_i$ random trees.

## Model Selection

Given a new data point $x$, each model predicts a response variable $y$ by averaging the predicted values $F_{D_{train}^{m_j}}(x)$, i.e.,

$$ \hat{y} = M_{D_{train}}^{t,m}(x) = \frac{1}{m} \sum_{k=1}^{m} F_{D_{train}^{k}}(x) $$

Then in a final step an estimate of the root-mean-square (RMSE) generalized error on the test set $D_{test}$ is computed, where RSME is defined as

$$ RMSE = \sum_{i=1}^{n} \sqrt{(y_i - \bar{y}_i)^2} $$

for $D_{test} = \{(x_1, y_1), (x_2, y_2), \dots, (x_n, y_n)\}$.

Finally, the model with the parameter pair $(t, m)$ having the minimal error (on the separate test set) is chosen in order to make the predictions.

## Data Availability

All data used in this paper (including the image and RSML libraries) are available at the address http://doi.org/10.5281/zenodo.208214

An archived version of the codes used in this paper is available at the address http://doi.org/10.5281/zenodo.208499

An archived version of the machine learning framework is available at the address https://github.com/FaustFrankenstein/RandomForestFramework/releases/tag/v1.0

# RESULTS AND DISCUSSIONS

## Production of a Large Library of Ground-Truth Root System Images

We combined existing tools into a single pipeline to produce a large library of ground-truth root system images. The pipeline combines a root model (ArchiSimple, Pagès et al., 2013), the Root System Markup Language (RSML) and the RSML Reader plugin from ImageJ (Lobet et al., 2015). In short, ArchiSimple was used to create a large number of root systems, based on random input parameter sets. Each output was stored as an RSML file (Figure 2A), which was then used by the RSML Reader plugin to create a graphical representation of the root system (as a jpeg file) and a ground-truth dataset (Figure 2B). Details about the different steps are presented in the Materials and Methods section.

We used the pipeline to create a library of 10,000 root system images, separated into fibrous (multiple first order roots and no secondary growth) and tap-root systems (one first order root and secondary growth). The ranges of the different ground-truth data are shown in Table 3 and their distribution is shown in the Supplemental Figure 1.

We started by evaluating whether fibrous and tap-root systems should be separated during the analysis. We performed a Principal Component Analysis on the ground-truth dataset to reduce its dimensionality and assess if the type grouping influenced the overall dataset structure (Figure 3A). Fibrous and tap-root systems formed distinct groups (MANOVA $p < 0.001$), with limited overlap. The first principal component, which represented 30.9% of the variation within the dataset, was mostly influenced by the number of first-order axes. The second principal component (19.1% of the variation) was influenced, in part, by the root diameters. These two effects were consistent with the clear root system type grouping, since they expressed the main difference between the two groups of root-system types. Therefore, since the type grouping had such a strong effect on the overall structure, we decided to separate them for the following analyses.

To demonstrate the utility of a synthetic library of ground-truth root systems, we analyzed every image of the library using a custom-built root image analysis tool, RIA-J. We decided to do so since our purpose was to test the usefulness of the synthetic analysis and not to
<table><caption>TABLE 3 | Ranges of the different ground-truth data from the root systems generated using ArchiSimple.</caption><thead><tr><th>Variable</th><th>Minimum value</th><th>Maximum value</th><th>Unit</th></tr></thead><tbody><tr><td colspan="4"><strong>FIBROUS</strong></td></tr><tr><td>tot_root_length</td><td>110</td><td>73971</td><td>mm</td></tr><tr><td>width</td><td>0.76</td><td>302</td><td>mm</td></tr><tr><td>depth</td><td>50</td><td>505</td><td>mm</td></tr><tr><td>n_1_order</td><td>1</td><td>20</td><td>-</td></tr><tr><td>tot_1_order_length</td><td>79</td><td>5409</td><td>mm</td></tr><tr><td>mean_1_order_length</td><td>27</td><td>470</td><td>mm</td></tr><tr><td>mean_1_order_diameter</td><td>0.2</td><td>0.4</td><td>mm</td></tr><tr><td>mean_2+_order_density</td><td>0</td><td>5</td><td>root/mm</td></tr><tr><td>n_2+_orders</td><td>0</td><td>4448</td><td>-</td></tr><tr><td>tot_2+_order_length</td><td>0</td><td>71556</td><td>mm</td></tr><tr><td>mean_2+_order_length</td><td>0</td><td>50</td><td>mm</td></tr><tr><td>mean_2+_order_diameter</td><td>0</td><td>0.3</td><td>mm</td></tr><tr><td>mean_2+_order_angle</td><td>0</td><td>88</td><td>°</td></tr><tr><td colspan="4"><strong>TAP ROOTED</strong></td></tr><tr><td>tot_root_length</td><td>78</td><td>41870</td><td>mm</td></tr><tr><td>width</td><td>0.1</td><td>173</td><td>mm</td></tr><tr><td>depth</td><td>55</td><td>505</td><td>mm</td></tr><tr><td>n_1_order</td><td>1</td><td>1</td><td>-</td></tr><tr><td>tot_1_order_length</td><td>59</td><td>509</td><td>mm</td></tr><tr><td>mean_1_order_length</td><td>59</td><td>509</td><td>mm</td></tr><tr><td>mean_1_order_diameter</td><td>0.2</td><td>16</td><td>mm</td></tr><tr><td>mean_2+_order_density</td><td>0</td><td>2.3</td><td>root/mm</td></tr><tr><td>n_2+_orders</td><td>0</td><td>3353</td><td>-</td></tr><tr><td>tot_2+_order_length</td><td>0</td><td>40225</td><td>mm</td></tr><tr><td>mean_2+_order_length</td><td>0</td><td>51</td><td>mm</td></tr><tr><td>mean_2+_order_diameter</td><td>0</td><td>2.2</td><td>mm</td></tr><tr><td>mean_2+_order_angle</td><td>0</td><td>97</td><td>°</td></tr></tbody></table>

assess the accuracy of existing tools. Nonetheless, RIA-J was designed using known and published algorithms, often used in root system quantification. A detailed description of RIA-J can be found in the Materials and Methods section and **Supplemental File 1**.

We extracted 10 descriptors from each root system image (Table 2) and compared them with the ground-truth data. For each pair of descriptor-data, we performed a linear regression and computed its r-squared value. Different types of information are highlighted in Figure 4. First, using a ground-truth image library allows for a quick and systematic analysis of all the descriptors extracted by the image analysis pipeline. Second, it allows researchers to identify which traits can be accurately evaluated (or not) and by which descriptors. Third, for some ground-truth data, such as the mean length of second order roots or the number of first order roots, it shows that none of the classical descriptors gave a good estimation (Figure 4, highlighted with arrows). Finally, the figure highlights that some correlations were different for fibrous- and tap-root systems. As an example, the correlation found between the *mean_2+_order_diameter* and *diam_mean* estimators was better for fibrous roots than within the tap-root dataset. Consequently, validation of the

<figure><img src="image_4.png" /></figure>

different image analysis algorithms should be performed, at least, for each group. An algorithm giving good results for a fibrous root system might fail when applied to tap-rooted ones.
<figure><img src="image_5.png" /><figcaption>**FIGURE 4** | Heatmap of the r-squared values between the different image descriptors and the ground-truth values, for the images without any noise. Black represents an r-squared value of 1; white represents a value of 0. Upper panel: fibrous root dataset. Lower panel: tap-root dataset. Arrows highlight the ground-truth data that cannot be accurately described with the different descriptors. The arrows were doubled when it was the case for both fibrous and tap-rooted root systems.</figcaption></figure>

## Errors from Image Descriptors Are Likely to Be Non-linear across Root System Sizes and Image Qualities

In addition to being related to the species of study, estimation errors are likely to increase with the root system size. As the root system grows and develops, the number of crossing and overlapping segments increases (Figure 5A), making the subsequent image analysis potentially more difficult and prone to error. However, a systematic analysis of such error is seldom performed.

Estimation errors are also likely to increase as the image quality decreases. Here we artificially added one type of noise (random “salt and pepper” particles) to the images, with two intensity levels. It should be noted that virtually any type of image degradation could be added to the original images using custom image filters (e.g., using ImageJ). Different types of degradation are expected to generate different levels of estimation errors.

**Figure 5** shows the relationship between the ground-truth and descriptor values for three parameters: the total root length (Figure 5B), the number of roots (Figure 5C), and the root system depth (Figure 5D). For each of these variables, we quantified the Mean Relative Error (see Materials and Methods for details) as a function of the overlap index. This was done for three levels of noise added to the images (“null,” “medium,” and “high”). We can observe that for the estimation of both the total root length and the number of lateral roots, the Mean Relative Error increased with the size of the root system (Figures 5B–C). As stated above, such increase of the error was somehow expected with increasing complexity. Moreover, depending on the metric of interest, such as the number of root tips, low image quality can result in high level of error. For other traits, such as the root system depth, no errors were expected (depth is supposedly an error-less variable) and the Mean Relative Error was close to 0 whatever the size of the root system and image quality.

The results presented here are tightly dependent on the specific algorithms used for image analysis and hence might be different for other published tools. However, they are a call for caution when analyzing root images: unexpected errors in ground-truth estimation can arise. Our image library can be used to better identify the errors generated by other analysis tools, current or future.

## Roadmap for Root Image Analysis Tools Calibration

To improve the calibration and validation of future root image analysis tools, we propose the following procedure:

1. Develop the new root image analysis pipeline;

2. Use it to analyse the images from the synthetic root library described here;

3. Compare the results from the new analysis with the corresponding ground-truth;

4. Identify, and clearly state, the type of root systems for which the pipeline works accurately;

5. When releasing the new pipeline, inform the users about the possible errors identified.

## Using the Synthetic Library to Train Machine Learning Algorithms

The main advantage of creating a synthetic library is to generate paired datasets of image descriptors and their corresponding ground-truth values. Having information on both can, in theory, be used to either calibrate the image analysis pipeline or to identify the best descriptors for the ground-truth traits of interest. Here, we explored the second approach and used a Random Forest algorithm to find which combination of descriptors would best describe each ground-truth data (see Material and Methods for details). In short, we randomly divided the whole dataset into training (3/4) and testing subsets (1/4). The training set was used to create a Random Forest model for each ground-truth data, which was then applied to the test set. The accuracy of these new predictions was then compared to the accuracy of the direct method (single descriptors) (Figure 2C).
<figure><img src="image_6.png" /><figcaption>**FIGURE 5** | Error estimation for three ground-truth parameters. (A) Evolution of the overlap index (proportion of root overlapping) with the root system size. (B–D) Left panel shows the relationship between the descriptors and the corresponding ground-truth variables. Right panels show the evolution of the Mean Relative Error (MRE) as a function of the overlap index. For the MRE calculations, the continuous variables were discretized in groups. (B) Number of lateral roots. (C) Total root length. (D) Root system depth. In the left panels, the gray line indicates the diagonal (1:1 relation).</figcaption></figure>
**Figure 6** shows the comparison of the accuracy (both the r-squared values from linear regressions and the Mean Relative Error, MRE) of both methods for each ground-truth data. We can clearly see that the Random Forest approach performed always better (sometimes substantially) than the direct approach, even for images with high level of noise. In addition, for most traits, the r-squared and MRE values were above 0.9 and below 0.1 respectively, which is very good, especially for such a wide range of images. In addition, the Random Forest approach allowed the correct estimation of traits that were difficult to estimate with the direct approach (such as the number of first-order axes or the mean second-order root density).

**Figure 7** shows the detailed comparison of both methods for the estimation of the total root length. Again, a clear improvement was visible with the Random Forest method, leading to small errors, even with large root systems and noisy images.

Here we presented how machine learning algorithms (Random Forest), could be used in combination with a

synthetic image library to improve the estimation of root system traits. Although both the training and test datasets used were made of synthetic images, we believe this approach presents an interesting perspective for the analysis of experimental images.

Indeed, a root architectural model can be used to build a custom library of synthetic images from a set of parameters evaluated on a small number of plants from the experimental dataset. Such library could then be used to train the machine learning model which, in turn, will enable the automatic evaluation of root traits from the remaining experimental images. Alternatively, the algorithm could be directly trained on a subset of experimental data obtained by manual or semi-automatic analyses to be then automatically applied to the rest of the dataset. One must keep in mind that the output of the machine learning strongly depends upon the quality of the dataset used for its training and hence must be analyzed carefully.

<figure><img src="image_7.png" /><figcaption>**FIGURE 6** | Comparison between the direct trait and the Random Forest approach, for the different root system types and the different levels of noise.
For each metric, we computed both the r-squared value from the linear regression between the estimation and the ground-truth (left panels), as well as the Mean Relative Error (right panel). The gray points represent the values obtained with the direct estimation (best descriptor, no noise). Color points represent the values obtained with the Random Forest approach, for different levels of noise. The dotted lines show the 0.9 (r-squared) and 0.1(MRE) thresholds.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**FIGURE 7** | Comparison between the direct trait estimation and the Random Forest approach, for the different root system types and the different levels of noise. (A) Comparison, for the total root length, of the accuracy of both approaches. The dotted line represents the diagonal. The plain line represents the linear regression. (B) Same, for the number of roots.</figcaption></figure>

## CONCLUSIONS

The automated analysis of root system images is routinely performed in many research projects. Here we used a library of 10,000 synthetic images to estimate the accuracy and usefulness

of different image descriptors extracted with a homemade root image analysis pipeline. Our study highlighted some limitations and biases of the image analysis process.

We found that the type of root system (fibrous vs. tap-rooted), its size and complexity, as well as the quality of the images had a strong influence on the accuracy of some commonly used image descriptors and their meaning and relevance for ground-truth extraction. So far, a large proportion of the root research has been focused on seedlings with small root systems and has *de facto* avoided such errors.

However, as the research questions are likely to focus more on mature root systems in the future, these limitations will become critical. We showed that synthetic datasets can be used for calibration or modeling (machine learning) steps that allow ground-truth extraction from comparable images. We then hope that our library will be helpful for the root research community to evaluate and improve other image analysis pipelines.

## AUTHOR CONTRIBUTIONS

GL, LP, PT, IK, MN, and CP designed the study. IK developed the image analysis pipeline RIA-J. MN and PM developed the Random Forest framework. GL generated the image library and did the data analysis. LP developed the Archisimple model. All authors have participated in the writing of the manuscript.

## FUNDING

This research was funded by the Interuniversity Attraction Poles Programme initiated by the Belgian Science Policy Office, P7/29. GL and MN are grateful to the F.R.S.-FNRS for a postdoctoral research grant (1.B.237.15F) and doctoral grant (1.A.320.16F), respectively.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00447/full#supplementary-material

<caption>**Supplemental Figure 1** | Distribution of the properties of the modeled root images.</caption>

<caption>**Supplemental Figure 2** | Distribution of the descriptors of the modeled root images.</caption>

<caption>**Supplemental File 1** | Definitions of the different descriptors extracted by RIA-J.</caption>

## REFERENCES

Armengaud, P., Zambaux, P., Hills, A., Sulpice, R., Pattison, R. J., Blatt, M. R., et al. (2009). EZ-Rhizo: integrated software for the fast and accurate measurement of root system architecture. *Plant J.* 57, 945–956. doi: 10.1111/j.1365-313X.2008.03739.x

Benoit, L., Rousseau, D., Belin, É., Demilly, D., and Chapeau-Blondeau, F. (2014). Simulation of image acquisition in machine vision dedicated to seedling elongation to validate image processing root segmentation

algorithms. *Comput. Electron. Agric.* 104, 84–92. doi: 10.1016/j.compag.2014.04.001

Breiman, L. (1996). Bagging predictors. *Mach. Learn.* 24, 123–140. doi: 10.1007/BF00058655

Breiman, L. (2001). Random Forests. *Mach. Learn.* 45, 5–32. doi: 10.1023/A:1010933404324

Bucksch, A., Burridge, J., York, L. M., Das, A., Nord, E., Weitz, J. S., et al. (2014). Image-based high-throughput field phenotyping of crop roots. *Plant Physiol.* 166, 470–486. doi: 10.1104/pp.114.243519
CRAN randomForest (2015). Available online at: https://cran.r-project.org/web/packages/randomForest/randomForest.pdf

Galkovskyi, T., Mileyko, Y., Bucksch, A., Moore, B., Symonova, O., Price, C. A., et al. (2012). *GiA Roots: software for the high throughput analysis of plant root system architecture*. *BMC Plant Biol.* 12:116. doi: 10.1186/1471-2229-12-116

Huynh-Thu, V. A., Wehenkel, L., and Geurts, P. (2013). *Gene Regulatory Network Inference from Systems Genetics Data Using Tree-Based Methods*. New York, NY: Springer.

Koevoets, I. T., Venema, J. H., Elzenga, J. T., and Testerink, C. (2016). Roots withstanding their environment: exploiting root system architecture responses to abiotic stress to improve crop tolerance. *Front. Plant Sci.* 7:1335. doi: 10.3389/fpls.2016.01335

Lobet, G., Pound, M. P., Diener, J., Pradal, C., Draye, X., Godin, C., et al. (2015). Root system markup language: toward a unified root architecture description language. *Plant Physiol.* 167, 617-627. doi: 10.1104/pp.114.253625

Marée, R., Geurts, P., and Wehenkel, L. (2016). Towards generic image classification using tree-based learning: an extensive empirical study. *Pattern Recognit. Lett.* 74, 17-23. doi: 10.1016/j.patrec.2016.01.006

Pagès, L., Bécel, C., Boukcim, H., Moreau, D., Nguyen, C., and Voisin, A.-S. (2013). Calibration and evaluation of ArchiSimple, a simple model of root system architecture. *Ecol. Modell.* 290, 76-84. doi: 10.1016/j.ecolmodel.2013.11.014

Pagès, L., Jordan, M. O., and Picard, D. (1989). A simulation model of the three-dimensional architecture of the maize root system. *Plant Soil* 119, 147-154. doi: 10.1007/BF02370279

Pagès, L., and Pellerin, S. (1996). Study of differences between vertical root maps observed in a maize crop and simulated maps obtained using a model for the

three-dimensional architecture of the root system. *Plant Soil* 182, 329-337. doi: 10.1007/BF00029063

Pagès, L., Vercambre, G., Drouet, J.-L., Lecompte, F., Collet, C., and LeBot, J. (2004). RootTyp: a generic model to depict and analyze the root system architecture. *Plant Soil* 258, 103-119. doi: 10.1023/B:PLSO.0000016540.47134.03

Pierret, A., Gonkhamdee, S., Jourdan, C., and Maeght, J.-L. (2013). IJ-Rhizo: an open-source software to measure scanned images of root samples. *Plant Soil* 373, 531-539. doi: 10.1007/s11104-013-1795-9

Rellán-Álvarez, R., Lobet, G., Lindner, H., Pradier, P.-L., Sebastian, J., Yee, M.-C., et al. (2015). GLO-Roots: an imaging platform enabling multidimensional characterization of soil-grown root systems. *Elife* 4:e07597. doi: 10.7554/eLife.07597

Sarkar, D. (2008). *Lattice: Multivariate Data Visualization with R*. New York, NY: Springer.

Wickham, H. (2009). *ggplot2*. New York, NY: Springer New York.

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2017 Lobet, Koevoets, Noll, Meyer, Tocquin, Pagès and Périlleux. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.