Desarrollo de Python
  1. Desarrollo de Python
  3. Interpolación RBF: LinAlgError: matriz singular
Intereting Posts
Asignar valores a varias columnas en Pandas Manera pythonica de evitar las declaraciones de “if x: return x” Convertir Array formateado como una cadena – Python Python Trailing L Problem Termine los hilos de python usando sys.exit () Python tkinter devuelve el valor de la función utilizada en el comando Matplotlib: agrega una barra de color a un objeto no asignable SQLAlchemy: selección de las columnas de un objeto en una consulta ¿Cómo firmo una solicitud POST utilizando HMAC-SHA512 y la biblioteca de solicitudes Python? En Python, ¿puedes llamar a un método de instancia de clase A, pero pasar en una instancia de clase B? python & suds “ImportError: no se puede importar el nombre getLogger” ¿Cómo fusionar listas en una lista de tuplas? Cómo asignar columnas mientras se ignora la alineación del índice ¿Cómo resaltar la línea actual de un widget de texto? ¿Cuáles son las diferencias entre Conda y Anaconda?

Interpolación RBF: LinAlgError: matriz singular

La siguiente convocatoria: 

rbf = Rbf(points[0], points[1], values,epsilon=2)

resulta en un error: 

 LinAlgError: singular matrix

con los siguientes valores: 

 In [3]: points Out[3]: (array([71, 50, 48, 84, 71, 74, 89, 76, 70, 77, 74, 79, 83, 71, 72, 78, 73, 84, 75, 65, 73, 82, 48, 86, 74, 86, 66, 74, 68, 74, 81, 74, 88, 66, 57, 50, 72, 86, 72, 92, 81, 67, 82, 78, 69, 70, 73, 71, 76, 72, 74, 75]), array([32, 34, 4, 35, 1, 7, 47, 16, 37, 14, 65, 18, 32, 4, 3, 27, 25, 34, 18, 25, 6, 25, 34, 41, 16, 35, 44, 2, 32, 2, 37, 60, 45, 32, 33, 42, 54, 31, 18, 38, 24, 18, 45, 48, 9, 63, 56, 45, 9, 59, 5, 12])) In [4]: values Out[4]: array([ 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.])

¿Qué puedo hacer para evitarlo y resolver el problema de interpolación?

Creo que lo que estás tratando de hacer es estimar la densidad del kernel . Puedes usar scipy.stats.gaussian_kde para esto: 

 import numpy as np from scipy.stats import gaussian_kde from matplotlib import pyplot as pp # kernel density estimate of the PDF kde = gaussian_kde(points) # evaluate the estimated PDF on a grid x,y = np.mgrid[40:101,-20:101] z = kde((x.ravel(),y.ravel())).reshape(*x.shape) # plot fig,ax = pp.subplots(1,1) ax.hold(True) pc = ax.pcolor(x,y,z) cb = pp.colorbar(pc) cb.ax.set_ylabel('Probability density') ax.plot(points[0],points[1],'o',mfc='w',mec='k') pp.show()

introduzca la descripción de la imagen aquí

El módulo statsmodels también tiene algunas herramientas más elaboradas para la estimación de la densidad del kernel.

Tengo el mismo error. Finalmente encontré por qué me sale el error.

Tienes 2 puntos con la misma coordenada. (74,2) verifique que los valores tengan la misma coordenada en 28, 30.

En mi opinión, incluso si tiene el mismo valor en el mismo punto, emite un error de matriz singular.