numpy argmin优雅的解决方案.

Jus*_*ime 3 python numpy scipy

在python中找到我使用的数组最小值的索引y = numpy.argmin(someMat)

我能找到这个矩阵的最小值,使它不能以一种整齐的方式位于指定的范围内吗?

DSM*_*DSM 7

"我可以找到这个矩阵的最小值,使它不会以一种整齐的方式位于指定的范围内吗?"

如果您只关心满足某些条件而不是位置的最小值,那么

>>> numpy.random.seed(1)
>>> m = numpy.random.randn(5.,5.)
>>> m
array([[ 1.62434536, -0.61175641, -0.52817175, -1.07296862,  0.86540763],
       [-2.3015387 ,  1.74481176, -0.7612069 ,  0.3190391 , -0.24937038],
       [ 1.46210794, -2.06014071, -0.3224172 , -0.38405435,  1.13376944],
       [-1.09989127, -0.17242821, -0.87785842,  0.04221375,  0.58281521],
       [-1.10061918,  1.14472371,  0.90159072,  0.50249434,  0.90085595]])
>>> m[~ ((m < 0.5) | (m > 0.8))].min()
0.50249433890186823
Run Code Online (Sandbox Code Playgroud)

如果你确实想通过argmin获取位置,那么这有点棘手,但一种方法是使用掩码数组:

>>> numpy.ma.array(m,mask=((m<0.5) | (m > 0.8))).argmin()
23
>>> m.flat[23]
0.50249433890186823
Run Code Online (Sandbox Code Playgroud)

请注意,此处的条件是翻转的,因为对于排除的值,掩码为True,而不是包含的值.


更新:似乎"在指定范围内"并不意味着最小值不在某些范围内,而是您希望基于x,y坐标从搜索中排除矩阵的部分.这是一种方式(与以前相同的矩阵):

>>> xx, yy = numpy.indices(m.shape)
>>> points = ((xx == 0) & (yy == 0)) | ((xx > 2) & (yy < 3))
>>> points
array([[ True, False, False, False, False],
       [False, False, False, False, False],
       [False, False, False, False, False],
       [ True,  True,  True, False, False],
       [ True,  True,  True, False, False]], dtype=bool)
>>> m[points]
array([ 1.62434536, -1.09989127, -0.17242821, -0.87785842, -1.10061918,
        1.14472371,  0.90159072])
>>> m[points].min()
-1.1006191772129212
Run Code Online (Sandbox Code Playgroud)

如果需要位置,请使用相应的蒙版数组变体.[编辑使用指数代替mgrid; 我实际上已经忘记了它,直到它今天用于另一个答案!]

如果我仍然错了:^)这也不是你想要的,请编辑你的问题以包括你想要的输入和输出的3x3例子.