img2grdgmt - Online w chmurze

PROGRAM:

IMIĘ

img2grd - Wyodrębnij podzbiór pliku img w formacie Mercator lub Geographic

STRESZCZENIE

img2grd plik obrazu plik grd region rodzaj [ ] [ [minlat/makslat] ] [ ] [ minuty ] [ ] [ nawigacja
] [ [skala] ] [ [poziom] ] [ maxlon ] [ -n]

Uwaga: Między flagą opcji a powiązanymi argumentami nie może być spacji.

OPIS

img2grd odczytuje plik w formacie img, wyodrębnia podzbiór i zapisuje go w pliku siatki. ten -M
opcja określa, czy projekcja Spherical Mercator pliku img jest
zachowane lub jeśli siatka geograficzna powinna zostać zapisana przez cofnięcie odwzorowania Mercatora. Gdyby
wybrano siatkę geograficzną, można również poprosić o ponowne próbkowanie na dokładny -R dany.

WYMAGANE ARGUMENTY

plik obrazu
Plik w formacie Mercator img, taki jak pola grawitacji morskiej lub topografii dna morskiego
oszacowane na podstawie danych z wysokościomierza satelitarnego Sandwella i Smitha. Jeśli użytkownik ustawił
zmienna środowiskowa $GMT_DATADIR, następnie img2grd spróbuje znaleźć plik obrazu in
$GMT_DATADIR; inaczej spróbuje się otworzyć plik obrazu bezpośrednio.

-Gplik grd
plik grd to nazwa wyjściowego pliku siatki.

-R[jednostka]zachód/wschód/południe/na północ[/z min/zmaks][R]
zachód, wschód, południe, na północ określ region zainteresowania i możesz określić
je w stopniach dziesiętnych lub w formacie [+-]dd:mm[:ss.xxx][W|E|S|N]. Dodać r jeśli niższy
zamiast w/e/s/n podano współrzędne lewego i prawego górnego rogu mapy. Dwójka
skróty -Rg i -R & D oznacza domenę globalną (0/360 i -180/+180 w długości geograficznej)
odpowiednio z -90/+90 na szerokości geograficznej). Alternatywnie do tworzenia siatki, daj
Rkodon/lat/nx/ny, Gdzie kod to dwuznakowa kombinacja L, C, R (dla lewej,
środek lub prawo) i T, M, B dla góry, środka lub dołu. np. BL dla dolnego lewego.
Wskazuje, który punkt na prostokątnym obszarze lon/lat odwołuje się współrzędne
do i wymiary siatki nx i ny z odstępami siatki przez -I służy do tworzenia
odpowiedni region. Alternatywnie podaj nazwę istniejącego pliku siatki
oraz -R ustawienia (i odstępy siatki, jeśli dotyczy) są kopiowane z siatki.
Korzystanie z -Rjednostka oczekuje rzutowanych (kartezjańskich) współrzędnych zgodnych z wybranym -J
i odwrotnie prognozujemy, aby określić rzeczywisty prostokątny region geograficzny. Do
widok perspektywiczny (-p), opcjonalnie dołącz /z min/zmaks. W przypadku widoku perspektywicznego
(-p), zakres z (z min, zmaks) można dołączyć w celu wskazania trzeciego wymiaru. Ten
należy to zrobić tylko podczas korzystania z -J z opcja, nie przy korzystaniu tylko z -p opcja.
W tym drugim przypadku wykreślany jest widok perspektywiczny samolotu, bez trzeciego
wymiar.

OPCJA ARGUMENTY

-C Ustaw współrzędne Mercatora x i y względem środka projekcji [Domyślnie
względem dolnego lewego rogu siatki]. Wymaga -M.

-D[minlat/makslat]
Użyj rozszerzonego zakresu szerokości geograficznej -80.738/+80.738. Alternatywnie dołącz
minlat/makslat jako zakres szerokości geograficznej wejściowego pliku img. [Domyślnie:
-72.006/72.006]. Zwykle nie jest to wymagane, ponieważ możemy określić zakres z
sprawdzenie rozmiaru pliku.

-E Może być używany, gdy -M nie jest ustawiony tak, aby wymuszał dokładnie taką samą siatkę końcową
region zgodnie z żądaniem z -R. Domyślnie końcowy region jest bezpośrednim odwzorowaniem
oryginalny region Mercator i zazwyczaj będzie rozciągał się nieco poza obszar
żądany zakres szerokości geograficznej, a ponadto przyrost siatki w szerokości geograficznej nie
dopasować przyrost długości geograficznej. Jednak dodatkowe resampling wprowadza małe
błędy interpolacji i powinny być używane tylko wtedy, gdy siatka wyjściowa musi być zgodna z
żądany region i mają x_inc = y_inc. W tym przypadku region ustawiony przez -R musi być
podane w wielokrotnościach przyrostu (np. -R0/45/45/72).

-I Wskazać minuty jako szerokość wejściowego piksela img w minutach długości geograficznej.
[Domyślnie 2.0]. Zwykle nie jest to wymagane, ponieważ możemy określić rozmiar piksela z
kontrola rozmiaru.

-M Wygeneruj siatkę Spherical Mercator [Domyślnie jest to geograficzna siatka lon/lat]. ten
Projekcja sferyczna Mercatora pliku img jest zachowana, dzięki czemu region -R
ustawiony przez użytkownika jest nieznacznie modyfikowany; zmodyfikowany region odpowiada krawędziom
pikseli [lub grup nawigacja piksele]. Nagłówek pliku siatki jest ustawiony tak, że x i
długości osi y reprezentują odległość od zachodniej i południowej krawędzi obrazu,
mierzone w domyślnych jednostkach użytkownika, z -Jm1 i skorygowany -R. Poprzez ustawienie
domyślnym PROJ_ ELIPSOIDA = Kula, użytkownik może tworzyć nakładki z dostosowanym -R
aby pasowały. Widzieć PRZYKŁADY poniżej. Dostosowany -R jest również napisane w siatce
uwaga nagłówka, aby można ją było znaleźć później. Widzieć -C ustawić współrzędne względem
centrum projekcji.

-Nnawigacja Uśrednij wartości w wejściowych pikselach img do nawigacja by nawigacja kwadraty i tworzenie
jeden piksel wyjściowy na każdy taki kwadrat. Jeśli jest używany z -T3 zgłosi średnią
ograniczenie od 0 do 1. Jeśli jest używane z -T2 wyjście będzie średnią wartością danych
lub NaN w zależności od tego, czy średnie ograniczenie jest > 0.5. nawigacja musi równomiernie dzielić
do wymiarów pliku imgfile w pikselach. [Domyślny 1 nie ma uśredniania].

-S[skala]
Pomnóż wartości pliku img przez skala przed zapisaniem w pliku siatki. [Domyślnie:
1.0]. Dla ostatnich plików img: pliki img topo są przechowywane w (poprawionych) miernikach [-S1];
pilniki grawitacyjne na wolnym powietrzu w mGal*10 [-S0.1, aby uzyskać mGal]; pilniki odchylania pionowego w
mikroradiany*10 [-S0.1, aby uzyskać mikro-radiany], pliki gradientu grawitacji pionowej w
Eotvos*50 [-S0.02, aby zdobyć Eotvos, lub -S0.002 do uzyskania mGal/km]). Jeśli nie skala jest podawany
próbujemy określić skalę, badając nazwę pliku w poszukiwaniu wskazówek.

-Trodzaj rodzaj obsługuje kodowanie informacji o ograniczeniach. rodzaj = 0 oznacza, że ​​nie
takie informacje są zakodowane w pliku img (używanym w wersjach sprzed 1995 r.
dane grawitacyjne) i pobiera wszystkie dane. rodzaj > 0 wskazuje, że informacja o ograniczeniu jest
zakodowane (1995 i późniejsze (aktualne) wersje plików img), dzięki czemu można:
utwórz plik siatki w następujący sposób: -T1 pobiera wartości danych we wszystkich punktach, -T2 pobiera dane
wartości w punktach ograniczonych i NaN w punktach interpolowanych; -T3 dostaje 1 w
punkty ograniczone i 0 w punktach interpolowanych [Domyślnie 1].

-V[poziom] (więcej ...)
Wybierz poziom szczegółowości [c]. Szczególnie polecane tutaj, ponieważ warto zobaczyć
jak dopasowywane są współrzędne.

-Wmaxlon
Wskazać maxlon jako maksymalny zakres długości wejściowego pliku img. Wersje
od 1995 roku maxlon = 360.0, podczas gdy niektóre wcześniejsze pliki miały maxlon = 390.0.
[Domyślnie 360.0].

-^ or właśnie -
Wydrukuj krótką wiadomość o składni polecenia, a następnie wyjdź (UWAGA: w systemie Windows
użyj tylko -).

-+ or właśnie +
Wydrukuj obszerny komunikat dotyczący użytkowania (pomocy), w tym wyjaśnienie wszelkich
opcja specyficzna dla modułu (ale nie typowe opcje GMT), a następnie kończy działanie.

-? or Nie argumenty
Wydrukuj pełny komunikat (pomoc) dotyczący użytkowania, w tym wyjaśnienie opcji, a następnie
wyjść.

--wersja
Wydrukuj wersję GMT i wyjdź.

--show-katalog danych
Wydrukuj pełną ścieżkę do katalogu współdzielonego GMT i wyjdź.

GEOGRAFICZNY PRZYKŁADY

Połączenia -M opcja powinna być wyłączona, jeśli chcesz, aby siatka wyjściowa była geograficzna
współrzędne. Aby wyodrębnić dane w regionie -R-40/40/-70/-30 od świat_grav.img.7.2 i
przeprojektuj, aby uzyskać współrzędne geograficzne, możesz spróbować

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

Ponieważ odstępy szerokości geograficznej w pliku img są równe odległości w jednostkach Mercator,
wynikowa siatka nie będzie pasować do podanej -R dokładnie, a odstępy szerokości geograficznej nie będą
równe odstępie długości geograficznej. Jeśli potrzebujesz dokładnego dopasowania do swojego -R i te same odstępy
w długości i szerokości geograficznej użyj -E opcja:

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

MERKATOR PRZYKŁADY

Ponieważ pliki img są w rzucie Mercator, NIE powinieneś wyodrębniać siatki geograficznej
jeśli twoim planem jest stworzenie mapy Mercator. Jeśli to zrobiłeś, kończysz projekcję i
reprojekcji siatki, tracąc szczegóły na krótkich falach. Lepiej używać -M i narysuj siatkę
przy użyciu projekcji liniowej o tej samej skali, co żądana projekcja Mercator (patrz GMT
Przykład 29). Aby wyodrębnić dane w regionie -R-40/40/-70/-30 od świat_grav.img.7.2, biegać

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

Zauważ, że -V opcja informuje nas, że zakres został dostosowany do
-R-40/40/-70.0004681551/-29.9945810754. Możemy również użyć informacje o grd aby znaleźć plik siatki
nagłówek pokazuje, jaki region ma być -R0/80/0/67.9666667 To jest przedział x,y który otrzymamy
z projekcji Spherical Mercator przy użyciu -R-40/40/-70.0004681551/-29.9945810754 i
-Jm1. Tak więc, aby pobrać plik ship.lonlatgrav i użyć go do spróbkowania pliku merc_grav.nc, możemy to zrobić:

gmt ustaw sferę PROJ_ELLIPSOID

gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i Ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | projekt mapy gmt \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > statek.lonlatgravsat

Zaleca się zastosowanie powyższej metody projekcji i odwzorowywania danych w takich
aplikacji, ponieważ wtedy jest tylko jeden krok interpolacji (w grdtrack). Gdyby jeden
najpierw próbuje przekonwertować plik siatki na lon,lat, a następnie próbkować, są dwa
kroki interpolacji (podczas konwersji i próbkowania).

Aby zrobić siatkę dł, łat z powyższej siatki, możemy użyć

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

W niektórych przypadkach nie będzie to łatwe, ponieważ -R w dwóch układach współrzędnych może się nie zgadzać
dobrze. Kiedy tak się stanie, możemy również użyć (w rzeczywistości zawsze może być lepiej użyć)

gmt grd2xyz merc_grav.nc | projekt mapy gmt \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
powierzchnia gmt -R-40/40/-70/70 -I2m -Ggrav.nc

Aby utworzyć mapę Mercator powyższego regionu, załóżmy, że nasza wartość gmt.conf dla
PROJ_LENGTH_UNIT to cal. Następnie, ponieważ powyższy plik merc_grav.nc jest wyświetlany z -Jm1i
ma 80 cali szerokości. Możemy zrobić mapę o szerokości 8 cali za pomocą -Jx0.1i w dowolnych programach map
zastosowane do tej siatki (np. grdkontur, grdobraz, widok grd), a następnie dla nakładek, które
praca w dł,lat (np. psx, pscoast) możemy użyć powyższego dostosowanego -R i -Jm0.1, aby uzyskać
dwa systemy do dopasowania.

Jednak możemy być mądrzejsi. Zdając sobie sprawę, że wejściowy plik img miał piksele 2.0
minut (lub sprawdzając nx i ny za pomocą grdinfo merc_grav.nc) zdajemy sobie sprawę, że
merc_grav.nc wykorzystał pełną rozdzielczość pliku img i ma 2400 na 2039 pikseli, oraz
przy szerokości 8 cali jest to 300 pikseli na cal. Postanawiamy, że nie potrzebujemy tak wielu i my
zadowoli się 100 pikselami na cal, więc chcemy uśrednić dane na 3 na 3
kwadraty. (Jeśli chcemy uzyskać wykres konturowy, prawdopodobnie zdecydujemy się na znacznie większe uśrednienie danych
(np. 6 na 6) aby uzyskać gładkie kontury. Ponieważ rok 2039 nie jest podzielny przez 3, otrzymamy
różne dostosowane -R tym razem:

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

Tym razem okazuje się, że dostosowany region to -R-40/40/-70.023256525/-29.9368261101 i
wyjście to 800 na 601 pikseli, dla nas lepszy rozmiar. Teraz możemy stworzyć sztuczną
plik oświetlenia do tego przy użyciu gradient:

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

a jeśli mamy również plik CPT o nazwie „grav.cpt”, możemy utworzyć mapę reliefową z cieniowaniem w kolorze
lubię to:

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

Załóżmy, że chcesz uzyskać tylko ograniczone wartości danych z pliku img, w formacie lat/lon
współrzędne. Następnie uruchomić img2grd z -T2 opcja, użyj grd2xyz zrzucić wartości, rura
przez grep -v NaN, aby wyeliminować NaN i przepuścić przez potok mapaprojekt z odwrotnością
projekcja jak wyżej.

Korzystaj z img2grdgmt online za pomocą usług onworks.net