mapprojectgmt - Online în cloud

PROGRAM:

NUME

mapproject - Efectuați transformări de hărți înainte și inversă, conversii de date și geodezie

REZUMAT

neaproape [ tabele ] parametrii regiune [ b|B|f|F|o|O[poate0/lat0] ] [ [dx/dy] ] [ c|i|p ] [
[datums] ] [ [unitate] ] [ [x0/y0/][[+|-]unitate][+|-] ] [ ] [ linie.xy[/[+|-]unitate][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]din[/la] ] [ [nivel] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o] ]

Notă: Nu este permis niciun spațiu între indicatorul opțiunii și argumentele asociate.

DESCRIERE

neaproape citește (longitudine, latitudine) poziții din tabele [sau intrare standard] și
calculează coordonatele (x,y) folosind proiecția hărții și scalele specificate. Opțional, acesta
poate citi (x,y) poziții și poate calcula valori (longitudine, latitudine) făcând inversul
transformare. Aceasta poate fi folosită pentru a transforma puncte liniare (x,y) obținute prin digitizare
o hartă a proiecției cunoscute la coordonatele geografice. Poate calcula, de asemenea, distanțe de-a lungul
traseu, la un punct fix sau cea mai apropiată apropiere de o linie. În cele din urmă, poate fi folosit pentru a efectua
diferite conversii de date. Câmpurile de date suplimentare sunt permise după primele 2 coloane
care trebuie să aibă (longitudine,latitudine) sau (x,y). Vezi opțiunea -: despre cum să citești
(latitudine, longitudine).

NECESARE ARGUMENTE

-Jparametrii (Mai mult ...)
Selectați proiecția hărții.

-R[unitate]xmin/XMAX/ymin/ymax[R] (Mai mult ...)
Specificați regiunea de interes. Caz special pentru proiecția UTM: Dacă -C este utilizat
și -R nu este dat, atunci regiunea este setată să coincidă cu zona UTM dată, deci
pentru a păstra soluția elipsoidală completă (vezi RESTRICȚII pentru mai multe
informație).

OPTIONAL ARGUMENTE

tabel Unul sau mai multe ASCII (sau binar, vezi -bi[ncols][tip]) fișiere de tabel de date care dețin a
numărul de coloane de date. Dacă nu sunt date tabele, citim din intrarea standard.

-Ab|B|f|F|o|O[poate0/lat0]
-Af calculează azimutul (înainte) din punctul fix lon/lat la fiecare punct de date.
Utilizare -Departe pentru a obține back-azimut de la punctele de date la punctul fix. Utilizare -Ao a obține
orientări (-90/90) mai degrabă decât azimuturi (0/360). Majuscule F, B or O voi
conversia de la latitudini geodezice la latitudini geocentrice și estimarea azimutului geodezicilor
(presupunând că elipsoidul actual nu este o sferă). Dacă nu este dat niciun punct fix, atunci
calculăm azimutul (sau back-azimuth) din punctul anterior.

-C[dx/dy]
Setați centrul coordonatelor proiectate să fie în centrul de proiecție a hărții [Implicit este
colțul din stânga jos]. Opțional, adăugați decalaje în unitățile proiectate care vor fi adăugate (sau
scăzut când -I este setat) la (de la) coordonatele proiectate, cum ar fi false
estări și norduri pentru anumite zone de proiecție [0/0]. Unitatea folosită pentru
offsets este unitatea de distanță a graficului în vigoare (a se vedea PROJ_LENGTH_UNIT) dacă nu -F is
utilizate, caz în care decalajele sunt în metri.

-Dc|i|p
Înlocuiți temporar PROJ_LENGTH_UNIT și utilizați c (cm), i (inch), sau p (puncte)
in schimb. Nu poate fi folosit cu -F.

-E[datums]
Convertiți de la geodezică (lon, lat, înălțime) în Earth Centered Earth Fix (ECEF)
(x,y,z) coordonate (adăugați -I pentru conversia inversă). Adăugați ID de referință (vezi -Qd)
sau da elipsoid:dx,dy,dz Unde elipsoid poate fi un ID elipsoid (vezi -Qe) Sau
dat ca a[,*inv_f*], unde a este semiaxa majoră și inv_f este inversul
aplatizare (0 dacă este omisă). Dacă datums este sau nu dat presupunem WGS-84.

-F[unitate]
Forțare scalare 1:1, adică ieșire (sau intrare, vezi -I) datele sunt proiectate efectiv
metri. Pentru a specifica alte unități, adăugați unitatea dorită (vezi UNITĂȚI). Fără -F,
ieșirea (sau intrarea, vezi -I) sunt în unitățile specificate de PROJ_LENGTH_UNIT (dar
vedea -D).

-G[x0/y0/][[+|-]unitate][+|-]
Calculați distanțele de-a lungul căii or la punctul opțional setat cu -Gx0/y0. Adăuga
unitatea de distanta (vezi UNITATE), inclusiv c (Distanța carteziană folosind intrare
coordonate) sau C (Distanța carteziană folosind coordonatele proiectate). The C unitate
Necesită -R și -J a fi setat. Fără nici un punct fix este dat, calculăm cumul
distante de-a lungul pistei. Adăuga - pentru a obține distanța incrementală între succesive
puncte. Adăuga + pentru a specifica al 2-lea punct prin intermediul a două coloane suplimentare în fișierul de intrare.

-I Faceți transformarea inversă, adică obțineți (longitudine, latitudine) din datele (x,y).

-Llinie.xy[/[+|-]unitate][+]
Determinați cea mai scurtă distanță de la punctele de date de intrare la liniile indicate în
fișierul multisegment ASCII linie.xy. Distanța și coordonatele lui
cel mai apropiat punct va fi adăugat la ieșire ca trei coloane noi. Adăugați
unitate de distanță (vezi UNITĂȚI), inclusiv c (Distanța carteziană folosind coordonatele de intrare)
or C (Distanța carteziană folosind coordonatele proiectate). The C unitatea necesită -R și
-J a fi setat. În cele din urmă, anexează + pentru a raporta id-ul segmentului de linie și fracționalul
numărul punctului în loc de lon/lat al celui mai apropiat punct.

-N[a|c|g|m]
Convertiți din latitudini geodezice (folosind elipsoidul curent; vedeți PROJ_ELLIPSOID)
la una dintre cele patru latitudini auxiliare diferite (longitudinile nu sunt afectate). Alege
din autalic, conformal, geocentric, și mlatitudini eridionale [geocentrice]. Utilizare -I
a converti de la latitudini auxiliare la latitudini geodezice.

-Q[d|e Listați toți parametrii de proiecție. Pentru a enumera doar datele de referință, utilizați -Qd. Doar la lista
elipsoizi, folosire -Qe.

-S Suprimați punctele care se încadrează în afara regiunii.

-T[h]din[/la]
Coordonați conversiile între date din și la folosind standardul Molodensky
transformare. Utilizare -T dacă a treia coloană de intrare are înălțimea deasupra elipsoidului [Implicit
presupune înălțimea = 0, adică pe elipsoid]. Specificați datele utilizând ID-ul de referință (vezi
-Qd) sau da elipsoid:dx,dy,dz Unde elipsoid poate fi un ID elipsoid (vezi -Qe) Sau
dat ca a[,*inv_f*], unde a este semiaxa majoră și inv_f este inversul
aplatizare (0 dacă este omisă). Dacă datums este sau nu dat presupunem WGS-84. -T poate fi
folosit împreună cu -R -J pentru a modifica data înainte de proiectarea coordonatelor
(adăuga -I pentru a aplica conversia datelor după proiecția inversă). Asigura-te ca
setarea PROJ_ELLIPSOID este corectă pentru cazul dvs.

-V[nivel] (Mai mult ...)
Selectați nivelul de verbozitate [c].

-W[w|h]
Imprimă lățimea și înălțimea hărții pe ieșire standard. Nu sunt citite fișiere de intrare. Numai la
scoateți lățimea sau înălțimea, adăugați w or h, respectiv. Unitățile de
dimensiunile pot fi modificate prin -D.

-bi[ncols][t] (Mai mult ...)
Selectați intrarea binară nativă. [Implicit este 2 coloane de intrare].

-bo[ncols][tip] (Mai mult ...)
Selectați ieșirea binară nativă. [Valoare implicită este aceeași cu intrarea].

-d[i|o]nu există date (Mai mult ...)
Înlocuiți coloanele de intrare egale nu există date cu NaN și faceți invers la ieșire.

-f[i|o]colinfo (Mai mult ...)
Specificați tipurile de date ale coloanelor de intrare și/sau de ieșire.

-g[a]x|y|d|X|Y|D|[cu]z[+|-]decalaj[U] (Mai mult ...)
Determinați golurile de date și întreruperile de linie.

-h[i|o][n][+c][+d][+rremarcă][+rtitlu] (Mai mult ...)
Omite sau produce înregistrări antet.

-icol[l][sscară][ocompensa][,...] (Mai mult ...)
Selectați coloanele de intrare (0 este prima coloană).

-ocol[,...] (Mai mult ...)
Selectați coloanele de ieșire (0 este prima coloană).

-p[x|y|z]azim/ridica[/zlevel][+wpoate0/lat0[/z0]][+vx0/y0] (Mai mult ...)
Selectați vizualizarea în perspectivă.

-s[col][a|r] (Mai mult ...)
Setați gestionarea înregistrărilor NaN.

-:[i|o] (Mai mult ...)
Schimbați prima și a doua coloană la intrare și/sau la ieșire.

-^ or doar -
Imprimă un mesaj scurt despre sintaxa comenzii, apoi iese (NOTĂ: pe Windows
foloseste doar -).

-+ or doar +
Imprimați un mesaj extins de utilizare (ajutor), inclusiv explicația oricăruia
opțiunea specifică modulului (dar nu opțiunile comune GMT), apoi iese.

-? or Nu. argumente
Apoi imprimați un mesaj complet de utilizare (ajutor), inclusiv explicația opțiunilor
iesirile.

--versiune
Tipăriți versiunea GMT și ieșiți.

--show-datadir
Imprimați calea completă către directorul de partajare GMT și ieșiți.

UNITĂȚI

Pentru unitatea de distanță pe hartă, anexați unitate d pentru gradul de arc, m pentru minut de arc și s pentru arc
al doilea, sau e pentru contor [Implicit], f pentru picior, k pentru km, M pentru mile statutare, n pentru nautice
milă, și u pentru piciorul sondajului SUA. În mod implicit, calculăm astfel de distanțe folosind o sferică
aproximare cu cercuri mari. Adăugați - la o distanță (sau unitatea nu este distanța este
dat) pentru a efectua calcule „Pământ plat” (mai rapide, dar mai puțin precise) sau pentru a adăuga înainte + la
efectuați calcule geodezice exacte (mai lente, dar mai precise).

ASCII FORMAT PRECIZIE

Formatele de ieșire ASCII ale datelor numerice sunt controlate de parametrii din dvs gmt.conf
fişier. Longitudinea și latitudinea sunt formatate conform FORMAT_GEO_OUT, în timp ce altele
valorile sunt formatate conform FORMAT_FLOAT_OUT. Fiți conștienți de faptul că formatul efectiv poate
duce la pierderea preciziei la ieșire, ceea ce poate duce la diverse probleme în aval. Dacă
descoperiți că rezultatul nu este scris cu suficientă precizie, luați în considerare trecerea la binar
ieșire (-bo dacă este disponibil) sau specificați mai multe zecimale utilizând setarea FORMAT_FLOAT_OUT.

EXEMPLE

Pentru a converti coordonatele UTM în metri în locații geografice, dat fiind un fișier utm.txt și
cunoscând zona UTM (și zona sau emisfera), încercați

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

Pentru a transforma un fișier cu (longitudine, latitudine) în poziții (x,y) în cm pe un Mercator
grilă pentru o scară dată de 0.5 cm pe grad, rulați

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

Pentru a transforma mai multe fișiere binare, cu două coloane, de precizie dublă cu (latitudine, longitudine)
în pozițiile (x,y) în inci pe o grilă Mercator transversală (longitudine centrală 75W) pentru
scară = 1:500000 și suprimați acele puncte care ar fi în afara zonei hărții, alergați

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

Pentru a converti coordonatele geodezice (lon, lat, height) în fișierul old.dat din NAD27
Date CONUS (ID de referință 131 care utilizează elipsoidul Clarke-1866) la WGS 84, rulați

gmt mapproject old.dat -Th131 > new.dat

Pentru a calcula cea mai apropiată distanță (în km) dintre fiecare punct din fișierul de intrare quakes.dat
iar segmentele de linie date în fișierul ASCII multisegment coastline.xy, rulează

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

RESTRICȚII

Regiunea de intrare dreptunghiulară setată cu -R va fi, în general, mapat într-o formă nedreptunghiulară
grilă. Dacă nu -C este setat, punctul cel mai din stânga al acestei grile are xvalue = 0.0 și
punctul cel mai de jos va avea valoarea y = 0.0. Astfel, înainte de a digitiza o hartă, rulați extrema
coordonatele hărții prin neaproape folosind scara adecvată și vedeți ce valori (x,y).
acestea sunt mapate pe. Utilizați aceste valori atunci când configurați pentru digitizare pentru a avea
transformarea inversă funcționează corect sau, alternativ, folosește Wow pentru a scala și a schimba
(x,y) valori înainte de transformare.

Pentru anumite proiecții, poate fi utilizată o soluție sferică, în ciuda faptului că utilizatorul a selectat un
elipsoid. Acest lucru se întâmplă atunci când utilizatorii -R setarea implică o regiune care depășește domeniul
în care sunt valabile expansiunile serie elipsoidale. Acestea sunt condițiile: (1)
Lambert conform conic (-JL) și Albers Equal-Area (-JB) va folosi soluția sferică
când scara hărții depășește 1.0E7. (2) Mercator transversal (-JT) și UTM (-JU) va va
utilizați soluția sferică atunci când limita de vest sau de est este dată în -R este mai mult decât
10 grade de meridianul central și (3) același lucru pentru Cassini (-JC) dar cu o limită de
doar 4 grade.

ELIPSOIZII AND SFEROIZII

GMT va folosi formule elipsoidale dacă acestea sunt implementate și utilizatorul a selectat un
elipsoid ca formă de referință (vezi PROJ_ELLIPSOID). Utilizatorul trebuie să fie conștient de câteva
potențiale capcane: (1) Pentru unele proiecții, cum ar fi Transverse Mercator, Albers și
Conica conformă a lui Lambert folosim expresiile elipsoidale atunci când zonele mapate sunt
mic și treceți la expresiile sferice (și înlocuind auxiliarul corespunzător
latitudini) pentru hărți mai mari. Formulele elipsoidale se folosesc după cum urmează: (a) Transversale
Mercator: Când toate punctele sunt la 10 grade de meridianul central, (b) proiecții conice
când intervalul longitudinal este mai mic de 90 de grade, (c) proiecția Cassini când toate punctele
sunt la 4 grade de meridianul central. (2) Când încercați să vă potriviți cu ceva istoric
date (de exemplu, coordonate obținute cu o anumită proiecție și o anumită referință
elipsoid) puteți descoperi că GMT oferă rezultate ușor diferite. Unul probabil
sursa acestei nepotriviri este că calculele mai vechi foloseau adesea cifre mai puțin semnificative. Pentru
De exemplu, exemplele lui Snyder folosesc adesea elipsoidul Clarke 1866 (definit de el ca având
o aplatizare f = 1/294.98). Din f obținem că excentricitatea la pătrat este 0.00676862818
(acesta este ceea ce folosește GMT), în timp ce Snyder rotunjește și folosește 0.00676866. Această diferență poate
dați discrepanțe de câteva zeci de cm. Dacă trebuie să reproduceți coordonatele proiectate
cu această excentricitate ușor diferită, ar trebui să specificați propriul dvs. elipsoid cu
aceiași parametri ca și Clarke 1866, dar cu f = 1/294.97861076. De asemenea, rețineți că mai în vârstă
datele pot fi referite la diferite date și, cu excepția cazului în care știți care date a fost utilizată și
convertiți toate datele într-o dată comună, este posibil să aveți nepotriviri de la zeci la sute
metri. (3) În cele din urmă, rețineți că PROJ_SCALE_FACTOR au anumite valori implicite pentru unii
proiecții, astfel încât este posibil să trebuiască să înlocuiți setarea pentru a se potrivi cu rezultatele produse
cu alte setări.

Utilizați mapprojectgmt online folosind serviciile onworks.net