Articles of Interpolation

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’ ?

Error de memoria spline cúbica

En una computadora con 4GB de memoria, esta simple interpolación conduce a un error de memoria: (basado en: http://docs.scipy.org/doc/scipy/reference/tutorial/interpolate.html ) import numpy as np from scipy.interpolate import interp1d x = np.linspace(0, 10, 80000) y = np.cos(-x**2/8.0) f2 = interp1d(x, y, kind=’cubic’) Pensé en cortar los datos en partes, pero ¿hay alguna manera de poder realizar […]

Interpolación intuitiva entre puntos espaciados de manera desigual.

Tengo el siguiente gráfico que quiero digitalizar para obtener una figura de grado de publicación de alta calidad utilizando Python y Matplotlib: Utilicé un progtwig digitalizador para tomar algunas muestras de uno de los 3 conjuntos de datos: x_data = np.array([ 1, 1.2371, 1.6809, 2.89151, 5.13304, 9.23238, ]) y_data = np.array([ 0.0688824, 0.0490012, 0.0332843, 0.0235889, […]

Cómo generar valores de interpolación equidistantes.

Tengo una lista de valores (x, y) que no están espaciados uniformemente. Aquí está el archivo utilizado en esta pregunta. Soy capaz de interpolar entre los valores, pero lo que obtengo no son puntos de interpolación separados por puntos. Esto es lo que hago: x_data = [0.613,0.615,0.615,…] y_data = [5.919,5.349,5.413,…] # Interpolate values for x […]

¿Qué debo hacer si quiero una interpolación 3D spline / suave de datos no estructurados aleatorios?

Me inspiré en esta respuesta de @James para ver cómo se pueden usar griddata y map_coordinates . En los ejemplos a continuación, muestro datos en 2D, pero mi interés es en 3D. Noté que griddata solo proporciona splines para 1D y 2D, y está limitado a la interpolación lineal para 3D y superior (probablemente por […]