mapprojectgmt — online w chmurze

PROGRAM:

IMIĘ

mapproject - Wykonuj przekształcenia map w przód i w tył, konwersje danych i geodezji

STRESZCZENIE

mapaprojekt [ stoły ] parametry region [ b|B|f|F|o|O[lon0/Lat0] ] [ [dx/dy] ] [ c|i|p ] [
[podstawa odniesienia] ] [ [jednostka] ] [ [x0/y0/][[+|-]jednostka][+|-] ] [ ] [ linia.xy[/[+|-]jednostka][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]od[/do] ] [ [poziom] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o] ]

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

OPIS

mapaprojekt odczytuje (długość, szerokość geograficzną) pozycje z stoły [lub standardowe wejście] i
oblicza współrzędne (x,y) przy użyciu określonego odwzorowania mapy i skal. Opcjonalnie to
potrafi odczytywać (x,y) pozycje i obliczać wartości (długość, szerokość geograficzną) wykonując odwrotność
transformacja. Można to wykorzystać do przekształcenia liniowych (x,y) punktów uzyskanych przez digitalizację
mapa znanej projekcji na współrzędne geograficzne. Może również obliczyć odległości wzdłuż
do ustalonego punktu lub najbliższego podejścia do linii. Wreszcie może być używany do wykonywania
różne konwersje danych. Dodatkowe pola danych są dozwolone po pierwszych 2 kolumnach
które muszą mieć (długość,szerokość geograficzna) lub (x,y). Zobacz opcję -: jak czytać
(szerokość, długość geograficzna).

WYMAGANE ARGUMENTY

-Jparametry (więcej ...)
Wybierz odwzorowanie mapy.

-R[jednostka]x min/xmaks/ymin/ymaks[R] (więcej ...)
Określ interesujący Cię region. Specjalny przypadek dla projekcji UTM: If -C służy
i -R nie jest podany, region jest ustawiony tak, aby pokrywał się z daną strefą UTM, więc
w celu zachowania pełnego elipsoidalnego rozwiązania (więcej informacji znajdziesz w OGRANICZENIACH
Informacja).

OPCJA ARGUMENTY

stół Jeden lub więcej ASCII (lub binarny, patrz -bi[podoficerowie][rodzaj]) plik(i) tabeli danych zawierający a
liczba kolumn danych. Jeśli nie podano żadnych tabel, czytamy ze standardowego wejścia.

-Ab|B|f|F|o|O[lon0/Lat0]
-Af oblicza (do przodu) azymut z punktu stałego dł./łac do każdego punktu danych.
Zastosowanie -Z dala aby uzyskać wsteczny azymut z punktów danych do punktu stałego. Posługiwać się -Ao dostać
orientacje (-90/90) zamiast azymutów (0/360). Duże litery F, B or O będzie
konwersja z szerokości geodezyjnych na geocentryczne i oszacowanie azymutu geodezyjnego
(zakładając, że aktualna elipsoida nie jest sferą). Jeśli nie podano punktu stałego, to
obliczamy azymut (lub azymut wsteczny) od poprzedniego punktu.

-C[dx/dy]
Ustaw środek rzutowanych współrzędnych tak, aby znajdował się w centrum odwzorowania mapy [Domyślnie:
lewy dolny róg]. Opcjonalnie dodaj odsunięcia w jednostkach rzutowania, które mają zostać dodane (lub
odejmowane, gdy -I jest ustawiona) na (z) rzutowanych współrzędnych, np. false
wschodnie i północne dla poszczególnych stref rzutu [0/0]. Jednostka używana do
offsets to obowiązująca jednostka odległości wydruku (patrz PROJ_LENGTH_UNIT), chyba że -F is
używane, w którym to przypadku przesunięcia są w metrach.

-Dc|i|p
Tymczasowo zastąp PROJ_LENGTH_UNIT i użyj c (cm), i (cale) lub p (punkty)
zamiast. Nie można używać z -F.

-MI[podstawa odniesienia]
Konwersja z geodezyjnego (lon, lat, height) na Earth Centered Earth Fixed (ECEF)
(x,y,z) współrzędne (dodaj -I dla konwersji odwrotnej). Dołącz identyfikator odniesienia (patrz -Qd)
lub dać elipsoida:dx,dy,dz gdzie elipsoida może być identyfikatorem elipsoidy (patrz -Q) lub
podane jako a[,*inv_f*], gdzie a jest półoś wielką i inv_f jest odwrotnością?
spłaszczenie (0, jeśli pominięto). Gdyby podstawa odniesienia jest - lub nie podano zakładamy WGS-84.

-F[jednostka]
Wymuś skalowanie 1:1, tj. wyjście (lub wejście, patrz -I) dane są w rzeczywistych rzutach
metrów. Aby określić inne jednostki, dodaj żądaną jednostkę (patrz JEDNOSTKI). Z pominięciem -F,
wyjście (lub wejście, patrz -I) są w jednostkach określonych przez PROJ_LENGTH_UNIT (ale
widzieć -D).

-G[x0/y0/][[+|-]jednostka][+|-]
Oblicz odległości wzdłuż toru or do opcjonalnego punktu ustawionego za pomocą -Gx0/rok0. Dodać
jednostka odległości (patrz JEDNOSTKI), w tym c (Odległość kartezjańska za pomocą wprowadzania
współrzędne) lub C (Odległość kartezjańska przy użyciu współrzędnych rzutowanych). ten C jednostka
Wymaga -R i -J do ustawienia. Bez podania stałego punktu obliczamy kumuluj
odległości wzdłuż toru. Dodać - aby uzyskać przyrostową odległość między kolejnymi
zwrotnica. Dodać + aby określić drugi punkt za pomocą dwóch dodatkowych kolumn w pliku wejściowym.

-I Wykonaj transformację odwrotną, tj. pobierz (długość, szerokość geograficzną) z danych (x,y).

-Llinia.xy[/[+|-]jednostka][+]
Określ najkrótszą odległość od punktów danych wejściowych do linii podanych w
plik wielosegmentowy ASCII linia.xy. Odległość i współrzędne
najbliższy punkt zostanie dodany do wyniku jako trzy nowe kolumny. Dołącz
jednostka odległości (patrz JEDNOSTKI), w tym c (Odległość kartezjańska przy użyciu współrzędnych wejściowych)
or C (Odległość kartezjańska przy użyciu współrzędnych rzutowanych). ten C jednostka wymaga -R i
-J do ustawienia. Na koniec dołącz + aby zgłosić identyfikator segmentu linii i ułamek
numer punktu zamiast lon/lat najbliższego punktu.

-N[a|c|g|m]
Konwersja z szerokości geodezyjnych (przy użyciu bieżącej elipsoidy; patrz PROJ_ELLIPSOID)
na jedną z czterech różnych szerokości geograficznych pomocniczych (długości geograficzne pozostają bez zmian). Wybierać
od autalny, cnieformalny, geocentryczny i mszerokości geograficzne [geocentryczne]. Posługiwać się -I
do konwersji z pomocniczych szerokości geograficznych na szerokości geodezyjne.

-Q[d|e Wymień wszystkie parametry projekcji. Aby wyświetlić tylko układy odniesienia, użyj -Qd. Aby tylko wymienić
elipsoidy, użyj -Q.

-S Pomijaj punkty, które wychodzą poza region.

-NS]od[/do]
Konwersje współrzędnych między układami odniesienia od i do używając standardowego Molodensky
transformacja. Posługiwać się -NS jeśli trzecia kolumna wejściowa ma wysokość powyżej elipsoidy [Domyślnie
przyjmuje wysokość = 0, czyli na elipsoidzie]. Określ punkty odniesienia używając ID punktu odniesienia (zobacz
-Qd) lub daj elipsoida:dx,dy,dz gdzie elipsoida może być identyfikatorem elipsoidy (patrz -Q) lub
podane jako a[,*inv_f*], gdzie a jest półoś wielką i inv_f jest odwrotnością?
spłaszczenie (0, jeśli pominięto). Gdyby podstawa odniesienia jest - lub nie podano zakładamy WGS-84. -T może być
używany w połączeniu z -R -J zmienić punkt odniesienia przed rzutowaniem współrzędnych
(Dodaj -I aby zastosować konwersję bazy pomiarowej po odwzorowaniu odwrotnym). Upewnij się, że
ustawienie PROJ_ELLIPSOID jest prawidłowe w Twoim przypadku.

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

-W[sz|h]
Drukuje szerokość i wysokość mapy na standardowym wyjściu. Żadne pliki wejściowe nie są odczytywane. Tylko
wypisz szerokość lub wysokość, dołącz w or h, odpowiednio. Jednostki
wymiary można zmienić za pomocą -D.

-bi[podoficerowie][T] (więcej ...)
Wybierz natywne wejście binarne. [Domyślnie 2 kolumny wejściowe].

-bo[podoficerowie][rodzaj] (więcej ...)
Wybierz natywne wyjście binarne. [Domyślnie jest takie samo jak wejście].

-d[i|o]brak danych (więcej ...)
Zastąp kolumny wejściowe, które są równe brak danych z NaN i zrób odwrotną stronę na wyjściu.

-f[i|o]kolinfo (więcej ...)
Określ typy danych kolumn wejściowych i/lub wyjściowych.

-g[a]x|y|d|X|Y|D|[col]z[+|-]szczelina[u] (więcej ...)
Określ luki w danych i podziały wierszy.

-h[i|o][n][+c][+d][+ruwagę][+rtytuł] (więcej ...)
Pomiń lub utwórz rekord(y) nagłówka.

-icols[l] [sskala][ooffset][,...] (więcej ...)
Wybierz kolumny wejściowe (0 to pierwsza kolumna).

-ocols[,...] (więcej ...)
Wybierz kolumny wyjściowe (0 to pierwsza kolumna).

-p[x|y|z]azim/podwyższenie[/zpoziom][+wlon0/Lat0[/z0]][+wx0/y0] (więcej ...)
Wybierz widok perspektywiczny.

-S[cols[a|r] (więcej ...)
Ustaw obsługę rekordów NaN.

-:[i|o] (więcej ...)
Zamień pierwszą i drugą kolumnę na wejściu i/lub wyjściu.

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

JEDNOSTKI

W przypadku jednostki odległości na mapie dołącz jednostka d dla stopnia łuku, m dla minuty kątowej, i s dla łuku
drugi, lub e na licznik [Domyślnie], f na stopę, k za km, M za milę ustawową, n dla żeglarstwa
mila i u dla amerykańskiej stopy pomiarowej. Domyślnie obliczamy takie odległości za pomocą sferycznego
zbliżenie z wielkimi kołami. Przedstaw - na odległość (lub jednostka nie ma odległości
podane), aby wykonać obliczenia „Płaskiej Ziemi” (szybsze, ale mniej dokładne) lub dołączyć + do
wykonać dokładne obliczenia geodezyjne (wolniej, ale dokładniej).

ASCII FORMAT PRECISION

Formaty wyjściowe ASCII danych liczbowych są kontrolowane przez parametry w twoim konf.gmt
plik. Długość i szerokość geograficzna są formatowane zgodnie z FORMAT_GEO_OUT, podczas gdy inne
wartości są sformatowane zgodnie z FORMAT_FLOAT_OUT. Należy pamiętać, że obowiązujący format może:
prowadzić do utraty precyzji na wyjściu, co może prowadzić do różnych problemów na dalszych etapach. Gdyby
okaże się, że dane wyjściowe nie są napisane z wystarczającą precyzją, rozważ przełączenie na tryb binarny
wyjście (-bo jeśli to możliwe) lub określ więcej miejsc dziesiętnych za pomocą ustawienia FORMAT_FLOAT_OUT.

PRZYKŁADY

Aby przekonwertować współrzędne UTM w metrach na lokalizacje geograficzne, biorąc pod uwagę plik utm.txt i
znając strefę UTM (i strefę lub półkulę), spróbuj

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

Aby przekształcić plik z (długość, szerokość geograficzna) na pozycje (x,y) w cm na Mercator
siatka dla danej skali 0.5 cm na stopień, bieg

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

Aby przekształcić kilka dwukolumnowych plików binarnych o podwójnej precyzji za pomocą (szerokość, długość geograficzna)
do (x,y) pozycji w calach na siatce poprzecznego Mercatora (centralna długość 75W) dla
skala = 1:500000 i stłumić te punkty, które wypadłyby poza obszar mapy, biegnij

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

Aby przekonwertować współrzędne geodezyjne (długość, szerokość, wysokość) w pliku old.dat z NAD27
Odniesienie CONUS (identyfikator odniesienia 131, który wykorzystuje elipsoidę Clarke-1866) do WGS 84, uruchom

gmt mapproject stary.dat -Th131 > nowy.dat

Aby obliczyć najbliższą odległość (w km) między każdym punktem w pliku wejściowym quakes.dat
oraz segmenty linii podane w wielosegmentowym pliku ASCII coastline.xy, run

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

OGRANICZENIA

Prostokątny obszar wejściowy ustawiony za pomocą -R generalnie zostanie odwzorowany na nieprostokątny
siatka. Chyba że -C jest ustawiony, skrajny lewy punkt na tej siatce ma wartość x = 0.0, a
najniższy punkt będzie miał wartość y = 0.0. Tak więc, zanim zdigitalizujesz mapę, przeprowadź ekstremum
współrzędne mapy przez mapaprojekt używając odpowiedniej skali i zobacz jakie wartości (x,y)
są mapowane. Użyj tych wartości podczas konfigurowania digitalizacji, aby mieć
transformacja odwrotna działa poprawnie lub alternatywnie użyj Awk skalować i przesuwać
(x,y) wartości przed przekształceniem.

W przypadku niektórych projekcji można zastosować rozwiązanie sferyczne, mimo że użytkownik wybrał
elipsoida. Dzieje się tak, gdy użytkownicy -R ustawienie oznacza region, który wykracza poza domenę
w którym obowiązują rozwinięcia w szereg elipsoidalny. Oto warunki: (1)
Stożkowy stożkowy Lamberta (-J L) i Albersa o równym obszarze (-JB) użyje rozwiązania sferycznego
gdy skala mapy przekracza 1.0E7. (2) Merkator poprzeczny (-JT) i UTM (-JU) będzie
użyj rozwiązania sferycznego, gdy zachodnia lub wschodnia granica podana w -R wiecej niż
10 stopni od południka centralnego i (3) to samo dla Cassini (-JC) ale z limitem
tylko 4 stopnie.

ELIPSOIDY ROLNICZE SFEROIDY

GMT użyje formuł elipsoidalnych, jeśli zostaną zaimplementowane, a użytkownik wybrał
elipsoida jako kształt odniesienia (patrz PROJ_ELLIPSOID). Użytkownik musi być świadomy kilku
potencjalne pułapki: (1) W przypadku niektórych prognoz, takich jak Transverse Mercator, Albers i
Stożek konforemny Lamberta używamy wyrażeń elipsoidalnych, gdy obszary mapowane są
małe i przejść do wyrażeń sferycznych (i zastępując odpowiednią pomocniczą)
szerokości geograficzne) dla większych map. Wzory elipsoidalne stosuje się w następujący sposób: (a) Poprzeczny
Mercator: gdy wszystkie punkty znajdują się w obrębie 10 stopni od południka centralnego, (b) rzuty stożkowe
gdy zakres podłużny jest mniejszy niż 90 stopni, (c) rzut Cassini, gdy wszystkie punkty
znajdują się w zakresie 4 stopni od południka środkowego. (2) Kiedy próbujesz dopasować jakiś historyczny
dane (np. współrzędne uzyskane z określoną projekcją i pewnym odniesieniem)
elipsoida) może się okazać, że GMT daje wyniki, które są nieco inne. Jeden prawdopodobny
źródłem tej niezgodności jest to, że w starszych obliczeniach często używano mniej znaczących cyfr. Do
na przykład przykłady Snydera często wykorzystują elipsoidę Clarke'a 1866 (określoną przez niego jako posiadającą
spłaszczenie f = 1/294.98). Z f otrzymujemy mimośród do kwadratu równy 0.00676862818
(tego używa GMT), podczas gdy Snyder zaokrągla i używa 0.00676866. Ta różnica może
podać rozbieżności rzędu kilkudziesięciu cm. Jeśli chcesz odtworzyć rzutowane współrzędne
przy tej nieco innej ekscentryczności należy określić własną elipsoidę za pomocą
takie same parametry jak Clarke 1866, ale z f = 1/294.97861076. Pamiętaj też, że starsze
dane mogą odnosić się do różnych układów odniesienia i chyba że wiesz, który układ odniesienia został użyty i
przekonwertować wszystkie dane na wspólny punkt odniesienia, mogą wystąpić niedopasowania od dziesiątek do setek
metrów. (3) Na koniec pamiętaj, że PROJ_SCALE_FACTOR ma pewne wartości domyślne dla niektórych
prognozy, więc być może będziesz musiał zmienić ustawienie w celu dopasowania uzyskanych wyników
z innymi ustawieniami.

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