grdgradientgmt - Online w chmurze

PROGRAM:

IMIĘ

grdgradient – ​​Oblicza pochodną kierunkową lub gradient z siatki

STRESZCZENIE

gradient plik_grd plik_grd [ azim[/azym2] ] [ [a][c][o][n] ] [
[s|p]azim/elev[/otoczenia/rozpraszać/wziernikowy/błyszczeć] ] [ flaga ] [ [e][t][amp][/sigma[/offset]]
] [ region ] [ plik nachylenia ] [ [poziom] ] [ -fg ] [ -n]

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

OPIS

gradient można użyć do obliczenia pochodnej kierunkowej w danym kierunku (-A),
lub kierunek (-S) [i wielkość (-D)] gradientu wektorowego danych.

Szacowane wartości w pierwszym/ostatnim wierszu/kolumnie wyników zależą od warunków brzegowych (patrz
-L).

WYMAGANE ARGUMENTY

plik_grd
Plik siatki 2D, z którego można obliczyć pochodną kierunkową. (Zobacz FORMATY PLIKÓW SIATKI)
poniżej).

-Gplik_grd
Nazwa pliku siatki wyjściowej dla pochodnej kierunkowej. (Zobacz FORMATY PLIKÓW SIATKI)
poniżej).

OPCJA ARGUMENTY

-Aazim[/azym2]
Kierunek azymutalny dla pochodnej kierunkowej; azim jest kątem w x,y
płaszczyzna mierzona w stopniach w kierunku dodatnim zgodnie z ruchem wskazówek zegara od północy (kierunek +y) w kierunku
wschód (kierunek +x). Negatywna wartość pochodnej kierunkowej,
-[dz/dx*sin(azim) + dz/dy*cos(azim)], jest znaleziony; negacja daje wartości dodatnie
gdy nachylenie z(x,y) jest spadkowe azim kierunek, właściwy sens
zacienianie oświetlenia obrazu (patrz grdobraz oraz widok grd) przez źródło światła
nad płaszczyzną x,y świecącą z azim kierunek. Opcjonalnie, podaj dwa
azymuty, -Aazim/azym2, w którym przypadku gradienty w każdym z tych kierunków są
obliczana i zachowywana jest ta większa; jest to przydatne w przypadku
oświetlanie danych dwoma kierunkami struktur liniowych, np. -A0/270
oświetla od północy (góra) i zachodu (lewa strona).

-D[a][c][o][n]
Znajdź kierunek dodatniego (wznoszącego się) gradientu danych. Zamiast tego
znajdź aspekt (kierunek spadku), użyj -TakDomyślnie kierunki są następujące:
mierzona zgodnie z ruchem wskazówek zegara od północy, jako azim in -A powyżej. Dodaj c używać konwencjonalnego
Kąty kartezjańskie mierzone przeciwnie do ruchu wskazówek zegara od dodatniego kierunku x (wschód).
Dodać o aby zgłosić orientacje (0-180) zamiast kierunków (0-360). Dołącz n
dodać 90 stopni do wszystkich kątów (np. w celu uzyskania lokalnych uderzeń powierzchni).

-E[s|p]azim/elev[/otoczenia/rozpraszać/wziernikowy/błyszczeć]
Oblicz radiancję Lamberta odpowiednią do użycia z grdobraz oraz widok grd,
Odbicie lambertowskie zakłada idealną powierzchnię odbijającą całe światło, które
Uderza w nią, a powierzchnia wydaje się jednakowo jasna z każdego kierunku. azim
oraz podwyższenie są azymutem i elewacją wektora światła. Opcjonalnie, podaj otoczenia
rozpraszać wziernikowy błyszczeć które są parametrami kontrolującymi właściwości odbicia
powierzchni. Wartości domyślne to: 0.55/0.6/0.4/10 Aby pozostawić niektóre wartości
nietknięty, określ = jako nową wartość. Na przykład -E60/30/=/0.5 ustawia azim podwyższenie
oraz rozpraszać do 60, 30 i 0.5, pozostawiając pozostałe parametry odbicia
nietknięty. Dodaj s aby użyć prostszego algorytmu Lamberta. Należy zauważyć, że w tej formie
musisz podać tylko parametry azymutu i elewacji. Dołącz p używać
Aproksymacja liniowa Peuckera kawałkami (prostszy, ale szybszy algorytm; w tym przypadku
dotychczasowy azim oraz podwyższenie są podłączone do 315 i 45 stopni. Oznacza to, że nawet jeśli
(podaj inne wartości – zostaną zignorowane.)

-Lflaga Warunek brzegowy flaga może być x or y or xy wskazujące, że dane są okresowe w zakresie
x lub y lub oba, lub flaga może być g wskazujące warunki geograficzne (x i y to
lon i lat). [Domyślnie używane są warunki „naturalne” (druga pochodna cząstkowa normalnej
do krawędzi wynosi zero).]

-Internet][amp][/sigma[/offset]]
Normalizacja. [Domyślnie: brak normalizacji.] Rzeczywiste gradienty g są przesunięte i
skalowane w celu wytworzenia znormalizowanych gradientów gn o maksymalnej wielkości wyjściowej amp.
If amp nie podano, domyślnie amp = 1. Jeśli offset nie jest podany, jest ustawiony na
Średnia z g. -N plony gn = amp * (g - offset)/maks.(bezwzględne(g - offset)). -Nie
normalizuje się za pomocą kumulatywnego rozkładu Laplace'a, co daje gn = amp * (1.0 -
exp (kwadrat(2) * (g - offset)/ sigma)) gdzie sigma szacuje się przy użyciu normy L1
(g - offset) jeśli nie jest podany. -Nt normalizuje się za pomocą kumulatywnego równania Cauchy’ego
dystrybucja dająca plony gn = (2 * amp /PI) * atan( (g - offset)/ sigma) gdzie sigma
szacuje się przy użyciu normy L2 (g - offset) jeśli nie jest podany.

-R[jednostka]x min/xmaks/ymin/ymaks[R] (więcej ...)
Określ interesujący Cię region. Używając -R opcja wybierze podsekcję
plik_grd siatka. Jeśli ta podsekcja wykracza poza granice siatki, tylko
wspólny region zostanie wyodrębniony.

-Splik nachylenia
Nazwa pliku siatki wyjściowej z wielkościami skalarnymi wektorów gradientu. Wymaga -D
ale sprawia -G opcjonalny.

-V[poziom] (więcej ...)
Wybierz poziom szczegółowości [c].

-fg Siatki geograficzne (wymiary długości i szerokości geograficznej) zostaną przeliczone na metry
poprzez przybliżenie „Płaskiej Ziemi” przy użyciu bieżących parametrów elipsoidy.

-n[b|c|l|n][+a][+bBC][+c][+tpróg] (więcej ...)
Wybierz tryb interpolacji dla siatek.

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

GRID DYSTANS JEDNOSTKI

Jeśli siatka nie ma licznika jako jednostki poziomej, dołącz +ujednostka do pliku wejściowego
nazwa do konwersji z określonej jednostki na metr. Jeśli twoja siatka jest geograficzna, przekonwertuj
odległości na metry przez dostarczanie -fg zamiast.

PORADY

Jeśli nie wiesz co -N opcje do wykorzystania przy tworzeniu pliku intensywności grdobraz or
widok grd, dobra pierwsza próba to -Nie0.6.

Zwykle 255 odcieni jest więcej niż wystarczające do celów wizualizacyjnych. Możesz zaoszczędzić 75% dysku
spacja poprzez dodanie =nb/a do nazwy pliku wyjściowego plik_grd.

Jeśli chcesz wykonać kilka podświetlanych map podregionów dużego zestawu danych i
jeśli chcesz, aby efekty oświetlenia były spójne na wszystkich mapach, użyj -N opcja i
dostarczyć taką samą wartość sigma oraz offset do gradient dla każdej mapy. Dobrym przypuszczeniem jest
offset = 0 i sigma znalezione przez informacje o grd -L2 or -L1 zastosowano do nienormalizowanego gradientu grd.

Jeśli po prostu potrzebujesz x- lub y-pochodne siatki, użyj grdmat.

GRID FILE FORMATY

Domyślnie GMT zapisuje siatkę jako zmiennoprzecinkową o pojedynczej precyzji w netCDF z reklamacją COARDS
format pliku. Jednak GMT jest w stanie tworzyć pliki grid w wielu innych powszechnie używanych siatkach
formatów plików, a także ułatwia tzw. „pakowanie” siatek, wypisywanie zmiennoprzecinkowych
dane jako 1- lub 2-bajtowe liczby całkowite. Aby określić precyzję, skalę i przesunięcie, użytkownik powinien:
dodaj przyrostek =id[/skala/offset[/nan]], gdzie id jest dwuliterowym identyfikatorem siatki
rodzaj i precyzja oraz skala oraz offset są opcjonalnym współczynnikiem skali i przesunięciem, które mają być
zastosowane do wszystkich wartości siatki, oraz nan to wartość używana do wskazania brakujących danych. W razie
dwie postacie id nie jest dostarczony, jak w =/skala niż a id=nf zakłada się. Kiedy
czytanie siatek, format jest zazwyczaj automatycznie rozpoznawany. Jeśli nie, ten sam przyrostek
można dodać do nazw plików siatki wejściowej. Widzieć grdconvert i format pliku siatki sekcji
Informacje techniczne GMT i książka kucharska, aby uzyskać więcej informacji.

Podczas odczytywania pliku netCDF, który zawiera wiele siatek, GMT domyślnie odczyta
pierwsza dwuwymiarowa siatka, którą można znaleźć w tym pliku. Aby nakłonić GMT do przeczytania innego
wielowymiarowa zmienna w pliku siatki, dołącz ?nazwa_zmiennej do nazwy pliku, gdzie
nazwa_zmiennej to nazwa zmiennej. Pamiętaj, że być może będziesz musiał uciec od specjalnego znaczenia
of ? w swoim programie powłoki, umieszczając przed nim ukośnik odwrotny lub umieszczając znak
nazwę pliku i przyrostek między cudzysłowami lub cudzysłowami. ten ?nazwa_zmiennej można również użyć sufiksu
dla siatek wyjściowych, aby określić nazwę zmiennej inną niż domyślna: „z”. Widzieć
grdconvert i sekcje modyfikatory-for-CF i format pliku siatki GMT Technical
Odnośnik i książka kucharska, aby uzyskać więcej informacji, w szczególności o tym, jak czytać sploty 3-,
Siatki 4- lub 5-wymiarowe.

PRZYKŁADY

Aby utworzyć plik do oświetlenia danych w geoid.nc przy użyciu znormalizowanych gradientów exp
zakres [-0.6,0.6] imitujący źródła światła w kierunku północnym i zachodnim:

gmt grdgradient geoid.nc -A0/270 -Ggradients.nc=nb/a -Ne0.6 -V

Aby znaleźć orientacje azymutalne tkaniny dna morskiego w pliku topo.nc:

gmt grdgradient topo.nc -Dno -Gazimuths.nc -V

LITERATURA

Horn, BKP, Cieniowanie wzniesień i mapa odbicia światła, Proceedings of the IEEE, tom 69, nr XNUMX (XNUMX).
1, styczeń 1981, s. 14-47. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

Użyj grdgradientgmt online za pomocą usług onworks.net