我有一个形状为40x40的测试矩阵(z),填充了零。
我需要在测试矩阵的特定位置添加4个形状子矩阵,称为c1,c2(5x5),c3(7x7)和c4(9x9)。
我想把子矩阵的中心放在相应的位置,然后简单地执行元素的加法。测试矩阵中的位置是:z(9,9),z(9,29),z(29,9),z(29,29)。
我试着看这些线程,但我不能得到一个明确的答案,如何解决我的问题。How to add different arrays from the center point of an array in Python/NumPyAdding different sized/shaped displaced NumPy matrices
我尝试的代码示例:
def zero_matrix(d):
matrix = np.zeros((d,d), dtype=np.float)
return matrix
z = zero_matrix(40)
c1 = np.genfromtxt('xxxxxx', dtype=None, delimiter = '\t')
c2 = np.genfromtxt('xxxxxx', dtype=None, delimiter = '\t')
c3 = np.genfromtxt('xxxxxx', dtype=None, delimiter = '\t')
c4 = np.genfromtxt('xxxxxx', dtype=None, delimiter = '\t')
def adding(z):
for i in range(z.shape[0]):
for j in range(z.shape[1]):
if i == 9 and j==9:
c1mid = c1.shape[0]//2
z[i,j] = c1[c1mid,c1mid]
print z
return z
字符串
但这只是增加了中心,而不是整个子矩阵。
它应该看起来像这样:
2条答案
按热度按时间dsf9zpds1#
numpy中数组切片的好处是你不需要使用for循环。它只放置中心元素的原因是因为你只放置一个元素(c1[c1mid,c1mid]是一个数字),下面是你可以做的:
字符串
uqjltbpv2#
内置的
np.ix_
工作完美字符串
输出是
型