Кто-нибудь знает о замене Python для функции Matlab / Octave bwdist ()? Эта функция возвращает эвклидово расстояние каждой ячейки до ближайшей ненулевой ячейки для данной матрицы. Я видел реализацию Octave C, чистую реализацию Matlab, и мне было интересно, нужно ли кому-либо реализовать это в ANSI C (который не содержит заголовков Matlab / Octave, поэтому я […]
У меня есть входной файл, который содержит числа с плавающей запятой до 4 знаков после запятой: ie 13359 0.0000 0.0000 0.0001 0.0001 0.0002` 0.0003 0.0007 … (первым является id). Мой класс использует метод loadVectorsFromFile который умножает его на 10000, а затем int() эти числа. Кроме того, я также просматриваю каждый вектор, чтобы убедиться, что внутри […]
Я пытаюсь понять 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, поэтому я попытаюсь объяснить свою проблему: Я пытаюсь сгруппировать документы (иерархическая агломеративная кластеризация) и создал вектор для каждого документа и создал симметричную матрицу расстояний. Vector_list содержит (очень длинные) векторы, представляющие каждый документ. Порядок этого списка векторов совпадает с моим списком входных […]
У меня есть горизонтальная полоса 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 […]
У меня есть изображение, которое я размываю и растягиваю так: 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 -> 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 […]
Я пытаюсь преобразовать следующий код MATLAB в Python, и у меня возникли проблемы с поиском решения, которое работает в любое разумное время. M = diag(sum(a)) – a; where = vertcat(in, out); M(where,:) = 0; M(where,where) = 1; Здесь a – разреженная матрица и где – вектор (как и в / out). Решение, которое я использую […]
Я импортирую: from scipy import misc, io Но я получаю эти ошибки: Traceback (most recent call last): File "C:\work_asaaki\code\generateProposals.py", line 20, in <module> from scipy import misc, io File "C:\Python27\lib\site-packages\scipy\misc\__init__.py", line 47, in <module> from scipy.special import comb, factorial, factorial2, factorialk File "C:\Python27\lib\site-packages\scipy\special\__init__.py", line 548, in <module> from .basic import * File "C:\Python27\lib\site-packages\scipy\special\basic.py", line 17, […]
Может ли кто-нибудь представить пример предоставления функции якобиана функции наименьших квадратов в scipy ? Я не могу понять, какой метод они хотят – они говорят, что это должна быть функция, но очень сложно понять, какие входные параметры в каком порядке должна принимать эта функция.