Articles of Interpolation

Diferencia entre interpolación spline cuadrática y de segundo orden en scipy

Estoy escribiendo funciones que calcularán las interpolaciones 1d en python usando la función scipy.interpolate. Usando la ayuda de la documentación, escribí 2 funciones diferentes para la interpolación spline cúbica y cúbica # calculate cubic interpolation def linear_interpolation(x): linear = interpolate.interp1d(support_x, support_y, ‘cubic’) return linear(x) # calculate cubic spline interpolation def cubic_spline_interpolation(x): tck = interpolate.splrep(support_x, support_y) […]

analogía a scipy.interpolate.griddata?

Quiero interpolar una nube de puntos 3D dada: Eché un vistazo a scipy.interpolate.griddata y el resultado es exactamente lo que necesito, pero según tengo entendido, debo ingresar “griddata” que significa algo como x = [[0,0,0],[1,1,1],[2,2,2]] . Pero mi nube de puntos 3D dada no tiene este aspecto de cuadrícula: los valores x, y no se […]

Interpolación lineal rápida en Numpy / Scipy “a lo largo de un camino”

Digamos que tengo datos de estaciones meteorológicas en 3 altitudes (conocidas) en una montaña. Específicamente, cada estación registra una medición de temperatura en su ubicación cada minuto. Tengo dos tipos de interpolación que me gustaría realizar. Y me gustaría poder realizar cada uno rápidamente. Así que vamos a configurar algunos datos: import numpy as np […]

Encuentra la distancia mínima de punto a curva complicada

Tengo una curva complicada definida como un conjunto de puntos en una tabla como tal (la tabla completa está aquí ): # xy 1.0577 12.0914 1.0501 11.9946 1.0465 11.9338 … Si trazo esta tabla con los comandos: plt.plot(x_data, y_data, c=’b’,lw=1.) plt.scatter(x_data, y_data, marker=’o’, color=’k’, s=10, lw=0.2) Me sale lo siguiente: donde he añadido los puntos […]

Interpolación de funciones con valores vectoriales usando NumPy / SciPy

¿Hay alguna manera de interpolar una función de vector con NumPy / SciPy ? Hay muchas ofertas que funcionan en funciones de valor escalar, y creo que puedo usar una de ellas para estimar cada componente del vector por separado, pero ¿hay alguna manera de hacerlo de manera más eficiente? Para elaborar, tengo una función […]

Interpolar los valores de Z en una superficie 3D, a partir de un conjunto irregular de puntos

Tengo una superficie que se parece a la Figura A, imagino que es una vista desde arriba. La superficie ha calculado el valor z. Ahora necesito encontrar todos los valores de Z en puntos nuevos como la figura B. ¿Cómo hacer esto? Intenté scipy.interpolate.interp2d pero da algunos resultados extraños como este: Solo quiero encontrar z […]

Scipy griddata no funciona dentro de un bucle / pérdida de memoria

Tengo un problema al utilizar los datos de cuadrícula de Scipy dentro de un bucle. Básicamente, lo que sucede es que la memoria crece sin límite mientras el bucle se está ejecutando. Para reproducir el problema solo pon el ejemplo en http://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.griddata.html dentro de un bucle: for i in range(100000): grid_z1 = griddata(points, values, (grid_x, […]

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 […]