Discussion :

[Shokola]

Bonjour,

Je tente de faire une interpolation 1d avec la bibliothèque: scipy.interpolate, et j'utilise interp1d. Mais j'obtiens une erreurs: ValueError: A value in x_new is above the interpolation range

Voici mon code :

Code :

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

1
2
3
4
5
6
7
8
 
 
x = np.arange(0,len(matrice_defaut_acq[0]),1) 
for col in matrice_defaut_acq:
    f1d = interpolate.interp1d(x,col,kind='cubic') 
    xnew=np.linspace(0,len(matrice_defaut_acq[0]),len(matrice_defaut_bdd[0]))
    plt.figure()
    plt.plot(x,col,'r-',x,col,'ro',xnew,f1d(xnew),'co',xnew,f1d(xnew),'c--')

matrice_defaut_acq est la matrice que je récupère de mon capteur, et matrice_defaut_bdd est la matrice de ma base de données correspondant à un même défaut. Le pas de la matrice de la base de données est 100x plus fin que celle du capteur.

Merci d'avance

[Shokola]

Pour ceux qui ont le même problème que moi, il suffit de mettre le paramètre bounds_error à false dans la fonction interp1d:

Code :

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

1
2
3
4
5
6
7
 
x = np.arange(0,len(matrice_defaut_acq[0]),1) 
for col in matrice_defaut_acq:
    f1d = interpolate.interp1d(x,col,bounds_error=False,kind='cubic') 
    xnew=np.linspace(0,len(matrice_defaut_acq[0]),len(matrice_defaut_bdd[0]))
    plt.figure()
    plt.plot(x,col,'r-',x,col,'ro',xnew,f1d(xnew),'co',xnew,f1d(xnew),'c--')