python画矢量场_Python中的图像渐变矢量场

我试图使用

Python获取

Gradient Vector Field的图像(类似于

this matlab question).

这是原始图片：

这是我的代码：

import numpy as np

import matplotlib.pyplot as plt

import Image

from PIL import ImageFilter

I = Image.open('test.png').transpose(Image.FLIP_TOP_BOTTOM)

I = I.filter(ImageFilter.BLUR)

p = np.asarray(I)

w,h = I.size

y, x = np.mgrid[0:h:500j, 0:w:500j]

dy, dx = np.gradient(p)

skip = (slice(None, None, 3), slice(None, None, 3))

fig, ax = plt.subplots()

im = ax.imshow(I, extent=[x.min(), x.max(), y.min(), y.max()])

ax.quiver(x[skip], y[skip], dx[skip], dy[skip])

ax.set(aspect=1, title='Quiver Plot')

plt.show()

这是结果：

问题是向量似乎是不正确的.放大图像时,此点变得更加清晰：

为什么有些向量按预期指向中心,而其他向量则没有？

也许调用np.gradient的结果有问题？