我的主管使用IDL，但我一直在使用Python，因为我对它更熟悉。我正在执行插值并保存了下限/上限值。有更快的方法吗？
变量
输入：
- sed = numpy array [6,221,6900]
- it0，it1，iz0，iz1 = numpy arrays [341499]
- 快照（其中38个）
输出：
- sed1，sed2，sed3，sed4 = numpy arrays [341499]
MWE
我想循环遍历38个快照，然后循环通过341499粒子，然后分配下面给出的结果numpy数组[6900]。
sed1 = sed[iz0, it0]
sed2 = sed[iz1, it0]
sed3 = sed[iz0, it1]
sed4 = sed[iz1, it1]
我尝试过什么
- 我无法初始化所需大小的数组，即numpy [38,341499,4,6900]，因为这会产生内存错误。意义不能使用向量[：]操作
进行分配
- 我尝试初始化一个大小为[38,341499]的numpy dtype对象数组，但这很慢