Question

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

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

Answer 1

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

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

>>> 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

如果你确实想要通过argmin获取位置，那么这有点棘手，但一种方法是使用屏蔽数组：

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

请注意，此处的条件是翻转的，因为对于排除的值，掩码为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

如果需要位置，请使用相应的蒙版数组变体

。 [编辑使用指数代替mgrid;我实际上已经忘记了它，直到它今天用于另一个答案！]

如果我仍然错了：^）这也不是你想要的，请编辑你的问题以包括你所需输入和输出的3x3示例。

Answer 2

我猜这是你想要实现的目标：

Argmin with arrays：

>>> from numpy import *
>>> a = array( [2,3,4] )
>>> argmin(a)
0
>>> print a[argmin(a)]
2

带有矩阵的Argmin：

>>> b=array( [[6,5,4],[3,2,1]] )
>>> argmin(b)
5
>>> print b[argmin(b)]
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
IndexError: index out of bounds

索引的相同方法对数组不起作用。原因是argmin（以及argmax）返回变量的索引 - 在矩阵的情况下，您需要将n维矩阵转换为1维索引数组。

要执行此操作，您需要致电ravel：

>>> print b
[[6 5 4]
 [3 2 1]]
>>> ravel(b)
array([6, 5, 4, 3, 2, 1])

当您将ravel与argmin结合使用时，您必须写下：

>>> print ravel(b)[argmin(b)]

numpy argmin优雅的解决方案。

2 个答案: