grdfiltergmt - Інтернет у хмарі

ПРОГРАМА:

ІМ'Я

grdfilter - фільтрує сітку в просторовій (або часовій) області

СИНТАКСИС

grdfilter Інгрід прапор_дистанції xширина[/ширина2][модифікатори] позасітка [ приріст ] [ i|p|r ]
[ регіон ] [ ] [ [рівень] ] [ -f]

Примітка: Між прапорцем параметра та пов’язаними аргументами не допускається пробіл.

ОПИС

grdfilter фільтруватиме файл сітки у часовій області за допомогою однієї з вибраних згорток
або ізотропні або прямокутні фільтри без згортки та обчислюйте відстані за допомогою декартової системи
або сферичні геометрії. Вихідний файл сітки можна за бажанням створити як підобласть
введення (через -R) та/або з новим приростом (через -I) або реєстрація (через -T). У цьому
Таким чином, у вхідних даних може бути "додатковий простір", тому краї не використовуватимуться
вивід може бути в межах однієї половини ширини вхідних країв. Якщо фільтр низьких частот,
тоді вихідні дані можуть дискретизуватися рідше, ніж вхідні.

ВИМАГАЄТЬСЯ АРГУМЕНТИ

Інгрід Файл сітки точок, які потрібно відфільтрувати. (Див. ФОРМАТИ ФАЙЛІВ СІТКИ нижче).

-Dпрапор_дистанції
відстань прапор розповідає, як сітка (x,y) пов’язана з фільтром ширина наступним чином:

прапор = p: сітка (px,py) з ширина непарна кількість пікселів; Декартові відстані.

прапор = 0: сітка (x,y) ті ж одиниці, що і ширина, декартові відстані.

прапор = 1: сітка (x,y) у градусах, ширина в кілометрах, декартові відстані.

прапор = 2: сітка (x,y) у градусах, ширина у км, dx у масштабі cos(середнє y), декартова система
відстані.

Наведені вище варіанти є найшвидшими, оскільки дозволяють обчислювати лише матрицю ваги
один раз. Наступні три варіанти повільніші, оскільки вони перераховують ваги для кожного
широта.

прапор = 3: сітка (x,y) у градусах, ширина в км, dx масштабований косинусом(y), декартовим
розрахунок відстані.

прапор = 4: сітка (x,y) у градусах, ширина в км, розрахунок сферичної відстані.

прапор = 5: сітка (x,y) у Mercator -Джм1 одиниця зображення, ширина в км, Сферична відстань
розрахунок

-Fxширина[/ширина2][модифікатори]
Встановлює тип фільтра. Виберіть між згортковими та неконволюційними фільтрами. Використовуйте будь-який
код фільтра x (перелічено нижче), а потім повний діаметр ширина. Це дає
ізотропний фільтр; додати /ширина2 для прямокутного фільтра (потрібно -Dp or -0 DD).
Деякі фільтри допускають необов’язкові аргументи та модифікатори.

Згорткові фільтри (і їхні коди) такі:

(b) Вагон: Усі ваги рівні.

(c) Арка косинуса: ваги слідують за кривою косинусної арки.

(g) Гаусс: ваги задано функцією Гаусса, де ширина становить 6 разів
звичайна сигма Гауса.

(f) Спеціальний: ваги надаються попередньо обчисленими значеннями в таблиці ваг фільтра
файл вага, які повинні мати непарні розміри; також вимагає -0 DD і вихідний інтервал
має відповідати інтервалу введення або бути цілим кратним.

(o) Оператор: Вагові коефіцієнти задаються попередньо обчисленими значеннями в таблиці вагових коефіцієнтів фільтра
файл вага, які повинні мати непарні розміри; також вимагає -0 DD і вихідний інтервал
має відповідати інтервалу введення або бути цілим кратним. Припускається, що сума ваги дорівнює
нуль, тому накопичення сум ваг і нормалізація не відбуватиметься.

Неконволюційні фільтри (і їхні коди) це:

(m) Медіана: повертає середнє значення. Щоб вибрати інший квантиль, додайте +qквантил in
діапазон 0-1 [за замовчуванням 0.5, тобто медіана].

(p) Імовірність максимальної правдоподібності (оцінка режиму): повертає модальне значення. Якщо більше
ніж один режим знайдено, ми повертаємо їх середнє значення. Додайте +l or +u якщо вам краще
хочете повернути найнижче або найвище з модальних значень.

(h) Режим гістограми (інший режим оцінки): повертає модальне значення як центр
домінуючого піку на гістограмі. Додати /binwidth щоб указати групування
інтервал. Використовуйте модифікатор +c щоб центрувати контейнери на кратних binwidth [За замовчуванням має
ребра контейнера, кратні binwidth]. Якщо знайдено більше ніж один режим ми
повернути їх середнє значення. Додайте +l or +u якщо ви бажаєте повернути
найнижче або найвище з модальних значень.

(l) Нижня: повертає мінімальне значення з усіх значень.

(L) Нижній: повертає лише мінімум усіх позитивних значень.

(u) Верхній: повертає максимум усіх значень.

(U) Верхній: повертає максимум або лише всі від’ємні значення.

У разі L|U можливо, що жодні дані не пройдуть початковий тест на ознаки; в
у цьому випадку фільтр поверне NaN.

-Gпозасітка
позасітка є вихідним файлом сітки фільтра. (Див. ФОРМАТИ ФАЙЛІВ СІТКИ нижче).

ДОДАТКОВО АРГУМЕНТИ

-Ixinc[блок][=|+][/yinc[блок][=|+]]
x_inc [і за бажанням y_inc] – вихідний приріст. Додати m для позначення дуги
хвилин, або s для позначення кутових секунд. Якщо новий x_inc, y_inc НЕ є цілими
кратні старим (у вхідних даних), фільтрація буде значною
повільніше. [За замовчуванням: те саме, що введено.]

-Ni|p|r
Визначте, як NaN-значення у вхідній сітці впливають на відфільтрований вихід: Додайте i до
ігнорувати всі NaN при обчисленні відфільтрованого значення [за замовчуванням], r це те саме, що i
за винятком випадків, коли вхідним вузлом було NaN, тоді вихідним вузлом буде встановлено NaN (тільки
застосовується, якщо обидві мережі зареєстровані спільно), і p що примусово відфільтровано значення
бути NaN, якщо всередині кола фільтра знайдено будь-які вузли сітки зі значеннями NaN.

-R захід, схід, південь та на північ визначає регіон вихідних точок. [За замовчуванням:
Те саме, що введення.]

-T Переключіть реєстрацію вузла для вихідної сітки, щоб стати протилежністю
сітка введення [за замовчуванням дає таку ж реєстрацію, як і сітка введення].

-V[рівень] (більше ...)
Виберіть рівень детальності [c].

-f [я | о]colinfo (більше ...)
Вкажіть типи даних вхідних та/або вихідних стовпців.

-^ or просто -
Надрукуйте коротке повідомлення про синтаксис команди, а потім завершує роботу (ПРИМІТКА: у Windows
використовувати тільки -).

-+ or просто +
Надрукуйте розширене повідомлення (довідка) про використання, включаючи пояснення будь-якого
параметр, специфічний для модуля (але не загальні параметри GMT), потім виходить.

-? or немає аргументація
Потім надрукуйте повне повідомлення про використання (довідку), включаючи пояснення параметрів
виходи.

-- версія
Роздрукуйте версію GMT і вийдіть.

--show-datadir
Роздрукуйте повний шлях до загального каталогу GMT та вийдіть.

GRID Фото ФОРМАТИ

За замовчуванням GMT виписує сітку, як з одинарної точністю поплавці в COARDS-скарзі NetCDF
формат файлу. Проте, GMT здатний створювати файли сітка під багато іншої часто респонденти користуються послугами сітці
формати файлів, а також сприяє так званої «упаковки» сіток, виписуючи з плаваючою точкою
дані як 1- або 2-байтові цілі числа. Для того, щоб визначити точність, масштаб і зміщення, користувач повинен
додати суфікс =id[/масштаб/зсув[/бабуся]], Де id є двобуквенним ідентифікатором сітки
тип і точність, і масштаб та зсув є необов'язковим масштабним коефіцієнтом і зміщення буде
застосовується до всіх значень сітки і бабуся це значення використовується для вказівки відсутніх даних. в разі
двох персонажів id не передбачено, як і в =/масштаб ніж a id=nf передбачається. Коли
читаючи сітки, формат зазвичай автоматично розпізнається. Якщо ні, той самий суфікс
можуть бути додані вхідні імена файлів сітки. Побачити grdconvert і розділ сітка-формат файлу
Технічний довідник GMT та Кулінарна книга для отримання додаткової інформації.

Під час читання файлу netCDF, який містить кілька сіток, GMT за замовчуванням прочитає файл
перших 2-мірної сітки, які можна знайти в цьому файлі. Щоб переконати GMT прочитати інше
багатовимірна змінна в файлі сітки, Append ?ім'я вар до імені файлу, де
ім'я вар це ім'я змінної. Зауважте, що вам може знадобитися уникнути особливого значення
of ? у вашій оболонки програми, поставивши зворотну косу риску перед ним, або шляхом розміщення
ім'я файлу та суфікс між лапками або подвійними лапками. The ?ім'я вар Суфікс може також використовуватися
для виведення сітки, щоб задати змінну інше ім'я від значення за замовчуванням: «г». Побачити
grdconvert і Розділи модифікатори-для-CF і сітка-файли формати GMT Технічний
Довідник та Кулінарна книга для отримання додаткової інформації, зокрема про те, як читати зрощення 3-,
4- або 5-вимірні сітки.

ГЕОГРАФІЧНИЙ І ЧАС КООРДИНАТИ

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

ПРИКЛАДИ

Припустимо, що north_pacific_etopo5.nc є файлом 5-хвилинної батиметрії від 140E до 260E
і від 0 до 50 пн., і ви хочете знайти медіани значень у радіусі 300 км (600 км
повна ширина) вихідних точок, які ви вибираєте як від 150E до 250E та від 10N до 40N,
і вам потрібні вихідні значення кожні 0.5 градуса. Використовуючи обчислення сферичної відстані,
тобі потрібно:

gmt grdfilter north_pacific_etopo5.nc -Gfiltered_pacific.nc -Fm600
-D4 -R150/250/10/40 -I0.5 -V

Якщо натомість ми хотіли отримати результат високої частоти, тоді можна виконати відповідну низьку частоту
фільтр із використанням грубого інтервалу сітки, оскільки grdfilter передискретизує результат до того самого
роздільна здатність як вхідну сітку, щоб ми могли обчислити залишки, наприклад,

gmt grdfilter north_pacific_etopo5.nc -Gresidual_pacific.nc -Fm-600
-D4 -R150/250/10/40 -I0.5 -V

Тут сітка residual_pacific.nc матиме таку ж 5-хвилинну роздільну здатність, що й оригінал.

Щоб відфільтрувати набір даних у ripples.nc за допомогою спеціального анізотропного фільтра Гауса exp
(-0.5*r^2), чия відстань r від центру визначається як (2x^2 + y^2 -2xy)/6, з великим
вісь під кутом 63 градуси до горизонталі, проб

gmt grdmath -R-10/10/-10/10 -I1 X 2 POW 2 MUL Y 2 POW ADD XY MUL 2 MUL
SUB 6 DIV NEG 2 DIV EXP DUP SUM DIV = gfilter.nc
gmt grdfilter ripples.nc -Ffgfilter.nc -D0 -Gsmooth.nc -V

ОБМЕЖЕННЯ

Для використання -D5 опція вхідна сітка Меркатора повинна бути створена img2mercgrd за допомогою -C
таким чином, початок значень y – екватор (тобто x = y = 0 відповідає lon =
lat = 0).

Використовуйте grdfiltergmt онлайн за допомогою сервісів onworks.net