Discussion :

[thecrazydonut]

Bonjour,

Voilà maintenant quelques heures que je me creuse la tête pour découper un vecteur d'une façon particulière.

J'ai par exemple un vecteur Qobs90s (il s'appelle comme ça chez moi) qui est défini comme tel :

Code :

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

Qobs90s=[31 22 71 -999 -999 12 31 11 -999 -999 17 -999 -999 21 41 -999]

Ce que je veux, c'est extraire l'emplacement des chiffres précédents et suivants les -999 mais aussi ajouter 1 dans le cas où la série de commence pas par un -999 et ajouter une length(Qobs90s) dans le cas où la série ne se termine pas par un -999. Donc ici la réponse désirée serait :

Code :

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

1
2
3
4
5
idx =
1 3
6 8
11 11
14 15

J'ai écris ce qui se trouve dessous mais je n'arrive pas à surmonter le problème où il y a deux -999 séparés seulement d'une seule valeur positive (ça fonctionne mais j'ai un nombre d'élément impair dans mon vecteur idx ce qui m’empêche d'utiliser le reshape). C'est un bout de code qui devrait préparer des bases de données et je voudrais donc qu'il soit robuste.

Le petit bout de code :

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
24
25
26
27
28
29
Qobs90s=Qobs90s';
idx=1;
for m=1:length(Qobs90s)-1
    if Qobs90s(m)<0
        if Qobs90s(m-1)>0
            idx=[idx,m-1];
            if Qobs90s(m+1)>0
                idx=[idx,m+1];
            end
        end
        if Qobs90s(m+1)>0
            idx=[idx,m+1];
        end
    end
end
 
idx=unique(idx);
 
if idx==1
    idx=[1;length(Qobs90s)];
else
    if Qobs90s(end)<0
        idx=reshape(idx,2,[]);
    else
        idx=[idx,length(Qobs90s)];
        idx=reshape(idx,2,[]);
    end
end
disp((idx'));

Merci beaucoup à ceux qui auront des idées (et à ceux qui n'en n'auront pas mais qui auront pris le temps de lire ce message )

[shaiHulud]

Bonjour,

Code :

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

1
2
3
4
5
6
7
 
i= find(Qobs90s == -999); %indices des -999
idx= [i(:)-1; i(:)+1]; % indices avant et après les -999
idx= idx(idx>0 & idx <= length(Qobs90s ));
idx= setdiff(idx,i); %enlève les -999 adjacents à un -999 (pas sur si c'est demandé)
 
data= Qobs90s (idx);

@FLB: ma méthode va effectivement doubler certaines valeurs, il faut rajouter

Code :

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

idx= unique(idx); % enlève les doublons

[FLB]

Salut,
je me permet d'ajouter une autre solution, moins sexy mais qui permet d'avoir le 11 en double.
Je me suis inspiré du post http://www.developpez.net/forums/d89...zeros-vecteur/

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
%On commence par s'assurer de l'état des bords du vecteur
if Qobs90s(1)>0
     Qobs90s = [-999 Qobs90s];
     enleverUn = 1;
else
    enleverUn = 0;
end
if Qobs90s(end)>0
     Qobs90s = [Qobs90s -999];
end
snz=diff(find(~[0 Qobs90s>0 0]))-1 ; 
numElem = snz(snz~=0);%on compte les suites de nombres positifs
liste = find(Qobs90s>0); %pour aller chercher l'indice du premier et dernier element de chaque suite de nombres positifs
idstart = cumsum([1 numElem(1:end-1)]);
idend = cumsum(numElem(1:end));
[liste(idstart)' liste(idend)']-enleverUn

[thecrazydonut]

Hehe ! Il semblerait en fait que ce ne soit pas si simple que ce que je pensais. Le code que FLB propose est en effet assez "sexy" mais il présente exactement le même problème que le mien, c'est à dire que pour une séquence -999 12 -999 il va me rendre [2] au lieu de [2,2]. Le deuxième code semble présenter quelque problèmes au niveau des indices...

+astucieux l'usage du setdiff

[thecrazydonut]

@FLB: ma méthode va effectivement doubler certaines valeurs, il faut rajouter

Code :

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

idx= unique(idx); % enlève les doublons

Tu es certain que ta fonction puisse doubler les valeurs qui sont entourées de -999 de part et d'autre ? Je reprend l'exemple que j'ai posté dans mon message précédent : si on prend Qobs90s = -999 12 -999, idx= 2 et non pas 2 2 (ce que je désire)

[FLB]

si je lance mon code sur

Code :

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

Qobs90s=[-999 12 -999];

J'obtiens

Code :

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

1
2
ans =
     2     2

C'est bien ce que tu veux, non?

[Invité]

Bonsoir,

Une autre solution :

Code :

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

1
2
3
4
5
Qobs90s = [31 22 71 -999 -999 12 31 11 -999 -999 17 -999 -999 21 41 -999];
idxBefore = strfind([Qobs90s==-999 true], [false true]);
idxAfter = strfind([true Qobs90s==-999], [true false]);
 
[idxAfter(:) idxBefore(:)]

[thecrazydonut]

si je lance mon code sur

Code :

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

Qobs90s=[-999 12 -999];

J'obtiens

Code :

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

1
2
ans =
     2     2

C'est bien ce que tu veux, non?

Oui c'est ce que je veux... mais pas ce que j'obtiens.

Code :

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

1
2
3
ans =
 
     2

[thecrazydonut]

Merci à tous. La solution de Winjerome me donne les résultats escomptés. (passer de la 20aine de ligne de mon code à 3 lignes qui fonctionnent mieux... je me dis que j'ai encore du chemin à faire ! )

@VV33D. Bha... je comprends pas ce qui se passe !

[shaiHulud]

Tu es certain que ta fonction puisse doubler les valeurs qui sont entourées de -999 de part et d'autre ?

Je pensais au cas d'indices uniques. setdiff(a,b) lorsqu'il y a des doublons dans a enlève ces doublons pour être cohérent avec la notion d'ensemble (il me semble). Cela me permet aussi d'ignorer ce qu'il se passe au bord : si 1 est dans i lorsque je fais

Code :

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

idx= [i(:)-1; i(:)+1]; et idx(idx>0 & idx <= length(Qobs90s ));

Pour garder les doublons, je préférerais

Code :

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

1
2
3
4
5
 
idx= find(values== -999);
idx= [idx-1, idx+1];
idx= idx(idx > 0 & idx<= length(values));
idx( values(idx) == -999 )= [];