Павел А. Дмитриев *, Анастасия А. Дмитриева, Борис Л. Козловский
Ботанический сад, Академия биологии и медицины, Южный федеральный университет, Ростов-на-Дону 344006, Россия; admit@sfedu.ru (А.А.Д.); blk@sfedu.ru (Б.Л.К.)
- Корреспонденция: pdmitriev@sfedu.ru
АннотацияПодсолнечник является одной из важнейших сельскохозяйственных культур в мире. Учитывая высокую значимость продуктов подсолнечника на мировом рынке и масштабы его возделывания, внедрение технологий точного земледелия в его культуру может дать значительный экономический и экологический эффект. В данном исследовании продемонстрирована принципиальная возможность разработки технологии быстрой, дистанционной и неинвазивной оценки влажности семян подсолнечника для определения оптимальных сроков десикации и уборки урожая. Показано, что влажность семян подсолнечника можно оценивать с высокой точностью на основе спектральных характеристик нижней стороны соцветий, полученных с помощью гиперспектральной камеры в видимом и ближнем инфракрасном диапазоне (VNIR) (от 450 до 950 нм). Для прогнозирования влажности семян подсолнечника использовалась регрессия случайного леса (RFR). Модель показала отличные результаты на обучающей выборке (R²_c = 1,00; MAE_c = 0,58; RMSE_c = 0,74; MAPE_c = 1,29) и высокую производительность на тестовой выборке (R²_t = 0,98; MAE_t = 2,99; RMSE_t = 3,28; MAPE_t = 12,22). Наиболее значимыми вегетационными индексами для определения влажности являются CCI, Booch, Datt3, Datt4, LSIRed, modPRI, SR5, TCARI и TCARI2.
Ключевые слова: Helianthus annuus; десикация; уборка урожая; гиперспектральная съёмка
Академический редактор: Хосе Антонио Эрнандес Кортес
Получено: 29 августа 2025 г.; Исправлено: 24 октября 2025 г.; Принято: 27 октября 2025 г.; Опубликовано: 29 октября 2025 г.
Цитирование: Дмитриев П.А., Дмитриева А.А., Козловский Б.Л. Оценка влажности семян подсолнечника по спектральным характеристикам соцветий в VNIR. Seeds 2025, 4, 55.
https://doi.org/10.3390/seeds4040055Авторское право: © 2025 авторы. Лицензиат MDPI, Базель, Швейцария. Эта статья распространяется на условиях открытого доступа в соответствии с лицензией Creative Commons Attribution (CC BY) (
https://creativecommons.org/licenses/by/4.0/).
1. ВведениеПодсолнечник (
Helianthus annuus L.), являясь экономически важной культурой, занимает четвертое место среди масличных культур в мире. По данным Министерства сельского хозяйства США (USDA), мировая площадь посевов подсолнечника в 2023–2024 годах составила более 27,3 млн гектаров, а общий урожай — 55,97 млн тонн (
https://www.fas.usda.gov/data/production/commodity/2224000, дата обращения 24 октября 2025 г.). Объём рынка подсолнечника в 2024 году составил 23,01 млрд долл. США. Ожидается, что к 2029 году он вырастет до 29,56 млрд долл. США при среднегодовом темпе роста 5,1% (
https://www.thebusinessresearchcompany.com/report/sunflower-global-market-report, дата обращения 14 августа 2025 г.).
Учитывая высокую значимость продуктов подсолнечника на мировом продовольственном рынке и масштабы его возделывания, внедрение технологий точного земледелия в его культуру может дать значительный экономический и экологический эффект [1–3]. Поэтому в настоящее время для определения сроков проведения агротехнических мероприятий востребованы точные, быстрые, крупномасштабные технологии, не требующие экспертного контроля и обеспечивающие свободный доступ к результатам [4]. Определение сроков десикации и уборки урожая являются важнейшими агротехническими мероприятиями для подсолнечника. Десикация сокращает сроки созревания подсолнечника и его уязвимость для птиц, помогая избежать проблем с грибковыми заболеваниями. Своевременная уборка подсолнечника помогает снизить потери при обмолоте и сэкономить на сушке семян. Учитывая масштабы возделывания этой культуры, точность и оперативное прогнозирование оптимальной даты десикации и уборки для каждого поля могут быть использованы в управлении логистикой, планировании операций с техникой, а также для прогнозирования урожайности. Влажность семян является предиктором для определения физиологической зрелости семян и, соответственно, даты десикации и уборки [5].
Применение десикантов (дибромид, глюфосинат-аммоний, карфентразон-этил) для подсолнечника рекомендуется при влажности семян 25–30%, что соответствует физиологической зрелости [6]. Оптимальная влажность семян для уборки составляет 10–12% при условии последующей сушки перед хранением до 8% и ниже [7–9].
Гиперспектральная съёмка (HSI) представляется перспективной для разработки технологий оперативной, дистанционной, неинвазивной, высокотехнологичной и надёжной оценки влажности семян. Прогресс в области гиперспектрального дистанционного зондирования (RS) стал предвестником революции в мониторинге сельского хозяйства [4,10]. Гиперспектральная съёмка может предоставить бесценные данные для высокопроизводительного фенотипирования культур [11], включая диагностику болезней [12], диагностику стресса растений [11], мониторинг и оценку урожайности [13], оценку фенологического состояния [14] и классификацию культур [15]. Использование методов машинного и глубокого обучения (ML) позволяет связать извлечённую гиперспектральную информацию и полученные фенотипы [11,16]. Согласно недавнему систематическому обзору по точному земледелию [17] применительно к подсолнечнику и роли HSI, исследования были сосредоточены на управлении посевами, мониторинге развития подсолнечника (определение стадий роста и вегетации), борьбе с сорняками и промышленном применении.
Анализ предыдущих результатов по влажности семян, исследованиям содержания воды в растениях и исследованиям засухоустойчивости подсолнечника и других сельскохозяйственных культур важен для достижения поставленной цели. Huang и др. [18] использовали гиперспектральные технологии и спектральные характеристики семян подсолнечника в диапазоне длин волн 384–1034 нм для прогнозирования жизнеспособности и влажности семян. Модель SNV-XGBoost-LightGBM продемонстрировала высокую точность определения влажности с коэффициентом детерминации (R²) 0,97 и среднеквадратичной ошибкой (RMSE) 0,84. Имеется ряд исследований других сельскохозяйственных культур, в которых лабораторная проксимальная HSI семян использовалась для определения их влажности с достаточной точностью [19–22].
В то же время только дистанционное зондирование может обеспечить оперативность и масштаб диагностики влажности семян сельскохозяйственных культур для определения сроков проведения агротехнических мероприятий. В связи с этим проблема многих культур (подсолнечник, кукуруза, бобовые) заключается в недоступности семян для съёмки, которая проводится с помощью беспилотных летательных аппаратов (БПЛА) или вышек. Поэтому большой интерес представляют исследования, в которых влажность семян определяется на основе взаимосвязи между семенами и спектральными характеристиками других органов растений (листьев, соцветий, плодов), которые могут быть дистанционно обнаружены в полях. Интересное исследование провели Litvishchenko и др. [23], в котором влажность семян подсолнечника определялась по спектральным характеристикам нижней стороны соцветия в микроволновом диапазоне. Xu и др. [24] прогнозировали сроки уборки кукурузы на основе взаимосвязи между содержанием хлорофилла в пологе и влажностью семян.
Исследование Steidle Neto и др. [25], в котором была установлена корреляция между содержанием хлорофилла и воды в пологе подсолнечника, указывает на возможность определения влажности семян подсолнечника в полевых условиях. Этот подход осуществим, поскольку определение фотосинтетических пигментов по спектральным характеристикам листьев и полога подсолнечника методически разработано [26,27]. Важно отметить, что недавний обзор методов определения содержания воды в растительности [28] указывает на то, что нормализованный разностный вегетационный индекс (NDVI), чувствительный к хлорофиллу, является более предпочтительным индикатором содержания воды, чем другие вегетационные индексы (VIs). Большинство VIs являются математическими функциями отражения в видимой и ближней инфракрасной областях и могут эффективно обнаруживать фотосинтетическую способность растительного покрова. В то же время разработаны специализированные VIs для определения содержания влаги, среди которых наиболее часто используются нормализованный разностный инфракрасный индекс (NDII), нормализованный разностный водный индекс (NDWI) и водный полосный индекс (WBI) [28]. Эти формулы используют спектральные диапазоны в ближней инфракрасной (NIR) и коротковолновой инфракрасной (SWIR) областях. Сравнение эффективности определения влажности растительности с использованием чувствительных к хлорофиллу VIs со специализированными «водными» VIs дало противоречивые результаты. По результатам Zhang и Zhou [29], среди десяти узкополосных VIs (NDVI, Red edge normalized ratio (NRred edge), Green chlorophyll index (CIgreen), Red edge chlorophyll index (CIred edge), Water index (WI), Moisture stress index (MSI), Simple ratio water index (SRWI), Normalized difference water index (NDWI), Normalized difference water index centered at 1640 nm (NDWI1640), Normalized difference water index centered at 2130 nm (NDWI2130)) CIgreen, NRred edge и CIred edge были наиболеечувствительны к изменениям влажности растительности. Zhou и др. [30] обнаружили, что WI и NDWI лучше отражают влажность растительности, чем NDII, NDVI и OSAVI.
Аналогичные результаты были получены при диагностике содержания воды с использованием спектральных диапазонов. В ряде исследований [19,31–34] успешно диагностировали содержание воды в растительности с использованием спектральных диапазонов в видимой и ближней инфракрасной областях (от 400 до 1000 нм).
Другие исследования показали, что высокая точность определения влажности растительности и семян может быть достигнута с использованием дополнительных спектральных диапазонов в SWIR-области [20,35].
Таким образом, в процессе созревания семена подсолнечника не могут быть сканированы спектральными сенсорами с БПЛА. Это ограничивает возможность быстрого и дистанционного определения их влажности. С другой стороны, ряд исследований показал тесную взаимосвязь между влажностью семян сельскохозяйственных культур и спектральными характеристиками их вегетативных органов.
На основе вышеуказанных пробелов мы выдвинули гипотезу о том, что влажность семян подсолнечника можно определять по спектральным характеристикам нижней стороны их соцветия. Для подтверждения гипотезы был проведён поисковый эксперимент на сорте «Джинн М». Целью эксперимента было доказать принципиальную возможность определения влажности семян подсолнечника на основе спектральных характеристик нижней стороны соцветия в VNIR-диапазоне. В ходе исследования использовалась лабораторная проксимальная HSI — метод гиперспектрального фенотипирования, который позволяет минимизировать действие случайных факторов и лучше подходит для исследовательских экспериментов, требующих принципиальной оценки гипотезы [36]. В результате были получены подробные временные ряды (TS) спектральных характеристик нижней стороны соцветия и соответствующие им временные ряды влажности семян. Было рассчитано 86 VIs, чувствительных к содержанию пигментов, и использовано для построения модели регрессии случайного леса (RFR).
Аббревиатуры, введённые в рукописи, и их определения перечислены в таблице 1.
Таблица 1. Список аббревиатур. Рисунок 1. Внешний вид поля подсолнечника (а). Верхняя и нижняя стороны соцветия подсолнечника (b).
2.2. Климатические и почвенные условия исследуемого районаСальский район расположен на юго-западе Ростовской области, которая имеет умеренно-континентальный климат.
Средняя сумма температур воздуха за многолетний период выше 10°C составляет 3304°C, среднегодовая температура воздуха — 8,9°C, максимальная температура в июне — от 22,9 до 23,8°C, минимальная в январе — от –4,7 до –5,5°C. Радиационный баланс составляет от 2641 до 2685 МДж/м² в год. Средняя сумма осадков за многолетний период составляет от 474 до 500 мм, в том числе за вегетационный период — от 290 до 300 мм; испаряемость за год — от 825 до 912 мм (
https://weatherspark.com/y/102167/Average-Weather-in-Sal%E2%80%99sk-Russia-Year-Round, дата обращения 14 августа 2025 г.).
Динамика среднесуточной температуры воздуха и количества осадков в период вегетации подсолнечника на экспериментальном поле представлена на рисунке 2.
Рисунок 2. Среднесуточная температура воздуха и количество осадков за период вегетации подсолнечника сорта «Джинн М» в 2024 г. (14 апреля 2024 г. – 7 сентября 2024 г.).
Почва — чернозём обыкновенный, очень тёплый и переходно-мерзлотный. Гранулометрический состав почвы в районе характеризуется как глинистый. Мощность гумусового горизонта составляет от 75 до 140 см, содержание гумуса — от 3,7 до 4,2%.
2.3. Методы отбора образцов
Для лабораторной HSI соцветия подсолнечника отбирались в период от развития плодов (BBCH 65) до конца созревания (BBCH 92). Для достижения максимального диапазона влажности семян и равномерного представления соцветий подсолнечника с различной влажностью семян на каждую дату эксперимента проводился стратифицированный отбор соцветий (рисунок 2).
Отбор соцветий в поле проводился с приблизительным определением их влажности путём визуальной оценки цвета нижней стороны соцветия подсолнечника. Следует отметить, что подсолнечник сильно различался по степени зрелости семян и их влажности соответственно (это хорошо видно на рисунке 3). Хотя такой отбор соцветий идеален для создания модели влажности семян на основе спектральных характеристик, он не является репрезентативным для поля, на котором проводилось исследование.
Рисунок 3. Соцветия подсолнечника (цифры обозначают порядковые номера) и соответствующая влажность семян (%). Влажность семян измерялась с помощью анализатора влажности Sartorius MA-35. Дата эксперимента: 5 августа 2024 г.
2.4. Метод гиперспектральной съёмки
Время от отбора образцов до доставки в лабораторию составляло три часа. Гиперспектральная съёмка нижней стороны соцветий подсолнечника проводилась при искусственном освещении с использованием гиперспектральной камеры Cubert UHD-185 (Cubert GmbH, Ульм, Германия) (рисунок 4). Основные характеристики гиперспектральной камеры приведены в таблице 2.
Рисунок 4. Гиперспектральная съёмка нижней стороны соцветия подсолнечника при искусственном освещении в лабораторных условиях.
Таблица 2. Ключевые характеристики гиперспектральной камеры Cubert UHD-185.
Соцветия помещались на чёрную пластину на расстоянии 40 см от объектива камеры. Размер пикселя на этом расстоянии составляет 0,25 см². Среднее количество пикселей на одно соцветие составило 1370.
2.5. Предварительная обработка гиперспектральных данных
Спектральные данные были предварительно обработаны с использованием фильтра Савицкого–Голая. Затем спектральные профили по всему гиперспектральному изображению были преобразованы в синтетические спектральные профили с использованием метода Random Reflectance [37] для снижения шума и повышения точности RFR. Для выбора области интереса (ROI) в гиперспектральном изображении использовался двухэтапный метод сегментации. Сначала ROI выбиралась по пороговому значению среднего Carter5 > 1,4, как описано в Dmitriev и др. [38]. Затем применялась морфологическая эрозия со структурирующим элементом 3×3.
2.6. Расчёт вегетационных индексов
Вегетационные индексы рассчитывались в среде R (R Core Team, версия 4.5.1) на основе значений спектральных каналов, полученных с помощью Cubert UHD-185 [39].
2.7. Определение влажности семян подсолнечника
Семена подсолнечника созревают неравномерно в одном соцветии — созревание идёт от края соцветия к его центру. Для баланса между пропускной способностью анализа и репрезентативностью из каждого соцветия отбиралась проба массой 20 г, поскольку предварительные тесты показали вариацию влажности менее 2% по трём подпробам. Определение влажности семян проводилось с использованием лабораторного анализатора влажности Sartorius MA-35 (Sartorius AG, Гёттинген, Германия). Значения влажности выражались в процентах.
2.8. Анализ данных
Исследовательский анализ данных (EDA) проводился с использованием стандартных статистических методов и анализа главных компонент (PCA). Мультиколлинеарность VIs оценивалась с использованием коэффициента корреляции Спирмена (r). Сила связи оценивалась по шкале Чеддока. Мультиколлинеарность переменных отрицательно влияет на качество моделей машинного обучения. Она приводит к переобучению модели, увеличению сложности модели и времени обучения, усложнению интерпретации параметров и снижению надёжности оценок параметров [40]. Среди алгоритмов машинного обучения RFR показал устойчивость к высокоразмерным мультиколлинеарным спектральным данным с минимальным риском переобучения [41]. Поэтому для прогнозирования влажности семян подсолнечника на основе спектральных характеристик нижней стороны соцветия использовался алгоритм RFR. Модель RFR обучалась на 70% отобранных соцветий и тестировалась на 30%. Для настройки гиперпараметров и оценки эффективности использовался метод 5-кратной перекрёстной проверки. Количество деревьев: 500; количество переменных, пробуемых на каждом разбиении: 3. Корректность и производительность модели оценивались с использованием средней абсолютной процентной ошибки (MAPE), среднеквадратичной ошибки (RMSE), средней абсолютной ошибки (MAE) и коэффициента детерминации (R²). Для обучающей выборки метрики обозначены как MAPE_c, RMSE_c, MAE_c и R²_c, а для тестовой — MAPE_t, RMSE_t, MAE_t и R²_t.