Se presentan resultados de ejemplo para la segmentación de raíces basada en imágenes RGB e imágenes HS. Para la imagen espectral, se proporciona un ejemplo de mapeo de agua de alta resolución. Finalmente, se muestran resultados que demuestran el contexto científico en el que se aplican los datos de raíces basados en imágenes.
Medición de raíces basada en RGB
Figura 8 muestra una serie temporal de imágenes RGB de la raíz de la variedad de remolacha azucarera Ferrara. Las imágenes revelan algunos artefactos debidos a una iluminación inhomogénea de las cajas radiculares, con áreas más brillantes a lo largo del lado izquierdo y diferente luminosidad en la zona de superposición entre las imágenes superior e inferior.

Figura 8: Serie temporal de crecimiento de raíces a partir de imágenes RGB.
Las imágenes muestran el cultivar de remolacha azucarera Ferrara en diferentes días posteriores a la siembra (DAS). Las imágenes muestran algunos artefactos debidos a una iluminación no homogénea en el lado izquierdo de la imagen y entre las imágenes superiores e inferiores. Barras de escala, 2 cm. Haga clic aquí para ver una versión más grande de esta figura.
Figura 9 proporciona detalles sobre la segmentación de raíces basada en el umbral de color para el cultivar Ferrara en el día después de la siembra (DAS) 35. Como referencia (Figura 9A), se utiliza una imagen binaria en la que todas las raíces fueron rastreadas manualmente con una tableta gráfica. El tiempo requerido para el rastreo manual de todo el sistema radical, completamente desarrollado y denso, de la remolacha azucarera fue de aproximadamente cuatro horas. Figura 9B ofrece una vista detallada de un área seleccionada en la parte superior de la imagen, donde están presentes raíces laterales viejas. Aquí, varios ejes radiculares no son clasificados por el umbral de color. En la parte inferior (Figura 9C), por el contrario, donde predominan las raíces jóvenes blancas, la segmentación basada en el color clasifica correctamente todos los ejes radiculares. El sistema radical binarizado (Figura 9D) muestra un área negra en el lado izquierdo debida a un artefacto de iluminación, que fue definido como región de exclusión antes de realizar el análisis cuantitativo. Figura 9E muestra los histogramas de píxeles correspondientes a características seleccionadas (raíces frente a suelo) para el canal rojo de la imagen RGB de Ferrara en el DAS 35. Los píxeles de las raíces (color azul) muestran claramente tres picos correspondientes a laterales jóvenes brillantes, laterales viejos oscuros y raíz principal. La superposición entre los laterales viejos y el fondo del suelo es muy fuerte, lo que lleva a ejes radiculares no clasificados (cf. Figura 9B).

Figura 9: Segmentación de raíces mediante umbral de color.
(A) Sistema de raíces segmentado manualmente utilizando una tableta gráfica, (B) área con ejes de raíces viejas mal segmentados en la parte superior y (C) ejes jóvenes correctamente segmentados en la parte inferior de la imagen, y (D) imagen binaria obtenida mediante umbralización basada en color. (E) Histogramas de píxeles para características seleccionadas de la imagen RGB. Las raíces están representadas por las barras azules con diferentes tipos de raíces indicados; el suelo está representado por las barras rojas. Barras de escala en A y D, 2 cm; barras de escala en B y C, 1 cm. Haga clic aquí para ver una versión más grande de esta figura.
La longitud total de raíces visibles resultante cuantificada para la imagen de referencia segmentada manualmente es de 1534,1 cm, mientras que la segmentación automatizada basada en color da una longitud total de raíces de 1427,6 cm.
Las imágenes en escala de grises obtenidas mediante iluminación UV no ofrecen una ventaja en el caso mostrado aquí y tuvieron un desempeño inferior en comparación con el umbralizado a color (longitud de raíz: 1679,7 cm). No fue posible segmentar las raíces viejas, y hubo más ruido en la imagen, probablemente debido a la menor intensidad lumínica de las lámparas UV. Sin embargo, en el caso de raíces jóvenes con alta auto-fluorescencia y un sustrato de fondo brillante, la iluminación UV aún puede ser una opción, como se muestra en una imagen obtenida de otro experimento en el que se utilizó arena como sustrato de fondo (Figura 10).

Figura 10: Iluminación con luz ultravioleta para visualizar raíces sobre un fondo claro.
Ejemplo de un sistema de raíces de trigo duro creciendo en una caja rizosférica llena de arena de cuarzo. La caja rizosférica se imagina con iluminación para (A) luz UV y (B) luz fluorescente (diurna). Barras de escala, 2 cm. Haga clic aquí para ver una versión más grande de esta figura.
Mediciones de raíces basadas en HSI
Figura 11 muestra los espectros promedio para tres regiones de interés en la raíz (laterales viejas y jóvenes, raíz principal) y dos regiones de interés en el suelo (parte superior e inferior del cajón radicular).

Figura 11: Espectros medios de raíz y suelo.
Espectros provenientes de regiones de interés (ROI) en la raíz (tres tipos de raíces) y en el suelo (parte superior e inferior del rizobox). Las ROI se seleccionan para determinar un criterio óptimo de segmentación entre raíz y suelo. Haga clic aquí para ver una versión más grande de esta figura.
Es evidente que las raíces principales y laterales jóvenes difieren considerablemente del fondo en la intensidad de la mayoría de las bandas espectrales. En el caso de las laterales viejas, las diferencias de intensidad son mucho menores. Una característica que puede inferirse visualmente es la pendiente diferente del espectro alrededor de la región de absorción de agua (1450 nm). Aquí, la pendiente de los espectros de las raíces es mayor en comparación con los espectros del suelo. Además, puede identificarse un cambio en los espectros de la raíz principal y de las laterales jóvenes en la región alrededor de 1100 nm, que no ocurre en las laterales viejas.
Figura 12A muestra el resultado del algoritmo de búsqueda que identifica una relación espectral con el mayor contraste entre primer plano y fondo. La relación de los espectros a 1476 nm y 1076 nm proporciona la mejor separación entre raíces y suelo. El histograma resultante de los píxeles de primer plano de las raíces y del fondo del suelo se muestra en la Figura 12B. Aunque existe cierta superposición, la mayoría de los píxeles están claramente separados del fondo del suelo. Al ajustar una curva gaussiana bimodal al histograma y utilizar la distancia de Bhattacharyya para la cuantificación, se obtiene un valor de 7,80. Un valor superior a 3,0 indica un fuerte contraste de imagen que permite una separación confiable28.

Figura 12: Diferencia en reflectancia entre el primer plano de la raíz y el fondo del suelo para diferentes relaciones de bandas espectrales e histograma de píxeles en la relación espectral utilizada para la segmentación.
(A) Los colores brillantes (amarillo) indican alto contraste entre el primer plano y el fondo, mientras que los colores oscuros (azul) indican bajo contraste. Las primeras 15 bandas se han eliminado debido al ruido. Las líneas rojas indican la relación de bandas con el mayor contraste. (B) Histograma de píxeles de las raíces (azul) y del suelo (rojo) en la relación espectral de segmentación. Las barras azules representan la raíz y las barras rojas representan el suelo. El valor de intensidad corresponde a la relación entre la banda espectral 160 y la banda espectral 49. Haga clic aquí para ver una versión más grande de esta figura.
La imagen binaria (Figura 13) se crea aplicando un umbral de intensidad global de la razón espectral identificada con un valor de 1,008 calculado a partir de la distancia del histograma27. El análisis de la longitud de las raíces de esta imagen arroja una longitud total de 1557,3 cm, lo que representa un error de solo el 1,5 % en comparación con la imagen de referencia trazada manualmente.

Figura 13: Imagen binaria del sistema radicular de la variedad de remolacha azucarera Ferrara.
La imagen se obtiene aplicando un umbral espectral global. Barra de escala (esquina inferior izquierda), 2 cm. Haga clic aquí para ver una versión más grande de esta figura.
Aunque la segmentación de raíces ha mejorado mediante el uso de información espectral en comparación con la información basada en color, la intención principal de la imagen hiperespectral es el análisis de las propiedades quimiométricas de la imagen. Esto se ejemplifica mediante la representación del contenido de agua de una imagen de una caja de rizosfera.
Figura 14A muestra los espectros medios de los compartimentos en la caja radicular de calibración (cf. Figura 7) llenos con suelo de diferente contenido de agua. La forma de los espectros es similar entre los compartimentos, es decir, en este caso una razón espectral no proporciona necesariamente un criterio de clasificación más estable. Por lo tanto, se identifica la intensidad en una única banda espectral (1680 nm), donde la diferencia media entre contenidos de agua adyacentes es máxima, como el mejor criterio de separación. Los histogramas resultantes de píxeles para esta longitud de onda espectral se muestran en la Figura 14B.

Figura 14: Características espectrales para la calibración del contenido de agua.
(A) Espectros promedio de nueve compartimentos de agua del rizobox de calibración con diferentes contenidos de agua; (B) Histogramas de píxeles para los compartimentos de agua en la banda 216, donde la distancia promedio entre compartimentos adyacentes es máxima. Haga clic aquí para ver una versión más grande de esta figura.
La relación entre la intensidad promedio de píxeles a 1680 nm y el contenido de agua medido se muestra en la Figura 15.

Figura 15: Relación entre reflectancia espectral y contenido volumétrico de agua.
La figura muestra pares de datos de contenido de agua medido y reflectancia espectral con curvas empíricas (lineal y exponencial) ajustadas a los datos, excluyendo los contenidos más altos de agua (triángulos rojos). Haga clic aquí para ver una versión más grande de esta figura.
La diferenciación de contenidos de agua más altos a partir de la intensidad espectral resulta difícil. Se puede ajustar una regresión significativa (ya sea lineal o exponencial) con un alto R2 para contenidos de agua hasta aproximadamente 0,30 cm3 cm-3. Las condiciones de suelo más húmedo no pueden predecirse de forma confiable mediante el valor de intensidad. Un comportamiento similar de una relación exponencial entre reflectividad y contenido de agua, con una respuesta decreciente para contenidos de agua superiores a 0,30 cm3 cm-3, también se ha observado en otros estudios30.
Se muestra una imagen de una caja radicular con una cartografía detallada del contenido de agua en la Figura 16. Se deben destacar cuatro aspectos. Primero, puede observarse una región con menor contenido de agua en las partes enraizadas de la caja radicular. Segundo, la mayor disminución se concentra en las proximidades de ejes radiculares individuales. Tercero, las zonas de agotamiento también ocurren donde no se ven ejes radiculares en la superficie, lo que indica regiones donde las raíces están ocultas en el suelo. Cuarto, la cartografía del agua sin procesamiento adicional de la imagen da como resultado un aspecto irregular debido a la agregación del suelo. Esto puede indicar una distribución heterogénea del contenido de agua a escala de agregado, pero también puede deberse a efectos de la morfología superficial sobre la calidad de la imagen. Las técnicas quimiométricas de procesamiento de imágenes son una opción para superar tales efectos morfológicos en imágenes espectrales31, pero hasta ahora no se han implementado en los scripts de Matlab utilizados aquí.

Figura 16: Mapeo del contenido de agua en una caja de raíces (rhizobox).
Los colores azul oscuro representan regiones con alto contenido de agua, mientras que las áreas de verde a rojo muestran regiones con bajo contenido de agua. La raíz de la planta aparece superpuesta en la imagen en color negro. Barra de escala (esquina inferior izquierda), 2 cm. Haga clic aquí para ver una versión más grande de esta figura.
Ejemplos de aplicación
Figura 17 relaciona rasgos cuantitativos de las raíces obtenidos mediante análisis de imágenes con mediciones de la parte aérea.

Figura 17: Ejemplos típicos de aplicaciones para datos de raíces.
(A) y (B) muestran información de raíces utilizada para la caracterización de partes aéreas y subterráneas de la planta en un contexto de fenotipado. (A) representa el crecimiento de la raíz de la variedad de remolacha azucarera Ferrara, (B) compara seis variedades de remolacha azucarera cultivadas en cajas de rizosfera utilizando la relación entre el área foliar y el área de raíz (datos de una réplica). (C) y (D) son relaciones funcionales entre caracteres tal como se observan en investigaciones de fisiología vegetal. (C) muestra la influencia de la relación entre el área foliar y el área de raíz en la producción de materia seca y (D) la relación entre el área superficial de la raíz y la conductancia estomática. Haga clic aquí para ver una versión más grande de esta figura.
Figura 17A y 17B son relevantes para la fenotipificación centrada en la caracterización completa de las partes aéreas y subterráneas de la planta. Figura 17A muestra el crecimiento radicular del cultivar de remolacha azucarera Ferrara (cf. Figura 8 para imágenes). La expansión del sistema radicular indica la capacidad de un cultivar para explorar el volumen del suelo en un período de tiempo determinado del ciclo vegetativo. Figura 17B muestra la relación entre el área superficial de las hojas y las raíces de seis cultivares de remolacha azucarera, proporcionando un descriptor del equilibrio entre el suministro (raíz) y la demanda (hoja) de la planta.
Figuras 17C y 17D muestran ejemplos de relaciones funcionales de interés en la investigación fisiológica. En la Figura 17C, la relación entre el área foliar y el área radicular se relaciona con la materia seca formada durante el experimento, lo que indica el papel predominante de la superficie asimiladora como factor limitante para la acumulación de materia seca. La falta de significancia, a pesar de un valor de R2 comparativamente alto, se debe al reducido número de datos apareados (n=6) utilizados aquí. La Figura 17D revela que los cultivares con mayor área superficial radicular (absorción mejorada) presentan una conductancia estomática promedio más alta a lo largo del experimento. El mayor área radicular aparentemente sostiene la extracción de agua, prolongando así la apertura estomática.