Surfer 17

Оказывается, в последней версии Surfer-а появился Кокригинг. Это довольно полезный в нашем деле алгоритм интерполяции, советую ознакомиться.

Перевод отметок рельефа из 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).

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

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

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