mapprojectgmt - Online di Cloud

PROGRAM:

NAMA

mapproject - Lakukan transformasi peta maju dan terbalik, konversi datum, dan geodesi

RINGKASAN

proyek peta [ tabel ] parameter wilayah [ b|B|f|F|o|O[kesepian0/lat0] ] [ [dx/dy] ] [ c|i|p ] [
[datums] ] [ [satuan] ] [ [x0/y0/] [[+|-]satuan][+|-] ] [ ] [ line.xy[/[+|-]satuan][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]dari[/untuk] ] [ [tingkat] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o] ]

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

DESKRIPSI

proyek peta membaca (bujur, lintang) posisi dari tabel [atau masukan standar] dan
menghitung (x,y) koordinat menggunakan proyeksi dan skala peta yang ditentukan. Opsional, itu
dapat membaca (x,y) posisi dan menghitung (bujur, lintang) nilai melakukan kebalikannya
transformasi. Ini dapat digunakan untuk mentransformasikan titik-titik linier (x,y) yang diperoleh dengan mendigitalkan
peta proyeksi yang diketahui ke koordinat geografis. Mungkin juga menghitung jarak sepanjang
trek, ke titik tetap, atau pendekatan terdekat ke garis. Akhirnya bisa digunakan untuk perform
berbagai konversi datum. Bidang data tambahan diizinkan setelah 2 kolom pertama
yang harus memiliki (bujur,lintang) atau (x,y). Lihat opsi -: tentang cara membaca
(lintang, bujur) file.

DIBUTUHKAN ARGUMEN

-Jparameter (lebih ...)
Pilih proyeksi peta.

-R[satuan]xmin/xmax/ymin/ymax[R] (lebih ...)
Tentukan wilayah yang diminati. Kasus khusus untuk proyeksi UTM: Jika -C digunakan
dan -R tidak diberikan maka wilayah diatur bertepatan dengan zona UTM yang diberikan sehingga
untuk mempertahankan solusi elipsoidal penuh (Lihat PEMBATASAN untuk lebih lanjut
informasi).

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.

-Ab|B|f|F|o|O[kesepian0/lat0]
-Af menghitung azimuth (maju) dari titik tetap panjang/lat ke setiap titik data.
penggunaan -Jauh untuk mendapatkan back-azimuth dari titik data ke titik tetap. Menggunakan -Ao untuk mendapatkan
orientasi (-90/90) daripada azimut (0/360). Huruf besar F, B or O akan
mengkonversi dari geodesi ke geosentris lintang dan memperkirakan azimut geodesi
(dengan asumsi ellipsoid saat ini bukan bola). Jika tidak ada titik tetap yang diberikan maka
kami menghitung azimuth (atau back-azimuth) dari titik sebelumnya.

-C[dx/dy]
Setel pusat koordinat yang diproyeksikan berada di pusat proyeksi peta [Awalnya adalah
pojok kiri bawah]. Secara opsional, tambahkan offset di unit yang diproyeksikan untuk ditambahkan (atau
dikurangi bila -I diatur) ke (dari) koordinat yang diproyeksikan, seperti false
timur dan utara untuk zona proyeksi tertentu [0/0]. Satuan yang digunakan untuk
offset adalah satuan jarak plot yang berlaku (lihat PROJ_LENGTH_UNIT) kecuali -F is
digunakan, dalam hal ini offset dalam meter.

-Dc|i|hal
Ganti sementara PROJ_LENGTH_UNIT dan gunakan c (cm), i (inci), atau p (poin)
sebagai gantinya. Tidak dapat digunakan dengan -F.

-E[datums]
Mengkonversi dari geodetik (lon, lat, tinggi) ke Earth Centered Earth Fixed (ECEF)
(x,y,z) koordinat (tambahkan -I untuk konversi terbalik). Tambahkan ID datum (lihat -Qd)
atau berikan elipsoid:dx,dy,dz dimana elipsoid mungkin ID ellipsoid (lihat -Qe) atau
diberikan sebagai a[,*inv_f*], di mana a adalah sumbu semi-mayor dan inv_f adalah kebalikannya
merata (0 jika dihilangkan). Jika datums adalah - atau tidak diberikan, kami menganggap WGS-84.

-F[satuan]
Penskalaan paksa 1:1, yaitu, keluaran (atau masukan, lihat -I) data dalam proyeksi aktual
meter. Untuk menentukan unit lain, tambahkan unit yang diinginkan (lihat UNIT). Tanpa -F,
keluaran (atau masukan, lihat -I) berada dalam unit yang ditentukan oleh PROJ_LENGTH_UNIT (tetapi
melihat -D).

-G[x0/y0/][[+|-]satuan][+|-]
Hitung jarak di sepanjang trek or ke titik opsional yang disetel dengan -Gx0/y0. Menambahkan
satuan jarak (lihat UNIT), termasuk c (Jarak Cartesian menggunakan input
koordinat) atau C (Jarak Cartesian menggunakan koordinat proyeksi). NS C satuan
membutuhkan -R dan -J untuk diatur. Tanpa titik tetap yang diberikan, kami menghitung cumulate
jarak di sepanjang lintasan. Menambahkan - untuk mendapatkan jarak tambahan antara berturut-turut
poin. Menambahkan + untuk menentukan titik ke-2 melalui dua kolom tambahan di file input.

-I Lakukan transformasi Invers, yaitu, dapatkan (bujur,lintang) dari data (x,y).

-Lline.xy[/[+|-]satuan] [+]
Tentukan jarak terpendek dari titik data input ke garis yang diberikan dalam
file multisegmen ASCII line.xy. Jarak dan koordinat
titik terdekat akan ditambahkan ke output sebagai tiga kolom baru. Tambahkan
satuan jarak (lihat UNIT), termasuk c (Jarak Cartesian menggunakan koordinat input)
or C (Jarak Cartesian menggunakan koordinat proyeksi). NS C unit membutuhkan -R dan
-J untuk diatur. Terakhir, tambahkan + untuk melaporkan id segmen garis dan pecahan
nomor titik bukan lon/lat dari titik terdekat.

-N[a|c|g|m]
Konversi dari garis lintang geodetik (menggunakan ellipsoid saat ini; lihat PROJ_ELLIPSOID)
ke salah satu dari empat garis lintang bantu yang berbeda (garis bujur tidak terpengaruh). Memilih
dari authalik, ctidak resmi, geosentris, dan mlintang eridional [geosentris]. Menggunakan -I
untuk mengkonversi dari garis lintang bantu ke garis lintang geodetik.

-Q[d|e Daftar semua parameter proyeksi. Untuk hanya membuat daftar datum, gunakan -Qd. Untuk hanya daftar
elipsoid, gunakan -Qe.

-S Menekan poin yang jatuh di luar wilayah tersebut.

-NS]dari[/untuk]
Konversi koordinat antar datum dari dan untuk menggunakan standar Molodensky
transformasi. Menggunakan -NS jika kolom input ke-3 memiliki ketinggian di atas ellipsoid [Default
mengasumsikan tinggi = 0, yaitu, pada ellipsoid]. Tentukan datum menggunakan ID datum (lihat
-Qd) atau memberikan elipsoid:dx,dy,dz dimana elipsoid mungkin ID ellipsoid (lihat -Qe) atau
diberikan sebagai a[,*inv_f*], di mana a adalah sumbu semi-mayor dan inv_f adalah kebalikannya
merata (0 jika dihilangkan). Jika datums adalah - atau tidak diberikan, kami menganggap WGS-84. -T mungkin
digunakan bersama dengan -R -J untuk mengubah datum sebelum koordinat proyeksi
(Menambahkan -I untuk menerapkan konversi datum setelah proyeksi terbalik). Pastikan bahwa
pengaturan PROJ_ELLIPSOID sudah benar untuk kasus Anda.

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

-W[w|h]
Mencetak lebar dan tinggi peta pada keluaran standar. Tidak ada file input yang dibaca. Hanya untuk
keluaran lebar atau tinggi, tambahkan w or h, masing-masing. Satuan dari
dimensi dapat diubah melalui -D.

-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).

-p[x|y|z]azim/tinggi[/tingkat z] [+wkesepian0/lat0[/z0]]] [+vx0/y0] (lebih ...)
Pilih tampilan perspektif.

-S[kerah][a|r] (lebih ...)
Mengatur penanganan catatan NaN.

-:[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 mengonversi koordinat UTM dalam meter ke lokasi geografis, diberikan file utm.txt dan
mengetahui zona UTM (dan zona atau belahan bumi), coba

proyek peta gmt utm.txt -Ju+11/1:1 -C -I -F

Untuk mengubah file dengan (bujur,lintang) menjadi (x,y) posisi dalam cm pada Mercator
grid untuk skala tertentu 0.5 cm per derajat, run

gmt mapproject lonlatfile -R20/50/12/25 -Jm0.5c > xyfile

Untuk mengubah beberapa file 2-kolom, biner, presisi ganda dengan (lintang, bujur)
ke dalam posisi (x,y) dalam inci pada kisi Transverse Mercator (bujur tengah 75W) untuk
skala = 1:500000 dan tekan titik-titik yang akan berada di luar area peta, jalankan

gmt mapproject track.* -R-80/-70/20/40 -Jt-75/1:500000 -: -S -Di -bo -bi2 > tmfile.b

Untuk mengonversi koordinat geodetik (lon, lat, height) di file old.dat dari NAD27
CONUS datum (Datum ID 131 yang menggunakan ellipsoid Clarke-1866) ke WGS 84, jalankan

gmt mapproject old.dat -Th131 > new.dat

Untuk menghitung jarak terdekat (dalam km) antara setiap titik dalam file input quakes.dat
dan segmen garis yang diberikan dalam file ASCII multisegmen, coast.xy, run

gmt mapproject quakes.dat -Lcoastline.xy/k > quake_dist.dat

BATASAN

Wilayah input persegi panjang diatur dengan -R secara umum akan dipetakan menjadi non-persegi panjang
kisi. Kecuali -C diatur, titik paling kiri pada kisi ini memiliki nilai x = 0.0, dan
titik paling bawah akan memiliki nilai y = 0.0. Jadi, sebelum Anda mendigitalkan peta, jalankan yang ekstrem
koordinat peta melalui proyek peta menggunakan skala yang sesuai dan lihat apa (x,y) nilai
mereka dipetakan. Gunakan nilai-nilai ini saat menyiapkan untuk digitalisasi agar
transformasi terbalik bekerja dengan benar, atau sebagai alternatif, gunakan Wow untuk menskala dan menggeser
(x,y) nilai sebelum transformasi.

Untuk beberapa proyeksi, solusi sferis dapat digunakan meskipun pengguna telah memilih
elips. Ini terjadi ketika pengguna -R pengaturan menyiratkan wilayah yang melebihi domain
di mana ekspansi deret ellipsoidal valid. Berikut syarat-syaratnya: (1)
Kerucut Konformal Lambert (-JL) dan Albers Equal-Area (-JB) akan menggunakan solusi sferis
ketika skala peta melebihi 1.0E7. (2) Mercator Transversal (-JT) dan UTM (-JU) akan akan
gunakan solusi bola ketika batas barat atau timur diberikan dalam -R lebih dari
10 derajat dari meridian pusat, dan (3) sama untuk Cassini (-JC) tetapi dengan batas
hanya 4 derajat.

ELLIPSOID DAN SPHEROID

GMT akan menggunakan rumus ellipsoidal jika diimplementasikan dan pengguna telah memilih satu
ellipsoid sebagai bentuk referensi (lihat PROJ_ELLIPSOID). Pengguna perlu menyadari beberapa
potensi jebakan: (1) Untuk beberapa proyeksi, seperti Transverse Mercator, Albers, dan
Kerucut konformal Lambert kami menggunakan ekspresi ellipsoidal ketika area yang dipetakan adalah
kecil, dan beralih ke ekspresi bola (dan mengganti alat bantu yang sesuai
lintang) untuk peta yang lebih besar. Rumus ellipsoidal digunakan sebagai berikut: (a) Transversal
Mercator: Ketika semua titik berada dalam 10 derajat meridian pusat, (b) Proyeksi kerucut
ketika rentang longitudinal kurang dari 90 derajat, (c) proyeksi Cassini ketika semua titik
berada dalam jarak 4 derajat dari meridian pusat. (2) Saat Anda mencoba mencocokkan beberapa sejarah
data (misalnya, koordinat yang diperoleh dengan proyeksi tertentu dan referensi tertentu
ellipsoid) Anda mungkin menemukan bahwa GMT memberikan hasil yang sedikit berbeda. Satu kemungkinan
sumber ketidakcocokan ini adalah bahwa perhitungan yang lebih lama sering menggunakan angka yang kurang signifikan. Untuk
contoh, contoh Snyder sering menggunakan ellipsoid Clarke 1866 (didefinisikan olehnya sebagai memiliki
a merata f = 1/294.98). Dari f kita mendapatkan kuadrat eksentrisitas menjadi 0.00676862818
(inilah yang digunakan GMT), sementara Snyder membulatkan dan menggunakan 0.00676866. Perbedaan ini dapat
memberikan perbedaan beberapa puluh cm. Jika Anda perlu mereproduksi koordinat yang diproyeksikan
dengan eksentrisitas yang sedikit berbeda ini, Anda harus menentukan ellipsoid Anda sendiri dengan
parameter yang sama seperti Clarke 1866, tetapi dengan f = 1/294.97861076. Juga, ketahuilah bahwa lebih tua
data dapat direferensikan ke datum yang berbeda, dan kecuali Anda tahu datum mana yang digunakan dan
mengonversi semua data ke datum umum, Anda mungkin mengalami ketidakcocokan puluhan hingga ratusan
meter. (3) Akhirnya, ketahuilah bahwa PROJ_SCALE_FACTOR memiliki nilai default tertentu untuk beberapa
proyeksi sehingga Anda mungkin harus mengganti pengaturan agar sesuai dengan hasil yang dihasilkan
dengan pengaturan lainnya.

Gunakan mapprojectgmt online menggunakan layanan onworks.net