CC BY
49
Вычислительные технологии
Том 12, Специальный выпуск 2, 2007
РАСПОЗНАВАНИЕ ФОРМАЦИЙ ЛЕСНОЙ РАСТИТЕЛЬНОСТИ С БЛИЗКИМИ СПЕКТРАЛЬНО-ЯРКОСТНЫМИ ХАРАКТЕРИСТИКАМИ ПО ДАННЫМ СЪЕМКИ СО СПУТНИКА LANDSAT 7 ЕТМ+
Н. В. Гопп
Институт почвоведения и агрохимии СО РАН, Новосибирск, Россия
e-mail: [email protected]
Е.А. Куликова, И. А. Пестунов, Ю.Н. Синявский,
В. В. Смирнов
Институт вычислительных технологий СО РАН, Новосибирск, Россия e-mail: pestunov@ict. nsc. ru, valentin@sbras. ru
Applicability of automatic classification algorithms (clusterization) for decomposition of the forest vegetation formations with close spectral-brightness characteristics using Landsat Т ETM+ data is considered.
Введение
При решении задач, связанных с анализом и распознаванием данных дистанционного зон пировали я (ДДЗ), в большинстве случаев отсутствует не только априорная информация о вероятностных характеристиках классов, но и классифицированная обучающая выборка. Естественный выход из этой ситуации — применение алгоритмов автоматической классификации [1, 2], Однако кластеризация ДДЗ без предварительной обработки зачастую приводит к неудовлетворительным результатам при разделении природных объектов с близкими спектрально-яркостными характеристиками, так как выбор параметров, позволяющих разделить эти классы, ведет к чрезмерной раздробленности полученной картосхемы,
В работе представлена методика, автоматизированной обработки данных спутниковой съемки, позволяющая успешно решать подобные задачи. Методика, прошла апробацию при исследовании лесной растительности Болотнинского района Новосибирской области по данным съемки со спутника І.аїкікаї 7 ЕТМ+,
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2007.
1. Методика исследования
Предлагаемая процедура обработки изображения начинается с пороговой сегментации по нормализованному вегетационному индексу
где Iir, Ir — значения сигналов, регистрируемые соответственно в 3-м и 4-м каналах. Для более наглядного представления гистограммы, значения вегетационного индекса преобразуются к интервалу [0, 255], Пороговые величины для отсечения непокрытых растительностью территорий (водная поверхность, минерализованные участки, постройки и т, п.) определяются пользователем в диалоговом режиме на основе визуального анализа гистограммы значений NDVI, Затем производится линейное растяжение динамических диапазонов значений спектральных признаков на весь допустимый диапазон [0, 255], После сегментации для всех возможных наборов признаков вычисляется значение OIF (Optimum Index Factor) [3]:
где S — набор каналов; a(i) — среднеквадратичное отклонение значений спектральных яркостей для i-го ran ала; r(i,j) — коэффициент корреляции значений спектральных яркостей между i-м и j-м каналами.
Затем определяются наборы признаков с наибольшими значениями OIF, Среди них на основе анализа спектральных откликов растительности определяется набор признаков, наиболее подходящий для решения конкретной задачи. Таким образом находится набор слабо коррелированных и информативных признаков.
Полученный набор признаков используется для автоматической классификации изображения при помощи непараметрического алгоритма ССА [4],
2. Объект исследования
Исследовались формации растительности Болотнинского района Новосибирской области, Этот район имеет сложное геологическое и тектоническое строение и находится в зоне сочленения мегаструктур разной геодинамической природы Западно-Сибирской плиты и Алтае-Саянской складчатой области [5], Изучаемый район ограничен 83°49/53// и 84°29/7// в.д. и 55°58/42// и 55°52/42// с.ш. Прямоугольниками выделены ключевые
участки с координатами 83°50/43//-83°55/20// в.д., 55°54/45//-55°53/15// с.ш. и 84°23/7//-° / // ° / // ° / //
Псследуемая территория находится на стыке двух геоморфологических районов. Большая ее часть относится к долине реки Оби, меньшая — к Сокурскому геоморфологическому району. Долина реки Оби выражена четырьмя надпойменными террасами, которые высоко подняты над современной поймой и тянутся полосой вдоль Оби, достигая на отдельных участках совокупной ширины в несколько десятков километров. Растительный покров изучаемой территории неоднороден и образован разными типами
NDVI
E*(i)
OIF
растительности, характерными для южтто-таежттой подзоны, лесостепной зотты, а также итттразотталытой растительности долины реки Оби.
Южтто-тасжттая подзона представлена незначительно и занимает северо-восточные участки изучаемой территории. В основном это тсмттохвойпые южтто-таежттые травятто-зслсттомотттттыс леса (пихта, кедр, ель).
На ттадпоймстшых террасах реки Оби расположены крупные болотные массивы, представленные травяными ттизиппыми болотами (тростниковыми, осоковыми). По периферии болот распространены заболоченные луга. Сосновые боры занимают дрсвттие террасы правого берега Оби тт ее крупных притоков. Среди ттих встречаются сосновые травяные тт травятто-кустарттттчковыс (брусничные, черничные) леса.
3. Результаты и обсуждение
Обработке подвергался фрагмент многозонального снимка Landsat 7 ЕТМ+ с пространственным разрешением 30 м, полученного 1 августа 1999 года тт прошедшего радиометрическую тт геометрическую коррекции (рис. 1). Для дешифрирования фрагмента использовались данные наземных мартттрутттых наблюдений па ключевых участках, собранные с помощью системы глобального позиционирования (GPS).
На рис. 2 представлены картосхемы, полученные методом визуалыто-ипструмспталь-ттого дешифрирования выделенных ключевых участков. С помощью этих картосхем по данным многозонального снимка были сформированы обучающие выборки для представленных типов лесттой растительности во всех семи спектральных диапазонах. В результате анализа спектралыто-яркостпых характеристик лесттой растительности, построенных по ттаидсттпым обучающим выборкам, установлено, что различные формации
Рис. 1. Фрагмент изображения, полученного по данным Landsat 7 ЕТМ+, с нанесенной географической координатной сеткой; прямоугольниками выделены ключевые участки, содержащие формации растительности с близкими спектралъно-яркостными характеристиками: 1 ключевой участок в долине реки Оби с сосновым лесом, эдификатором служит сосна обыкновенная (Pinus silvestris); 2 ключевой участок в южно-таежной подзоне с темнохвойными лесами, эдификаторы: пихта (Abies sibirica), ель (Picea obovata), кедр (Pinus sibirica)
растительности имеют высокую спектральную селективность. Исключением являются темпохвойпые (пихта, кедр, ель) и сосновые (травяные и травяпо-кустарпичковые) леса, имеющие близкие характеристики (рис, 3), Указанная особенность существенно усложняет разделение этих типов леса без предварительной обработки снимка и выбора информативных признаков.
а
О 0.5 1.0 1.5 2.0 км
Сосновые травяные; и травяно-кустарни чковыо леса Березовые и березово-осиновые травяные леса
Луговая и кустарниковая растительность в нойме реки Пес чаныо отложения Река Обь
б
84°24'00" 84°25'12" 84°26'24
I-----Ь----ь-----Ь-----I
0 0.5 1.0 1.5 2.0 км
Томнохвойно-борозово-осиновыо леса Темнохвойные (пихтовые, кедровые, еловые) леса
Луговая разнотравная растительность
Березовые травяные леса Сельскохозяйственные угодья Озеро Большое Черное
Рис. 2. Картосхемы участков 1 и 2 (а и б соответственно), полученные на основе визуальноинструментального дешифрирования снимка Ьапска! 7 ЕТМ—
60
О
1 2 3 4 5 7 6
Номер канапа
Сосновые травяные и травяно-кустарни чковыо леса Темнохвойные (пихтовые, кедровые, еловые) леса
Березовые и березово-осиновые травяные леса
Томнохвойно-борозово-осиновыо леса
Рис. 3. Характеристики спектрального отклика лесной растительности
О 50 100 150 200 250
Рис. 4. Гистограмма значений вегетационного индекса
Рис. 5. Результаты пороговой сегментации но вегетационному индексу
В соответствии с предложенной методикой была построена гистограмма значений вегетационного индекса (рис. 4) и произведена сегментация изображения с пороховыми значениями 140 и 240 (рис. 5). Это позволило снизить влияние не покрытых растительностью территорий (па рис. 5 изображены белым цветом) при определении информативных признаков с помощью величины OIF.
В таблице представлены значения OIF для различных наборов признаков, вычисленные по сегментированному изображению. Анализ спектральных характеристик (см. рис. 3) показывает, что для разделения темпохвойпых и сосновых травяных лесов каналы 1 и 2 малоипформативпы, а капал 3 значительно информативнее, чем капал 7. По этой причине среди наборов с наибольшим значением OIF для классификации выбран набор признаков {3, 4, 5}.
По выбранным признакам проводилась обработка исследуемого фрагмента с помощью алгоритмов ССА (рис. 6) и ISODATA, реализованного в пакете ERDAS Imagine (рис. 7). На рис. 8 и 9 ноказапо, что результаты автоматической классификации хорошо согласуются с результатами визуально-инструментального дешифрирования.
Значения OIF для различных наборов из трех признаков
Набор признаков Величина OIF Набор признаков Величина OIF Набор признаков Величина OIF Набор признаков Величина OIF
{1,4,5} 29.7 {1,2,4} 25.4 {1,2,5} 22.6 {1,2,7} 16.2
{1,4,7} 27.1 {1,3,4} 25.1 {1,5,7} 19.5 {3,5,7} 15.2
{2,4,5} 26.8 {2,4,7} 24.1 {1,3,5} 18.7 {1,3,7} 14.9
{4,5,7} 26.2 {3,4,7} 23.0 {2,5,7} 16.7 {1,2,3} 14.7
{3,4,5} 26.1 {2,3,4} 22.7 {2,3,5} 16.2 {2,3,7} 12.6
Рис. 6. Результаты классификации алгоритмом CCA: 1 территории, исключенные из обработки; 2 сосновые травяные и травяно-кустарниковые леса; 3 темнохвойные (пихта, ель, кедр) леса; 4 смешанные томнохвойно-борозово-осиновыо леса; 5 редкостойные томнохвойно-борозовыо местами заболоченные леса; 6 разнотравные заболоченные луга; 7 березовые травяные леса; 8 травяные низинные болота; 9 травяно-моховые низинные болота; 10 редкостойные березовые леса на переходных болотах; 11 пойменная луговокустарниковая растительность; 12 лугово-разнотравная растительность; 13 мелкие классы (суммарной площадью менее 0.4 % от изображения)
Рис. 7. Результаты классификации алгоритмом ШСШАТА: 1 территории, исключенные из обработки; 2 сосновые травяные и травяно-кустарниковые леса; 3 темнохвойные (пихта, ель, кедр) леса; 4 смешанные томнохвойно-борозово-осиновыо леса; 5 местами заболоченные томнохвойно-борозовыо леса; 6 редкостойные томнохвойно-борозовыо леса; 7 разнотравные заболоченные луга; 8* — травяные низинные болота; 9* — травяно-моховые низинные болота; 10* — березовые травяные леса; 11* — редкостойные березовые леса на переходных болотах; 12 пойменная лугово-кустарниковая растительность
* Формации растительности, для которых не удалось получить удовлетворительное разделение
Рис. 8. Ключевой участок 1: фрагмент исходного изображения после пороговой сегментации (а), результаты визуально-инструментального дешифрирования (б) и результаты классификации с помощью алгоритмов ССА (в) и ISODATA (г)
Рис. 9. Ключевой участок 2: фрагмент исходного изображения после пороговой сегментации (а), результаты визуально-инструментального дешифрирования (б) и результаты классификации с помощью алгоритмов ССА (в) и ISODATA (г)
Заключение
Проведенные исследования показывают, что применение распространенного алгоритма ISODATA, а также непараметрического алгоритма ССА непосредственно к изображению не позволяет разделить формации растительности, близкие по спектрально-яркостным характеристикам. Предварительная обработка изображения позволила улучшить результаты классификации с помощью алгоритма ISODATA, Однако и в этом случае не произошло полного разделения изучаемых формаций растительности. Предложенная методика на основе алгоритма ССА позволила разделить темнохвойные (пихта, кедр, ель) и сосновые (травяные и травяно-кустарничковые) леса. Результаты этого алгоритма характеризуются хорошей проработкой мелких деталей изображения, которые зачастую не выделяются при использовании алгоритма ISODATA,
Исследование формаций растительности Болотнинского района Новосибирской области позволяет заключить, что результаты обработки хорошо согласуются с результатами визуально-инструментального дешифрирования.
Список литературы
[1] Mercer D.P. Clustering Large Datasets. Linacre College, 2003. http://www.stats.ox.ac.uk/ mercer/documents/Transfer.pdf
[2] SCHOWENGERDT R.A. Remote sensing, models and methods for image processing / 2nd ed. Acad. Press, 1997. 522 p.
[3] Chavez P.S., Berlin G.L., Sowers L.B. Statistical method for selecting Landsat MSS ratios // J. of Applied Photographic Engineering. 1984. Vol. 8. P. 23-30.
[4] Пестунов И.А., Будкина E.A., Синявский Ю.Н. Алгоритм кластеризации многоспектральных изображений на основе формирования сеточной структуры в пространстве признаков // Труды Международной конференции “Вычислительные и информационные технологии в науке, технике и образовании”. Павлодар, 2006. Т. 2. С. 124-131.
[5] Свиридов В.Г., Краснов В.И., Сурков B.C. и др. Геологическое строение и полезные ископаемые Западной Сибири. Т. 1: Геологическое строение. МПР РФ. Новосибирск: Нзд-во СО РАН, НИЦ ОНГГМ, 1999. 228 с.
Поступила в редакцию 30 августа 2007 г.
