Полигоны и лесополосы

Разделение лесов, кустов и полей на старых ландсатах при помощи растра высот

Вот смотрите. Берем седьмой Ландсат за июль 1999 года на территорию Игнатьевского сельского поселения в Адыгее и строим по нему растр вегетационного индекса. Он же — стандартный NDVI, который считается по формуле (nir-red)/(nir+red), где nir — инфракрасный канал, а red — красный. В результате имеем вот такую хренотень:
NDVI-index

А теперь берем снимки Alos, которые улучшенный Астер, который, в свою очередь, отснял эту же территорию несколько месяцев спустя:
Alos

Теперь открываем JOSM и рисуем осевые линии лесополос:
JOSM

Которые затем буферизуем в QGis-e до ширины 100 м (что-бы ухватить разные пикселы Ландсата) и режем на стометровые отрезки поперек. Это в QGis требует особой акробатики, но мы справляемся. Самое главное, что-бы квадраты обязательно выходили за границы лесополосы:
лесополосы

Теперь считаем зональную статистику по двум растрам. С данных Alos снимаем размах, который соответствует разности максимальной и минимальной высоты в каждом полигоне (перепад высот). С растра NDVI снимаем среднее значение индекса в каждом полигоне. Экспортируем все это во внешнюю таблицу и строим примитивный график:
график

Благодаря тому, что полигоны превосходят по ширине лесные полосы, они захватывают области перепада высот между поверхностью земли и вершиной деревьев. Таким образом, мы делим область внутри полигонов на три нечеткие группы:
график с пояснениями

Остается только повторить процедуру в обратном порядке и получим простейшую карту растительности.

Уверен, что эти размышления уже много раз описаны, но как-то все описания мимо меня прошли. Но самое главное в другом. Даже если отмести всякий лженаучный бред про информационные поля, все-равно выходит, что информация очень тесно связана между собой. Остаточная информация хранит свойства объекта также, как ящик с песком из которого вытащили отлитую деталь. Собственно, из-за размышлений на эти философские темы, я и проебал заявку на тендер.

Добавить комментарий

Ваш e-mail не будет опубликован.

2 comments

  1. freeExec:

    Десяти метровые кусты, нормально чё.