pygmt.grdsample

pygmt.grdsample(grid, \*, outgrid=None, spacing=None, region=None,
translate=None, verbose=None, coltypes=None, interpolation=None,
registration=None, cores=None, \**kwargs)

Modifier l’enregistrement, l’espacement ou les nœuds dans un fichier de grille.

Cela lit un fichier de grille et l’interpole pour créer un nouveau fichier de grille. Il peut modifier l’enregistrement avec translate ou registration, modifier l’espacement de la grille ou le nombre de nœuds avec spacing et définir une nouvelle sous-région à l’aide de region. Une interpolation bicubique [par défaut], bilinéaire, B-spline ou du plus proche voisin est définie avec interpolation.

Lorsque regionest omis, la grille de sortie couvrira la même région que la grille d’entrée. Lorsque spacingest omis, l’espacement de grille de la grille de sortie sera le même que celui de la grille d’entrée. Soit registrationou translatepeut être utilisé pour modifier l’enregistrement de la grille. Lorsqu’il est omis, la grille de sortie aura le même enregistrement que la grille d’entrée.

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.

  • outgrid ( str ou None ) – Le nom du fichier netCDF de sortie avec l’extension .nc dans lequel stocker la grille.

  • 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émentation. Choisissez parmi m pour indiquer les arc-minutes ou s pour indiquer les arc-secondes. 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, km, mile, mile nautique ou US survey foot, respectivement, et sera converti en les 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é égal à x_inc; sinon, il sera converti en degrés de latitude.

    • Toutes les coordonnées : si +e est ajouté, le max correspondant x ( est ) ou y ( nord ) peut être légèrement ajusté pour s’adapter exactement à l’incrément donné [par défaut, l’incrément peut être légèrement ajusté pour s’adapter 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 alors recalculé à partir du nombre de nœuds, du registration et du domaine. La valeur d’incrément résultante varie selon que vous avez sélectionné une grille enregistrée sur la grille ou enregistrée sur les pixels ; voir Formats de fichier GMT pour plus de détails.

    Remarque : 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.

  • région ( str ou liste ) – xmin/xmax/ymin/ymax [ +r ][ +u unité ]. Spécifiez la région d’intérêt. L’ajout region sélectionnera une sous-section de la grille. Si cette sous-section dépasse les limites de la grille, seule la région commune sera sortie.

  • transalte ( bool ) - Traduire entre l’enregistrement de la grille et des pixels ; si l’entrée est enregistrée sur la grille, la sortie sera enregistrée sur les pixels et vice-versa.

  • registration ( str ) – g | p . Forcer l’enregistrement du nœud quadrillage ( g ) ou pixel ( p ) [la valeur par défaut est g (ridline)].

  • 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

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

  • 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

  • cores ( bool ou int ) – [[ - ] n ]. Limitez le nombre de cœurs à utiliser dans tous les algorithmes multithreads compatibles OpenMP. Par défaut, nous essayons d’utiliser tous les cœurs disponibles. Définissez un nombre n pour n’utiliser que n cœurs (s’il est trop grand, il sera tronqué au maximum de cœurs disponibles). Enfin, donnez un nombre négatif -n pour sélectionner (tous - n ) cœurs (ou au moins 1 si n est égal ou supérieur à tous).

Retour

  • ret ( xarray.DataArray ou None ) – Le type de retour dépend si le outgrid paramètre est défini :

    • xarray.DataArray si outgrid n’est pas défini

    • Aucun si outgrid est défini (la sortie de la grille sera stockée dans le fichier défini par outgrid)

Exemple

import pygmt
# Load a grid of @earth_relief_30m data, with an x-range of 10 to 30,
# and a y-range of 15 to 25
grid = pygmt.datasets.load_earth_relief(
    resolution="30m", region=[10, 30, 15, 25]
)
# Create a new grid from an input grid, change the registration,
# and set both x- and y-spacing to 0.5 degrees
new_grid = pygmt.grdsample(
    grid=grid, translate=True, spacing=[0.5, 0.5]
)