Oblicz gradientowe pole wektorowe obrazu

Question

Chcę odczytać obraz - obrazek koła i obliczyć gradientowe pole wektorowe tego obrazu (tj. Wektory wskazujące równomiernie i normalnie na okrąg). Moja logika zawodzi mnie trochę, ale mam:

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) 

gdybym po prostu nie wyżej, mam pole wektorowe, ale to jest po prostu gradient pustej siatki (czyli tylko pola wektorowego gradientu y = x). Co faktycznie chcę jest użycie

[dx,dy] = gradient(im); 

do wykrywania krawędzi okręgu na obrazie, a następnie obliczyć pole wektorowe gradientu ze względu na kole w obrazie. oczywiście, przypisanie u = x i v = y da mi tylko pole wektorowe linii prostej - tak w zasadzie chcę osadzić gradient obrazu w wektorach u i v. Jak to zrobić?

Answer 1

Zrobiłeś błąd w kodzie (innego niż, że to działa dobrze). należy wymienić następujące:

u = dx; 
v = dy; 

nie

u = x; 
v = y; 

Współpracuje z this obrazu jak urok!

EDIT: Jeśli chcesz super nałożyć wektory na obrazie, a następnie wykonaj następujące czynności:

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) 

Zauważ, że nie przekonwertować im podwoić, ponieważ nie wydaje się poprawnie imshow (wymaga uint8). W zależności od wymiarów obrazu możesz powiększyć widok, aby zobaczyć wektory gradacji.

można zobaczyć powiększony obszar wektorów nałożony na obraz poniżej:

Lepsza jakość obrazu jest na http://i.stack.imgur.com/fQbwI.jpg