pygmt.grdtrack

pygmt.grdtrack(grid, points=None, newcolname=None, outfile=None, \*,
resample=None, crossprofile=None, dfile=None, profile=None, critical=None,
region=None, no_skip=None, stack=None, radius=None, verbose=None, z_only=None,
aspatial=None, binary=None, nodata=None, find=None, coltypes=None, gap=None,
header=None, incols=None, distcalc=None, interpolation=None, outcols=None,
skiprows=None, wrap=None, \**kwargs)

Exemples de grilles à des emplacements (x,y) spécifiés.

Lit un ou plusieurs fichiers de grille et une table (à partir d’une entrée de fichier ou de tableau ; mais voir profile pour les exceptions) avec des positions (x, y) [ou (lon, lat)] dans les deux premières colonnes (plus de colonnes peuvent être présentes). Il interpole la grille (s) aux positions dans la table et écrit la table avec les valeurs interpolées ajoutées comme (une ou plusieurs) nouvelles colonnes. Alternativement ( crossprofile ), l’entrée est considérée comme des segments de ligne et nous créons des profils croisés orthogonaux à chaque point de données ou avec une séparation équidistante et échantillonnons la ou les grilles le long de ces profils. Une interpolation bicubique [par défaut], bilinéaire, B-spline ou au plus proche voisin est utilisée, nécessitant des conditions aux limites aux limites de la région (voir interpolation; La valeur par défaut utilise des conditions « naturelles » (la dérivée partielle seconde normale au bord est nulle) à moins que la grille ne soit automatiquement reconnue comme périodique.)

Paramètre

  • grid ( str ou xarray.DataArray ) – Le nom de fichier de la grille d’entrée ou de la grille chargée en tant que xarray.DataArray. C’est le seul paramètre obligatoire.

  • points ( str ou numpy.ndarray ou pandas.DataFrame ou xarray.Dataset ou geopandas.GeoDataFrame ) – Transmettez soit un nom de fichier à une table de données ASCII, un 2-D numpy.ndarray, un pandas.DataFrame, un xarray.Dataset composé de xarray.DataArray variables de données 1-D, ou a geopandas.GeoDataFrame contenant les données tabulaires.

  • newcolname ( str ) – Obligatoire si pointsest un pandas.DataFrame. Le nom de la nouvelle colonne dans la table de suivi pandas.DataFrame où les valeurs échantillonnées seront placées.

  • outfile ( str ) – Le nom de fichier pour le fichier ASCII de sortie.

  • resample ( str ) – f | p | m | r | R [ +l ] Pour le rééchantillonnage de piste (si crossprofileou profilesont définis), nous pouvons sélectionner la manière dont cela doit être effectué. Ajoutez f pour conserver les points d’origine, mais ajoutez des points intermédiaires si nécessaire [Par défaut], m comme f , mais suivez d’abord le méridien (le long de y) puis parallèle (le long de x), p comme f , mais suivez d’abord le parallèle (le long de y) puis le méridien (le long de x), r pour rééchantillonner à des emplacements équidistants ; les points d’entrée ne sont pas nécessairement inclus dans la sortie, et R commer , mais ajustez l’espacement donné pour s’adapter exactement à la longueur de la piste. Enfin, ajoutez +l si les distances géographiques doivent être mesurées le long des lignes de rhumb (loxodromes) au lieu des grands cercles. Ignoré sauf si crossprofile est utilisé.

  • corssprofile ( str ) – longueur / ds [ / espacement ][ +a | +v ][ l | r ]. Utilisez des segments de ligne d’entrée pour créer un ensemble de profils de croisement équidistants et (éventuellement) équidistants le long desquels nous échantillonnons la ou les grilles [par défaut, échantillonne simplement la ou les grilles aux emplacements d’entrée]. Spécifiez deux échelles de longueur qui contrôlent la façon dont l’échantillonnage est effectué : longueur définit la longueur totale de chaque profil transversal, tandis que ds est l’espacement d’échantillonnage le long de chaque profil transversal. Facultativement, ajouter / espacerpour un espacement équidistant entre les profils transversaux [La valeur par défaut érige les profils transversaux aux coordonnées d’entrée] ; voyez resample comment le rééchantillonnage de la piste d’entrée est contrôlé. Par défaut, tous les profils transversaux ont la même direction (de gauche à droite lorsque nous regardons dans la direction du segment de ligne d’entrée). Ajoutez +a pour alterner la direction des profils en travers, ou v pour appliquer une vue « ouest-est » ou « sud-nord ». Par défaut, tous les profils sont sortis. Choisissez de ne sortir que les moitiés gauche ou droite des profils en ajoutant +l ou +r , respectivement. Ajouter les unités appropriées à la longueur ; il définit l’unité utilisée pour ds [et l’espacement ] (voirUnités ). L’unité par défaut pour les grilles géographiques est le mètre tandis que les grilles cartésiennes impliquent l’unité utilisateur. Les colonnes de sortie seront lon , lat , dist , azimuth , z1 , z2 , …, zn (les zi sont les valeurs échantillonnées pour chacune des n grilles).

  • dfile ( str ) – De concert avec crossprofile, nous pouvons enregistrer les lignes d’origine (éventuellement rééchantillonnées) dans dfile [par défaut, enregistre uniquement les profils croisés]. Les colonnes seront lon , lat , dist , azimut , z1 , z2 , … (valeur échantillonnée pour chaque grille).

  • profile ( str ) – ligne [, ligne ,…][ +a az ][ +c ][ +d ][ +g ][ +i inc ][ +l longueur ][ +n np ][ +o az ][ +r rayon ]. Au lieu de lire les coordonnées des pistes d’entrée, spécifiez les profils via les coordonnées et les modificateurs. Le format de chaque ligne est start / stop , où start ou stop sont soit lon / lat (x / y pour les données cartésiennes) ou une clé XY à 2 caractères qui utilise le format de justification de style texte pour spécifier un point sur la carte comme [LCR][BMT]. Chaque ligne sera un segment séparé sauf si +c est utilisé, ce qui connectera les segments avec des joints partagés en un seul segment. En plus des coordonnées de ligne, vous pouvez utiliser Z-, Z+ pour désigner les emplacements globaux minimum et maximum dans la grille (uniquement disponible si une seule grille est donnée via outfile ). Vous pouvez ajouter +i inc pour définir l’intervalle d’échantillonnage ; s’il n’est pas donné, nous utilisons par défaut la moitié de l’intervalle de grille minimum. Pour une ligne le long de parallèles ou de méridiens, vous pouvez ajouter +g pour signaler des degrés de longitude ou de latitude au lieu de distances orthodromiques commençant à zéro. Au lieu de deux coordonnées, vous pouvez spécifier une origine et l’une de +a , +o ou +r . Le +a définit l’azimut d’un profil de longueur donnée à partir de l’origine donnée, tandis que +o centre le profil sur l’origine ; les deux nécessitent +l . Pour un échantillonnage circulaire, spécifiez +r pour définir un cercle de rayon donné centré sur l’origine ; cette option nécessite soit +n soit +i . Le modificateur +n np définit le nombre de points souhaité, tandis que +l longueurdonne la longueur totale du profil. Ajoutez +d pour sortir les distances le long de la trace après les coordonnées.

    Remarque : Aucun fichier de piste ne sera lu. Notez également qu’une seule unité de distance peut être choisie. Donner des unités différentes entraînera une erreur. Si aucune unité n’est spécifiée, nous utiliserons par défaut les distances orthodromiques en km (si géographiques). Si vous travaillez avec des données géographiques, vous pouvez utiliser distcalc pour contrôler le mode de calcul de distance [La valeur par défaut est Great Circle]. Remarque : Si crossprofile est défini et que l’espacement est donné, ce schéma d’échantillonnage annule tout modificateur défini dans profile.

  • critical ( str ) – [ +b ][ +n ][ +r ][ +z z0 ]. Trouvez les points critiques le long de chaque profil transversal en fonction de la distance le long de la trace. Nécessite crossprofile et une seule grille d’entrée ( z ). Nous examinons chaque profil croisé généré et rapportons ( dist , lonc , latc , distc , azimuthc , zc ) au pic central de la valeur z maximale , ( lonl , latl , distl ) et ( lonr , latr , distr ) au premier et au dernier point non NaN dont la valeur z dépasse z0 , respectivement, et la largeur basée sur les deux points extrêmes trouvés. Ici, dist est la distance le long de l’entrée d’origine points et les 12 autres colonnes de sortie sont fonction de cette distance. Lors de la recherche du pic central et des premières et dernières valeurs extrêmes qui dépassent le seuil, nous supposons que le profil est positif vers le haut. Si nous recherchons plutôt un creux, vous devez utiliser +n pour retourner temporairement le profil en positif. La valeur du seuil z0 est toujours donnée comme >= 0 ; utilisez +z pour le modifier [la valeur par défaut est 0]. Sinon, utilisez +b pour déterminer le point d’équilibre et l’écart type du profil ; il s’agit de la moyenne pondérée et de l’écart type pondéré des distances, z servant de poids. Enfin, utilisez +r pour obtenir la valeur efficace pondérée autour du centre de la voie transversale ( distc == 0). Remarque : Nous arrondissons les résultats exacts aux nœuds de distance les plus proches le long des profils transversaux. Nous écrivons 13 colonnes de sortie par piste : dist, lonc, latc, distc, azimuthc, zc, lonl, latl, distl, lonr, latr, distr, width .

  • région ( str ou liste ) – xmin/xmax/ymin/ymax [ +r ][ +u unité ]. Spécifiez la région d’intérêt.

  • no_skip ( bool ) – Ne saute pas les points qui se trouvent en dehors du domaine de la ou des grilles [par défaut, seuls les points de sortie dans le domaine de la grille].

  • stack ( str ou liste ) – méthode / modificateurs . En conjonction avec crossprofile, calculez un seul profil empilé à partir de tous les profils de chaque segment. Choisissez comment l’empilement doit être calculé [la méthode par défaut est a ] :

    • a = moyenne (moyenne)

    • m = médiane

    • p = mode (vraisemblance maximale)

    • l = inférieur

    • L = inférieur mais ne considère que les valeurs positives

    • u = supérieur

    • U = supérieur mais ne considère que les valeurs négatives.

    Les modificateurs contrôlent la sortie ; choisissez un ou plusieurs parmi ces choix : - +a : ajoute des valeurs empilées à tous les profils croisés. - +d : Ajoute des déviations de pile à tous les profils croisés. - +r : Ajoute des résidus de données (données - pile) à tous les profils croisés. - +s [ fichier ] : enregistre le profil empilé dans un fichier [le nom de fichier par défaut est grdtrack_stacked_profile.txt]. - +c fact : calcule l’enveloppe sur le profil empilé sous la forme ± fact * déviation [la valeur factuelle par défaut est 2].

    Remarques:

    1. Les écarts dépendent de la méthode et sont st.dev ( a ), échelle L1, c’est-à-dire 1,4826 * écart absolu médian (MAD) (pour m et p ), ou demi-intervalle (supérieur-inférieur)/2.

    2. Le fichier de profil empilé contient une colonne principale plus des groupes de 4 à 6 colonnes, avec un groupe pour chaque grille échantillonnée. La colonne principale contient la distance transversale, tandis que les quatre premières colonnes d’un groupe contiennent respectivement la valeur empilée, l’écart, la valeur minimale et la valeur maximale. Si la méthode fait partie d” a | m | p puis on écrit aussi les bornes de confiance inférieure et supérieure (voir +c ). Lorsqu’un ou plusieurs de +a , +d et +r sont utilisés, nous ajoutons également les résultats d’empilement à la fin de chaque ligne, pour tous les profils croisés. L’ordre est toujours la valeur empilée ( +a ), suivie des écarts ( +d ) et enfin des résidus (+r ). Lorsque plus d’une grille est échantillonnée, cette séquence de 1 à 3 colonnes est répétée pour chaque grille.

  • radius ( bool ou int ou float ou str ) – [ rayon ][ +e | p ]. À utiliser avec un échantillonnage de grille normal et limité à une seule grille non IMG. Si le nœud le plus proche du point d’entrée est NaN, recherchez vers l’extérieur jusqu’à ce que nous trouvions le nœud non-NaN le plus proche et rapportons cette valeur à la place. Spécifiez éventuellement un rayon de recherche qui limite la prise en compte aux points situés à cette distance du point d’entrée. Pour indiquer l’emplacement du nœud le plus proche et sa distance par rapport au point d’entrée, ajoutez +e . L’unité par défaut pour les distances de grille géographique est le degré sphérique. Utilisez le rayon pour changer l’unité et donnerradius = 0 si vous ne souhaitez pas limiter la recherche de rayon. Pour remplacer à la place le point d’entrée par les coordonnées du nœud le plus proche, ajoutez +p .

  • 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

  • z_only ( bool ) – N’écrivez que les valeurs z échantillonnées [la valeur par défaut écrit toutes les colonnes].

  • aspatial ( bool ou str ) – [ col =] nom [,…]. Contrôlez la manière dont les données aspatiales sont gérées lors de l’entrée et de la sortie.

  • binary ( 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’une des suivantes :

    • c - int8_t (caractère signé 1 octet)

    • u - uint8_t (caractère non signé de 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é de 4 octets)

    • l - int64_t (entier signé sur 8 octets)

    • L - uint64_t (entier non signé de 8 octets)

    • f - Flottant simple précision de 4 octets

    • d - Flottant double précision de 8 octets

    • x - utilisé pour ignorer les ncols n’importe où dans l’enregistrement

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

    • w après n’importe quel élément pour forcer l’échange d’octets.

    • +l | b pour indiquer que l’ensemble du fichier de données doit être lu respectivement en petit ou gros boutien.

  • 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.

  • 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).

  • espace ( str ou liste ) – x | y | z | d | X | Y | D écart [ u ][ +a ][ +c col ][ +n | p ]. Examinez l’espacement entre les points de données consécutifs afin d’imposer des ruptures de ligne. Pour spécifier plusieurs critères, fournissez une liste avec chaque élément contenant une chaîne décrivant un ensemble de critères.

    • x | X - définit un écart lorsqu’il y a un changement suffisamment important dans les coordonnées x (majuscules pour utiliser les coordonnées projetées).

    • y | Y - définit un écart lorsqu’il y a un changement suffisamment important dans les coordonnées y (majuscules pour utiliser les coordonnées projetées).

    • d | D - définit un écart lorsqu’il y a une distance suffisamment grande entre les coordonnées (majuscules pour utiliser les coordonnées projetées).

    • z - définit un espace lorsqu’il y a un changement suffisamment important dans les données z. Utilisez +c col pour modifier la colonne de données z [ la colonne par défaut est 2 (c’est-à-dire la 3e colonne)].

    Une unité u peut être ajoutée à l” écart spécifié : - Pour les données géographiques ( x | y | d ), l’unité peut être arc- d (egrees), m (inutes) et s (econds) , ou (m) e (ters), f (eet), k (ilomètres) ), M (iles) ou n (miles automatiques) [La valeur par défaut est (m) e (ters)]. - Pour les données projetées ( X | Y | D ), l’unité peut être i (inches), c (entimètres) ou p (points). Ajoutez le modificateur +a pour spécifier que tous les critères doivent être remplis [la valeur par défaut impose des pauses si l’un des critères est rempli].

    L’un des modificateurs suivants peut être ajouté :

    • +n - spécifie que la valeur précédente moins la valeur de la colonne actuelle doit dépasser l’écart pour qu’une rupture soit imposée.

    • +p - spécifie que la valeur actuelle moins la valeur précédente doit dépasser l’écart pour qu’une pause soit imposée.

  • 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].

  • distcalc ( str ) – e | f | g . Déterminez comment les distances sphériques sont calculées.

    • e - Mode ellipsoïdal (ou géodésique)

    • f - Mode Terre plate

    • g - Distance orthodromique [Défaut]

    Tous les calculs de distance sphérique dépendent de l’ellipsoïde courant ( PROJ_ELLIPSOID ), de la définition du rayon moyen ( PROJ_MEAN_RADIUS ) et de la spécification du type de latitude ( PROJ_AUX_LATITUDE ). Les calculs de distance géodésique sont également pilotés par la méthode ( PROJ_GEODESIC ).

  • interpolation ( str ) – [b|c|l|n][+a][+bBC][+c][+t threshold]. Sélectionnez le mode d’interpolation pour les grilles. Vous pouvez sélectionner le type de spline utilisé :

    • b pour B-spline

    • c pour bicubique [Défaut]

    • l pour bilinéaire

    • n pour le plus proche voisin

  • outcols ( str ou tableau 1-D ) – cols [,…][, t [ mot ]]. Spécifiez les colonnes de données pour la sortie 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 écrites 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 de sortie (par exemple, outcols=[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 est par défaut égal à 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, outcols= »0:2,4 » pour afficher les trois premières colonnes suivies du 5ème colonne). Pour écrire à 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 écrire du texte de fin, ajoutez la colonne t . Ajoutez le numéro de mot à t pour n’écrire qu’un seul mot à partir du texte de fin. Au lieu de spécifier des colonnes, utilisez outcols= »n » pour lire simplement l’entrée numérique et ignorer le texte de fin.

    Remarque : Si incols est également utilisé alors les colonnes indiquées pour outcols correspondre à la commande après que la incols sélection a eu lieu.

  • skiprows ( bool ou str ) – [ colonnes ][ +a ][ +r ]. Supprimer la sortie pour les enregistrements dont la valeur z est égale à NaN [la valeur par défaut génère tous les enregistrements]. Éventuellement, fournissez une liste séparée par des virgules de toutes les colonnes ou plages de colonnes à prendre en compte pour ce test NaN [la valeur par défaut ne prend en compte que la troisième colonne de données (c’est-à-dire, cols = 2 )]. Les plages de colonnes doivent être données au format start [: inc ]: stop , où inc est par défaut égal à 1 s’il n’est pas spécifié. Les modificateurs suivants sont pris en charge :

    • +r pour inverser la suppression, c’est-à-dire ne sortir que les enregistrements dont la valeur z est égale à NaN.

    • +a pour supprimer la sortie de l’enregistrement si une ou plusieurs des colonnes sont égales à NaN [la valeur par défaut ignore l’enregistrement uniquement si les valeurs de toutes les colonnes spécifiées sont égales à NaN].

  • wrap ( str ) – y | un | w | d | h | m | s | c période [/ phase ][ +c col ]. Convertissez la coordonnée x d’entrée en une coordonnée cyclique, ou une colonne différente si elle est sélectionnée via +c col . 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 - deuxième cycle (seconde)

Retour

  • track ( pandas.DataFrame ou None ) – Le type de retour dépend si le outfile paramètre est défini :

    • pandas.DataFrame table avec (x, y, …, newcolname) si outfile n’est pas défini

Aucun si outfile est défini (la sortie de la piste sera stockée dans le fichier défini par outfile)

Exemple

import pygmt
# Load a grid of @earth_relief_30m data, with an x-range of -118 to
# -107, and a y-range of -49 to -42
grid = pygmt.datasets.load_earth_relief(
    resolution="30m", region=[-118, -107, -49, -42]
)
# Load a pandas dataframe with ocean ridge points
points = pygmt.datasets.load_sample_data(name="ocean_ridge_points")
# Create a pandas dataframe from an input grid and set of points
# The output dataframe adds a column named "bathymetry"
output_dataframe = pygmt.grdtrack(
    points=points, grid=grid, newcolname="bathymetry"
)