Tag: scipy

Оцените стандартное отклонение установленных параметров в scipy.odr?

( Несколько связанный с этим вопросом Linear fit, включая все ошибки с NumPy / SciPy и код заимствования из этого одного Linear фитинга в python с неопределенностью как в координатах x, так и y ) Я подхожу к линейной модели ( y=a*x+b ), используя фиксированные ошибки в x,y используя scipy.odr (код ниже), и я получаю: […]

Numpy: Row Wise Уникальные элементы

Кто-нибудь знает, как получить уникальные элементы в мудрую матрицу. Например, входная матрица может быть такой: a = [[1,2,1,3,4,1,3], [5,5,3,1,5,1,2], [1,2,3,4,5,6,7], [9,3,8,2,9,8,4], [4,6,7,4,2,3,5]] Он должен вернуть следующее: b = rowWiseUnique(a) => b = [[1,2,3,4,0,0,0], [5,3,1,2,0,0,0], [1,2,3,4,5,6,7], [9,3,8,2,4,0,0], [4,6,7,2,3,5,0]] Каков наиболее эффективный способ сделать это в numpy? Я попробовал следующий код, есть ли лучший и более короткий […]

Удалите все нулевые строки в 2D-матрице

Есть ли эффективная и / или встроенная функция для удаления всех нулевых строк массива 2d? Я смотрю на документацию numpy, но я ее не нашел.

Полосовой фильтр в питоне

Я пытаюсь получить полосовой фильтр с 128-точечным окном Хэмминга с частотой среза 0,7-4 Гц в питоне. Я получаю свои образцы для моего сигнала из изображений. (1 образец = 1 изображение). Часто меняется частота кадров. Как это можно сделать в python? Я читал это: http://mpastell.com/2010/01/18/fir-with-scipy/, но я считаю, что firwin довольно запутанным. Как это можно сделать […]

Привязка векторной функции с кривой_fit в Scipy

Я хочу, чтобы функция с векторным выводом использовалась с использованием кривой curve_fit (или что-то более подходящее, если доступно). Например, рассмотрим следующую функцию: import numpy as np def fmodel(x, a, b): return np.vstack([a*np.sin(b*x), a*x**2 – b*x, a*np.exp(b/x)]) Каждый компонент является другой функцией, но они разделяют параметры, которые я хочу поместить. В идеале я бы сделал что-то […]

Как работает функция scipy distance_transform_edt?

https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.ndimage.morphology.distance_transform_edt.html Мне трудно понять, как работает функция евклидовой дистанционной трансляции в Scipy. Из того, что я понимаю, он отличается от функции Matlab (bwdist). В качестве примера для ввода: [[ 0. 0. 0. 0. 0.] [ 0. 1. 0. 0. 0.] [ 0. 0. 0. 0. 0.] [ 0. 0. 0. 1. 0.] [ 0. 0. […]

Эффективное сопоставление двух массивов (как использовать KDTree)

У меня два 2d массива, obs1 и obs2 . Они представляют собой две независимые серии измерений, и оба имеют dim0 = 2 и немного разные dim1 , скажем, obs1.shape = (2, 250000) и obs2.shape = (2, 250050) . obs1[0] и obs2[0] означают время, а obs1[1] и obs2[1] означают некоторую пространственную координату. Оба массива (более или […]

найти средние значения бинов, используя python histogram2d

Как вы вычисляете средние значения для бинов с 2D-гистограммой в python? У меня есть диапазоны температур для оси x и y, и я пытаюсь построить вероятность молнии с использованием бункеров для соответствующих температур. Я читаю данные из файла csv, и мой код таков: filename = 'Random_Events_All_Sorted_85GHz.csv' df = pd.read_csv(filename) min37 = df.min37 min85 = df.min85 […]

python scipy.stats.powerlaw отрицательный показатель

Я хочу предоставить отрицательный показатель для процедуры scipy.stats.powerlaw, например a = -1.5, для того, чтобы нарисовать случайные выборки: """ powerlaw.pdf(x, a) = a * x**(a-1) """ from scipy.stats import powerlaw R = powerlaw.rvs(a, size=100) Почему требуется> 0, как я могу предоставить отрицательный а для генерации случайных выборок и как я могу обеспечить коэффициент / преобразование […]

Интерполяция 3d-массива в Python. Как избежать циклов?

У меня есть массив, который я хочу интерполировать по 1-й оси. На данный момент я делаю это, как этот пример: import numpy as np from scipy.interpolate import interp1d array = np.random.randint(0, 9, size=(100, 100, 100)) new_array = np.zeros((1000, 100, 100)) x = np.arange(0, 100, 1) x_new = np.arange(0, 100, 0.1) for i in x: for […]

Python - лучший язык программирования в мире.