способ идентификации типов растительности
Классы МПК: | G01C11/00 Фотограмметрия или видеограмметрия, например стереограмметрия; топографическая съемка местности с помощью фотографирования G01C11/04 расшифровка изображений G01C11/10 с использованием вычислительных устройств для контроля положения изображения G03B37/00 Панорамная фотосъемка и съемка широкоэкранных кинофильмов; фотографирование обширных поверхностей, например для топографической съемки; фотографирование полостей, например внутренней поверхности труб A01C23/00 Устройства, специально предназначенные для распределения навозной жижи или других жидких удобрений, в том числе аммиачной воды, например резервуары для перевозки, повозки для разбрызгивания G06F17/14 преобразования Фурье, Уолша или аналогичные преобразования G06F17/15 вычисление корреляционных функций G06F17/16 матричные или векторные вычисления G06F17/17 вычисление функций приближенными методами, например интерполяцией или экстраполяцией, сглаживанием, методом наименьших квадратов |
Автор(ы): | Давыдов В.Ф., Григорьева О.Ю., Щербаков А.А., Васильев Н.И. |
Патентообладатель(и): | Московский государственный университет леса |
Приоритеты: |
подача заявки:
1994-09-27 публикация патента:
20.07.1998 |
Формула изобретения
Способ идентификации типов растительности путем определения ее состояния и тематической оценки на контрольных площадках, дистанционного зондирования подстилающей поверхности, содержащей контрольные площадки, с орбиты космического аппарата средствами, установленными на его автономной платформе, отслеживания заданных районов наблюдения вводом уставок баллистических данных в систему управления автономной платформы, отличающийся тем, что осуществляют дистанционную, покадровую видеосъемку объектов, передают получаемые изображения и регистрируют их на наземных средствах, разбивают видеокадр на последовательность единичных участков и преобразуют изображение каждого участка в цифровую функцию интенсивности, вычисляют двумерный пространственный спектр функции интенсивности и физических размеров элементарного участка, получают из пространственного спектра энергетический спектр сигнала изображения, вычисляют обратным Фурье-преобразованием автокорреляционную функцию сигнала изображения, рассчитывают эталонную регрессионную зависимость между параметрами автокоррреляционных функций для различных типов растительности контрольных площадок, синтезируют из последовательно проанализированных единичных участков мозаичную картину растительности на отснятых видеокадрах по всей площади наблюдения.Описание изобретения к патенту
Предлагаемое изобретение относится к лесному хозяйству, в частности к оперативной диагностике растительных покровов на вырубках, оценке текущих и последующих стадий возобновления леса. При тематическом картировании групп лесов и условий произрастания выделяют следующие типы растительности: лиственничное мелколесье, лиственничные леса высоко-, средне- и низкобонитетные, сосновые леса, гари, аласы (кустарниковые, травянистые). В лесоводстве традиционный метод установления типа вырубки и стадии возобновления леса заключается в их визуально-инструментальном обследовании и картировании. За минимальную величину выдела принимается 1 га, масштаб 1:10000, или линейные размеры элементарной площадки 100 х 100 м. Известна оценка состояния вырубки путем стационарного наблюдения учебно-опытного участка на интервале 20 лет (см., например, И.С. Мелехов. Лесоводство. М.: Агропромиздат, 1989, с.8, табл. 2 - аналог). При известном способе формирование древостоя на вырубках на всех стадиях - густого молодняка, крупного жердняка, приспевающего и спелого леса оценивают количеством деревьев на 1 га и их толщиной более 6 см на высоте 1,3 м. Известны также различные методы (ИЮФРО, по Крафту, Жилкина Б.Д.) классификации деревьев по качеству кроны, диаметру стволов, жизненности, положению в древостое (см. там же, аналог, с.245 - 247). Недостатками известных аналогов являются большая трудоемкость, связанная с необходимостью индивидуального пересчета каждого дерева, неоперативность, недоступность горных и отдаленных районов. Ближайшим аналогом по технической сущности к заявляемому способу является способ, представленный в кн. Исаева А.С. и др. Аэрокосмический мониторинг лесов. М.: Наука, 1991, с. 28 - 30. В ближайшем аналоге оценку (в условных баллах) категорий состояния лесов проводят дистанционным методом, путем измерения коэффициентов спектральной яркости зондируемых площадок в синем, зеленом и красном участках видимого диапазона, расчета значений хроматических коэффициентов и сопоставления измеренных значений с тарированными значениями функции регрессии хроматических коэффициентов, полученной зондированием эталонных (пробных) площадок. Известный способ реализуется следующей последовательностью операций:выбирают контрольные площадки и определяют визуально-инструментальным методом состояние древостоя в баллах;
размещают спектрометрический модуль на автономной поворотной платформе искусственного спутника Земли;
отслеживают заложенные площадки вводом уставок баллистических данных в систему управления автономной платформы;
измеряют коэффициенты спектральной яркости (
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115048/961.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115048/961.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115048/961.gif)
квантуют непрерывные значения амплитуд сигналов измеренной функции с шагом квантования не более 1/256;
рассчитывают хроматические коэффициенты жизненности (g) и поражения (r):
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-2t.gif)
вычисляют функцию корреляции хроматических коэффициентов и тарируют ее по эталонным измерениям контрольных площадок с известными категориями состояния в баллах;
получают оценку состояния лесных массивов по всей площади наблюдения. Известный способ имеет следующие недостатки. 1. Способ не может быть непосредственно использован для идентификации типов растительности, поскольку принимаемый сигнал несет интегральную характеристику подстилающей поверхности, из которого уже нельзя отдельно выделить составляющую структуры полога. 2. Коэффициенты спектральной яркости различных типов растительности имеют один и тот же вид зависимостей и отличаются друг от друга на единицы процентов, что не обеспечивает необходимую точность и достоверность идентификации. В принципе, по индексу цветности в вегетационный период трудно отличить молодую зелень от зелени луга. Задача, решаемая данным изобретением, заключается в регистрации и преобразовании текстурных признаков изображений, получаемых дистанционными методами в "шероховатость" полога растительности, достаточную для идентификации ее типа, а также использовании преимуществ, вытекающих из космических технологий в виде повышения оперативности, производительности и глобальности получения целевой информации. Поставленная цель достигается введением в ближайший аналог следующих технологических операций:
осуществляют высокодетальную съемку подстилающей поверхности, содержащей контрольные площадки, и записывают полученные изображения на бортовой магнитофон;
передают в сеансах связи по радиоканалу в Центр управления полетом полученные изображения и регистрируют их на кассеты;
разбивают видеокадр на последовательность единичных участков и преобразуют изображение каждого участка в цифровую функцию интенсивности I (x, y);
вычисляют двумерный пространственный спектр от функции интенсивности и физических размеров элементарного участка прямым Фурье-преобразованием (на основе программных методов быстрого Фурье-преобразования) в соответствии с формулой
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-3t.gif)
получают из пространственного спектра G (fx, fy) энергетический спектр сигнала изображения S(F);
вычисляют обратным Фурье-преобразованием автокорреляционную функцию B(
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-4t.gif)
рассчитывают эталонную регрессионную зависимость между параметрами автокорреляционных функций B(
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
синтезируют из последовательно проанализированных единичных участков мозаичную картину растительности на отснятых видеокадрах по всей площади наблюдения. Сопоставительный анализ заявляемого решения с ближайшим аналогом показывает, что заявляемый способ отличается от известного введением новых технологических операций, обеспечивающих достижение свойств, закономерности которых не были известны и проявились в заявляемом объекте впервые. Действительно, основным селектируемым параметром ближайшего аналога является индекс цветности, в то время как в заявляемом способе информационным параметром является шероховатость растительного покрова. В свою очередь "шероховатость" зондируемого участка "извлекается" из текстуры его изображения. Это позволяет утверждать, что заявляемый способ удовлетворяет критерию изобретения "новизна". Наличие таких признаков, как введенные операции, связанные с функциональными преобразованиями текстуры получаемого изображения, позволяющими извлекать информацию о "шероховатости" зондируемых участков, расчет функции регрессии между типом растительности эталонных участков и параметров автокорреляции отраженного от данного типа растительности сигнала, позволяет сделать вывод о соответствии заявляемого технического решения критерию "Существенные отличия". Техническая сущность изобретения заключается в следующем. На настоящий момент основным селектируемым параметром при тематической обработке аэрокосмических снимков растительности является индекс цветности. Кривые спектральных коэффициентов яркости у всех древесных пород имеют одни и те же закономерности и отличаются друг от друга на единицы процентов. Это затрудняет тематическую обработку изображений по индексу цветности. Кроме физиологических факторов на величину КСЯ оказывает влияние и архитектура растений - высота, форма кроны, соотношение ярусов. Следовательно, "шероховатость" зондируемой поверхности также может служить селектируемым признаком растительности, т.е. целевую информацию о подстилающей поверхности содержит не только тон изображения, но и его текстура. Разреженность или густота полога, диаметры крон, их геометрическая повторяемость участвуют в формировании изображения и его крупных дискретных контрастных элементов, что отражается в качестве текстуры. Анализ известных определений понятия текстуры в области дешифрирования изображений показал, что общепринятого, формализованного определения текстуры как дешифровочного признака не существует. Известные определения понятия текстуры сводятся, как правило, к ее качественному описанию: мелко-, крупнозернистая, гладкая, линейчатая, нерегулярная и т.д. Авторами разработан способ получения количественных характеристик текстуры путем вычисления пространственного спектра анализируемого участка изображения с последующим уменьшением размерности его вектора признаков для идентификации до двух параметров: ширины автокорреляционной функции и ее максимальной амплитуды. Методы оценки "шероховатости" поверхности основаны на вычислении пространственного спектра соответствующего участка снимка. Из математики известно (см., например, Пискунов Н.С. Дифференциальное и интегральное исчисления, т. 2, М. : Наука, 1964, с.240 - 242), что любая функция может быть разложена в интеграл Фурье. Причем чем больше скорость изменения функции по координате, тем шире ее комплексный спектр. Следовательно, "шероховатость" или изрезанность полога определяют частотно-спектральный образ растительности. По определению, пространственный спектр вычисляют как двумерное Фурье-преобразование от функции интенсивности I (x, y) участка изображения:
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-5t.gif)
Однако непосредственно по виду огибающей пространственного спектра трудно идентифицировать тип растительности ввиду ее изрезанности, многомодности и неявно выраженных пиков. Другой характеристикой скорости флуктуации каких-либо процессов являются их автокорреляционные функции. Но непосредственно из функции интенсивности изображения автокорреляционная функция вычислена быть не может из-за потерь временной координаты при получении изображения (моментальности снимка). Авторами разработана последовательность операций, реализуемая посредством комплекса программно-аппаратных средств, позволяющая вычислять функцию автокорреляции сигнала изображения, включающая следующие промежуточные процедуры:
вычисление двумерного пространственного спектра G (fx, fy);
ввод физических размеров анализируемого участка;
нормирование и получение энергетического спектра сигнала изображения S (F);
вычисление автокорреляционной функции B(
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-6t.gif)
Оцениваемыми параметрами для идентификации типа растительности служат два значения функции автокорреляции - ее максимальное значение при
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
выбирают контрольные площадки с характерными типами растительности и проводят их тематическую оценку;
размещают видеокомплекс на автоматической поворотной платформе орбитальной станции;
отслеживают заданные районы наблюдения вводом уставок баллистических данных в систему управления автономной платформы;
осуществляют высокодетальную дистанционную съемку подстилающей поверхности, содержащей контрольные площадки, и записывают получаемые изображения на бортовой видеомагнитофон;
передают в сеансах связи по радиолинии в ЦУП полученные изображения и осуществляют их перепись на кассеты;
последовательно отбирают для анализа участки видеокадра и преобразуют их изображение в цифровую функцию интенсивности I (x, y);
вычисляют на основе программных методов быстрого Фурье-преобразования комплексный Фурье-спектр G (fx, fy) от функции I (x, y) в соответствии с формулой:
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-7t.gif)
преобразуют двумерный Фурье-спектр G (fx, fx) в энергетический спектр S(F) сигнала с учетом физических размеров введенного участка изображения;
вычисляют автокорреляционную функцию B(
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-8t.gif)
рассчитывают эталонную функцию регрессии между параметрами автокорреляционных функций контрольных площадок и типом растительности на них;
синтезируют из последовательно проанализированных участков мозаичную картину растительности на отснятых видеокадрах по всей площади наблюдения. Пример конкретного осуществления способа. Оценка метрических и информационных характеристик частотно-спектральных образов различных типов растительности и производных от них автокорреляционных функций проводилась на базе средств орбитальной станции "Мир". Реализация заявляемого способа нашла отражение в "Методике дистанционного решения лесоводческих и лесоэкологических задач видеоспектральными средствами орбитального комплекса "Мир", НПО "Энергия", М., 1994, 31 с. Заявляемый способ может быть реализован на базе устройства по схеме фиг.1. Фиг. 2-5 также поясняют предлагаемый способ. На автономной поворотной платформе 1, управляемой как из Центра управления полетом, так и от бортовой системы управления 2 орбитальной станции "Мир", установлен видеокомплекс 3 высокого пространственного разрешения типа "Астра". По программе или в режиме непосредственного телеуправления осуществляются съемка и запись изображений объектов 4 подстилающей поверхности, отслеживаемых видеокомплексом, на бортовой видеомагнитофон 5 типа "Betta Kam". Накопленная информация видеокадров сбрасывается в сеансах связи по радиолинии 6 типа "Ромашка" в ЦУП, где регистрируется на кассетах устройства 7. Наземный комплекс программно-аппаратных средств содержит ПЭВМ типа IBM PC/386/387 с набором периферийных элементов в составе устройства телеввода 8 изображений в компьютер, набор дискет 9 с эталонными видеоснимками ключевых участков, набор дискет 10 для записи файлов цифровой информации единичных участков. Последовательность цифровых файлов изображений единичных участков видеокадра вводится в оперативное запоминающее устройство центрального процессора (11) ПЭВМ с возможностью их визуализации на мониторе 12 типа Super VGA. Комплекс специализированных программ обработки записан на ПЗУ 13. Полученная в результате обработки снимков целевая информация в виде графиков расчетных функций выводится на принтер 14. Процедуру идентификации типа растительности по операциям заявляемого способа с помощью средств фиг.1 проиллюстрируем на примере отдельного видеокадра, представленного на фиг. 2. На фиг. 2 изображен участок подстилающей поверхности, попадающий в поле зрения бортового видеокомплекса "Астра". Для анализа типа растительности поле кадра последовательно разбивается на элементарные участки (штрихпунктирные линии на фиг. 2), определяемые требованиями точности идентификации и возможностями пространственного разрешения. Единичный элемент выбирается из условия соизмеримости с геометрическими размерами контрольной площадки (100 х 100 м). Каждый элементарный участок (единичный элемент) видеокадра телевводом 8 преобразуют в файл цифровой информации функции интенсивности I (x, y) участка изображения в формате 515 х 512 элементов. Операция преобразования изображения в файлы цифровой информации функции интенсивности реализуется на средствах телеввода (см., например, Система телеввода Panasonik. Инструкция использователю, техническое описание). Файлы цифровой информации с дисков 10 последовательно считываются в ОЗУ процессора 11. От каждого массива цифровой информации файла вычисляется двумерный пространственный спектр, используя комплекс специализированных программ, содержащихся в ПЗУ 13. О реализуемости данной операции быстрого преобразования Фурье см. , например, Марпл С.А. (мл), перевод с английского. Цифровой спектральный анализ. М.: Мир, 1990, с.77 - 79, алгоритм БПФ. Для структур, не имеющих преимущественных периодичностей, так называемых псевдоструктур, в Фурье-спектре будет присутствовать множество компонент. Обычно двумерный Фурье-спектр псевдослучайных структур выглядит в виде расплывчатой эллиптической фигуры. Дисперсность (шероховатость) структуры отражается в размахе Фурье-спектра. Следующий шаг анализа - оценка вклада пространственных структур различного размера в вид огибающей. Интегрированием двумерного Фурье-спектра по кольцевым сегментам получают одномерную функцию пространственного спектра S(1/
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115110/955.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115110/955.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115110/955.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-9t.gif)
где
v - скорость распространения физических волн. Поскольку съемка осуществляется с орбиты КА, а объекты подстилающей поверхности неподвижны, то за скорость распространения пространственных волн можно принять относительную скорость движения объектов и измерителя. Таким образом, пространственный спектр и спектр частот связаны зависимостью G(f) = v/S(1/
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115110/955.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-10t.gif)
т. е. энергетический спектр равен квадрату амплитудного спектра, усредненного на интервале наблюдения единичной реализации. Вычисление автокорреляционных функций B(
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-11t.gif)
В рассматриваемом примере конкретной реализации энергетические спектры S(F) двух типов растительности аппроксимированы функциями
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115887/2115887-12t.gif)
и вычислялись аналитически. Рассчитанные значения автокорреляционных функций представлены графиками фиг.4. Данная операция может быть автоматизирована на основе программных методов расчета, специализированных программ, записанных в ПЗУ 13. О реализуемости программных расчетов автокорреляционных функций см., например, Марпл С.А. Цифровой спектральный анализ и его приложение. Перевод с англ. М.: Мир, 1990, с.254. Из представленных зависимостей фиг. 4 следует, что автокорреляционные функции типов растительности существенно отличаются как амплитудой (мощностью переменных составляющих), так и интервалом корреляции (скоростью флуктуаций или шероховатостью). При наличии тарировочной характеристики единичные участки растительности могут быть достоверно дешифрованы. Тарирование тракта для получения дешифровочной характеристики осуществляется путем ввода в ПЭВМ и обработки снимков с дискет 9 ключевых участков. В результате осуществления данной процедуры получают совокупность реализаций автокорреляционных функций, соответствующих типам растительности ключевых участков. По пиковым значениям амплитуд B(
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
![способ идентификации типов растительности, патент № 2115887](/images/patents/356/2115084/964.gif)
Класс G01C11/00 Фотограмметрия или видеограмметрия, например стереограмметрия; топографическая съемка местности с помощью фотографирования
Класс G01C11/04 расшифровка изображений
Класс G01C11/10 с использованием вычислительных устройств для контроля положения изображения
Класс G03B37/00 Панорамная фотосъемка и съемка широкоэкранных кинофильмов; фотографирование обширных поверхностей, например для топографической съемки; фотографирование полостей, например внутренней поверхности труб
Класс A01C23/00 Устройства, специально предназначенные для распределения навозной жижи или других жидких удобрений, в том числе аммиачной воды, например резервуары для перевозки, повозки для разбрызгивания
Класс G06F17/14 преобразования Фурье, Уолша или аналогичные преобразования
Класс G06F17/15 вычисление корреляционных функций
Класс G06F17/16 матричные или векторные вычисления
Класс G06F17/17 вычисление функций приближенными методами, например интерполяцией или экстраполяцией, сглаживанием, методом наименьших квадратов