gmtspatialgmt - Online di Cloud

PROGRAM:

NAMA

gmtspatial - Lakukan operasi geospasial pada garis dan poligon

RINGKASAN

gmtspasial [ tabel ] [ [amin_dist][satuan]] [ ] [
[+ffillet][+aamax][+dd maks][+c|Ccmax][+sfakta] ] [ +|- ] [ [l] ] [ -I[e|i] ] [
file p[+a][+pawal][+r][+z] ] [ [[-|+]*satuan*][+h][+l][+p] ] [ wilayah ] [ i|u|s|j ] [
[klippoligon] ] [ [tingkat] ] [ -b] [ -d] [ -f] [ -g] [
-h] [ -i] [ -o] [ -:[i|o] ]

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

DESKRIPSI

gmtspasial membaca satu atau lebih file data (yang mungkin berupa file multisegmen) yang berisi:
poligon tertutup dan mengoperasikan poligon ini dengan cara yang ditentukan. Operasi termasuk
perhitungan luas, pembalikan wenangan, dan persimpangan poligon.

DIBUTUHKAN ARGUMEN

Tidak ada.

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.

-A Amin_dist][satuan]
Lakukan analisis tetangga terdekat (NN) spasial: Tentukan tetangga terdekat dari
setiap titik dan laporkan jarak NN dan ID titik yang terlibat dalam setiap pasangan (ID
adalah nomor catatan input mulai dari 0). Menggunakan -A A untuk memusnahkan kumpulan data jadi
bahwa tidak ada jarak NN yang lebih rendah dari ambang batas min_dist. Dalam hal ini kita menulis
keluar (mungkin rata-rata) koordinat dan jarak NN dan ID titik yang diperbarui.
Angka poin negatif berarti poin asli diganti dengan rata-rata tertimbang
(nilai ID absolut memberikan ID dari ID titik asli pertama yang dimasukkan
rata-rata.)

-C Klip poligon ke wilayah peta, termasuk batas peta ke poligon sesuai kebutuhan.
Hasilnya adalah poligon tertutup (lihat -T untuk pemotongan sebagai gantinya). Memerlukan -R.

-D[+ffillet] [+sebuahamax] [+ dd maks] [+c|Ccmax] [+ dtkfakta]
Periksa duplikat di antara garis input atau poligon, atau, jika fillet diberikan melalui
+f, periksa apakah fitur input sudah ada di antara fitur-fitur di fillet. Kami
pertimbangkan kasus-kasus kecocokan eksak (jumlah dan koordinat yang sama) dan perkiraan
(jarak rata-rata antara titik terdekat dari dua fitur kurang dari ambang batas).
Kami juga mempertimbangkan bahwa beberapa fitur mungkin telah dibalik. Fitur dipertimbangkan
perkiraan kecocokan jika jarak minimumnya kurang dari d maks [0] (lihat UNIT) dan
kedekatan mereka (didefinisikan sebagai rasio antara jarak rata-rata antara
fitur dibagi dengan panjang rata-ratanya) kurang dari cmax [0.01]. Untuk setiap
duplikat ditemukan, catatan keluaran dimulai dengan satu huruf Y (sama persis) atau
~ (perkiraan cocok). Jika dua segmen yang cocok berbeda panjangnya lebih dari a
faktor 2 maka kita anggap duplikatnya sebagai subset (-) atau superset
(+). Akhirnya, kami juga mencatat jika dua garis adalah hasil dari pemisahan kontinu
garis melintasi Dateline (|). Untuk poligon kami juga mempertimbangkan pecahan
perbedaan wilayah; duplikat harus berbeda kurang dari amax [0.01]. Secara default,
kami menghitung pemisahan garis rata-rata. Menggunakan +Ccmin sebagai gantinya menghitung garis median
pemisahan dan karena itu nilai kedekatan yang kuat. Juga secara default kami mempertimbangkan semua
jarak antara titik-titik pada satu garis dengan garis lainnya. Menambahkan +p untuk membatasi perbandingan
ke titik-titik yang menonjol secara tegak lurus ke titik-titik pada garis lain (dan bukan garisnya)
perpanjangan).

-E+|- ]
Setel ulang kehandalan semua poligon agar sesuai dengan yang diberikan + (berlawanan arah jarum jam) atau -
(searah jarum jam). menyiratkan -Q+.

-F[l] Paksa data input menjadi poligon pada output, yaitu, tutup secara eksplisit jika tidak
sudah ditutup. Secara opsional, tambahkan l untuk memaksa geometri garis.

-Saya[e|i]
Tentukan lokasi perpotongan antara semua pasangan poligon. Menambahkan i untuk
hanya menghitung persilangan internal (yaitu, poligon berpotongan sendiri) atau e hanya untuk
menghitung eksternal (yaitu, antara paris poligon) crossover [Default adalah keduanya].

-Nfile p[+a] [+ halawal][+r][+z]
Tentukan apakah satu (atau semua, dengan +a) poin dari setiap fitur dalam data input adalah
di dalam salah satu poligon yang diberikan dalam file p. Jika di dalam, laporkan poligon mana
ini; ID poligon diambil dari nilai aspatial yang ditetapkan ke Z,
tajuk segmen (pertama -Z, kemudian -L dipindai), atau ditugaskan menjalankan
nomor yang diinisialisasi ke awal [0]. Secara default segmen input yang
ditemukan berada di dalam poligon ditulis ke stdout dengan ID poligon yang dikodekan dalam
tajuk segmen sebagai -ZID. Atau, tambahkan +r untuk hanya melaporkan poligon mana
berisi fitur atau +z agar ID ditambahkan sebagai kolom data tambahan pada output.
Segmen yang gagal berada di dalam poligon tidak akan ditulis. Jika lebih dari satu
poligon berisi segmen yang sama kita lewati skenario kedua (dan selanjutnya).

-Q[[-|+]*satuan*][+j][+l][+p]
Ukur luas semua poligon atau panjang segmen garis. Menggunakan -Q+jam untuk menambahkan
area ke setiap header segmen poligon [Default cukup tulis area ke stdout].
Untuk poligon kami juga menghitung lokasi centroid sedangkan untuk data garis kami menghitung
posisi titik tengah (setengah panjang). Tambahkan satuan jarak untuk memilih satuan
digunakan (lihat UNIT). Perhatikan bahwa area akan tergantung pada pengaturan saat ini dari
PROJ_ELLIPSOID; ini harus menjadi ellipsoid terbaru untuk mendapatkan hasil yang akurat. NS
centroid dihitung menggunakan rata-rata dari vektor Cartesian 3-D yang membentuk
simpul poligon, sedangkan luas diperoleh melalui proyeksi luas yang sama. Untuk
panjang garis yang dapat Anda tambahkan -|+ ke unit dan perhitungan akan menggunakan Flat
Algoritma Bumi atau Geodesi, masing-masing [Default adalah jarak lingkaran besar]. Oleh
default, kami menganggap poligon terbuka sebagai garis. Menambahkan +p untuk menutup poligon terbuka dan
jadi pertimbangkan semua input sebagai poligon, atau tambahkan +l untuk mempertimbangkan semua input sebagai garis,
bahkan jika ditutup.

-R[satuan]barat/timur/selatan/utara[/zmin/zmax][R]
barat, timur, selatan, dan utara tentukan wilayah yang diminati, dan Anda dapat menentukan
dalam derajat desimal atau dalam format [+-]dd:mm[:ss.xxx][W|E|S|N]. Menambahkan r jika lebih rendah
koordinat peta kiri dan kanan atas diberikan sebagai ganti w/e/s/n. keduanya
singkatan -Rg dan -Rd singkatan dari domain global (0/360 dan -180/+180 dalam bujur
masing-masing, dengan -90/+90 di garis lintang). Atau untuk pembuatan kisi, berikan
Rkodelon/lat/nx/ny, Di mana kode adalah kombinasi 2 karakter dari L, C, R (untuk kiri,
tengah, atau kanan) dan T, M, B untuk atas, tengah, atau bawah. misalnya, BL untuk kiri bawah.
Ini menunjukkan titik mana pada daerah persegi panjang yang lon/lat koordinat mengacu
ke, dan dimensi kisi nx dan ny dengan jarak kisi melalui -I digunakan untuk membuat
wilayah yang sesuai. Atau, tentukan nama file kisi yang ada
dan -R pengaturan (dan spasi grid, jika berlaku) disalin dari grid.
Menggunakan -Rsatuan mengharapkan koordinat yang diproyeksikan (Cartesian) yang kompatibel dengan yang dipilih -J
dan kami memproyeksikan terbalik untuk menentukan wilayah geografis persegi panjang yang sebenarnya. Untuk
pandangan perspektif (-p), opsional tambahkan /zmin/zmax. Dalam hal tampilan perspektif
(-p), rentang-z (zmin, zmax) dapat ditambahkan untuk menunjukkan dimensi ketiga. Ini
hanya perlu dilakukan saat menggunakan -Jz pilihan, bukan saat hanya menggunakan -p .
Dalam kasus terakhir, pandangan perspektif pesawat diplot, tanpa ketiga
dimensi. Klip poligon ke wilayah peta, termasuk batas peta ke poligon
sesuai kebutuhan. Hasilnya adalah poligon tertutup.

-Si|j|s|u
Pemrosesan spasial poligon. Pilih dari -Iya yang mengembalikan persimpangan
poligon (tertutup), -Su yang mengembalikan penyatuan poligon (tertutup), -H.H yang akan
membagi poligon yang mengangkangi Dateline, dan -Sj yang akan bergabung dengan poligon itu
dipisahkan oleh Dateline. Catatan: Hanya -H.H telah dilaksanakan.

-T[klippoligon]
Memotong poligon terhadap poligon tertentu yang diberikan, mungkin menghasilkan open
poligon. Jika tidak ada argumen yang diberikan kepada -T kami membuat poligon kliping dari -R yang
maka diperlukan. Perhatikan bahwa ketika -R kliping berlaku kami juga akan mencari
poligon dengan panjang 4 atau 5 yang sama persis dengan -R kliping poligon.

-V[tingkat] (lebih ...)
Pilih tingkat verbositas [c].

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

-bo[ncol][mengetik] (lebih ...)
Pilih keluaran biner asli. [Default sama dengan input].

-d[i|o]tidak ada data (lebih ...)
Ganti kolom input yang sama tidak ada data dengan NaN dan lakukan kebalikannya pada output.

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

-g[a]x|y|d|X|Y|D|[col]z[+|-]celah[atau] (lebih ...)
Tentukan kesenjangan data dan jeda baris.

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

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

-okerah[,...] (lebih ...)
Pilih kolom keluaran (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.

UNIT

Untuk satuan jarak peta, tambahkan satuan d untuk derajat busur, m untuk menit busur, dan s untuk busur
kedua, atau e untuk meteran [Awal], f untuk kaki, k untuk km, M untuk mil undang-undang, n untuk bahari
mil, dan u untuk survei kaki AS. Secara default kami menghitung jarak tersebut menggunakan bola
pendekatan dengan lingkaran besar. tambahkan - untuk jarak (atau unit tidak ada jarak adalah
diberikan) untuk melakukan perhitungan "Bumi Datar" (lebih cepat tetapi kurang akurat) atau tambahkan + untuk
melakukan perhitungan geodesik yang tepat (lebih lambat tetapi lebih akurat).

ASCII FORMAT KETEPATAN

Format output ASCII dari data numerik dikendalikan oleh parameter di gmt.conf
mengajukan. Bujur dan lintang diformat menurut FORMAT_GEO_OUT, sedangkan lainnya
nilai diformat menurut FORMAT_FLOAT_OUT. Ketahuilah bahwa format yang berlaku dapat
menyebabkan hilangnya presisi dalam output, yang dapat menyebabkan berbagai masalah hilir. Jika
Anda menemukan output tidak ditulis dengan presisi yang cukup, pertimbangkan untuk beralih ke biner
keluaran (-bo jika tersedia) atau tentukan lebih banyak desimal menggunakan pengaturan FORMAT_FLOAT_OUT.

CONTOH

Untuk mengubah semua baris dalam file multisegmen lines.txt menjadi poligon tertutup, jalankan

gmt spasial lines.txt -F > polygons.txt

Untuk menghitung luas semua poligon geografis dalam file multisegmen polygons.txt, jalankan

gmt spasial polygons.txt -Q > area.txt

Data yang sama, tetapi sekarang arahkan semua poligon berlawanan arah jarum jam dan tuliskan areanya ke
header segmen, jalankan

gmt spasial polygons.txt -Q+h -E+ > area.txt

Untuk menentukan perpotongan antara poligon A.txt dan B.txt, jalankan

gmt spasial A.txt B.txt -Ce > crossovers.txt

Gunakan gmtspatialgmt online menggunakan layanan onworks.net