numpy.nditer#

class numpy.nditer(op, flags=None, op_flags=None, op_dtypes=None, order='K', casting='safe', op_axes=None, itershape=None, buffersize=0)[源代码]#

用于迭代数组的高效多维迭代器对象。要开始使用此对象,请参阅数组迭代入门指南

参数:
opndarray 或 array_like 序列

要迭代的数组。

flags字符串序列,可选

控制迭代器行为的标志。

  • buffered 在需要时启用缓冲。

  • c_index 导致跟踪 C 顺序的索引。

  • f_index 导致跟踪 Fortran 顺序的索引。

  • multi_index 导致跟踪多重索引,或每个迭代维度一个索引的元组。

  • common_dtype 导致所有操作数转换为公共数据类型,并在必要时进行复制或缓冲。

  • copy_if_overlap 导致迭代器确定读取操作数是否与写入操作数重叠,并在必要时进行临时复制以避免重叠。在某些情况下,可能会出现假阳性(不必要的复制)。

  • delay_bufalloc 延迟缓冲区的分配,直到调用 reset() 函数。允许在将值复制到缓冲区之前初始化 allocate 操作数。

  • external_loop 导致给定的 values 成为具有多个值的一维数组,而不是零维数组。

  • grow_inner 当同时使用 bufferedexternal_loop 时,允许 value 数组的大小大于缓冲区大小。

  • ranged 允许将迭代器限制为 iterindex 值的子范围。

  • refs_ok 启用引用类型的迭代,例如对象数组。

  • reduce_ok 启用广播的 readwrite 操作数的迭代,也称为归约操作数。

  • zerosize_ok 允许 itersize 为零。

op_flags字符串列表的列表,可选

这是每个操作数的标志列表。至少必须指定 readonlyreadwritewriteonly 中的一个。

  • readonly 表示操作数只能从中读取。

  • readwrite 表示操作数将从中读取和写入。

  • writeonly 表示操作数只能写入。

  • no_broadcast 防止操作数被广播。

  • contig 强制操作数数据是连续的。

  • aligned 强制操作数数据是对齐的。

  • nbo 强制操作数数据为本机字节顺序。

  • copy 如果需要,允许进行临时的只读复制。

  • updateifcopy 如果需要,允许进行临时的读写复制。

  • allocate 如果 op 参数中为 None,则导致分配数组。

  • no_subtype 防止 allocate 操作数使用子类型。

  • arraymask 表示此操作数是用于在写入设置了 “writemasked” 标志的操作数时选择元素的掩码。迭代器不强制执行此操作,但是当从缓冲区写回数组时,它只复制此掩码指示的元素。

  • writemasked 表示只写入所选 arraymask 操作数为 True 的元素。

  • overlap_assume_elementwise 可用于标记仅按迭代器顺序访问的操作数,以便在存在 copy_if_overlap 时允许进行不太保守的复制。

op_dtypesdtype 或 dtype 元组,可选

操作数所需的数据类型。如果启用了复制或缓冲,则数据将转换为/从其原始类型转换。

order{‘C’, ‘F’, ‘A’, ‘K’},可选

控制迭代顺序。“C” 表示 C 顺序,“F” 表示 Fortran 顺序,“A” 表示如果所有数组都是 Fortran 连续的则为 ‘F’ 顺序,否则为 ‘C’ 顺序,而 ‘K’ 表示尽可能接近数组元素在内存中出现的顺序。这也影响 allocate 操作数的元素内存顺序,因为它们被分配为与迭代顺序兼容。默认为 ‘K’。

casting{‘no’, ‘equiv’, ‘safe’, ‘same_kind’, ‘unsafe’},可选

控制在复制或缓冲时可能发生的数据类型转换类型。不建议将其设置为 “unsafe”,因为它会对累积产生不利影响。

  • “no” 表示根本不应该转换数据类型。

  • “equiv” 表示只允许字节顺序更改。

  • “safe” 表示只允许可以保留值的转换。

  • “same_kind” 表示只允许安全转换或同类型内的转换,如 float64 到 float32。

  • “unsafe” 表示可以进行任何数据转换。

op_axes整数列表的列表,可选

如果提供,则为每个操作数的整数列表或 None。操作数的轴列表是从迭代器的维度到操作数的维度的映射。可以为条目放置 -1 值,导致将该维度视为 newaxis

itershape整数元组,可选

迭代器所需的形状。这允许通过 op_axes 映射的维度与不同操作数的维度不对应的 allocate 操作数,在该维度上的值不等于 1。

buffersizeint, 可选

当启用缓冲时,控制临时缓冲区的大小。设置为 0 表示默认值。

说明

nditer 取代了 flatiternditer 背后的迭代器实现也通过 NumPy C API 公开。

Python 公开提供了两个迭代接口,一个遵循 Python 迭代器协议,另一个镜像 C 风格的 do-while 模式。在大多数情况下,原生的 Python 方法更好,但是如果您需要迭代器的坐标或索引,请使用 C 风格的模式。

示例

下面是如何使用 Python 迭代器协议编写一个 iter_add 函数的示例

>>> import numpy as np
>>> def iter_add_py(x, y, out=None):
...     addop = np.add
...     it = np.nditer([x, y, out], [],
...                 [['readonly'], ['readonly'], ['writeonly','allocate']])
...     with it:
...         for (a, b, c) in it:
...             addop(a, b, out=c)
...         return it.operands[2]

这是相同的函数,但遵循 C 风格的模式

>>> def iter_add(x, y, out=None):
...    addop = np.add
...    it = np.nditer([x, y, out], [],
...                [['readonly'], ['readonly'], ['writeonly','allocate']])
...    with it:
...        while not it.finished:
...            addop(it[0], it[1], out=it[2])
...            it.iternext()
...        return it.operands[2]

这是一个外部乘积函数的示例

>>> def outer_it(x, y, out=None):
...     mulop = np.multiply
...     it = np.nditer([x, y, out], ['external_loop'],
...             [['readonly'], ['readonly'], ['writeonly', 'allocate']],
...             op_axes=[list(range(x.ndim)) + [-1] * y.ndim,
...                      [-1] * x.ndim + list(range(y.ndim)),
...                      None])
...     with it:
...         for (a, b, c) in it:
...             mulop(a, b, out=c)
...         return it.operands[2]
>>> a = np.arange(2)+1
>>> b = np.arange(3)+1
>>> outer_it(a,b)
array([[1, 2, 3],
       [2, 4, 6]])

这是一个类似于 “lambda” ufunc 的函数示例

>>> def luf(lamdaexpr, *args, **kwargs):
...    '''luf(lambdaexpr, op1, ..., opn, out=None, order='K', casting='safe', buffersize=0)'''
...    nargs = len(args)
...    op = (kwargs.get('out',None),) + args
...    it = np.nditer(op, ['buffered','external_loop'],
...            [['writeonly','allocate','no_broadcast']] +
...                            [['readonly','nbo','aligned']]*nargs,
...            order=kwargs.get('order','K'),
...            casting=kwargs.get('casting','safe'),
...            buffersize=kwargs.get('buffersize',0))
...    while not it.finished:
...        it[0] = lamdaexpr(*it[1:])
...        it.iternext()
...    return it.operands[0]
>>> a = np.arange(5)
>>> b = np.ones(5)
>>> luf(lambda i,j:i*i + j/2, a, b)
array([  0.5,   1.5,   4.5,   9.5,  16.5])

如果使用操作数标志 "writeonly""readwrite",则操作数可能是具有 WRITEBACKIFCOPY 标志的原始数据的视图。在这种情况下,必须将 nditer 用作上下文管理器,或者必须在使用结果之前调用 nditer.close 方法。当调用 __exit__ 函数时,临时数据将写回到原始数据,而不是在此之前。

>>> a = np.arange(6, dtype='i4')[::-2]
>>> with np.nditer(a, [],
...        [['writeonly', 'updateifcopy']],
...        casting='unsafe',
...        op_dtypes=[np.dtype('f4')]) as i:
...    x = i.operands[0]
...    x[:] = [-1, -2, -3]
...    # a still unchanged here
>>> a, x
(array([-1, -2, -3], dtype=int32), array([-1., -2., -3.], dtype=float32))

需要注意的是,一旦退出迭代器,悬空引用(如示例中的 x)可能与也可能不与原始数据 a 共享数据。如果回写语义处于活动状态,即如果 x.base.flags.writebackifcopyTrue,则退出迭代器将断开 xa 之间的连接,写入 x 将不再写入 a。如果回写语义未处于活动状态,则 x.data 仍将指向 a.data 的某个部分,写入其中一个将影响另一个。

上下文管理和 close 方法出现在 1.15.0 版本中。

属性:
dtypesdtype(s) 元组

value 中提供的值的数据类型。如果启用了缓冲,则可能与操作数数据类型不同。仅在迭代器关闭之前有效。

finishedbool

是否完成了对操作数的迭代。

has_delayed_bufallocbool

如果为 True,则迭代器是使用 delay_bufalloc 标志创建的,并且尚未在其上调用 reset() 函数。

has_indexbool

如果为 True,则迭代器是使用 c_indexf_index 标志创建的,并且可以使用属性 index 来检索它。

has_multi_indexbool

如果为 True,则迭代器是使用 multi_index 标志创建的,并且可以使用属性 multi_index 来检索它。

index

当使用 c_indexf_index 标志时,此属性提供对索引的访问。如果访问且 has_index 为 False,则引发 ValueError。

iterationneedsapibool

迭代是否需要访问 Python API,例如,如果其中一个操作数是对象数组。

iterindexint

与迭代顺序匹配的索引。

itersizeint

迭代器的大小。

itviews

operands 在内存中的结构化视图,与重新排序和优化的迭代器访问模式匹配。仅在迭代器关闭之前有效。

multi_index

当使用 multi_index 标志时,此属性提供对索引的访问。如果访问且 has_multi_index 为 False,则引发 ValueError。

ndimint

迭代器的维度。

nopint

迭代器操作数的数量。

operands操作数元组

operands[切片]

shape整数元组

形状元组,迭代器的形状。

value

当前迭代中 operands 的值。通常,这是一个数组标量的元组,但是如果使用标志 external_loop,则它是一个一维数组的元组。

方法

close()

解析可写操作数中的所有回写语义。

copy()

获取迭代器当前状态的副本。

debug_print()

nditer 实例的当前状态和调试信息打印到标准输出。

enable_external_loop()

当在构造期间未使用 "external_loop",但需要使用时,此方法会修改迭代器的行为,就像指定了该标志一样。

iternext()

检查是否还剩余迭代,并执行单次内部迭代而不返回结果。

remove_axis(i, /)

从迭代器中移除轴 i

remove_multi_index()

当指定 "multi_index" 标志时,此方法会将其删除,从而允许进一步优化内部迭代结构。

reset()

将迭代器重置为其初始状态。