Tag: numpy

Нежелательные дополнительные размеры в массиве numpy

Я открыл изображение .fits: scaled_flat1 = pyfits.open('scaled_flat1.fit') scaled_flat1a = scaled_flat1[0].data и когда я печатаю его форму: print scaled_flat1a.shape Я получаю следующее: (1, 1, 510, 765) Я хочу, чтобы это прочитало: (510,765) Как мне избавиться от этих двух?

Сплайны с Python (с использованием узлов управления и конечных точек)

Я пытаюсь сделать что-то вроде следующего (изображение, извлеченное из википедии) #!/usr/bin/env python from scipy import interpolate import numpy as np import matplotlib.pyplot as plt # sampling x = np.linspace(0, 10, 10) y = np.sin(x) # spline trough all the sampled points tck = interpolate.splrep(x, y) x2 = np.linspace(0, 10, 200) y2 = interpolate.splev(x2, tck) # […]

Как получить матричное умножение по типу элемента (произведение Адамара) в numpy?

У меня две матрицы a = [[1,2], [3,4]] b = [[5,6], [7,8]] Я хочу получить элементарную продеку, которая будет [[1*5,2*6], [3*7,4*8]] равна [[5,12], [21,32]] Я пробовал с numpy print(np.dot(x,y)) а также print(x*y) результат: [[19 22], [43 50]] который дает мне результат умножения матрицы, а не элементный продукт. Как я могу получить элементарный продукт. Продукт Hadamard […]

Подмножество pandas dataframe

У меня есть DataFrame, загруженный из файла tsv. Я хотел создать несколько поисковых участков. Проблема в том, что набор данных большой (~ 1 миллион строк), поэтому на графике слишком много точек, чтобы увидеть тренд. Кроме того, на графику требуется некоторое время. Я хотел подмножить 10000 случайно распределенных строк. Кроме того, это должно быть воспроизводимым, поэтому […]

Pythonic и эффективный способ сделать элемент «in» с помощью numpy

Я ищу способ эффективного получения массива булевых, где заданы два массива с равными размерами a и b , каждый элемент является истинным, если соответствующий элемент a появляется in соответствующем элементе b . Например, следующая программа: a = numpy.array([1, 2, 3, 4]) b = numpy.array([[1, 2, 13], [2, 8, 9], [5, 6], [7]]) print(numpy.magic_function(a, b)) Должен […]

Написание более быстрого симулятора физики Python

Я играю с написанием своего собственного движка физики в Python как упражнение в физике и программировании. Я начал, следуя учебнику, расположенному здесь . Это получилось хорошо, но затем я нашел статью «Продвинутая физика персонажей» от thomas jakobsen, которая охватила интеграцию Верле для симуляций, которую я нашел захватывающей. Я пытаюсь написать свой собственный базовый физический симулятор, […]

Методы ввода уравнений при программировании в C / C ++, Python или Fortran

Я пишу код, который имел длинные математические уравнения со многими тригонометрическими и другими тождествами. Есть ли способ визуализации одного и того же выражения в латексе и создания выражения C или python из него или наоборот. Как ввести и проверить математические выражения, чтобы увидеть, находятся ли скобки и т. Д. В правильном положении и использовать их […]

Сортировка строк массива другим массивом в python

Я пытаюсь сортировать строки одного массива значениями другого. Например: import numpy as np arr1 = np.random.normal(1, 1, 80) arr2 = np.random.normal(1,1, (80,100)) Я хочу сортировать arr1 в порядке убывания и поддерживать текущую связь между arr1 и arr2 (т. Е. После сортировки обеих строк строки arr1[0] и arr2[0, :] совпадают).

Найти все n-мерные линии и диагонали с NumPy

Используя NumPy, я хотел бы создать список всех строк и диагонали n-мерного массива с длиной k. Возьмем случай следующего трехмерного массива длиной 3. array([[[ 0, 1, 2], [ 3, 4, 5], [ 6, 7, 8]], [[ 9, 10, 11], [12, 13, 14], [15, 16, 17]], [[18, 19, 20], [21, 22, 23], [24, 25, 26]]]) Для […]

Numpy: Почему не работает «a + = aT»?

Как сказано в лекционных лекциях, это не будет работать так, как ожидалось: a = np.random.randint(0, 10, (1000, 1000)) a += aT assert np.allclose(a, aT) Но почему? Как восприятие влияет на это поведение?

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