Перевод отметок рельефа из Autocad в табличный вид

При подготовке данных для моделирования весьма часто возникает необходимость в конвертации отметок рельефа из геоподосновы, сохраненной в формате Autocad DWG, в табличный вид (типа X, Y, Z) для того, чтоб скормить эти данные какому-нибудь интерполятору (Surfer или встроенный в PmWin “Field Interpolator”).
Я почти уверен точно знаю, что эта задача может быть легко и быстро решена с помощью самого Autocad — достаточно запустить соответствующую программу на LISP-е и радоваться жизни. К сожалению, я LISP-а не знаю, да и вообще не являюсь большим специалистом в автокаде.
В современных версиях автокада эта проблема решается еще проще: через инструмент, расположенный в пункте меню Tools\Data Extraction. Инструмент довольно мощный, но в нестандартных случаях возможно придется повозиться.
Но проблему как-то надо решать. Я предлагаю использовать для этого MapInfo (согласен, для кого-то это выглядит сменой шила на мыло). Далее по пунктам:

1. Конвертируем слой с отметками (желательно, чтоб слой содержал только отметки в текстовом виде, без самих точек) из формата Autocad DWG  в формат MapInfo TAB (с помощью встроенного в MapInfo мини-приложения Universal Translator).
2. Подчищаем полученную таблицу от нетекстовых элементов: это можно сделать несколькими способами, наиболее удобный и быстрый — мини-приложение MapCad, но можно и с помощью Query Select и функции ObjectInfo(obj, 1), но там придется сначала создать дополнительный столбец в таблицу, занести него результат выполнения функции ObjectInfo(obj, 1), а уж потом делать Query Select по этому столбцу, выбирая значения, отличные от 10 (а 10 — это как раз текстовые).
3. Еще разок запускаем Update Column (создайте новый столбец с типом float или смените тип существующего столбца) с той же функцией, но с другими параметрами: ObjectInfo(obj, 3). Если все сделано правильно, то в вашей таблице появится столбец типа float с отметками рельефа.
4. Запускаем мини-приложение Coordinate Extractor: в таблице теперь будут столбцы с координатами центра текстовой подписи отметки рельефа. Вот тут важно отметить явный недостаток рассматриваемого метода: наши точки будут немного смещены относительно реальных отметок — ровно на столько, на сколько отличаются координаты середины метки от координат точки замера. Если вы страдаете перфекционизмом, то эту проблему можно решить с помощью простейших математических операций со свежеполученными координатами.
5. Запускаем Table/Create Points, если хотим заменить тектовые метки на точки (а уж сами метки пусть MapInfo своими силами рисует, благо соответствующий столбец в таблице уже есть).
6. Экспортируем полученную таблицу в нужный текстовый формат (txt или csv).

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

Сколько метров в градусе

У моего блога небольшой юбилей — 100-я запись. Чуть более назад, когда я начинал это дело, мне и в голову не могло придти, что рубеж в сотню постов вообще будет преодолен. В голове тогда было идей на пару десятков записей, а вот поди ж ты. А поделиться в этот «зднаменательный» пост я хотел небольшим “know how”, которое будет очень полезно тем, кто столкнулся с реальными задачами, а не тестовыми примерами из обучающих наборов.
Большая часть программ для гидрогеологического моделирования с географическими координатами, выраженными в градусах, работать либо не умеют, либо умеют, но так, что лучше бы не умели вовсе. В то же время, общедоступные подборки векторных карт от сообщества OpenStreetMap, которые очень удобно использовать в качестве подложек для моделей, на беду представлены именно в градусной координатной сетке (Longitude/Latitude WGS84, если «выражаться» в терминах MapInfo). Так что же делать? Надо поменять градусную координатную сетку на метровую. В том же MapInfo это делается путем пересохранения нужного слоя в проекции Universal Transverse Mercator WGS84. Главное тут — не ошибиться сильно с зоной (Москва, в частности, находится в 37-й зоне), но в этом вопросе поможет разграфка зон UTM (via: Часто задаваемые вопросы по координатам, проекциям, системам координат).

Как ускорить работу MapInfo

Думаю, для тех, кто много работает с этой программой, мой следующий совет новостью не станет, хотя лично я о этой функции частенько забываю. Так вот, если вы заметили, что MapInfo начал ни с того ни с сего заметно подтормаживать в процессе отрисовки карты при зуммировании или перетаскивании, то вам поможет функция Table\Maintenance\Pack Table, «натравленная» на самые большие открытые таблицы. Скорость работы возрастает на порядки. Не забудьте предварительно сохранить все таблицы и рабочее пространство (особенно его, т.к. в процессе «упаковки» таблицы будут закрыты).

Импорт WMF в MapInfo

Возникла у меня не совсем типичная задачка: прислали мне векторную карту, сохраненную в формате WMF (Windows Media File). Мне для работы желательно загнать эту карту в MapInfo. Решение «в лоб»: привязать wmf, как растровую подложку. Не сработало. MapInfo в принципе это умеет делать, но с большими файлами (а у меня он весит почти 5 мегабайт) не дружит — показывает от силы половину карты, а то и того меньше. К тому же, как-то некомильфо получается — векторный исходник подкладывать растром в векторную же программу.
Пришлось идти по длинному пути: сначала открываем наш WMF в AutoCAD-е с помощью команды _wmfin, не меняя масштаба по осям, потом командой _explode разбиваем блок (на всякий случай, т.к. у MapInfo случаются затыки с импортом блоков), сохраняем файл в виде DXF. Транслируем этот DXF в MapInfo. С помощью MapBasic-утилиты Register Vector перемещаем и масштабируем карту. Увы, если у вас был на карте текст, то он скорее всего потеряется. Особенно, если он написан на кириллице. Для решения этой проблемы, по всей видимости, придется перемещать и масштабировать карту сразу в AutoCAD-е, еще до импорта ее в MapInfo.

Единая государственная картографическая основа г. Москвы, предоставляемая ГУП "Мосгоргеотрест".

Ух-ты, смотрите что на торрентах выложили. Там на первом этапе возникнет необходимость всосать файлы в формате mif в MapInfo. Если это делать в лоб, то можно сильно утомиться. Для решения этой проблемы существует удобный набор скриптов: BatchTools.

Формат: MapInfo, Access, Oracle, DWG, DGN, DXF
Описание: Единая государственная картографическая основа г. Москвы, предоставляемая ГУП "Мосгоргеотрест".
"В соответствии с пунктом 1.1 Постановления Правительства Москвы от 19.01.1999г. №24 "О внедрении Единой государственной картографической основы г.Москвы для решения задач управления городским хозяйством с использованием автоматизированных технологий" Мосгоргеотрест Москомархитектуры является единственной уполномоченной Правительством Москвы организацией по созданию и ведению Единой государственной картографической основы г. Москвы на базе данных цифровых топографических планов, карт и схем масштабного ряда 1:2000 - 1:25000, формируемой в Мосгоргеотресте в рамках Геофонда г. Москвы, и несет ответственность за достоверность и юридическую чистоту информации, содержащейся в Единой государственной картографической основе г. Москвы."

Pitney Bowes Mapinfo Professional v11

Вышла очередная версия Mapinfo Professional:

Work with data you can truly trust with MapInfo Professional

* Quickly and easily find and open the data you want
* Obtain its source, vintage, resolution and owner with a click of a button.
* Add your map data to MapInfo Manager's library right from MapInfo Professional

New Browser Window means easier data creation, analysis and manipulation

* Speed up tabular data entry with flexible browser navigation options
* Hide fields, add calculated fields or sort your table right from the right-click menu.
Data analysis was never this easy
* Visibility enhancements such as alternating row colours, full row selections, and quick access to font size and style setting are sure to speed up any data analysis

Enhanced support of 64 bit operating environments

* Leverage up to double the application memory space (4Gb) than before
* More memory means ability to perform larger operations

Инструментарий — визуализация

Третий этап работы над геофильтрационной моделью: оформление отчета и визуализация результатов. Отчет я предпочитаю писать в Word. Безусловно есть и более удобные инструменты, а главное — более гибкие, но клиенты не особо жалуют pdf, а альтернативы еще из детских штанишек не вылезли.
Подготовка табличных данных — ясное дело — Excel. А вот графики предпочитаю доверять Golden Software Grapher. Опять же, её есть гораздо более мощные программы. Тот же Origin, к примеру, так что если есть время их их изучать, то никто не запрещает. Коллеги утверждают, что оно того стоит.
Это все была присказка, главное будет дальше. Ведь самое важное, сложное и интересное — визуализация модельных данных. Для простеньких задач, не требующих прохождения экспертизы, вполне бывает достаточно встроенных средств программ-препроцессоров (тут они выступают уже в качестве «постпроцессоров»). Результат будет выглядеть бледненько, но для понимания картины вполне достаточно, особенно если не полениться и дорисовать по-человечески условные обозначения. Увы, для уважаемых экспертов понимание вовсе не главное, на первом месте — «чтобы было красиво». А для этого встроенных средств как правило недостаточно. Хотя GMS или тот же Visual Modflow достаточно продвинуты в этом плане, но для получения удовлетворительного результата придется изрядно повозиться. Впрочем, некоторых «фишек», вроде «свечек» с невязкой по наблюдательным скважинам в GMS, иногда очень не хватает.
Так вот. Все известные мне программы для моделирования умеют экспортировать результаты расчетов в виде таблицы x, y, H, где x, y — координаты центра ячеек, а H — значение напора (или другого расчетного параметра). Полученная таблица отлично поддается интерполяции с помощью того-же GS Surfer (чего-то мой бложик становится все больше и больше похож на рекламную компанию продукции Golden Software). Полученные изолинии из Surfer можно этим экспортировать в любую GIS. Мне нравится Mapinfo. Существует мнение, что ArcGIS существенно мощнее, но я как-то сломался в нем разбираться, MI существенно проще для новичка, а главное - для поставленных задач больше и не надо.
В результате получаются карты типа вот такой (на карте показано влияния «стены в грунте» на водоносный горизонт):

Инструментарий

Мне очень хочется рассказать про все имеющееся в арсенале специалиста гидрогеолога-модельера программные средства, но это займет не один пост. И, пожалуй, даже и не два — посмотрим, как дело пойдет.
Поэтому начну с того, с чем регулярно приходится работать лично мне. Опишу свой, как говорится, workflow.
Создание геофильтрационной модели как правило происходит в несколько этапов:
1) Сбор и анализ имеющихся материалов. Если со сбором всё более или менее понятно — поднимаются связи в региональных фондах, достаются полулегально (о причинах этого маразма я обещаю как-нибудь рассказать) колонки скважин и отчеты по изысканиям, детальные карты и т.п.; то с анализом всего этого, нарытого праведно и неправедно добра, часто возникают затруднения. Лично мне удобнее работать с цифровыми данными, поэтому в первую очередь я перевожу все бумажные материалы в цифровой вид. Сканирую текст отчетов, карты (о них, кстати, будет еще много записей в этом блоге, ибо тема весьма благодатная), вбиваю в табличном виде содержимое колонок и так далее. Инструменты на данном этапе использую следующие: программа для сканирования и распознавания текста; Golden Software Surfer — для интерполяции точечных данных; GIS Mapinfo — привязка и векторизация карт и схем, а также создание новых (об этой замечательной программе я тоже еще буду много рассказывать); Microsoft Excel — ну, это понятно: работа с любыми табличными данными. Много еще каких программ используются на этом этапе по мере необходимости, но эти четыре поменяются мною всегда.
2) Второй этап — самый муторный, заключается в подготовке полученных на первом этапе данных к экспорту в программу моделирования. Тут используются все те же Excel и Surfer. Ну, и встроенные средства моделирующих программ, о которых позже.
Про следующие этапы — непосредственно моделирование и визуализацию результатов я расскажу в следующий раз.