回答此问题可获得 20 贡献值,回答如果被采纳可获得 50 分。
<p>我需要与Matlab(八进制)离散Laplacian运算符(函数)del2()等价的Python/Numpy。我尝试了一些Python解决方案,但没有一个与del2的输出相匹配。我有八度音阶</p>
<pre><code>image = [3 4 6 7; 8 9 10 11; 12 13 14 15;16 17 18 19]
del2(image)
</code></pre>
<p>这就产生了结果</p>
<pre><code> 0.25000 -0.25000 -0.25000 -0.75000
-0.25000 -0.25000 0.00000 0.00000
0.00000 0.00000 0.00000 0.00000
0.25000 0.25000 0.00000 0.00000
</code></pre>
<p>在Python上我试过</p>
<pre><code>import numpy as np
from scipy import ndimage
import scipy.ndimage.filters
image = np.array([[3, 4, 6, 7],[8, 9, 10, 11],[12, 13, 14, 15],[16, 17, 18, 19]])
stencil = np.array([[0, 1, 0],[1, -4, 1], [0, 1, 0]])
print ndimage.convolve(image, stencil, mode='wrap')
</code></pre>
<p>结果是</p>
<pre><code>[[ 23 19 15 11]
[ 3 -1 0 -4]
[ 4 0 0 -4]
[-13 -17 -16 -20]]
</code></pre>
<p>我也试过了</p>
<pre><code>scipy.ndimage.filters.laplace(image)
</code></pre>
<p>结果就是这样</p>
<pre><code>[[ 6 6 3 3]
[ 0 -1 0 -1]
[ 1 0 0 -1]
[-3 -4 -4 -5]]
</code></pre>
<p>所以所有的输出似乎都不匹配。八度码del2.m表明它是拉普拉斯算子。我遗漏了什么吗?</p>