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