Articles of Interpolation

Representación de spline con scipy.interpolate: Interpolación deficiente para funciones de baja amplitud, que oscilan rápidamente

Necesito (numéricamente) calcular la primera y segunda derivada de una función para la que he intentado usar splrep y UnivariateSpline para crear splines con el fin de interpolar la función para tomar las derivadas. Sin embargo, parece que hay un problema inherente en la representación de spline en sí misma para las funciones cuya magnitud […]

Invertir interpolación para dar la variable asociada con un valor de función de interpolación deseado

Estoy intentando invertir una función interpolada utilizando la función de interpolación de scipy. Digamos que creo una función interpolada, import scipy.interpolate as interpolate interpolatedfunction = interpolated.interp1d(xvariable,data,kind=’cubic’) ¿Hay alguna función que pueda encontrar x cuando especifico un: interpolatedfunction(x) == a En otras palabras, “Quiero que mi función interpolada sea igual a a; ¿cuál es el valor […]

Alternativas a scipy.interpolate.griddata que no cuelgan en puntos alineados

Tengo un conjunto de datos de puntos que estoy tratando de interpolar en una cuadrícula. Estos puntos están alineados en forma de cuadrícula con algunos puntos faltantes, ver a continuación: Para complicarlo, es posible que otros conjuntos de puntos de entrada no se alineen en una cuadrícula, por lo que estoy tratando de usar scipy.interpolate.griddata […]

Error de tiempo de ejecución de griddata – Python / SciPy (Interpolación)

Uso la función griddate de scipy para la interpolación. ¿Qué significa el siguiente mensaje de error que aparece cuando python está ejecutando la función griddata? File “C:\Python25\lib\site-packages\scipy\interpolate\ndgriddata.py”, line 182, in griddata ip = LinearNDInterpolator(points, values, fill_value=fill_value) File “interpnd.pyx”, line 192, in interpnd.LinearNDInterpolator.__init__ (scipy\interpolate\interpnd.c:2524) File “qhull.pyx”, line 917, in scipy.spatial.qhull.Delaunay.__init__ (scipy\spatial\qhull.c:4030) File “qhull.pyx”, line 170, in […]

Interpolación ciencia de matriz grande.

Tengo un ndarray (Z) con unos 500000 elementos en una cuadrícula rectangular (X, Y). Ahora quiero interpolar valores en unas 100 ubicaciones en x, y que no están necesariamente en la cuadrícula. Tengo un código trabajando en Matlab: data = interp2(X,Y,Z, x,y); Sin embargo, cuando trato de usar el mismo enfoque con scipy.interpolate, obtengo varios […]

Python – Arreglo 2D de interpolación para arreglos enormes

Me gustaría interpolar la “prueba” de matriz 2D cuyas dimensiones son 4×4 (como ejemplo, en realidad cerca de 1000×1000) con una cuadrícula de forma 8×8. import numpy as np X = np.arange(0,4,1) Y = np.arange(0,4,1) points = np.vstack((X,Y)) points = points.T #my coordinates #my values as a 2D array test = np.array([[ 1.2514318 , 1.25145821, […]

Dado un conjunto de puntos definidos en (X, Y, Z) coordenadas, interpolar el valor Z en arbitrario (X, Y)

Dado un conjunto de puntos en coordenadas (X, Y, Z) que son puntos en una superficie, me gustaría poder interpolar valores Z en coordenadas arbitrarias (X, Y). He encontrado cierto éxito al usar mlab.griddata para interpolar valores en una cuadrícula, pero quiero poder llamar a una función de uso general para cualquier coordenada (X, Y). […]

Las pandas se interpolan dentro de un grupo

Tengo un dataframe con la siguiente información: filename val1 val2 t 1 file1.csv 5 10 2 file1.csv NaN NaN 3 file1.csv 15 20 6 file2.csv NaN NaN 7 file2.csv 10 20 8 file2.csv 12 15 Me gustaría interpolar los valores en el dataframe según los índices, pero solo dentro de cada grupo de archivos . […]

Interpolación circular acotada en python.

Mi pregunta es algo similar a la pregunta aquí . En el término simple tengo una serie de datos de ángulo de tiempo que está delimitada en [0, 360]. Necesito interpolar valores entre medidas. Actualmente estoy usando scipy.interpolate.interp1d . Para hacer mi pregunta aún más clara aquí hay un ejemplo de escenario, import numpy as […]

¿Qué significan los diferentes valores del argumento kind en scipy.interpolate.interp1d?

La documentación de SciPy explica que el kind argumento de interp1d puede tomar los valores ‘linear’ ‘zero’ , ‘slinear’ , ‘quadratic’ ‘zero’ , ‘slinear’ , ‘quadratic’ , ‘cubic’ . Los últimos tres son órdenes spline y ‘linear’ se explica por sí mismo. ¿Qué hacen ‘nearest’ y ‘zero’ ?