Interpolation 3D et extrapolation

Bonjour,

Travaillant sur le reformatage d'images, je dois utiliser des interpolations 3D et faire également une extrapolation à partir de mes données.

Mes données sont éparses, le problème en 2D a été résolu grâce à la fonction griddata qui permet l'interpolation sur grille non régulière et également l'extrapolation des mes données grâce à la méthode 'v4' .

Devant désormais passer en 3D, la fonction griddata3 quant à elle ne propose que les méthodes 'nearest' ou 'linear' pour l'interpolation (des plus perfectionnées sont disponible avec griddata), de plus griddata3 ne permet pas d'extrapolation (absence de la méthode 'v4').

Connaissez-vous une solution à mon problème ou une fonction développée qui permettrait cela ? (j'ai regardé dans le file exchange mais rien trouvé...)

Merci de votre aide

cette solution ?

D'après le lien envoyé :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 in1=[1:2]; a(:,:,1,1)=ones(2); a(:,:,1,2)=2*ones(2); a(:,:,2,1)=2*ones(2); a(:,:,2,2)=4*ones(2); testin=[10*rand(1,4)]; %interpn returns a numerical value if the method 'spline' is chosen and the %point lies outside the range defined by the look-up table extspl=interpn(in1,in1,in1,in1,a,... testin(1),testin(2),testin(3),testin(4),'spline') %interpn returns NaN if the method 'linear' is chosen and the %point lies outside the range defined by the look-up table extlin=interpn(in1,in1,in1,in1,a,... testin(1),testin(2),testin(3),testin(4),'linear')
L'extrapolation se réalise donc toute seule suivant la méthode d'interpolation choisie, ici dans ce cas l'interpolation spline ?

Deux questions se posent alors:

- cette extrapolation se base sur quoi et est-elle valable ?
- dans mon cas j'ai des données non régulières, donc j'utilise griddata3...comment réaliser donc de même avec cette fonction ?

Merci