Ini adalah perintah dimfiltergmt yang dapat dijalankan di penyedia hosting gratis OnWorks menggunakan salah satu dari beberapa workstation online gratis kami seperti Ubuntu Online, Fedora Online, emulator online Windows atau emulator online MAC OS
PROGRAM:
NAMA
dimfilter - Pemfilteran terarah dari file grid 2-D dalam domain ruang (atau waktu)
RINGKASAN
filter redup masukan_file.nc jarak_bendera [mode] keluaran_file.nc
[ kerah ] [ kenaikan ] [ wilayah ] [ ] [ [tingkat] ] [ -f]
Catatan: Tidak ada spasi yang diperbolehkan antara flag opsi dan argumen terkait.
DESKRIPSI
filter redup akan menyaring .nc file dalam ruang (atau waktu) domain dengan membagi yang diberikan
saring lingkaran menjadi n_sektor, menerapkan salah satu konvolusi primer yang dipilih atau
filter non-konvolusi untuk setiap sektor, dan memilih hasil akhir sesuai dengan
filter sekunder yang dipilih. Ini menghitung jarak menggunakan geometri Cartesian atau Spherical.
Hasil .nc file opsional dapat dihasilkan sebagai subregion dari input dan/atau dengan a
yang baru -Imeningkat. Dengan cara ini, seseorang mungkin memiliki "ruang ekstra" dalam data input sehingga ada
tidak akan ada efek tepi untuk kisi keluaran. Jika filter low-pass, maka outputnya
sampel mungkin lebih jarang daripada input. -Q adalah untuk mode analisis kesalahan dan hanya
membutuhkan jumlah total kolom dalam file input, yang berisi filter
kedalaman. Akhirnya, orang harus tahu itu filter redup tidak akan menghasilkan output yang mulus seperti yang lain
filter spasial melakukannya karena mengembalikan median minimum dari N median dari N sektor. NS
output bisa kasar kecuali data input bebas noise. Dengan demikian, penyaringan tambahan
(misalnya, Gaussian melalui filter grd) dari data yang difilter DiM umumnya direkomendasikan.
DIBUTUHKAN ARGUMEN
masukan_file.nc
Kisi data yang akan difilter.
-Djarak_bendera
Jarak bendera memberitahu bagaimana grid (x,y) berhubungan dengan filter lebar, sebagai berikut:
bendera = 0: kisi (x,y) satuan yang sama dengan lebar, jarak Cartesian. bendera = 1: kisi
(x,y) dalam derajat, lebar dalam kilometer, jarak Cartesian. bendera = 2: kisi (x,y)
dalam derajat, lebar dalam km, dx diskalakan oleh cos(y tengah), jarak Cartesian.
Opsi di atas adalah yang tercepat karena hanya memungkinkan matriks bobot dihitung
satu kali. Tiga opsi berikutnya lebih lambat karena mereka menghitung ulang bobot untuk masing-masing
Garis Lintang.
bendera = 3: kisi (x,y) dalam derajat, lebar dalam km, dx diskalakan oleh cosinus(y), Cartesian
perhitungan jarak.
bendera = 4: kisi (x,y) dalam derajat, lebar dalam km, perhitungan jarak Bulat.
-F[mode]
Menyetel jenis filter utama. Pilih di antara filter konvolusi dan non-konvolusi.
Tambahkan kode filter diikuti dengan diameter penuh lebar. Konvolusi yang tersedia
filter adalah:
(b) Boxcar: Semua bobot sama.
(c) Lengkungan Cosinus: Bobot mengikuti kurva lengkung kosinus.
(g) Gaussian: Bobot diberikan oleh fungsi Gaussian.
Filter non-konvolusi adalah:
(m) Median: Mengembalikan nilai median.
(p) Probabilitas kemungkinan maksimum (penaksir mode): Mengembalikan nilai modal. Jika lebih
dari satu mode ditemukan, kami mengembalikan nilai rata-ratanya. Tambahkan - atau + ke filter
lebar jika Anda lebih suka mengembalikan nilai modal terkecil atau terbesar.
-N
Menyetel jenis filter sekunder dan jumlah sektor dasi kupu-kupu. n_sektor harus
menjadi bilangan bulat dan lebih besar dari 0. Ketika n_sektor diatur ke 1, filter sekunder adalah
tidak efektif. Filter sekunder yang tersedia adalah:
(l) Bawah: Mengembalikan nilai minimum dari semua nilai yang difilter.
(u) Atas: Mengembalikan nilai maksimum dari semua nilai yang difilter.
(a) Rata-rata: Mengembalikan rata-rata dari semua nilai yang difilter.
(m) Median: Mengembalikan median dari semua nilai yang difilter.
(p) Mode: Mengembalikan mode dari semua nilai yang difilter.
-Gkeluaran_file.nc
keluaran_file.nc adalah keluaran filter.
OPSIONAL ARGUMEN
-I x_inc [dan opsional y_inc] adalah peningkatan keluaran. Menambahkan m untuk menunjukkan menit,
or c untuk menunjukkan detik. Jika baru x_inc, y_inc BUKAN kelipatan bilangan bulat dari
yang lama (dalam data input), penyaringan akan jauh lebih lambat. [Awal: Sama
sebagai masukan.]
-R barat, timur, selatan, dan utara mendefinisikan Wilayah titik output. [Bawaan:
Sama seperti yang dimasukkan.]
-T Alihkan pendaftaran simpul untuk kisi keluaran sehingga menjadi kebalikan dari
grid input [Default memberikan registrasi yang sama dengan grid input].
-Qkerah kerah adalah jumlah total kolom dalam file tabel teks input. Untuk modus ini,
itu mengharapkan untuk membaca kedalaman terdiri dari beberapa kolom. Setiap kolom mewakili a
kisi yang difilter dengan lebar filter, yang dapat diperoleh dengan grd2xyz -Z. Hasilnya
akan menjadi median, MAD, dan mean. Jadi, kolom dengan median digunakan untuk menghasilkan
komponen regional dan kolom dengan MAD digunakan untuk melakukan kesalahan
analisis.
-V[tingkat] (lebih ...)
Pilih tingkat verbositas [c].
-f[i|o]com.colinfo (lebih ...)
Tentukan tipe data kolom input dan/atau output.
-^ or hanya -
Cetak pesan singkat tentang sintaks perintah, lalu keluar (CATATAN: di Windows
gunakan saja -).
-+ or hanya +
Cetak pesan penggunaan (bantuan) ekstensif, termasuk penjelasan tentang apa pun
opsi khusus modul (tetapi bukan opsi umum GMT), lalu keluar.
-? or tidak argumen
Cetak pesan penggunaan (bantuan) lengkap, termasuk penjelasan opsi, lalu
keluar.
--Versi: kapan
Cetak versi GMT dan keluar.
--tunjukkan-datadir
Cetak jalur lengkap ke direktori berbagi GMT dan keluar.
GRID FILE FORMAT
Secara default GMT menulis grid sebagai presisi tunggal mengapung di netCDF keluhan COARDS
format berkas. Namun, GMT mampu menghasilkan file grid di banyak grid lain yang umum digunakan
format file dan juga memfasilitasi apa yang disebut "pengemasan" grid, menulis floating point
data sebagai bilangan bulat 1- atau 2-byte. Untuk menentukan presisi, skala, dan offset, pengguna harus
tambahkan akhiran =id[/skala/mengimbangi[/nan]], di mana id adalah pengidentifikasi dua huruf dari kisi
jenis dan presisi, dan skala dan mengimbangi adalah faktor skala opsional dan offset menjadi
diterapkan ke semua nilai kisi, dan nan adalah nilai yang digunakan untuk menunjukkan data yang hilang. Dalam hal
dua karakter id tidak disediakan, seperti pada =/skala dari a id=nf diasumsikan. Kapan
membaca grid, format umumnya secara otomatis dikenali. Jika tidak, sufiks yang sama
dapat ditambahkan ke input nama file grid. Lihat konversi dan Format file-grid Bagian dari
Referensi Teknis GMT dan Buku Masak untuk informasi lebih lanjut.
Saat membaca file netCDF yang berisi beberapa kisi, GMT akan membaca, secara default,:
grid 2 dimensi pertama yang dapat ditemukan di file itu. Untuk membujuk GMT agar membaca yang lain
variabel multi-dimensi dalam file grid, tambahkan ?nama belakang ke nama file, di mana
nama belakang adalah nama variabelnya. Perhatikan bahwa Anda mungkin perlu menghindari makna khusus
of ? di program shell Anda dengan meletakkan garis miring terbalik di depannya, atau dengan menempatkan
nama file dan akhiran antara tanda kutip atau tanda kutip ganda. NS ?nama belakang akhiran juga dapat digunakan
untuk kisi keluaran untuk menentukan nama variabel yang berbeda dari default: "z". Lihat
konversi dan Pengubah Bagian-untuk-CF dan format file-grid dari Teknis GMT
Referensi dan Buku Masak untuk informasi lebih lanjut, terutama tentang cara membaca splices dari 3-,
4-, atau grid 5 dimensi.
GEOGRAFIS DAN WAKTU KOORDINASI
Ketika tipe grid keluaran adalah netCDF, koordinat akan diberi label "bujur",
"lintang", atau "waktu" berdasarkan atribut input data atau grid (jika ada) atau pada
-f or -R pilihan. Misalnya, keduanya -f0x -f1t dan -R90w/90e/0t/3t akan menghasilkan a
garis bujur/waktu. Ketika koordinat x, y, atau z adalah waktu, itu akan disimpan di grid
sebagai waktu relatif sejak zaman sebagaimana ditentukan oleh TIME_UNIT dan TIME_EPOCH di gmt.conf fillet
atau di baris perintah. Selain itu, satuan atribut variabel waktu akan menunjukkan
baik unit ini dan zaman.
CONTOH
Misalkan north_pacific_dbdb5.nc adalah file batimetri 5 menit dari 140E ke 260E dan
0N hingga 50N, dan Anda ingin mencari median nilai dalam radius 300km (600km penuh
lebar) dari titik keluaran, yang Anda pilih dari 150E hingga 250E dan 10N hingga 40N, dan
Anda ingin nilai output setiap 0.5 derajat. Untuk mencegah median menjadi bias oleh
bidang miring, Anda ingin membagi lingkaran filter menjadi 6 sektor dan untuk memilih
nilai terendah di antara 6 median. Menggunakan perhitungan jarak bola, Anda membutuhkan:
gmt dimfilter north_pacific_dbdb5.nc -Gfiltered_pacific.nc -Fm600 -D4 \
-Nl6 -R150/250/10/40 -I0.5 -V
Misalkan cape_verde.nc adalah file batimetri 0.5 menit dari 32W ke 15W dan 8N ke
25N, dan Anda ingin menghapus fitur skala panjang kecil untuk mendefinisikan swell dalam
area yang membentang dari 27.5W hingga 20.5W dan 12.5N hingga 19.5N, dan Anda menginginkan nilai keluaran setiap
2 menit. Dengan menggunakan perhitungan jarak kartesius, Anda memerlukan:
gmt dimfilter 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
Misalkan Anda menemukan rentang lebar filter untuk area tertentu, dan Anda memfilter
diberikan data batimetri menggunakan kisaran lebar filter (misalnya, f100.nc f110.nc f120.nc
f130.nc), dan Anda ingin menentukan tren regional menggunakan rentang lebar filter, dan
Anda ingin mendapatkan perkiraan deviasi absolut median (MAD) pada setiap titik data. Maka kamu
perlu melakukan:
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
rekatkan f100.d f110.d f120.d f130.d > depths.d
gmt dimfilter depths.d -Q4 > output.z
PEMBATASAN
Saat bekerja dengan grid geografis (lat, lon), ketiga filter konvolusi (boxcar,
lengkung kosinus, dan gaussian) akan menormalkan dengan benar bobot filter untuk variasi dalam
ukuran kotak kisi dengan garis lintang, dan tentukan dengan benar simpul mana yang diperlukan untuk
konvolusi ketika "lingkaran" filter melintasi batas periodik (0-360) atau mengandung a
kutub geografis. Namun, filter spasial, seperti filter median dan mode, tidak menggunakan
bobot dan dengan demikian hanya boleh digunakan pada kisi Cartesian (atau pada garis lintang yang sangat rendah) saja.
Jika Anda ingin menerapkan filter spasial seperti itu, Anda harus memproyeksikan data Anda ke area yang sama
proyeksi dan jalankan dimfilter pada kisi Cartesian yang dihasilkan.
SCRIPT TEMPLAT
dim.template.sh adalah skrip kerangka kerangka yang dapat digunakan untuk menyiapkan DiM . lengkap
analisis, termasuk analisis MAD.
REFERENSI
Kim, S.-S., dan Wessel, P. (2008), Directional Median Filtering untuk Regional-Residual
Pemisahan Batimetri, Geokimia. Geofis. Geosistem., 9, Q03005,
doi:10.1029/2007GC001850.
Gunakan dimfiltergmt online menggunakan layanan onworks.net