Articles of scipy

Producto puntual de un vector en SciPy / NumPy (obteniendo ValueError: los objetos no están alineados)

Acabo de empezar a aprender SciPy y estoy luchando con las funciones más básicas. Considere el siguiente vector estándar: In [6]: W=array([[1],[2]]) In [7]: print W [[1] [2]] Si lo comprendo correctamente, esta debería ser la representación en SciPy de un vector matemático estándar de 2×1, como este: (1) (2) El producto puntual de este […]

Consumo de memoria de la función NumPy por desviación estándar.

Actualmente estoy usando los enlaces Python de GDAL para trabajar en conjuntos de datos raster bastante grandes (> 4 GB). Ya que cargarlos en la memoria a la vez no es una solución viable para mí, los leo en bloques más pequeños y hago los cálculos pieza por pieza. Para evitar una nueva asignación para […]

¿Cómo visualizar la triangulación de delaunay 3D en Python?

Tengo un conjunto de puntos 3D que he usado scipy.spatial.Delaunay para hacer la triangulación / tetraédrica. Ahora tengo un conjunto de caras únicas de todos los tetraedros, y me gustaría visualizarlas en 3D. ¿Hay alguna biblioteca de Python (o bibliotecas con un envoltorio de Python) que puedan hacer esto?

Log-computaciones en Python

Estoy buscando para calcular algo como: Donde f(i) es una función que devuelve un número real en [-1,1] para cualquier i en {1,2,…,5000} . Obviamente, el resultado de la sum está en algún lugar en [-1,1] , pero cuando parece que no puedo calcularlo en Python usando encoding directa, como 0.5 5000 convierte en 0 […]

El ajuste de la curva falla con exponencial pero zunzun lo hace bien

Estoy tratando de calcular el mejor ajuste de dos formas de un exponencial a algunos datos x, y (el archivo de datos se puede descargar desde aquí ) Aquí está el código: from scipy.optimize import curve_fit import numpy as np # Get x,y data data = np.loadtxt(‘data.txt’, unpack=True) xdata, ydata = data[0], data[1] # Define […]

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

scipy curve_fit falla en ajuste exponencial

Cuando trato de hacer un ajuste exponencial usando curve_fit, scipy devuelve un error. ¿Estoy haciendo algo mal? Eliminar el signo negativo de np.exp (-b * t) permite que curve_fit funcione, pero los valores que devuelve están muy lejos. #!/usr/bin/python import numpy as np import scipy as sp from scipy.optimize import curve_fit import scipy.optimize as opt […]

Dendrogtwig u otro gráfico de la matriz de distancia

Tengo tres matrices para comparar. Cada uno de ellos son 5×6. Originalmente quería usar agrupamiento jerárquico para agrupar las matrices, de modo que las matrices más similares se agrupen, dado un umbral de similitud. No pude encontrar ninguna de estas funciones en python, así que implementé la medida de distancia a mano, (p-norma donde p […]

Definir evaluación numérica de una derivada de una función sympy.

¿Cómo puedo definir la evaluación numérica de una derivada de una función en sympy? Tengo algunas funciones que puedo describir con splines para la función y es derivativo usando scipy.interpolate. Quiero manipular algunas expresiones con esta función y luego evaluar las expresiones con las splines. Puedo usar lambdify para hacer que una función sympy se […]

¿Cómo leer / recorrer / cortar matrices dispersas Scipy (LIL, CSR, COO, DOK) más rápido?

Para manipular matrices Scipy, normalmente, se utilizan los métodos incorporados. Pero a veces es necesario leer los datos de la matriz para asignarlos a tipos de datos no dispersos. Por el bien de la demostración, creé una matriz dispersa aleatoria de LIL y la convertí a una matriz Numpy (¡los tipos de datos de python […]