gmtspatialgmt - En ligne dans le Cloud

PROGRAMME:

Nom

gmtspatial - Effectuer des opérations géospatiales sur des lignes et des polygones

SYNOPSIS

gmtspatial [ table ] [ [amin_dist][unité]] [ ] [
[+ffilet][+aamax][+ddmax][+c|Ccmax][+sfait] ] [ +|- ] [ [l] ] [ -I[e|i] ] [
pfichier[+a][+pCommencer][+r][+z] ] [ [[-|+]*unité*][+h][+l][+p] ] [ région ] [ i|u|s|j ] [
[clippolygone] ] [ [niveau] ] [ -b] [ -d] [ -f] [ -g] [
-h] [ -i] [ -o] [ -:[i|o] ]

Remarque: Aucun espace n'est autorisé entre l'indicateur d'option et les arguments associés.

DESCRIPTION

gmtspatial lit un ou plusieurs fichiers de données (qui peuvent être des fichiers multisegments) qui contiennent
polygones fermés et exploite ces polygones de la manière spécifiée. Les opérations comprennent
calcul de surface, inversions manuelles et intersections de polygones.

REQUIS ARGUMENTS

Aucun.

EN OPTION ARGUMENTS

table Un ou plusieurs ASCII (ou binaire, voir -bi[ncols][type]) fichier(s) de table de données contenant un
nombre de colonnes de données. Si aucun tableau n'est fourni, nous lisons à partir de l'entrée standard.

-UNE[unemin_dist][unité]
Effectuez une analyse spatiale du voisin le plus proche (NN) : déterminez le voisin le plus proche de
chaque point et rapporter les distances NN et les identifiants de points impliqués dans chaque paire (identifiants
sont les numéros d'enregistrement d'entrée commençant à 0). Utilisation -Aa décimer un ensemble de données afin
qu'aucune distance NN n'est inférieure au seuil min_dist. Dans ce cas on écrit
les coordonnées (éventuellement moyennées) et les distances NN mises à jour et les identifiants de point.
Un nombre de points négatif signifie que le point d'origine a été remplacé par une moyenne pondérée
(la valeur d'ID absolue donne l'ID du premier ID de point d'origine à inclure
dans la moyenne.)

-C Découpe les polygones à la région de la carte, y compris la limite de la carte au polygone si nécessaire.
Le résultat est un polygone fermé (voir -T pour la troncature à la place). A besoin -R.

-D[+ffilet][+unamax][+ddmax][+c|Ccmax][+sfait]
Vérifiez les doublons parmi les lignes ou les polygones d'entrée, ou, si filet est donné via
+f, vérifiez si les entités en entrée existent déjà parmi les entités dans filet. Nous
considérer les cas de correspondances exactes (même nombre et coordonnées) et approximatives
(la distance moyenne entre les points les plus proches de deux entités est inférieure à un seuil).
Nous considérons également que certaines fonctionnalités peuvent avoir été inversées. Les caractéristiques sont considérées
correspondances approximatives si leur distance minimale est inférieure à dmax [0] (voir UNITÉS) et
leur proximité (définie comme le rapport entre la distance moyenne entre les
éléments divisés par leur longueur moyenne) est inférieur à cmax [0.01]. Pour chaque
doublon trouvé, l'enregistrement de sortie commence par la lettre unique Y (correspondance exacte) ou
~ (correspondance approximative). Si les deux segments correspondants diffèrent en longueur de plus d'un
facteur de 2 alors nous considérons que le doublon est soit un sous-ensemble (-) soit un sur-ensemble
(+). Enfin, nous notons également si deux lignes sont le résultat de la division d'un continu
ligne à travers la ligne de date (|). Pour les polygones, nous considérons également le fractionnaire
différence dans les zones; les doublons doivent différer de moins de amax [0.01]. Par défaut,
nous calculons la séparation moyenne des lignes. Utilisation +Ccmmin à la place calculer la ligne médiane
séparation et donc une valeur de proximité robuste. Par défaut également, nous considérons tous
distances entre les points d'une ligne et d'une autre. Ajouter +p limiter la comparaison
aux points qui se projettent perpendiculairement aux points de l'autre droite (et non à ses
extension).

-E+|- ]
Réinitialiser la latéralité de tous les polygones pour correspondre à la donnée + (sens antihoraire) ou -
(dans le sens des aiguilles d'une montre). Implique -Q+.

-F[l] Forcer les données d'entrée à devenir des polygones en sortie, c'est-à-dire les fermer explicitement sinon
déjà fermé. Ajoutez éventuellement l pour forcer la géométrie de la ligne.

-je[e|je]
Déterminer les emplacements d'intersection entre toutes les paires de polygones. Ajouter i à
calculer uniquement les croisements internes (c'est-à-dire les polygones auto-sécants) ou e à seulement
calculer les croisements externes (c'est-à-dire entre les paris de polygones) [la valeur par défaut est les deux].

-Npfichier[+a][+pCommencer][+r][+z]
Déterminez si un (ou tous, avec +a) les points de chaque entité dans les données d'entrée sont
à l'intérieur de l'un des polygones donnés dans le pfichier. Si à l'intérieur, indiquez quel polygone
il est; l'ID du polygone est soit tiré de la valeur aspatiale attribuée à Z, le
en-tête de segment (premier -Z, puis -L sont scannés), ou il se voit attribuer le
numéro initialisé à Commencer [0]. Par défaut, les segments d'entrée qui sont
trouvés à l'intérieur d'un polygone sont écrits sur stdout avec l'ID de polygone encodé dans
l'en-tête de segment comme -ZID. Sinon, ajoutez +r pour simplement signaler quel polygone
contient une caractéristique ou +z pour que les ID soient ajoutés en tant que colonne de données supplémentaire à la sortie.
Les segments qui ne se trouvent pas à l'intérieur d'un polygone ne sont pas écrits. Si plus d'un
polygone contient le même segment que nous sautons le deuxième (et plus loin) scénario.

-Q[[-|+]*unité*][+h][+l][+p]
Mesurez l'aire de tous les polygones ou la longueur des segments de ligne. Utilisation -Q+h d'ajouter le
zone à chaque en-tête de segment de polygones [La valeur par défaut écrit simplement la zone dans la sortie standard].
Pour les polygones, nous calculons également l'emplacement du centre de gravité tandis que pour les données de ligne, nous calculons
la position médiane (demi-longueur). Ajouter une unité de distance pour sélectionner l'unité
utilisé (voir UNITÉS). Notez que la zone dépendra du réglage actuel de
PROJ_ELLIPSOID ; cela devrait être un ellipsoïde récent pour obtenir des résultats précis. Les
le centroïde est calculé en utilisant la moyenne des vecteurs cartésiens 3-D constituant le
sommets des polygones, tandis que la surface est obtenue via une projection à surface égale. Pour
longueurs de ligne que vous pouvez ajouter -|+ à l'unité et le calcul utilisera Flat
Algorithmes terrestres ou géodésiques, respectivement [La valeur par défaut est les distances des grands cercles]. Par
par défaut, nous considérons les polygones ouverts comme des lignes. Ajouter +p pour fermer des polygones ouverts et
donc considérer toutes les entrées comme des polygones, ou ajouter +l considérer toutes les entrées comme des lignes,
même si fermé.

-R[unité]ouest/est/sud/nord[/zmin/zmax][r]
ouest, est, sudet nord spécifiez la région d'intérêt, et vous pouvez spécifier
en degrés décimaux ou au format [+-]dd:mm[:ss.xxx][W|E|S|N]. Ajouter r si inférieur
les coordonnées de gauche et en haut à droite de la carte sont données au lieu de w/e/s/n. Les deux
raccourcis -Rg et -Rd signifie domaine global (0/360 et -180/+180 en longitude
respectivement, avec -90/+90 en latitude). Alternativement pour la création de grille, donnez
Rcodelon/lat/nx/ny, Où code est une combinaison de 2 caractères de L, C, R (pour gauche,
au centre ou à droite) et T, M, B pour le haut, le milieu ou le bas. par exemple, BL pour le coin inférieur gauche.
Ceci indique quel point sur une région rectangulaire le lon/lat la coordonnée fait référence
à, et les dimensions de la grille nx et ny avec des espacements de grille via -I est utilisé pour créer
la région correspondante. Vous pouvez également spécifier le nom d'un fichier de grille existant
et la -R les paramètres (et l'espacement de la grille, le cas échéant) sont copiés à partir de la grille.
En utilisant -Runité attend des coordonnées projetées (cartésiennes) compatibles avec celles choisies -J
et nous projetons inversement de déterminer la région géographique rectangulaire réelle. Pour
vue de perspective (-p), éventuellement ajouter /zmin/zmax. En cas de vue en perspective
(-p), une plage z (zmin, zmax) peut être ajouté pour indiquer la troisième dimension. Cette
doit être fait uniquement lors de l'utilisation du -Jz option, pas lorsque vous utilisez uniquement l'option -p option.
Dans ce dernier cas, une vue en perspective de l'avion est tracée, sans troisième
dimension. Délimite les polygones à la région de la carte, y compris la limite de la carte au polygone
comme requis. Le résultat est un polygone fermé.

-Si|j|s|u
Traitement spatial des polygones. Choisissez parmi -Si qui renvoie l'intersection de
polygones (fermés), -Le sien qui renvoie l'union des polygones (fermé), -H.H Qui va
diviser les polygones qui chevauchent la ligne de date, et -Sj qui rejoindra des polygones qui
ont été divisés par la ligne de date. Remarque : uniquement -H.H a été mis en place.

-T[clippolygone]
Tronquer les polygones par rapport au polygone spécifié donné, ce qui peut entraîner une ouverture
polygones. Si aucun argument n'est donné à -T nous créons un polygone de découpage à partir de -R qui
alors est requis. Notez que lorsque le -R l'écrêtage est en vigueur, nous allons également rechercher
polygones de longueur 4 ou 5 qui correspondent exactement à la -R polygone de découpage.

-V[niveau] (plus ...)
Sélectionnez le niveau de verbosité [c].

-bi[ncols][t] (plus ...)
Sélectionnez l'entrée binaire native. [La valeur par défaut est de 2 colonnes d'entrée].

-bo[ncols][type] (plus ...)
Sélectionnez la sortie binaire native. [La valeur par défaut est identique à l'entrée].

-d[i|o]pas de données (plus ...)
Remplacer les colonnes d'entrée égales pas de données avec NaN et faire l'inverse en sortie.

-f[je|o]colinfo (plus ...)
Spécifiez les types de données des colonnes d'entrée et/ou de sortie.

-g[a]x|y|d|X|Y|D|[avec]z[+|-]écart[u] (plus ...)
Déterminez les lacunes dans les données et les sauts de ligne.

-h[je|o][n][+c][+d][+rremarque][+rtitre] (plus ...)
Ignorer ou produire des enregistrements d'en-tête.

-icols[l][sen échelon][ocompenser][,...] (plus ...)
Sélectionnez les colonnes d'entrée (0 est la première colonne).

-ocols[,...] (plus ...)
Sélectionnez les colonnes de sortie (0 est la première colonne).

- :[i|o] (plus ...)
Permutez la 1ère et la 2ème colonne en entrée et/ou en sortie.

-^ or juste -
Imprime un court message sur la syntaxe de la commande, puis quitte (REMARQUE : sous Windows
utiliser juste -).

-+ or juste +
Imprimez un message d'utilisation détaillé (aide), y compris l'explication de tout
option spécifique au module (mais pas les options communes GMT), puis se ferme.

-? or aucune arguments
Imprimez un message d'utilisation (aide) complet, y compris l'explication des options, puis
sorties.

--version
Imprimer la version GMT et quitter.

--show-datadir
Affichez le chemin complet vers le répertoire de partage GMT et quittez.

UNITÉS

Pour l'unité de distance de la carte, ajoutez unité d pour le degré d'arc, m pour la minute d'arc, et s pour arc
deuxième, ou e pour le compteur [Par défaut], f pour pied, k pour km, M pour mille terrestre, n pour nautique
mille, et u pour pied d'arpentage américain. Par défaut, nous calculons de telles distances en utilisant un
approximation avec les grands cercles. Préfixer - à une distance (ou l'unité est aucune distance n'est
donné) pour effectuer des calculs « Terre plate » (plus rapide mais moins précis) ou ajouter + à
effectuer des calculs géodésiques exacts (plus lents mais plus précis).

ASCII Format PRÉCISION

Les formats de sortie ASCII des données numériques sont contrôlés par les paramètres de votre gmt.conf
déposer. La longitude et la latitude sont formatées selon FORMAT_GEO_OUT, tandis que d'autres
les valeurs sont formatées selon FORMAT_FLOAT_OUT. Sachez que le format en vigueur peut
conduire à une perte de précision dans la sortie, ce qui peut conduire à divers problèmes en aval. Si
vous trouvez que la sortie n'est pas écrite avec suffisamment de précision, envisagez de passer en binaire
sortir (-bo si disponible) ou spécifiez plus de décimales à l'aide du paramètre FORMAT_FLOAT_OUT.

EXEMPLE

Pour transformer toutes les lignes du fichier multisegment lines.txt en polygones fermés, exécutez

gmt lignes spatiales.txt -F > polygones.txt

Pour calculer l'aire de tous les polygones géographiques dans le fichier multisegment polygons.txt, exécutez

gmt spatial polygons.txt -Q > secteurs.txt

Mêmes données, mais maintenant orientez tous les polygones dans le sens inverse des aiguilles d'une montre et écrivez leurs aires sur
les en-têtes de segment, exécutez

gmt spatial polygons.txt -Q+h -E+ > area.txt

Pour déterminer les intersections entre les polygones A.txt et B.txt, exécutez

gmt spatial A.txt B.txt -Ce > crossovers.txt

Utiliser gmtspatialgmt en ligne en utilisant les services onworks.net