我在2D阵列上的恒定大小的移动窗口上应用操作.是否有一个有效的类似矢量化的操作,我可以实现这样做而无需在Python中循环?我目前的结构看起来像这样
for i in range(1,xmax-1):
for j in range(1,ymax-1):
out[i][j] = f(in[i][j],in[i+1][j],in[i-1][j],in[i][j+1],in[i][j-1],...)
Run Code Online (Sandbox Code Playgroud)
这些评论可以吃留在这个问题暗指矢量化这种操作这种可能性,但没有进一步的细节矢量索引/切片在numpy的/ SciPy的?
我有一个包含东距 (x)、北距 (y) 和高程数据 (z) 的文本文件,如下所示:
x y z
241736.69 3841916.11 132.05
241736.69 3841877.89 138.76
241736.69 3841839.67 142.89
241736.69 3841801.45 148.24
241736.69 3841763.23 157.92
241736.69 3841725.02 165.01
241736.69 3841686.80 171.86
241736.69 3841648.58 178.80
241736.69 3841610.36 185.26
241736.69 3841572.14 189.06
241736.69 3841533.92 191.28
241736.69 3841495.71 193.27
241736.69 3841457.49 193.15
241736.69 3841419.27 194.85
241736.69 3841381.05 192.31
241736.69 3841342.83 188.73
241736.69 3841304.61 183.68
241736.69 3841266.39 176.97
241736.69 3841228.18 160.83
241736.69 3841189.96 145.69
241736.69 3841151.74 129.09
241736.69 3841113.52 120.03
241736.69 3841075.30 111.84
241736.69 3841037.08 …Run Code Online (Sandbox Code Playgroud)