mapprojectgmt - Online in de cloud

PROGRAMMA:

NAAM

mapproject - Voer voorwaartse en inverse kaarttransformaties, datumconversies en geodesie uit

KORTE INHOUD

kaartproject [ tafels ] parameters regio [ b|B|f|F|o|O[kan0/lat0] ] [ [dx/dy] ] [ c|i|p ] [
[datum] ] [ [eenheid] ] [ [x0/y0/][[+|-]eenheid][+|-] ] [ ] [ lijn.xy[/[+|-]eenheid][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]oppompen van[/naar] ] [ [niveau] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o] ]

Opmerking: Er is geen spatie toegestaan ​​tussen de optievlag en de bijbehorende argumenten.

PRODUCTBESCHRIJVING

kaartproject leest (lengtegraad, breedtegraad) posities van tafels [of standaard invoer] en
berekent (x,y)-coördinaten met behulp van de gespecificeerde kaartprojectie en schalen. Optioneel kan het
kan (x,y) posities lezen en waarden (lengtegraad, breedtegraad) berekenen door het omgekeerde te doen
transformatie. Dit kan worden gebruikt om lineaire (x,y) punten te transformeren die zijn verkregen door digitalisering
een kaart met bekende projectie naar geografische coördinaten. Kan ook afstanden mee berekenen
spoor, naar een vast punt of de dichtste nadering van een lijn. Ten slotte kan worden gebruikt om uit te voeren
verschillende datumconversies. Extra gegevensvelden zijn toegestaan ​​na de eerste 2 kolommen
die (lengtegraad,breedtegraad) of (x,y) moet hebben. Zie optie -: over hoe te lezen
(breedtegraad, lengtegraad) bestanden.

VERPLICHT ARGUMENTEN

-Jparameters (meer ...)
Selecteer kaartprojectie.

-R[eenheid]x min/xmax/ymin/ymax[R] (meer ...)
Geef de regio van interesse op. Speciaal geval voor de UTM-projectie: If -C is gebruikt
en -R is niet gegeven, dan wordt de regio ingesteld om samen te vallen met de gegeven UTM-zone dus
om de volledige ellipsoïde oplossing te behouden (zie BEPERKINGEN voor meer)
informatie).

OPTIONELE ARGUMENTEN

tafel Een of meer ASCII (of binaire, zie -bi[ncols][ type dan: ]) gegevenstabelbestand(en) met a
aantal gegevenskolommen. Als er geen tabellen worden gegeven, lezen we uit de standaardinvoer.

-Ab|B|f|F|o|O[kan0/lat0]
-Af berekent de (voorwaartse) azimut vanaf een vast punt lang/latijn naar elk datapunt.
Te gebruiken -Weg om back-azimuth van datapunten naar een vast punt te krijgen. Gebruik maken van -oa te krijgen
oriëntaties (-90/90) in plaats van azimuts (0/360). hoofdletters F, B or O wil
converteer van geodetische naar geocentrische breedtegraden en schat het azimut van geodeten
(ervan uitgaande dat de huidige ellipsoïde geen bol is). Als er geen vast punt is gegeven dan
we berekenen de azimut (of back-azimut) van het vorige punt.

-C[dx/dy]
Stel het middelpunt van de geprojecteerde coördinaten in op het projectiecentrum van de kaart [Standaard is
linksonderin]. Voeg optioneel offsets toe in de geprojecteerde eenheden die moeten worden toegevoegd (of
afgetrokken wanneer -I is ingesteld) naar (van) de geprojecteerde coördinaten, zoals false
oosten en noorden voor bepaalde projectiezones [0/0]. De eenheid die wordt gebruikt voor de
offsets is de eenheid van de plotafstand die van kracht is (zie PROJ_LENGTH_UNIT) tenzij -F is
gebruikt, in welk geval de offsets in meters zijn.

-Dc|i|p
Tijdelijk overschrijven van PROJ_LENGTH_UNIT en gebruiken c (cm), i (inch), of p (punten)
in plaats daarvan. Kan niet worden gebruikt met -F.

-E[datum]
Converteren van geodetisch (lon, lat, height) naar Earth Centered Earth Fixed (ECEF)
(x,y,z) coördinaten (add -I voor de inverse conversie). Datum ID toevoegen (zie -Qd)
of geef ellipsoïde:dx,dy,dz WAAR ellipsoïde kan een ellipsoïde ID zijn (zie -Qe) Of
gegeven als a[,*inv_f*], waar a is de halve lange as en inv_f is het omgekeerde
afvlakking (0 indien weggelaten). Indien datum is - of niet gegeven, we gaan uit van WGS-84.

-F[eenheid]
Forceer 1:1 schaling, dwz output (of input, zie -I) gegevens zijn in werkelijkheid geprojecteerd
meter. Om andere eenheden te specificeren, voegt u de gewenste eenheid toe (zie EENHEDEN). Zonder -F,
de output (of input, zie -I) zijn in de eenheden gespecificeerd door PROJ_LENGTH_UNIT (maar
zien -D).

-G[x0/y0/][[+|-]eenheid][+|-]
Bereken afstanden langs het spoor or naar het optionele punt ingesteld met -Gx0/j0. toevoegen
de afstandseenheid (zie EENHEDEN), inclusief c (Cartesiaanse afstand met invoer)
coördinaten) of C (Cartesiaanse afstand met behulp van geprojecteerde coördinaten). De C eenheid
vereist -R en -J er klaar voor zijn. Zonder vast punt wordt gegeven berekenen we cumuleren
afstanden langs het spoor. toevoegen - om een ​​incrementele afstand te verkrijgen tussen opeenvolgende
punten. toevoegen + om het 2e punt te specificeren via twee extra kolommen in het invoerbestand.

-I Voer de inverse transformatie uit, dwz haal (lengtegraad, breedtegraad) uit (x,y) gegevens.

-Llijn.xy[/[+|-]eenheid][+]
Bepaal de kortste afstand van de invoergegevenspunten tot de lijn(en) gegeven in
het ASCII-multisegmentbestand lijn.xy. De afstand en de coördinaten van de
het dichtstbijzijnde punt wordt als drie nieuwe kolommen aan de uitvoer toegevoegd. Voeg de toe
afstandseenheid (zie EENHEDEN), inclusief c (Cartesiaanse afstand met invoercoördinaten)
or C (Cartesiaanse afstand met behulp van geprojecteerde coördinaten). De C eenheid vereist: -R en
-J er klaar voor zijn. Voeg ten slotte toe + om het lijnsegment-ID en de fractionele . te rapporteren
puntnummer in plaats van lon/lat van het dichtstbijzijnde punt.

-N[a|c|g|m]
Converteren van geodetische breedtegraden (met behulp van de huidige ellipsoïde; zie PROJ_ELLIPSOID)
naar een van de vier verschillende hulpbreedtegraden (lengtegraden worden niet beïnvloed). Kiezen
oppompen van authalisch, conformeel, geocentrisch, en meridionale breedtegraden [geocentrisch]. Gebruik maken van -I
om te zetten van hulpbreedten naar geodetische breedtegraden.

-Q[d|e Maak een lijst van alle projectieparameters. Gebruik . om alleen datums weer te geven -Qd. Alleen lijst
ellipsoïden, gebruik -Qe.

-S Onderdruk punten die buiten de regio vallen.

-NS]oppompen van[/naar]
Coördinatenconversies tussen datums oppompen van en naar met behulp van de standaard Molodensky
transformatie. Gebruik maken van -NS als de 3e invoerkolom een ​​hoogte heeft boven de ellipsoïde [Standaard
neemt hoogte = 0 aan, dwz op de ellipsoïde]. Specificeer datums met behulp van de datum ID (zie
-Qd) of geef ellipsoïde:dx,dy,dz WAAR ellipsoïde kan een ellipsoïde ID zijn (zie -Qe) Of
gegeven als a[,*inv_f*], waar a is de halve lange as en inv_f is het omgekeerde
afvlakking (0 indien weggelaten). Indien datum is - of niet gegeven, we gaan uit van WGS-84. -T kan zijn
gebruikt in combinatie met -R -J om de datum vóór coördinaatprojectie te wijzigen
(toevoegen -I om de datumconversie toe te passen na de inverse projectie). Zeker weten dat
de PROJ_ELLIPSOID-instelling is correct voor uw geval.

-V[niveau] (meer ...)
Selecteer breedsprakigheidsniveau [c].

-W[w|h]
Drukt kaartbreedte en hoogte af op standaarduitvoer. Er worden geen invoerbestanden gelezen. om alleen
voer de breedte of de hoogte uit, voeg toe w or h, respectievelijk. De eenheden van de
afmetingen kunnen gewijzigd worden via -D.

-bi[ncols][T] (meer ...)
Selecteer native binaire invoer. [Standaard is 2 invoerkolommen].

-bo[ncols][ type dan: ] (meer ...)
Selecteer native binaire uitvoer. [Standaard is hetzelfde als invoer].

-d[i|o]geen informatie (meer ...)
Vervang invoerkolommen die gelijk zijn aan geen informatie met NaN en doe het omgekeerde bij uitvoer.

-f[i|o]colinfo (meer ...)
Specificeer gegevenstypen van invoer- en/of uitvoerkolommen.

-g[a]x|y|d|X|Y|D|[col]z[+|-]kloof[U] (meer ...)
Bepaal hiaten in de gegevens en regeleinden.

-h[i|o][n][+c][+d][+ropmerking][+rtitel] (meer ...)
Koprecord(s) overslaan of produceren.

-icols[l][sschaal][Ocompenseren][,...] (meer ...)
Selecteer invoerkolommen (0 is de eerste kolom).

-ocols[,...] (meer ...)
Selecteer uitvoerkolommen (0 is de eerste kolom).

-p[x|y|z]Azim/leerling[/z-niveau][+wkan0/lat0[/z0]][+vx0/y0] (meer ...)
Selecteer perspectiefweergave.

-S[cols][a|r] (meer ...)
Afhandeling van NaN-records instellen.

-:[i|o] (meer ...)
Wissel 1e en 2e kolom om op input en/of output.

-^ or voor slechts -
Druk een kort bericht af over de syntaxis van de opdracht en sluit vervolgens af (OPMERKING: op Windows
gebruik gewoon -).

-+ or voor slechts +
Druk een uitgebreid gebruiks(help)bericht af, inclusief de uitleg van eventuele
module-specifieke optie (maar niet de algemene GMT-opties), en wordt vervolgens afgesloten.

-? or geen argumenten
Druk een volledig gebruiks(help)bericht af, inclusief de uitleg van de opties, dan
uitgangen.

--versie
GMT-versie afdrukken en afsluiten.

--show-datadir
Druk het volledige pad af naar de GMT-sharedirectory en sluit af.

UNITS

Voor kaartafstandseenheid, voeg toe eenheid d voor booggraad, m voor boogminuut, en s voor boog
tweede, of e voor meter [Standaard], f voor voet, k voor kilometer, M voor statuut mijl, n voor nautische
mijl, en u voor Amerikaanse onderzoeksvoet. Standaard berekenen we dergelijke afstanden met behulp van een sferische
benadering met grote cirkels. voeg toe - tot een afstand (of de eenheid is geen afstand is
gegeven) om "Flat Earth"-berekeningen uit te voeren (sneller maar minder nauwkeurig) of prepend + naar
exacte geodetische berekeningen uitvoeren (langzamer maar nauwkeuriger).

ASCII FORMAT PRECISION

De ASCII-uitvoerformaten van numerieke gegevens worden bestuurd door parameters in uw gmt.conf
het dossier. Lengte- en breedtegraad zijn opgemaakt volgens FORMAT_GEO_OUT, terwijl andere
waarden zijn opgemaakt volgens FORMAT_FLOAT_OUT. Houd er rekening mee dat het huidige formaat kan
leiden tot verlies van precisie in de uitvoer, wat stroomafwaarts tot verschillende problemen kan leiden. Indien
je vindt dat de uitvoer niet met voldoende precisie is geschreven, overweeg dan om over te schakelen naar binair
uitvoer (-bo indien beschikbaar) of geef meer decimalen op met de instelling FORMAT_FLOAT_OUT.

Voorbeelden

Om UTM-coördinaten in meters naar geografische locaties te converteren, geeft u een bestand utm.txt en
ken de UTM-zone (en zone of halfrond), probeer

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

Om een ​​bestand met (lengtegraad,breedtegraad) om te zetten in (x,y) posities in cm op een Mercator
raster voor een gegeven schaal van 0.5 cm per graad, run

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

Om verschillende 2-koloms, binaire, dubbele precisiebestanden te transformeren met (breedtegraad, lengtegraad)
in (x,y) posities in inch op een transversaal Mercatorraster (centrale lengte 75W) voor
schaal = 1:500000 en onderdruk die punten die buiten het kaartgebied zouden vallen, run

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

Om de geodetische coördinaten (lon, lat, height) in het bestand old.dat van de NAD27 om te zetten
CONUS-datum (Datum ID 131 die de Clarke-1866-ellipsoïde gebruikt) naar WGS 84, run

gmt mapproject old.dat -Th131 > new.dat

Om de dichtstbijzijnde afstand (in km) tussen elk punt in het invoerbestand quakes.dat te berekenen
en de lijnsegmenten gegeven in het multisegment ASCII-bestand kustlijn.xy, run

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

BEPERKINGEN

Het rechthoekige invoergebied ingesteld met -R zal in het algemeen worden afgebeeld in een niet-rechthoekige
rooster. Behalve -C is ingesteld, heeft het meest linkse punt op dit raster xwaarde = 0.0, en de
onderste punt heeft een y-waarde = 0.0. Dus, voordat u een kaart digitaliseert, voer het uiterste uit
kaart coördinaten door kaartproject gebruik de juiste schaal en kijk welke (x,y) waarden
ze zijn in kaart gebracht. Gebruik deze waarden bij het instellen van digitalisering om de
inverse transformatie correct werken, of als alternatief, gebruik awk opschalen en verschuiven van de
(x,y) waarden vóór transformatie.

Voor sommige projecties kan een bolvormige oplossing worden gebruikt ondanks dat de gebruiker een
ellipsoïde. Dit gebeurt wanneer de gebruikers -R instelling impliceert een regio die het domein overschrijdt
waarin de ellipsvormige reeksuitbreidingen geldig zijn. Dit zijn de voorwaarden: (1)
Lambert Conform Conisch (-JL) en Albers Equal-Area (-JB) zal de sferische oplossing gebruiken
wanneer de kaartschaal groter is dan 1.0E7. (2) Transversale Mercator (-JT) en UTM (-JU) zal zal
gebruik de sferische oplossing wanneer ofwel de west- of oostgrens wordt gegeven in -R is meer dan
10 graden vanaf de centrale meridiaan, en (3) hetzelfde voor Cassini (-JC) maar met een limiet van
slechts 4 graden.

ELLIPSODEN EN SFERODEN

GMT zal ellipsvormige formules gebruiken als deze zijn geïmplementeerd en de gebruiker een
ellipsoïde als de referentievorm (zie PROJ_ELLIPSOID). De gebruiker moet op de hoogte zijn van een paar
mogelijke valkuilen: (1) Voor sommige projecties, zoals Transversale Mercator, Albers en
Volgens de conforme kegelsnede van Lambert gebruiken we de ellipsvormige uitdrukkingen wanneer de in kaart gebrachte gebieden zijn
klein, en schakel over naar de sferische uitdrukkingen (en vervang de juiste hulp
breedtegraden) voor grotere kaarten. De ellipsoïde formules worden als volgt gebruikt: (a) Transversaal
Mercator: Als alle punten binnen 10 graden van de centrale meridiaan liggen, (b) Kegelvormige projecties
wanneer het longitudinale bereik kleiner is dan 90 graden, (c) Cassini-projectie wanneer alle punten
zijn binnen 4 graden van de centrale meridiaan. (2) Wanneer u probeert een historisch te matchen
gegevens (bijv. coördinaten verkregen met een bepaalde projectie en een bepaalde referentie)
ellipsoïde) zult u merken dat GMT resultaten geeft die iets anders zijn. een waarschijnlijke
bron van deze mismatch is dat oudere berekeningen vaak minder significante cijfers gebruikten. Voor
Snyder's voorbeelden gebruiken bijvoorbeeld vaak de Clarke 1866-ellipsoïde (door hem gedefinieerd als hebbende)
een afvlakking f = 1/294.98). Van f krijgen we de excentriciteit in het kwadraat van 0.00676862818
(dit is wat GMT gebruikt), terwijl Snyder afrondt en 0.00676866 gebruikt. Dit verschil kan
afwijkingen van enkele tientallen cm geven. Als u geprojecteerde coördinaten moet reproduceren
met deze iets andere excentriciteit, zou je je eigen ellipsoïde moeten specificeren met de
dezelfde parameters als Clarke 1866, maar met f = 1/294.97861076. Houd er ook rekening mee dat oudere
gegevens kunnen naar verschillende datums worden verwezen, en tenzij u weet welke datum is gebruikt en
converteer alle gegevens naar een gemeenschappelijk datum, u kunt mismatches van tientallen tot honderden ervaren
meter. (3) Houd er ten slotte rekening mee dat PROJ_SCALE_FACTOR bepaalde standaardwaarden heeft voor sommigen
projecties, dus het kan zijn dat u de instelling moet overschrijven om de geproduceerde resultaten te matchen
met andere instellingen.

Gebruik mapprojectgmt online met onworks.net-services