Tag: scipy

Значение по умолчанию scipy.sparse

В разреженном формате матрицы (dok) предполагается, что значения ключей, не содержащихся в словаре, равны нулю. Есть ли способ заставить его использовать значение по умолчанию, отличное от нуля? Кроме того, существует ли способ вычисления логарифма разреженной матрицы (аналогичной np.log в регулярной матрице numpy)

Ограничение значений для curve_fit (scipy.optimize)

Я пытаюсь подогнать кривую логистического роста к моим данным, используя curve_fit, используя в качестве входной функции следующую функцию. def logistic(x, y0, k, d, a, b): if b > 0 and a > 0: y = (k * pow(1 + np.exp(d – (a * b * x) ), (-1/b) )) + y0 elif b >= -1 […]

Matlab / Octave bwdist () в Python или C

Кто-нибудь знает о замене Python для функции Matlab / Octave bwdist ()? Эта функция возвращает эвклидово расстояние каждой ячейки до ближайшей ненулевой ячейки для данной матрицы. Я видел реализацию Octave C, чистую реализацию Matlab, и мне было интересно, нужно ли кому-либо реализовать это в ANSI C (который не содержит заголовков Matlab / Octave, поэтому я […]

Scipy Отрицательное расстояние? Какие?

У меня есть входной файл, который содержит числа с плавающей запятой до 4 знаков после запятой: ie 13359 0.0000 0.0000 0.0001 0.0001 0.0002` 0.0003 0.0007 … (первым является id). Мой класс использует метод loadVectorsFromFile который умножает его на 10000, а затем int() эти числа. Кроме того, я также просматриваю каждый вектор, чтобы убедиться, что внутри […]

Как использовать scipy.weave.inline вместе с внешними библиотеками C?

Я пытаюсь понять weave.inline, чтобы обернуть C-код в моих программах Python. В приведенном ниже коде просто берется массив Numpy и умножается на все его элементы на 2. inl.py import numpy import scipy.weave a = numpy.array([1.0, 2.0, 3.0]) N = a.shape[0] print a code = \ """ int i; for(i = 0; i < N; i++) […]

Кластеризация с scipy-кластерами через матрицу расстояний, как вернуть исходные объекты

Я не могу шов, чтобы найти достаточно простые учебники или описания кластеризации в scipy, поэтому я попытаюсь объяснить свою проблему: Я пытаюсь сгруппировать документы (иерархическая агломеративная кластеризация) и создал вектор для каждого документа и создал симметричную матрицу расстояний. Vector_list содержит (очень длинные) векторы, представляющие каждый документ. Порядок этого списка векторов совпадает с моим списком входных […]

настройка только одной стороны меток в matplotlib с использованием шипов

У меня есть горизонтальная полоса matplotlib: import matplotlib.pyplot as plt from numpy import * from scipy import * bars = arange(5) + 0.1 vals = rand(5) print bars, vals plt.figure(figsize=(5,5), dpi=100) spines = ["bottom"] ax = plt.subplot(1, 1, 1) for loc, spine in ax.spines.iteritems(): if loc not in spines: spine.set_color('none') # don't draw ytick marks […]

Как удалить расширение на blob, вызванное морфологией

У меня есть изображение, которое я размываю и растягиваю так: kernel = np.ones((5,5),np.float32)/1 eroded_img = cv2.erode(self.inpainted_adjusted_image, kernel, iterations=10) dilated_img = cv2.dilate(eroded_img, kernel, iterations=10) Вот результат эрозии и дилатации: и тогда я беру порог его так: self.thresh = cv2.threshold(dilated_img, 0, 255, cv2.THRESH_BINARY + cv2.THRESH_OTSU)[1] Но порог дает мне нежелательное расширение, которое я отметил на изображении ниже […]

Эрозия изображения и расширение с помощью Scipy

Я пытаюсь использовать scipy, чтобы сделать эрозию и расширение изображения. Это выглядит довольно просто, используя scipy -> binary_erosion / dialation . Однако выход не является тем, что ожидается. Вот мой базовый код: import scipy from scipy import ndimage import matplotlib.pyplot as plt import numpy as np import Image #im = Image.open('flower.png') im = ndimage.imread('flower.png') im […]

Каков наиболее эффективный способ установки строк в нули для разреженной матрицы scipy?

Я пытаюсь преобразовать следующий код MATLAB в Python, и у меня возникли проблемы с поиском решения, которое работает в любое разумное время. M = diag(sum(a)) – a; where = vertcat(in, out); M(where,:) = 0; M(where,where) = 1; Здесь a – разреженная матрица и где – вектор (как и в / out). Решение, которое я использую […]

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