Discussion :

[rom3478]

Bonjour,
Je pense avoir un soucis dans l'interpolation au sein de ma fonction. Je vais essayer d'être clair:
- je possede des trajectoires définies par une latitude et une longitude pour chaque point
- je possede une map avec des vecteurs vents définis pour certaines positions
- je souhaite avoir le vecteur vent associé à chaque point de ma trajectoire

Problème:
- quand je regarde les valeurs de mes composantes vent en x et y, elles tournent autour de 0,1 à 0,9 m/s
- apres interpolation sur ma trajectoire, les composantes interpolés en x et y tournent sont toutes supérieures à 1 alors que je m'attendais à des valeurs toujours inférieure à 1

Je ne comprends pas pourquoi? Si quelqu'un a une idée sur mon probléme. Je vous joins le code de mon interpolation:

U est ma matrice où se trouve la valeur de mon vent (composante x)
X et Y correspondent à mon meshgrid (longitude et latitude) construit sur la position de chaque vecteur vent
x(i) et y(i) correspondent aux longitudes et latitudes pour chaque point de ma trajectoire

Code :

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

 uw=[uw;interp2(X,Y,U,x(i),y(i))];

Merci pour votre aide

[Jerome Briot]

Serait-il possible d'avoir un jeu de données réelles dans un fichier mat (mis dans une archive zip ou rar) ?

[ol9245]

Tu as du te gaufrer dans les données d'entrée parceque ta formule interp2 a l'air propre.

Sinon, apparemment tu fais une boucle sur l'indice i, ce qui est inutile : tu peux balancer ta trajectoire complète. cf la doc de interp2 :

Code :

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

1
2
XI and YI can be matrices, in which case interp2 returns the values of Z 
corresponding to the points (XI(i,j),YI(i,j)).