Articles of Interpolation

Uso correcto de scipy.interpolate.RegularGridInterpolator

Estoy un poco confundido por la documentación de scipy.interpolate.RegularGridInterpolator . Por ejemplo, tengo una función f: R ^ 3 => R que se muestrea en los vértices del cubo unitario. Me gustaría interpolar para encontrar valores dentro del cubo. import numpy as np # Grid points / sample locations X = np.array([[0,0,0], [0,0,1], [0,1,0], [0,1,1], […]

Scipy map_coordinates interpolación bilineal en comparación con interpolación interpola y IDL

Estoy en el proceso de reescribir un código IDL de compañeros de trabajo en python y estoy encontrando algunas diferencias con las que estoy confundido. De acuerdo con otras preguntas de SO y los hilos de la lista de correo que he encontrado si usa scipy.ndimage.interpolation.map_coordinates y especifica order=1 , se supone que debe hacer […]

Cómo mejorar el rendimiento al interpolar en datos 3D con SciPy

Tengo datos 3D que representan la atmósfera. Ahora quiero interpolar estos datos a una coordenada Z común (lo que quiero decir con esto debería quedar claro en la doctrine de la función). El siguiente código funciona bien, pero me preguntaba si habría alguna manera de mejorar el rendimiento … def interpLevel(grid,value,data,interp=’linear’): “”” Interpolate 3d data […]

Problemas con la interpolación 2D en Scipy

En mi aplicación, los datos de datos se muestrean en una cuadrícula distorsionada, y me gustaría volver a muestrearla en una cuadrícula no distorsionada. Para probar esto, escribí este progtwig con distorsiones ejemplares y una función simple como datos: from __future__ import division import numpy as np import scipy.interpolate as intp import pylab as plt […]

Uso de funciones de base radial para interpolar una función en una esfera

Primero, un poco de fondo: Estoy usando armónicos esféricos como ejemplo de una función en la superficie de una esfera como las esferas delanteras en esta imagen: Produje una de estas esferas, coloreadas según el valor de la función armónica en los puntos de su superficie. Primero hago esto para una gran cantidad de puntos, […]

¿Cómo realizar la interpolación spline cúbica en python?

Tengo dos listas para describir la función y (x): x = [0,1,2,3,4,5] y = [12,14,22,39,58,77] Me gustaría realizar una interpolación spline cúbica para que, dado un cierto valor u en el dominio de x, por ejemplo u = 1.25 Puedo encontrarte (u). Encontré esto en SciPy pero no estoy seguro de cómo usarlo.

Interpolación Python / Scipy (map_coordinates)

Estoy tratando de hacer alguna interpolación con scipy. He pasado por muchos ejemplos, pero no encuentro exactamente lo que quiero. Digamos que tengo algunos datos donde la variable de fila y columna puede variar de 0 a 1. El delta cambia entre cada fila y columna no siempre es la misma (ver más abajo). | […]

interpolar a tiempo específico

Digamos que tengo este código: import numpy as np import time from datetime import datetime class Measurements(): def __init__(self, time_var, value): self.time_var = time_var self.value = value a = np.array([ Measurements(’30-01-2017 12:02:15.880922′, 100), Measurements(’30-01-2017 12:02:16.880922′, 100), Measurements(’30-01-2017 12:02:17.880922′, 110), Measurements(’30-01-2017 12:02:18.880922′, 99), Measurements(’30-01-2017 12:02:19.880922′, 96)]) b = np.array([ Measurements(’30-01-2017 12:02:15.123444′, 10), Measurements(’30-01-2017 12:02:18.880919′, 12), ]) […]

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