【Python numpy / scipy特征向量似乎对马尔可夫链模型不正确】教程文章相关的互联网学习教程文章

python – PIL Image构造来自numpy数组的奇怪图像 – 为什么?【代码】

我想要一种方法来生成红色,绿色或蓝色的小RGB方形图像.它应该产生坚固的色块,但PIL输出的图像非常奇怪.为什么?import numpy as np from PIL import Imageclass MakeSquares():def __init__(self):self.num_rows = 3self.num_cols = 3self.colourmap = {'red': [255, 0, 0],'green': [0, 255, 0],'blue': [0, 0, 255]}def generateExample(self, label):arr = []colour = labelcolour_array = self.colourmap[colour]for i in rang...

python – numpy的百分位函数究竟做了什么?【代码】

根据我的理解,numpy’s percentile计算数据的第q百分位数. 但它是如何做到的呢? 比如说,给定x = np.array([1.3,1.7,2.4,2.8,3.5,5.6,6.6,7.7,8.8,9.9])(里面有10个浮点数). 如果我做np.percentile(x,100),它会返回9.9000000000000004. 如果我做np.percentile(x,90),它应该返回8.8,对吗?但它回馈了8.9100000000000001. 为什么会有这样的差异?这些差异是否可以接受?解决方法:从版本1.9.0开始,Numpy的百分位函数有一个插值参数,如...

python – 相当于熊猫的numpy【代码】

我有以下numpy数组:x = np.arange(9.).reshape(3, 3)从而:>>>> x array([[ 0., 1., 2.],[ 3., 4., 5.],[ 6., 7., 8.]])我想选择第三列大于2的所有行. 在熊猫中,我会这样做:import pandas as pd d = DataFrame(x)>>>> d[d.iloc[:,2]>2]输出将是:0 1 2 1 3 4 5 2 6 7 8如何使用numpy获得相同的输出?我尝试了numpy,但它返回索引,而不是值.谢谢!解决方法: In [120]: x[x[: , 2] > 2] Out[120]: array([[ 3., 4., ...

python – 使用Numpy手动反转FFT【代码】

我有一个用于计算方波的傅立叶变换的小脚本,当我使用numpy.fft.ifft()反转fft时,该脚本运行良好并正确返回方波.但是,我无法通过在将它们乘以我从numpy.fft.fft获得的各自系数后手动累加谐波来反转变换()下面是我的脚本,我相信你会看到我的意图.from numpy import zeros, concatenate, sin, pi, linspace from numpy.fft import fft, fftfreq, ifft import numpy as np import matplotlib.pyplot as pltN = 1024 # samples T = 1 # ...

python – numpy中的随机值【代码】

当条件为True时,我试图获得一系列随机选择的数字. 我有一个充满布尔数的NumPy数组:In [1]: multiples Out[1]: array([False, False, False, True, True, False, False, False, False,False, False, False, False, False, False, False, False, False,False, False, False, False, False, False, False, False, False,False, False, True, True, True, True, False, False, False,False, False, False, True, True, False,...

python – X * Y项的平均值并保持numpy数组的大小【代码】

你如何在二维阵列上取平均值例如四个附近的项目(2 * 2)?我的意见是:[[1,1,1,1],[1,1,0,0],[0,0,1,1],[0,0,0,0]]应该导致:[[1, 0.5],[0, 0.5]]numpy.mean(x.reshape(-1,4),1)将展平数组,并按错误的顺序略微平均四个项目. 附加信息 例如,通过此方法生成数组:n = 10 l = 100 A = np.zeros((l, l)) points = l*np.random.random((2, n**2)) A[(points[0]).astype(np.int), (points[1]).astype(np.int)] = 1 A = ndimage.gaussian_f...

如何打印列表每次在python中跳过一个元素而没有numpy?【代码】

特定samplelist = [100,101,102,103,104,105,106,107,108,109]然后我想要输出如下:[100,[101,102,103,104,105,106,107,108,109]] [101,[100,102,103,104,105,106,107,108,109]] [102,[100,101,103,104,105,106,107,108,109]]注意:在第一行输出100被跳过,即列表[0],然后在第二行中跳过列表[1]. 我尝试弹出并显示,但下次没有显示弹出的所有元素.所以我尝试再次将它们插入循环结束但没有运气.解决方法:一个简单的列表理解应该这样做...

python – 一个numpy数组在更改另一个时意外更改,尽管它是分开的【代码】

我在我的大代码中发现了一个错误,我将问题简化为下面的案例. 虽然在每个步骤中我只改变w2,但是当我在每一步打印出w1时,它也会被改变,因为第一个循环的结束我将它们分配为相等.我读到了这个,但有写的以防万一我w1 = W2 [:]它会解决问题,但事实并非如此import numpy as np import mathw1=np.array([[1,2,3],[4,5,6],[7,8,9]]) w2=np.zeros_like(w1) print 'w1=',w1 for n in range(0,3):for i in range(0,3):for j in range(0,3):pr...

python – 将两个numpy.arrays之间的比较应用于一列,但检索整行【代码】

我有两个numpy数组,每个数组有两列.import numpy as npa = np.array([[1131, 1],[4131, 2],[421, 1],[41, 1]])b = np.array([[5881, 2],[637, 2],[742, 2],[36, 2]])我想创建一个包含2列的第三个数组,第一列将包含a和b的第一列之间的最小值,第二列将包含包含最小值的数组的第二列. 第三个数组应该是c = np.array([[1131, 1],[637, 2],[421, 1],[36, 2]])我怎样才能有效地做到这一点?解决方法:看起来像是在np.where内使...

python – 在块矩阵中排列numpy数组【代码】

我有3个numpy数组A,B和C.为简单起见,我们假设它们都是形状[n,n].我想将它们安排为块矩阵A B B^t C其中B ^ t表示B的转置.当然,我可以通过一系列连接来做到这一点top_row = np.concatenate([A, B], axis=1) bottom_row = np.concatenate([B.transpose(), C], axis=1) result = np.concatenate([top_row, bottom_row], axis=0)有更简单,更易读的方式吗?解决方法:从NumPy 1.13开始,有np.block.这是从嵌套的块列表构建矩阵,但它也更...

python – Matlab的缓冲区是否等同于numpy?【代码】

我看到有一个array_split并拆分methods但是当你必须拆分一个长度不是块大小的整数倍的数组时,这些都不是很方便.而且,这些方法输入是切片的数量而不是切片大小.我需要更像Matlab的buffer方法,它更适合信号处理. 例如,如果我想将信号缓冲到大小为60的块,我需要这样做:np.vstack(np.hsplit(x.iloc [0:((len(x)// 60)* 60)],len (x)// 60))这很麻烦.解决方法:我编写了以下例程来处理我需要的用例,但是我还没有实现/测试“underlap”....

Numpy-如何在OpenCV python中迭代创建ROI?【代码】

我试图在较小图像的网格中分割图像,以便我可以分别处理每个小图像.为此,我意识到我必须将每个小图像定义为ROI,我可以从那里轻松使用它. 现在,我的网格大小不固定.即,如果用户输入5,我必须制作55的网格. 逐像素迭代图像会很慢,所以我决定使用Numpy通过使用这个构造来创建ROI:#Assuming user entered grid size =5 roiwidth=w/5 roiheight=h/5 roi0=img[0:roiheight,0:roiwidth]这将是我的第一片. h和w分别是图像的高度和宽度.对...

将Python列表转换为Numpy数组InPlace【代码】

我有一个巨大的python列表(16 GB),我想将它转换为numpy数组,inplace.我承担不起这个说法huge_array = np.array(huge_list).astype(np.float16)我正在寻找一些有效的方法将这个huge_list转换为numpy数组,而不是复制它. 有谁能建议一个有效的方法来做到这一点?这可能涉及首先将列表保存到磁盘然后加载为numpy数组,我很好. 我非常感谢任何帮助. 编辑1:huge_list是一个内存python列表,它是在运行时创建的,所以它已经占用了16GB.我需要...

python – Numpy np.where多重条件【代码】

我需要使用numpy处理多个条件. 我正在尝试这个似乎有用的代码. 我的问题是:还有另一种可以做同样工作的替代方案吗?Mur=np.array([200,246,372])*pq.kN*pq.m Mumax=np.array([1400,600,700])*pq.kN*pq.m Mu=np.array([100,500,2000])*pq.kN*pq.m Acreq=np.where(Mu<Mur,0,"zero") Acreq=np.where(((Mur<Mu)&(Mu<Mumax)),45,Acreq) Acreq=np.where(Mu>Mumax,60,Acreq) Print(Acreq) ['0' '45' '60']解决方法:从这开始:Mur = np...

我可以在Android上运行Numpy(或其他Python包)吗?

参见英文答案 > Using python and matplotlib on android 12个我已经实现了一个python脚本,它导入Numpy和Pandas,我想在Android上运行这个脚本.更确切地说,我想将这个脚本嵌入到一个应用程序中.我想知道它是否可能?如果是这样,实施它的最佳做法是什么?我非常感谢任何帮助!解决方法:如果您不想构建网站或应用程序并将Python / Pandas作为后端运行.您可以在Android上使用Kivy作为packager to r...