São apresentados resultados de exemplo para segmentação de raízes com base em imagens RGB e HS. Para a imagem espectral, é fornecido um exemplo de mapeamento de água em alta resolução. Por fim, são apresentados resultados que demonstram o contexto científico em que os dados de raízes baseados em imagens são aplicados.
Medição da raiz baseada em RGB
Figura 8 mostra uma série temporal de imagens de raízes em cores RGB da cultivar de beterraba sacarina Ferrara. As imagens revelam alguns artefatos decorrentes da iluminação não homogênea das caixas radiculares, com áreas mais claras ao longo do lado esquerdo e diferentes níveis de brilho na região de sobreposição entre as imagens superiores e inferiores.

Figura 8: Série temporal de crescimento da raiz a partir da imagem RGB.
As imagens mostram o cultivar de beterraba Ferrara em diferentes dias após a semeadura (DAS). As imagens apresentam alguns artefatos devido à iluminação não homogênea no lado esquerdo da imagem e entre as imagens superiores e inferiores. Barras de escala, 2 cm. Clique aqui para visualizar uma versão maior desta figura.
Figura 9 apresenta detalhes sobre a segmentação de raízes com base no limiar de cor para a cultivar Ferrara no dia após a semeadura (DAS) 35. Como referência (Figura 9A), utiliza-se uma imagem binária na qual todas as raízes foram rastreadas manualmente com um Tablet Gráfico. O tempo necessário para o rastreamento manual de todo o sistema radicular denso e completamente desenvolvido da beterraba foi de aproximadamente quatro horas. Figura 9B mostra uma visão detalhada de uma área selecionada na parte superior da imagem, onde estão presentes raízes laterais antigas. Nessa região, vários eixos radiculares não são classificados pelo limiar de cor. Na parte inferior (Figura 9C), ao contrário, onde predominam raízes jovens brancas, a segmentação baseada na cor classifica corretamente todos os eixos radiculares. O sistema radicular binarizado (Figura 9D) exibe uma área preta no lado esquerdo, decorrente de um artefato de iluminação, que foi definida como região de exclusão antes da realização da análise quantitativa. Figura 9E mostra os histogramas de pixels correspondentes às características selecionadas (raízes versus solo) para o canal vermelho da imagem RGB da Ferrara no DAS 35. Os pixels das raízes (cor azul) mostram claramente três picos correspondentes às laterais jovens claras, laterais antigas escuras e raiz principal. A sobreposição entre as laterais antigas e o fundo do solo é muito acentuada, resultando em eixos radiculares não classificados (cf. Figura 9B).

Figura 9: Segmentação de raízes utilizando um limiar de cor.
(A) Sistema radicular segmentado manualmente utilizando uma mesa digitalizadora, (B) área com eixos de raízes antigas mal segmentados na parte superior e (C) eixos jovens adequadamente segmentados na parte inferior da imagem, e (D) imagem binária obtida a partir do limiar baseado em cor. (E) Histogramas de pixels para características selecionadas da imagem RGB. As raízes são representadas pelas barras azuis com diferentes tipos de raízes indicados; o solo é representado pelas barras vermelhas. Barras de escala em A e D, 2 cm; barras de escala em B e C, 1 cm. Clique aqui para visualizar uma versão maior desta figura.
O comprimento total visível das raízes quantificado para a imagem de referência segmentada manualmente é de 1534,1 cm, enquanto a segmentação automatizada baseada em cor fornece um comprimento total de raiz de 1427,6 cm.
Imagens em escala de cinza obtidas por iluminação UV não apresentam vantagem no caso mostrado aqui e tiveram desempenho inferior em comparação com o limiar de cor (comprimento das raízes: 1679,7 cm). As raízes antigas não puderam ser segmentadas, e houve mais ruído na imagem, provavelmente devido à menor intensidade luminosa das lâmpadas UV. No entanto, no caso de raízes jovens com alta autofluorescência e um substrato de fundo claro, a iluminação UV ainda pode ser uma opção, conforme mostrado em uma imagem obtida de outro experimento em que areia foi usada como substrato de fundo (Figura 10).

Figura 10: Iluminação UV para visualizar raízes sobre fundo claro.
Exemplo de um sistema radicular de trigo duro crescendo em uma caixa de rizosfera preenchida com areia de quartzo. A caixa de rizosfera é fotografada com iluminação para (A) luz UV e (B) luz fluorescente (diurna). Barras de escala, 2 cm. Clique aqui para visualizar uma versão maior desta figura.
Medições da raiz baseadas em HSI
Figura 11 apresenta os espectros médios de três áreas de interesse nas raízes (laterais velhas e jovens, raiz principal) e duas áreas de interesse no solo (parte superior e inferior da caixa radicular).

Figura 11: Espectros médios de raiz e solo.
Espectros de regiões de interesse (ROIs) na raiz (três tipos de raiz) e no solo (parte superior e inferior da caixa de rizosfera). As ROIs são selecionadas para determinar um critério ótimo de segmentação entre raiz e solo. Clique aqui para visualizar uma versão maior desta figura.
É evidente que as raízes principais e laterais jovens diferem substancialmente do fundo quanto à intensidade na maioria das bandas espectrais. Para as laterais antigas, as diferenças de intensidade são muito menores. Uma característica que pode ser inferida visualmente é a inclinação distinta do espectro na região de absorção da água (1450 nm). Nesta região, a inclinação dos espectros das raízes é maior em comparação com os espectros do solo. Além disso, pode-se identificar uma alteração nos espectros da raiz principal e das laterais jovens na região próxima a 1100 nm, que não ocorre nas laterais antigas.
Figura 12A mostra o resultado do algoritmo de busca que identifica uma razão espectral com o maior contraste entre primeiro plano e plano de fundo. A razão entre os espectros em 1476 nm e 1076 nm fornece a melhor separação entre raízes e solo. O histograma resultante dos pixels de primeiro plano das raízes e do plano de fundo do solo é apresentado na Figura 12B. Embora haja alguma sobreposição, a maioria dos pixels está claramente separada do plano de fundo do solo. Ajustando uma curva gaussiana bimodal ao histograma e utilizando a distância de Bhattacharyya para quantificação, obtém-se um valor de 7,80. Um valor superior a 3,0 indica um alto contraste na imagem, permitindo uma separação confiável28.

Figura 12: Diferença na refletância entre o primeiro plano da raiz e o plano de fundo do solo para diferentes razões de bandas espectrais e histograma de pixels na razão espectral utilizada para segmentação.
(A) Cores claras (amarelo) indicam alto contraste entre primeiro plano e plano de fundo, cores escuras (azul) indicam baixo contraste. As primeiras 15 bandas foram removidas devido ao ruído. As linhas vermelhas indicam a razão de bandas com maior contraste. (B) Histograma de pixels de raízes (azul) e solo (vermelho) na razão espectral de segmentação. As barras azuis representam a raiz e as barras vermelhas representam o solo. O valor de intensidade corresponde à razão entre a banda espectral 160 e a banda espectral 49. Clique aqui para visualizar uma versão maior desta figura.
A imagem binária (Figura 13) é criada aplicando-se um limiar global de intensidade da razão espectral identificada em um valor de 1,008 calculado a partir da distância do histograma27. A análise do comprimento da raiz desta imagem fornece um comprimento total de 1557,3 cm, o que representa um erro de apenas 1,5% em comparação com a imagem de referência rastreada manualmente.

Figura 13: Imagem binária do sistema radicular da cultivar de beterraba açucareira Ferrara.
A imagem foi obtida mediante a aplicação de um limiar espectral global. Barra de escala (canto inferior esquerdo), 2 cm. Clique aqui para visualizar uma versão maior desta figura.
Embora a segmentação da raiz tenha melhorado com o uso de informações espectrais em comparação com informações baseadas em cor, a principal intenção da imagem hiperspectral é a análise de propriedades quimiométricas da imagem. Isso é exemplificado por meio do mapeamento do teor de água de uma imagem de caixa de raiz.
Figura 14A mostra os espectros médios dos compartimentos na rizocaixa de calibração (cf. Figura 7) preenchida com solo de diferentes conteúdos de água. A forma dos espectros é semelhante entre os compartimentos, ou seja, aqui uma razão espectral não fornece necessariamente um critério de classificação mais estável. Assim, a intensidade em uma única banda espectral (1680 nm), onde a diferença média entre conteúdos de água adjacentes é maximizada, é identificada como o melhor critério de separação. Os histogramas de pixels resultantes para este comprimento de onda espectral são mostrados na Figura 14B.

Figura 14: Características espectrais para calibração do teor de água.
(A) Espectros médios de nove compartimentos de água da caixa-radicular de calibração com diferentes conteúdos de água; (B) Histogramas de pixels para os compartimentos de água na banda 216, onde a distância média entre compartimentos vizinhos é máxima. Clique aqui para visualizar uma versão maior desta figura.
A relação entre a intensidade média do pixel em 1680 nm e o teor de água medido é mostrada na Figura 15.

Figura 15: Relação entre reflectância espectral e teor volumétrico de água.
A figura mostra pares de dados de teor de água medido e reflectância espectral com curvas empíricas (linear e exponencial) ajustadas aos dados, excluindo os teores mais altos de água (triângulos vermelhos). Clique aqui para visualizar uma versão maior desta figura.
A diferenciação de teores de água mais elevados com base na intensidade espectral torna-se difícil. Uma regressão significativa (linear ou exponencial) com alto valor de R2 pode ser ajustada para teores de água até cerca de 0,30 cm3 cm-3. Condições de solo mais úmido não podem ser previstas com confiabilidade pelo valor de intensidade. Comportamento semelhante, com uma relação exponencial entre refletividade e teor de água e com resposta decrescente para teores de água superiores a 0,30 cm3 cm-3, também foi observado em outros estudos30.
Uma imagem de rizocaixa com mapeamento detalhado do teor de água é mostrada na Figura 16. Quatro aspectos devem ser destacados. Primeiro, uma região com menor teor de água pode ser observada nas partes com raízes da rizocaixa. Segundo, a maior depleção está concentrada nas proximidades de eixos radiculares individuais. Terceiro, zonas de depleção também ocorrem onde não há eixos radiculares visíveis na superfície, indicando regiões onde as raízes estão ocultas no solo. Quarto, o mapeamento de água sem processamento adicional de imagem resulta em uma aparência irregular devido à agregação do solo. Isso pode indicar uma distribuição heterogênea do teor de água na escala do agregado, mas também pode refletir o efeito da morfologia da superfície na qualidade da imagem. Técnicas quimiométricas de processamento de imagens são uma opção para superar esses efeitos morfológicos em imagens espectrais31, mas ainda não foram implementadas nos scripts Matlab utilizados aqui.

Figura 16: Mapeamento do teor de água em uma caixa radicular.
As cores azul-escuras representam regiões com alto teor de água, enquanto as áreas de verde a vermelho indicam regiões com baixo teor de água. A raiz da planta é sobreposta à imagem em preto. Barra de escala (canto inferior esquerdo), 2 cm. Clique aqui para visualizar uma versão maior desta figura.
Exemplos de aplicação
Figura 17 relaciona características quantitativas da raiz obtidas por análise de imagem com medições da parte aérea.

Figura 17: Exemplos típicos de aplicações para dados radiculares.
(A) e (B) mostram informações sobre raízes utilizadas para a caracterização da parte aérea e subterrânea da planta em um contexto de fenotipagem. (A) representa o crescimento radicular do cultivar de beterraba sacarina Ferrara, (B) compara seis cultivares de beterraba sacarina cultivados em caixas de rizosfera utilizando a razão entre área foliar e radicular (dados de uma repetição). (C) e (D) são relações funcionais entre características conforme observadas em pesquisas fisiológicas de plantas. (C) mostra a influência da razão entre área foliar e radicular na produção de matéria seca e (D) a relação entre área superficial da raiz e condutância estomática. Clique aqui para visualizar uma versão ampliada desta figura.
Figura 17A e 17B são relevantes para fenotipagem com foco na caracterização abrangente de partes aéreas e subterrâneas da planta. Figura 17A mostra o crescimento radicular da cultivar de beterraba sacarina Ferrara (cf. Figura 8 para imagens). A expansão do sistema radicular indica a capacidade de uma cultivar de explorar o volume do solo em um determinado período do ciclo de vegetação. Figura 17B mostra a razão entre a área foliar e a área radicular de seis cultivares de beterraba sacarina, fornecendo um descritor para o equilíbrio entre o suprimento (raiz) e a demanda (folha) da planta.
Figuras 17C e 17D apresentam exemplos de relações funcionais de interesse na pesquisa fisiológica. Na Figura 17C, a relação entre área foliar e radicular é relacionada à matéria seca formada durante o experimento, indicando o papel predominante da superfície assimiladora como fator limitante para o acúmulo de matéria seca. A falta de significância, apesar de um valor comparativamente alto de R2, está relacionada ao baixo número de pares de dados (n=6) utilizados aqui. A Figura 17D revela que cultivares com maior área radicular (absorção aprimorada) apresentam condutância estomática média mais elevada ao longo do experimento. A área radicular maior aparentemente sustenta a extração de água, prolongando assim a abertura dos estômatos.