Tag: scipy

Как получить значения из массива NumPy с использованием нескольких индексов

У меня есть массив NumPy, который выглядит так: arr = [100.10, 200.42, 4.14, 89.00, 34.55, 1.12] Как я могу получить несколько значений из этого массива по индексу: Например, как я могу получить значения в позициях индекса 1, 4 и 5? Я пытался что-то вроде этого, что неверно: arr[1, 4, 5]

Scipy – 3d griddata – Почему нужно вводить аргумент griddata xi для кортежа?

Почему нижний вызов griddata терпит неудачу? import scipy.interpolate import numpy as np grid_vals = np.meshgrid(*([np.linspace(-1,1,200)] * 3)) interp_vals = scipy.interpolate.griddata(np.random.randn(50,3), np.random.randn(50), grid_vals, 'linear') Возникает следующее исключение: ValueError: количество измерений в xi не соответствует x Если я передал аргумент xi (grid_vals) для кортежа: interp_vals = scipy.interpolate.griddata(np.random.randn(50,3), np.random.randn(50), tuple(grid_vals), 'linear') Ошибка исчезнет. Зачем?

правильное использование scipy.optimize.fmin_bfgs

Я играю с логистической регрессией в Python. Я реализовал версию, в которой минимизация функции стоимости выполняется с помощью градиентного спуска, и теперь я бы хотел использовать алгоритм BFGS из scipy ( scipy.optimize.fmin_bfgs ). У меня есть набор данных (функции в матрице X с одним образцом в каждой строке X и корреляционные метки в вертикальном векторе […]

Использование pytables, которое более эффективно: scipy.sparse или numpy плотная матрица?

При использовании pytables нет поддержки (насколько я могу судить) для матричных форматов scipy.sparse , поэтому для хранения матрицы мне нужно сделать некоторое преобразование, например def store_sparse_matrix(self): grp1 = self.getFileHandle().createGroup(self.getGroup(), 'M') self.getFileHandle().createArray(grp1, 'data', M.tocsr().data) self.getFileHandle().createArray(grp1, 'indptr', M.tocsr().indptr) self.getFileHandle().createArray(grp1, 'indices', M.tocsr().indices) def get_sparse_matrix(self): return sparse.csr_matrix((self.getGroup().M.data, self.getGroup().M.indices, self.getGroup().M.indptr)) Проблема в том, что функция get_sparse занимает некоторое время (чтение […]

иерархическая кластеризация на корреляциях в Python scipy / numpy?

Как запустить иерархическую кластеризацию на корреляционной матрице в scipy / numpy? У меня есть матрица из 100 строк по 9 столбцов, и я хотел бы иерархически кластеризовать по корреляциям каждой записи в 9 условиях. Я бы хотел использовать корреляцию 1-pearson как расстояния для кластеризации. Предполагая, что у меня есть массив numpy «X», который содержит матрицу […]

Как я могу получить быструю оценку расстояния между точкой и бикубической поверхностью сплайна в Python?

Как я могу получить быструю оценку расстояния между точкой и бикубической поверхностью сплайна в Python? Есть ли существующее решение, которое я мог бы использовать в SciPy, NumPy или другом пакете? У меня есть поверхность, определяемая бикубической интерполяцией: import numpy as np import scipy.interpolate # Define regular grid surface xmin,xmax,ymin,ymax = 25, 125, -50, 50 x […]

Открытое квантовое моделирование системы

Я долгое время работал над моделированием открытой квантовой системы с использованием уравнения Линдблада . Гамильтониан следующий: Однако к гамильтониану добавляются две другие матрицы. Один из них имеет все диагональные члены, равные -33.3333i, и все остальное ноль. Другая – матрица с третьим диагональным членом, равным -0.033333i. Уравнение Линдблада таково: где L_i – матрицы (в списке: [L1, […]

Ускорение линейной интерполяции многих местоположений пикселей в NumPy

Я попытался воспроизвести основное узкое место в одной из моих программ. Я хочу получить линейно (или, скорее, билинейно) интерполированные значения нескольких нецелых значений пикселей одновременно. Это не так, что каждая координата пикселя возмущена одинаково. Ниже приведен полный / минимальный сценарий вместе с комментариями, которые демонстрируют проблему. Как ускорить вычисление result ? import numpy as np […]

Изменение размера изображения SciPy – ожидаемое поведение или ошибка?

Я заметил что-то странное с scipy.misc.resize – кажется, что любой другой метод интерполяции, отличный от «ближайших» результатов, смещается примерно на 1×1 пиксель от (0,0) в полученном изображении. Вот полностью синтетический пример взятия 3×3 изображения в 6×6: >>> src array([[ 0., 0., 0.], [ 0., 64., 0.], [ 0., 0., 0.]]) >>> imresize(src, (6, 6), interp='bicubic',mode='F') […]

Вычисление расстояний между уникальными областями массива Python?

У меня есть растр с набором уникальных патчей / регионов ID, которые я преобразовал в двумерный массив numpy Python. Я хотел бы вычислить парные евклидовы расстояния между всеми областями, чтобы получить минимальное расстояние, разделяющее ближайшие края каждого растрового патча. Поскольку массив был первоначально растром, решение должно учитывать диагональные расстояния по ячейкам (я всегда могу преобразовать […]

 
Interesting Posts for Van-Lav

Как генерировать случайные числа, которые отличаются?

Python ORM, который автоматически генерирует / обновляет таблицы и использует SQLite?

Перемещение неперекрывающегося окна в Numpy

Исключения модуля sqlite3 на Python – где находится документация?

сумма матричных столбцов в python

Appengine – возможно ли получить сущность, использующую только ключевую строку без имени модели?

Как реализовать аутентификацию для REST API?

Импортированный модуль, но все же необходимо использовать полное имя

Запись в CSV, получение «Ошибка: необходимость выхода» для пустой строки

Использование True / False в качестве ключей – как / почему это работает?

Загрузка pysftp не удалась из-за того, что клиент превысил внутренние буферы сервера

Как добавить врагов в pygame?

Как передать загруженное изображение в template.html в Flask

Автозаполнение в динамических языковых IDE, в частности Python в PyDev

Чтение Python Wiegand drop of Zeros

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