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 dovuti a 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 radicale ottenuta dall'imaging RGB.
Le immagini mostrano la cultivar di barbabietola da zucchero Ferrara a diversi giorni dopo la semina (DAS). Le immagini mostrano 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 della segmentazione delle radici basata sulla soglia di colore per la varietà Ferrara al giorno dalla 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 denso, della barbabietola da zucchero è stato di circa quattro ore. Figura 9B fornisce una vista 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 prevalgono 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, definita come regione di esclusione prima dell'analisi quantitativa. Figura 9E mostra gli istogrammi corrispondenti dei pixel relativi a 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 a laterali giovani luminose, laterali vecchie scure e radice principale. L'overlap tra le laterali vecchie e lo sfondo del suolo è molto intenso, portando a radici non classificate (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 scarsamente 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. Barre della scala in A e D, 2 cm; 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 e hanno prestato peggio rispetto alla sogliatura a colori (lunghezza delle radici: 1679,7 cm). Le radici vecchie non sono state segmentate e nell'immagine era presente un rumore maggiore, probabilmente a causa dell'intensità luminosa inferiore delle lampade UV. Tuttavia, nel caso di radici giovani con elevata autofluorescenza e un substrato di fondo luminoso, 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 rizocassa riempita con sabbia di quarzo. La rizocassa è stata immaginata con illuminazione per (A) luce UV e (B) luce fluorescente (diurna). Le barre della scala corrispondono a 2 cm. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Misurazioni della radice basate sull'HSI
Figura 11 mostra gli spettri medi di tre aree di interesse radicale (laterale vecchia e giovane, radice principale) e di due aree di interesse del suolo (parte superiore e inferiore della rizobox).

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 sono 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 dedotta visivamente è la diversa pendenza dello spettro nella regione di assorbimento dell'acqua (1450 nm). In questa zona, 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 delle radici e allo sfondo del suolo è mostrato in Figura 12B. Sebbene vi sia un certo sovrapposizione, la maggior parte dei pixel è chiaramente separata 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 la porzione di radice in primo piano e lo sfondo del suolo per diversi rapporti spettrali e istogramma dei pixel al rapporto spettrale utilizzato per la segmentazione.
(A) I colori chiari (giallo) indicano un alto contrasto tra primo piano e sfondo, mentre i colori scuri (blu) indicano un basso contrasto. Le prime 15 bande sono state rimosse a causa del rumore. Le linee rosse indicano il rapporto tra bande 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 la radice e le barre rosse il suolo. Il valore di intensità corrisponde al rapporto tra la banda spettrale 160 e la banda 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 varietà di barbabietola da zucchero Ferrara.
L'immagine è stata 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 si sia migliorata utilizzando informazioni spettrali rispetto a quelle basate sul colore, l'obiettivo principale dell'immagine iperspettrale è l'analisi delle proprietà chimicometriche dell'immagine. Questo viene illustrato tramite la mappatura del contenuto di acqua di un'immagine di una rhizobox.
Figura 14A mostra gli spettri medi dei compartimenti nella rizocassa di calibrazione (cf. Figura 7) riempita con suolo a diverso contenuto idrico. La forma degli spettri è simile tra i compartimenti, pertanto in questo caso 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 dei pixel risultanti 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 idrici; (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 d'acqua.
La figura mostra coppie di dati di contenuto d'acqua misurato e riflettanza spettrale, con curve empiriche (lineare ed esponenziale) adattate ai dati escludendo i contenuti d'acqua più elevati (triangoli rossi). Cliccare qui per visualizzare una versione ingrandita di questa figura.
La differenziazione di contenuti idrici più elevati in base all'intensità spettrale diventa difficile. Una regressione significativa (lineare o esponenziale) con alto valore di R2 può essere adattata ai contenuti idrici fino a circa 0,30 cm3 cm-3. Le condizioni di suolo più umido non possono essere previste in modo affidabile tramite il valore di intensità. Un comportamento simile, caratterizzato da una relazione esponenziale tra riflettività e contenuto idrico con 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 aree 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 eterogenea 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 finora non sono state implementate negli script Matlab utilizzati in questo contesto.

Figura 16: Mappatura del contenuto idrico in un 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 varietà di barbabietola da zucchero Ferrara, (B) confronta sei varietà di barbabietola da zucchero coltivate in rizobox utilizzando il rapporto tra area fogliare e area radicale (dati da un'unica replica). (C) e (D) sono relazioni funzionali tra caratteri come si riscontrano nella ricerca fisiologica vegetale. (C) mostra l'influenza del rapporto tra area fogliare e area radicale sulla produzione di materia secca e (D) la relazione tra area superficiale radicale e conduttanza stomatica. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Figura 17A e 17B sono pertinenti per la fenotipizzazione incentrata sulla caratterizzazione completa della parte aerea e della parte sotterranea della pianta. Figura 17A mostra la crescita radicale del cultivar di barbabietola da zucchero Ferrara (cf. Figura 8 per le immagini). L'espansione del sistema radicale indica la capacità di un cultivar di esplorare il volume del suolo in un determinato lasso di tempo del periodo vegetativo. Figura 17B mostra il rapporto tra superficie fogliare e superficie radicale di sei cultivar di barbabietola da zucchero, fornendo un parametro descrittivo dell'equilibrio tra fornitura (radice) e domanda (foglia) della pianta.
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 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 stomi.