【Numpy,Python中的“拉伸”直方图(级别)】教程文章相关的互联网学习教程文章

python – 从NumPy数组到Mat(OpenCV)的C转换【代码】

我正在围绕ArUco增强现实库(基于OpenCV)编写一个薄的包装器.我想要构建的接口非常简单: > Python将图像传递给C代码;> C代码检测标记并将其位置和其他信息作为dicts元组返回给Python. 但是,我无法弄清楚如何在Python中表示图像以将其传递给C.对于GUI和相机管理,我将使用PyQt,所以最初它将是QImage,但我不能简单地将它传递给OpenCV(或者我可以?).起初,我尝试使用嵌套元组来表示每个像素的行,列和颜色,所以我最终得到了这个示例代码...

python – Py_initialize / Py_Finalize没有使用numpy两次【代码】

在下面的代码的第二次调用,我的应用程序段错误,所以我想我错过了一些东西:Py_Initialize(); pName = PyString_FromString("comp_macbeth"); pModule = PyImport_Import(pName); Py_DECREF(pName);if(pModule == NULL) {PyErr_Print();Py_Finalize();return; }pFunc = PyObject_GetAttrString(pModule, "compute"); /* pFunc is a new reference */if (!pFunc || !PyCallable_Check(pFunc) ) {PyErr_Print();Py_Finalize();return;...

使用Python和NumPy从矩阵中获取最小/最大n值和索引的有效方法【代码】

在给定NumPy矩阵(2D数组)的情况下,返回数组中最小/最大n值(及其索引)的有效方法是什么? 目前我有:def n_max(arr, n):res = [(0,(0,0))]*nfor y in xrange(len(arr)):for x in xrange(len(arr[y])):val = float(arr[y,x])el = (val,(y,x))i = bisect.bisect(res, el)if i > 0:res.insert(i, el)del res[0]return res这比pyopencv生成我要运行它的数组的图像模板匹配算法要长三倍,我认为这很愚蠢.解决方法:从另一个答案的时间开始,...

python – 将NumPy数组转换为set需要太长时间【代码】

我正在尝试执行以下操作from numpy import * x = array([[3,2,3],[711,4,104],.........,[4,4,782,7845]]) # large nparray for item in x:set(item)与以下相比需要很长时间:x = array([[3,2,3],[711,4,104],.........,[4,4,782,7845]]) # large nparray for item in x:item.tolist()为什么将NumPy数组转换为集合而不是列表需要更长的时间?我的意思是基本上都有复杂度O(n)?解决方法:TL; DR:set()函数使用Pythons迭代协议创建...

python – 解释为什么不应该从源目录导入numpy【代码】

免责声明: 我已经检查了以下其他StackOverflow问题: > How to import numpy in python shell> How can I use numpy without installing it?> Import a module from a relative path 也许对某些人来说,这些可能会回答我的问题,但据我所知,我仍然不了解情况. 我试图导入numpy以便matplotlib可以工作,但是在numpy文件夹中执行__init__.py文件时,会显示以下错误消息:ImportError: Error importing numpy: you should not try to imp...

连续的,重叠的数组子集(NumPy,Python)【代码】

我有一个NumPy阵列[1,2,3,4,5,6,7,8,9,10,11,12,13,14],并希望有一个像[[1,2,3, 4],[2,3,4,5],[3,4,5,6],…,[11,12,13,14]]. 当然这可以通过循环大数组并将长度为4的数组添加到新数组中,但我很好奇是否有一些秘密的’魔术’Python方法正在做这个:)解决方法:最快的方法似乎是预分配数组,在本答案的底部作为选项7给出.>>> import numpy as np >>> A=np.array([1,2,3,4,5,6,7,8,9,10,11,12,13,14]) >>> A array([ 1, 2, 3, 4, 5, ...

为什么在Python列表上的“for”比在Numpy数组上更快?【代码】

所以我没有讲一个很长的故事,我正在研究一些代码,我从二进制文件中读取一些数据,然后使用for循环遍历每一个点.所以我完成了代码,它的运行速度非常慢.我从大约128个数据通道循环了大约60,000个点,这需要一分钟或更长时间来处理.这比我预期的Python运行速度要慢.所以我通过使用Numpy使整个事情变得更有效率但是在试图弄清楚为什么原始进程运行得如此之慢以至于我们进行了一些类型检查并发现我在循环Numpy数组而不是Python列表.好吧没...

numpy数组数组的pythonic方法(带行索引)【代码】

我想在表中找到与特定索引相对应的值.例如,这是我的表:import numpy as np my_array = np.array([[0,1,0,1,0,1,0],[1,2,1,2,1,2,1],[4,5,4,3,3,4,5]])#--------------------------------------------------------------------- # my_array : [[0, 1, 0, 1, 0, 1, 0], # [1, 2, 1, 2, 1, 2, 1], # [4, 5, 4, 3, 3, 4, 5]])以下是一系列索引.此数组中的值是my_array的行. (列未编制索...

python – 沿动态指定的轴切割numpy数组【代码】

我想动态切片沿特定轴的numpy数组.鉴于这种:axis = 2 start = 5 end = 10我希望得到与此相同的结果:# m is some matrix m[:,:,5:10]使用这样的东西:slc = tuple(:,) * len(m.shape) slc[axis] = slice(start,end) m[slc]但是:值不能放在元组中,所以我无法弄清楚如何构建切片.解决方法:我认为一种方法是使用切片(无):>>> m = np.arange(2*3*5).reshape((2,3,5)) >>> axis, start, end = 2, 1, 3 >>> target = m[:, :, 1:3] >>>...

python – 使用最后的非零值填充1d numpy数组的零值【代码】

假设我们有一个填充了一些int值的1d numpy数组.让我们说其中一些是0. 有没有办法,使用numpy数组的幂,用找到的最后一个非零值填充所有0值? 例如:arr = np.array([1, 0, 0, 2, 0, 4, 6, 8, 0, 0, 0, 0, 2]) fill_zeros_with_last(arr) print arr[1 1 1 2 2 4 6 8 8 8 8 8 2]一种方法是使用此功能:def fill_zeros_with_last(arr):last_val = None # I don't really care about the initial valuefor i in range(arr.size):if arr[i...

python – Numpy和16位PGM【代码】

使用numpy在Python中读取16位PGM图像的有效且清晰的方法是什么? 我无法使用PIL加载16位PGM图像due to a PIL bug.我可以使用以下代码读取标题:dt = np.dtype([('type', 'a2'),('space_0', 'a1', ),('x', 'a3', ),('space_1', 'a1', ),('y', 'a3', ),('space_2', 'a1', ),('maxval', 'a5')]) header = np.fromfile( 'img.pgm', dtype=dt ) print header这打印出正确的数据:(‘P5′,”,’640′,”,’480′,”,’65535’)但我感觉这不...

python – NumPy的axis参数如何工作?【代码】

有人可以准确解释NumPy中轴参数的作用吗? 我非常困惑. 我正在尝试使用myArray.sum函数(axis = num) 起初我认为如果数组本身是3维,则axis = 0将返回三个元素,由同一位置中所有嵌套项的总和组成.如果每个维度包含五个维度,我希望axis = 1返回五个项目的结果,依此类推. 但事实并非如此,文档并没有很好地帮助我(they use a 3x3x3 array so it’s hard to tell what’s happening) 这是我做的:>>> e array([[[1, 0],[0, 0]],[[1, 1],[...

python – 在NumPy数组中查找等于零的元素索引【代码】

NumPy具有有效的函数/方法nonzero()来识别ndarray对象中的非零元素的索引.获取值为零的元素的索引的最有效方法是什么?解决方法:numpy.where()是我的最爱.>>> x = numpy.array([1,0,2,0,3,0,4,5,6,7,8]) >>> numpy.where(x == 0)[0] array([1, 3, 5])

python – numpy中的快速外张量积【代码】

我有两个numpy数组:x of shape ((d1,...,d_m)) y of shape ((e_1,...e_n)) 我想形成外张量产品,即numpy数组z of shape ((d1,...,d_m,e_1,...,e_n))这样的z[i_1,...,i_n,i_{n+1}...,i_{m+n}] == x[i_1,...i_m]*y[i_{m+1},...,i_{m+n}]我必须多次执行上面的外部乘法,所以我想尽可能加快速度.解决方法:外部的替代方法是明确扩展尺寸.对于1d阵列,这将是x[:,None]*y # y[None,:] is automatic.对于1010阵列,并推广尺寸扩展,我得到相...

python – Numpy将数组乘以矩阵(外积)【代码】

我有2个numpy形状阵列(5,1)说:A = [1,2,3,4,5]B = [2,4,2,3,6] 如何使矩阵乘以每个第i个元素与每个第j个?喜欢:..a = [1,2,3,4,5] b 2 2, 4, 6, 8,10 4 4, 8,12,16,20 2 2, 4, 6, 8,10 3 3, 6, 9,12,15 6 6,12,18,24,30不使用forloops?我可以使用重塑,缩减或乘法的任何组合吗? 现在我沿着行和沿着colums创建每个数组的a * b平铺,然后乘以元素,但在我看来必须有一个更简单的方法.解决方法:使用numpy.outer()和n...