Discussion :

[C_Rongieras]

Bonjour,

Je ne suis pas vraiment experte sur Scilab…
Je dois analyser des données de position enregistrées avec le système VICON.
Dans plusieurs fichiers je me retrouve avec des «*valeurs*» NaN dans certaines colonnes.
Cela n’arrive pas à des moment de l’enregistrement qui m’intéressent en revanche ça bloque le filtre que j’applique puisque dés qu’il y a un NaN dans une colonne le filtre me transforme toute la colonne en NaN.

J’ai vu la fonction interp1 mais malgré les aides je ne sais pas comment l’appliquer dans le cas d’enregistrement de positions

Merci pour votre aide.


CR

[Jerome Briot]

Cela n’arrive pas à des moment de l’enregistrement qui m’intéressent

Dans ce cas, tu peux remplacer ces NaN par n'importe quelle valeur.

Si les NaN sont isolés tu peux par exemple faire ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
X = [1 2 3 %nan 5 6 %nan 8 9]
 
disp(X)
 
for i = 1:length(X)-1
    if isnan(X(i)) then
        X(i) = (X(i-1)+X(i+1))/2
    end
end
 
disp(X)

Ce qui donne :

Code :

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

1
2
3
   1.   2.   3.   Nan   5.   6.   Nan   8.   9.
 
   1.   2.   3.   4.   5.   6.   7.   8.   9.

Si les NaN sont consécutifs, tu peux par exemple faire ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
X = [1 2 %nan %nan 5 %nan %nan 8 9]
 
disp(X)
 
i = length(X)-1
 
while i > 0
    if isnan(X(i)) then
        X(i) = X(i+1)
    end    
    i = i-1
end
 
disp(X)

Ce qui donne :

Code :

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

1
2
3
   1.   2.   Nan   Nan   5.   Nan   Nan   8.   9.
 
   1.   2.   5.   5.   5.   8.   8.   8.   9.

Si tu veux vraiment utiliser une fonction d'interpolation, voici comment faire :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
X = [11 12 %nan %nan 15 %nan %nan 18 19]
 
disp(X)
 
i = find(~isnan(X))
ii = 1:length(X)
X = X(~isnan(X))
 
Xi = interp1(i, X, ii)
 
disp(X)
disp(i)
disp(ii)
disp(Xi)

Ce qui donne :

Code :

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

1
2
3
4
5
6
7
8
9
   11.   12.   Nan   Nan   15.   Nan   Nan   18.   19.
 
   11.   12.   15.   18.   19.
 
   1.   2.   5.   8.   9.
 
   1.   2.   3.   4.   5.   6.   7.   8.   9.
 
   11.   12.   13.   14.   15.   16.   17.   18.   19.

[C_Rongieras]

Merci beaucoup pour votre réponse je vais essayer ça