Discussion :

[Gooby]

Bonjour,

J'ai de nombreux vecteurs colonne d'environ 2000 points. Ces vecteurs colonnes prennent des valeurs comprises entre 0 et 40. Je cherche à déterminer les points de mon vecteur colonne inférieur à, par exemple, la valeur 5.

Pour ce faire je peux utiliser simplement find(monVecteur < 5) et j'obtiens, sur un exemple:

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
         677
         678
        1102
        1104
        1105
        1107
        1108
        1109
        1110
        1111
        1112
        1113
        1114
        1120
        1121
        1456
        1460
        1461
        1462
        1463

Cependant, je ne m'intéresse qu'aux larges plages de points < 5, par exemple 30 points en dessous de la valeur 5 (il n'y en a donc pas dans cet exemple ci )

Je pourrais vérifier ça avec des for, des if, des while, mais ça donnerait pas un code super propre et ça prendrait du temps de calcul.

Ma question est : connaissez vous un moyen pour vérifier que j'ai bien 30 points successifs < 5 dans mon vecteur?

En vous remerciant par avance,
Gooby

[Jean Dumoncel]

Bonjour,

tu peux utiliser un truc du genre :

Code :

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

diff(find(monVecteur < 5))==1

pour trouver les points successifs (attention, le résultat possède un élément de moins que le résultat de find(monVecteur < 5) à cause du diff.

Ensuite tu peux t'inspirer de Taille et nombre d'éléments non zéros dans un vecteur

[Gooby]

Merci pour la réponse:

En prenant <=2 pour ne pas prendre en compte le fait qu'un point isolé peut se retrouver >5 dans une plage, j'obtiens:

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
diff(find(distance_contour<5))<=2
 
ans =
 
     1
     0
     1
     1
     1
     1
     1
     1
     1
     1
     1
     1
     0
     1
     0
     0
     1
     1
     1

La question maintenant est: Comment puis-je faire pour identifier les indices de ma plage < 5 ?

[Jean Dumoncel]

Je n'ai pas compris pourquoi tu prends <=2, il ne faut pas que les points se suivent?

Pour la suite, pas facile à expliquer, il y a plusieurs façons de le faire, en voici une (dans mon test, distance_contour est un vecteur ligne):

Code :

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

1
2
3
4
5
6
7
8
9
d= find(distance_contour<5);
e =  (diff([d(1) d]))==1;
 
snz = diff(find(~[0 e 0]))-1 ;
snz = snz(snz~=0) ;
 
debut = find(diff([0 e])==1);
d(debut)-1
snz+1

d(debut)-1 contient les indices des premiers éléments de chaque suite et snz+1 contient les nombres d'éléments consécutifs.

[Gooby]

J'ai pris <= 2 dans le cas où par exemple j'aurai 15 valeurs successives <5 sauf une se trouvant au milieu qui serait par exemple égale à 5.1. Cela ferait que mon intervalle passerait inaperçu pour une valeur ne remplissant pas la condition. Le <= 2 aurait permis de "sauter" cette valeur. Je ne sais pas si c'est très clair.

En tout cas merci pour ta contribution, cela fonctionne très bien. Je dois cependant dire qu'après un bon moment d'analyse, je me perds à cette ligne snz = diff(find(~[0 e 0]))-1 ;

[Jean Dumoncel]

Les 2 premières lignes permettent de repérer les éléments successifs (sauf les premiers éléments, ce qui pourrait être amélioré). Il faut alors repérer les indices des premiers élément et les longueurs de chaque séquence de 1. Le rôle des lignes 4 et 5 est de compter le nombre d'éléments (ces lignes viennent de la contribution cité dans mon premier message). Je te conseille de décomposer la ligne diff(find(~[0 e 0]))-1 et de regarder le résultat de chaque opération pour comprendre l'ensemble.