Articles of arrays

Python numpy o pandas equivalente de la función R de barrido ()

¿Cuál es el equivalente numpy o pandas del sweep() de la función R sweep() ? Para elaborar: en R digamos que tenemos un vector de coeficiente (por ejemplo, tipo beta – numérico) y una matriz (por ejemplo, datos – tipo numérico 20×5). Quiero superponer el vector en cada fila de la matriz y multiplicar los […]

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

Python, utiliza el multiprocesamiento para acelerar aún más una función de cython

el código que se muestra aquí se simplifica pero activa el mismo PicklingError. Sé que hay mucha discusión sobre lo que puede y no puede ser decapado, pero encontré la solución de ellos. Escribo un script de cython simple con la siguiente función: def pow2(int a) : return a**2 La comstackción está funcionando, puedo llamar […]

¿Guardar un dict con np.savez da un resultado inesperado?

¿Puedo almacenar un diccionario usando np.savez? Los resultados son sorprendentes (al menos para mí) y no puedo encontrar una manera de recuperar mis datos por clave. In [1]: a = {‘0’: {‘A’: array([1,2,3]), ‘B’: array([4,5,6])}} In [2]: a Out[2]: {‘0’: {‘A’: array([1, 2, 3]), ‘B’: array([4, 5, 6])}} In [3]: np.savez(‘model.npz’, **a) In [4]: a […]

crear pequeños arreglos en cython toma una enorme cantidad de tiempo

Estaba escribiendo un nuevo generador de números aleatorios para números que produce números aleatorios según una distribución arbitraria cuando me encontré con este comportamiento realmente extraño: esto es test.pyx #cython: boundscheck=False #cython: wraparound=False import numpy as np cimport numpy as np cimport cython def BareBones(np.ndarray[double, ndim=1] a,np.ndarray[double, ndim=1] u,r): return u def UntypedWithLoop(a,u,r): cdef int […]

Función de Python que maneja escalar o arrays

¿Cuál es la mejor forma de escribir una función que pueda aceptar ya sea flotantes escalares o vectores numpy (matriz 1-d), y devolver una matriz escalar, 1-d o 2-d, según la entrada? La función es costosa y se llama a menudo, y no quiero poner una carga en la persona que llama para hacer conversiones […]

Mapa de color de Python pero con todos los valores cero asignados a negro

Tengo una matriz cuadrada de elementos que corresponden a los sitios de celosía. Algunos de los elementos son cero y el rest varía entre 1 y aproximadamente 2700. Con imshow y el mapa de color OrRd, quiero que todos los sitios de celosía mayores a 0 muestren el color correspondiente, pero lo más importante es […]

Función vectorizada de pandas cumsum versus numpy

Al responder a la pregunta Vectorizar el cálculo de un Marco de datos de Pandas , noté un problema interesante con respecto al rendimiento. Tenía la impresión de que funciones df.min() como df.min() , df.mean() , df.cumsum() , etc. Sin embargo, estoy viendo una gran discrepancia en el rendimiento entre df.cumsum() y una alternativa numpy […]

Cargue datos comprimidos (.npz) desde el archivo usando numpy.load

Tengo una matriz: >>> data = np.ones((1,3,128)) Lo savez_compressed en un archivo usando savez_compressed : >>> with open(‘afile’,’w’) as f: np.savez_compressed(f,data=data) Cuando bash cargarlo, parece que no puedo acceder a los datos: >>> with open(‘afile’,’r’) as f: b=np.load(f) >>> b.files [‘data’] >>> b[‘data’] Traceback (most recent call last): File “”, line 1, in b[‘data’] File […]

¿Cómo crear una matriz numpy desde un AudioSegment pydub?

Soy consciente de la siguiente pregunta: ¿Cómo crear un AudioSegment pydub usando una matriz numpy? Mi pregunta es lo contrario. Si tengo un Pydub AudioSegment, ¿cómo puedo convertirlo en una matriz numpy? Me gustaría usar filtros scipy y así sucesivamente. No me queda muy claro cuál es la estructura interna de los datos en bruto […]