Discussion :

[Titened]

Bonjour à tous,

Je n'ai pas trouvé (ou peut être pas compris les réponses...) la solution à mon problème sur le forum.
J'ai un nuage de points type [lon, lat, z].
Et j'ai une grille régulière format meshgrid pour des latitudes (XX) et des longitudes (YY).

Le but du jeu, c'est d'avoir une valeur de z sur les points de ma grille et non un nuage de point. Donc d'interpoler mes données sur cette grille prédéfinie.

Comme lon, lat, z sont des vecteurs, je voulais utiliser triscatterinterp, sauf que XX et YY de mon maillages sont des matrices....
Est-ce que quelqu'un a une petite idée de comment procéder ??

Autre petit soucis, z comporte des NaN, comment triscatterinterp se comporte avec les NaN ?

Merci beaucoup à mon(mes) futur(s) sauveur(s) !!

[Jean Dumoncel]

Bonjour,

triscatterinterp permet de définir la méthode d'interpolation, pas d'interpoler directement, tu peux faire par exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

F = TriScatteredInterp(lat,lon,z);

puis :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

qz = F(XX,YY);

Si tu as une version récente de matlab, il est préférable d'utiliser scatteredInterpolant car triscatterinterp est voué à être supprimé.

Pour les NaN je ne sais pas, il sera peut-être nécessaire de les retirer avant de faire l'interpolation.

[Titened]

Super ! Ça marche très bien. Merci beaucoup !!

Pour ceux qui sont dans le même cas que moi, pour éviter les NaN, j'ai extrait dans un vecteur la zone où je veux que l'interpolation se fasse. Etant donné que la-dite zone est de forme irrégulière mais que j'avais les coordonnées du contour de cette zone, j'ai utilisé la fonction :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

[in, on] = inpolygon (XX, YY, lat_polygone, lon_polygone)