Лекция "Цифровая картография"

Открытая лекция «Цифровая картография»

24 июня в 13.00 я приглашаю всех в бар «Бакалавриат» (Петербург, Улица Маяковского, дом 10. (метро Маяковская). Во-первых, зайти воскресным днем в бар — сама по себе неплохая идея, а во-вторых, как раз в это время я буду выступать там с открытой лекцией о современной цифровой картографии.

Это будет рассказ о том, как устроены современные картографические сервисы (Google maps, Яндекс-карты, OpenStreetMap и др.). Почему спутниковые снимки искажают действительность, что такое «геоданные» и «картостиль». Обсудим, как картографы сделали большой крюк, прежде чем вернулись к привычным способам наблюдения за миром. Рассмотрим практическую значимость геоинформационных систем и научимся правильно понимать значение географических координат. Лекция рассчитана на неподготовленных слушателей, поэтому страшных и непонятных вещей там не будет.

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

Еще раз: Бар «Бакалавриат» на улице Маяковского 10. В час дня 24 июня. Вход бесплатный — приходите!

P.S. При репосте этого сообщения постгрес начинает работать в 2.3 раза быстрее — проверенный факт.
P.P.S. Дабы избежать организационных трудностей, рекомендую любым способом сообщить мне о том, что вы придете. Это желательно, но не обязательно.

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

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

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

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

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

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

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

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

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

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

Точки на карте

Основы панка. Бетмен с яйцами

«Вон геолог с картой идет. Сейчас дорогу спрашивать будет»
Анекдот от опытного геолога

С геологической картой района работ я познакомился на следующий день после приезда. Геологи усмехались, глядя как я держу ее вверх-ногами, но позже оказалось, что как ее ни крути, один хрен содержимое карты не сходилось с реальностью.

Как создается геологическая карта? Очень просто. Берете картину любого абстракциониста, например Девида Духовны, который, как утверждал один из номеров газеты «Спидинфо», творит свои полотна, размазывая краски на холсте голой жопой. Расставляете на этой картине индексы комплексов и свит, добавляете легенду и карта готова:

Геологическая карта
Я не утрирую – геологическая карта как никакая другая далека от истины и представляет собой скорее продукт воображения, чем реальную интерпретацию полевых данных. Смотришь на карту — вот он, Сумий-Сариолий, представленный в виде двух яиц Селецкой свиты. Или вот голова Бетмена в подробностях, сложенная нюкозерскими гранитами:

Голова Бетмена
Все так детально, красиво прорисовано. Но приходишь на место и весь день пытаешься найти хоть какое-нибудь обнажение. Какой там, нахуй Бетмен с яйцами. Порой даже просто подтвердить наличие указанной структуры невозможно без бурения — кругом одни сплошные морены, озы, камы, болота и непросматриваемые вырубки, обросшие густой ольхово-березовой хуетой.

Болотина

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

И странно еще, что какие-то полевые данные вообще удается собрать, слишком уж много времени в маршруте уходит на откровенную хуету – разведку дорог и последующую прокладку маршрутов. Казалось бы, кто мешает перед полем подготовить абрисы дорог? Хрен с ним, нет денег на свежие спутниковые снимки, но хоть в том-же OpenStreetMap-е дороги-то можно было обрисовать? Чем так возбуждают карты генштаба, которые наполовину не соответствовали действительности при издании, а на вторую половину устарели еще в прошлом тысячелетии?

Геологи нихуя не картографы. От слова совсем. Даже постоянная работа в гисах не может привить им современной парадигмы картографии. Иногда это выглядит достаточно мило и винтажно: помнить номенклатуры, стандартные ряды масштабов, ходить по азимуту, периодически прикладывая транспортир к карте и навигатору на телефоне (я сам охуел когда первый раз увидел) и раз за разом оказываться в местах, в которых не планировал оказаться.

Однажды я перепил пива и подготовил демку на район. Хочешь горизонтали вытаскивай, хочешь DTM-фильтром склоны выделяй. Вместо этого последовал вопрос: «Демка – это что, демо-версия?». Да, блядь. Я подготовил демо-версию нашего района. Триал – главное все успеть, а то, как пизданет через тридцать дней по району работ очередным метеоритом.

Но хуй-то со всеми этими демками, фильтрами и парадигмами картографии. Хуже всего то, что многие (замечу, не все) из геологов нихуя не умеют ходить по навигатору. Зато все берут азимут с точностью до градуса — похуй, что трек по форме больше похож на кардиограмму при аритмии.

— Нам по азимуту триста четырнадцать градусов, вон туда идти
— Обратно через болото, что-ли? Оттуда же пришли
— Эээ, ну значит тогда вот туда…
— Бля, отдай мне навигатор — сам поведу.

Просто пиздец какой-то. Идешь с геологом и гадаешь — куда он тебя сейчас заведет и как его потом из этой жопы вытаскивать. Прямо как мой пес по кличке Лишай. Тот тоже всегда выбирал ебанутые маршруты, но в отличие от геологов всегда знал как вернуться. 

Хождение по азимуту — это как алхимия в средние века — все верят, но кто-бы не пробовал, получает в результате какое-то говно. Причем, так было есть и будет во все времена. Именно поэтому глобальная навигация и возникла.