Discussion :

[Falcosama]

Bonjour a tous,
est-ce que quelqu'un sait s'il existe une fonction Matlab permettant d'interpoler un matrice creuse? "interp2" ne va pas et pour l'instant je n'y arrive pas avec "griddata". Plus precisemment j'ai une matrice 17*6 mais seuleument 15 valeurs. La matrice doit etre "lisse". J'aimerais pouvoir trouver les valeurs manquantes (cubic splines...).

Merci d'avance!

[Jerome Briot]

Tu pourrais mettre cette matrice dans un fichier mat et l'attacher à un message (dans une archive zip ou rar) ?

[Falcosama]

Bien sur. J'ai mis la matrice creuse dans le m file. Les zero sont en fait des valeurs inconnues (c'est le probleme quand je charge la marice a partir de Excel, ca interprete le vide comme des zero et ensuite ca fait rater toute l'interpolation...). Pour etre plus clair, j'ai aussi colle dans le fichier excel le type de resultat que j'aimerais obtenir.
Merci

[Jerome Briot]

Essaie avec INPAINT_NANS :

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
30
A=[    0.5925         0         0         0         0         0
         0         0         0         0         0         0
         0         0    0.7695         0         0         0
         0         0         0         0         0         0
         0         0         0         0    0.8225         0
         0         0         0         0         0         0
    0.5458         0         0         0         0         0
         0         0    0.6405         0         0         0
         0         0         0         0         0         0
         0         0         0         0         0         0
         0         0         0         0         0         0
         0    0.4393         0         0         0         0
         0         0    0.4675         0         0    0.4775
         0         0         0         0         0         0
         0         0         0         0    0.4370         0
    0.3668         0         0         0         0         0]
 
 
idx = A==0;
 
A(idx) = nan;
 
Ai = inpaint_nans(A);
 
figure
surf(Ai)
hold on
 
[y,x] = find(~idx);
plot3(x(:),y(:),A(~idx),'r*')

Cette fonction est disponible sur le FEX : inpaint_nans - Interpolates (& extrapolates) NaN elements in a 2d array

[Falcosama]

Top! Je regarderai plus en details demain mais c'est exactement ce que je cherchais.
Merci pour l'aide, ca me depanne bien.

[Falcosama]

En y regardant de plus pres, j'ai une autre question.

Dans la matrice a interpoler, l'espacement des donnees n'est pas constant (cf. le fichier excel, l'increment des abscisses et ordonnees n'est pas contant). Or dans INPAINT_NANS(), on ne peut passer que "l'interieur" Z de la matrice (contrairement a INTERP2() ou GRIDDATA() ou on peut passer les X et Y aussi).

Alors je peux rajouter des NaNs pour m'aligner sur le plus petit increment et ainsi obtenir des abscisses/ordonnees a increments constants mais suivant les matrices, ca peut n'etre qu'approximatif et peut etre qu'il y a plus intelligent ou une autre fonction...

Merci