Skip to main content

Root phenotypic detection of different vigorous maize seeds based on Progressive Corrosion Joining algorithm of image



The root phenotypes of different vigorous maize seeds vary a lot. Imaging roots of growing maize is a non-invasive, affordable and high throughput approach. However, it’s difficult to get integral root images because of the block of the soil. The paper proposed an algorithm to repair incomplete root images for maize root fast non-invasive phenotyping detection.


A two-layer transparent stress growth device with two concentric cylinders was developed as mesocosms and the maize seeds were planted in the annulus of it. The maize roots grow in soil against two acrylic plastic surfaces due to the press of the small growing area to acquire more root details during roots visualization and imaging. Even though, parts of the roots are occluded which means that it’s tough to extract the information of root general physical construction. For recovering gaps from disconnected root segments, Progressive Corrosion Joining (PCJ) algorithm was proposed based on the physiological characteristics of hydrotropism, geostrophic and continuity with three steps which are root image thinning, progressive corrosion and joining processing respectively. The experiments indicate that maize phenotyping parameters are negative correlation with seed aging days. And specifically, Root Number (RTN), Root Length (RTL), Root Width (RTW) and Root Extension Length (REL) of unaged and 14-day-aged maize seeds are decreased from 15.40, 82.40 mm, 1.53 mm and 82.20 mm to 4.58, 38.6 mm, 1.35 mm and 55.20 mm, and the growing speed of them are changed from 1.68 per day, 8.80 mm/d, 0.06 mm/d, 9.0 mm/d to 0.70 per day, 4.3 mm/d, 0.05 mm/d and 5.70 mm/d respectively. Whereas Root Extension Angle (REA) is basically irrelevant with the level of maize seed aging.


The developed double-layer Annular Root Phenotyping Container (ARPC) can satisfy the general physical construction of maize as well as push each root growing along the inner wall of the container which help to acquire more root information. The presented novel PCJ algorithm can recover the missing parts, even for big gaps, of maize roots effectively according to root morphological properties. The experiments show that the proposed method can be applied to evaluate the vigor of maize seeds which has vast application prospect in high throughput root phenotyping area.


Plant root is an important organ for water, nutrients and minerals uptake. Root morphology research is an important content of plant nutrition, physiology, breeding and ecology [1, 2]. As the largest production crop in the world, maize can be used as both food and fodder whose root, as the plant holder, greatly affects the ability of press resistance such as water and wind which significantly results in final yield reduction. Therefore, root phenotyping parameters of maize seeding such as root spatial distribution, root length, root width, root extension characteristics, root hair density and bending angle are applied to evaluate maize seed vigor.

For root character observation clearly and directly, germination paper method was proposed to seed screening such as soybean [3], corn [4], wheat [5], rapeseed [6] and pearl millet [7]. Computer vision technology was introduced to speed up germination paper detection but the root can only grow on the surface of paper [8]. To extend plant root grow in 3D space from 2D surface, hydroponics [9] and atomization cultivation [10] patterns were used but root morphology is affected by gravity. So, the gel transparent matrix which can keep the plant roots in a consistent spatial structure has been developed, whereas the chemical properties of the matrix are different from that of soil which may affect root development and configuration [11]. The research mentioned above can overcome the problem of soil interference, control the homogeneity of the culture medium and minimize the inherent variability of the observed root traits as well as enable clear imaging. However, plant root system growing environment in the field is essentially different from that of air, water and matrix.

Therefore, it becomes a key problem to obtain as much root information as possible on the premise of ensuring the root growth in the soil. The most direct and effective traditional approach such as soil core, soil column and root canal [12, 13] were applied to study a root system by digging out the root and observe it by manual inspecting, machine vision [14], spectral imaging [15] and laser scanning imaging [16] technologies after removing the soil. However, it will cause irreversible damage and configurational change to the root system. To achieve non-destructive in situ plant root phenotypic detection, X-ray computed tomography (CT) [17, 18] and nuclear magnetic resonance (NMR) imaging technology [19] used in medicine science popularly were introduced to obtain clear stereoscopic images of roots, but they hurt the roots to some extent after repeated measurements and the instruments are also very expensive.

Minirhizotron [20, 21] and its improvement methods [22], also as a non-destructive way, were provided to monitor plant roots over a long period of time in the field for forests, orchards and agro-ecosystems research [23,24,25,26] but are labor intensive and time consuming, even more serious is that the methods can only sample local pictures of roots and lead to excessive root reproduction because the set can disturb the balance of roots and soil.

For obtaining root general and detailed information in the soil rapidly, non-destructively and simultaneously, root phenotype pipeline system [27] was developed to image plant roots in cylindrical or flat-plate type transparent culture containers which is popularly used in breeding industry nowadays, whereas parts of root are occluded by soil which results in incomplete root image even under perfect root segmentation. So, root gap correcting algorithm was proposed for the extraction of finely grained root system architecture traits using deep learning CNN and obtained good effects but the training set and testing set are idealized manual drawn pictures which need a lot of manual work [28,29,30].

In this paper, we developed a double-layer Annular Root Phenotyping Container (ARPC) to acquire more root information considering the spatial configuration of maize root. Five rules were given for maize root gaps recovering considering the physiological characteristics of hydrotropism, geostrophic and continuity to propose a novel maize root gaps recovering algorithm, Progressive Corrosion Joining (PCJ) algorithm. Finally, the proposed algorithm was applied to improve maize root images which can be used to evaluate the vigor of different aging-day seeds effectively. The proposed method, compared with other existed root phenotyping approaches, can be used for high-throughput, low-cost, non-destructive phenotypic detection of maize roots, and can obtain more phenotypic information of maize roots. So, the proposed ARPC combing with PCJ algorithm is a good method which has vast application prospect in high throughput root phenotyping area (Tables 1, 2).

Table 1 The captures of the cylinder
Table 2 The captures of the ring

Materials and methods

Plant culture device and seed growth conditions

Because of the root growing in straight lines and its tendency such as hydrotropism, geostrophic, a double-layer ARPC device (Fig. 1) was developed to get more root information. The APRC device consists of two concentric hollow acrylic cylinders of different diameters with the same height and a ring with 8 leaking holes. The hollow annular between two acrylic layers are filled with homogeneous cultivated soil (organic matter 478 g/1000 g, N + P+K: 6.75 mg/1000 g, probiotics: 90 μg/1000 g, trace element 1.7 μg/1000 g, humic acid: 130 μg/1000 g, pH 6.5–6.8). Heinuo-1 (Black Corn, Shangdu Seed Co., Ltd, Henan, China) maize seeds were sown in the annular in the greenhouse in Nanjing Agricultural University, College of Engineering (32° 18′ N, 118° 46′ E). The parameters of the device are shown in Fig. 1.

Fig. 1
figure 1

The double-layer Annular Root Phenotyping Container. a is the vertical view of ARPC. b is the front view of ARPC. c is the simulation view, and the green part of the root can be seen by the naked eye

Image acquisition device

The major factor influencing image processing quality is the reflection from the device. To obtain high quality pictures, the following devices are used.

  1. 1.

    Camera: Canon D600 digital SLR camera (18–135 mm medium focal length lens, Canon). The camera was set to macro mode. The shutter speed, aperture and sensor sensitivity were set as 1/3 s, 5 mm and 800 ISO respectively.

  2. 2.

    Polarizer: A set of polarizers were placed in order along the axis of light to eliminate reflection from the plant container. The angle of polarization between the polarizer and analyzer was adjusted manually to optimum condition.

  3. 3.

    The centers of camera lens, polarizer, analyzer and root were set along the same light axis. The lens of camera focused on the root center to acquire clear root images (Fig. 2).

    Fig. 2
    figure 2

    Shooting scene pictures

Experimental procedure

  1. 1.

    Seed aging: The screened seeds were divided into 8 parts after dormancy treatment and were artificially aged for 0, 2, 4, 6, 8, 10, 12 and 14 days respectively. During aging time,seeds were placed in glass trays and enclosed in intelligent artificial climate chamber (RXZ type, Ningbo Jiangnan Instrument Factory) under the condition of 45 °C and 90% humidity.

  2. 2.

    Seed disinfection: Black corn Heinuo-1 seeds with the germination not less than 85% were used for root phenotyping experiment. The seeds were immersed in 70% (v/v) ethanol and 5% sodium hypochlorite solution for 30 s and 10 min respectively before transferred into ARPC device.

  3. 3.

    Seeding culture: Two seeds were sowed in each ARPC device and the seedlings were poured once every 5 days with 5 ml volume at 9 am after seeds sowing. The two sides of ARPC device were covered with black plastic film to simulate root growing luminous environment in the field (Fig. 3).

    Fig. 3
    figure 3

    The picture of shaded device

  4. 4.

    Maize root images acquisition: Six seeds of every aging day were used for experiment, as a total 48 seedings were used for root imaging. Root images were collected indoor twice every day at 9 am and 4 pm for monitoring the development of seeding roots dynamically, and 4 pictures with different perspectives (covering 360°) were acquired every time. Each culture container was fixed on the same marked position and rotated manually while imaging to get the sequence of images cover 360° views.

Root repair


  1. a.

    Image stitching: The paper judges the continuity of maize root by detecting whether there are more than 4 consecutive white pixels at the image edge. And the feature points between the adjacent images are extracted for stitching a complete root according to the L-ORB algorithm if the edge of root image is incomplete [31].

  2. b.

    Image calibration: Considering that the images taken by the camera are distorted, a piece of checkerboard paper was pasted to the surface of the ARPC for camera calibration with reference to the literature [32] (Fig. 4).

    Fig. 4
    figure 4

    The device pasted with the checkered paper for correcting

  3. c.

    The transformational approaches such as gray, binary, median filtering, gaussian filtering and flooding algorithms were used successively as preprocessing to obtain clear images from image noise made by light reflection and soil color difference.

Feature point extraction

Two rules were proposed to determine root feature points as below before image processing algorithm.

  1. (a)

    Endpoint decision rule: One endpoint only connects and belongs to one root. Define the image area as Ω, P is one of the point on it. The minimized square D centered at point P with 9 pixels is selected to judge whether point P connects with no more than 2 connected points. The mathematical expression of the rule is as follows.

    $$ e\_p = \left\{ {\begin{aligned} & {Ture,\;\;n = 1} \\ &{Ture,\;\;(n = 2)\;and\;(Q\;{\text{are}}\;connected)} \\ &{False,\;\;else} \\ \end{aligned} } \right. $$

    m is 8 points in square D except point P, point Q \( \in \) m, n is the number of points in m. e_p is the end point.

  2. (b)

    Bifurcation point rule: The branching point is connected to multiple roots. Define the image area as Ω, P is one of the point on it. The minimized square D centered at point P with 9 pixels is selected to judge whether point P connects with not less than 2 disconnected points. The mathematical expression of the rule is as follows.

    $$ b\_p = \left\{ \begin{aligned} & Ture,\;\;(n > 1)\;and\;(Q\;{\text{are}}\;connected) \hfill \\ & False,\;\;else \hfill \\ \end{aligned} \right. $$

    m is 8 points in square D except point P, point Q \( \in \) m, n is the number of points in m. b_p is the bifurcation point.

Based on the above two rules, the Progressive Corrosion Joining algorithm is presented as follow steps:

  • Step 1: Image thinning

    Zhang-Suen thinning algorithm was used to refine the root trunk [33]. Firstly, all the foreground pixels are traversed and removed if they meet the following conditions (Fig. 5).

    Fig. 5
    figure 5

    Foreground pixel illustration

    $$ \left\{ \begin{array}{ll} 2 \le N(P1) \le 6 \hfill \\ S(P1) = 1 \hfill & \\ P2*P4*P6 = 0 \hfill & \\ P4*P6*P8 = 0 \hfill \end{array} \right. $$

    where N(P1) represents the number of foreground pixel points among the 8 pixels adjacent to P1; S(P1) represents the cumulative numbers of 0–1 model, the adjacent two pixels are 0 and 1, counting along P2-P9-P2 line, where 0 represents the background and 1 indicates the prospect.

    Then, P1 is also traversed and removed if it satisfies the conditions below.

    $$ \left\{ \begin{array}{ll} 2 \le N(P1) \le 6 \hfill \\ S(P1) = 1 \hfill \\ P2*P4*P8 = 0 \hfill \\ P2*P6*P8 = 0 \hfill \\ \end{array} \right. $$

    The refined root will be obtained after the above two steps.

  • Step 2: Endpoint extraction

    The endpoint set M is created by using the endpoint decision rule. The coordinates of all the elements in M are stored in the queue “a”, the elements in M are numbered and the serial numbers are stored in queue “b”.

    To obtain the endpoint point set M of the refined image, and numbering the element M[i] in M, storing the element M[i] in the queue “a”, and storing the sequence number in the queue “b”. There are only two endpoints for each line, and after subsequent processing, the point set numbers on each line will be the same (Fig. 6).

    Fig. 6
    figure 6

    A queue for dynamically storing endpoints information

  • Step 3: Image mapping

    Create an empty image “c” with the same size of thinning image, then save the elements of “b” in the corresponding coordinates of “c”, that is, the coordinates of the elements in “c” are equal to the values of “a”.

    $$ c(a[i]) = b[index(a[i])] $$

    where index is the coordinate of the index value of a queue.

  • Step 4: Progressive corrosion

The steps of the proposed progressive corrosion algorithm are as follow and loop until “a” is empty.

  1. a.

    A point P pops off from queue “a”, then the value of the point in the thinning image which coordinate equal to the value of P is set to 0.

  2. b.

    If one of the surrounding 8 pixels, defined as point q, of P in “c” has been assigned a serial number, then assign the value of P to q.

  3. c.

    If P is an endpoint, the surrounding point connected to P is named as P1. And the value of the point with the same coordinate with P1 in “c” is assigned to the value of P.

  4. d.

    If P is a bifurcation point, the coordinates of the surrounding points of P are pushed into “a”. Then the surrounding points are numbered in “b” and the corresponding points in “c” are assigned the corresponding value.

The flowchart of the PCJ algorithm can be seen at Fig. 7.

Fig. 7
figure 7

Progressive Corrosion algorithm Flowchart

Root connection

After the operation of above steps, each continuous segment of a root has the same serial number. Three rules were proposed for root gap connection.

Internal continuity rule

The pixels of each root segment in the thinning image are step-forced but a root should be smooth actually. So, the least squares method is used to fit the pixels as shown in the formula below.

$$ f(x) = a_{0} + a_{1} x + a_{2} x + \cdots + a_{k} x_{k} $$

where x and f(x) are the abscissa and ordinate respectively.

The loss value is expressed as below.

$$ l = \sum\limits_{i = 1}^{n} {[y_{i} - (a_{0} + a_{1} x + \cdots + a_{k} x_{k} )]^{2} } $$

l is loss.

The partial derivatives of ai (i = 1…k) are absolved and let be 0, we can get the Vandermonde matrix,

$$ \left[ {\begin{array}{*{20}c} n & {\sum\limits_{{i = 1}}^{n} {x_{i} } } & \cdots & {\sum\limits_{{i = 1}}^{n} {x_{i}^{k} } } \\ {\sum\limits_{{i = 1}}^{n} {x_{i} } } & {\sum\limits_{{i = 1}}^{n} {x_{i}^{2} } } & \cdots & {\sum\limits_{{i = 1}}^{n} {x_{i}^{{k + 1}} } } \\ \vdots & \vdots & \ddots & \vdots \\ {\sum\limits_{{i = 1}}^{n} {x_{i}^{k} } } & {\sum\limits_{{i = 1}}^{n} {x_{i}^{{k + 1}} } } & \cdots & {\sum\limits_{{i = 1}}^{n} {x_{i}^{{2k}} } } \\ \end{array} } \right] = {\text{ }}\left[ {\begin{array}{*{20}c} {a_{0} } \\ {a_{1} } \\ \vdots \\ {a_{k} } \\ \end{array} } \right] = \left[ {\begin{array}{*{20}c} {\sum\limits_{{i = 1}}^{n} {y_{i} } } \\ {\sum\limits_{{i = 1}}^{n} {x_{i} y_{i} } } \\ \vdots \\ {\sum\limits_{{i = 1}}^{n} {x_{k} y_{k} } } \\ \end{array} } \right] $$

The matrix can be simplified as,

$$ \left[ {\begin{array}{*{20}c} 1 & {x_{1} } & \cdots & {x_{1}^{k} } \\ 1 & {x_{2} } & \cdots & {x_{2}^{k} } \\ \vdots & \vdots & \ddots & \vdots \\ 1 & {x_{n} } & \cdots & {x_{n}^{k} } \\ \end{array} } \right]\left[ \begin{aligned} a_{0} \hfill \\ a_{1} \hfill \\ \vdots \hfill \\ a_{k} \hfill \\ \end{aligned} \right] = \left[ \begin{aligned} y_{1} \hfill \\ y_{2} \hfill \\ \vdots \hfill \\ y_{n} \hfill \\ \end{aligned} \right] $$


$$ x = \left[ {\begin{array}{*{20}c} 1 & {x_{1} } & \cdots & {x_{1}^{k} } \\ 1 & {x_{2} } & \cdots & {x_{2}^{k} } \\ \vdots & \vdots & \ddots & \vdots \\ 1 & {x_{n} } & \cdots & {x_{n}^{k} } \\ \end{array} } \right]\,\,A = \left[ \begin{aligned} a_{0} \hfill \\ a_{1} \hfill \\ \vdots \hfill \\ a_{k} \hfill \\ \end{aligned} \right]\,\,Y = \left[ \begin{aligned} y_{1} \hfill \\ y_{2} \hfill \\ \vdots \hfill \\ y_{n} \hfill \\ \end{aligned} \right] $$

So, the coefficient matrix A is obtained as A = (XTX)1 XTY, and at the same time, we can get the fitting curve.

Segmental continuity rule

The plant root grows continuously which means that the derivative difference of a root changes slowly. Therefore, the derivative difference values of adjacent root segments are very close. And the endpoints derivatives of the adjacent root segments of the same root vary within a narrow range.

Proximity rule

On the basis of the above rules, the endpoints of discrete root segments connect to that of nearest segments.

The detailed steps for root segments connection is as follow.

  • Step 1: Root segment characteristics extraction

    The derivative differences of two endpoints and the straight line slope of a root segment calculated according to the method of least squares are obtained as root segment characteristics for root segments matching.

  • Step 2: Root segments matching

    The root segments match well as the same root if they meet the following conditions.

    1. a.

      The difference of the straight line slopes of the two segments is less than a certain value d1 and at the same time is more than a certain value d2.

    2. b.

      The derivative difference values of adjacent endpoints of the two root segments for intended connection is not greater than a set value d3.

    3. c.

      The loss value of the two segments calculated according to formula (7) is not more than a special value d4.

    4. d.

      Based on the root growing trait of divergence, a root segment only connects to one near root segment but can be concatenated by more segments.

    The root segments matching method need 4 parameters which were adjusted manually for obtaining good results which optimization values are d1 = 4, d2 = 10, d3 = 50, d4 = 0.2.

  • Step 3: Root segments connecting

    The adjacent endpoints of two root segments are linked together by using the quantic polynomial interpolation and minimum mean square error fitting method if they meet the above conditions.

    The tiny segments (not more than 10 pixels) are filled with black. These segments can be classified to two types: one is tiny root fragments, and the other one is noise caused by reflection light and impurities in the soil. The tiny fragments will be covered by the connection line created using the proposed PCJ algorithm and the noise will be removed as well. In addition, the reduction of plenty of tiny root segments can also speed up the computer processing.

    The whole flowchart of PCJ algorithm is as shown in Fig. 8.

    Fig. 8
    figure 8

    PCJ algorithm flowchart

Root phenotypic parameters for analysis

Five root physiological parameters such as Root Number (RTN, obtained by calculating the number of end points in the repaired image), Root Length (RTL, acquired by calculating the total number of pixel points of each fitting line in the repaired image according to the L-ORB algorithm), Root Width (RTW, got by calculating the mean number of pixels perpendicular to the extension direction on each root in the binary image based on the L-ORB algorithm), Root Extension Length (REL, obtained by calculating the horizontal distance between the leftmost white pixel(s) and the rightmost pixel(s)) and Root Extension Angle (REA, calculated by measuring the angle between the two fitting lines of roots) are selected for root phenotyping analysis because they are important indexes of maize growing condition which can ensure the efficient absorption of moisture and nutrient uptake [34]. Before the measurement of the above five parameters, the relationship between the physical size and root image of a root is calibrated. Finally, the root phenotypic parameters are used to evaluate maize seed vigor to verify the validity of the method.


Image processing results of PCJ algorithm

In our experiments, as shown in Fig. 9a, there are some obvious gaps in the original maize root image after denoising processing, and the length and width of different root segments vary a lot. After binary processing, a clear boundary between root system and the soil is obtained and the width of different root segments can be calculated as shown in Fig. 9b.

Fig. 9
figure 9

The result of the algorithm processing the picture of the maize root. a The image after denoising; b the image after binarization; c the image after thinning; d the end point of extracted maize root. Every endpoint is zoomed up by 5 times for the sake of observation; e The image processed by PCJ algorithm, and the points belonging to the same root are in the same color; f The completed image

Then the skeletal image of the root is got according to thinning algorithm, at this point the width of all root segments is one pixel in Fig. 9c. In Fig. 9d, all the endpoints of root segments are extracted as the beginning points of corrosion step. After progressive corrosion treatment, the different segments of the same root are marked as the same ID with the same color in Fig. 9e. Finally, the disconnected root segments are filled with green curves along the fitting curves for root recovering in Fig. 9f.

Root phenotypic parameters of different vigor maize seeds

Static phenotypic parameters of maize seed

Root growing condition of a seed at the early stage is a highly effective way to evaluate seed vigor. In this paper, five root physiological parameters such as Root Number (RTN), Root Length (RTL), Root Width (RTW), Root Extension Length (REL) and Root Extension Angle (REA) are used as root phenotypic parameters for seed vigor indexes. As shown in Tables 3, 4 and 5, the experimental results of the root phenotypic parameters above calculated based on the refined root images using the presented PCJ algorithm show that there is remarkable correlation between maize vigor and root phenotypic parameters mentioned above. For the performance evaluation of the proposed PCJ algorithm, the maize seeds with 0, 3 and 7 aged days respectively were used for comparison test which results show that the error between the measured value and actual value is not more than 5%. Among them, the error of RTL is maximum because some of roots are blocked by the soil.

Table 3 Data record of the 7th day of the unaged root system
Table 4 Data of the 7th day of aging 3 days of seed
Table 5 Day 7 data records of aging 7 days of seeds

The experiments indicate that, in Fig. 10, Root Number (RTN), Root Length (RTL), Root Width (RTW) and Root Extension Length (REL) of unaged and 14-day-aged maize seeds are decreased from 15.40, 82.40 mm, 1.53 mm and 82.20 mm to 4.58, 38.6 mm, 1.35 mm and 55.20 mm, respectively. After analyzing the data, it can be found that there is an evident negative linear between seed vigor and RTN, RTL, RTW, REL, but there is uncorrelation between REA and maize seed vigor.

Fig. 10
figure 10

Comparison of parameters of roots of different vigorous maize seeds after 7 days of growth. a The comparison table of aging days and root numbers on the 7th day of maize root growth. b The comparison table of aging days and root length on the 7th day of maize root growth. c The comparison table of aging days and root width on the 7th day of maize root growth. d The comparison table of aging days and root extension length on the 7th day of maize root growth. e The comparison table of aging days and root extension angle on the 7th day of maize root growth

The regression relationships between aging day and root growth parameters are as follows:

$$ RTN \, = \, 0.72d^{2} - \, 7.118d \, + \, 21.888 $$
$$ RTW = \, 0.0033d^{3} - \, 0.0307d^{2} + \, 0.036d \, + \, 1.52 $$
$$ RTL = \, 2.0667d^{3} - \, 16.721d^{2} + \, 25.312d \, + \, 72.7 $$
$$ RTET = 0.075d^{3} - 0.7493d^{2} + 1.4957d \, + 7.348 $$

d is the aging day.

By comparison, it shows that the relation curves between RTL, REL and maize seed aging degree are very similar to that of seed vigor and deterioration [35].

Dynamic phenotypic parameters of maize seed

To obtain maize seed vigor information quickly, the dynamic phenotypic parameters of maize seed root system during early stage were used to mine the information about seed vigor.

As shown in Fig. 11, the differences between RTN, RTL, RTW, REL and ageing degree are even more remarkable, moreover, the gaps of RTN, RTL and REL are increasing during 2 to 12 growth days from 1.5, 4.2, 0.2, 0.3 to 14.4, 52.9, 0.4, 4.7 respectively. The average growing speed of RTN, RTL, RTW and REL change from 1.68 per day, 8.80 mm/d, 0.06 mm/d, 9.0 mm/d to 0.70 per day, 4.3 mm/d, 0.05 mm/d and 5.70 mm/d respectively. Based on the analysis, it can be found that there is a remarkable correlation between the above four dynamic phenotypic parameters and maize seed aged day.

Fig. 11
figure 11

Comparison of the growth rate of root parameters of different vigorous maize seeds after 7 days of growth. a The comparison table of grow days and root number of different vigor maize seeds. b The comparison table of grow days and root length of different vigor maize seeds. c The comparison table of grow days and root width of different vigor maize seeds. d The comparison table of grow days and root extension length of different vigor maize seeds

In summary, the static as well as dynamic maize root phenotyping parameters can be used to evaluate the vigor of seeds.


In order to obtain more root image information directly and rapidly as well as acquire the complete root image from the root segments blocked by the soil, the paper developed an ARPC device which can coax root system growing along the inner wall of the container to expose more parts due to narrow growing annulus space and capillary effect on the boundary between the soil and the container. Moreover, the PCJ algorithm was proposed to recover gaps from disconnected root segments for reconstructing a complete root system architecture.

It should be noted that the toroidal surfaces of the container need be covered by black items to avoid sunlight because the root system will grow into the soil due to root negative phototropism trait which will result in root information reduction in the image. Besides, polarizers and pure black soil were used to overcome the reflection from the transparent container and the disturbance from the soil respectively during image sampling. For obtaining high quality image, the professional SLR camera was used in macro mode in a dark environment indoors with artificial light and the parameters such as polarized angle, aperture, ISO and shutter speed need to be adjusted according to ambient luminous intensity.

For improving the processing speed of the algorithm in computer, the short root segments, not more than 30 pixels in our experiment, can be neglected and covered by the calculated joining root.

However, the proposed ARPC device and PCJ algorithm are focusing on the root phenotyping detection of the early growing stage of fibrous root plant owing to the limitation of the narrow growing annulus space and the increasing proportion of root blocked by the soil, whereas CT X-ray can be applied to acquire the complete image of plant root in any growing stage accurately with huge cost.


In order to realize maize root fast non-invasive phenotyping detection, a novel ARPC device was developed to push each root growing along the inner wall of the container to acquire more root information. Moreover, the PCJ algorithm was proposed to recover various gaps among root segments effectively to generate complete root system images. Finally, five physiological parameters such as RTN, RTL, RTW, REL and REA obtained from the recovered root image by using the presented ARPC device combined with PCJ algorithm were applied to evaluate maize seed vigor which results show that there is an evident negative linear between seed vigor and RTN, RTL, RTW, REL. In addition, the correlation between the above four root phenotypic parameters and maize seed aged day is also remarkable. The proposed ARPC device combined with PCJ algorithm has vast application prospect in high throughput root phenotyping area.

Availability of data and materials

Original procedure and data for root system architecture and analysis is available on request.



the number of roots


root length: the total length of all roots


root width: the average width of each root


root extension length, the horizontal length of the root system


root extension angle, the average angle of each root bend


Progressive Corrosion Joining


Annular Root Phenotyping Container


  1. Rogers ED, Benfey PN. Regulation of plant root system architecture: implications for crop advancement. Curr Opin Biotechnol. 2015;32(32C):93–8.

    Article  CAS  Google Scholar 

  2. Iyer-Pascuzzi AS, Benfey PN. Imaging and analysis platform for automatic phenotyping and trait ranking of plant root systems. Plant Physiol. 2010;152(3):1148.

    Article  CAS  Google Scholar 

  3. Bonser AM, Lynch J, Snapp S. Effect of phosphorus deficiency on growth angle of basal roots in Phaseolus vulgaris. New Phytol. 1996;132:281–8.

    Article  CAS  Google Scholar 

  4. Hund A, Trachsel S, Stamp P. Growth of axile and lateral roots of maize: i development of a phenotying platform. Plant Soil. 2009;325:335–49.

    Article  CAS  Google Scholar 

  5. Atkinson JA, Wingen LU, Griffiths M, Pound MP, Gaju O, Foulkes MJ, Gouis JL, Griffiths S, Bennett MJ, King J, et al. Phenotyping pipeline reveals major seedling root growth QTL in hexaploid wheat. J Exp Bot. 2015;66:2283–92.

    Article  CAS  Google Scholar 

  6. Gioia T, Galinski A, Lenz H, Muller C, Lentz J, Heinz K, Briese C, Putz A, Fiorani F, Watt M, et al. GrowScreen-PaGe, a noninvasive, high-throughput phenotyping system based on germination paper to quantify crop phenotypic diversity and plasticity of root traits under varying nutrient supply. Funct Plant Biol. 2017;44:76–93.

    Article  Google Scholar 

  7. Passot S, Gnacko F, Moukouanga D, Lucas M, Guyomarch S, Ortega BM, Atkinson JA, Belko MN, Bennett MJ, Gantet P, et al. Characterization of pearl millet root architecture and anatomy reveals three types of lateral roots. Front Plant Sci. 2016;7:829.

    Article  Google Scholar 

  8. French A, Ubeda-Tomás S, Holman TJ, Bennett MJ, Pridmore T. Highthroughput quantification of root growth using a novel image-analysis tool. Plant Physiol. 2009;150:1784–95.

    Article  CAS  Google Scholar 

  9. Clark RT, MacCurdy RB, Jung JK, Shaff JE, McCouch SR, Aneshansley DJ, Kochian LV. Three-dimensional root phenotyping with a novel imaging and software platform. Plant Physiol. 2011;156:455–65.

    Article  CAS  Google Scholar 

  10. Topp CN, Iyer-Pascuzzi AS, Anderson JT, Lee C-R, Zurek PR, Symonova O, Zheng Y, Bucksch A, Mileyko Y, Galkovskyi T, et al. 3D phenotyping and quantitative trait locus mapping identifycore regions of the rice genome controlling root architecture. Proc Natl Acad Sci USA. 2013;110:E1695–704.

    Article  CAS  Google Scholar 

  11. Iyer-Pascuzzi AS, Symonova O, Mileyko Y, Hao Y, Belcher H, Harer J, Weitz JS, Benfey PN. Imaging and analysis platform for automatic phenotyping and trait ranking of plant root systems. Plant Physiol. 2010;152:1148–57.

    Article  CAS  Google Scholar 

  12. Taylor BN, Beidler KV, Strand AE, et al. Improved scaling of minirhizotron data using an empirically-derived depth of field and correcting for the underestimation of root diameters. Plant Soil. 2014;374(1–2):941–8.

    Article  CAS  Google Scholar 

  13. Joslin JD, Wolfe MH. Disturbances during minirhizotron installation can affect root observation data. Soil Sci Soc Am J. 1999;63(1):218–21.

    Article  CAS  Google Scholar 

  14. Majdi H, Smucker AJM, Persson H. A comparison between minirhizotron and monolith sampling methods for measuring root growth of maize (Zea mays L.). Plant Soil. 1992;147:127–34.

    Article  Google Scholar 

  15. Bodner G, Nakhforoosh A, Arnold T, Leitner D. Hyperspectral imaging: a novel approach for plant root phenotyping. Plant methods. 2018;14:84.

    Article  Google Scholar 

  16. Lobet G, Draye X. Novel scanning procedure enabling the vectorization of entire rhizotron-grown root systems. Plant Methods. 2013;9(1):1.

    Article  Google Scholar 

  17. Koenig C, Wey H, Binkley T. Precision of the XCT 3000 and comparison of densitometric measurements in distal radius scans between XCT 3000 and XCT 2000 peripheral quantitative computed tomography scanners. J Clin Densitom. 2008;11(4):575.

    Article  Google Scholar 

  18. Metzner R, Eggert A, Van Dusschoten D, et al. Direct comparison of MRI and X-ray CT technologies for 3D imaging of root systems in soil: potential and challenges for root trait quantification. Plant Methods. 2015;11(1):17.

    Article  Google Scholar 

  19. Pflugfelder D, Metzner R, Dusschoten D, et al. Non-invasive imaging of plant roots in different soils using magnetic resonance imaging (MRI). Plant Methods. 2017;13(1):102.

    Article  Google Scholar 

  20. Bates GH. A device for the observation of root growth in the soil. Nature. 1937;139(3527):966–7.

    Article  Google Scholar 

  21. van Dusschoten D, et al. Quantitative 3D analysis of plant roots growing in soil using magnetic resonance imaging. Plant Physiol. 2016;170:1176–88.

    Article  Google Scholar 

  22. Mohamed A, Monnier Y, Mao Z, et al. An evaluation of inexpensive methods for root image acquisition when using rhizotrons. Plant Methods. 2017;13(1):11.

    Article  CAS  Google Scholar 

  23. Johnson MG, Tingey DT, Phillips DL, et al. Advancing fine root research with minirhizotrons. Environ Exp Bot. 2001;45(3):263.

    Article  Google Scholar 

  24. Tierney GL, Fahey TJ. Evaluating minirhizotron estimates of fine root longevity and production in the forest floor of a temperate broadleaf forest. Plant Soil. 2001;229(2):167–76.

    Article  CAS  Google Scholar 

  25. Zhang Z, Li X, Zhang J, et al. Root growth dynamics of Caragana korshinskii using minirhizotrons. J Plant Ecol. 2006;30(3):457–64.

    Article  CAS  Google Scholar 

  26. Nadelhoffer KJ, Aber JD, Melillo JM. Fine roots, net primary production, and soil nitrogen availability: a new hypothesis. Ecology. 1985;66(4):1377–90.

    Article  Google Scholar 

  27. Lu W, Wang X, Wang F. Adaptive minirhizotron for pepper roots observation and its installation based on root system architecture traits. Plant Methods. 2019;15(1):29.

    Article  Google Scholar 

  28. Lobet G, Koevoets IT, Noll M, et al. Using a structural root system model to evaluate and improve the accuracy of root image analysis pipelines. Front Plant Sci. 2017;8:447.

    Article  Google Scholar 

  29. Lin P, Li D, Zou Z, et al. Deep convolutional neural network for automatic discrimination between Fragaria × Ananassa flowers and other similar white wild flowers in fields. Plant Methods. 2018;14(1):64.

    Article  Google Scholar 

  30. Chen H, Giuffrida MV, Tsaftaris SA, et al. Root gap correction with a deep Inpainting Model. In: BMVC. 2018. p. 325.

  31. Ghosh D, Kaabouch N. A survey on image mosaicing techniques. J Vis Commun Image Represent. 2016;34:1–11.

    Article  Google Scholar 

  32. Zhang Z. Camera calibration with one-dimensional objects. IEEE Trans Pattern Anal Mach Intell. 2004;26(7):892–9.

    Article  Google Scholar 

  33. Zhang TY, Suen CY. A fast parallel algorithm for thinning digital patterns. Comm ACM. 1984;27:236–9.

    Article  Google Scholar 

  34. Djanaguiraman M, Prasad PVV, Kumari J, et al. Root length and root lipid composition contribute to drought tolerance of winter and spring wheat. Plant Soil. 2019;439(1–2):57–73.

    Article  CAS  Google Scholar 

  35. Marcos-Fiho Julio. Seed vigor testing: an overview of the past, present and future perspective. Scientia Agricola. 2015;72(4):363–74.

    Article  Google Scholar 

Download references


We would like to express our sincere thanks to Dr. Ling Wang and Hui Luo who have supported us in image processing algorithm.


This work was supported by the National Natural Science Foundation of China (No. 11604154), the Natural Science Foundation of Jiangsu Province (No. BK20181315), the Agricultural Machinery Three New Project (No. SZ120170036), the Asia hub on WEF and Agriculture, and the NAU-MSU Joint Project (No. 2017-H-11), and the Key Research Plan of Yangzhou (No. YZ2018038).

Author information

Authors and Affiliations



This method was proposed by WL, YL and YD. WL developed these devices and proposed a progressive patching algorithm. YL and YD completed the collection and analysis of experimental data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Wei Lu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors have given consent for the publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

16 original root images and the corresponding 16 processed images using the presented algorithm in the paper.

Additional file 2.

The key codes of the Progressive Corrosion algorithm.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lu, W., Li, Y. & Deng, Y. Root phenotypic detection of different vigorous maize seeds based on Progressive Corrosion Joining algorithm of image. Plant Methods 15, 137 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: