pygmt.Figure.velo¶
Figure.velo(data=None, *, vector=None, frame=None, cmap=None,
rescale=None, uncertaintyfill=None, fill=None, scale=None,
shading=None, projection=None, line=None, no_clip=None,
region=None, spec=None, verbose=None, pen=None, zvalue=None,
panel=None, nodata=None, find=None, header=None, incols=None,
perspective=None, transparency=None, **kwargs)
Tracez les vecteurs de vitesse, les croix, les barres d’anisotropie et les coins.
Lit les valeurs de données à partir de fichiers numpy.ndarray ou pandas.DataFrame trace le symbole de géodésie sélectionné sur une carte. Vous pouvez choisir parmi les vecteurs de vitesse et leurs incertitudes, les coins de rotation et leurs incertitudes, les barres d’anisotropie ou les croix de déformation. Les remplissages de symboles ou leurs contours peuvent être colorés en fonction de paramètres constants ou via des tables de correspondance de couleurs.
Doit fournir data et spec.
Paramètre¶
data ( str ou list ou numpy.ndarray ou pandas.DataFrame ou xarray.Dataset ou geopandas.GeoDataFrame ) – Transmettez soit un nom de fichier à une table de données ASCII, une liste Python, un 2-D numpy.ndarray, un pandas.DataFrame, un xarray.Dataset composé de Variables de données 1D xarray.DataArray ou a geopandas.GeoDataFrame contenant les données tabulaires.
spec ( str ) – Sélectionnez la signification des colonnes dans le fichier de données et la figure à tracer. Dans tous les cas, les échelles sont en unités de données par unité de longueur et les tailles sont en unités de longueur (l’unité de longueur par défaut est contrôlée par PROJ_LENGTH_UNIT sauf si c , i ou p est ajouté).
- e [ velscale /] confiance [ +f police ] Ellipses de vitesse en convention (N,E). Le velscale définit l’échelle des flèches de vélocité. Si velscale n’est pas donné, nous le lisons à partir du fichier de données en tant que colonne supplémentaire. La confiance définit la limite de confiance bidimensionnelle pour l’ellipse, par exemple 0,95 pour une ellipse de confiance à 95 %. Utilisez +f pour définir la police et la taille du texte [la valeur par défaut est 9p, Helvetica, black] ; donner +f 0 pour désactiver l’étiquetage. La flèche sera dessinée avec les attributs de stylo spécifiés par le pen paramètre et la tête de flèche peut être colorée via fill. L’ellipse sera remplie avec la couleur ou la nuance spécifiée par le uncertaintyfill paramètre [La valeur par défaut est transparente], et son contour sera dessiné si line est sélectionné à l’aide du stylo sélectionné (par pen s’il n’est pas donné par line). Les paramètres doivent se trouver dans les colonnes suivantes :
1 , 2 : longitude, latitude de la station
3 , 4 : vitesse vers l’est, vers le nord
5 , 6 : incertitude des vitesses vers l’est, vers le nord (1-sigma)
7 : corrélation entre les composantes est et nord
Texte de fin : nom de la station (optionnel)
n [ échelle à barres ] Barres d’anisotropie. barscale définit l’échelle des barres. Si barscale n’est pas donné, nous le lisons à partir du fichier de données comme une colonne supplémentaire. Les paramètres doivent se trouver dans les colonnes suivantes :
1 , 2 : longitude, latitude de la station
3 , 4 : composantes vers l’est et vers le nord du vecteur d’anisotropie
- r [ velscale /] confiance [ +f police ] Ellipses de vitesse en convention de rotation. Le velscale définit l’échelle des flèches de vélocité. Si velscale n’est pas donné, nous le lisons à partir du fichier de données en tant que colonne supplémentaire. La confiance définit la limite de confiance bidimensionnelle pour l’ellipse, par exemple 0,95 pour une ellipse de confiance à 95 %. Utilisez +f pour définir la police et la taille du texte [la valeur par défaut est 9p, Helvetica, black] ; donner +f 0 pour désactiver l’étiquetage. La flèche sera dessinée avec les attributs de stylo spécifiés par le pen paramètre et la tête de flèche peut être colorée via fill. L’ellipse sera remplie avec la couleur ou la nuance spécifiée par le uncertaintyfill paramètre [La valeur par défaut est transparente], et son contour sera dessiné si line est sélectionné à l’aide du stylo sélectionné (par pen s’il n’est pas donné par line). Les paramètres doivent se trouver dans les colonnes suivantes :
1 , 2 : longitude, latitude de la station
3 , 4 : vitesse vers l’est, vers le nord
5 , 6 : axes semi-majeur, semi-mineur
7 : angle anti-horaire, en degrés, de l’axe horizontal au grand axe de l’ellipse.
Texte de fin : nom de la station (optionnel)
w [ échelle de coin /] coinmag
Coins de rotation. L” échelle de coin définit la taille des coins. Si wedgescale n’est pas donné, nous le lisons à partir du fichier de données comme une colonne supplémentaire. Les valeurs de rotation sont multipliées par wedgemag avant le tracé. Par exemple, régler wedgemag sur 1.e7 fonctionne bien pour des rotations de l’ordre de 100 nanoradians/an. Utilisez fill pour définir la couleur ou l’ombre de remplissage du biseau et uncertaintyfill pour définir la couleur ou l’ombre de l’incertitude. Les paramètres doivent se trouver dans les colonnes suivantes :
1 , 2 : longitude, latitude de la station
3 : rotation en radians
4 : incertitude de rotation en radians
- x [ échelle_croisée ] La souche traverse. Le cross_scale définit la taille de la croix. Si cross_scale n’est pas donné, nous le lisons à partir du fichier de données en tant que colonne supplémentaire. Les paramètres doivent se trouver dans les colonnes suivantes :
1 , 2 : longitude, latitude de la station
3 : eps1, la valeur propre la plus extensionnelle du tenseur de déformation, avec extension prise positive.
4 : eps2, la valeur propre la plus compressive du tenseur de déformation, avec extension prise positive.
5 : azimut de eps2 en degrés CW depuis le Nord.
projection ( str ) – Sélectionnez la projection cartographique .
region ( chaîne ou liste ) – xmin/xmax/ymin/ymax [ +r ][ +u unité ]. Spécifiez la région d’intérêt.
vector ( bool ou str ) – Modifie les paramètres vectoriels. Pour les têtes vectorielles, ajoutez la taille de la tête vectorielle [la valeur par défaut est 9p].
frame ( bool ou str ou list ) – Définit le cadre de délimitation de la carte et les attributs des axes .
cmap ( str ) – Nom de fichier d’un fichier CPT ou d’une série de couleurs séparées par des virgules (par exemple, color1 , color2 , color3 ) pour créer automatiquement un CPT linéaire continu à partir de ces couleurs.
rescale ( str ) – Peut être utilisé pour redimensionner les incertitudes des vitesses ( spec= »e » et spec= »r ») et des rotations ( spec= »w »). Peut être combiné avec la confidence variable.
uncertaintyfill ( str ) – Définissez la couleur ou le motif pour remplir les coins d’incertitude ( spec= »w ») ou les ellipses d’erreur de vitesse ( spec= »e » ou spec= »r »). Si uncertaintyfill n’est pas spécifié, les régions d’incertitude seront transparentes. Remarque : L’utilisation cmap de et zvalue= »+e » mettra à jour la couleur de remplissage de l’incertitude en fonction de la mesure sélectionnée dans zvalue [La valeur par défaut est l’erreur d’amplitude].
fill ( str ) – Définit la couleur ou le motif des symboles de remplissage [par défaut, pas de remplissage]. Remarque : L’utilisation cmap de (et éventuellement de zvalue) mettra à jour la couleur de remplissage du symbole en fonction de la mesure sélectionnée dans zvalue [La valeur par défaut est la magnitude].
scale ( float ou bool ) – [ échelle ]. Mettre à l’échelle les tailles des symboles et les largeurs de plumes pour chaque enregistrement en utilisant l” échelle lue à partir de l’ensemble de données, indiquée dans la première colonne après les colonnes (facultatives) z et taille [la valeur par défaut est aucune mise à l’échelle]. La taille du symbole est fournie par spec ou via la colonne de taille d’entrée . Vous pouvez également ajouter une échelle constante qui doit être utilisée au lieu de lire une colonne d’échelle.
line ( str ) – [ stylo [ +c [ f | l ]]]. Tracer des lignes. Les ellipses et les biseaux rotatifs auront leurs contours dessinés à l’aide du stylet actuel (voir pen). Vous pouvez également ajouter un stylo séparé à utiliser pour les contours d’erreur. Si le modificateur +cl est ajouté, la couleur du stylo est mise à jour à partir du CPT (voir cmap). Si à la place le modificateur +cf est ajouté, la couleur du fichier cpt est appliquée uniquement au remplissage d’erreur [par défaut]. Utilisez simplement +c pour définir à la fois la plume et la couleur de remplissage.
no_clip ( bool ) – Ne saute pas les symboles qui tombent en dehors des limites du cadre [la valeur par défaut est False, c’est-à-dire, les symboles de tracé à l’intérieur des limites du cadre uniquement].
- verbose ( bool ou str ) – Sélectionnez le niveau de verbosité [la valeur par défaut est w ], qui module les messages écrits sur stderr. Choisissez parmi 7 niveaux de verbosité :
q - Des messages d’erreur silencieux, même fatals, ne sont pas générés
e - Messages d’erreur uniquement
w - Avertissements [par défaut]
t - Timings (temps d’exécution des rapports pour les algorithmes chronophages)
i - Messages d’information (comme verbose=True)
c - Avertissements de compatibilité
d - Messages de débogage
pen ( str ) – [ pen ][ +c [ f | l ]]. Définissez les attributs du stylet pour les flèches de vitesse, la circonférence de l’ellipse et les bords du plan de faille [la valeur par défaut est « 0.25p,black,solid »]. Si le modificateur +cl est ajouté, la couleur du stylo est mise à jour à partir du CPT (voir cmap). Si à la place le modificateur +cf est ajouté, la couleur du fichier cpt est appliquée uniquement au remplissage du symbole [par défaut]. Utilisez simplement +c pour définir à la fois la plume et la couleur de remplissage.
zvalue ( str ) – [ m | e | n | u][ +e ]. Sélectionnez la quantité qui sera utilisée avec le CPT donné via cmap pour définir la couleur de remplissage. Choisissez parmi la m agnitude (magnitude du vecteur ou magnitude de rotation), la vitesse est -ouest, la vitesse n ord -sud ou la colonne de données fournie par l’utilisateur (fournie après les colonnes requises). Pour utiliser à la place les estimations d’erreur correspondantes (c’est-à-dire vecteur ou incertitude de rotation) pour rechercher la couleur et peindre l’ellipse ou le coin d’erreur à la place, ajoutez +e .
panel ( bool ou int ou list ) – [ row,col | indice ]. Sélectionnez un panneau de sous-parcelle spécifique. Autorisé uniquement en mode sous-parcelle. Utilisez panel=True pour passer au panneau suivant dans l’ordre sélectionné. Au lieu de row,col, vous pouvez également donner un index de valeur scalaire qui dépend de l’ordre que vous autolabel avez défini lors de la définition de la sous-parcelle. Remarque : row , col et index commencent tous à 0.
nodata ( str ) – i | o aucune donnée . 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. Ajoutez i à la valeur nodata pour les colonnes d’entrée uniquement. Ajoutez o à la valeur nodata pour les colonnes de sortie uniquement.
find ( chaîne ) – [ ~ ] « motif » | [ ~ ]/ expression régulière /[ i ]. Ne transmettez que les enregistrements qui correspondent au modèle ou aux expressions régulières donnés [par défaut, traite tous les enregistrements]. Ajoutez ~ 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.
- header (str) – [i|o][n][+c][+d][+msegheader][+rremark][+ttitle]. Spécifiez que les fichiers d’entrée et/ou de sortie ont n enregistrements d’en-tête [la valeur par défaut est 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 des noms de colonne à 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, pas d’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, pas de titre]. La chaîne de titre peut contenir n pour indiquer des sauts de ligne.
- incols ( str ou tableau 1-D ) – 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, 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 2ème colonne suivie de la 1ère colonne).
- Pour str: spécifiez des colonnes individuelles ou des plages de colonnes au format start [: inc ]: stop , où inc prend par défaut la valeur 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 saisir les trois premières colonnes suivies du 5e colonne transformée en log). Pour lire à partir d’une colonne donnée jusqu’à la fin de l’enregistrement, laissez off stop lors de la spécification de la plage de colonnes. Pour lire le texte de fin, ajoutez la colonne t . Ajoutez le numéro de mot à t pour ingérer un seul mot du texte de fin. Au lieu de spécifier des colonnes, utilisez incols= »n » pour lire simplement l’entrée numérique et ignorer le texte de fin. Ajoutez éventuellement l’un des modificateurs suivants à n’importe quelle colonne ou plage de colonnes pour transformer les colonnes d’entrée :
+l pour prendre le log10 des valeurs d’entrée.
+d pour diviser les valeurs d’entrée par le facteur diviseur [la valeur par défaut est 1].
+s pour multiplier les valeurs d’entrée par l” échelle des facteurs [la valeur par défaut est 1].
+o pour ajouter le décalage donné aux valeurs d’entrée [la valeur par défaut est 0].
perspective ( liste ou chaîne ) – [ x | y | z ] azim [/ elev [/ zlevel ]][ +w lon0 / lat0 [/ z0 ]][ +v x0 / y0 ]. Sélectionnez la vue en perspective et définissez l’azimut et l’angle d’élévation du point de vue [la valeur par défaut est [180,90]].
transparence ( int ou float ) – Définissez le niveau de transparence, dans la plage [0-100] pour cent [la valeur par défaut est 0, c’est-à-dire opaque]. Uniquement visible lorsque la sortie au format PDF ou raster est sélectionnée. Seule la sélection du format PNG ajoute une couche de transparence dans l’image (pour un traitement ultérieur).
Exemple¶
import pygmt
# Load sample compilation of fracture lengths and azimuth as
# hypothetically digitized from geological maps
data = pygmt.datasets.load_sample_data(name="fractures")
fig = pygmt.Figure()
fig.rose(
# use columns of the sample dataset as input for the length and azimuth
# parameters
length=data.length,
azimuth=data.azimuth,
# specify the "region" of interest in the (r,azimuth) space
# [r0, r1, az0, az1], here, r0 is 0 and r1 is 1, for azimuth, az0 is 0 and
# az1 is 360 which means we plot a full circle between 0 and 360 degrees
region=[0, 1, 0, 360],
# set the diameter of the rose diagram to 7.5 cm
diameter="7.5c",
# define the sector width in degrees, we append +r here to draw a rose
# diagram instead of a sector diagram
sector="10+r",
# normalize bin counts by the largest value so all bin counts range from
# 0 to 1
norm=True,
# use red3 as color fill for the sectors
fill="red3",
# define the frame with ticks and gridlines every 0.2
# length unit in radial direction and every 30 degrees
# in azimuthal direction, set background color to
# lightgray
frame=["x0.2g0.2", "y30g30", "+glightgray"],
# use a pen size of 1p to draw the outlines
pen="1p",
)
fig.show()