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

ПРОГРАММА:

ИМЯ

dimfilter - Направленная фильтрация файлов с двумерной сеткой в ​​пространственной (или временной) области

СИНТАКСИС

димфильтр input_file.nc Distance_flag [Режим] выходной_файл.nc
[ смещ_по_столбцам ] [ увеличить ] [ область ] [] [[уровень]] [ -f]

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

ОПИСАНИЕ

димфильтр отфильтрует .NC файл в пространственной (или временной) области путем разделения заданного
фильтровать круг в n_секторов, применяя одну из выбранных первичной свертки или
несверточные фильтры для каждого сектора и выбор окончательного результата в соответствии с
выбранный вторичный фильтр. Он вычисляет расстояния с использованием декартовой или сферической геометрии.
Выход .NC файл может быть необязательно сгенерирован как подобласть ввода и / или с
new -Increment. Таким образом, во входных данных может быть «лишнее пространство», так что там
не будет краевых эффектов для выходной сетки. Если фильтр низкочастотный, то на выходе
выборка может выполняться реже, чем входная. -Q предназначен для режима анализа ошибок и только
требуется общее количество столбцов во входном файле, который содержит отфильтрованные
глубины. Наконец, следует знать, что димфильтр не даст плавного вывода, как другие
пространственные фильтры работают, потому что они возвращают минимальную медиану из N медианы N сектора.
вывод может быть грубым, если входные данные не содержат шумов. Таким образом, дополнительная фильтрация
(например, гауссовский через грдфильтр) данных, отфильтрованных с помощью DiM, обычно рекомендуется.

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

input_file.nc
Сетка данных для фильтрации.

-DDistance_flag
Расстояние флажок сообщает, как сетка (x, y) относится к фильтру ширина, следующим образом:

флажок = 0: сетка (x, y) те же единицы, что и ширина, Декартовы расстояния. флажок = 1: сетка
(x, y) в градусах, ширина в километрах, декартовых расстояниях. флажок = 2: сетка (x, y)
в градусах, ширина в км, dx в масштабе cos (середина y), декартовы расстояния.

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

флажок = 3: сетка (x, y) в градусах, ширина в км, dx с косинусом (y), декартово
расчет расстояния.

флажок = 4: сетка (x, y) в градусах, ширина в км, расчет сферического расстояния.

-F[Режим]
Устанавливает тип основного фильтра. Выберите фильтр свертки или фильтр без свертки.
Добавьте код фильтра, после которого укажите полный диаметр. ширина. Доступная свертка
фильтры:

(bКрытый вагон: все веса одинаковы.

(c) Косинусная дуга: веса следуют кривой косинусоидальной дуги.

(g) Гауссовский: веса задаются функцией Гаусса.

Несверточные фильтры:

(m) Median: возвращает среднее значение.

(p) Максимальная вероятность правдоподобия (оценка режима): возврат модального значения. Если больше
чем один режим найден, возвращаем их среднее значение. Добавить - или + к фильтру
width, если вы предпочитаете возвращать наименьшее или наибольшее из модальных значений.

-N
Устанавливает тип вторичного фильтра и количество секторов галстука-бабочки. n_секторов должен
быть целым и больше 0. Когда n_секторов установлен на 1, вторичный фильтр
не эффективны. Доступные вторичные фильтры:

(l) Нижний: возвращает минимум всех отфильтрованных значений.

(u) Верхний: возвращает максимум всех отфильтрованных значений.

(a) Среднее: возвращает среднее всех отфильтрованных значений.

(m) Медиана: возвращает медиану всех отфильтрованных значений.

(p) Mode: возврат режима всех отфильтрованных значений.

-Gвыходной_файл.nc
выходной_файл.nc это выход фильтра.

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

-I x_inc [и, возможно, y_inc] - выходное приращение. Добавить m для обозначения минут,
or c для обозначения секунд. Если новый x_inc, y_inc НЕ являются целыми кратными
старые (во входных данных) фильтрация будет значительно медленнее. [По умолчанию: То же
как ввод.]

-R запад, восток, и к северу определяет область выходных точек. [Дефолт:
Так же, как вход.]

-T Переключите регистрацию узла для выходной сетки, чтобы она стала противоположностью
сетка ввода [по умолчанию дает такую ​​же регистрацию, как и сетка ввода].

-Qсмещ_по_столбцам смещ_по_столбцам - общее количество столбцов во входном текстовом табличном файле. Для этого режима
он рассчитывает прочитать глубины, состоящие из нескольких столбцов. Каждый столбец представляет собой
фильтрованная сетка с шириной фильтра, которая может быть получена грд2xyz -Z. Результат
будет медианным, безумным и средним. Итак, столбец с медианами используется для генерации
региональная составляющая и столбец с MAD используются для проведения ошибки
анализа.

-V [уровень] (Больше ...)
Выберите уровень детализации [c].

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

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

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

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

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

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

GRID ФАЙЛОВ ФОРМАТЫ

По умолчанию GMT записывает сетку как числа с плавающей запятой одинарной точности в netCDF жалобы COARDS.
формат файла. Тем не менее, GMT может создавать файлы сетки во многих других часто используемых сетках.
форматы файлов, а также облегчает так называемую «упаковку» сеток с записью с плавающей запятой
данные в виде 1- или 2-байтовых целых чисел. Чтобы указать точность, масштаб и смещение, пользователь должен
добавить суффикс =id[/лестница/смещение[/бабушка]], где id двухбуквенный идентификатор сетки
тип и точность, и лестница и смещение являются необязательными масштабным коэффициентом и смещением, которые должны быть
применяется ко всем значениям сетки, и бабушка это значение, используемое для обозначения отсутствующих данных. В случае
два персонажа id не предусмотрено, как в =/лестница чем id=nf предполагается. Когда
считывая сетки, формат обычно распознается автоматически. Если нет, то тот же суффикс
могут быть добавлены к именам файлов входной сетки. Видеть грдконвертировать и Секция grid-file-format из
Технический справочник и поваренная книга GMT для получения дополнительной информации.

При чтении файла netCDF, содержащего несколько сеток, GMT по умолчанию считывает
первая двумерная сетка, которую можно найти в этом файле. Чтобы уговорить GMT прочитать другой
многомерная переменная в файле сетки, добавить ?VarName к имени файла, где
VarName это имя переменной. Обратите внимание, что вам может потребоваться избежать специального значения
of ? в вашей программе оболочки, поставив перед ней обратную косую черту или поместив
имя файла и суффикс в кавычках или двойных кавычках. В ?VarName суффикс также может быть использован
для выходных сеток указать имя переменной, отличное от значения по умолчанию: «z». Видеть
грдконвертировать и модификаторы разделов для CF и формата файла сетки технических средств GMT.
Справочник и поваренная книга для получения дополнительной информации, особенно о том, как читать сращивания 3-,
4- или 5-мерные сетки.

ГЕОГРАФИЧЕСКИЙ И ВРЕМЯ КООРДИНАТЫ

Если тип выходной сетки - netCDF, координаты будут помечены как «долгота»,
"широта" или "время" на основе атрибутов входных данных или сетки (если есть) или
-f or -R параметры. Например, оба -f0x -f1t и -R90w / 90e / 0t / 3t приведет к
сетка долготы / времени. Когда координаты x, y или z - время, они будут сохранены в сетке.
как относительное время с эпохи, как указано TIME_UNIT и TIME_EPOCH в gmt.conf файл
или в командной строке. В дополнение Ед. изм атрибут временной переменной укажет
и эта единица и эпоха.

ПРИМЕРЫ

Предположим, что north_pacific_dbdb5.nc - это файл 5-минутной батиметрии от 140E до 260E и
От 0N до 50N, и вы хотите найти медианы значений в радиусе 300 км (600 км
width) выходных точек, которые вы выбираете от 150E до 250E и от 10N до 40N, и
вам нужны выходные значения каждые 0.5 градуса. Чтобы не допустить смещения медианы
наклонной плоскости, вы хотите разделить круг фильтра на 6 секторов и выбрать
наименьшее значение среди 6 медиан. Используя расчет сферических расстояний, вам необходимо:

gmt dimfilter North_pacific_dbdb5.nc -Gfiltered_pacific.nc -Fm600 -D4 \
-Nl6 -R150 / 250/10/40 -I0.5 -V

Предположим, что cape_verde.nc - это файл 0.5-минутной батиметрии от 32 до 15 Вт и от 8 до
25N, и вы хотите удалить мелкомасштабные элементы, чтобы определить волнение в
область, простирающаяся от 27.5 до 20.5 Вт и от 12.5 до 19.5 Н, и вы хотите, чтобы выходное значение выходило каждые
2 мин. Используя вычисления декартовых расстояний, вам необходимо:

Фильтр затемнения gmt cape_verde.nc -Gt.nc -Fm220 -Nl8 -D2 -R-27.5 / -20.5 / 12.5 / 19.5 -I2m -V
gmt grdfilter t.nc -Gcape_swell.nc -Fg50 -D2 -V

Предположим, вы нашли диапазон ширины фильтра для данной области и отфильтровали
заданные батиметрические данные с использованием диапазона ширины фильтра (например, f100.nc f110.nc f120.nc
f130.nc), и вы хотите определить региональный тренд, используя диапазон ширины фильтра, и
вы хотите получить оценки среднего абсолютного отклонения (MAD) в каждой точке данных. Затем вы
нужно будет сделать:

gmt grd2xyz f100.nc -Z> f100.d
gmt grd2xyz f110.nc -Z> f110.d
gmt grd2xyz f120.nc -Z> f120.d
gmt grd2xyz f130.nc -Z> f130.d
паста f100.d f110.d f120.d f130.d> глубины.d
gmt dimfilter depths.d -Q4> output.z

ОГРАНИЧЕНИЯ

При работе с географическими сетками (широта, долгота) все три фильтра свертки (boxcar,
косинусная дуга и гауссовский) правильно нормализует веса фильтра для изменения
размер сетки с широтой и правильно определить, какие узлы необходимы для
свертка, когда фильтр "круг" пересекает периодическую (0–360) границу или содержит
географический полюс. Однако пространственные фильтры, такие как медианные и модовые фильтры, не используют
веса и, следовательно, должны использоваться только в декартовых сетках (или на очень низких широтах).
Если вы хотите применить такие пространственные фильтры, вы должны проецировать данные на равновеликий
projection и запустите dimfilter на получившейся декартовой сетке.

SCRIPT ШАБЛОН

Dim.template.sh - это скелетный сценарий оболочки, который можно использовать для настройки полного DiM
анализ, в том числе анализ MAD.

СПРАВКА

Ким, С.-С., и Вессель, П. (2008), Направленная медианная фильтрация для регионально-остаточной
Разделение батиметрии, Геохим. Геофиз. Геосист., 9, Q03005,
DOI: 10.1029 / 2007GC001850.

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