我想读入一幅图像——一张圆的图片,并计算该图像的梯度矢量场(即矢量均匀地指向圆并垂直于圆)。我的逻辑有点让我失望,但我有:
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 中。我该怎么做?
最佳答案
你在代码中犯了一个错误(除此之外,它工作正常)。您应该替换以下内容:
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)。根据您的图像尺寸,您可能想要放大以查看渐变向量。
您可以看到叠加在图像上的矢量放大区域,如下所示:
质量更好的图片位于 http://i.stack.imgur.com/fQbwI.jpg
关于image - 计算图像的梯度矢量场,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/10416918/