是否存在与numpy.pad()函数相反/相反的函数?

2024-05-15 14:49:51 发布

您现在位置:Python中文网/ 问答频道 /正文

是否有一个函数做了与numpy.pad()相反的事情?

我要找的是一个函数(均匀地)在每个方向上减少numpy数组(矩阵)的维数。我试图用负值调用numpy.pad(),但它给出了一个错误:

import numpy as np

A_flat = np.array([0,1,2,3,4,5,6,7,8,9,10,11])
A = np.reshape(A_flat, (3,2,-1))

#this WORKS:
B = np.pad(A, ((1,1),(1,1),(1,1)), mode='constant')

# this DOES NOT WORK:
C = np.pad(B, ((-1,1),(1,1),(1,1)), mode='constant')

错误:ValueError: ((-1, 1), (1, 1), (1, 1)) cannot contain negative values.

我知道这个函数numpy.pad()不接受负值,但是否有numpy.unpad()或类似的值?


Tags: 函数numpymode错误np矩阵数组this
3条回答

正如mdurant所建议的,只需使用切片索引:

In [59]: B[1:-1, 1:-1, 1:-1]
Out[59]: 
array([[[ 0,  1],
        [ 2,  3]],

       [[ 4,  5],
        [ 6,  7]],

       [[ 8,  9],
        [10, 11]]])

您想要的操作:

C = np.pad(B, ((-1,1),(1,1),(1,1)), mode='constant')

可替换为pad和常规切片的组合:

C = np.pad(B, ((0,1),(1,1),(1,1)), mode='constant')[1:,...]

一般的解决办法是:

def unpad(x, pad_width):
    slices = []
    for c in pad_width:
        e = None if c[1] == 0 else -c[1]
        slices.append(slice(c[0], e))
    return x[tuple(slices)]

# Test
import numpy as np
pad_width = ((0, 0), (1, 0), (3, 4))

a = np.random.rand(10, 10, 10)
b = np.pad(a, pad_width, mode='constant')
c = unpad(b, pad_width)
np.testing.assert_allclose(a, c)

相关问题 更多 >