binstats¶
Ici vous aurez des explications sur la commande et ses différents paramètres.
pygmt.binstats(data, \*, statistic=None, empty=None, outgrid=None,
spacing=None, normalize=None, region=None, search_radius=None,
rttyytrtghgffverbose=None, weight=None, aspatial=None, binary=None,
header=None, incols=None, registration=None, \**kwargs)
Utilité :¶
Regroupe les données spatiales en bacs et détermine des statistiques par bac.
Lit des points situés de manière arbitraire (x, y[, z][, w]) à partir des données, avec 2 à 4 colonnes, et pour chaque nœud dans la disposition de grille spécifiée, détermine quels points se trouvent dans le rayon donné. Ces points sont ensuite utilisés dans le calcul de la statistique spécifiée. Les résultats peuvent être présentés tels quels ou normalisés par la superficie du cercle pour obtenir éventuellement des estimations de densité.
Paramètres :¶
data (str or numpy.ndarray or pandas.DataFrame or xarray.Dataset or geopandas.GeoDataFrame) - Nom du fichier d’une table de données ASCII ou d’un tableau numpy.ndarray à 2 dimensions, d’un pandas.DataFrame, d’un xarray.Dataset composé de variables de données xarray.DataArray à 1 dimension, ou d’un geopandas.GeoDataFrame contenant les données tabulaires.
outgrid (str ou None) - Nom du fichier de sortie netCDF avec l’extension .nc pour stocker la grille.
statistic (str) -a|d|g|i|l|L|m|n|o|p|q[quant]|r|s|u|U|z. Choisissez la statistique qui sera calculée par nœud en fonction des points situés à une distance de rayon du nœud. Choisissez parmi :
a pour la moyenne (average)
d pour la déviation absolue médiane (MAD)
g pour l’étendue complète (max-min)
i pour la plage interquartile 25-75%
l pour le minimum (low)
L pour le minimum des valeurs positives uniquement
m pour la médiane
n pour le nombre de valeurs
o pour l’échelle LMS
p pour le mode (maximum de vraisemblance)
q pour un quantile sélectionné (ajoutez le quantile souhaité dans la plage de 0 à 100% [50])
r pour la racine carrée de la moyenne des carrés (r.m.s.)
s pour l’écart type
u pour le maximum (upper)
U pour le maximum des valeurs négatives uniquement
z pour la somme
empty (float or int) - Définit la valeur assignée aux nœuds vides [Par défaut, NaN].
normalize (bool) - Normalise les valeurs de grille résultantes par la surface représentée par le rayon de recherche [pas de normalisation].
search_radius (float or str) - Définit le rayon de recherche qui détermine quels points de données sont considérés comme proches d’un nœud. Ajoutez l’unité de distance. Non compatible avec le découpage en tuiles (tiling).
weight (str) - Les données d’entrée contiennent une colonne supplémentaire contenant le poids du point d’observation. Si des poids sont donnés, des quantités statistiques pondérées seront calculées, tandis que le nombre total sera la somme des poids au lieu du nombre de points. Si les poids sont en fait des incertitudes (une sigma), ajoutez +s et weight = 1/sigma.
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) : En option, ajoutez 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, l’incrément est supposé être donné en mètre, pied, kilomètre, mille, mille marin ou pied d’étalonnage 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 à 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 la valeur maximale correspondante de x (est) ou y (nord) peut être légèrement ajustée 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 spécifier un incrément, vous pouvez spécifier le nombre de nœuds souhaité en ajoutant +n à l’argument entier fourni ; l’incrément est alors recalculé à partir du nombre de nœuds, de l’enregistrement et du domaine. La valeur d’incrément résultante dépend de la grille enregistrée sur les lignes de grille ou enregistrée sur les pixels ; consultez la documentation sur les formats de fichiers GMT pour plus de détails.
Remarque : Si region=grdfile est utilisé, l’espacement de la grille et l’enregistrement ont déjà été initialisés ; utilisez spacing et registration pour remplacer ces valeurs.
region (str ou liste) - xmin/xmax/ymin/ymax[+r][+uunit]. Spécifiez la région d’intérêt.
verbose (bool ou str) - Sélectionnez le niveau de verbosité [par défaut w], qui module les messages écrits sur stderr. Choisissez parmi 7 niveaux de verbosité :
q - Silencieux, même les messages d’erreur fatale ne sont pas produits
e - Messages d’erreur uniquement
w - Avertissements [par défaut]
t - Chronométrage (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
aspatial (bool ou str) - [col=]name[,…]. Contrôle la façon dont les données aspatiales sont gérées lors de l’entrée et de la sortie. La documentation complète se trouve à l’adresse https://docs.generic-mapping-tools.org/latest/gmt.html#aspatial-full.
binaire (bool ou str) - i|o[ncols][type][w][+l|b]. Sélectionnez l’entrée binaire native (en utilisant binary= »i ») ou la sortie (en utilisant binary= »o »), où ncols est le nombre de colonnes de données de type, qui doit être l’un des suivants :
c - int8_t (char signé sur 1 octet)
u - uint8_t (char 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 (entier signé sur 8 octets)
L - uint64_t (entier non signé sur 8 octets)
f - flottant simple précision sur 4 octets
d - flottant double précision sur 8 octets
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 (pas d’espace). Les modificateurs suivants sont pris en charge :
w après n’importe quel élément pour forcer le renversement des octets.
+l|b pour indiquer que le fichier de données entier doit être lu en tant que petit ou gros-boutien, respectivement.
La documentation complète se trouve à l’adresse https://docs.generic-mapping-tools.org/latest/gmt.html#bi-full.
header (str) - [i|o][n][+c][+d][+msegheader][+rremark][+ttitle]. Spécifie que le(s) fichier(s) d’entrée et/ou de sortie a(n) n en-têtes [par défaut, n=0]. Préfixez i si seule l’entrée principale doit avoir des en-têtes. Préfixez o pour contrôler l’écriture des en-têtes, avec les modificateurs suivants pris en charge :
+d pour supprimer les en-têtes existants.
+c pour ajouter un commentaire d’en-tête avec les noms de colonnes à la sortie [par défaut, aucun nom de colonne].
+m pour ajouter un en-tête de segment segheader à la sortie après le bloc d’en-tête [par défaut, aucun en-tête de segment].
+r pour ajouter un commentaire de remarque à la sortie [par défaut, aucun 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, aucun 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.
ncols (str ou tableau 1-D) -
Spécifie 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, toutes les colonnes sont lues dans l’ordre, en commençant par la première (c’est-à-dire la colonne 0)].
Pour un tableau 1-D : 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 de caractères : spécifiez les colonnes individuelles ou les plages de colonnes au format start[:inc]:stop, où inc est par défaut à 1 s’il n’est pas spécifié, avec des colonnes et/ou des 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 5ème avec transformation logarithmique). 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 supplémentaire, ajoutez la colonne t. Ajoutez le mot number à t pour n’ingérer qu’un seul mot du texte supplémentaire. Au lieu de spécifier des colonnes, utilisez incols= »n » pour simplement lire une entrée numérique et ignorer le texte supplémentaire. Facultativement, 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 facteur divisor [par défaut, 1].
+s pour multiplier les valeurs d’entrée par le facteur scale [par défaut, 1].
+o pour ajouter l’offset donné aux valeurs d’entrée [par défaut, 0].
registration (str) - g|p. Force l’enregistrement des nœuds de grille (g) ou des pixels (p) [par défaut, g (grille)].