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

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

Математические операции с сетками можно осуществлять, используя пункт меню «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).

 

 

Рис. 8.20

 

В Surfer эти расчеты реализованы в пункте главного меню «Grid», подпункте «Matrix Smooth...». После выбора этого пункта меню и задания файла с исходной моделью появляется диалоговое окно «Matrix Smooth» (рис. 8.21).

Традиционно в верхней части окна расположена строка «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» задается имя файла с результирующей моделью.