Articles of interpolación

Cómo realizar la interpolación bilineal en Python

Me gustaría realizar la interpolación blinear usando python. Ejemplo de punto gps para el que quiero interpolar altura es: B = 54.4786674627 L = 17.0470721369 utilizando cuatro puntos adyacentes con coordenadas conocidas y valores de altura: n = [(54.5, 17.041667, 31.993), (54.5, 17.083333, 31.911), (54.458333, 17.041667, 31.945), (54.458333, 17.083333, 31.866)] z01 z11 z z00 z10 […]

¿Cómo puedo obtener la salida de un gráfico de matplotlib como SVG?

Necesito tomar la salida de un gráfico de matplotlib y convertirlo en una ruta SVG que pueda usar en un cortador láser. import matplotlib.pyplot as plt import numpy as np x = np.arange(0,100,0.00001) y = x*np.sin(2*pi*x) plt.plot(y) plt.show() Por ejemplo, a continuación ves una forma de onda. Me gustaría poder emitir o guardar esta forma […]

Tamaño de voxel de escaneo ct de unificación mediante interpolación en Python

He utilizado interp2 en Matlab, como el siguiente código, que forma parte de la respuesta de @ rayryeng en: Interpolación de matriz tridimensional (3D) en Matlab : d = size(volume_image) [X,Y] = meshgrid(1:1/scaleCoeff(2):d(2), 1:1/scaleCoeff(1):d(1)); for ind = z %Interpolate each slice via interp2 M2D(:,:,ind) = interp2(volume_image(:,:,ind), X, Y); end Ejemplo de dimensiones: The image size […]

Volver a muestrear series de tiempo en pandas a un intervalo semanal

¿Cómo vuelvo a muestrear una serie de tiempo en pandas a una frecuencia semanal donde las semanas comienzan en un día arbitrario? Veo que hay una base de palabras clave opcional, pero solo funciona para intervalos más cortos que un día.

interpolar valores de un dataframe basado en un valor de columna

Suponiendo que tengo el siguiente problema: import pandas as pd import numpy as np xp = [0.0, 0.5, 1.0] np.random.seed(100) df = pd.DataFrame(np.random.rand(10, 4), columns=[‘x0’, ‘y1’, ‘y2’, ‘y3’]) df x0 y1 y2 y3 0 0.5434 0.2784 0.4245 0.8448 1 0.0047 0.1216 0.6707 0.8259 2 0.1367 0.5751 0.8913 0.2092 3 0.1853 0.1084 0.2197 0.9786 4 0.8117 […]

Speedup scipy griddata para múltiples interpolaciones entre dos grillas irregulares

Tengo varios valores que están definidos en la misma cuadrícula irregular (x, y, z) que quiero interpolar en una nueva cuadrícula (x1, y1, z1) . es decir, tengo f(x, y, z), g(x, y, z), h(x, y, z) y quiero calcular f(x1, y1, z1), g(x1, y1, z1), h(x1, y1, z1) . En este momento estoy haciendo […]

Interpolar los valores de NaN en una matriz numpy

¿Hay una manera rápida de reemplazar todos los valores de NaN en una matriz de números por (digamos) los valores interpolados linealmente? Por ejemplo, [1 1 1 nan nan 2 2 nan 0] sería convertido en [1 1 1 1.3 1.6 2 2 1 0]

interpolar el volumen 3D con adormilado y o scipy

Estoy extremadamente frustrado porque después de varias horas parece que no puedo hacer una interpolación 3D aparentemente fácil en python. En Matlab todo lo que tenía que hacer era Vi = interp3(x,y,z,V,xi,yi,zi) ¿Cuál es el equivalente exacto de esto usando ndimage.map_coordinate u otros métodos numpy de scipy‘s? Gracias

La línea suave con objetos spline + datetime no funciona

He estado tratando de hacer una ttwig más suave, como se hace aquí , pero mis X son objetos de fecha y hora que no son compatibles con linspace … Convierto las X a fechas de matplotlib: Xnew = matplotlib.dates.date2num(X) X_smooth = np.linspace(Xnew.min(), Xnew.max(), 10) Y_smooth = spline(Xnew, Y, X_smooth) Pero luego me sale una […]

pandas multiindex dataframe, interpolación ND para valores perdidos

¿Es posible en los pandas interpolar los valores faltantes en el dataframe multiindex? Este ejemplo a continuación no funciona como se esperaba: arr1=np.array(np.arange(1.,10.,1.)) arr2=np.array(np.arange(2.,20.,2.)) df1=pd.DataFrame(zip(arr1,arr2,arr1+arr2,arr1*arr2),columns=[‘x’,’y’,’xplusy’,’xtimesy’]) df1.set_index([‘x’,’y’],inplace=True) df2=df1.reindex(index=zip(*df1.index.levels)+[(2,2),(3,2),(5,5)]) df2.sortlevel([0,1],inplace=True) df2.interpolate(method=’linear’,inplace=True) no muestra lo que esperaba en las columnas xplusy y xtimesy para los índices agregados. ———– —- — (1.0, 2.0) 3 2 (2.0, 2.0) 4.5 5 […]