Articles of Interpolation

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

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