xyz2grd

Ici vous aurez des explications sur la commande et ses différents paramètres.

pygmt.xyz2grd(data=None, x=None, y=None, z=None, *, duplicate=None,
outgrid=None, spacing=None, projection=None, region=None,
verbose=None, convention=None, binary=None, nodata=None,
find=None, coltypes=None, header=None, incols=None,
registration=None, wrap=None, **kwargs)

Utilité :

Crée un fichier de grille à partir de données de table.

Lit une ou plusieurs tables avec des colonnes x, y, z et crée un fichier de grille binaire. xyz2grd signalera si certains des nœuds ne sont pas remplis avec des données. Ces nœuds non contraints sont définis sur une valeur spécifiée par l’utilisateur [Par défaut, NaN]. Les nœuds avec plusieurs valeurs seront définis sur la valeur moyenne.

Paramètres :

data (str ou numpy.ndarray ou pandas.DataFrame ou xarray.Dataset ou geopandas.GeoDataFrame) : Fournissez les valeurs (x, y, z) ou (longitude, latitude, élévation) en spécifiant le nom d’un fichier de tableau ASCII, un tableau numpy 2D, un DataFrame pandas, un Dataset xarray composé de variables de données xarray.DataArray 1D, ou un GeoDataFrame geopandas contenant les données tabulaires.

x/y/z (tableaux 1D) : Tableaux des coordonnées x et y et des valeurs z des points de données.

outgrid (str ou None) - Le nom du fichier de sortie netCDF avec l’extension .nc pour stocker la grille.

duplicate (str) - [d|f|l|m|n|r|S|s|u|z]. Par défaut, nous calculerons les valeurs moyennes si plusieurs entrées tombent sur le même nœud. Utilisez -A pour modifier ce comportement, sauf si -Z est spécifié. Ajoutez f ou s pour simplement conserver le premier ou dernier point de données attribué à chaque nœud. Ajoutez l ou u ou d pour trouver la valeur la plus basse (minimum) ou supérieure (maximum) ou la différence entre la valeur maximale et minimale de chaque nœud, respectivement. Ajoutez m ou r ou S pour calculer la valeur moyenne ou la valeur RMS ou l’écart-type de chaque nœud, respectivement. Ajoutez n pour simplement compter le nombre de points de données attribués à chaque nœud (cela ne nécessite que deux colonnes d’entrée x et y car z n’est pas consulté). Ajoutez z pour faire la somme des valeurs multiples qui appartiennent au même nœud.

spacing (str) - x_inc[+e|n][/y_inc[+e|n]]. x_inc [et éventuellement y_inc] est l’espacement de la grille.

Coordonnées géographiques (degrés) : Ajoutez éventuellement une unité d’incrément. Choisissez parmi m pour indiquer les minutes d’arc ou s pour indiquer les secondes d’arc. Si l’une des unités e, f, k, M, n ou u est ajoutée à la place, l’incrément est supposé être donné en mètre, pied, kilomètre, mile, mille marin ou pied d’enquête américain, respectivement, et sera converti en degrés de longitude équivalents à la latitude moyenne de la région (la conversion dépend de PROJ_ELLIPSOID). Si y_inc est donné mais défini sur 0, il sera réinitialisé à x_inc ; sinon, il sera converti en degrés de latitude.

Toutes les coordonnées : Si +e est ajouté, alors le max correspondant x (est) ou y (nord) peut être légèrement ajusté pour correspondre exactement à l’incrément donné [par défaut, l’incrément peut être légèrement ajusté pour correspondre au domaine donné]. Enfin, au lieu de donner un incrément, vous pouvez spécifier le nombre de nœuds souhaités en ajoutant +n à l’argument entier fourni ; l’incrément est ensuite recalculé à partir du nombre de nœuds, de l’enregistrement et du domaine. La valeur d’incrément résultante dépend de la sélection d’une grille enregistrée par lignes de grille ou enregistrée par pixels ; voir les détails dans les formats de fichier GMT.

Note : Si region=grdfile est utilisé, alors l’espacement de la grille et l’enregistrement ont déjà été initialisés ; utilisez spacing et registration pour remplacer ces valeurs.

projection (str) – codeproj[paramsproj/]largeur. Sélectionnez la projection cartographique.

region (str ou liste) - xmin/xmax/ymin/ymax[+r][+uunit]. Spécifie la région d’intérêt.

verbose (bool or str) - Sélectionnez le niveau de verbosité [Par défaut, w], qui module les messages écrits dans stderr. Choisissez parmi 7 niveaux de verbosité :

q - Silencieux, aucun message d’erreur fatale n’est produit

e - Messages d’erreur uniquement

w - Avertissements [Par défaut]

t - Durées (rapporte les temps d’exécution pour les algorithmes gourmands en temps)

i - Messages informatifs (identique à verbose=True)

c - Avertissements de compatibilité

d - Messages de débogage

convention (str) - [drapeaux]. Lit une table ASCII [ou binaire] à 1 colonne. Cela suppose que tous les nœuds sont présents et triés selon la convention d’ordonnancement spécifiée dans les drapeaux. Si les données entrantes représentent des lignes, commencez les drapeaux par T(op) si la première ligne est y = ymax ou B(ottom) si la première ligne est y = ymin. Ensuite, ajoutez L ou R pour indiquer que le premier élément est à gauche ou à droite de la ligne. De même pour les formats de colonnes : commencez par L ou R pour positionner la première colonne, puis ajoutez T ou B pour positionner le premier élément dans une ligne. Note : Ces deux indicateurs de ligne/colonne ne sont nécessaires que pour les grilles ; pour les autres tables, ils ne s’appliquent pas. Pour les grilles enregistrées par lignes de grille : Si les données sont périodiques en x mais que les données entrantes ne contiennent pas la colonne (redondante) à x = xmax, ajoutez x. Pour les données périodiques en y sans ligne redondante à y = ymax, ajoutez y. Ajoutez sn pour ignorer les n premiers octets (probablement un en-tête). Si l’ordre des octets ou les mots doivent être inversés, ajoutez w. Sélectionnez un des plusieurs types de données (tous binaires sauf a) :

A Représentation ASCII de une ou plusieurs valeurs à virgule flottante par enregistrement

a Représentation ASCII d’un seul élément par enregistrement

c int8_t, caractère signé sur 1 octet

u uint8_t, caractère non signé sur 1 octet

h int16_t, entier signé sur 2 octets

H uint16_t, entier non signé sur 2 octets

i int32_t, entier signé sur 4 octets

I uint32_t, entier non signé sur 4 octets

l int64_t, long (entier sur 8 octets)

L uint64_t, unsigned long (entier non signé sur 8 octets)

f Flottant simple précision sur 4 octets

d Flottant double précision sur 8 octets

[Le format par défaut est l’orientation des lignes de balayage des nombres ASCII : La]. La différence entre A et a est que ce dernier peut décoder à la fois les formats dateTclock et ddd:mm:ss[.xx], mais il s’attend à ce que chaque enregistrement d’entrée ait une seule valeur, tandis que le premier peut traiter plusieurs valeurs par enregistrement mais ne peut analyser que des valeurs à virgule flottante régulières. Traduire les valeurs z entrantes via le paramètre incols.

binary (bool or str) - i|o[ncols][type][w][+l|b]. Sélectionnez une entrée binaire native (en utilisant binary= »i ») ou une sortie binaire (en utilisant binary= »o »), où ncols est le nombre de colonnes de données du type spécifié, qui peut être l’un des suivants :

c - int8_t (1 octet signé)

u - uint8_t (1 octet non signé)

h - int16_t (2 octets signés)

H - uint16_t (2 octets non signés)

i - int32_t (4 octets signés)

I - uint32_t (4 octets non signés)

l - int64_t (8 octets signés)

L - uint64_t (8 octets non signés)

f - flottant sur 4 octets (simple précision)

d - flottant sur 8 octets (double précision)

x - utiliser pour sauter ncols n’importe où dans l’enregistrement

Pour les enregistrements avec des types mixtes, ajoutez des combinaisons supplémentaires de ncols type séparées par des virgules (sans espace). Les modificateurs suivants sont pris en charge :

w après un élément pour forcer l’inversion des octets.

+l|b pour indiquer que le fichier de données complet doit être lu en tant que petit ou grand boutiste, respectivement.

La documentation complète se trouve à l’adresse https://docs.generic-mapping-tools.org/latest/gmt.html#bi-full.

nodata (str) - i|onodata. Remplacez les valeurs spécifiques par NaN (pour les données tabulaires). Par exemple, nodata= »-9999 » remplacera toutes les valeurs égales à -9999 par NaN lors de l’entrée et toutes les valeurs NaN par -9999 lors de la sortie. Préfixez i à la valeur nodata pour les colonnes d’entrée uniquement. Préfixez o à la valeur nodata pour les colonnes de sortie uniquement.

find (str) - [~]“pattern” | [~]/regexp/[i]. Ne transmettre que les enregistrements correspondant au motif ou à l’expression régulière donnée [Par défaut, tous les enregistrements sont traités]. Préfixez ~ au motif ou à l’expression régulière pour ne transmettre que les expressions de données qui ne correspondent pas au motif. Ajoutez i pour une correspondance insensible à la casse. Cela ne s’applique pas aux en-têtes ou aux en-têtes de segment.

coltypes (str) - [i|o]colinfo. Spécifiez les types de données des colonnes d’entrée et/ou de sortie (données temporelles ou géographiques). La documentation complète se trouve à l’adresse https://docs.generic-mapping-tools.org/latest/gmt.html#f-full.

header (str) - [i|o][n][+c][+d][+msegheader][+rremark][+ttitle]. Spécifiez que le(s) fichier(s) d’entrée et/ou de sortie ont n enregistrements d’en-tête [Par défaut, 0]. Préfixez i si seule l’entrée principale doit avoir des enregistrements d’en-tête. Préfixez o pour contrôler l’écriture des enregistrements d’en-tête, avec les modificateurs suivants pris en charge :

+d pour supprimer les enregistrements d’en-tête existants.

+c pour ajouter un commentaire d’en-tête avec les noms des colonnes à la sortie [Par défaut, pas de noms de colonnes].

+m pour ajouter un en-tête de segment segheader à la sortie après le bloc d’en-tête [Par défaut, pas d’en-tête de segment].

+r pour ajouter un commentaire remarque à la sortie [Par défaut, pas de commentaire]. La chaîne de remarque peut contenir n pour indiquer des sauts de ligne.

+t pour ajouter un commentaire de titre à la sortie [Par défaut, pas de titre]. La chaîne de titre peut contenir n pour indiquer des sauts de ligne.

Les lignes vides et les lignes commençant par # sont toujours ignorées.

incols (str ou tableau 1D) - Spécifiez les colonnes de données pour l’entrée principale dans un ordre arbitraire. Les colonnes peuvent être répétées et les colonnes non répertoriées seront ignorées [Par défaut, lit toutes les colonnes dans l’ordre, en commençant par la première (c’est-à-dire, colonne 0)].

Pour un tableau 1D : spécifiez les colonnes individuelles dans l’ordre d’entrée (par exemple, incols=[1,0] pour la 2e colonne suivie de la 1ère colonne).

Pour une chaîne : spécifiez les colonnes individuelles ou les plages de colonnes selon le format start[:inc]:stop, où inc est par défaut à 1 s’il n’est pas spécifié, avec les colonnes et/ou les plages de colonnes séparées par des virgules (par exemple, incols= »0:2,4+l » pour lire les trois premières colonnes suivies de la colonne 5e transformée en logarithme). Pour lire à partir d’une colonne donnée jusqu’à la fin de l’enregistrement, laissez stop lorsque vous spécifiez la plage de colonnes. Pour lire du texte de fin, ajoutez la colonne t. Ajoutez le mot number à t pour n’ingérer qu’un seul mot du texte de fin. Au lieu de spécifier des colonnes, utilisez incols= »n » pour simplement lire une entrée numérique et ignorer le texte de fin. En option, ajoutez l’un des modificateurs suivants à n’importe quelle colonne ou plage de colonnes pour transformer les colonnes d’entrée :

+l pour prendre le logarithme décimal des valeurs d’entrée.

+d pour diviser les valeurs d’entrée par le diviseur [Par défaut, 1].

+s pour multiplier les valeurs d’entrée par l’échelle [Par défaut, 1].

+o pour ajouter le décalage donné aux valeurs d’entrée [Par défaut, 0].

registration (str) - g|p. Force l’enregistrement des nœuds sur les lignes de quadrillage (g) ou les pixels (p) [Par défaut, g (lignes de quadrillage)].

wrap (str) - y|a|w|d|h|m|s|cperiod[/phase][+ccol]. Convertir la coordonnée x d’entrée en une coordonnée cyclique, ou une colonne différente si elle est sélectionnée via +ccol. Les transformations de coordonnées cycliques suivantes sont prises en charge :

y - cycle annuel (normalisé)

a - cycle annuel (mensuel)

w - cycle hebdomadaire (jour)

d - cycle quotidien (heure)

h - cycle horaire (minute)

m - cycle minute (seconde)

s - cycle seconde (seconde)

c - cycle personnalisé (normalisé)