Представлены примеры результатов сегментации корней на основе RGB- и HS-визуализации. Для спектральной визуализации приведен пример картирования воды с высоким разрешением. В завершение показаны результаты, демонстрирующие научный контекст применения данных о корнях, полученных с помощью анализа изображений.
Измерение корней на основе RGB-изображений
Рисунок 8 представляет собой временную серию RGB-изображений корня сахарной свеклы сорта Ferrara. На изображениях видны артефакты, вызванные неоднородным освещением ризобоксов: более светлые области вдоль левой стороны и разная яркость в зоне перекрытия верхнего и нижнего изображений.

Рисунок 8: Временные ряды роста корней, полученные с помощью RGB-визуализации.
На снимках представлена сахарная свекла сорта Ferrara в разные дни после посева (DAS). На изображениях видны некоторые артефакты, вызванные неоднородным освещением в левой части снимка, а также между верхним и нижним изображениями. Масштабные отрезки — 2 cm. Пожалуйста, нажмите здесь, чтобы просмотреть этот рисунок в большем размере.
Рисунок 9 содержит подробную информацию о сегментации корней на основе цветового порога для сорта Ferrara на 35-й день после посева (DAS). В качестве эталона (Рисунок 9A) используется бинарное изображение, на котором все корни были отслежены вручную с помощью графического планшета. Время, затраченное на ручное отслеживание всей полностью развитой густой корневой системы сахарной свеклы, составило около четырех часов. На Рисунке 9B представлен детальный вид выбранного участка в верхней части изображения, где присутствуют старые боковые корни. Здесь несколько осевых корней не классифицированы с помощью цветового порога. Напротив, в нижней части (Рисунок 9C), где преобладают белые молодые корни, сегментация по цвету правильно классифицирует все осевые корни. Бинаризированная корневая система (Рисунок 9D) имеет черную область с левой стороны из-за артефакта освещения, которая была определена как область исключения перед проведением количественного анализа. На Рисунке 9E показаны соответствующие пиксельные гистограммы выбранных признаков (корни против почвы) для красного канала RGB-изображения сорта Ferrara на 35-й DAS. Пиксели корней (синий цвет) четко демонстрируют три пика, соответствующих светлым молодым боковым корням, темным старым боковым корням и стержневому корню. Перекрытие между старыми боковыми корнями и фоном почвы очень сильное, что приводит к появлению неклассифицированных осевых корней (ср. Рисунок 9B).

Рисунок 9: Сегментация корней с использованием цветового порога.
(A) Корневая система, сегментированная вручную с помощью графического планшета; (B) область с плохо сегментированными старыми осями корня в верхней части изображения и (C) правильно сегментированные молодые оси в нижней части изображения; (D) бинарное изображение, полученное путем пороговой обработки по цвету. (E) Гистограммы пикселей для выбранных признаков RGB-изображения. Корни представлены синими столбцами с указанием различных типов корней; почва представлена красными столбцами. Масштабные отрезки на A и D — 2 cm; масштабные отрезки на B и C — 1 cm. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Итоговая общая длина видимых корней, определенная для эталонного изображения при ручной сегментации, составляет 1534.1 cm, тогда как автоматизированная сегментация по цвету дает общую длину корней 1427.6 cm.
Изображения в градациях серого, полученные при УФ-освещении, не дали преимущества в данном случае и показали худшие результаты по сравнению с цветовым пороговым анализом (длина корня: 1679,7 cm). Старые корни не удалось сегментировать, а на изображении наблюдалось больше шума, что, вероятно, связано с более низкой интенсивностью света УФ-ламп. Тем не менее, в случае молодых корней с высокой автофлуоресценцией и светлым фоном субстрата УФ-освещение все еще может быть подходящим вариантом, что демонстрирует изображение из другого эксперимента, где в качестве фонового субстрата использовался песок (Рисунок 10).

Рисунок 10: УФ-освещение для визуализации корней на светлом фоне.
Пример корневой системы твердой пшеницы, растущей в ризобоксе, заполненном кварцевым песком. Ризобокс сфотографирован при освещении (A) УФ-светом и (B) флуоресцентным (дневным) светом. Масштабные отрезки — 2 cm. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Измерения корней на основе HSI
Рисунок 11 представляет средние спектры для трех областей интереса (ROI) корня (старый и молодой боковой корень, стержневой корень) и двух областей интереса почвы (верх и низ ризобокса).

Рисунок 11: Средние спектры корня и почвы.
Спектры из областей интереса (ROI) на корне (три типа корней) и в почве (верхняя и нижняя части ризобокса). Области ROI выбраны для определения оптимального критерия сегментации между корнем и почвой. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Очевидно, что главный и молодые боковые корни существенно отличаются от фона по интенсивности большинства спектральных зон. Для старых боковых корней разница в интенсивности гораздо ниже. Характерной особенностью, которую можно определить визуально, является различный наклон спектра в области поглощения воды (1450 nm). Здесь наклон спектров корня выше по сравнению со спектрами почвы. Кроме того, в области около 1100 nm можно обнаружить изменение спектров главного и молодых боковых корней, которое отсутствует у старых боковых корней.
Рисунок 12A демонстрирует результат работы алгоритма поиска, определяющего спектральное отношение с наиболее сильным контрастом между объектом и фоном. Отношение спектров на 1476 nm к 1076 nm обеспечивает наилучшее разделение корней и почвы. Результирующая гистограмма пикселей переднего плана (корней) и фона (почвы) представлена на Рисунке 12B>. Несмотря на некоторое перекрытие, большинство пикселей четко отделены от почвенного фона. При аппроксимации гистограммы бимодальной кривой Гаусса и использовании расстояния Бхаттачарйи для количественной оценки получено значение 7,80. Значение выше 3,0 указывает на сильный контраст изображения, что позволяет проводить надежное разделение28.

Рисунок 12: Разница в отражательной способности между передним планом (корнями) и фоном (почвой) для различных соотношений спектральных зон и гистограмма пикселей при спектральном соотношении, используемом для сегментации.
(A) Яркие цвета (желтый) указывают на высокий контраст между передним планом и фоном, темные цвета (синий) — на низкий контраст. Первые 15 зон были удалены из-за шума. Красные линии обозначают соотношение зон с наивысшим контрастом. (B) Гистограмма пикселей корней (синий) и почвы (красный) при спектральном соотношении сегментации. Синие столбцы представляют корни, красные — почву. Значение интенсивности соответствует отношению спектральной зоны 160 к спектральной зоне 49. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Бинарное изображение (Рисунок 13) создается путем применения глобального порога интенсивности определенного спектрального отношения со значением 1.008, рассчитанным на основе расстояния по гистограмме27. Анализ длины корня на этом изображении дает общую длину 1557.3 cm, что соответствует погрешности всего 1.5% по сравнению с эталонным изображением, обработанным вручную.

Рисунок 13: Бинарное изображение корневой системы сахарной свеклы сорта Ferrara.
Изображение получено путем применения глобального спектрального порога. Масштабная линейка (нижний левый угол), 2 cm. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Хотя сегментация корней при использовании спектральной информации улучшилась по сравнению с использованием цветовой информации, основной целью гиперспектральной (HS) визуализации является анализ хемометрических свойств изображения. Это продемонстрировано на примере картирования содержания воды на изображении ризобокса.
Рисунок 14A демонстрирует средние спектры отсеков в калибровочном ризобоксе (см. Рисунок 7), заполненном почвой с различным содержанием воды. Форма спектров в разных отсеках схожа, то есть в данном случае спектральное отношение не обязательно обеспечивает более стабильный критерий классификации. Таким образом, в качестве наилучшего разделяющего критерия определена интенсивность на одной спектральной полосе (1680 nm), где средняя разница между соседними значениями содержания воды максимальна. Полученные гистограммы пикселей для данной длины волны показаны на Рисунке 14B.

Рисунок 14: Спектральные характеристики для калибровки содержания воды.
(A) Средние спектры девяти водных отсеков из калибровочного ризобокса с разным содержанием воды; (B) Пиксельные гистограммы для водных отсеков на полосе 216, где среднее расстояние между соседними отсеками максимально. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Зависимость средней интенсивности пикселей на длине волны 1680 nm от измеренного содержания воды показана на рисунке 15.

Рисунок 15: Зависимость спектрального коэффициента отражения от объемного содержания воды.
На рисунке представлены пары данных измеренного содержания воды и спектрального коэффициента отражения с эмпирическими кривыми (линейной и экспоненциальной), построенными по данным за исключением значений наибольшего содержания воды (красные треугольники). Пожалуйста, нажмите здесь, чтобы просмотреть этот рисунок в большем размере.
Дифференциация более высокого содержания воды по спектральной интенсивности становится затруднительной. Значимая регрессия (линейная или экспоненциальная) с высоким R2 может быть аппроксимирована для содержания воды примерно до 0.30 cm3 cm-3. Более влажные условия почвы не могут быть надежно предсказаны по значению интенсивности. Аналогичное поведение экспоненциальной зависимости между отражательной способностью и содержанием воды с уменьшением отклика при содержании воды выше 0.30 cm3 cm-3 также было обнаружено в других исследованиях30.
Изображение ризобокса с детальным картированием содержания воды показано на рисунке 16. Следует отметить четыре аспекта. Во-первых, в укорененных частях ризобокса видна область с более низким содержанием воды. Во-вторых, наиболее сильное истощение сосредоточено вблизи отдельных осей корней. В-третьих, зоны истощения также наблюдаются там, где оси корней не видны на поверхности, что указывает на области, где корни скрыты в почве. В-четвертых, картирование воды без дополнительной обработки изображений приводит к пятнистому виду из-за агрегатов почвы. Это может указывать как на неоднородное распределение содержания воды на уровне агрегатов, так и на влияние морфологии поверхности на качество изображения. Хемометрические методы обработки изображений позволяют преодолеть такие морфологические эффекты в спектральных изображениях31, однако они пока не реализованы в используемых здесь скриптах Matlab.

Рисунок 16: Карта содержания воды в ризобоксе.
Темно-синий цвет соответствует областям с высоким содержанием воды, области от зеленого до красного цвета показывают регионы с низким содержанием воды. Корень растения наложен на изображение черным цветом. Масштабная линейка (в нижнем левом углу) — 2 cm. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Примеры применения
Рисунок 17 сопоставляет количественные признаки корней, полученные в результате анализа изображений, с наземными измерениями.

Рисунок 17: Типичные примеры применения данных о корневой системе.
(A) и (B) демонстрируют информацию о корнях, используемую для характеризации надземных и подземных частей растения в контексте фенотипирования. (A) представляет рост корней сорта сахарной свеклы Ferrara, (B) сравнивает шесть сортов сахарной свеклы, выращенных в ризобоксах, с использованием соотношения площади поверхности листьев и корней (данные одного повтора). (C) и (D) представляют функциональные зависимости между признаками, выявленные в ходе физиологических исследований растений. (C) показывает влияние соотношения площади поверхности листьев и корней на производство сухого вещества, а (D) — зависимость площади поверхности корня от устьичной проводимости. Щелкните здесь, чтобы просмотреть увеличенную версию этого рисунка.
Рисунок 17A и 17B относятся к фенотипированию, направленному на комплексную характеристику надземной и подземной частей растений. На рисунке 17A показан рост корней сахарной свеклы сорта Ferrara (см. рисунок 8 для просмотра изображений). Расширение корневой системы указывает на способность сорта осваивать объем почвы за определенный промежуток времени вегетационного периода. На рисунке 17B показано соотношение площади поверхности листьев к площади поверхности корней у шести сортов сахарной свеклы, что служит дескриптором баланса между обеспечением растения (корень) и его потребностями (лист).
Рисунки 17C и 17D приводят примеры функциональных зависимостей, представляющих интерес в физиологических исследованиях. На рисунке 17C отношение площади поверхности листьев к площади поверхности корней соотносится с количеством сухого вещества, образовавшегося в ходе эксперимента, что указывает на определяющую роль ассимиляционной поверхности как лимитирующего фактора накопления сухого вещества. Отсутствие статистической значимости, несмотря на сравнительно высокий R2, связано с малым количеством парных данных (n=6), использованных в данном случае. Рисунок 17D показывает, что сорта с большей площадью поверхности корней (улучшенным поглощением) имеют в среднем более высокую устьичную проводимость на протяжении всего эксперимента. Большая площадь корней, по всей видимости, поддерживает экстракцию воды, тем самым продлевая открытие устьиц.