Des résultats d'exemple sont présentés pour la segmentation des racines basée sur l'imagerie RVB et HS. Pour l'imagerie spectrale, un exemple de cartographie de l'eau à haute résolution est fourni. Enfin, des résultats sont montrés afin de démontrer le contexte scientifique dans lequel les données racinaires issues de l'imagerie sont appliquées.
Mesure racinaire basée sur RGB
Figure 8 montre une série temporelle d'images racinaires RVB du cultivar de betterave sucrière Ferrara. Les images révèlent certains artefacts dus à une illumination hétérogène des rhizoboxes, avec des zones plus lumineuses le long du côté gauche et des différences de luminosité dans la zone de recouvrement entre les images supérieures et inférieures.

Figure 8 : Série chronologique de la croissance racinaire issue de l'imagerie RVB.
Les images montrent le cultivar de betterave sucrière Ferrara à différents jours après le semis (JAS). On observe certains artefacts dus à un éclairage non homogène du côté gauche de l'image ainsi qu'entre les images du haut et du bas. Les échelles correspondent à 2 cm. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Figure 9 fournit des détails sur la segmentation des racines basée sur le seuil de couleur pour le cultivar Ferrara au 35e jour après semis (JAS). À titre de référence (Figure 9A), une image binaire est utilisée où toutes les racines ont été suivies manuellement à l’aide d’une tablette graphique. Le temps nécessaire pour le suivi manuel de l’ensemble du système racinaire dense de betterave sucrière, complètement développé, était d’environ quatre heures. Figure 9B donne une vue détaillée d'une zone sélectionnée en haut de l'image où sont présentes d'anciennes racines latérales. Ici, plusieurs axes racinaires ne sont pas classifiés par le seuil de couleur. En bas (Figure 9C) au contraire, lorsque les jeunes racines blanches sont prédominantes, la segmentation fondée sur la couleur classe correctement tous les axes racinaires. Le système racinaire binarisé (Figure 9D) montre une zone noire sur le côté gauche due à un artefact d'illumination, définie comme région d'exclusion avant l'analyse quantitative. Figure 9E montre les histogrammes de pixels correspondants des caractéristiques sélectionnées (racines vs. sol) pour le canal rouge de l'image RVB de Ferrare à 35 JDS. Les pixels racinaires (en bleu) présentent clairement trois pics correspondant aux jeunes latérales claires, aux anciennes latérales foncées et à la racine pivot. Le recouvrement entre les anciennes latérales et le fond du sol est très important, conduisant à des axes racinaires non classifiéscf. Figure 9B).

Figure 9 : Segmentation des racines à l'aide d'un seuil de couleur.
(A) Système racinaire segmenté manuellement à l’aide d’une tablette graphique, (B) zone présentant une mauvaise segmentation des axes racinaires anciens en haut de l’image et (C) axes jeunes correctement segmentés en bas de l’image, et (D) image binaire obtenue par seuillage basé sur la couleur. (E) Histogrammes des pixels pour les caractéristiques sélectionnées de l’image RVB. Les racines sont représentées par les barres bleues, avec différents types de racines indiqués ; le sol est représenté par les barres rouges. Les barres d’échelle dans A et D, 2 cm ; les barres d’échelle dans B et C, 1 cm. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
La longueur totale des racines visibles obtenue pour l'image de référence segmentée manuellement est de 1534,1 cm, tandis que la segmentation automatisée basée sur la couleur donne une longueur totale des racines de 1427,6 cm.
Les images en niveaux de gris provenant de l'illumination UV n'offrent pas d'avantage dans le cas présenté ici et se sont révélées moins performantes par rapport au seuillage en couleur (longueur des racines : 1679,7 cm). Les vieilles racines n'ont pas pu être segmentées, et l'image comportait davantage de bruit, probablement en raison de l'intensité lumineuse plus faible des lampes UV. Toutefois, dans le cas de jeunes racines présentant une forte auto-fluorescence et un substrat d'arrière-plan lumineux, l'illumination UV peut néanmoins constituer une option, comme le montre une image obtenue lors d'une autre expérience utilisant du sable comme substrat d'arrière-plan (Figure 10).

Figure 10 : Éclairage UV pour visualiser les racines sur un fond clair.
Exemple d'un système racinaire de blé dur en cours de croissance dans une rhizoboîte remplie de sable de quartz. La rhizoboîte est imagée avec un éclairage en lumière (A) UV et (B) fluorescente (lumière du jour). Les barres d'échelle représentent 2 cm. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Mesures racinaires basées sur l'HSI
Figure 11 présente les spectres moyens pour trois régions d'intérêt racinaires (latérale ancienne et jeune, racine pivotante) et deux régions d'intérêt du sol (partie supérieure et inférieure du rhizobox).

Figure 11 : Spectres moyens des racines et du sol.
Spectres provenant de régions d'intérêt (ROIs) sur les racines (trois types de racines) et dans le sol (partie supérieure et inférieure de la rhizoboîte). Les ROIs sont sélectionnées afin de déterminer un critère optimal de segmentation entre racine et sol. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Il est évident que les racines pivotantes et latérales jeunes diffèrent sensiblement du fond par l'intensité de la plupart des bandes spectrales. Pour les latérales anciennes, les différences d'intensité sont beaucoup plus faibles. Une caractéristique pouvant être déduite visuellement est la pente différente du spectre autour de la région d'absorption de l'eau (1450 nm). Ici, la pente des spectres racinaires est plus élevée par rapport aux spectres du sol. De plus, on peut identifier une modification des spectres de la racine pivotante et des latérales jeunes dans la région autour de 1100 nm, qui ne se produit pas chez les latérales anciennes.
Figure 12A montre le résultat de l'algorithme de recherche identifiant un rapport spectral présentant le meilleur contraste entre premier plan et arrière-plan. Le rapport des spectres à 1476 nm et 1076 nm assure la meilleure séparation entre les racines et le sol. L'histogramme résultant des pixels correspondant aux racines (premier plan) et au sol (arrière-plan) est présenté dans Figure 12B. Bien qu'il existe un certain recouvrement, la majorité des pixels sont clairement distincts du fond correspondant au sol. En ajustant une courbe gaussienne bimodale à l'histogramme et en utilisant la distance de Bhattacharyya pour la quantification, une valeur de 7,80 est obtenue. Une valeur supérieure à 3,0 indique un fort contraste d'image permettant une séparation fiable28.

Figure 12 : Différence de réflectance entre le premier plan racinaire et le fond de sol pour différents rapports de bandes spectrales et histogramme des pixels au rapport spectral utilisé pour la segmentation.
(A) Les couleurs vives (jaune) indiquent un fort contraste entre le premier plan et l'arrière-plan, tandis que les couleurs sombres (bleu) indiquent un faible contraste. Les 15 premières bandes ont été supprimées en raison du bruit. Les lignes rouges indiquent le rapport de bande présentant le contraste le plus élevé. (B) Histogramme des pixels des racines (bleu) et du sol (rouge) au rapport spectral de segmentation. Les barres bleues représentent les racines et les barres rouges représentent le sol. La valeur d'intensité correspond au rapport entre la bande spectrale 160 et la bande spectrale 49. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
L'image binaire (Figure 13) est créée en appliquant un seuil d'intensité global du rapport spectral identifié à une valeur de 1,008 calculée à partir de la distance de l'histogramme27. L'analyse de la longueur des racines de cette image donne une longueur totale de 1557,3 cm, ce qui représente une erreur de seulement 1,5 % par rapport à l'image de référence tracée manuellement.

Figure 13 : Image binaire du système racinaire du cultivar de betterave sucrière Ferrara.
L'image est obtenue en appliquant un seuil spectral global. Échelle (coin inférieur gauche), 2 cm. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Bien que la segmentation des racines se soit améliorée grâce aux informations spectrales par rapport aux informations basées sur la couleur, l'objectif principal de l'imagerie hyperspectrale est l'analyse des propriétés chimiométriques de l'image. Cela est illustré par la cartographie de la teneur en eau d'une image de rhizobox.
Figure 14A montre les spectres moyens des compartiments dans la rhizobox de calibration (cf. Figure 7) remplie avec un sol présentant différentes teneurs en eau. La forme des spectres est similaire entre les compartiments, c’est-à-dire qu’un rapport spectral ne fournit pas nécessairement un critère de classification plus stable. Ainsi, l’intensité à une seule bande spectrale (1680 nm), où la différence moyenne entre des teneurs en eau adjacentes est maximisée, est identifiée comme le meilleur critère de séparation. Les histogrammes de pixels résultants pour cette longueur d’onde spectrale sont présentés dans Figure 14B.

Figure 14 : Caractéristiques spectrales pour l'étalonnage de la teneur en eau.
(A) Spectres moyens de neuf compartiments d'eau du rhizobox d'étalonnage, ayant différentes teneurs en eau ; (B) Histogrammes des pixels pour les compartiments d'eau à la bande 216, où la distance moyenne entre les compartiments adjacents est maximale. Veuillez cliquer ici pour consulter une version agrandie de cette figure.
La relation entre l'intensité moyenne des pixels à 1680 nm et la teneur en eau mesurée est illustrée dans la Figure 15.

Figure 15 : Relation entre la réflectance spectrale et la teneur en eau volumétrique.
La figure présente des paires de données de teneur en eau mesurée et de réflectance spectrale, avec des courbes empiriques (linéaire et exponentielle) ajustées aux données, à l'exception des teneurs en eau les plus élevées (triangles rouges). Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
La différenciation de teneurs en eau élevées à partir de l'intensité spectrale devient difficile. Une régression significative (linéaire ou exponentielle) avec un R2 élevé peut être ajustée aux teneurs en eau jusqu'à environ 0,30 cm3 cm-3. Les conditions de sol plus humides ne peuvent pas être prédites de façon fiable par la valeur d'intensité. Un comportement similaire, avec une relation exponentielle entre la réflectivité et la teneur en eau et une réponse décroissante pour des teneurs en eau supérieures à 0,30 cm3 cm-3, a également été observé dans d'autres études30.
Une image de rhizobox avec une cartographie fine de la teneur en eau est présentée dans la Figure 16. Quatre aspects doivent être soulignés. Premièrement, une zone de teneur en eau plus faible est visible dans les parties racinées de la rhizobox. Deuxièmement, la déplétion la plus marquée est concentrée à proximité des axes racinaires individuels. Troisièmement, des zones de déplétion apparaissent également là où aucun axe racinaire n'est visible en surface, indiquant des régions où les racines sont enfouies dans le sol. Quatrièmement, la cartographie de l'eau sans traitement d'image supplémentaire donne une apparence tachetée en raison de l'agrégation du sol. Cela peut refléter une distribution hétérogène de la teneur en eau à l'échelle de l'agrégat, mais aussi l'effet de la morphologie de surface sur la qualité de l'image. Des techniques de traitement d'image chimiométriques permettent de surmonter de tels effets morphologiques dans les images spectrales31, mais n'ont pas encore été intégrées aux scripts Matlab utilisés ici.

Figure 16 : Cartographie de la teneur en eau dans une rhizobox.
Les couleurs bleu foncé représentent les régions à forte teneur en eau, tandis que les zones allant du vert au rouge indiquent une faible teneur en eau. La racine de la plante est superposée en noir sur l'image. Échelle (en bas à gauche), 2 cm. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Exemples d'application
Figure 17 met en relation les caractéristiques racinaires quantitatives issues de l'analyse d'image avec les mesures aériennes.

Figure 17 : Exemples d'applications typiques pour les données racinaires.
(A) et (B) montrent des informations racinaires utilisées pour la caractérisation des parties aériennes et souterraines de la plante dans un contexte de phénotypage. (A) représente la croissance racinaire du cultivar de betterave sucrière Ferrara, (B) compare six cultivars de betterave sucrière cultivés en rhizobox à l’aide du rapport entre la surface foliaire et la surface racinaire (données provenant d’un seul réplicat). (C) et (D) illustrent des relations fonctionnelles entre des caractères telles qu’observées en recherche physiologique végétale. (C) montre l’influence du rapport entre la surface foliaire et la surface racinaire sur la production de matière sèche et (D) la relation entre la surface racinaire et la conductance stomatique. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Figure 17A et 17B sont pertinentes pour le phénotypage axé sur une caractérisation complète des parties aériennes et souterraines de la plante. Figure 17A montre la croissance racinaire du cultivar de betterave sucrière Ferrara (cf. Figure 8 pour les images). L'expansion du système racinaire indique la capacité d'un cultivar à explorer le volume du sol sur une période donnée du cycle végétatif. Figure 17B montre le rapport entre la surface foliaire et la surface racinaire de six cultivars de betterave sucrière, fournissant un descripteur de l'équilibre entre l'approvisionnement (racine) et la demande (feuille) de la plante.
Figures 17C et 17D présentent des exemples de relations fonctionnelles d'intérêt en recherche physiologique. Dans la figure 17C, le rapport entre la surface foliaire et la surface racinaire est mis en relation avec la matière sèche formée au cours de l'expérience, indiquant le rôle prédominant de la surface d'assimilation comme facteur limitant l'accumulation de matière sèche. L'absence de significativité, malgré un coefficient de détermination R2 relativement élevé, est liée au faible nombre de données appariées (n=6) utilisé ici. La figure 17D révèle que les cultivars ayant une surface racinaire plus importante (absorption améliorée) présentent en moyenne une conductance stomatique plus élevée tout au long de l'expérience. La surface racinaire plus étendue permet apparemment une extraction d'eau accrue, prolongeant ainsi l'ouverture des stomates.