Articles of numpy

Python, por que eu = + 1 não resulta em um loop infinito?

Então no outro dia eu estava preso em um problema por causa de um erro de digitação. Em vez de percorrer o loop nested com i + = 1, usei i = + 1. Eu não percebi isso até que comecei a imprimir o número de etapas e vi que estava imprimindo a etapa 1 […]

Transmissão de elementos para comparar dois arrays NumPy?

Digamos que eu tenha uma matriz como esta: import numpy as np base_array = np.array([-13, -9, -11, -3, -3, -4, 2, 2, 2, 5, 7, 7, 8, 7, 12, 11]) Suponha que eu queira saber: “quantos elementos em base_array são maiores que 4?” Isso pode ser feito simplesmente explorando a transmissão: np.sum(4 < base_array) Para […]

Binarize dados de imagem

Eu tenho 10 scans de ressonância magnética em escala de cinzentos do BrainWeb . Eles são armazenados como um array numpy 4d, brains , com forma (10, 181, 217, 181) . Cada um dos 10 cérebros é composto de 181 fatias ao longo do plano z (passando pelo topo da cabeça até o pescoço), onde […]

Localização da documentação sobre methods especiais reconhecidos por numpy

Uma das diferenças entre math.exp e numpy.exp é que, se você tiver uma class C personalizada que tenha um método numpy.exp , numpy.exp notará e delegará a este método, enquanto math.exp não: class C: def exp(self): return ‘hey!’ import math math.exp(C()) # raises TypeError import numpy numpy.exp(C()) # evaluates to ‘hey!’ No entanto, se você […]

Encontre todos os pares de pontos mais próximos do que uma dada distância máxima

Eu quero encontrar (eficientemente) todos os pares de pontos que estão mais próximos do que alguma distância max_d . Meu método atual, usando o cdist , é: import numpy as np from scipy.spatial.distance import cdist def close_pairs(X,max_d): d = cdist(X,X) I,J = (d<max_d).nonzero() IJ = np.sort(np.vstack((I,J)), axis=0) # remove diagonal element IJ = IJ[:,np.diff(IJ,axis=0).ravel()0] # […]

Conversão de código Matlab para Python: as matrizes não estão alinhadas

Abaixo está um exemplo de código do MATLAB e seu código Python eqv usando o pacote Numpy. O código MATLAB funciona bem, mas o código Python está dando problemas: MATLAB / OCTAVE N=1200 YDFA_P0 = double([1;2;3;4;5]) P0=YDFA_P0 *ones(1, N) octave:27> whos P0 Variables in the current scope: Attr Name Size Bytes Class ==== ==== ==== […]

Python Localizar max na coluna dataframe para fazer um loop para encontrar todos os valores

Então eu tenho um grande dataframe, usando pandas. Quando faço max(df[‘A’]) ele reporta um máximo de 9999 quando deveria ser 396450 por observação. import numpy as numpy import pandas as pd f = open(“20170901.as-rel2.txt”, ‘r’) #read file into array, ignore first 6 lines lines = loadtxt(“20170901.as-rel2.txt”, dtype=’str’, comments=”#”, delimiter=”|”, unpack=False) #ignore col 4 lines=lines[:, :3] […]

interpolar para hora específica

Digamos que eu tenha este código: import numpy as np import time from datetime import datetime class Measurements(): def __init__(self, time_var, value): self.time_var = time_var self.value = value a = np.array([ Measurements(’30-01-2017 12:02:15.880922′, 100), Measurements(’30-01-2017 12:02:16.880922′, 100), Measurements(’30-01-2017 12:02:17.880922′, 110), Measurements(’30-01-2017 12:02:18.880922′, 99), Measurements(’30-01-2017 12:02:19.880922′, 96)]) b = np.array([ Measurements(’30-01-2017 12:02:15.123444′, 10), Measurements(’30-01-2017 12:02:18.880919′, 12), […]

Erro numpy / CAPI com import_array () ao compilar vários módulos

Eu estou tentando compilar um módulo de C ++ para usar em scipy.weave que é composto de vários headers e arquivos de origem C ++. Esses arquivos contêm classs e methods que usam extensivamente a interface Numpy / C-API. Mas eu não estou conseguindo descobrir como include import_array() sucesso. Eu tenho lutado com isso durante […]

Plotar eficientemente uma tabela no formato csv usando Python

Eu estou tentando plotar uma tabela formatada csv usando Python. Até agora, consegui o resultado que queria lendo perguntas semelhantes no site, mas a minha solução não parece muito “pythonic”, nem encontrei uma maneira muito simples de fazer isso. Tenho certeza de que há uma maneira mais eficiente de plotar uma tabela, por isso estou […]