计算图像的梯度矢量场

时间:2012-05-02 15:31:31

标签: image matlab image-processing

我想在图像中读取 - 圆形图片,并计算该图像的渐变矢量场(即向量均匀指向并垂直于圆圈的向量)。我的逻辑让我失望了一点,但我有:

clear all;
im = im2double(imread('littlecircle.png'));
im = double(im);
[nr,nc]=size(im);
[dx,dy] = gradient(im);
[x y] = meshgrid(1:nc,1:nr);
u = x;
v = y;
quiver(x,y,u,v)

如果我只是简单地执行上述操作,我会得到一个矢量场,但它只是空网格的渐变(即只是渐变的矢量场y = x)。我真正想要的是使用

[dx,dy] = gradient(im);

检测图像中圆的边缘,然后根据图像中的圆计算梯度矢量场。很明显,分配u = x和v = y只会给我一条直线的矢量场 - 所以基本上,我想把图像的渐变嵌入到矢量u和v中。我该怎么做?

my result

image that i am getting error with

1 个答案:

答案 0 :(得分:14)

你在代码中犯了一个错误(除此之外,它工作正常)。您应该替换以下内容:

u = dx;
v = dy;

u = x;
v = y;

它可以像this图片一样使用魅力!

编辑: 如果要在图像上超级施加矢量,请执行以下操作:

clear all;
im = imread('littlecircle.png');
[nr,nc]=size(im);
[dx,dy] = gradient(double(im));
[x y] = meshgrid(1:nc,1:nr);
u = dx;
v = dy;
imshow(im);
hold on
quiver(x,y,u,v)

请注意,我没有将im转换为double,因为imshow无法正确显示(需要uint8)。根据图像尺寸的不同,您可能需要放大才能看到渐变矢量。

您可以看到叠加在图像上的矢量的放大区域,如下所示:

Gradient vectors of a circle in an image

质量更好的图片位于http://i.stack.imgur.com/fQbwI.jpg