numpy

    0热度

    1回答

    欲索引通过多维阵列这样numpy的索引: a = range(12).reshape(3, 2, 2) def fun(axis, state): # if axis=0 return a[state, :, :] # if axis=1 it should return a[:, state, :] 样本输出: fun(0, 1) array([[4, 5]

    0热度

    1回答

    我无法使用brentq函数迭代数组中的每个元素。下面定义的函数中的q是一个FITS文件数组,我们使用该数组中的每个元素作为输入来运行brentq函数,以便解决T。 本质上,我的问题在于不是特别知道在何处或如何实现合适的for循环来迭代函数遍历每个元素q。 有关如何解决此问题的任何建议? def f(T,q,coeff1,coeff2,coeff3): return q*const3 -

    1热度

    1回答

    我写了一个函数,它接受一组随机笛卡尔坐标并返回保留在某个空间域内的子集。为了说明: grid = np.ones((5,5)) grid = np.lib.pad(grid, ((10,10), (10,10)), 'constant') >> np.shape(grid) (25, 25) random_pts = np.random.random(size=(100, 2)) *

    2热度

    3回答

    我有一个图像我加载到python使用matplotlib.pyplot.imread其结果作为numpy数组包含rgb值的数组。这里是除了两个像素的白色哑片段: >>> test array([[[255, 255, 255], [255, 255, 255]], [[ 1, 255, 255], [255, 255, 255]], [[255

    4热度

    2回答

    我在学习numpy,对广播有点困惑,这里是我的设置。我有两个矩阵 >>> y=np.array([1,2,3]) >>> v = np.array([1,2,3]) >>> r=np.reshape(v, (3, 1)) 因此r是(3 * 1)矩阵,而y是形状为(3,)的秩为1的矩阵。如果我做了y.dot(r),我得到了14,假设numpy在y上应用广播,使它成为(1 * 3),然后用r(

    1热度

    2回答

    什么是numpy的最好的方式,以产生与以下形式的n个值的随机阵列 arr = [1,-1,-1,1,1,1,...] 该平均掉尽可能接近到预定的值m,以便 print 1/n*np.sum(arr) >>> #value that is as close as possible to m 我一直在试验 numpy.random.choice([-1,1], size=n) 但

    1热度

    3回答

    我对Python和Numpy相当陌生,当我将一个MATLAB程序翻译成Python时遇到了这个问题。 据我所知,下面的代码通过修改一个全局变量来表现异常,即使它不应该。 import numpy as np A = np.matrix([[0, 1, 2], [3, 4, 5], [6, 7, 8]]) def function(B): B[1,1] = B[1,1]/2

    -2热度

    1回答

    我的数据集看起来 1 2 3 4 5 ... 我有一个中间步骤,应该做如下因素 1 1,2 1,2,3 1,2,3,4 1,2,3,4,5 .... 最后计算出其平均 1 1.5 2 2.5 3 ... 问题 一)有没有办法在python/py-spark中实现这个? b)有没有开箱即用的方法/ api? c)我搜索这种解决方案最接近的是移动平均/滚动平均

    4热度

    1回答

    我试图运行这个setup.py文件,该文件位于.tat.gz文件中。我要去的目录它在运行蟒蛇setup.py但收到此错误信息: [email protected]:~/Downloads/numpy-1.9.2$ python setup.py Running from numpy source directory. /usr/lib/python2.7/distutils/dist.py:2

    2热度

    2回答

    我有一个巨大的一组目录有不同的列和每个列不同的标题名称,每个标题名称的描述作为注释给出我的ascii文件连续开始。将它们读入pandas.DataFrame的最佳方法是什么,同时它可以设置列的名称,而无需从头开始定义它。以下是我的目录的例子: # 1 MAG_AUTO Kron-like elliptical aperture magnitude [mag] # 2 rh hal