Vengono presentati risultati di esempio per la segmentazione delle radici basata su immagini RGB e HS. Per l'immagine spettrale, viene fornito un esempio di mappatura ad alta risoluzione dell'acqua. Infine, vengono mostrati risultati che dimostrano il contesto scientifico in cui i dati delle radici basati su immagini vengono applicati.
Misurazione della radice basata su RGB
Figura 8 mostra una serie temporale di immagini RGB delle radici della coltivar di barbabietola da zucchero Ferrara. Le immagini rivelano alcuni artefatti derivanti da un'illuminazione non omogenea delle rizocasse, con aree più luminose lungo il lato sinistro e diverse intensità luminose nella zona di sovrapposizione tra le immagini superiori e inferiori.

Figura 8: Serie temporale della crescita delle radici ottenuta dall'imaging RGB.
Le immagini mostrano la coltivar di barbabietola da zucchero Ferrara in diversi giorni dopo la semina (DAS). Le immagini evidenziano alcuni artefatti dovuti a un'illuminazione non omogenea sul lato sinistro dell'immagine e tra le immagini superiori e inferiori. Le barre della scala corrispondono a 2 cm. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Figura 9 mostra i dettagli sulla segmentazione delle radici basata sulla soglia di colore per la varietà Ferrara al giorno dopo la semina (DAS) 35. Come riferimento (Figura 9A), viene utilizzata un'immagine binaria in cui tutte le radici sono state tracciate manualmente mediante una tavoletta grafica. Il tempo necessario per il tracciamento manuale dell'intero sistema radicale, completamente sviluppato e fitto, della barbabietola da zucchero è stato di circa quattro ore. Figura 9B fornisce un'immagine dettagliata di un'area selezionata nella parte superiore dell'immagine, dove sono presenti radici laterali vecchie. In questo caso, diversi assi radicali non vengono classificati sulla base della soglia di colore. Al contrario, nella parte inferiore (Figura 9C), dove predominano le radici giovani di colore bianco, la segmentazione basata sul colore classifica correttamente tutti gli assi radicali. Il sistema radicale binarizzato (Figura 9D) mostra un'area nera sul lato sinistro, dovuta a un artefatto di illuminazione, che è stata definita come regione di esclusione prima dell'analisi quantitativa. Figura 9E mostra gli istogrammi dei pixel corrispondenti alle caratteristiche selezionate (radici vs. suolo) per il canale rosso dell'immagine RGB di Ferrara al DAS 35. I pixel delle radici (colore blu) mostrano chiaramente tre picchi corrispondenti alle laterali giovani luminose, alle laterali vecchie scure e alla radice principale. L'overlap tra le laterali vecchie e lo sfondo del suolo è molto intenso, causando assi radicali non classificati (cf. Figura 9B).

Figura 9: Segmentazione delle radici mediante soglia cromatica.
(A) Sistema radicale segmentato manualmente utilizzando una tavoletta grafica, (B) area con assi radicali vecchi mal segmentati nella parte superiore e (C) assi giovani correttamente segmentati nella parte inferiore dell'immagine, e (D) immagine binaria ottenuta mediante sogliatura basata sul colore. (E) Istogrammi dei pixel per le caratteristiche selezionate dell'immagine RGB. Le radici sono rappresentate dalle barre blu con indicazione dei diversi tipi radicali; il suolo è rappresentato dalle barre rosse. Le barre della scala in A e D, 2 cm; le barre della scala in B e C, 1 cm. Cliccare qui per visualizzare una versione ingrandita di questa figura.
La lunghezza totale delle radici visibili risultante quantificata per l'immagine di riferimento segmentata manualmente è di 1534,1 cm, mentre la segmentazione automatizzata basata sul colore fornisce una lunghezza totale delle radici di 1427,6 cm.
Le immagini in scala di grigi ottenute con illuminazione UV non offrono vantaggi nel caso qui mostrato, risultando peggiori rispetto alla sogliatura del colore (lunghezza delle radici: 1679,7 cm). Le radici vecchie non sono state segmentate correttamente e nell'immagine era presente un rumore maggiore, probabilmente dovuto all'intensità luminosa inferiore delle lampade UV. Tuttavia, nel caso di radici giovani con elevata autofluorescenza e un substrato di fondo chiaro, l'illuminazione UV può comunque rappresentare un'opzione, come mostrato da un'immagine ottenuta da un altro esperimento in cui la sabbia è stata utilizzata come substrato di fondo (Figura 10).

Figura 10: Illuminazione UV per visualizzare le radici su sfondo chiaro.
Esempio di un sistema radicale di grano duro in crescita in una rhizobox riempita con sabbia di quarzo. La rhizobox è stata immaginata utilizzando un'illuminazione a (A) luce UV e (B) luce fluorescente (giorno). Barre della scala, 2 cm. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Misurazioni della radice basate sull'HSI
Figura 11 mostra gli spettri medi per tre aree di interesse (ROI) delle radici (laterali vecchie e giovani, radice principale) e due ROI del suolo (parte superiore e inferiore della rizocassetta).

Figura 11: Spettri medi di radice e suolo.
Spettri provenienti da regioni di interesse (ROI) sulla radice (tre tipi di radice) e nel suolo (parte superiore e inferiore del rizobox). Le ROI vengono selezionate per determinare un criterio ottimale di segmentazione tra radice e suolo. Cliccare qui per visualizzare una versione ingrandita di questa figura.
È evidente che le radici principali e le radici laterali giovani differiscono sostanzialmente dal fondo per l'intensità della maggior parte delle bande spettrali. Per le radici laterali vecchie, le differenze di intensità sono molto più basse. Una caratteristica che può essere rilevata visivamente è la diversa pendenza dello spettro nella regione di assorbimento dell'acqua (1450 nm). In questo caso, la pendenza degli spettri delle radici è maggiore rispetto a quella degli spettri del suolo. Inoltre, si può osservare un cambiamento negli spettri della radice principale e delle radici laterali giovani nella regione intorno a 1100 nm, che non si verifica nelle radici laterali vecchie.
Figura 12A mostra il risultato dell'algoritmo di ricerca che identifica un rapporto spettrale con il massimo contrasto tra primo piano e sfondo. Il rapporto tra gli spettri a 1476 nm e 1076 nm fornisce la migliore separazione tra radici e suolo. L'istogramma risultante dei pixel relativi al primo piano (radici) e allo sfondo (suolo) è mostrato in Figura 12B. Sebbene vi sia un certo sovrapposizione, la maggior parte dei pixel è chiaramente distinta dallo sfondo del suolo. Adattando una curva gaussiana bimodale all'istogramma e utilizzando la distanza di Bhattacharyya per la quantificazione, si ottiene un valore di 7,80. Un valore superiore a 3,0 indica un forte contrasto dell'immagine, che permette una separazione affidabile28.

Figura 12: Differenza di riflettanza tra primo piano radicale e sfondo del suolo per diversi rapporti spettrali e istogramma dei pixel al rapporto spettrale utilizzato per la segmentazione.
(A) Colori chiari (giallo) indicano un alto contrasto tra primo piano e sfondo, colori scuri (blu) indicano un basso contrasto. I primi 15 canali sono stati rimossi a causa del rumore. Le linee rosse indicano il rapporto spettrale con il contrasto più elevato. (B) Istogramma dei pixel delle radici (blu) e del suolo (rosso) al rapporto spettrale di segmentazione. Le barre blu rappresentano le radici e le barre rosse il suolo. Il valore di intensità corrisponde al rapporto tra il canale spettrale 160 e il canale spettrale 49. Cliccare qui per visualizzare una versione ingrandita di questa figura.
L'immagine binaria (Figura 13) viene creata applicando una soglia di intensità globale al rapporto spettrale identificato, con un valore di 1,008 calcolato dalla distanza dell'istogramma27. L'analisi della lunghezza delle radici di questa immagine fornisce una lunghezza totale di 1557,3 cm, che corrisponde a un errore del solo 1,5% rispetto all'immagine di riferimento tracciata manualmente.

Figura 13: Immagine binaria del sistema radicale della cultivar di barbabietola da zucchero Ferrara.
L'immagine è ottenuta applicando una soglia spettrale globale. Barra della scala (angolo in basso a sinistra), 2 cm. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Sebbene la segmentazione delle radici sia migliorata grazie alle informazioni spettrali rispetto a quelle basate sul colore, l'obiettivo principale dell'imaging iperspettrale è l'analisi delle proprietà chimio metriche delle immagini. Questo viene illustrato tramite la mappatura del contenuto di acqua di un'immagine di una rhizobox.
Figura 14A mostra gli spettri medi delle compartimentazioni nella rizocassetta di calibrazione (cf. Figura 7) riempita con suolo a diverso contenuto idrico. La forma degli spettri è simile tra le diverse compartimentazioni, pertanto un rapporto spettrale non fornisce necessariamente un criterio di classificazione più stabile. Di conseguenza, l'intensità in una singola banda spettrale (1680 nm), in corrispondenza della quale la differenza media tra contenuti idrici adiacenti è massimizzata, viene identificata come il criterio migliore per la separazione. Gli istogrammi risultanti dei pixel per questa lunghezza d'onda spettrale sono mostrati in Figura 14B.

Figura 14: Caratteristiche spettrali per la calibrazione del contenuto di acqua.
(A) Spettri medi di nove compartimenti di acqua dalla scatola radicale di calibrazione con diversi contenuti di acqua; (B) Istogrammi dei pixel per i compartimenti di acqua alla banda 216, dove la distanza media tra compartimenti adiacenti è massima. Cliccare qui per visualizzare una versione ingrandita di questa figura.
La relazione tra l'intensità media del pixel a 1680 nm e il contenuto di acqua misurato è mostrata nella Figura 15.

Figura 15: Relazione tra riflettanza spettrale e contenuto volumetrico di acqua.
La figura mostra coppie di dati relative al contenuto di acqua misurato e alla riflettanza spettrale, con curve empiriche (lineari ed esponenziali) adattate ai dati escludendo i contenuti di acqua più elevati (triangoli rossi). Cliccare qui per visualizzare una versione ingrandita di questa figura.
La differenziazione di contenuti idrici più elevati a partire dall'intensità spettrale diventa difficile. Si può adattare una regressione significativa (lineare o esponenziale) con un alto valore di R2 ai contenuti idrici fino a circa 0,30 cm3 cm-3. Le condizioni di suolo più umido non possono essere previste in modo affidabile a partire dal valore di intensità. Un comportamento simile, con una relazione esponenziale tra riflettività e contenuto idrico e una risposta decrescente per contenuti idrici superiori a 0,30 cm3 cm-3, è stato riscontrato anche in altri studi30.
Un'immagine di una rizocassa con una mappatura dettagliata del contenuto idrico è mostrata in Figura 16. Vanno evidenziati quattro aspetti. Primo, si può osservare una regione con contenuto idrico inferiore nelle parti della rizocassa con radici. Secondo, il maggiore esaurimento è concentrato nelle vicinanze di assi radicali singoli. Terzo, si verificano zone di deplezione anche dove non sono visibili assi radicali sulla superficie, indicando regioni in cui le radici sono nascoste nel suolo. Quarto, la mappatura dell'acqua senza ulteriore elaborazione dell'immagine produce un aspetto irregolare a causa dell'aggregazione del suolo. Ciò può indicare una distribuzione non omogenea del contenuto idrico alla scala dell'aggregato, ma anche effetti della morfologia superficiale sulla qualità dell'immagine. Tecniche chemiometriche di elaborazione delle immagini rappresentano un'opzione per superare tali effetti morfologici nelle immagini spettrali31, ma non sono ancora state implementate negli script Matlab utilizzati in questo caso.

Figura 16: Mappatura del contenuto idrico in una rhizobox.
Le tonalità blu scuro rappresentano aree con elevato contenuto idrico, mentre le zone verdi fino al rosso indicano regioni con basso contenuto idrico. La radice della pianta è sovrapposta all'immagine in nero. Scala grafica (angolo in basso a sinistra), 2 cm. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Esempi di applicazione
Figura 17 mette in relazione tratti radicali quantitativi ottenuti dall'analisi delle immagini con misurazioni della parte aerea.

Figura 17: Esempi tipici di applicazione per i dati radicolari.
(A) e (B) mostrano informazioni sulle radici utilizzate per la caratterizzazione della parte aerea e sotterranea della pianta in un contesto di fenotipizzazione. (A) rappresenta la crescita radicale della cultivar di barbabietola da zucchero Ferrara, (B) confronta sei cultivar di barbabietola da zucchero coltivate in rizobox utilizzando il rapporto tra area fogliare e area radicale (dati provenienti da un replicato). (C) e (D) mostrano relazioni funzionali tra caratteri come quelle riscontrate nella ricerca fisiologica vegetale. (C) illustra l'influenza del rapporto tra area fogliare e area radicale sulla produzione di materia secca e (D) la relazione tra superficie radicale e conduttanza stomatica. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Figura 17A e 17B sono pertinenti per la fenotipizzazione, con un focus sulla caratterizzazione completa della parte aerea e della parte sotterranea della pianta. Figura 17A mostra l'accrescimento radicale del cultivar di barbabietola da zucchero Ferrara (cf. Figura 8 per le immagini). L'espansione dell'apparato radicale indica la capacità di un cultivar di esplorare il volume del suolo in un determinato arco temporale del periodo vegetativo. Figura 17B mostra il rapporto tra superficie fogliare e superficie radicale di sei cultivar di barbabietola da zucchero, fornendo un indice dell'equilibrio tra fornitura (radice) e domanda (foglia).
Figure 17C e 17D mostrano esempi di relazioni funzionali di interesse nella ricerca fisiologica. In Figura 17C il rapporto tra superficie fogliare e radicale è correlato alla massa secca formata durante l'esperimento, indicando il ruolo predominante della superficie assimilante come fattore limitante per l'accumulo di massa secca. La mancanza di significatività, nonostante un valore di R2 relativamente alto, è legata al ridotto numero di dati accoppiati (n=6) utilizzati in questo caso. Figura 17D rivela che i cultivar con una maggiore superficie radicale (assorbimento migliorato) presentano una conduttanza stomatica media più elevata nel corso dell'esperimento. L'area radicale maggiore sostiene evidentemente l'estrazione d'acqua, prolungando così l'apertura degli stomati.