Il s'agit de la commande gmtselectgmt qui peut être exécutée dans le fournisseur d'hébergement gratuit OnWorks en utilisant l'un de nos nombreux postes de travail en ligne gratuits tels que Ubuntu Online, Fedora Online, l'émulateur en ligne Windows ou l'émulateur en ligne MAC OS
PROGRAMME:
Nom
gmtselect - Sélectionner des sous-ensembles de tables de données en fonction de plusieurs critères spatiaux
SYNOPSIS
gmtselect [ table ] [ zone_min[/niveau_min/Niveau maximum][+ag|i|s|S][+r|l][ppour cent] ] [
dist[unité]/fichier pt ] [ RAPIDE[+] ] [ [fn] ] [ fichier polygone ] [ -I[cflrsz] ] [ paramètres
] [ [p]dist[unité]/fichier de lignes ] [ valeurs de masque ] [ région ] [ m.[/max][+cavec] ] [ [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
gmtselect est un filtre qui lit les positions (longitude, latitude) des 2 premières colonnes
of fichiers [ou entrée standard] et utilise une combinaison de 1 à 6 critères pour réussir ou rejeter
les enregistrements. Les enregistrements peuvent être sélectionnés selon qu'ils sont ou non 1) à l'intérieur d'un
région rectangulaire (-R [et -J]), 2) dans dist km de n'importe quel point dans fichier pt, 3) dans dist
km de n'importe quelle ligne dans fichier de lignes, 4) à l'intérieur d'un des polygones du fichier polygone, 5) à l'intérieur
caractéristiques géographiques (basées sur les côtes), ou 6) a des valeurs z dans une plage donnée.
le sens des tests peut être inversé pour chacun de ces 6 critères en utilisant le -I option.
Voir l'option -: sur la façon de lire les fichiers (latitude, longitude).
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.
-Azone_min[/niveau_min/Niveau maximum][+ag|i|s|S][+r|l][+ppour cent]
Caractéristiques d'une superficie inférieure à zone_min en km^2 ou de niveau hiérarchique qui
est inférieur à niveau_min ou supérieur à Niveau maximum ne sera pas tracé [La valeur par défaut est
0/0/4 (toutes les fonctionnalités)]. Le niveau 2 (lacs) contient des lacs réguliers et une large rivière
des corps que nous incluons normalement comme des lacs ; ajouter +r pour simplement obtenir des lacs fluviaux ou +l
pour obtenir juste des lacs réguliers. Par défaut (+ai) nous sélectionnons la limite de la banquise comme
le littoral pour l'Antarctique ; ajouter +ag sélectionner à la place la ligne de mise à la terre
comme littoral. Pour les utilisateurs experts qui souhaitent imprimer leur propre littoral de l'Antarctique
et les îles via psxy vous pouvez utiliser +comme pour ignorer toutes les fonctionnalités GSHHG en dessous de 60S ou +aS à
à la place, ignorez toutes les fonctionnalités au nord de 60S. Enfin, ajoutez +ppour cent exclure
polygones dont le pourcentage de surface de l'entité pleine résolution correspondante est inférieur
que pour centVoir les informations GSHHG ci-dessous pour plus de détails. Ignoré sauf si -N is
défini.
-Cdist[unité]/fichier pt
Transmettre tous les enregistrements dont l'emplacement est à l'intérieur dist de l'un des points de l'ASCII
filet fichier pt. Si dist est nul alors la 3ème colonne de fichier pt doit avoir chaque point
Rayon d'influence individuel. Les distances sont cartésiennes et exprimées en unités utilisateur ; précisez
-fg pour indiquer les distances sphériques et ajouter une unité de distance (voir UNITÉS).
Sinon, si -R et -J sont utilisées, puis les coordonnées géographiques sont projetées sur
coordonnées cartographiques (en cm, pouces ou points, telles que déterminées par PROJ_LENGTH_UNIT) avant
Les distances cartésiennes sont comparées à dist.
-DRAPIDE[+]
Ignoré à moins que -N est défini. Sélectionne la résolution de l'ensemble de données du littoral à utiliser.
((f)très, (h)haut, (i)intermédiaire, (l)ow, ou (c)impoli). La résolution diminue de
~80 % entre les ensembles de données. [La valeur par défaut est l]. Ajoutez (+) pour sélectionner automatiquement une valeur inférieure
résolution si celle demandée n'est pas disponible [abandonner si non trouvée]. Remarque
que parce que les côtes diffèrent dans les détails, il n'est pas garanti qu'un point
restera à l'intérieur [ou à l'extérieur] lorsqu'une résolution différente est sélectionnée.
-E[fn] Spécifiez précisément comment les points situés sur la limite d'un polygone doivent être pris en compte. Par défaut,
ces points sont considérés comme étant à l'intérieur du polygone. n et/ou f pour changer
ce comportement pour le -F et -N options, respectivement, de sorte que les points limites soient
considéré comme étant à l'extérieur.
-Ffichier polygone
Transmettez tous les enregistrements dont l'emplacement se trouve dans l'un des polygones fermés du
fichier à segments multiples fichier polygonePour les polygones sphériques (lon, lat), assurez-vous qu'aucun
les points consécutifs sont séparés par 180 degrés ou plus en longitude. Notez que
fichier polygone doit être en ASCII, que -bi est utilisé.
-Je[cflrsz]
Inverse le sens du test pour chacun des critères spécifiés :
c sélectionner les enregistrements qui ne se trouvent PAS dans le cercle d'influence d'un point.
f sélectionnez les enregistrements qui ne se trouvent PAS à l'intérieur d'un des polygones.
l sélectionnez les enregistrements qui ne se trouvent PAS dans la distance spécifiée d'une ligne.
r sélectionner les enregistrements NON à l'intérieur de la région rectangulaire spécifiée.
s sélectionner les enregistrements NON pris en compte à l'intérieur comme spécifié par -N (Et -A, -D).
z sélectionner les enregistrements NON compris dans la plage spécifiée par -Z.
-Jparamètres (plus ...)
Sélectionnez la projection cartographique.
-L[p]dist[unité]/fichier de lignes
Transmettre tous les enregistrements dont l'emplacement est à l'intérieur dist de l'un des segments de ligne dans le
Fichier ASCII à segments multiples fichier de lignes. Si dist est nul, alors nous allons scanner chaque
sous-en-tête dans le fichier pt pour un intégré -Ddist paramètre qui définit chaque ligne
Valeur de distance individuelle. Les distances sont cartésiennes et exprimées en unités utilisateur ; précisez -fg
Pour indiquer les distances sphériques, ajoutez une unité de distance (voir UNITÉS). Alternativement,
if -R et -J sont utilisées puis les coordonnées géographiques sont projetées sur les coordonnées cartographiques
(en cm, pouces, m ou points, comme déterminé par PROJ_LENGTH_UNIT) avant cartésien
les distances sont comparées à dist. utilisation -LP pour garantir uniquement les points dont l'orthogonalité
les projections sur le segment de ligne le plus proche se situent dans les extrémités des segments
[Par défaut, les points sont considérés comme « au-delà » des extrémités de la ligne.
-Nvaleurs de masque
Transmettre tous les enregistrements dont l'emplacement se trouve à l'intérieur des entités géographiques spécifiées. Spécifier
si les enregistrements doivent être ignorés (s) ou conservés (k) en utilisant l'un des 1 formats :
-Nhumide / sec.
-Nocéan/terre/lac/île/étang.
[La valeur par défaut est s/k/s/k/s (c'est-à-dire s/k), qui passe tous les points sur la terre ferme].
-R[unité]x min/xmax/ymin/ymax[r] (plus ...)
Spécifiez la région d'intérêt. Si aucune projection cartographique n'est fournie, nous définissons implicitement
-Jx1.
-V[niveau] (plus ...)
Sélectionnez le niveau de verbosité [c].
-Zm.[/max][+cavec]
Transmettre tous les enregistrements dont la 3ème colonne (z; avec = 2) se situe dans la plage donnée ou est
NaN (utiliser -s pour ignorer les enregistrements NaN). Si max est omis, alors nous testons si z équivaut à m.
Le fichier d'entrée doit comporter au moins trois colonnes. Pour indiquer qu'il n'y a pas de limite minimale,
ou max, spécifiez un tiret (-). Si votre 3e colonne est une heure absolue, n'oubliez pas de
la quantité -f2T. Pour spécifier une autre colonne, ajoutez +cavec, et de spécifier plusieurs tests
il suffit de répéter le Z autant de fois que nécessaire pour tester les colonnes. Remarque :
plus d'un Z l'option est donnée alors le Iz l'option ne peut pas être utilisée.
-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).
-s[cols][a|r] (plus ...)
Définir la gestion des enregistrements NaN.
- :[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.
Cette note s'applique uniquement à la sortie ASCII en combinaison avec une entrée binaire ou netCDF ou la
-: option. Voir également la note ci-dessous.
REMARQUE ON TRAITEMENT DES SEMENCES ASCII CONTRIBUTION RECORDS
À moins que vous n'utilisiez le -: option, les enregistrements d'entrée ASCII sélectionnés sont copiés mot pour mot dans
sortie. Cela signifie que des options comme -foT et des paramètres comme FORMAT_FLOAT_OUT et
FORMAT_GEO_OUT n'aura aucun effet sur le résultat. En revanche, il permet
sélectionner des enregistrements au contenu diversifié, y compris des chaînes de caractères, entre guillemets ou non,
commentaires et autres contenus non numériques.
REMARQUE ON DISTANCES
Si les options -C or -L sont sélectionnées, les distances sont alors cartésiennes et en unités utilisateur ; utilisez -fg
pour impliquer des distances sphériques en km et des coordonnées géographiques (lon, lat). Alternativement,
spécifier -R et -J pour mesurer les distances cartésiennes projetées en unités cartographiques (cm, pouce ou
points, tels que déterminés par PROJ_LENGTH_UNIT).
Ce programme a évolué au fil des ans. À l'origine, le -R et -J étaient obligatoires pour
pour gérer les données géographiques, mais il existe désormais une prise en charge complète des calculs sphériques. Ainsi,
-J ne doit être utilisé que si vous souhaitez que les tests soient appliqués sur des données projetées et non sur les
coordonnées d'origine. Si -J sont utilisées les distances données via -C et -L sont projetés
les distances.
REMARQUE ON SEGMENTS
Les en-têtes de segment dans les fichiers d'entrée sont copiés dans la sortie si un ou plusieurs enregistrements d'un
Le segment réussit le test. La sélection s'effectue toujours point par point, et non segment par segment.
signifie que seuls les points d'un segment qui réussissent le test seront inclus dans le résultat. Si
vous souhaitez découper les lignes et inclure les nouveaux points de délimitation aux extrémités du segment que vous
il faut plutôt utiliser gmtspatial.
EXEMPLES
Pour extraire le sous-ensemble de données qui se trouve à moins de 300 km de l'un des points de pts.d mais
à plus de 100 km des lignes dans lines.d, courir
gmt select lonlatfile -fg -C300k/pts.d -L100/lines.d -Il > sous-ensemble
Ici, vous devez préciser -fg afin que le programme sache que vous traitez des données géographiques.
Pour conserver tous les points de data.d dans la région spécifiée, à l'exception des points sur terre (comme
déterminé par les côtes à haute résolution), utiliser
gmt select data.d -R120/121/22/24 -Dh -Nk/s > sous-ensemble
Pour renvoyer tous les points dans quakes.d qui sont à l'intérieur ou sur le polygone sphérique lonlatpath.d,
Essai
gmt select quakes.d -Flonlatpath.d -fg > subset1
Pour renvoyer tous les points dans stations.d qui sont à moins de 5 cm du point dans origin.d pour un
certaine projection, essayez
gmt sélectionne les stations.d -C5/origine.d -R20/50/-10/20 -JM20c \
--PROJ_LENGTH_UNIT=cm > sous-ensemble 2
GSHHS INFORMATION
La base de données du littoral est GSHHG (anciennement GSHHS) qui est compilée à partir de trois sources :
World Vector Shorelines (WVS), CIA World Data Bank II (WDBII) et Atlas de la cryosphère
(AC, pour l'Antarctique uniquement). Hormis l'Antarctique, tous les polygones de niveau 1 (océan-terre
frontière) sont dérivés du WVS plus précis tandis que tous les polygones de niveau supérieur (niveau
2-4, représentant terre/lac, lac/île-dans-lac, et
limites d'île dans un lac/lac dans une île dans un lac) sont tirées de WDBII. L'Antarctique
les côtes se déclinent en deux saveurs : front de glace ou ligne d'échouage, sélectionnable via le -A option.
Une grande partie du traitement a eu lieu pour convertir les données WVS, WDBII et AC en une forme utilisable pour
GMT : assemblage de polygones fermés à partir de segments de ligne, vérification des doublons et
correction des croisements entre les polygones. L'aire de chaque polygone a été déterminée
afin que l'utilisateur puisse choisir de ne pas dessiner d'entités plus petites qu'une zone minimale (voir -A); une
peut également limiter le niveau hiérarchique le plus élevé de polygones à inclure (4 est le
maximum). Les 4 bases de données à basse résolution ont été dérivées de la base de données à pleine résolution
en utilisant l'algorithme de simplification de ligne Douglas-Peucker. La classification des rivières et
les frontières suivent celle du WDBII. Voir le livre de recettes GMT et l'annexe K de référence technique
pour plus de détails
Utilisez gmtselectgmt en ligne à l'aide des services onworks.net