surfacegmt - Online di Cloud

PROGRAM:

NAMA

permukaan - Data tabel kisi menggunakan splines kelengkungan kontinu tegangan yang dapat disesuaikan

RINGKASAN

permukaan [ tabel ] file keluaran.nc kenaikan wilayah [ aspek_rasio ] [ batas_konvergensi ] [
lmenurunkan ] [ -Luatas ] [ max_iterations ] [ ] [ radius_pencarian[m|s] ] [ [i|b]faktor_ketegangan
] [ [tingkat] ] [ faktor_relaksasi berlebihan ] [ -a] [ -dua] [ -du] [
-f] [ -h] [ -i] [ -:[i|o] ]

Catatan: Tidak ada spasi yang diperbolehkan antara flag opsi dan argumen terkait.

DESKRIPSI

permukaan membaca secara acak (x,y,z) tiga kali lipat dari input standar [atau tabel] dan menghasilkan
file grid biner dari nilai grid z(x,y) dengan memecahkan:
(1 - T) * L (L (z)) + T * L (z) = 0

di mana T adalah faktor tegangan antara 0 dan 1, dan L menunjukkan operator Laplacian. T = 0
memberikan solusi "kelengkungan minimum" yang setara dengan SuperMISP dan ISM
paket. Kelengkungan minimum dapat menyebabkan osilasi yang tidak diinginkan dan maksima lokal palsu atau
minima (Lihat Smith dan Wessel, 1990), dan Anda mungkin ingin menggunakan T > 0 untuk menekannya
efek. Pengalaman menyarankan T ~ 0.25 biasanya terlihat bagus untuk data lapangan potensial dan T
harus lebih besar (T ~ 0.35) untuk data topografi curam. T = 1 memberikan permukaan harmonik (tidak
maxima atau minima dimungkinkan kecuali pada titik data kontrol). Disarankan agar
pengguna melakukan pra-proses data dengan blok berarti, median blok, atau mode blok untuk menghindari spasial
aliasing dan menghilangkan data yang berlebihan. Anda dapat memberlakukan batas bawah dan/atau batas atas pada
larutan. Ini dapat dimasukkan dalam bentuk nilai tetap, kisi dengan nilai, atau hanya
menjadi nilai input data minimum/maksimum. Kondisi batas alami diterapkan pada
tepi, kecuali untuk data geografis dengan rentang 360 derajat di mana kami menerapkan batas periodik
kondisi dalam arah bujur.

DIBUTUHKAN ARGUMEN

-Gfile keluaran.nc
Nama file keluaran. Outputnya adalah biner 2-D .nc mengajukan. Perhatikan bahwa kisi terkecil
dimensi harus minimal 4.

-Ixinci[satuan][=|+][/yin[satuan][=|+]]
x_inc [dan opsional y_inc] adalah jarak kisi. Secara opsional, tambahkan sufiks
pengubah. Geografis (derajat) koordinat: Menambahkan m untuk menunjukkan menit busur atau s
untuk menunjukkan detik busur. Jika salah satu unit e, f, k, M, n or u ditambahkan
sebagai gantinya, kenaikan diasumsikan diberikan dalam meter, kaki, km, Mil, bahari
mil atau kaki survei AS, masing-masing, dan akan dikonversi ke yang setara
derajat bujur di garis lintang tengah wilayah (konversi tergantung pada
PROJ_ELLIPSOID). Jika /y_inc diberikan tetapi diatur ke 0 itu akan diatur ulang sama dengan x_inc;
jika tidak maka akan dikonversi ke derajat lintang. Semua koordinat: Jika = is
ditambahkan maka maks yang sesuai x (timur) atau y (utara) mungkin sedikit disesuaikan
agar sesuai dengan kenaikan yang diberikan [secara default kenaikan dapat disesuaikan
sedikit agar sesuai dengan domain yang diberikan]. Akhirnya, alih-alih memberikan kenaikan, Anda mungkin
tentukan jumlah of node diinginkan dengan menambahkan + ke bilangan bulat yang disediakan
argumen; kenaikan tersebut kemudian dihitung ulang dari jumlah node dan
domain. Nilai kenaikan yang dihasilkan tergantung pada apakah Anda telah memilih a
kisi-kisi-terdaftar atau grid-terdaftar-piksel; lihat Format file aplikasi untuk detailnya.
Catatan: jika -Rfile grd digunakan maka jarak grid telah diinisialisasi; menggunakan
-I untuk mengesampingkan nilai-nilai.

-R[satuan]xmin/xmax/ymin/ymax[R] (lebih ...)
Tentukan wilayah yang diminati.

OPSIONAL ARGUMEN

tabel Satu atau lebih ASCII (atau biner, lihat -dua[ncol][mengetik]) file tabel data yang menyimpan
jumlah kolom data. Jika tidak ada tabel yang diberikan maka kita membaca dari input standar.

-Aaspek_rasio
Rasio aspek. Jika diinginkan, anisotropi grid dapat ditambahkan ke persamaan. Memasuki
aspek_rasio, dimana dy = dx / aspek_rasio menghubungkan dimensi grid. [Awal =
1 mengasumsikan kisi isotropik.]

-Cbatas_konvergensi[%]
Batas konvergensi. Iterasi diasumsikan telah konvergen ketika maksimum absolut
perubahan nilai kisi apa pun kurang dari batas_konvergensi. (Satuan sama dengan data z
unit). Sebagai alternatif, berikan batas dalam persentase deviasi rms dengan menambahkan %.
[Default diskalakan ke 1e-4 dari deviasi akar-rata-rata-kuadrat data dari a
bidang yang paling cocok (kuadrat terkecil).]. Ini adalah batas konvergensi akhir di
jarak grid yang diinginkan; untuk kisi menengah (lebih kasar) konvergensi efektif
batas diskalakan oleh pengali spasi grid.

-IImenurunkan dan -Luatas
Memaksakan batasan pada solusi keluaran. lmenurunkan menetapkan batas bawah. menurunkan bisa menjadi
nama file grid dengan nilai batas bawah, nilai tetap, d untuk mengatur ke minimum
nilai masukan, atau u untuk tidak dibatasi [Default]. uatas menetapkan batas atas dan dapat
menjadi nama file grid dengan nilai batas atas, nilai tetap, d untuk mengatur ke
nilai input maksimum, atau u untuk tidak dibatasi [Default]. File kisi digunakan untuk mengatur
batas mungkin berisi NaN. Di hadapan NaNs, batas simpul yang ditutupi dengan
NaN tidak dibatasi.

-Nmax_iterations
Jumlah iterasi. Iterasi akan berhenti ketika batas_konvergensi tercapai atau
ketika jumlah iterasi mencapai max_iterations. Ini adalah iterasi terakhir
batas pada jarak grid yang diinginkan; untuk grid menengah (kasar) yang efektif
batas iterasi diskalakan oleh pengali spasi grid. [Defaultnya adalah 500.]

-Q Sarankan dimensi kisi yang memiliki faktor persekutuan terbesar yang sangat komposit. Ini
memungkinkan permukaan untuk menggunakan beberapa langkah perantara dalam solusi, menghasilkan lebih cepat
waktu berjalan dan hasil yang lebih baik. Ukuran yang disarankan oleh -Q dapat dicapai dengan mengubah
-R dan / atau -I. Anda dapat memulihkan -R dan -I Anda inginkan nanti dengan menggunakan contoh grd or
jalan pintas pada keluaran permukaan.

-Sradius_pencarian[m|dtk]
radius pencarian. Memasuki radius_pencarian dalam satuan yang sama dengan data x,y; menambahkan m untuk menunjukkan
menit busur atau s untuk detik busur. Ini digunakan untuk menginisialisasi grid sebelum
iterasi pertama; itu tidak sepadan dengan waktu kecuali kisi kisi prima dan
tidak dapat memiliki tahapan regional. [Default = 0.0 dan tidak ada pencarian yang dilakukan.]

-T[i|b]faktor_ketegangan
Faktor ketegangan[s]. Ini harus antara 0 dan 1. Ketegangan dapat digunakan dalam
solusi interior (di atas persamaan, di mana ia menekan osilasi palsu) dan
dalam kondisi batas (di mana ia cenderung meratakan solusi mendekati
tepi). Menggunakan nol untuk kedua nilai menghasilkan permukaan kelengkungan minimum dengan free
tepi, yaitu, spline bicubic alami. Menggunakan -Tifaktor_ketegangan untuk mengatur interior
ketegangan, dan -Tbfaktor_ketegangan untuk mengatur tegangan batas. Jika Anda tidak menambahkan i or
b, keduanya akan disetel ke nilai yang sama. [Default = 0 untuk keduanya memberikan minimum
solusi kelengkungan.]

-V[tingkat] (lebih ...)
Pilih tingkat verbositas [c]. -V3 akan melaporkan konvergensi setelah setiap iterasi;
-V akan melaporkan hanya setelah setiap grid regional terkonvergensi.

-Zfaktor_relaksasi berlebihan
Faktor over-relaksasi. Parameter ini digunakan untuk mempercepat konvergensi; ini
angka antara 1 dan 2. Nilai 1 mengulangi persamaan dengan tepat, dan akan
selalu menjamin konvergensi yang stabil. Nilai yang lebih besar melebih-lebihkan inkremental
berubah selama konvergensi, dan akan mencapai solusi lebih cepat tetapi dapat menjadi
tidak stabil. Jika Anda menggunakan nilai besar untuk faktor ini, ada baiknya untuk memantau
setiap iterasi dengan -Vl pilihan. [Default = 1.4 konvergen dengan cepat dan hampir
selalu stabil.]

-acol=nama[...] (lebih ...)
Tetapkan asosiasi kolom aspatial col=nama.

-dua[ncol][T] (lebih ...)
Pilih masukan biner asli. [Default adalah 3 kolom input].

-dutidak ada data (lebih ...)
Ganti kolom input yang sama tidak ada data dengan NaN.

-f[i|o]com.colinfo (lebih ...)
Tentukan tipe data kolom input dan/atau output.

-h[i|o]n][+c][+d][+rkomentar] [+ rjudul] (lebih ...)
Lewati atau buat rekaman header. Tidak digunakan dengan data biner.

-ikerah[aku] [sskala][Haimengimbangi] [,...] (lebih ...)
Pilih kolom input (0 adalah kolom pertama).

-:[i|o] (lebih ...)
Tukar kolom 1 dan 2 pada 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 NILAI KETEPATAN

Terlepas dari ketepatan input data, program GMT yang membuat file grid akan
secara internal menahan grid dalam array floating point 4-byte. Ini dilakukan untuk menghemat memori
dan selanjutnya sebagian besar jika tidak semua data nyata dapat disimpan menggunakan floating point 4-byte
nilai-nilai. Data dengan presisi yang lebih tinggi (yaitu, nilai presisi ganda) akan kehilangan itu
presisi setelah GMT beroperasi di grid atau menulis grid baru. Untuk membatasi kehilangan
presisi saat memproses data, Anda harus selalu mempertimbangkan untuk menormalkan data sebelum
pengolahan.

CONTOH

Untuk grid 5 oleh 5 menit berarti blok gravitasi dari data ASCII di hawaii_5x5.xyg, menggunakan
faktor_ketegangan = 0.25, a batas_konvergensi = 0.1 miligal, menulis hasilnya ke file
disebut hawaii_grd.nc, dan memantau setiap iterasi, coba:

gmt permukaan hawaii_5x5.xyg -R198/208/18/25 -I5m -Ghawaii_grd.nc -T0.25 -C0.1 -Vl

Gunakan surfacegmt online menggunakan layanan onworks.net