Программные средства и технологии, используемые для обработки геолого-геофизической информации:стандартные программы MSOffice;
программы статистической обработки информации
(Statistica, Коскад);
программы компьютерной графики:
стандартные программы (CorelDraw, Photoshop…);
программы инженерной графики (Surfer, Grapher, Voxler,
Strater);
системы автоматизированного проектирования
(AutoCAD и др.);
специализированные системы обработки и
интерпретации геолого-геофизической информации;
системы комплексного анализа и интерпретации
геолого-геофизических данных;
геоинформационные системы.

План прохождения дисциплины
Содержание курса:
Баллы
1. Основы картопостроения в программном пакете
Surfer (Golden Software).
40 (16)
2. Создание трехмерных моделей полей в программе
Voxler (Golden Software).
20 (8)
3. Основы проектирования в системе Autocad (Autodesk)
40 (17)
4. Решение геологических задач в геоинформационной
системе ArcGIS (ESRI)
30 (12)
5. Создание 3D модели залежи и подсчет запасов в
системе Micromine (Micromine).
30 (12)
Итоговая аттестация
40 (17)

ТЕМА №1.

Основы картопостроения в
программном пакете Surfer

Программа Surfer (Golden Software, США)

Основным назначением пакета является построение
карт поверхностей z = f(x, y).
3D проекция

Интерфейс программы

Панели
инструментов
Меню
программы
Окно Plot
Окно Worksheet
Менеджер
объектов

Структура системы

Программа включает 3 основных
функциональных блока:
1. построение
цифровой модели
поверхности;
2. вспомогательные операции с цифровыми
моделями поверхности;
3. визуализация поверхности.

Построение цифровой модели поверхности
Цифровая модель поверхности Z(x, y) представляется
в виде значений в узлах прямоугольной регулярной сетки, дискретность
которой определяется в зависимости от конкретной решаемой задачи.
y
x ≠ y
x
y
z1
z5
z9
z13
z17 узел
z2
z6
z10
z14
z18
z3
z7
z11
z15
z19
z4
z8
z12
z16
z20
x

Для хранения применяются файлы типа [.GRD] (двоичного или
текстового формата).
количество ячеек по осям X и Y
min и max значения X, Y, Z
линия y
(Y=const)
линия x
(X=const)
Программа Surfer позволяет использовать готовые цифровые модели
поверхностей в форматах других систем USGS [.DEM], GTopo30 [.HDR],
SDTS [.DDF], Digital Terrain Elevation Model (DTED) [.DT*] .

В пакете реализованы 3 варианта
получения значений в узлах сетки:
по исходным данным, заданным в произвольных точках области (в
узлах нерегулярной сетки), с использованием алгоритмов
интерполяции двухмерных функций;
вычисление значений функции, заданной пользователем в явном виде;
переход от одной регулярной сетки к другой.

Создание grid по нерегулярному набору данных
Исходные данные:
Таблицы форматов [.BLN], [.BNA], [.CSV], [.DAT], [.DBF], [.MDB], [.SLK],
[.TXT], [.WKx], [.WRx], [.XLS], [.XLSX]
Данные XYZ

Выбор
данных
Пункт меню Grid>Data
Выбор метода
интерполяции
Определение геометрии сетки

Выбор размера ячейки сетки
Выбор плотности сети следует производить в соответствии с
исходными данными или требуемым масштабом карты.
Если известен масштаб, в котором надо изобразить карту, то шаг
между линиями сетки надо задать равным тому количеству единиц
карты, которые помещаются в 1 мм изображения.
Например, при масштабе 1:50 000 это 50 м.
Если требуемый масштаб заранее не известен, то шаг между линиями
сетки можно задать равным половине среднего расстояния
между точками данных.

Gridding - методы

Inverse Distance (Обратно взвешенные расстояния),
Kriging (Крикинг),
Minimum Curvature (Минимальная кривизна),
Polynomial Regression (Полиноминальная регрессия),
Triangulation with Linear Interpolation (Триангуляция с
линейной интерполяцией),
Nearest Neighbor (Ближайший сосед),
Shepard"s Method (Метод Шепарда),
Radial Basis Functions (Радиальные вазисные функции),
Moving Average (Скользящее среднее) и т.д.

ИНТЕРПОЛЯЦИЯ:
метод Triangulation with Linear
Interpolation
Метод Триангуляция с линейной интерполяцией (Triangulation with
Linear Interpolation) базируется на триангуляции Делоне по входным точкам и
линейной интерполяцией отметок поверхности в пределах плоских граней.
z
точка с неизвестным
значениям (узел)
x
y
Триангуляция Делоне
точки с известными
значениями

ИНТЕРПОЛЯЦИЯ: метод Inverse Distance to a Power (IDW)
Метод Обратно взвешенные расстояния (Inverse Distance to a Power)
рассчитывает значения ячеек путем усреднения значений в опорных точках,
находящихся в окрестности каждой ячейки. Чем ближе точка к центру ячейки,
значение которой вычисляется, тем большее влияние, или вес, она имеет в
процессе усреднения
7,5
11,8
,
100 м
где
150 м
60 м
3,0
i – вес измеренного значения;
k – показатель степени
?
70 м
21,6
точки с известными
значениями
?
точки с неизвестными
значениями
Радиус
интерполяции

ИНТЕРПОЛЯЦИЯ: метод Minimum Curvature
Метод Минимальная кривизна (Minimum Curvature) рассчитывает значения с
использованием математической функции, которая минимизирует общую
кривизну поверхности и строит сглаженную поверхность проходящую через
опорные точки

Интерполяция: метод Polynomial Regression
Метод Полиноминальная регрессия (Polynomial Regression) основан на
аппроксимации поверхности полиномом определенного порядка:
z(х)=a0+a1x1+a2x2+…..+anxn - полином n-го порядка
Метод наименьших квадратов минимизирует сумму
- рассчитанное (оценочное) значение параметра z
- наблюденное значение параметра z

первого порядка
Аппроксимация поверхности полиномом
второго порядка

Интерполяция: метод Kriging
Метод Кригинг (Kriging) базируется на статистических моделях, которые
учитывают пространственную автокорреляцию (статистическую взаимосвязь
между опорными точками)
Случайные, но пространственнокоррелированные флуктуации
высот
Случайный шум
(валуны)
Дрейф (общий тренд
изменения высоты)
Иллюстрация элементов кригинга. Дрейф (общая тенденция), случайные, но
пространственно коррелированные колебания высоты (небольшие отклонения от общей
тенденции), и случайный шум.

Вариограмма
Полудисперсия (расстояние h) = 0.5 * среднее[ (значение в точке i – значение в точке j)2]
для всех пар точек, разделенных расстоянием h
Полудисперсия
h
h
Расстояние (лаг)
Полудисперсия
Образование пар точек:
красная точка образует пары со всеми
другими точками измерений
Остаточная
дисперсия
(nugget)
Предельный
радиус
корреляции
(range)
Расстояние (лаг)

Моделирование вариограммы
Полудисперсия
Полудисперсия
Расстояние (лаг)
Сферическая модель
Расстояние (лаг)
Полудисперсия
Экспоненциальная модель
Расстояние (лаг)
Линейная модель

Вычисление значений в узлах сети
7,5
11,8
точки с известными
значениями
100 м
150 м
60 м
3,0
?
точки с неизвестными
значениями
?
70 м
21,6
i – вес измеренного значения,
вычисляется
на
основе
модели
вариограммы
и
пространственного
распределения точек замеров вокруг
оцениваемой точки
Радиус
интерполяции

Сравнение методов интерполяции
Обратно
взвешенные
расстояния
Триангуляция с
линейной
интерполяцией
Минимальная
кривизна
Кригинг

Дополнительные опции
IV
R2
1. Определение области исходных данных для расчета значений в узлах
grid файла
I
R1
III
II

2. Учет «следов разломов» (Breaklines) и разломов (Faults)
Faults
С помощью задания Faults имитируется положение
разрывных нарушений типа сброс/взброс.
Структура файла [.BLN]
Количество точек
задания объекта
Код
(0-обнуление грида вне
контура,
1- обнуление грида
внутри контура)
X1
Y1
X2
Y2
X3
Y3
Xn
Yn
Задание Fault
Учет Faults поддерживают методы интерполяции: Inverse Distance to a
Power, Minimum Curvature, Nearest Neighbor, and Data Metrics.

Breaklines
Структура файла [.BLN]
Количество
точек
задания
объекта
Код
(0-обнуление грида
вне контура,
1- обнуление грида
внутри
контура)
X1
Y1
Z1
X2
Y2
Z2
X3
Y3
Z3
Xn
Yn
Zn
Задание Breakline
Учет Breakline поддерживают методы интерполяции:
Inverse Distance to a Power, Kriging, Minimum Curvature,
Nearest Neighbor, Radial Basis Function, Moving Average, Local
Polynomial

Учет разрывных нарушений

Учет
Breaklines
Контурная карта без
учета разломов
Учет
Faults

Визуализация изображений поверхности

Контурная карта
Карта основы
Карта точечных данных
Растр
Затененный рельеф
Векторная карта
Трехмерная сетка
Трехмерная поверхность
Результат построения сохраняется в виде векторной
графики в файле [.srf].

Контурные карты
Contour Maps

Трехмерные
изображения
поверхности
3D Surface Maps

Трехмерные сетки
3D Wireframe Maps

Векторные карты
Vector Maps

Растры
Image Maps

Карта
затененного рельефа
Shaded Relief Maps

Карты основы
Base Maps
Импортируемые форматы:
AN?, BLN, BMP, BNA, BW, DCM, DIC,
DDF, DLG, DXF, E00, ECW, EMF, GIF,
GSB, GSI, JPEG, JPG, LGO, LGS, MIF,
PCX, PLT, PLY, PNG,
PNM/PPM/PGM/PBM, RAS, RGB,
RGBA, SHP, SID, SUN, TGA, TIF, TIFF,
VTK, WMF, X, XIMG

Карты водоразделов
Watershed Maps
депрессии
водные потоки
бассейны
Карты отражают дренажные системы

Моделирование дискретных объектов

Данные XYZ
(BLN, BNA, CSV, DAT, DBF, MDB, SLK, TXT, WKx, WRx, XLS, XLSX)

Карты точечных данных (Post Maps)

Карты классифицированных точечных данных
Classed Post Maps

Boundary Files [.bln]
Количество точек
задания объекта
Код
(0-обнуление грида вне контура,
1- обнуление грида внутри контура)
X1
Y1
X2
Y2
X3
Y3
Полигон (замкнутый)
X5 ,Y5
X3 ,Y3
X4 ,Y4
X2 ,Y2
Xn
X6 ,Y6
Yn
X10 ,Y10
X1 ,Y1
Линия
X6 ,Y6
X7 ,Y7
X4 ,Y4
X2,Y2
X5 ,Y5
X3 ,Y3
X1 ,Y1
X7 ,Y7
X8 ,Y8
X9 ,Y9
X1=X10
Y1=Y10

Вычисление ошибок интерполяции,
Графическое редактирование грида.

Ручная коррекция сетки (Grid Node Editor)

Графический редактор для ввода и коррекции значений данных
сеточной области

Оценка точности интерполяции (Residuals)

Пункт меню Grid

Математические операции над гридами (Math)
Входной грид 1
Позволяет осуществлять
вычисления над одним или
двумя гридами
Входной грид 2
Выходной грид
Формула вычислений
-
Кровля
=
Подошва
Мощность

Анализ поверхностей (Calculus)
Методы
Позволяет осуществлять анализ
формы поверхности
Входной грид
Выходной грид
Углы
наклона
Terrain
Slope
Ориентация
склонов
Terrain Aspect

Фильтрация (Filter)
Входной грид
Выходной грид
Размер
оператора
Методы
Позволяет выделять
разночастотные составляющие
моделей поверхностей
Оператор
Низкочастотная
фильтрация
41 41

Бланкирование (Blank)
Позволяет обнулять участки карты, определенные файлом [.bln]
Входной грид
+ Файл [.bln] = Выходной грид
Бланкировка
Blank
Границы полигона

Построение разрезов (Slice)
Позволяет рассечь поверхность вдоль линии, положение
которой определено файлом [.bln]
Входной грид
+ Файл [.bln] = Выходной файл [.dat]
X
Y
Z
Расстояние
по профилю
Линия профиля
64
Разрез по профилю
Z
56
48
40
0
20000
40000
Расстояние по профилю
60000
80000

Михаил Владимирович Морозов:
персональный сайт

Мат.модели (занятие, карта-1): Построение геохимических карт в Golden Software Surfer (общий подход, этапы и содержание работы, форма отчета)

Курс "Математические методы моделирования в геологии "

Карты-1. Построение геохимических карт в Golden Software Surfer: общий подход, этапы и содержание работы. Форма отчета.
Карты-2. Принципы работы с Golden Software Surfer.

Чтобы найти место скопления полезного металла в земной коре, требуется геохимическая карта. Как ее построить? Для этого необходимы хорошее программное обеспечение и системный подход. Познакомимся с принципами и основными этапами этой работы.

ТЕОРИЯ

Построение геохимической карты в программе Golden Software Surfer.

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

Координаты : в программе Surfer мы используем прямоугольные координаты (в метрах) , хотя в свойствах карты можно выбрать среди возможных систем координат также и различные полярные координаты (в градусах-минутах-секундах). На практике при работе с изображениями на плоском листе бумаге удобнее работать в системе прямоугольных координат в пользовательском формате.

Откуда берутся координаты:
1. При документации точки на месте координаты берутся из топопривязчика GPS или ГЛОНАСС в виде полярных координат (например, в системе координат WGS 84 ). Топопривязчик может нынче иметь вид смартфона, но удобнее и надежнее использовать специальный прибор, который ласково называют "джипиэской".
2. При переносе данных на компьютер из топопривязчика, координаты преобразуются из полярных в используемую систему прямоугольных координат (например, в системах UTM , Пулково-1942 , но можно использовать и местную геодезическую систему, принятую на конкретном предприятии). Для преобразования полярных координат в прямоугольные удобно использовать программу Ozi Explorer .
3. В столбцах электронной таблицы, подготовленной для работы с Surfer, должны располагаться прямоугольные координаты в метрах.

Картируемая величина : для построения учебной карты в изолиниях мы будем использовать логарифм содержания какого-либо химического элемента. Почему логарифм? Потому что закон распределения содержаний микроэлементов почти всегда логарифмический. Разумеется, в реальной работе сперва нужно проверить закон распределения , чтобы выбрать вид величины: исходное значение или его логарифм.

Виды карт, используемые в геохимии . Помимо карты в изолиниях геохимики часто используют некоторые другие типы карт, но не все великое разнообразие типов карт, которые умеет строить Surfer, а только строго определенные. Они перечислены ниже.

1. Карта фактов. Представляет собой набор точек, показывающих места опробования на местности. Около точек можно выводить метки - номера пикетов, но при геохимических поисках точек так много, что обычно метки лишь "засоряют" пространство карты и не приводятся. Для построения карты фактов используем функцию Post Map .

2. Точечная карта содержаний химического элемента. На ней кружками (или другими символами) разных размеров обозначаются разные содержания химического элемента в точках опробования. Если мы используем такую карту, то отдельная карта фактов уже не нужна - точки обеих карт наложатся друг на друга. Точечная карта (или "карта-разноска") строится так, чтобы высокие содержания искомого элемента бросались в глаза. На легенде обозначается соответствие между размером кружка и содержанием элемента в г/т. Помимо размера может изменяться цвет кружка. Каждому типу (размеру, цвету) кружка соответствует назначенный вручную диапазон содержаний. Т.е. разные типы кружков - это разные классы точек по содержаниям элемента. Поэтому инструмент для создания такой карты называется Classed Post Map . Удобно строить карту-разноску поверх карты в изолиниях, чтобы видеть, как последняя (которая является расчетной картой, т.е. построена по результатам интерполяции данных) сочетается с исходными, полученными из лаборатории, т.е. "истинными" содержаниями. Удобно наносить разноску одного важного элемента (например, золота) на карту в изолиниях другого поискового параметра (элемента-спутника, статистического фактора, геофизического параметра и т.п.). Важно: после построения карту типа Classed Post Map нельзя преобразовать в Post Map, наоборот тоже нельзя.

3. Карта в изолиниях. Собственно карта искомого параметра, где разные градации содержаний отображены разными цветными заливками. Также требует легенды, которая связывает цвет заливки с уровнем содержаний. Градации заливок настраиваются вручную. Инструмент - Contour Map . Помимо собственно содержаний элементов (или их логарифмов) в геохимии широко используются карты многоэлементных показателей. Это могут быть мультипликативные коэффициенты (где содержания нескольких элементов перемножаются), карты значений фактора (главной компоненты) и т.п. Собственно, задача геохимика - найти показатель, который позволяет решить геологическую задачу. Коль скоро такие показатели, как правило, выражаются в коллективном поведении элементов, вполне естественно, что моноэлементные карты (т.е. карты одного отдельно взятого элемента) часто менее информативны, чем полиэлементные. Поэтому этап построения карт обычно предваряется этапом статистической обработки данных с получением результатов многомерного статистического анализа, например, МГК (метода главных компонент).

4. Обводка карты. По умолчанию Surfer создает прямоугольную карту. В том случае, если точки опробования не образуют прямоугольник, получается, что область опробования вписана в искусственно созданный прямоугольник, в котором часть площади в реальности не опробовалась. Карта в изолиниях будет построена на всю площадь, поэтому неопробованные участки карты будут содержать фиктивные данные. Чтобы избежать этого, нужно ограничить область построения карты той частью площади, на которую имеются данные опробования. Для этого область опробования нужно оконтурить специальной линией, которая может быть построена вручную. Вывод контура обводки осуществляется посредством функции Base Map .

Этапы построения карты.

3. Построение карты фактов [карта-3]. 5. Построение точечной карты ("карты-разноски") [карта-5]. 9. Построение карты поверхности и ее оформление для достижения оптимальной информативности [карта-6, продолжение].

ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

Дано : таблица содержаний химического элемента и его логарифмов с координатами точек опробования.

Задание :

1. Построить карту фактов .

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

3. Самостоятельно создать контур площади картирования и построить его .

4. Совместить контур площади, точечную карту элемента и карту фактов в данном порядке в менеджере объектов. Вывести легенду для точечной карты.

5. Построить сеточный файл ("грид") для логарифмов содержаний элемента методом триангуляции , проверить его . Повторить другими методами.

6. Построить вариограмму для построения сеточного файла методом крайгинга , проверить его .

7. Построить сеточный файл ("грид") для логарифмов содержаний элемента методом крайгинга с использованием параметров вариограммы.

8. Сгладить полученный сеточный файл простым фильтром .

9. Восстановить сеточный файл из логарифмов в содержания .

10. Обрезать сеточный файл по созданному ранее контуру .

11. Построить карты поверхности в изолиниях и градиентной заливке по созданным сеточным файлам, добавить легенды.

12. Экспортировать построенные карты как файлы JPG, вставить в отчет в формате Word (DOC).

Форма отчета.

Построение регулярной ЦМР в Golden Software Surfer 8

8.1 Построения регулярных координатных сеток

Основой представления пространственной информации как в векторной, так и в растровой модели данных, являются дискретные операционно- территориальные единицы (ОТЕ). Каждая ОТЕ является пространственным объектом, для которого предполагается однородность имеющейся о нем атрибутивной информации с точки зрения изучаемого явления.

В растровой модели данных ОТЕ соответствуют ячейкам регулярной или нерегулярной сетки, которая полностью покрывает всю территорию характерного размера представляемых пространственных объектов и явлений. Наиболее популярным является использование регулярных прямоугольных сеток, ячейки которых в картографической проекции являются квадратами или прямоугольниками.

В пакете геостатистического моделирования Surfer создание модели, в качестве ОТЕ которой используется регулярная сетка, осуществляется через пункт меню «Grid», подпункт «Data». После нажатия на него появляется диалоговое окно, предлагающее выбрать файл с исходными данными.

После выбора исходного файла появляется диалоговое окно, где мы можем задать все необходимые для построения растровой модели параметры (рис. 8.1). Данное окно имеет 4 закладки. Закладка «General» является главной и имеет несколько строк ввода, кнопок и переключателей, объединенных в различные группы. Первая группа «Grid Line Geometry» позволяет задать координаты левого нижнего угла в столбце «Minimum», координаты правого верхнего угла в столбце «Maximum», шаг сетки (на каком расстоянии будут отстоять друг от друга узлы растровой модели по Х и по Y) в столбце «Spacing». Здесь же в столбце «# of Lines» вычисляется количество узлов.

Ниже расположена группа «Gridding Method», состоящая из выпадающего списка в левой части, где приведены названия методов, которыми можно вычислить значения в узлах нашей сетки. Все приведенные методы можно разделить на интерполяционные, аппроксимационные и методы геостатистики. При нажатии на кнопку «Options» отображаются параметры выбранного метода расчета. Более подробно о методах и параметрах, которые в них используются, можно прочитать в дополнителььной литературе .

На второй закладке «Data» (рис. 8.3) можно задать, какие данные из исходного файла будут использованы для расчета значений в узлах растровой модели (в группе выпадающих строк «Data Columns» необходимо задать, в какой колонке находятся координаты Х, Y и переменная Z).

Рис. 8.3

На этой же закладке (рис. 8.3) в строке «Data Exclusion Filter» можно задать условие, необходимое для исключения части данных из расчетов (исключены будут данные, удовлетворяющие условию). На той же закладке данного окна (рис. 8.3) приведены некоторые статистические показатели трех переменных X, Y, Z (общее количество записей, количество исключенных, среднее, минимум, максимум и т.д.). При нажатии на кнопку «Update Statistics» список статистик обновляется с учетом того условия, которое было введено в строке «Data Exclusion Filter» (если ничего не было введено, то список останется тем же.)

На закладке «Data» в верхнем правом углу есть группа операторов «Duplicates», которые определяют порядок действий при обнаружении дублирующих данных, а также определяют, что является дублированием. Выпадающий список «To Keep» определяет, какая из одинаковых по положению в пространстве точек будет взята (по умолчанию первая в исходном файле). Строки «X Tolerance» и «Y Tolerance» задают расстояние, на которое может отстоять точка от точки, чтобы считаться дублирующими.

Рис. 4

Закладка «Search» (рис.8.4) доступна не всегда. Данная закладка доступна, если на закладке «Data» в списке «Gridding Method» выбраны такие методы расчета как: «Radial Basic Function», «Inverse Distance to Power», «Kriging», «Nearest Neighbor». Эта закладка позволяет задать, какие исходные данные будут участвовать в вычислении значения в текущем узле сетки. Для определения данных, которые будут использованы, все пространство вокруг узла разделяется на 4 равных сектора и ограничено эллипсом. По умолчанию для вычисления значения в текущем узле сетки используются все исходные данные, при этом в верхнем левом углу стоит галка в переключателе «No Search (use all of the data)», а все остальные операторы данного окна недоступны. Если эту галочку убрать, то становятся доступными остальные элементы закладки, задающие поиск. Строка ввода «Number of sectors to search» необходима для задания количества секторов, данные из которых будут участвовать в расчете значения в узле сетки. Строка «Maximum number of data to use from each sectors» определяет максимальное количество точек в каждом секторе, значения из которых будут использоваться для расчетов в узлах сетки. Строка «Minimum number of data in all sectors (node is blanked if fewer)» задает минимальное количество точек во всех секторах, которое необходимо для расчета значения в текущем узле модели (узел будет бланкирован, если этих точек меньше). Строка «Blank node if more than this many sectors are empty» нужна для определения минимально необходимого количества непустых секторов для расчета значений в узлах модели. Строки «Radius 1» и «Radius 2» задают радиусы поискового эллипса, а строка «Angle» - его угол наклона.

где n - количество точек, образующих ломаную линию, X и Y - это координаты этих точек в той же системе координат, что и исходные точки. При использовании линий типа «Fault» расстояние от точки с исходными данными до узла сетки вычисляется, как показано на рис. 8.6.

Рис. 8.6

где n - количество точек, образующих ломаную линию, X, Y и Z - это координаты этих точек в той же системе координат, что и исходные точки. При использовании линий типа «Breaklines» расстояние от точки с исходными данными до узла сетки вычисляется так, как показано на рис. 8.7. Забегая вперед, заметим, что эти линии можно подготовить непосредственно в программе «Surfer», используя пункт главного меню «Map», подпункт «Digitze». «Breaklines» , в отличие от «Fault», не являются барьерами, это линии с дополнительной информацией, которая задействуется при использовании точек данных, лежащих по отношению к узлу сетки за этими линиями. Значение, которое при этом используется с «Breaklines» , рассчитывается методом линейной интерполяции между двумя вершинами линии.

После того, как мы задали все необходимые параметры, на данных четырех закладках необходимо нажать кнопку «Ок», после чего начнется создание модели. Здесь стоит отметить, что по тем или иным причинам значение в узле модели может быть не рассчитано, тогда туда записывается код незаданного значения = 1,70141 * 1038.

Рис. 7

8.2. Картографирование полученных моделей

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

Создание карт реализовано в пункте меню «Мар».

Карты изолиний реализованы в подпункте меню «Contour Map». Данный пункт меню содержит три подпункта «New Contour Label», «Edit Labels...», «Export Contours.». Непосредственно создание карт изолиний осуществляется через пункт «New Contour Label». Пункт меню «Edit Labels.» отвечает за ручную корректировку положения подписей изолиний. Через «Export Contours.» осуществляется экспорт изолиний в *.dxf формат.

В группе операторов «Filled Contours» есть два переключателя: 1 - «Fill Contours» отвечает за включение заливки между изолиниями, 2 - «Color Scale» переключатель, доступ к которому появляется при включенном «Fill Contours». Второй переключатель разрешает или запрещает вывод цветовой шкалы заливки на экран. Группа переключателей «Smothing» отвечает за сглаживание изолиний (сама сетка не сглаживается). При включенном «Smooth Contours» сглаживание начинает работать. Степень сглаживания задается чуть ниже в выпадающем списке «Amount». Не следует особенно увлекаться использованием подобного сглаживания, так как при большом количестве ОТЕ в модели данная функция для своей работы требует больших аппаратных ресурсов. Группа «Blanked

Regions» используется для определения того, каким образом будут показаны на карте узлы растровой модели, в которых записан код не заданного значения. Кнопка «Fault Line Properties» отвечает за задание типа, толщины и цвета линий типа «Fault» (если сетка строилась с их учетом).

На закладке «Levels» (рис. 8.9) задаются следующие параметры изолинейной карты: Level - шаг горизонталей; Line - тип и цвет изолинии определенного уровня; Fill - цвет и тип заливки пространств между изолиниями; Label - свойства подписей изолиний (интервал подписи изолиний, фонт и т.д.); Hach - свойства бергштрихов (штрихи - перпендикулярные изолинии, показывающие направление убывания или возрастания функции).

Рис. 8.9

Кроме всего прочего, на данной закладке справа есть четыре кнопки: 1 - «Add» позволяет задать дополнительную горизонталь с отметкой, заданной исследователем; «Delete» - кнопка, удаляющая любую горизонталь; «Save» - позволяет сохранить пользовательский набор горизонталей в файл формата *.lvl и загрузить его, используя клавишу «Load» для другой сетки.

Метод цветовых шкал реализован в подпункте меню «Image map». При выборе этого пункта меню появляется окно, предлагающее выбрать *.grd файл, по которому будет построена карта. После выбора нужного файла появляется диалоговое окно «Image Map Properties» (рис. 8.10). В его верхней части отображается строка с именем того *.grd файла, который мы выбрали.

Ниже слева есть два доступных переключателя: 1 - «Interpolate Pixels» осуществляет сглаживание пиксельной структуры растровой картинки (если включен, то пиксельная структура меньше просматривается); 2 - «Show Color Scale» отображает цветовую шкалу на карту.

Рис. 8.10

По центру окна «Image Map Properties» под названием «Colors» расположена кнопка, нажав на которую можно получить окошко «Color Spectrum» (рис. 8.11). Используя это окно, можно задать цвета для всего диапазона значений. Используя указатель курсора мыши и клавишу «Ctrl», можно добавить в шкалу необходимые для отображения цвета. Созданный цветовой набор можно сохранять, а затем загружать для отображения других сеток. Для этого используются клавиши «Load.», «Save...». В нижней левой части экрана есть группа операторов «Data to Color Mapping», которая предназначена для определения диапазона значений, отображаемых заданными цветами. После определения всех необходимых параметров необходимо нажать на кнопку «Ok», и карта будет построена.

Рис. 8.11

Метод аналитической отмывки реализован в том же пункте меню «Map», подпункте «Shaded relief Map». При выборе данного пункта появляется окно, предлагающее определиться с *.grd файлом, на основе которого будет построена карта. Далее появляется окно «Shaded Relief Map Properties» (рис. 8.12), где предлагается задать параметры данного вида отображения. В верхней части экрана традиционно расположено название файла сеточной модели. Ниже и левее есть блок «Light Position Angles», задающий в градусах положение источника света по вертикали - «Vertical» и горизонтали - «Horizontal». В блоке «Relief Parameters» расположены различные отображения рельефа. «Gradient Method» задает степень сглаженности картинки: «Central Difference» - более сглаженная картина; «Midpoint Difference» - менее сглаженная картина. «Shading Method» задает алгоритм расчета отраженного света от каждой ячейки. В строке «Z Scale Factor» задается коэффициент масштабирования по Z (чем он больше, тем длиннее тени). В блоке «Colors», как нетрудно догадаться, задаются цвета, которыми будут окрашиваться затененные и освещенные участки поверхности.

В блоке «Missing Data» задается цвет узлов сетки с кодами незаданных значений. После задания всех параметров необходимо нажать на кнопку «Ок» и посмотреть на результат (полученную карту).

Кроме трех вышеперечисленных карт поверхностей, можно получить карту постов исходных наблюдений. Такая возможность реализована в пункте меню «Post Map», подпункте «New Post Map». Выбирая данный пункт меню, получаем диалоговое окно «Post Map Properties» (рис. 8.13), которое имеет две закладки. Первая - «General» (рис.8.13), на которой указываются общие свойства карты постов. В верхнем левом углу название файла постов, ниже в группе «Worksheet Columns» расположены четыре выпадающих списка, которые определяют, в какой колонке исходного файла располагаются координаты Х, Y (X Coord и Y Coord), где расположен номер символа, которым будет отображена точка (Symbol) и угол (Angle) наклона данного символа. В блоке «Default Symbol» задается то, каким образом будут отображены посты, если в строках «Symbol» и «Angle» ничего не задано (значение «None»). Снизу справа в блоке «Symbol Size» задается размер поста. Возможны два варианта: 1 - «Fixed Size» - фиксированный размер (и указывается его размер); 2 - «Proportional» - размер поста на карте выставляется пропорционально значению какой-либо колонки из исходного файла.

Закладка «Labels» (рис. 8.14) отвечает за оформление подписей к карте постов. Здесь в выпадающем списке «Worksheet Column for Labels» приведен список названий колонок имен файлов, из которых будут браться подписи. Необходимо выбрать одно из них. Если ничего не выбрано, то подписи точки не будет.

Рис. 8.12

Рис. 8.13

Кроме всего прочего, в этом окне в блоке «3D Label Lines» можно задать, каким образом будет указываться точка на трехмерной поверхности (например, на карте, созданной через пункт меню «Wareframe...»). Группа операторов, объединенная под названием «Position Relative to Symbol», необходима для определения положения подписи относительно точки поста: Above - выше; Below - ниже; Center - по центру; Left - слева; Right - справа. При выборе из этого списка пункта «User Defined» появляется доступ к строкам «X Offset», «Y Offset», где можно задать пользовательское положение подписи. Строка «Label Angle» - угол наклона подписи. В правом верхнем углу есть еще две кнопки: «Font.» - задает свойства шрифтов подписи; «Format...» - задает количество знаков, отводимых на написание подписи, в том числе после запятой.

Рис. 8.14

Кроме отображения карт по отдельности, их можно накладывать друг на друга и показывать совмещенно (операция «Overlay»). Для этого необходимо выбрать те карты, которые мы хотим наложить (используя Ctrl и указывая нужную карту или нажав клавишу F2 на клавиатуре, при этом выберутся все карты), выбрать пункт меню «Map» ^ «Overlay Maps». Для карт в Surfer можно задавать масштаб. Для этого необходимо выбрать нужную карту указателем курсора мыши и нажать левую клавишу мыши, а затем выбрать пункт главного меню «Map», подпункт «Scale». После этих действий появляется диалоговое окно «Scale» (рис. 8.15), где можно задать масштабы для каждой координаты.

В данном окне можно установить пропорциональные масштабы для координат Х и Y, если включить опцию «Proportional XY Scaling». Когда мы выставляем масштаб для трехмерных отображений поверхности, то нам становится доступной последняя третья строка «Z Scale», в которой можно задать масштаб по координате Z.

Рис. 8.15

8.3. Операции с поверхностями

Помимо создания поверхностей с использованием вышеперечисленных методов с поверхностями можно совершать различные математические операции (сложение, вычитание, умножение и т.д.). Можно изменять координаты регулярных сеток, рассчитывать различные показатели, основанные на вычислении производной первого и второго порядка и использовании «скользящего окна» разной формы и размера (сглаживание).

Математические операции с сетками можно осуществлять, используя пункт меню «Grid», подпункт «Math...». После выбора этого пункта меню появляется диалоговое окно, предлагающее выбрать *.grd файл, с которым будут совершены определенные действия. Затем появляется окно «Grid Math» (рис. 8.16). Здесь в первой строчке под названием «Input Grid File A» указывается путь и имя файла, который мы указали в предыдущем диалоговом окне. Стоит заметить, что можно осуществить операцию с двумя моделями (например, сложить), тогда в строке «Input Grid File В», задаем еще одну поверхность. Модели А и В должны иметь равное количество строк и колонок (количество узлов по Х и Y в обеих моделях должно быть одинаковое). Вторая модель может и не использоваться, если мы хотим задействовать константу для вычисления новой модели. Тогда в последней строке вместо литеры В используется значение этой константы. В третьей строке «Output Grid File С» задается имя результирующей сетки.

Рис. 8.16

Рис. 8.17

Общий смысл операции над двумя моделями представлен на рис.8.17. Фактически при этом складываются значения Z1 и Z2 из узлов моделей А и В, а результат Z3 записывается в узел модели С, которая аналогична по размерам.

Используя пункт главного меню «Grid», подпункт «Calculus», мы можем рассчитать модели некоторых показателей, в основу которых заложено вычисление производных растровой модели. Более подробно будет рассмотрено вычисление производных рельефа (угла наклона, экспозиции, профильной и плановой кривизны) или морфометрических показателей.

После выбора соответствующего пункта меню и файла модели появляется диалоговое окно «Grid Calculus» (рис. 8.18).

Рис. 8.18

В строке «Input Grid File» приведено имя файла с исходной моделью. В последней строке «Output Grid File» указывается имя файла, в котором будут записаны результаты расчетов. Мы можем изменить имя файла и место его расположения. В выпадающем списке «Grid- to-Grid Operation» приведено 4 строки. Для расчета морфометрических характеристик необходимо выбрать строку «Terrain Modeling». После этого нужно нажать на кнопку «Options...», в результате появится новое диалоговое окно «Terrain Modeling Options» (рис. 8.19), где мы можем выбрать, какой морфометрический показатель необходимо рассчитать («Terrain Slope» - уклон, «Terrain Aspect» - экспозиция, «Profile Curvature» - профильная кривизна, «Plan Curvature» - плановая кривизна, «Tangential Curvature» - тангенциальная кривизна).

Рис. 8.19

Кроме расчета показателей, основанных на вычислении производной, в Surfer есть возможность произвести так называемое матричное сглаживание. Матричное сглаживание в Surfer представляет собой операцию по вычислению нового значения в узле модели, используя расчет среднего значения или вычисление средневзвешенного значения (вес при этом обратно пропорционален расстоянию). Для расчета среднего и средневзвешенного используются значения узлов растровой модели из локальной окрестности. Размер этой окрестности и количество узлов, участвующих в вычислениях, задается скользящим окном. Размер стороны скользящего окна в Surfer может изменяться от 3 до n, где n - размер одной из сторон модели. В самом простом случае узлы модели для вычисления нового значения используются те, что указаны на рис. 8.20. Новое значение вычисляется по формуле (1) .

Традиционно в верхней части окна расположена строка «Input Grid File», в которой записано имя файла с исходной моделью. Ниже слева в блоке «Method» можно выбрать, какое из значений (Average - среднее, Weighted - средневзвешенное) будет рассчитано в локальном окне. Правее ниже расположены 4 строки. Строка «Weight of Matrix Center» предназначена для определения веса центрального значения в матрице. При вычислении по формуле (1) - это параметр n. Строка «Rows on Either Side of Center» задает параметр n, который необходим для вычисления размера скользящего окна по Y. Строка «Rows on Either Side of Center» задает параметр n, который необходим для вычисления размера скользящего окна по X. Размер одной из сторон вычисляется по формуле (2 * n) + 1. Строка «Distance Weighting Power» становится доступной только в том случае, если выбран метод «Weighted». Эта строка задает степень влияния на вычисляемое значение расстояния от узла сетки до центра скользящего окна.

В строке «Output Grid File» задается имя файла с результирующей моделью.

В Surfer есть возможность производить трансформацию системы координат растровой модели. Данная возможность реализована в пункте меню «Grid» - «Transform». Выбирая различные строки в выпадающем списке «Operation», мы можем произвести смещение координат по Х и Y (Offset), масштабировать по координатам осям (Scale), повернуть на определенный угол (Rotate), зеркально отобразить (Mirror X, Mirror Y) (рис. 8.22). Правее этого списка приведены строки, где мы можем указать то, на сколько смещать или поворачивать координаты. В строках «Input Grid File» и «Output Grid File» приведены имена файлов, где записаны исходные данные и результирующая модель.

Рис.8.22