img2grdgmt - Bulutta Çevrimiçi

Program:

ADI

img2grd - img dosyasının alt kümesini Mercator veya Coğrafi biçimde çıkar

SİNOPSİS

img2grd img dosyası grd dosyası bölge tip [ ] [ [minlat/makslat] ] [ ] [ dakika ] [ ] [ gezinme
] [ [ölçek] ] [seviye] ] [ makslon ] [ -n]

Not: Seçenek bayrağı ve ilişkili bağımsız değişkenler arasında boşluk bırakılamaz.

TANIM

img2grd img biçimindeki bir dosyayı okur, bir alt kümeyi ayıklar ve bir ızgara dosyasına yazar. NS -M
seçeneği, img dosyasının Spherical Mercator projeksiyonunun olup olmadığını belirler.
korunmuşsa veya Mercator projeksiyonu geri alınarak bir Coğrafi ızgara yazılmalıdır. Eğer
coğrafi ızgara seçildi, ayrıca kesin olarak yeniden örnekleme talep edebilirsiniz. -R Verilen.

GEREKLİ ARGÜMANLAR

img dosyası
Deniz yerçekimi veya deniz tabanı topografya alanları gibi bir Mercator img formatı dosyası
Sandwell ve Smith tarafından uydu altimetre verilerinden tahmin edilmiştir. Kullanıcı ayarladıysa
bir ortam değişkeni $GMT_DATADIR, Daha sonra img2grd bulmaya çalışacak img dosyası in
$GMT_DATADIR; yoksa açmaya çalışacak img dosyası direkt olarak.

-Ggrd dosyası
grd dosyası çıktı ızgara dosyasının adıdır.

-R[birim]batı/doğu/güney/kuzey[/zmin/zmaks][R]
batı, doğu, güney, ve kuzey ilgilendiğiniz bölgeyi belirtin ve belirtebilirsiniz
ondalık derece veya [+-]dd:mm[:ss.xxx][W|E|S|N] biçiminde. Ekle r daha düşükse
w/e/s/n yerine sol ve sağ üst harita koordinatları verilmiştir. İki
stenografi -Rg ve -Rd küresel alanı temsil eder (0/360 ve -180/+180 boylamda)
sırasıyla, enlemde -90/+90 ile). Alternatif olarak ızgara oluşturma için,
Rkodlon/lat/nx/ny, Burada kod L, C, R'nin 2 karakterlik bir birleşimidir (sol için,
orta veya sağ) ve üst, orta veya alt için T, M, B. örneğin, sol alt için BL.
Bu, dikdörtgen bir bölge üzerindeki hangi noktanın lon/lat koordinat ifade eder
ve ızgara boyutları nx ve ny ızgara aralıkları ile -I oluşturmak için kullanılır
karşılık gelen bölge. Alternatif olarak, mevcut bir ızgara dosyasının adını belirtin
ve -R ayarlar (ve varsa ızgara aralığı) ızgaradan kopyalanır.
kullanma -Rbirim seçilen ile uyumlu öngörülen (Kartezyen) koordinatları bekler -J
ve gerçek dikdörtgen coğrafi bölgeyi belirlemek için ters olarak projeksiyon yapıyoruz. İçin
perspektif görünüş (-p), isteğe bağlı olarak ekle /zmin/zmaks. Perspektif görünüm durumunda
(-p), bir z aralığı (zmin, zmaks) üçüncü boyutu belirtmek için eklenebilir. Bu
yalnızca kullanılırken yapılmalıdır. -Jz seçeneği, yalnızca -p seçeneği.
İkinci durumda, üçüncü bir nokta olmadan, düzlemin bir perspektif görünümü çizilir.
boyut.

İSTEĞE BAĞLI ARGÜMANLAR

-C x ve y Mercator koordinatlarını projeksiyon merkezine göre ayarlayın [Varsayılan
ızgaranın sol alt köşesine göre]. Gereklilikler -M.

-NS[minlat/makslat]
Genişletilmiş enlem aralığı -80.738/+80.738'i kullanın. Alternatif olarak, ekle
minlat/makslat giriş img dosyasının enlem kapsamı olarak. [Varsayılan
-72.006/72.006]. kapsamını belirleyebileceğimiz için genellikle gerekli değildir.
dosya boyutunun incelenmesi.

-E Ne zaman kullanılabilir -M son ızgarayı tam olarak aynı olmaya zorlamak için ayarlanmadı
ile talep edilen bölge -R. Varsayılan olarak, son bölge,
orijinal Mercator bölgesidir ve tipik olarak bölgenin biraz ötesine uzanır.
istenen enlem aralığı ve ayrıca enlemdeki ızgara artışı
boylam artışını eşleştirin. Ancak, ekstra yeniden örnekleme küçük
enterpolasyon hataları ve yalnızca çıktı ızgarasının aşağıdakilerle eşleşmesi gerekiyorsa kullanılmalıdır.
istenen bölge ve x_inc = y_inc var. Bu durumda tarafından belirlenen bölge -R olmalı
artışın katları olarak verilir (örn. -R0/45/45/72).

-I belirtmek dakika dakika boylam cinsinden bir girdi img pikselinin genişliği olarak.
[Varsayılan 2.0'dır]. Piksel boyutunu şuradan belirleyebileceğimiz için genellikle gerekli değildir.
boyutunun incelenmesi.

-M Bir Küresel Merkatör ızgarası çıktısı alın [Varsayılan bir coğrafi boylam/enlem ızgarasıdır]. NS
img dosyasının küresel Mercator projeksiyonu korunur, böylece bölge -R
kullanıcı tarafından ayarlanan biraz değiştirilir; değiştirilmiş bölge kenarlara karşılık gelir
piksel sayısı [veya gezinme piksel]. Izgara dosyası başlığı, x ve
y ekseni uzunlukları görüntünün batı ve güney kenarlarından olan mesafeyi temsil eder,
kullanıcı varsayılan birimlerinde ölçülür, -Jm1 ve ayarlanmış -R. ayarlayarak
varsayılan PROJ_ ELİPSOİD = Küre, kullanıcı ayarlanan ile bindirmeler yapabilir -R
eşleşsinler diye. Görmek ÖRNEKLER aşağıda. ayarlanmış -R ızgaraya da yazılır
başlık açıklaması, böylece daha sonra bulunabilir. Görmek -C göre koordinatları ayarlamak için
projeksiyon merkezi.

-Ngezinme Giriş img piksellerindeki değerlerin ortalamasını alın gezinme by gezinme kareler ve oluşturmak
bu tür her kare için bir çıktı pikseli. ile kullanılırsa -T3 bir ortalama rapor edecek
0 ile 1 arasında kısıtlama -T2 çıktı ortalama veri değeri olacaktır
veya ortalama kısıtlamanın > 0.5 olup olmadığına göre NaN. gezinme eşit olarak bölmeli
piksel cinsinden imgfile boyutlarına. [Varsayılan 1 ortalama yok].

-S[ölçek]
img dosyası değerlerini şu şekilde çarpın: ölçek ızgara dosyasına kaydetmeden önce. [Varsayılan
1.0]. En son img dosyaları için: img topo dosyaları (düzeltilmiş) metrelerde saklanır [-S1];
mGal*10'da serbest hava yerçekimi dosyaları [-S0.1 mGal almak için]; dikey sapma dosyaları
mikro radyan*10 [-S0.1 mikro radyan elde etmek için], dikey yerçekimi gradyan dosyaları
Eotvolar*50 [-SEotvos almak için 0.02 veya -SmGal/km elde etmek için 0.002]). Eğer hayırsa ölçek verilmiş
ipuçları için dosya adını inceleyerek ölçeği belirlemeye çalışıyoruz.

-Ttip tip kısıtlama bilgilerinin kodlamasını yönetir. tip = 0 hayır olduğunu gösterir
bu tür bilgiler img dosyasında kodlanmıştır (1995 öncesi sürümleri için kullanılır).
yerçekimi verileri) ve tüm verileri alır. tip > 0, kısıtlama bilgisinin
kodlanmış (1995 ve img dosyalarının sonraki (geçerli) sürümleri), böylece
aşağıdaki gibi bir ızgara dosyası oluşturun: -T1 tüm noktalarda veri değerlerini alır, -T2 veri alır
kısıtlı noktalarda değerler ve enterpolasyonlu noktalarda NaN; -T3 1 alır
sınırlı noktalar ve enterpolasyonlu noktalarda 0 [Varsayılan 1'dir].

-V[seviye] (Daha ...)
Ayrıntı düzeyini seçin [c]. Görmek faydalı olduğu için özellikle burada tavsiye edilir.
koordinatlar nasıl ayarlanır.

-Wmakslon
belirtmek makslon giriş img dosyasının maksimum boylam kapsamı olarak. Sürümler
1995'ten beri var makslon = 360.0, daha önceki bazı dosyalar makslon = 390.0.
[Varsayılan 360.0'dır].

-^ or sadece -
Komutun sözdizimi hakkında kısa bir mesaj yazdırın, ardından çıkar (NOT: Windows'ta
sadece kullan -).

-+ or sadece +
Açıklama da dahil olmak üzere kapsamlı bir kullanım (yardım) mesajı yazdırın.
modüle özgü seçenek (ancak GMT ortak seçenekleri değil), ardından çıkar.

-? or yok hayır argümanlar
Seçeneklerin açıklamasını da içeren eksiksiz bir kullanım (yardım) mesajı yazdırın, ardından
çıkışlar.

--versiyon
GMT sürümünü yazdırın ve çıkın.

--gösteri-veri dizini
GMT paylaşım dizininin tam yolunu yazdırın ve çıkın.

COĞRAFİ ÖRNEKLER

The -M çıktı ızgarasının coğrafi olarak olması gerekiyorsa, seçenek hariç tutulmalıdır.
koordinatlar. Bölgedeki verileri çıkarmak için -R-40/40/-70/-30 world_grav.img.7.2 ve
coğrafi koordinatları elde etmek için yeniden projelendirin, deneyebilirsiniz

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

img dosyasındaki enlem aralığı Mercator birimlerinde eşit uzaklıkta olduğundan,
sonuçtaki ızgara belirtilenle eşleşmeyecek -R tam olarak ve enlem aralığı
boylam aralığına eşittir. Sizinle tam bir eşleşmeye ihtiyacınız varsa -R ve aynı aralık
boylam ve enlemde, -E seçenek:

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

MERKATOR ÖRNEKLER

İmg dosyaları bir Mercator projeksiyonunda olduğundan, bir coğrafi ızgara ÇIKARMAMALISINIZ
planınız bir Mercator haritası yapmaksa. Bunu yaptıysanız, projelendirmeyi sonlandırdınız ve
kısa dalga boyu detayını kaybederek ızgarayı yeniden yansıtın. kullanmak daha iyi -M ve ızgarayı çiz
istenen Mercator projeksiyonu ile aynı ölçekte doğrusal bir projeksiyon kullanarak (bkz. GMT
Örnek 29). Bölgedeki verileri çıkarmak için -R-40/40/-70/-30 world_grav.img.7.2, koşmak

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

Unutmayın -V seçeneği bize aralığın şu şekilde ayarlandığını söyler:
-R-40/40/-70.0004681551/-29.9945810754. Biz de kullanabiliriz grdin bilgisi ızgara dosyasını bulmak için
başlık bölgesini gösterir -R0/80/0/67.9666667 Bu x,y aralığını alacağız
kullanarak bir Küresel Mercator projeksiyonundan -R-40/40/-70.0004681551/-29.9945810754 ve
-Jm1. Böylece, ship.lonlatgrav'ı almak ve onu merc_grav.nc'yi örneklemek için kullanmak için şunu yapabiliriz:

gmt seti PROJ_ELLIPSOID Küre

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

Bu tür durumlarda verileri yansıtmak ve projeden çıkarmak için yukarıdaki yöntemin kullanılması tavsiye edilir.
bir uygulamadır, çünkü o zaman yalnızca bir enterpolasyon adımı vardır (içinde grdtrack). Eğer biri
önce grid dosyasını lon,lat'a dönüştürmeye çalışır ve sonra onu örnekleyin, iki tane var
enterpolasyon adımları (dönüşümde ve örneklemede).

Yukarıdaki ızgaradan bir lon,lat ızgarası yapmak için kullanabiliriz

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

Bazı durumlarda bu kolay olmayacak, çünkü -R iki koordinat sisteminde hizalanmayabilir
kuyu. Bu olduğunda biz de kullanabiliriz (aslında kullanmak her zaman daha iyi olabilir)

gmt grd2xyz merc_grav.nc | gmt harita projesi \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
gmt yüzey -R-40/40/-70/70 -I2m -Ggrav.nc

Yukarıdaki bölgenin bir Mercator haritasını yapmak için, gmt.conf değerimizi varsayalım.
PROJ_LENGTH_UNIT inçtir. O zaman yukarıdaki merc_grav.nc dosyası ile yansıtıldığından beri -Jm1i
80 inç genişliğindedir. kullanarak 8 inç genişliğinde bir harita yapabiliriz. -Jx0.1i herhangi bir harita programında
bu ızgaraya uygulanan (örn. grdcontour, küçük resim, grd görünümü) ve ardından bindirmeler için
lon,lat'ta çalışın (örn. PSxy, pcoast) yukarıda ayarlanmış olanı kullanabiliriz -R ve -Jm0.1 almak
eşleştirmek için iki sistem.

Ancak, bundan daha akıllı olabiliriz. Girdi img dosyasının piksel 2.0 olduğunu fark etmek
dakika genişliğinde (veya grdinfo merc_grav.nc ile nx ve ny'yi kontrol ederek) fark ederiz ki
merc_grav.nc, img dosyasının tam çözünürlüğünü kullandı ve 2400 x 2039 piksele sahip ve
8 inç genişliğinde bu, inç başına 300 pikseldir. O kadar çok şeye ihtiyacımız olmadığına karar veriyoruz ve
inç başına 100 piksel ile tatmin olacağından, verilerin ortalamasını 3'e 3 olarak almak istiyoruz
kareler. (Bir kontur grafiği istiyorsak, muhtemelen verilerin ortalamasını çok daha fazla seçeceğiz.
(örneğin, 6'ya 6) düzgün konturlar elde etmek için.) 2039 3'e tam bölünemediğinden
farklı ayarlanmış -R bu zaman:

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

Bu sefer ayarlanan bölgeyi buluyoruz: -R-40/40/-70.023256525/-29.9368261101 ve
çıktı 800 x 601 piksel, bizim için daha iyi bir boyut. Şimdi yapay bir tane oluşturabiliriz.
Bunun için aydınlatma dosyası kullanarak grdgradyan:

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

ve "grav.cpt" adında bir CPT dosyamız varsa, renk gölgeli bir kabartma haritası oluşturabiliriz
bunun gibi:

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

Enlem/boylam olarak bir img dosyasından yalnızca kısıtlı veri değerlerini almak istediğinizi varsayalım.
koordinatlar. O zaman koş img2grd ile -T2 seçenek, kullanın grd2xyz değerleri boşaltmak için, boru
NaN'leri ortadan kaldırmak için grep -v NaN aracılığıyla ve boru yoluyla harita projesi tersi ile
projeksiyon yukarıdaki gibidir.

onworks.net hizmetlerini kullanarak img2grdgmt'yi çevrimiçi kullanın