img2grdgmt - Online di Cloud

PROGRAM:

NAMA

img2grd - Ekstrak subset file img dalam format Mercator atau Geografis

RINGKASAN

img2grd file img file grd wilayah mengetik [ ] [ [minlat/makslat] ] [ ] [ menit ] [ ] [ navigasi
] [ [skala] ] [ [tingkat] ] [ maxlon ] [ -n]

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

DESKRIPSI

img2grd membaca file format img, mengekstrak subset, dan menulisnya ke file grid. NS -M
opsi menentukan apakah proyeksi Spherical Mercator dari file img adalah
dipertahankan atau jika kisi Geografis harus ditulis dengan membatalkan proyeksi Mercator. Jika
grid geografis dipilih, Anda juga dapat meminta pengambilan sampel ulang ke yang tepat -R diberikan.

DIBUTUHKAN ARGUMEN

file img
File format img Mercator seperti gravitasi laut atau bidang topografi dasar laut
diperkirakan dari data altimeter satelit oleh Sandwell dan Smith. Jika pengguna telah mengatur
variabel lingkungan $GMT_DATADIR, kemudian img2grd akan mencoba mencari file img in
$GMT_DATADIR; lain itu akan mencoba untuk membuka file img langsung.

-Gfile grd
file grd adalah nama file grid keluaran.

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

OPSIONAL ARGUMEN

-C Atur koordinat Mercator x dan y relatif terhadap pusat proyeksi [Awalnya adalah
relatif terhadap sudut kiri bawah grid]. Memerlukan -M.

-D[minlat/makslat]
Gunakan rentang garis lintang yang diperluas -80.738/+80.738. Atau, tambahkan
minlat/makslat sebagai tingkat lintang dari file img input. [Defaultnya adalah
-72.006/72.006]. Biasanya tidak diperlukan karena kita dapat menentukan sejauh mana dari
pemeriksaan ukuran file.

-E Dapat digunakan saat -M tidak diatur untuk memaksa kisi terakhir memiliki hal yang sama persis
wilayah seperti yang diminta dengan -R. Secara default, wilayah akhir adalah proyeksi langsung dari
wilayah Mercator asli dan biasanya akan sedikit melampaui
rentang garis lintang yang diminta, dan selanjutnya kenaikan kisi di garis lintang tidak
cocok dengan kenaikan bujur. Namun, pengambilan sampel ulang tambahan menghasilkan sedikit
kesalahan interpolasi dan hanya boleh digunakan jika kisi keluaran harus cocok dengan
wilayah yang diminta dan memiliki x_inc = y_inc. Dalam hal ini wilayah yang ditetapkan oleh -R harus
diberikan dalam kelipatan dari kenaikan (.eg, -R0/45/45/72).

-I Menunjukkan menit sebagai lebar piksel img masukan dalam bujur menit.
[Awalnya adalah 2.0]. Biasanya tidak diperlukan karena kami dapat menentukan ukuran piksel dari
pemeriksaan ukuran.

-M Keluarkan kisi Spherical Mercator [Default adalah kisi lon/lat geografis]. NS
Proyeksi Mercator Bulat dari file img dipertahankan, sehingga wilayah -R
ditetapkan oleh pengguna sedikit dimodifikasi; wilayah yang dimodifikasi sesuai dengan tepi
piksel [atau kelompok navigasi piksel]. Header file grid diatur sehingga x dan
panjang sumbu y mewakili jarak dari tepi barat dan selatan gambar,
diukur dalam unit default pengguna, dengan -Jm1 dan disesuaikan -R. Dengan mengatur
kegagalan PROJ_ ELIPSOID = Sphere, pengguna dapat membuat overlay dengan yang disesuaikan -R
sehingga mereka cocok. Lihat CONTOH di bawah. yang disesuaikan -R juga ditulis di grid
komentar header, sehingga dapat ditemukan nanti. Lihat -C untuk mengatur koordinat relatif terhadap
pusat proyeksi.

-Nnavigasi Rata-rata nilai dalam piksel img masukan menjadi navigasi by navigasi persegi, dan buat
satu piksel keluaran untuk setiap kotak tersebut. Jika digunakan dengan -T3 itu akan melaporkan rata-rata
batasan antara 0 dan 1. Jika digunakan dengan -T2 outputnya akan menjadi nilai data rata-rata
atau NaN menurut apakah kendala rata-rata > 0.5. navigasi harus dibagi rata
ke dalam dimensi imgfile dalam piksel. [Bawaan 1 tidak rata-rata].

-S[skala]
Kalikan nilai file img dengan skala sebelum disimpan dalam file grid. [Defaultnya adalah
1.0]. Untuk file img terbaru: file img topo disimpan dalam meter (dikoreksi) [-S1];
file gravitasi udara bebas dalam mGal*10 [-S0.1 untuk mendapatkan mGal]; file defleksi vertikal di
mikro-radian*10 [-S0.1 untuk mendapatkan radian mikro], file gradien gravitasi vertikal di
Eotvo*50 [-S0.02 untuk mendapatkan Eotvos, atau -S0.002 untuk mendapatkan mGal/km]). Jika tidak skala diberikan
kami mencoba menentukan skala dengan memeriksa nama file untuk petunjuk.

-Tmengetik mengetik menangani pengkodean informasi kendala. mengetik = 0 menunjukkan bahwa tidak
informasi tersebut dikodekan dalam file img (digunakan untuk versi pra-1995 dari
data gravitasi) dan mendapatkan semua data. mengetik > 0 menunjukkan bahwa informasi kendala adalah
versi file img yang dikodekan (1995 dan yang lebih baru (saat ini)) sehingga seseorang dapat
menghasilkan file grid sebagai berikut: -T1 mendapatkan nilai data di semua titik, -T2 mendapat data
nilai pada titik terbatas dan NaN pada titik interpolasi; -T3 mendapat 1 pada
titik dibatasi dan 0 pada titik interpolasi [Default adalah 1].

-V[tingkat] (lebih ...)
Pilih tingkat verbositas [c]. Sangat direkomendasikan di sini, karena sangat membantu untuk dilihat
bagaimana koordinat disesuaikan.

-Wmaxlon
Menunjukkan maxlon sebagai batas bujur maksimum dari file img input. Versi
sejak tahun 1995 telah maxlon = 360.0, sementara beberapa file sebelumnya memiliki maxlon = 390.0.
[Awalnya adalah 360.0].

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

GEOGRAFIS CONTOH

Grafik -M opsi harus dikecualikan jika Anda membutuhkan kisi keluaran berada di geografis
koordinat. Untuk mengekstrak data di wilayah -R-40/40/-70/-30 dari world_grav.img.7.2 dan
memproyeksikan ulang untuk menghasilkan koordinat geografis, Anda dapat mencoba

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Karena jarak lintang dalam file img sama dalam unit Mercator,
kisi yang dihasilkan tidak akan cocok dengan yang ditentukan -R persis, dan jarak lintang tidak akan
sama dengan jarak bujur. Jika Anda membutuhkan kecocokan yang tepat dengan Anda -R dan jarak yang sama
dalam bujur dan lintang, gunakan -E opsi:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

MERKATOR CONTOH

Karena file img berada dalam proyeksi Mercator, Anda TIDAK boleh mengekstrak kisi geografis
jika rencana Anda adalah membuat peta Mercator. Jika Anda melakukan itu, Anda mengakhiri proyeksi dan
memproyeksikan ulang grid, kehilangan detail panjang gelombang pendek. Lebih baik digunakan -M dan plot grid
menggunakan proyeksi linier dengan skala yang sama dengan proyeksi Mercator yang diinginkan (lihat GMT
Contoh 29). Untuk mengekstrak data di wilayah -R-40/40/-70/-30 dari world_grav.img.7.2, lari

gmt img2grd -M dunia_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Perhatikan bahwa -V opsi memberi tahu kami bahwa rentangnya disesuaikan dengan
-R-40/40/-70.0004681551/-29.9945810754. Kami juga dapat menggunakan info grd untuk menemukan bahwa file grid
header menunjukkan wilayahnya menjadi -R0/80/0/67.9666667 Ini adalah kisaran x,y yang akan kita dapatkan
dari proyeksi Spherical Mercator menggunakan -R-40/40/-70.0004681551/-29.9945810754 dan
-Jm1. Jadi, untuk mengambil ship.lonlatgrav dan menggunakannya untuk sampel merc_grav.nc, kita dapat melakukan ini:

gmt setel PROJ_ELLIPSOID Sphere

proyek peta gmt -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | proyek peta gmt \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > ship.lonlatgravsat

Disarankan untuk menggunakan metode proyeksi dan unproyeksi data di atas sedemikian rupa
aplikasi, karena hanya ada satu langkah interpolasi (dalam jalur grd). Jika satu
pertama mencoba mengonversi file grid ke lon,lat dan kemudian mengambil sampelnya, ada dua
langkah-langkah interpolasi (dalam konversi dan dalam pengambilan sampel).

Untuk membuat lon,lat grid dari grid di atas bisa kita gunakan

gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

Dalam beberapa kasus ini tidak akan mudah karena -R dalam dua sistem koordinat mungkin tidak sejajar
dengan baik. Ketika ini terjadi, kami juga dapat menggunakan (sebenarnya, mungkin selalu lebih baik untuk digunakan)

gmt grd2xyz merc_grav.nc | proyek peta gmt \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
permukaan gmt -R-40/40/-70/70 -I2m -Ggrav.nc

Untuk membuat peta Mercator dari wilayah di atas, misalkan nilai gmt.conf kami untuk
PROJ_LENGTH_UNIT adalah inci. Kemudian karena file merc_grav.nc di atas diproyeksikan dengan -Jm1i
itu adalah 80 inci lebar. Kita dapat membuat peta dengan lebar 8 inci dengan menggunakan -Jx0.1i pada program peta apa pun
diterapkan ke kisi ini (misalnya, kontur, gambar, tampilan), dan kemudian untuk overlay yang
bekerja di lon, lat (misalnya, psxy, pscoast) kita dapat menggunakan penyesuaian di atas -R dan -Jm0.1 untuk mendapatkan
dua sistem untuk mencocokkan.

Namun, kita bisa lebih pintar dari ini. Menyadari bahwa file img input memiliki piksel 2.0
menit lebar (atau memeriksa nx dan ny dengan grdinfo merc_grav.nc) kami menyadari bahwa
merc_grav.nc menggunakan resolusi penuh dari file img dan memiliki 2400 x 2039 piksel, dan
pada lebar 8 inci ini adalah 300 piksel per inci. Kami memutuskan kami tidak membutuhkan sebanyak itu dan kami
akan puas dengan 100 piksel per inci, jadi kami ingin rata-rata data menjadi 3 kali 3
kotak. (Jika kita menginginkan plot kontur, kita mungkin akan memilih rata-rata data lebih banyak
(misalnya, 6 dengan 6) untuk mendapatkan kontur yang halus.) Karena 2039 tidak habis dibagi 3, kita akan mendapatkan a
berbeda disesuaikan -R kali ini:

gmt img2grd -M dunia_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

Kali ini kami menemukan wilayah yang disesuaikan adalah -R-40/40/-70.023256525/-29.9368261101 dan
outputnya adalah 800 x 601 piksel, ukuran yang lebih baik bagi kami. Sekarang kita bisa membuat buatan
file iluminasi untuk penggunaan ini gradien:

gmt grdgradien merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

dan jika kita juga memiliki file CPT bernama "grav.cpt" kita dapat membuat peta relief berbayang warna
seperti ini:

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

Misalkan Anda hanya ingin mendapatkan nilai data yang dibatasi dari file img, dalam lat/lon
koordinat. Lalu lari img2grd pada pengatur terkenal. Pengatur ini menawarkan bantuan hukum kepada traderapabila trader berselisih dengan broker yang terdaftar dengan mereka. -T2 opsi, gunakan grd2xyz untuk membuang nilainya, pipa
melalui grep -v NaN untuk menghilangkan NaN, dan pipa melalui proyek peta dengan kebalikannya
proyeksi seperti di atas.

Gunakan img2grdgmt online menggunakan layanan onworks.net