Discussion :

[Kokiett]

Bonjour bonjour,

Je m'arrache les cheveux.... Tout ce que je trouve sur internet c'est pour supprimer les doublons, ahah, moi je suis bête, je veux en créer!!

Je veux identifier les flotteurs présents dans les tourbillons d'une trajectoire.
Pour chaque trajectoire, je récupère la date de tous les pas de temps de cette trajectoire, ce qui me donne par exemple un vecteur monotone

Code :

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

map = [1021:1:1023];

(exemple minuscule pour pouvoir mieux visualiser)
Je veux ensuite tous les flotteurs qui sont sur chaque pas de temps, mais sans faire une boucle for parce que j'ai beaucoup de données à traiter...
Du coup j'utilise:

Code :

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

[aa,bb] = ismember(date_ALM22042,map);

Comme length(date_flotteurs) = 1499, je vais avoir aa et bb de taille [1 1499]. aa contient beaucoup de 0 parce que le vecteur map est petit (du coup je vire toutes les lignes avec des 0). Mais si j'ai 2 fois une occurence d'une date de la trajectoire dans mon set date_flotteur, j'ai par exemple

Code :

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

1
2
3
4
bb = [ 1021
         1023
         1022
         2021]

J'ai un autre vecteur qui contient le numéro des tourbillons, chaque ligne correspond à chaque valeur de date

Code :

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

1
2
3
eddy = [65
           74
           22]

And the question is: comment je peux créer un autre vecteur de tourbillon, eddy2, qui contient les valeurs de eddy en répétant cas valeurs quand la date correspondante a été répétée?? Et tout ça bien sur sans faire de boucle for...

Code :

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

1
2
3
4
eddy2 = [65
             22
             74
             65]

Si quelqu'un a une idée (ou sait où trouver de la doc...), il me sauvera de la chauvitude ^^ !
(et recevra une bonne de reconnaissance)

[Jerome Briot]

Comme ceci :

Code :

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

1
2
3
[ubb, a, b] = unique(bb);
 
eddy2 = eddy(b);

[Kokiett]

Merci beaucoup!
En effet avec unique ça marche (j'en étais restée à la première fonctionnalité sans aller voir plus loin...)

Cette méthode est aussi plus facile que l'autre solution que j'ai trouvée finalement...

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
map_traj = [1020 1021 1022 1023]; % le petit vecteur des maps connues
eddy_traj = [17 55 42 39] % le numéro des tourbillons associés
 
map_float = [1022 1023 1020 1020 1021 1023] % le grand vecteur des maps
eddy_float = [ 42 39 17 18  55 39] % le numéro des tourbillons associés
 
a = dataset(map_traj',eddy_traj','varnames',{'key1','eddy'})
 
b = dataset(map_float',eddy_float','varnames',{'key1','eddyf'})
 
 c = join(b,a,'key1')
 
 c = 
 
    key1    eddyf    eddy
    1022    42       42  
    1023    39       39  
    1020    17       17  
    1020    18       17  
    1021    55       55  
    1023    39       39  
 
% la colonne eddy répète les valeurs de eddy_traj comme je voulais et je peux identifier les tourbillons communs!!

Donc merci Dut, comme promis je t'envoie une dose de reconnaissance (j'espère que tu vas la recevoir!)