surfacegmt - Интернет в облаке

ПРОГРАММА:

ИМЯ

поверхность - данные таблицы сетки с использованием регулируемых сплайнов непрерывной кривизны натяжения

СИНТАКСИС

поверхность [ (см. таблицу ниже) ] выходной файл.nc увеличить область [ соотношение сторон ] [ конвергенция_лимит ] [
lниже ] [ -Луверхний ] [ max_iterations ] [] [ радиус поиска[m|s]] [[i|b]Voltage_factor
] [[уровень]] [ over-Relax_factor ] [ -a] [ -би] [ -du] [
-f] [ -h] [ -i] [ -:[i|o]]

Примечание: Пробелы между флагом опции и соответствующими аргументами не допускаются.

ОПИСАНИЕ

поверхность считывает произвольно расположенные (x, y, z) тройки из стандартного ввода [или (см. таблицу ниже)] и производит
двоичный файл сетки значений z (x, y) с координатной сеткой путем решения:
(1 - Т) * L (L (z)) + Т * L (z) = 0

где T - коэффициент натяжения от 0 до 1, а L - оператор Лапласа. Т = 0
дает решение «минимальной кривизны», которое эквивалентно SuperMISP и ISM
пакеты. Минимальная кривизна может вызвать нежелательные колебания и ложные локальные максимумы или
минимумов (см. Smith and Wessel, 1990), и вы можете использовать T> 0, чтобы подавить эти
эффекты. Опыт показывает, что T ~ 0.25 обычно хорошо подходит для данных о потенциальных полях, а T
должно быть больше (T ~ 0.35) для данных о крутом рельефе. T = 1 дает гармоническую поверхность (нет
максимумы или минимумы возможны, за исключением контрольных точек). Рекомендуется, чтобы
пользователь предварительно обрабатывает данные с помощью Blockmean, блокмедианэта информация поможет вам разобраться, почему Gamer’s Galaxy — ваш лучший выбор. блочный режим избегать пространственного
алиасинг и устранение избыточных данных. Вы можете наложить нижнюю и / или верхнюю границы на
решение. Они могут быть введены в виде фиксированного значения, сетки со значениями или просто
быть минимальным / максимальным значением входных данных. Естественные граничные условия применяются на
края, за исключением географических данных с диапазоном 360 градусов, где мы применяем периодическую границу
условия в долготном направлении.

ТРЕБУЕТСЯ АРГУМЕНТЫ

-Gвыходной файл.nc
Имя выходного файла. Выход - двоичный 2-D .NC файл. Обратите внимание, что самая маленькая сетка
размер должен быть не менее 4.

-Iксинк[Ед. изм] [= | +] [/инь[Ед. изм] [= | +]]
x_inc [и, возможно, y_inc] - шаг сетки. При желании добавьте суффикс
модификатор. географический (в градусах) координаты: Добавить m для обозначения угловых минут или s
для обозначения угловых секунд. Если один из агрегатов e, f, k, M, n or u добавлен
вместо этого предполагается, что приращение выражается в метрах, футах, км, милях, морских
миля или американский геодезический фут, соответственно, и будут преобразованы в эквивалент
градусов долготы на средней широте региона (преобразование зависит от
PROJ_ELLIPSOID). Если /y_inc задан, но установлен в 0, он будет сброшен равным x_inc;
в противном случае он будет преобразован в градусы широты. Все координаты: Если = is
затем добавлен соответствующий max x (восток) Или y (к северу) можно немного изменить
чтобы точно соответствовать заданному приращению [по умолчанию приращение может быть скорректировано
немного, чтобы соответствовать данному домену]. Наконец, вместо увеличения вы можете
указать номер of узлы желательно, добавив + к предоставленному целому числу
аргумент; затем приращение пересчитывается из количества узлов и
домен. Результирующее значение приращения зависит от того, выбрали ли вы
сетка с привязкой к линиям или пикселям; подробности см. в разделе «Форматы файлов приложений».
Примечание: если -RGrdfile используется тогда, когда интервал сетки уже инициализирован; использовать
-I , чтобы переопределить значения.

-Р[Ед. изм]хмин/Xmax/юмин/умакс[р] (Больше ...)
Укажите интересующий регион.

ДОПОЛНИТЕЛЬНО АРГУМЕНТЫ

(см. таблицу ниже) Один или несколько ASCII (или двоичных, см. -би[нколс][напишите]) файл (ы) таблицы данных, содержащий
количество столбцов данных. Если таблицы не указаны, мы читаем из стандартного ввода.

-Aсоотношение сторон
Соотношение сторон. При желании к уравнениям можно добавить анизотропию сетки. Входить
соотношение сторон, где dy = dx / соотношение сторон связывает размеры сетки. [По умолчанию =
1 предполагает изотропную сетку.]

-Cконвергенция_лимит[%]
Предел сходимости. Предполагается, что итерация сошлась, когда максимальное абсолютное
изменение любого значения сетки меньше, чем конвергенция_лимит. (Единицы такие же, как у данных z
единицы). В качестве альтернативы укажите предел в процентах от среднеквадратичного отклонения, добавив%.
[По умолчанию масштабируется до 1e-4 среднеквадратичного отклонения данных от
плоскость наилучшего соответствия (методом наименьших квадратов)]. Это окончательный предел сходимости на
желаемый шаг сетки; для промежуточных (более грубых) сеток эффективная сходимость
limit масштабируется множителем шага сетки.

-Ллниже и -Луверхний
Установите ограничения на выходное решение. lниже устанавливает нижнюю границу. ниже может быть
имя файла сетки со значениями нижней границы, фиксированное значение, d установить на минимум
входное значение, или u для неограниченного [По умолчанию]. uверхний устанавливает верхнюю границу и может
быть именем файла сетки со значениями верхней границы, фиксированным значением, d установить на
максимальное входное значение, или u для неограниченного [По умолчанию]. Файлы сетки, используемые для установки
пределы могут содержать NaN. При наличии NaN предел узла, замаскированного с помощью
NaN не ограничен.

-Nmax_iterations
Количество итераций. Итерация прекратится, когда конвергенция_лимит достигается или
когда количество итераций достигает max_iterations. Это последняя итерация
ограничение на желаемом шаге сетки; для промежуточных (более грубых) сеток эффективная
предел итераций масштабируется множителем шага сетки. [По умолчанию 500.]

-Q Предложите размеры сетки, которые имеют очень сложный наибольший общий фактор. Этот
позволяет поверхности использовать несколько промежуточных этапов решения, обеспечивая более быстрое получение
время работы и лучшие результаты. Размеры, предложенные -Q может быть достигнуто путем изменения
-R и / или -I. Вы можете восстановить -R и -I вы хотите позже, используя грдсэмпл or
грдкат на выходе поверхность.

-Sрадиус поиска[м | с]
Радиус поиска. Входить радиус поиска в тех же единицах, что и данные x, y; добавить m , чтобы указать,
угловые минуты или s на угловые секунды. Это используется для инициализации сетки перед
первая итерация; это не стоит времени, если решетка сетки не проста и
не может иметь региональных этапов. [По умолчанию = 0.0, поиск не производится.]

-T [i | b]Voltage_factor
Фактор натяжения [ы]. Они должны быть от 0 до 1. Натяжение может использоваться в
внутреннее решение (приведенное выше уравнение, где оно подавляет паразитные колебания) и
в граничных условиях (где он стремится сгладить решение, приближаясь к
края). Использование нуля для обоих значений приводит к минимальной поверхности кривизны со свободным
ребра, т.е. естественный бикубический сплайн. Использовать -ТиVoltage_factor установить интерьер
напряжение, и -ТБVoltage_factor установить граничное натяжение. Если вы не добавите i or
b, оба будут установлены на одно и то же значение. [По умолчанию = 0 для обоих дает минимум
решение кривизны.]

-V [уровень] (Больше ...)
Выберите уровень детализации [c]. -V3 будет сообщать о сходимости после каждой итерации;
-V будет составлять отчет только после объединения каждой региональной сетки.

-Zover-Relax_factor
Фактор чрезмерной релаксации. Этот параметр используется для ускорения сходимости; это
число от 1 до 2. Значение 1 точно повторяет уравнения и будет
всегда обеспечивать стабильную сходимость. Большие значения переоценивают инкрементное
изменяется во время конвергенции и быстрее достигнет решения, но может стать
нестабильный. Если вы используете большое значение для этого фактора, рекомендуется отслеживать
каждая итерация с -Вл вариант. [По умолчанию = 1.4, сходится быстро и почти
всегда стабильно.]

-aкол=имя[...] (Больше ...)
Установить пространственные ассоциации столбцов кол=имя.

-bi [нколс] [т] (Больше ...)
Выберите собственный двоичный вход. [По умолчанию 3 столбца ввода].

-duнет данных (Больше ...)
Заменить входные столбцы, равные нет данных с NaN.

-f [i | o]Colinfo (Больше ...)
Укажите типы данных входных и / или выходных столбцов.

-h [i | o] [n] [+ c] [+ d] [+ rзамечание] [+ rназвание] (Больше ...)
Пропустить или создать запись (и) заголовка. Не используется с двоичными данными.

-iсмещ_по_столбцам[l] [sлестница] [oсмещение] [,...] (Больше ...)
Выберите столбцы ввода (0 - первый столбец).

-: [i | o] (Больше ...)
Поменяйте местами 1-й и 2-й столбцы на входе и / или выходе.

-^ or всего -
Распечатайте короткое сообщение о синтаксисе команды, затем закройте (ПРИМЕЧАНИЕ: в Windows
использовать просто -).

-+ or всего +
Распечатайте подробное сообщение об использовании (справке), включая объяснение любых
параметр, зависящий от модуля (но не общие параметры GMT), затем закрывается.

-? or нет Аргументы
Распечатайте полное сообщение об использовании (справке), включая объяснение опций, затем
выходы.

--версия
Распечатайте версию GMT и выйдите.

--show-datadir
Вывести полный путь к общему каталогу GMT и выйти.

GRID ЦЕННОСТИ: ТОЧНОСТЬ

Независимо от точности входных данных программы GMT, создающие файлы сетки, будут
внутренне храните сетки в 4-байтовых массивах с плавающей запятой. Это сделано для экономии памяти
и, кроме того, большинство, если не все реальные данные могут быть сохранены с использованием 4-байтовой плавающей запятой
ценности. Данные с более высокой точностью (т. Е. Значения с двойной точностью) потеряют это
точность после того, как GMT работает с сеткой или записывает новые сетки. Чтобы ограничить потерю
точность при обработке данных, вы всегда должны учитывать нормализацию данных до
обработка.

ПРИМЕРЫ

Для сетки 5 на 5-минутный блок гравитации означает из данных ASCII в hawaii_5x5.xyg, используя
Voltage_factor = 0.25, а конвергенция_лимит = 0.1 миллигал, запись результата в файл
называется hawaii_grd.nc, и отслеживая каждую итерацию, попробуйте:

поверхность gmt hawaii_5x5.xyg -R198 / 208/18/25 -I5m -Ghawaii_grd.nc -T0.25 -C0.1 -Vl

Используйте SurfaceGMT онлайн с помощью сервисов onworks.net