Discussion :

[mylia04]

Bonjour,

Je suis toute nouvelle sur le forum et ravie de faire partie de l'aventure!

Voila mon soucis:

J'aimerais remplacer les NaN d'une matrice A=[0 0 0 0 0;0 0 45 48 NaN;0 43 42 NaN 0;50 NaN NaN NaN 45]; par la moyenne des plus proches voisins.

Auriez-vous des suggestions?

Merci d'avance

[Jerome Briot]

Il nous faudrait une définition exacte de ce que entends pas "plus proches voisins"

[Jerome Briot]

Une première idée :

Code :

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

1
2
3
4
5
6
7
8
9
A = [0 0 0 0 0;0 0 45 48 NaN;0 43 42 NaN 0;50 NaN NaN NaN 45]
 
idx = find(~isnan(A(:)));
 
[i,j] = ind2sub(size(A),idx);
 
[yi,xi] = meshgrid(1:size(A,2),1:size(A,1));
 
Ai = griddata(i,j,A(idx),xi,yi)

donne :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
A =
 
     0     0     0     0     0
     0     0    45    48   NaN
     0    43    42   NaN     0
    50   NaN   NaN   NaN    45
 
 
Ai =
 
         0         0         0         0         0
         0         0   45.0000   48.0000         0
         0   43.0000   42.0000   21.0000         0
   50.0000   48.7500   47.5000   46.2500   45.0000

[mylia04]

Bonjour,

C'est un peu complique :

Tout dépend de la position du NaN. Dans l’idéal, les plus proches voisins sont les valeurs formant un 'carre' (3 lignes * 3 colonnes) autour de la valeur a laquelle on s’intéresse: Si le focus est sur A(3,2)=43, alors le 'carre' qui l'entoure serait :

0    0     45
0    43    42
50  NaN   NaN

Mais, aux bords de la matrice, le 'carre' de ne vaudrait que (2 lignes * 2 colonnes)

Exemples : Si i et j sont les indices du NaN se trouvant à A(4,2); j'aimerais remplacer ce NaN par la moyenne des valeurs non nulles a i-1 et j-1:

Pour :

0    43 
50  NaN

Ici le NaN serait valeurnan1= (50+43)/2= 46.5.

Connaissant à présent valeurnan1; on pourrait de manière récursive connaitre le NaN a A(4,3) par:

43  +  42         
    +
46.5  NaN

Ici valeurnan2 = (43+42+46.5)/3= 43.833 .

Je ne sais pas du tout si c'est clair !

[mylia04]

Merci beaucoup pour cette premiere reponse qui interpole ligne par ligne.
Y aurait-il egalement un moyen pour interpoler dans toutes les directions, comme le montre l'exemple ci-dessus?

Merci par avance